GM_Julia
is a zero-dimensional chemistry reaction solver for low-temperature
plasmas. It takes a list of reactions, some details on the plasma species, and
evolves them in time. These systems can be computationally intensive due to the
number of reactions and their disparate timescales: typical chemistry sets in
GM_Julia
are around 600 reaction equations, with rate coefficients spanning
many orders of magnitude. This makes for very numerically stiff systems,
requiring specialised solvers. The original author made the decision to use
Julia in order to take advantage of OrdinaryDiffEq.jl’s Rosenbrock
integrators, designed for stiff systems.