Comparison of explicit and implicit difference schemes of calculations for chemical non-equilibrium processes in nozzles

Fluid, gas and plasma mechanics


Krioukov V. G.*, Abdullin A. L.**, Demin A. V.***, Safiullin I. I.****

Kazan National Research Technical University named after A.N. Tupolev, KNRTU-KAI, 10, Karl Marks str., Kazan, 420111, Russia



The calculation of chemically non-equilibrium flows in the nozzles of rocket engines is normally performed using numerical implicit schemes (due to the stiffness of chemical kinetics equations). But the progress in developing of the stable explicit methods creates the possibility to use these simple methods instead of the implicit schemes. In this study, we introduce a calculating method of potential number of integration steps for explicit schemes, and their effectiveness is evaluated. The method includes:

-calculation of a chemically non-equilibrium flow along the nozzle length using the implicit scheme;

-parallel determination of the Jacobi matrix eigenvalues;

-calculation of the number of potential steps for the explicit integration scheme based on the stability bound.

The calculation of flows was carried out within the framework of the inverse nozzle problem, using the implicitly differential scheme of Pirumov U. G. Numerical research was carried out for the combustion products:

- liquid propellants (O2 + kerosene; N2O4 + C2H8N2) for Laval nozzles at: excess oxidant ratio; pressure Poc = 20…100 atm; minimum radius rm = 0.006…0.06 m and geometric expansion fa = 53.0;

- solid propellants (metalized fuel - C10.8760H46.546O25.806AL9.665CL1.517N6.781; nitrocellulose fuel - C23.498H30.259O34.190N10.011) for Laval nozzles at pressure Poc = 20…70 atm; minimum radius rm = 0.005…0.05 m and geometric expansion fa = 33.9.

The number of steps of the explicit scheme (K1) was calculated for the Runge-Kutta method of the 4th order. For the reactive media of liquid propellants, a huge number of potential steps (К1 ≈ 109) were obtained at high Poc and rm values. However, with a decrease of the Poc and rm parameters, K1 also decreased (to about К1 ≈ 107). In the subsonic part of the nozzle, the K1 values were approximately 10 times higher than in the supersonic part. For the reactive media of solid propellants, the results show the same trend, but at a lower level of K1 values, especially in the case of nitrocellulose fuel, when max (K1) ≈ 106.


chemical non-equilibrium flows, engine nozzles, mathematical model, eigenvalues


  1. Gidaspov V.YU. Trudy MAI, 2015, no. 83:

  2. Gidaspov V.YU. Trudy MAI, 2013, no. 66:

  3. Kartovitskii L.L., Levin V.M., Yakovlev A.A. Aviatsionnaya tekhnika, 2015, no. 4, pp. 67-72.

  4. Dregalin A.F., Barysheva O.B, Cherenkov A.S. Aviatsionnaya tekhnika, 2007, no. 3, pp. 46-49.

  5. L. Boccaletto, J.P. Dussauge. High-Performance Rocket Nozzle Concept. Journal of Propulsion and Power, 2010, vol. 26, no. 5, pp. 969-979.

  6. S.Keshav, Y.G. Utkin, M. Nishihara, J. W. Rich, I. V. Adamovich, A Bao. Studies of Chemi-Ionization and Chemiluminescence in Supersonic Flows of Combustion Products. Journal of Thermophysics and Heat Transfer, 2008, vol. 22, no. 2, pp. 157-167.

  7. Pirumov U.G., Roslyakov G.S. Gazovaya dinamika sopel (Gas Dynamics of nozzles), Moscow, Nauka, 1990, 368 p.

  8. Khairer E., Vanner G. Reshenie obyknovennykh differentsial’nykh uravneniy, zhestkie i differentsial’no-algebraicheskie zadachi (Solution of ordinary differential equations, stiff and differential-algebraic problems), Moscow, Mir, 1999, 685 p.

  9. Skvortsov L.M. Matematicheskoe modelirovanie, 2009, vol. 21, no. 9, pp. 54-65.

  10. Lebedev V.I., Medovikov A.A. Izvestiya vuzov. Matematika, 1998, no. 9, pp. 55-63.

  11. Naumov V. I., Krioukov V.G., Abdullin A.L., Demin A.V., Iskhakova R.L. Chemical non-equilibrium model for simulation of combustion and flow in propulsion and power generation systems. Proceedings of ASME-International Mechanical Engineering Congress and Exposition, Florida, USA, 2005, vol.1, 11 pp.

  12. Abdullin A.L., Durigon A., Kriukov V.G., Primenenie metoda splain-funktsii dlya resheniya zadach khimicheskoi kinetiki. Vestnik KGTU, 2004, vol. 3, pp. 31-34.

  13. Press W. H., Flinnery B. P. and Vetterling W. T., et al. Numerical Recipes in C: The art of scientific equation models by polynomial approximation. New Jersey: Prentice-Hall, 1988, 735p.

  14. Glarborg P., Miller J.A., Kee R.J. Kinetic Modeling and Sensitivity Analysis of Nitrogen Oxide Formation in Well-Stirred Reactors. Combustion and Flame, 1986, vol.65, pp. 177-202.

  15. Kondrat’ev V.N. Konstanty skorosti gazofaznykh reaktsii (Rate constants of gas phase reactions), Moscow, Nauka, 1974, 512 p.

Download — informational site MAI

Copyright © 2000-2022 by MAI