|
Article information
2016 , Volume 21, ¹ 5, p.54-76
Kozelkov A.S., Meleshkina D.P., Kurkin A.A., Tarasova N.V., Lashkin S.V., Kurulin V.V.
Fully implicit method for solution of Navier - Stokes equations for simulation of multiphase flows with free surface
This paper presents a method for calculation of multiphase flows with a free surface. According to the current experience for calculations of a variety of different flows, one of the most universal methods is the SIMPLE algorithm, based on the modified semiimplicit scheme that uses splitting into physical processes . Along with the ability to generalize this method to the case of strongly compressible flows, its major advantage is an efficient way of implicit approximation of the equations, allowing arbitrarily large calculation time for its application. Such a modification is a combined algorithm, which increases the rate of convergence of the total calculation algorithm compared to the classical SIMPLE method due to an implicit relation between the speed and the pressure. Such implicit term in the equation of conservation of momentum is the pressure gradient, whereas in the continuity equation it is attributed to the mass flow. This paper presents a numerical scheme for modelling of multiphase flows with a free surface. The procedure to the solution relies on the classical approach - the SIMPLE method with splitting into physical processes. With the help of some modification it leads to the combined method based on a fully implicit relation between the pressure and velocity fields. For the simulation of multiphase flows, the considered approach employs one-speed model to describe the mix. This model allows simulating any number of phases, including any combination of liquid, solid or gaseous media. Only one system of equations for all phases is required. For the calculation of free surface flows the numerical scheme is closed by a transport equation for the volume fractions of the phases, which simulates the motion of a free surface using the special M-CICSAM scheme. The correct implementation of the model is tested on a series of numerical experiments for problems of the collapse of a dam, the fall of a box into the water and a hydraulic shock.
[full text]
Keywords: multiphase flows, free surface, numerical simulation, equations Navier-Stokes
Author(s):
Kozelkov Andrey Sergeevich
Dr.
Position: Professor
Office: Russian Federal Nuclear Center, All-Russian Research Institute OF Experimental Physics
Address: 607189, Russia, Sarov, Mira ave., 37
Phone Office: (83130) 2 73 43
E-mail: askozelkov@mail.ru
SPIN-code: 6563-1107Meleshkina Dariya Pavlovna
Position: Junior Research Scientist
Office: Russian Federal Nuclear Center, All-Russian Research Institute OF Experimental Physics
Address: 607189, Russia, Sarov, Mira ave., 37
Phone Office: (83130) 28017
E-mail: dasha2008k@mail.ru
Kurkin Andrey Aleksandrovich
Dr. , Professor
Position: General Scientist
Office: Nizhny Novgorod State Technical University n.a. R.E. Alekseev
Address: 603950, Russia, Nizhny Novgorod, 24 Minin str.
Phone Office: (831) 436 04 89
E-mail: aakurkin@gmail.com
Tarasova Nataliya Vladimirovna
Position: Senior Research Scientist
Office: Russian Federal Nuclear Center, All-Russian Research Institute OF Experimental Physics
Address: 607189, Russia, Sarov, Mira ave., 37
Phone Office: (83130) 2 75 50
E-mail: tara271978@mail.ru
Lashkin Sergey Victorovich
Position: Head of Laboratory
Office: Russian Federal Nuclear Center, All-Russian Research Institute OF Experimental Physics
Address: 607189, Russia, Sarov, Mira ave., 37
Phone Office: (83130) 27506
E-mail: lashkinsv@gmail.com
Kurulin Vadim Viktorovich
PhD.
Position: Head of department
Office: Russian Federal Nuclear Center All-Russian research institute of experimental physics
Address: 607189, Russia, Sarov, Mira ave., 37
Phone Office: (83130) 27567
E-mail: kurulin@mail.ru
SPIN-code: 1984-4399 References:
[1] Leonov, À.À., Chudanov, V.V., Aksenova, À.Å. Metody pryamogo chislennogo modelirovaniya v dvukhfaznykh sredakh [Methods of direct numerical simulation in two-phase media]. Moscow: Nauka; 2013: 197. (In Russ.)
[2] Volkov, K.N., Emelyanov, V.N. Techeniya gaza s chastitsami [Flows of gas with particles]. Moscow: Fismatlit; 2008: 600. (In Russ.)
[3] Kolev, N.I. Multiphase flow dynamics. Berlin, Heidelberg: Springer-Verlag, 2007: 751.
[4] Strelets, M.Kh., Shur, M.L. Methods of scaling the compressibility in calculating stationary flows of a viscous gas at arbitrary mach numbers. USSR Computational Mathematics and Mathematical Physics. 1988; 28(1):165–173.
[5] Tarasova, N.V., Kozelkov, A.S., Meleshkina, D.P., Lashkin, S.V., Denisova, O.V., Sizova, M.A. Specifics of the simple algorithm application for the calculation of incompressible flows. Voprosy Atomnoy Nauki i Tekhniki. Series: Mathematical Modelling of Physical Processes. 2015; (3):20–34. (In Russ.)
[6] Turkel, E. Preconditioned methods for solving the incompressible and low speed compressible Equations. Journal of Computational Physics. 1987; 72(2):277–298.
[7] Najm, H.N., Knio, O.M. Modelling low Mach number reacting flow with detailed chemistry and transport. Journal of Scientific Computing. 2005; 25(1-2):263–287.
[8] Ferziger, J.H., Peric, M. Computational methods fluid dynamics. Berlin: Springer; 2001: 423.
[9] Fletcher, C.A.J. Computational techniques for fluid dynamics. Berlin: Springer-Verlag; 1988: 893.
[10] Chen, Z.J., Przekwas, A.J. A coupled pressure-based computational method for incompressible/compressible flows. Journal of Computational Physics. 2010; (229):9150–9165.
[11] Darwish, M., Sraj, I., Moukalled, F. A coupled finite volume solver for the solution of incompressible flows on unstructured grids. Journal of Computational Physics. 2009; (228):180–201.
[12] Moukalled, F., Darwish, M. Pressure-based algorithms for multifluid flow at all speeds. Part I: Mass conservation formulation. Numerical heat transfer. Part B: Fundamentals. 2004; 45(6):495–522.
[13] Kozelkov, A.S., Deryugin, Y.N., Lashkin, S.V., Silaev, D.P., Simonov, P.G., Tyatyushkina, E.S. Implementation in logos software of a computational scheme for a viscous incompressible fluid using the multigrid method based on algorithm simple. Voprosy Atomnoy Nauki i Tekhniki. Series: Mathematical Modelling of Physical Processes. 2013; (4):44–56. (In Russ.)
[14] Kozelkov, A.V., Volkov, K.N., Deryugin, Yu.N., Emelyanov, V.N., Karpenko, A.G. Metody uskoreniya gazodinamicheskikh raschetov na nestrukturirovannykh setkakh [The acceleration of gas dynamic calculations on unstructured grids]. Moscow: Fizmatlit; 2013: 536. (In Russ.)
[15] Rhie, C.M., Chow, W.L. A numerical study of the turbulent flow past an isolated airfoil with trailing edge separation. American Institute of Aeronautics and Astronautics. 1983; (21):1525–1532.
[16] Jasak, H. Error Analysis and Estimation for the finite volume method with applications to fluid flows. PhD, Department of Mechanical Engineering, Imperial College of Science. London; 1996: 394.
[17] Kozelkov, A.S., Shagaliev, R.M., Dmitriev, S.M., Kurkin, A.A., Volkov, K.N., Deryugin, Yu.N., Emelyanov, V.N., Pelinovsky, E.N., Legchanov, M.A. Matematicheskie modeli i algoritmy dlya chislennogo modelirovaniya zadach gidrodinamiki I aerodinamiki: Uchebnoe posobie [ Mathematical models and algorithms for numerical modeling of hydrodynamics and aerodynamics: A tutorial]. Nizhniy Novgorod: NGTU; 2014: 166. (In Russ.)
[18] Kozelkov, A.S., Kurulin, V.V., Tyatyushkina, E.S., Puchkova, O.L. Application of the detached eddy simulation model for viscous incompressible turbulent flow simulations on unstructured grids. Matematicheskoe modelirovanie. 2014; 26(8):81–96. (In Russ.)
[19] Kozelkov, A.S., Kurulin, V.V. Eddy-resolving numerical scheme for simulation of turbulent incompressible flows. USSR Computational Mathematics and Mathematical Physics. 2015; 55(7):1232-1241.
[20] Roache, P.J. Computational Fluid Dynamics. Hermosa: Albuquerque; 1982: 446.
[21] Waclawczyk, T., Koronowicz, T. Remarks on prediction of wave drag using VOF method with interface capturing àpproach. Archives of Civil and Mechanical Engineering. 2008; (8):5–14.
[22] Bakhvalov, N.S. On the convergence of a relaxation method with natural constraints on the elliptic operator. USSR Computational Mathematics and Mathematical Physics.1966; 6(5):101–135
[23] Fedorenko, R.P. A relaxation method for solving elliptic difference equations. USSR Computational Mathematics and Mathematical Physics. 1961; 1(5):1092-1096.
[24] Petrov, I.B., Lobanov, A.I. Lektsii po vychisl. matematike: Uchebnoe posobie [Lectures on computational mathematics: Tutorial]. Moscow: BINOM; 2006: 523.(In Russ.)
[25] Saad, Y. Iterative methods for sparse linear systems. Minneapolis: SIAM; 2003: 528.
[26] Brandt, A. Guide to multigrid development. In book: Eds. Dold A., Eckmann B. Lecture notes in mathematics. Berlin: Springer-Verlag; 1982: 651.
[27] Stuben, K., Trottenberg, U. Multigrid methods: fundamental algorithms, model problem analysis and applications. Berlin: Springer; 1982: 176.
[28] Vanek, P., Mandel, J., Brezina, M. Algebraic multigrid based on smoothed aggregation for second and fourth order problems. Computing. 1996; 56(3):179–196.
[29] Volkov, K.N., Deriugin, Yu.N., Emelyanov, V.N., Kozelkov, A.S., Teterina, I.V. An algebraic multigrid method in problems of computational. Numerical Methods and Programming. 2014; 15(2):183–200. (In Russ.)
[30] Kozelkov, A.S., Lashkin, S.V., Yalozo, A.V. The research of the supercomputers potential for scalable numerical simulation of problems of hydrodynamics in industrial applications. Proceedings of the XI All-russian Congress on Fundamental Problems of the Theoretical and Applied Mechanics, Kazan; 2015: 1853–1855. (In Russ.)
[31] Kozelkov, A.S., Kurulin, V.V., Puchkova, O.L., Lashkin, S.V. Simulation of turbulent flows using an algebraic Reynolds stress model with universal wall functions. Computational Continuum Mechanics. 2014; 7(1):40–51. (In Russ.)
[32] Kozelkov, A.S., Kurkin, A. A., Pelinovskii, E.N. Effect of the angle of water entry of a body on the generated wave heights. Fluid Dynamics. 2016; 51(2):288–289.
[33] Kozelkov, A.S., Kurkin, A.A., Pelinovskii, E.N., Kurulin, V.V., Tyatyushkina, E.S. Modelling the disturbances in the lake Chebarkul caused by the fall of the meteorite in 2013. Fluid Dynamics. 2015; 50(6):828–840.
[34] Kozelkov, A.S. Methods of numerical modelling of landslide type tsunami on the basis of the Navier-Stokes equations. Computational continuum mechanics. 2016; 9(2):218-236. (In Russ.)
[35] Ubbink, O. Numerical prediction of two fluid systems with sharp interfaces. London: Imperial College of Science, Technology and Medicine; 1997: 193.
[36] Bukreev, V.I., Gusev, A.V. Gravity waves generated by a body falling onto shallow water. Journal of Applied Mechanics and Technical Physics. 1996; 37(2):224-231.
[37] Hirt, C.W., Nichols, B.D. Volume of fluid (VOF) method for the dynamics of free boundaries. Journal of Computational Physics. 1981; (39):201–225.
Bibliography link:
Kozelkov A.S., Meleshkina D.P., Kurkin A.A., Tarasova N.V., Lashkin S.V., Kurulin V.V. Fully implicit method for solution of Navier - Stokes equations for simulation of multiphase flows with free surface // Computational technologies. 2016. V. 21. ¹ 5. P. 54-76
|
