Semi-analytical implicit direct time integration method for 1-D gas dynamic problem


  • Igor Orynyak Department of Applied Mathematics at National Technical University Kiev Polytechnic Institute, Ukraine
  • Iryna Kostyushko Igor Sikorsky Kyiv Polytechnic Institute, Kyiv, Ukraine
  • Roman Mazuryk Ph D student, Department of Applied Mathematics at National Technical University Kiev Polytechnic Institute, Ukraine



stability, Lagrangian formulation, implicit method, Sod’s task, ideal gas, transfer matrix method


Sharp wave treatment for 1-D gas dynamic problem is still a chellenge for modern numerical methods. They often require too many space and time steps, produce spurious oscillation of solution, exhibit a strong numerical dissipation or divergence of results. This paper is further extension of authors’ idea of employment the analytical solution for space coordinate, where time step is a parameter which used in the space solution. Its peculiarity consists in development of additional procedure of linearization of dependence between the pressure and density. It is performed in premise that actual pressure for each space element is close to the basic pressure, attained at previous moment of time. The efficiency of method is tested on the very popular task of Sod, where two different ideal gases in a tube are separated by diaphragm, which is suddenly broken. The problem considered in Lagrangian coordinates formulation. The results obtained show the very good efficiency of method, which requires the essentially lesser time and space steps, leads to no spurious oscillation and give consistent and predictable results with respect to meshing. The accuracy of method is mostly controlled by time step, which should be larger than clearly stated theoretical lower limit. Other advantage of method is that it can calculate the process to any desired moment of time, and space meshing can be variable in time and space and can be easily adapted during the process of calculation.


  1. Eleuterio F. Toro. Riemann Solvers and Numerical Methods for Fluid Dynamics, A Practical Introduction. Third Edition Springer-Verlag. Berlin Heidelberg, 2009. DOI:
  2. J. Nеumann and R. Richtmyer, “A method for the numerical calculation of hydrodynamic shocks”, Journ. Appl. Physics, Vol. 21, No. 3, pp. 232–237, 1950. DOI:
  3. P.D. Lax, “Weak solutions of nonli near hyperbolic equations and their numerical computation”, Communic. on Pure and Appl. Math., VII, No. 1, pp. 159–193, 1954. DOI:
  4. S.K. Godunov, “Difference method for numerical calculation of discontinuous solutions of hydrodynamic equations”, Mat. Sbornik, 47, pp. 271, 1959.
  5. G.A. Sod, “Asurvey of several finite difference methods for systems of nonlinear hyperbolic conservation laws”, Journal of computational physics, 27(1), pp. 1–31, 1978. DOI:
  6. T. Kiuchi, “An implicit method for transient gas flows in pipe networks”, International Journal of Heat and Fluid Flow, 15(5), pp. 378–383, 1994. DOI:
  7. G.P. Greyvenstein, “An implicit method for the analysis of transient flows in pipe networks”, International journal for numer- ical methods in engineering, 53(5), pp. 1127–1143, 2002. DOI:
  8. L. Zhang, “Simulation of the transient flow in a natural gas compression system using a high-order upwind scheme considering the real-gas behaviors”, J. Nat. Gas Sci. Eng., 28, pp. 479–490, 2016. DOI:
  9. P.Wang et al., “Adaptive implicit finite difference method for natural gas pipeline transient flow”, Oil & Gas Sciences and Technology–Revue d’IFP Energies nouvelles, 73, 21. 2018. DOI:
  10. M. Witek and F. Uilhoorn, “Impact of hydrogen blended natural gas on linepack energy for existing high pressure pipelines”, Archives of Thermodynamics, pp. 111–124, 2022. DOI:
  11. J.F. Helgaker et al., “Validation of 1D flow model for high pressure offshore natural gas pipelines”, Journal of Natural Gas Science and Engineering, 16, pp. 44–56, 2014. DOI:
  12. V.F. Chekurin and О.М. Khymko, “Transient Processes in a Gas Flow Through a Pipeline Caused by a Local Leak”, J Math Sci, 263, pp. 166–184, 2022. DOI:
  13. M. Popescu, “Modeling of fluid dynamics interacting with ductile fraction propagation in high pressure pipeline”, Acta Me- chanica Sinica, 25(3), pp. 311–318, 2009. DOI:
  14. Liu Fayue, Chen Huandong and Chen Jiayue, “Numerical simulation of shock wave problems with the two-phase two-fluid model”, Progress in Nuclear Energy, 121, 2020.DOI:
  15. T. Feng et al., “Numerical research on water hammer phenomenon of parallel pump-valve system by coupling FLUENT with RELAP5”, Ann. Nucl. Energy, 109, pp. 318–326, 2017. DOI:
  16. S.J. Yi, M.S. Chung and I. Toth, “Pressure wave propagation phenomena in the pipe system of a nuclear power plant”, Ann. Nucl. Energy, 87, pp. 157–161, 2016. DOI:
  17. R. Mahmoodi, A. Zolfaghari and A. Minuchehr, “Simulation of pressure waves propagation following LOCA in piping systems using Laplace Transform Finite Volume”, Ann. Nucl. Energy, 124, pp. 164–171, 2019. DOI:
  18. Marc O. Delchini, Jean C. Ragusa and Ray A. Berry, “Simulations of single- and two-phase shock tubes and gravity-driven wave problems with the RELAP-7 nuclear reactor system analysis code”, Nuclear Engineering and Design, Vol. 319, pp. 106–116, 2017. DOI:
  19. F. Daude and P. Galon, “A Finite-Volume approach for compressible single-and two-phase flows in flexible pipelines with fluid-structure interaction”, Journal of Computational Physics, Vol. 362, pp. 375–408, 2018. DOI:
  20. F. Daude and P. Galon, “Simulations of single- and two-phase shock tubes across abrupt changes of area and branched junc- tions”, Nuclear Engineering and Design, Vol. 365, 2020. DOI:
  21. I. Orynyak, R. Mazuryk and V. Tsybulskyi, “Semi-analytical implicit direct time integration scheme on example of 1-D wave propagation problem”, Mechanics and Advanced Technologies, Vol. 6, No. 2, pp. 115–123, 2022. DOI:
  22. Peng Wang, Bo Yu, Yajun Deng and Yu Zhao, “Comparison study on the accuracy and efficiency of the four forms of hydraulic equation of a natural gas pipeline based on linearized solution”, Journal of Natural Gas Science and Engineering, Vol. 22, pp. 235–244, 2015. DOI:
  23. E.B. Wylie, V.L. Streeter, and M.A. Stoner, “Unsteady-state natural-gas calculations in complex pipe systems”, Soc. Petrol. Eng. J., Vol. 14 (01), pp. 35–43, 1974. DOI:
  24. Barley, J., “Thermal decoupling: an investigation”,in PSIG Annual Meeting, Paper Number: PSIG-1213, May. 2012. pp. 15–18.
  25. J.F. Helgaker and T. Ytrehus, “Coupling between continuity/momentum and energy equation in 1D gas flow”, Energy Procedia, Vol. 26, pp. 82–89, 2012. DOI:
  26. F. Leckie and E. Pestel, “Transfer-matrix fundamentals”, International Journal of Mechanical Sciences, Vol. 2, Issue 3, pp. 137–167, 1960. DOI:
  27. I. Orynyak, R. Mazuryk and A. Orynyak, “Basic (discontinuous) and smoothing up (conjugated) solutions in transfer matrix method for static geometrically nonlinear beam and cable in plane”, Journal of Engineering Mechanics, Vol. 46, 5, 2020. DOI:




How to Cite

I. Orynyak, I. Kostyushko, and R. Mazuryk, “Semi-analytical implicit direct time integration method for 1-D gas dynamic problem”, Mech. Adv. Technol., vol. 7, no. 1 (97), pp. 91–99, Apr. 2023.