Article information
2019 , Volume 24, ¹ 1, p.106-119
Tyatyushkina E.S., Kozelkov A.S., Kurkin A.A., Kurulin V.V., Efremov V.R., Utkin D.A.
Evaluation of numerical diffusion of the finite volume method when modelling surface waves
The application of numerical simulation methods based on the solution of the full three-dimensional Navier-Stokes equations for modelling of wave propagation on the water surface requires the construction of a grid model containing countable nodes throughout the entire volume of water medium. Insufficient grid resolution leads to insufficient detailing of the fields of velocity and pressure, as well as volume fraction of the liquid, which increases the numerical diffusion of the method and, ultimately, leads to an underestimation of oscillation amplitudes of the medium. A large time step also results in a “blurring” of the solution and significantly reduces its stability, especially when using the schemes which compress the front of the media interface. This paper presents the results of an assessment of acceptable grid sizes and time steps expressed in dimensionless parameters with respect to the wave parameters necessary to ensure accuracy of the solution sufficient for geophysical applications. The estimate is given for the method of calculating three-dimensional multiphase flows with a free surface based on solving the Navier-Stokes equations in a one-velocity approximation based on a completely implicit connection between velocity and pressure using the finite volume method. The finite volume method for the numerical solution of the Navier-Stokes equations is implemented for use on arbitrary unstructured grids. The methodology for estimation of numerical diffusion of the calculation method is proposed. This estimation is expressed as a percentage of the wave amplitude decrease at the distance equal to the one wavelength. In turn the methodology is based on the parameters entered to estimate the acceptable grid sizes and time step for the calculation method. Based on the described methodology, the results of the estimation of the grid resolution in the horizontal and vertical directions, the estimation of the time step, and the evaluation of the influence of the discretization scheme of the convective term are presented.
[full text] [link to elibrary.ru]
Keywords: surface waves, volume of fluid, method of finite volume, difference scheme, grid resolution
doi: 10.25743/ICT.2019.24.1.008
Author(s): Tyatyushkina Elena Sergeevna Position: Junior Research Scientist Office: Russian Federal Nuclear Center All-Russian research institute of experimental physics, Nizhny Novgorod State Technical University n.a. R.E. Alekseev Address: 607189, Russia, Sarov, Mira ave., 37
Phone Office: (83130) 27372 E-mail: leno4ka-07@mail.ru 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-1107Kurkin 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 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-4399Efremov Valentin Robertovich Position: Councilor Office: Joint-Stock Company KBP n.a. Academician A. Shipunov Address: 300001, Russia, Tula, 59 Shcheglovskaya Zaseka Str.
Phone Office: (4872)25-26-27 E-mail: valentin_e@mail.ru Utkin Dmitriy Alexandrovich Office: Russian Federal Nuclear Center All-Russian research institute of experimental physics Address: 607189, Russia, Sarov, Mira ave., 37
Phone Office: (83130) 27144 E-mail: dimitryavich@yandex.ru
References: [1] Loytsyanskiy, L.G. Mekhanika zhidkosti i gaza [Mechanics of liquid and gas]. Moscow: Nauka; 1973: 847. (In Russ.)
[2] Landau, L.D., Lifshits, E.M. Gidrodinamika [Hydrodynamics]. Moscow: Nauka; 1986: 736. (In Russ.)
[3] Ferziger, J.H., Peric, M. Computational Method for Fluid Dynamics. New York: Springer-Verlag; 2002: 423.
[4] Kozelkov, A.S. Numerical technology for landslide tsunami simulations based on Navier–Stokes equations. Computational Continuum Mechanics. 2016; 9(2):218-236.
[5] Kozelkov, A.S., Kurkin, A.A., Pelinovsky, E.N., Tyatyushkina, E.S., Kurulin, V.V., Tarasova, N.V. Landslide-type tsunami modelling based on the Navier-Stokes Equations. Science of Tsunami Hazards. 2016; 35(3):106-144.
[6] Ubbink, O. Numerical prediction of two fluid systems with sharp interfaces. Thesis for the degree of PhD. London: Imperial College of Science, Technology & Medicine; 1997:138.
[7] Khrabryy, A.I. Chislennoe modelirovanie nestatsionarnykh turbulentnykh techeniy zhidkosti so svobodnoy poverkhnost'yu [Numerical simulation of unsteady turbulent fluid flow with free surface. Ph.D. Thesis]. Sankt-Peterburg; 2014: 154. (In Russ.)
[8] 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. Mathematical Models and Computer Simulations. 2014: 26(8):81-96. (In Russ.)
[9] Kozelkov, A., Kurulin, V., Emelyanov, V., Tyatyushkina, E., Volkov, K. Comparison of convective flux discretization schemes in detached-eddy simulation of turbulent flows on unstructured meshes. Journal of Scientific Computing. 2016; (67):176–191.
[10] Kozelkov, A.S., Kurulin, V.V., Lashkin, S.V., Shagaliev, R.M., Yalozo, A.V. Investigation of supercomputer capabilities for the scalable numerical simulation of computational fluid dynamics problems in industrial applications. Computational Mathematics and Mathematical Physics. 2016; 56(8):1506-1516.
[11] 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; 21(5):54-76. ( In Russ.)
[12] Fenton, J.D. A Fifth-Order Stokes Theory for Steady Waves. Coastal & Ocean Engineering. 1985; 111(2):216-234.
[13] Safronov, A.V., Deryugin, Yu.N., Zhuchkov, R.N., Zelenckii, D.K., Kozelkov, A.S., Sarazov, A.S., Kudimov, N.F., Lipnickii, Ju.M., Panasenko, A.V. Validation results for the LOGOS multifunction software package in solving problems of aerodynamics and gas dynamics for the lift-off and injection of launch vehicles. Mathematical Models and Computer Simulations. 2015; 7(2):144-153.
[14] Weinman, K. A., Valentino, M. Comparison of Hybrid RANS-LES Calculations within the Framework of Compressible and Incompressible Unstructured Solvers. Editors Shia-Hui Peng, Piotr Doerffer, Werner Haase Progress in Hybrid RANS-LES Modelling. Springer; 2010: 329-338.
[15] Kozelkov, A.S., Krutyakova, O.L., Kurulin, V.V., Lashkin, S.V., Tyatyushkina, E.S. Application of numerical schemes with singling out the boundary layer for the computation of turbulent flows using eddy-resolving approaches on unstructured grids. Computational Mathematics and Mathematical Physics. 2017; 57(6):1036-1047.
Bibliography link: Tyatyushkina E.S., Kozelkov A.S., Kurkin A.A., Kurulin V.V., Efremov V.R., Utkin D.A. Evaluation of numerical diffusion of the finite volume method when modelling surface waves // Computational technologies. 2019. V. 24. ¹ 1. P. 106-119
|