A parameter-free dynamic alternative to hyper-viscosity for coupled transport equations: application to the simulation of 3D squall lines using spectral elements. (English) Zbl 1351.86003
Summary: The stabilization of high order spectral elements to solve the transport equations for tracers in the atmosphere remains an active topic of research among atmospheric modelers. This paper builds on our previous work on variational multiscale stabilization (VMS) and discontinuity capturing (DC) [ibid. 231, No. 21, 7187–7213 (2012; Zbl 1284.65119)] and shows the applicability of VMS+DC to realistic atmospheric problems that involve physics coupling with phase change in the simulation of 3D deep convection. We show that the VMS+DC approach is a robust technique that can damp the high order modes characterizing the spectral element solution of complex coupled transport problems. The method has important properties that techniques of more common use often lack: 1) it is free of a user-defined parameter, 2) it is anisotropic in that it only acts along the flow direction, 3) it is numerically consistent, and 4) it can improve the monotonicity of high-order spectral elements. The proposed method is assessed by comparing the results against those obtained with a fourth-order hyper-viscosity programmed in the same code. The main conclusion that arises is that tuning can be fully avoided without loss of accuracy if the dissipative scheme is properly designed. Finally, the cost of parallel communication is that of a second order operator which means that fewer communications are required by VMS+DC than by a hyper-viscosity method; fewer communications translate into a faster and more scalable code, which is of vital importance as we approach the exascale range of computing.
MSC:
| 86-08 | Computational methods for problems pertaining to geophysics |
| 65M70 | Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs |
| 86A05 | Hydrology, hydrography, oceanography |
Keywords:
residual-based stabilization; variational multiscale method; crosswind discontinuity capturing; dynamic artificial diffusion; hyper-viscosity; spectral element method; monotonicity-preserving high-order methods; stabilization of spectral elements; deep convection; Kessler microphysicsCitations:
Zbl 1284.65119Software:
AlyaReferences:
| [1] | Jablonowski, C.; Williamson, D. L., The pros and cons of diffusion, filters and fixers in atmospheric general circulation models, (Lauritzen, P. H.; Jablonowski, C.; Taylor, M. A.; Nair, R. D., Numerical Techniques for Global Atmospheric Models. Numerical Techniques for Global Atmospheric Models, Lecture Notes in Computational Science and Engineering, vol. 80 (2011), Springer), 381-482 · Zbl 1215.86001 |
| [2] | Baruzzi, G.; Habashi, W.; Hefez, N., A second order accurate finite element method for the solutions of the Euler and Navier-Stokes equations, (Proceedings of the 13th International Conference on Numerical Methods in Fluid Dynamics (1992), Springer: Springer Rome), 509-513 |
| [3] | Codina, R., Stabilized finite element approximation of transient incompressible flows using orthogonal subscales, Comput. Methods Appl. Mech. Eng., 191, 4295-4321 (2002) · Zbl 1015.76045 |
| [4] | Hughes, T. J.R.; Brooks, A. N., A multidimensional upwind scheme with no crosswind diffusion, (Hughes, T. J.R., Finite Element Methods for Convection Dominated Flows, vol. 32 (1979), ASME), 19-35 · Zbl 0423.76067 |
| [5] | Hughes, T. J.R.; Mallet, M.; Mizukami, A., A new finite element formulation for computational fluid dynamics: II. Beyond SUPG, Comput. Methods Appl. Mech. Eng., 54, 341-355 (1986) · Zbl 0622.76074 |
| [6] | Codina, R., A discontinuity-capturing crosswind-dissipation for the finite element solution of the convection-diffusion equation, Comput. Methods Appl. Mech. Eng., 110, 325-342 (1993) · Zbl 0844.76048 |
| [7] | Marras, S.; Kelly, J. F.; Giraldo, F. X.; Vázquez, M., Variational multiscale stabilization of high-order spectral elements for the advection-diffusion equation, J. Comput. Phys., 231, 7187-7213 (2012) · Zbl 1284.65119 |
| [8] | Vandeven, H., Family of spectral filters for discontinuous problems, J. Sci. Comp., 6, 159-192 (1991) · Zbl 0752.35003 |
| [9] | Boyd, J. P., Two comments on filtering for Chebyshev and Legendre spectral and spectral element methods, J. Comput. Phys., 143, 283-288 (1998) · Zbl 0920.65046 |
| [10] | Malm, J.; Schlatter, P.; Fischer, P. F.; Henningson, D. S., Stabilization of the spectral element method in convection dominated flows by recovery of skew-symmetry, J. Sci. Comput., 57, 254-277 (2013) · Zbl 1282.76145 |
| [11] | Hughes, T., Multiscale phenomena: Green’s functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods, Comput. Methods Appl. Mech. Eng., 127, 387-401 (1995) · Zbl 0866.76044 |
| [12] | Gjesdal, T.; Wasberg, C. E.; Reif, B. A.P.; Andreassen, Ø., Variational multiscale turbulence modelling in a high order spectral element method, J. Comput. Phys., 228, 7333-7356 (2009) · Zbl 1172.76021 |
| [13] | Marras, S.; Moragues, M.; Vázquez, M. R.; Jorba, O.; Houzeaux, G., Simulations of moist convection by a variational multiscale stabilized finite element method, J. Comput. Phys., 252, 195-218 (2013) · Zbl 1349.86037 |
| [14] | Kelly, J. F.; Giraldo, F. X., Continuous and discontinuous Galerkin methods for a scalable three-dimensional nonhydrostatic atmospheric model: limited-area mode, J. Comput. Phys., 231, 7988-8008 (2012) · Zbl 1284.65134 |
| [15] | Giraldo, F. X.; Kelly, J. F.; Constantinescu, E., Implicit-explicit formulations of a three-dimensional Nonhydrostatic Unified Model of the Atmosphere (NUMA), SIAM J. Sci. Comput., 35, 1162-1194 (2013) · Zbl 1280.86008 |
| [16] | Marchuk, G. I., Numerical Methods in Weather Prediction (1974), Academic Press |
| [17] | Kessler, E., On the distribution and continuity of water substance in atmospheric circulation, Meteorol. Monogr., 10, 32 (1969) |
| [18] | Soong, S.; Ogura, Y., A comparison between axisymmetric and slab-symmetric cumulus cloud models, J. Atmos. Sci., 30, 879-893 (1973) |
| [19] | Knobloch, P., On the definition of the SUPG parameter, Electron. Trans. Numer. Anal., 32, 76-89 (2008) · Zbl 1171.65079 |
| [20] | Owen, H.; Houzeaux, G.; Samaniego, C.; Lesage, A. C.; Vázquez, M., Recent ship hydrodynamics developments in the parallel two-fluid flow solver Alya, Comput. Fluids, 80, 168-177 (2013) · Zbl 1284.76073 |
| [21] | Houzeaux, G.; Eguzkitza, B.; Vázquez, M., A variational multiscale model for the advection-diffusion-reaction equation, Commun. Numer. Methods Eng., 25, 787-809 (2009) · Zbl 1168.65413 |
| [22] | Guba, O.; Taylor, M. A.; Ullrich, P. A.; Overfelt, J. R.; Levy, M. N., The spectral element method on variable resolution grids: evaluating grid sensitivity and resolution-aware numerical viscosity, Geosci. Model Dev., 7, 4081-4117 (2014) |
| [23] | Weisman, M. L.; Klemp, J. B.; Rotunno, R., Structure and evolution of numerically simulated squall lines, J. Atmos. Sci., 45, 1990-2013 (1988) |
| [24] | Frame, J.; Markowski, P., The interaction of simulated squall lines with idealized mountain ridges, Mon. Weather Rev., 134, 1919-1941 (2006) |
| [25] | Weisman, M. L.; Klemp, J. B., The dependence of numerically simulated convective storms on vertical wind shear and buoyancy, Mon. Weather Rev., 110, 504-520 (1982) |
| [26] | Weisman, M. L.; Rotunno, R., “A theory for strong long-lived squall lines” revisited, J. Atmos. Sci., 61, 361-382 (2004) |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
