High-fidelity simulations of gravity currents using a high-order finite-difference spectral vanishing viscosity approach
File(s)2021_CF_LAIZET.pdf (4.39 MB)
Accepted version
Author(s)
Frantz, Ricardo
Deskos, Georgios
Laizet, Sylvain
Silvestrini, Jorge
Type
Journal Article
Abstract
This numerical work investigates the potential of a high-order finite-difference spectral
vanishing viscosity approach to simulate gravity currents at high Reynolds numbers.
The method introduces targeted numerical dissipation at small scales through altering
the discretisation of the second derivatives of the viscous terms in the incompressible
Navier-Stokes equations to mimic the spectral vanishing viscosity (SVV) operator, originally designed for the regularisation of spectral element method (SEM) solutions of pure
advection problems. Using a sixth-order accurate finite-difference scheme, the adoption of
the SVV method is straightforward and comes with a negligible additional computational
cost. In order to assess the ability of this high-order finite-difference spectral vanishing
viscosity approach, we performed large-eddy simulations (LES) of a gravity current in
a channelised lock-exchange set-up with our SVV model and with the well-known explicit static and dynamic Smagorinsky sub-grid scale (SGS) models. The obtained data
are compared with a direct numerical simulation (DNS) based on more than 800 million
mesh nodes, and with experimental measurements. A framework for the energy budget
is introduced to investigate the behaviour of the gravity current. First, it is found that
the DNS is in good agreement with the experimental data for the evolution of the front
location and velocity field as well as for the stirring and mixing inside the gravity current. Secondly, the LES performed with less than 0.4% of the total number of mesh nodes
compared to the DNS, can reproduce the main features of the gravity currents, with the
SVV model yielding slightly more accurate results. It is also found that the dynamic
Smagorinsky model performs better than its static version. For the present study, the
static and dynamic Smagorinsky models are 1.8 and 2.5 times more expensive than the
SVV model, because the latter does not require the calculation of explicit SGS terms in
the Navier-Stokes equations nor spatial filtering operations.
vanishing viscosity approach to simulate gravity currents at high Reynolds numbers.
The method introduces targeted numerical dissipation at small scales through altering
the discretisation of the second derivatives of the viscous terms in the incompressible
Navier-Stokes equations to mimic the spectral vanishing viscosity (SVV) operator, originally designed for the regularisation of spectral element method (SEM) solutions of pure
advection problems. Using a sixth-order accurate finite-difference scheme, the adoption of
the SVV method is straightforward and comes with a negligible additional computational
cost. In order to assess the ability of this high-order finite-difference spectral vanishing
viscosity approach, we performed large-eddy simulations (LES) of a gravity current in
a channelised lock-exchange set-up with our SVV model and with the well-known explicit static and dynamic Smagorinsky sub-grid scale (SGS) models. The obtained data
are compared with a direct numerical simulation (DNS) based on more than 800 million
mesh nodes, and with experimental measurements. A framework for the energy budget
is introduced to investigate the behaviour of the gravity current. First, it is found that
the DNS is in good agreement with the experimental data for the evolution of the front
location and velocity field as well as for the stirring and mixing inside the gravity current. Secondly, the LES performed with less than 0.4% of the total number of mesh nodes
compared to the DNS, can reproduce the main features of the gravity currents, with the
SVV model yielding slightly more accurate results. It is also found that the dynamic
Smagorinsky model performs better than its static version. For the present study, the
static and dynamic Smagorinsky models are 1.8 and 2.5 times more expensive than the
SVV model, because the latter does not require the calculation of explicit SGS terms in
the Navier-Stokes equations nor spatial filtering operations.
Date Issued
2021-05-15
Date Acceptance
2021-02-24
Citation
Computers and Fluids, 2021, 221, pp.1-18
ISSN
0045-7930
Publisher
Elsevier
Start Page
1
End Page
18
Journal / Book Title
Computers and Fluids
Volume
221
Copyright Statement
© 2021 Elsevier Ltd. All rights reserved. This manuscript is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International Licence http://creativecommons.org/licenses/by-nc-nd/4.0/
Identifier
https://www.sciencedirect.com/science/article/pii/S0045793021000682?via%3Dihub
Subjects
Applied Mathematics
0102 Applied Mathematics
0913 Mechanical Engineering
0915 Interdisciplinary Engineering
Publication Status
Published
Date Publish Online
2021-03-03