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, 738 p.
  2. J. Nеumann , R. Richtmyer , A method for the numerical calculation of hydrodynamic shocks, Journ. Appl. Physics, 21, № 3 (1950), 232—237
  3. P. D. Lax, Weak solutions of nonlinear hyperbolic equations and their numerical computation. Communic. on Pure and Appl. Math., VII, N 1 (1954), 159—193.
  4. С. К. Годунов, Разностный метод численного расчета разрывных решений уравнений гидродинамики, Матем. сб., 1959, том 47(89), номер 3, 271–306. S. K. GODUNOV, Difference method for numerical calculation of discontinuous solutions of hydrodynamic equations Mat. Sbornik 47 (1959), 271.
  5. Sod, G. A. (1978). A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws. Journal of computational physics, 27(1), 1-31,
  6. Kiuchi, T. (1994). An implicit method for transient gas flows in pipe networks. International Journal of Heat and Fluid Flow, 15(5), 378-383.
  7. Greyvenstein, G. P. (2002). An implicit method for the analysis of transient flows in pipe networks. International journal for numerical methods in engineering, 53(5), 1127-1143.
  8. Zhang L. (2016) 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, 479–490.
  9. Wang, P., Yu, B., Han, D., Li, J., Sun, D., Xiang, Y., & Wang, L. (2018). Adaptive implicit finite difference method for natural gas pipeline transient flow. Oil & Gas Sciences and Technology–Revue d’IFP Energies nouvelles, 73, 21.
  10. Witek, M., & Uilhoorn, F. (2022). Impact of hydrogen blended natural gas on linepack energy for existing high pressure pipelines. Archives of Thermodynamics, 111-124.
  11. Helgaker, J. F., Oosterkamp, A., Langelandsvik, L. I., & Ytrehus, T. (2014). Validation of 1D flow model for high pressure offshore natural gas pipelines. Journal of Natural Gas Science and Engineering, 16, 44-56.
  12. Chekurin, V.F., Khymko, О.М. Transient Processes in a Gas Flow Through a Pipeline Caused by a Local Leak. J Math Sci 263, 166–184 (2022).
  13. Popescu, M. (2009). Modeling of fluid dynamics interacting with ductile fraction propagation in high pressure pipeline. Acta Mechanica Sinica, 25(3), 311-318.
  14. Fayue Liu, Huandong Chen, Jiayue Chen. Numerical simulation of shock wave problems with the two-phase two-fluid model. Progress in Nuclear Energy, 121 (2020), 103259
  15. Feng, T., Zhang, D., Song, P., et al., 2017. Numerical research on water hammer phenomenon of parallel pump-valve system by coupling FLUENT with RELAP5. Ann. Nucl. Energy 109, 318–326.
  16. Yi, S.J., Chung, M.S., Toth, I., 2016. Pressure wave propagation phenomena in the pipe system of a nuclear power plant. Ann. Nucl. Energy 87, 157–161.
  17. Mahmoodi, R., Zolfaghari, A., Minuchehr, A., 2019. Simulation of pressure waves propagation following LOCA in piping systems using Laplace Transform Finite Volume. Ann. Nucl. Energy 124, 164–171.
  18. Marc O. Delchini, Jean C. Ragusa, 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 319 (2017) 106–116.
  19. Daude, F., & Galon, P. (2018). A Finite-Volume approach for compressible single-and two-phase flows in flexible pipelines with fluid-structure interaction. Journal of Computational Physics, 362, 375-408.
  20. F. Daude, P. Galon. Simulations of single- and two-phase shock tubes across abrupt changes of area and branched junctions. Nuclear Engineering and Design 365 (2020) 110734.
  21. Orynyak, I., Mazuryk, R., & Tsybulskyi, V. (2022). Semi-analytical implicit direct time integration scheme on example of 1-D wave propagation problem. Mechanics and Advanced Technologies, 6(2), 115–123.
  22. Peng Wang, Bo Yu , Yajun Deng, 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, 22, (2015), 235-244
  23. Wylie, E.B., Streeter, V.L., Stoner, M.A., 1974. Unsteady-state natural-gas calculations in complex pipe systems. Soc. Petrol. Eng. J. 14 (01), 35-43.
  24. Barley, J., May. 2012. Thermal decoupling: an investigation. In: PSIG Annual Meeting, Paper Number: PSIG-1213, pp. 15-18
  25. Helgaker, J.F., Ytrehus, T., 2012. Coupling between continuity/momentum and energy equation in 1D gas flow. Energy Procedia 26, 82-89.
  26. F. Leckie, E. Pestel. Transfer-matrix fundamentals. International Journal of Mechanical Sciences, Volume 2, Issue 3, 1960, Pages 137-167.
  27. Orynyak I., Mazuryk R., Orynyak A. (2020). “Basic (discontinuous) and smoothing up (conjugated) solutions in transfer matrix method for static geometrically nonlinear beam and cable in plane”. Journal of Engineering Mechanics, V46,




How to Cite

Orynyak, I., Kostyushko, I., & Mazuryk, R. (2023). Semi-analytical implicit direct time integration method for 1-D gas dynamic problem. Mechanics and Advanced Technologies, 7(1).