4.7 Article

A parameter-free dynamic alternative to hyper-viscosity for coupled transport equations: Application to the simulation of 3D squall lines using spectral elements

Journal

JOURNAL OF COMPUTATIONAL PHYSICS
Volume 283, Issue -, Pages 360-373

Publisher

ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2014.11.046

Keywords

Residual-based stabilization; Variational multiscalemethod; VMS; VMM; Crosswind discontinuity capturing; Dynamic artificial diffusion; Hyper-viscosity; Spectral element method; Monotonicity-preserving high-order methods; Stabilization of spectral elements; Deep convection; Kessler microphysics

Funding

  1. Office of Naval Research [PE-0602435N]
  2. National Science Foundation (Division of Mathematical Sciences) [121670]

Ask authors/readers for more resources

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) (Marras et al. (2012) [7]) 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. (C) 2014 Elsevier Inc. All rights reserved.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available