New Numerical Algorithm for the Multi-Layer Shallow Water Equations Based on the Hyperbolic Decomposition and the CABARET Scheme

V. M. Goloviznin✉, Pavel A. Maiorov, Petr A. Maiorov, A. V. Solovjov

Nuclear Safety Institute, Russian Academy of Sciences, Moscow, Russian Federation

✉ e-mail: gol@ibrae.ac.ru

Abstract

Purpose. The present article is devoted to describing a new method of numerical solution for hydrostatic approximation of incompressible hydrodynamic problems with free surfaces and variable density.

Methods and Results. The algorithm is based on the hyperbolic decomposition method, i. e. representation of a multilayer model as a sum of the one-layer models interacting by means of the reaction forces through the layers’ interfaces. The forces acting on the upper and lower interfaces of each layer are interpreted as the external ones which do not break hyperbolicity of the equations system for each layer. The explicit CABARET scheme is used to solve a system of hyperbolic equations with variable density in each layer. The scheme is of the second approximation order and the time reversibility. Its feature consists in the increased number of freedom degrees: along with the conservative-type variables referred to the centers of the calculated cells, applied are the flux-type variables related to the middle of the vertical edges of these cells. The system of the multilayer shallow water equations is not unconditionally hyperbolic, and in case hyperbolicity is lost, it becomes ill-posed. Hyperbolic decomposition does not remove incorrectness of the original system of the multilayer shallow water equations. To regularize the numerical solution, the following set of tools is propose: filtration of the flow variables at each time step; super-implicit approximation of the pressure gradient; linear artificial viscosity and transition to the Euler-Lagrangian (SEL) variables that leads to the mass and momentum exchange between the layers. Such transition to the SEL variables is the basic tool for stabilizing numerical solution at large times. The rest of the tricks are the auxiliary ones and used for fine tuning.

Conclusions. It is shown that regularizing and guaranteeing the problems’ stability requires not only reconstruction of the computational grid at each time step, but also application of the flow-type variables’ filtering and the artificial viscosity simulating turbulent mixing.

Keywords

non-stationary hydrodynamics, free surface, hydrostatic approximation, ill-posed problems, variable density, numerical algorithm, hyperbolic decomposition, CABARET scheme

Acknowledgements

The work was supported by the Russian Science Foundation, project No 18-11-00163. The authors are grateful to Zalesny V.B. and Semenov E.V. for productive discussions and constructive comments.

Original russian text

Original Russian Text © V.M. Goloviznin, Pavel A. Maiorov, Petr A. Maiorov, A.V. Solovjov, 2019, published in MORSKOY GIDROFIZICHESKIY ZHURNAL, Vol. 35, Iss. 6, pp. 600–620 (2019)

For citation

Goloviznin, V.M., Maiorov, Pavel A., Maiorov, Petr A. and Solovjov A.V., 2019. New Numerical Algorithm for the Multi-Layer Shallow Water Equations Based on the Hyperbolic Decomposition and the CABARET Scheme. Physical Oceanography, 26(6), pp. 528-546. doi:10.22449/1573-160X-2019-6-528-546

DOI

10.22449/1573-160X-2019-6-528-546

References

  1. Marchuk, G.I., Dymnikov, V.P. and Zalesny, V.B., 1987. [Mathematical Models in Geophysical Hydrodynamics and Methods of their Numerical Implementation]. Leningrad: Gidrometeoizdat, 296 p. (in Russian).
  2. Zalesny, V.B., Marchuk, G.I., Agoshkov, V.I., Bagno, A.V., Gusev, A.V., Diansky, N.A., Moshonkin, S.N., Tamsalu, R. and Volodin, E.M., 2010. Numerical Simulation of Large- Scale Ocean Circulation Based on the Multicomponent Splitting Method. Russian Journal of Numerical Analysis and Mathematical Modelling, [e-journal] 25(6), pp. 581-609. https://doi.org/10.1515/RJNAMM.2010.036
  3. Marchuk, G.I., Paton, B.E., Korotaev, G.K. and Zalesny, V.B., 2013. Data-Computing Technologies: A New Stage in the Development of Operational Oceanography. Izvestiya, Atmospheric and Oceanic Physics, [e-journal] 49(6), pp. 579-591. doi:10.1134/S000143381306011X
  4. Dotsenko, S.F., Zalesny, V.B. and Sannikova, N.K.V., 2016. Block Approach to the Simulation of Circulation and Tides in the Black Sea. Physical Oceanography, [e-journal] (1), pp. 3-19. doi:10.22449/1573-160X-2016-1-3-19
  5. Zalesny, V.B., Gusev, A.V. and Fomin, V.V., 2016. Numerical Model of Nonhydrostatic Ocean Dynamics Based on Methods of Artificial Compressibility and Multicomponent Splitting. Oceanology, [e-journal] 56(6), pp. 876-887. doi:10.1134/S0001437016050167
  6. Goloviznin, V.M., Zaytsev, M.A., Karabasov, S.A. and Korotkin, I.A., 2013. [New Algorythms of Numerical Hydrodynamic for Multiprocessor Computational System]. Moscow: MSU, 467 p. (in Russian).
  7. Audusse, E., 2005. A Multilayer Saint-Venant Model: Derivation and Numerical Validation. Discrete & Continuous Dynamical Systems - B, [e-journal] 5(2), pp. 189-214. https://doi.org/10.3934/dcdsb.2005.5.189
  8. Audusse, E. and Bristeau, M.-O., 2007. Finite-Volume Solvers for a Multilayer Saint-Venant System. International Journal of Applied Mathematics and Computer Science, [e-journal] 17(3), pp. 311-320. https://doi.org/10.2478/v10006-007-0025-0
  9. Ovsyannikov, L.V., 1979. Two-Layer “Shallow Water” Model. Journal of Applied Mechanics and Technical Physics, [e-journal] 20(2), pp. 127-135. https://doi.org/10.1007/BF00910010
  10. Duchêne, V., 2013. A Note on the Well-Posedness of the One-Dimensional Multilayer Shallow Water Model. Available at: https://hal.archives-ouvertes.fr/hal-00922045/document [Accessed: 12 July 2019].
  11. Chandrasekhar, S., 1961. Hydrodynamic and Hydromagnetic Stability. London: Oxford University Press, 652 p.
  12. Chukharev, A.M., Runovsky, K.V. and Kulsha, O.E., 2017. Modeling of Turbulent Patches Statistical Distribution in the Stratified Ocean Layers. Physical Oceanography, [e-journal] (5), pp. 31-41. doi:10.22449/1573-160X-2017-5-31-41
  13. Kollman, V., ed., 1980. Prediction Methods for Turbulent Flows. Washington, D.C.: Hemisphere Pub., 468 p.
  14. Audusse, E., Bristeau, M.-O., Perthame, B. and Sainte-Marie, J., 2011. A Multilayer Saint- Venant System with Mass Exchanges for Shallow Water Flows. Derivation and Numerical Validation. ESAIM: M2AN, [e-journal] 45(1), pp. 169-200. doi:10.1051/m2an/2010036
  15. Audusse, E., Benkhaldoun, F., Sari, S., Seaid, M. and Tassi, P., 2014. A Fast Finite Volume Solver for Multi-Layered Shallow Water Flows with Mass Exchange. Journal of Computational Physics, [e-journal] 272, pp. 23-45. https://doi.org/10.1016/j.jcp.2014.04.026
  16. Hirt, C.W., Amsden, A.A. and Cook, J.L., 1974. An Arbitrary Lagrangian Eulerian Computing Method for All Flow Speeds. Journal of Computational Physics, [e-journal] 14(3), pp. 227-253. https://doi.org/10.1016/0021-9991(74)90051-5
  17. Ringler, T.D. and Randall, D.A., 2002. The ZM Grid: An Alternative to the Z Grid. Monthly Weather Review, [e-journal] 130(5), pp. 1411-1422. doi:10.1175/1520- 0493(2002)130<1411:TZGAAT<2.0.CO;2
  18. Audusse, E., Bristeau, M.-O., Pelanti, M. and Sainte-Marie, J., 2011. Approximation of the Hydrostatic Navier-Stokes System for Density Stratified Flows by a Multilayer Model: Kinetic Interpretation and Numerical Solution. Journal of Computational Physics, [e-journal] 230(9), pp. 3453-3478. https://doi.org/10.1016/j.jcp.2011.01.042
  19. Stewart, A.L. and Dellar, P.J., 2010. Multilayer Shallow Water Equations with Complete Coriolis Force. Part 1. Derivation on a Non-Traditional Beta-Plane. Journal of Fluid Mechanics, [e-journal] 651, pp. 387-413. https://doi.org/10.1017/S0022112009993922
  20. Bermudez, A. and Vazquez, M.E., 1994. Upwind Methods for Hyperbolic Conservation Laws with Source Terms. Computers & Fluids, [e-journal] 23(8), pp. 1049-1071. https://doi.org/10.1016/0045-7930(94)90004-3
  21. Toro, E.F., 2009. Riemann Solvers and Numerical Methods for Fluid Dynamics. Berlin: Springer, 724 p. doi:10.1007/b79761
  22. Karabasov, S.A. and Goloviznin, V.M., 2009. Compact Accurately Boundary-Adjusting High-Resolution Technique for Fluid Dynamics. Journal of Computational Physics, [e-journal] 228(19), pp. 7426-7451. https://doi.org/10.1016/j.jcp.2009.06.037
  23. Goloviznin, V.M., Kanyukova, V.D. and Samarskaya, E.A., 1983. Super-Implicit Difference Schemes of Gas Dynamics. Differential Equations, 19(7), pp. 1186-1197.
  24. Goloviznin, V.M. and Isakov, V.A., 2017. Balance-Characteristic Scheme as Applied to the Shallow Water Equations over a Rough Bottom. Computational Mathematics and Mathematical Physics, [e-journal] 57(7), pp. 1140-1157. https://doi.org/10.1134/S0965542517070089
  25. Richtmyer, R.D. and Morton, K.W., 1967. Difference Methods for Initial-Value Problems. New York: Interscience Publishers, 405 p.
  26. Goloviznin, V.M., 1982. On a Method of Introducing Artificial Dissipation into the Variational-Difference Schemes of Magnetohydrodynamics. USSR Computational Mathematics and Mathematical Physics, [e-journal] 22(1), pp. 150-156. https://doi.org/10.1016/0041-5553(82)90172-0
  27. Goloviznin, V.M., Solovjov, A.V. and Zalesny, V.B., 2018. A New Algorithm for Solving the Shallow Water Equations on the Sphere Based on the Cabaret Scheme. Journal of Physics: Conference Series, [e-journal] 1128(1), 012091. https://doi.org/10.1088/1742-6596/1128/1/012091
  28. Evtushenko, Y.G., Gorchakov, A.Y. and Goloviznin, V.M., 2018. Fast Automatic Differentiation in Problems Variations Four-Dimensional Data Assimilation (4Dvar). Journal of Physics: Conference Series, [e-journal] 1128(1), 012001. https://doi.org/10.1088/1742-6596/1128/1/012001
  29. Vedernikov, A.B. and Kholodov, A.S., 1990. Numerical Investigations of Flows of the Double-and Threelayer Fluid in Framework of “Shallow Water” Model. Matematicheskoe Modelirovanie, 2(6), pp. 9-18 (in Russian).
  30. Kurganov, A. and Petrova, G., 2009. Central-Upwind Schemes for Two-Layer Shallow Water Equations. SIAM Journal on Scientific Computing, [e-journal] 31(3), pp. 1742-1773. https://doi.org/10.1137/080719091
  31. Elizarova, T.G. and Ivanov, A.V., 2016. Quasi-Gasdynamic Algorithm for Numerical Solution of Two-Layer Shallow Water Equations. KIAM Preprint No. 69. 27 p. (in Russian). doi:10.20948/prepr-2016-69
  32. Shankar, N.J., Cheong, H.F. and Sankaranarayanan, S., 1997. Multilevel Finite-Difference Model for Three-Dimensional Hydrodynamic Circulation. Ocean Engineering, [e-journal] 24(9), pp. 785-816. https://doi.org/10.1016/S0029-8018(96)00036-4
  33. Bouchut, F. and Zeitlin, V., 2010. A Robust Well-Balanced Scheme for Multi-Layer Shallow Water Equations. Discrete & Continuous Dynamical Systems - B, [e-journal] 3(4), pp. 739- 758. http://doi.org/10.3934/dcdsb.2010.13.739
  34. Couderc, F., Duran, A. and Vila, J.-P., 2017. An Explicit Asymptotic Preserving low Froude Scheme for the Multilayer Shallow Water Model with Density Stratification. Journal of Computational Physics, [e-journal] 343, pp. 235-270. https://doi.org/10.1016/j.jcp.2017.04.018

Download the article (PDF)

We use Yandex Metrics. Read more.

Мы используем Яндекс Метрику. Подробнее.