Shallow water equations are widely used in the simulation of those geophysical flows for which the flow horizontal length scale is much greater than the vertical one. Inspired by the example of lava flows, we consider here a modified model with an additional transport equation for a scalar quantity (e.g., temperature), and the derivation of the shallow water equations from depth-averaging the Navier-Stokes equations is presented. The assumption of constant vertical profiles for some of the model variables is relaxed allowing the pres- ence of vertical profiles, and it follows that the non-linearity of the flux terms results in the introduction of appropriate shape coefficients. The space discretization of the resulting system of hyperbolic partial differential equations is obtained with a modified version of the finite volume central-upwind scheme intro- duced by Kurganov and Petrova in 2007. The time discretization is based on an implicit- explicit Runge-Kutta method which couples properly the hyperbolic part and the stiffsource terms, avoiding the use of a very small time step; the use of complex arithmetic increases accuracy in the implicit treatment of stiffterms. The whole scheme is proved to preserve the positivity of flow thickness and the stationary steady-states. Some numerical experiments are performed to validate the proposed method and to show the incidence on the numerical solutions of shape coefficients introduced in the model.

Modified shallow water model for viscous fluids and positivity preserving numerical approximation

E. Biagioli;F. Di Benedetto
2021

Abstract

Shallow water equations are widely used in the simulation of those geophysical flows for which the flow horizontal length scale is much greater than the vertical one. Inspired by the example of lava flows, we consider here a modified model with an additional transport equation for a scalar quantity (e.g., temperature), and the derivation of the shallow water equations from depth-averaging the Navier-Stokes equations is presented. The assumption of constant vertical profiles for some of the model variables is relaxed allowing the pres- ence of vertical profiles, and it follows that the non-linearity of the flux terms results in the introduction of appropriate shape coefficients. The space discretization of the resulting system of hyperbolic partial differential equations is obtained with a modified version of the finite volume central-upwind scheme intro- duced by Kurganov and Petrova in 2007. The time discretization is based on an implicit- explicit Runge-Kutta method which couples properly the hyperbolic part and the stiffsource terms, avoiding the use of a very small time step; the use of complex arithmetic increases accuracy in the implicit treatment of stiffterms. The whole scheme is proved to preserve the positivity of flow thickness and the stationary steady-states. Some numerical experiments are performed to validate the proposed method and to show the incidence on the numerical solutions of shape coefficients introduced in the model.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11567/1065662
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 0
  • ???jsp.display-item.citation.isi??? 0
social impact