A Parallel Algorithm for a Two-Phase Gas-Solid-Particle Model with Chemical Reactions and Laser Radiation

The article is devoted to the development of a parallel computational algorithm for the model of a two-phase chemically active medium with laser radiation. The model consists of a system of low-Mach-number approximations of Navier–Stokes equations for multicomponent reacting mixtures, supplemented by equations of chemical kinetics, an equation for the radiation intensity, and equations for the solid phase. The solid phase is described in terms of particle concentrations. We describe them employing systems of convection-diffusion-reaction equations supplemented with thermal terms from laser radiation and reaction terms for compounds on the surface. The general solution algorithm is based on the splitting scheme for physical processes: chemical reactions, convection-diffusion, and stationary process for pressure correction. We use the WENO scheme to construct the computational algorithm for approximating the convective terms in the Navier–Stokes equations. Moreover, we use the RADAU5 plug-in module to calculate the equations of chemical kinetics and the equations for the radiation intensity and the temperature of the solid phase. The time derivatives are approximated explicitly. This approach guarantees the efficiency of using the suggested difference scheme for parallel implementation. The parallel algorithm is based on the principles of geometric parallelism; the MPI standard is used for interprocessor interaction. To obtain information about the efficiency of the parallel algorithm, we performed calculations on a sequence of shredding grids with different numbers of computing nodes. We applied the parallel code to compute the flow of a two-phase gas-solid medium with chemical reactions of light hydrocarbons in the gas phase and on the particle surface under the influence of laser radiation.