A numerical algorithm for studying subsonic chemically reacting flow in the presence of laser radiation

The article presents a numerical algorithm for studying subsonic viscous chemically active flows in the presence of laser radiation. The process model is described in the Navier – Stokes approximation adjusted for the subsonic flow regime, with addition of source terms corresponding to chemical transformations. An additional ordinary differential equation describing the propagation of laser radiation along the length of the region under study is introduced as well. The computational algorithm is based on splitting by physical processes. This makes it possible to calculate separately changes in concentrations during chemical transformations, convective fluxes, dissipative terms, dynamic pressure deviation and propagation of laser radiation. To account for the dissipative terms (diffusion, viscosity, and thermal conductivity), the local iteration method based on Chebyshev polynomials’ ordering. Due to the possible use of a larger total calculation time step, the software implementation of the constructed algorithm reveals shorter calculation times using the local iteration method for calculating dissipative terms in comparison with the algorithm calculating them based on a scheme with central differences. The algorithm was verified using the example of methane conversion by comparing it with the calculation of the stoichiometric balance of the brutto-reaction, as well as by studying the convergence of the solution on a sequence of thickening grids. Based on the developed algorithm, a numerical study of non-oxidative conversion of methane under the influence of laser radiation in a circular tube was carried out, and graphs of the distribution of the main characteristics of the mixture were obtained.