Depth-averaged unsteady RANS simulation of resonant shallow flows in lateral cavities using augmented WENO-ADER schemes

被引:27
作者
Navas-Montilla, A. [1 ]
Juez, C. [2 ]
Franca, M. J. [3 ,4 ]
Murillo, J. [1 ]
机构
[1] Univ Zaragoza, CSIC, Fluid Mech Dept LIFTEC, Zaragoza, Spain
[2] Univ Politecn Madrid, ETSIDI, Fluid Mech Dept, Madrid, Spain
[3] Delft Univ Technol, IHE Delft Inst Water Educ, Water Sci & Engn Dept, Delft, Netherlands
[4] Delft Univ Technol, Dept Hydraul Engn, Delft, Netherlands
关键词
Shallow flows; Turbulence; Shock capturing; Lateral cavities; Seiching; DISCONTINUOUS GALERKIN METHOD; HYPERBOLIC CONSERVATION-LAWS; GENERALIZED RIEMANN PROBLEM; SOURCE TERMS; EFFICIENT IMPLEMENTATION; WATER EQUATIONS; FINE SEDIMENT; OPEN-CHANNEL; SOLVERS; SYSTEMS;
D O I
10.1016/j.jcp.2019.06.037
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Turbulent shallow flows are characterized by the presence of horizontal large-scale vortices, caused by local variations of the velocity field. Apart from these 2D large vortices, small scale 3D turbulence, mainly produced by the interaction of the flowing water with the solid boundaries, is also present. The energy spectrum of turbulent shallow flows shows the presence of a 2D energy cascade at low wave numbers and a 3D energy cascade at high wave numbers, with a well-defined separation region between them. Horizontal flow movements (e.g. 2D large-scale vortical structures) at low wave numbers mostly determine the hydrodynamic behavior of these flows. Moreover, the generation of standing waves often occurs closely associated to the interaction of 2D horizontal flows with lateral boundaries, this is the case of seiches. To adequately reproduce these phenomena, a mathematical and numerical model able to resolve 2D turbulence is required. We herein show that depth-averaged (DA) unsteady Reynolds averaged Navier Stokes (URANS) models based on the Shallow Water Equations (SWE) are a suitable choice for the resolution of turbulent shallow flows with sufficient accuracy in an affordable computational time. The 3D small-scale vortices are modeled by means of diffusion terms, whereas the 2D large-scales are resolved. A high order numerical scheme is required for the resolution of 2D large eddies. In this work, we design a DA-URANS model based on a high order augmented WENO-ADER scheme. The mathematical model and numerical scheme are validated against observation of complex experiments in an open channel with lateral cavities that involve the presence of resonant phenomena (seiching). The numerical results evidence that the model accurately reproduces both longitudinal and transversal resonant waves and provides an accurate description of the flow field. The high order WENO-ADER scheme combined with a SWE model allows to obtain a powerful, reliable and efficient URANS simulation tool. (C) 2019 Elsevier Inc. All rights reserved.
引用
收藏
页码:511 / 536
页数:26
相关论文
共 52 条
[1]  
Akutina Y., 2015, THESIS
[2]  
[Anonymous], 2013, IAHR MONOGRAPH
[3]  
[Anonymous], 2004, Arch. Hydro-Eng. Environ. Mech.
[4]   Large eddy simulations for quasi-2D turbulence in shallow flows: A comparison between different subgrid scale models [J].
Awad, Esam ;
Toorman, Erik ;
Lacor, Chris .
JOURNAL OF MARINE SYSTEMS, 2009, 77 (04) :511-528
[5]   Efficient implementation of ADER schemes for Euler and magnetohydrodynamical flows on structured meshes - Speed comparisons with Runge-Kutta methods [J].
Balsara, Dinshaw S. ;
Meyer, Chad ;
Dumbser, Michael ;
Du, Huijing ;
Xu, Zhiliang .
JOURNAL OF COMPUTATIONAL PHYSICS, 2013, 235 :934-969
[6]  
Bouffard D, 2019, ANNU REV FLUID MECH, V51, P189, DOI [10.1146/annurev-fluid-010518-040506, 10.1146/annurev-fluid-010518040506]
[7]   Fourth-order balanced source term treatment in central WENO schemes for shallow water equations [J].
Caleffi, V. ;
Valiani, A. ;
Bernini, A. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2006, 218 (01) :228-245
[8]   A 2D local discontinuous Galerkin method for contaminant transport in channel bends [J].
Caleffi, Valerio ;
Valiani, Alessandro .
COMPUTERS & FLUIDS, 2013, 88 :629-642
[9]   Well-balanced high-order centred schemes for non-conservative hyperbolic systems. Applications to shallow water equations with fixed and mobile bed [J].
Canestrelli, Alberto ;
Siviglia, Annunziato ;
Dumbser, Michael ;
Toro, Eleuterio F. .
ADVANCES IN WATER RESOURCES, 2009, 32 (06) :834-844
[10]   Solvers for the high-order Riemann problem for hyperbolic balance laws [J].
Castro, C. E. ;
Toro, E. F. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2008, 227 (04) :2481-2513