scholarly journals DYNAMICO, an icosahedral hydrostatic dynamical core designed for consistency and versatility

2015 ◽  
Vol 8 (2) ◽  
pp. 1749-1800 ◽  
Author(s):  
T. Dubos ◽  
S. Dubey ◽  
M. Tort ◽  
R. Mittal ◽  
Y. Meurdesoif ◽  
...  
Keyword(s):  
Hamiltonian Structure ◽  
Potential Temperature ◽  
Time Integration ◽  
Hamiltonian Formulation ◽  
Three Dimensional ◽  
Primitive Equations ◽  
Finite Volume Scheme ◽  
Dynamical Core ◽  
Vertical Coordinate ◽  
Model Code

Abstract. The design of the icosahedral dynamical core DYNAMICO is presented. DYNAMICO solves the multi-layer rotating shallow-water equations, a compressible variant of the same equivalent to a discretization of the hydrostatic primitive equations in a Lagrangian vertical coordinate, and the primitive equations in a hybrid mass-based vertical coordinate. The common Hamiltonian structure of these sets of equations is exploited to formulate energy-conserving spatial discretizations in a unified way. The horizontal mesh is a quasi-uniform icosahedral C-grid obtained by subdivision of a regular icosahedron. Control volumes for mass, tracers and entropy/potential temperature are the hexagonal cells of the Voronoi mesh to avoid the fast numerical modes of the triangular C-grid. The horizontal discretization is that of Ringler et al. (2010), whose discrete quasi-Hamiltonian structure is identified. The prognostic variables are arranged vertically on a Lorenz grid with all thermodynamical variables collocated with mass. The vertical discretization is obtained from the three-dimensional Hamiltonian formulation. Tracers are transported using a second-order finite volume scheme with slope limiting for positivity. Explicit Runge–Kutta time integration is used for dynamics and forward-in-time integration with horizontal/vertical splitting is used for tracers. Most of the model code is common to the three sets of equations solved, making it easier to develop and validate each piece of the model separately. Representative three-dimensional test cases are run and analyzed, showing correctness of the model. The design permits to consider several extensions in the near future, from higher-order transport to more general dynamics, especially deep-atmosphere and non-hydrostatic equations.

scholarly journals DYNAMICO-1.0, an icosahedral hydrostatic dynamical core designed for consistency and versatility

2015 ◽  
Vol 8 (10) ◽  
pp. 3131-3150 ◽  
Author(s):  
T. Dubos ◽  
S. Dubey ◽  
M. Tort ◽  
R. Mittal ◽  
Y. Meurdesoif ◽  
...  
Keyword(s):  
Hamiltonian Structure ◽  
Potential Temperature ◽  
Time Integration ◽  
Hamiltonian Formulation ◽  
Three Dimensional ◽  
Primitive Equations ◽  
Finite Volume Scheme ◽  
Dynamical Core ◽  
Vertical Coordinate ◽  
Model Code

Abstract. The design of the icosahedral dynamical core DYNAMICO is presented. DYNAMICO solves the multi-layer rotating shallow-water equations, a compressible variant of the same equivalent to a discretization of the hydrostatic primitive equations in a Lagrangian vertical coordinate, and the primitive equations in a hybrid mass-based vertical coordinate. The common Hamiltonian structure of these sets of equations is exploited to formulate energy-conserving spatial discretizations in a unified way. The horizontal mesh is a quasi-uniform icosahedral C-grid obtained by subdivision of a regular icosahedron. Control volumes for mass, tracers and entropy/potential temperature are the hexagonal cells of the Voronoi mesh to avoid the fast numerical modes of the triangular C-grid. The horizontal discretization is that of Ringler et al. (2010), whose discrete quasi-Hamiltonian structure is identified. The prognostic variables are arranged vertically on a Lorenz grid with all thermodynamical variables collocated with mass. The vertical discretization is obtained from the three-dimensional Hamiltonian formulation. Tracers are transported using a second-order finite-volume scheme with slope limiting for positivity. Explicit Runge–Kutta time integration is used for dynamics, and forward-in-time integration with horizontal/vertical splitting is used for tracers. Most of the model code is common to the three sets of equations solved, making it easier to develop and validate each piece of the model separately. Representative three-dimensional test cases are run and analyzed, showing correctness of the model. The design permits to consider several extensions in the near future, from higher-order transport to more general dynamics, especially deep-atmosphere and non-hydrostatic equations.

scholarly journals A Three-Dimensional Wave-Activity Relation for Pseudomomentum

2007 ◽  
Vol 64 (6) ◽  
pp. 2126-2134 ◽  
Author(s):  
Lingkun Ran ◽  
Shouting Gao
Keyword(s):  
Potential Temperature ◽  
Three Dimensional ◽  
Wave Activity ◽  
Primitive Equations ◽  
Casimir Function ◽  
Activity Density ◽  
Atmospheric Data ◽  
Local Wave ◽  
Basic States ◽  
Wave Activity Flux

A three-dimensional, nonhydrostatic local wave-activity relation for pseudomomentum is derived from the nonhydrostatic primitive equations in Cartesian coordinates by using an extension of the momentum–Casimir method. The stationary and zonally symmetric basic states are chosen and a Casimir function, which is the single-valued function of potential vorticity and potential temperature, is introduced in the derivation. The wave-activity density and wave-activity flux of the local wave-activity relation for pseudomomentum are expressed entirely in terms of Eulerian quantities so that they are easily calculated with atmospheric data and do not require the knowledge of particle placements. Constructed in the ageostrophic and nonhydrostatic dynamical framework, the local wave-activity relation for pseudomomentum is applicable to diagnosing the evolution and propagation of mesoscale weather systems.

scholarly journals WAVETRISK-1.0: an adaptive wavelet hydrostatic dynamical core

2019 ◽  
Author(s):  
Nicholas K.-R. Kevlahan ◽  
Thomas Dubos
Keyword(s):  
Shallow Water ◽  
General Circulation ◽  
Wave Train ◽  
Shallow Water Equation ◽  
Three Dimensional ◽  
Circulation Model ◽  
Time Step ◽  
Adaptive Wavelet ◽  
Dynamical Core ◽  
Vertical Coordinate

Abstract. This paper presents the new adaptive dynamical core wavetrisk. The fundamental features of the wavelet-based adaptivity were developed for the shallow water equation on the β-plane in Dubos and Kevlahan (2013) and extended to the icosahedral grid on the sphere in Aechtner et al. (2015). The three-dimensional dynamical core solves the compressible hydrostatic multilayer rotating shallow water equations on a multiscale dynamically adapted grid. The equations are discretized using a Lagrangian vertical coordinate version of dynamico introduced in Dubos et al. (2015). The horizontal computational grid is adapted at each time step to ensure a user-specified relative error in either the tendencies or the solution. The Lagrangian vertical grid is remapped using an adaptive Lagrangian-Eulerian (ALE) algorithm onto the initial hybrid σ pressure-based coordinates as necessary. The resulting grid is adapted horizontally, but uniform over all vertical layers. Thus, the three-dimensional grid is a set of columns of varying sizes. The code is parallelized by domain decomposition using mpi and the variables are stored in a hybrid data structure of dyadic quad trees and patches. A low storage explicit fourth order Runge-Kutta scheme is used for time integration. Validation results are presented for three standard dynamical core test cases: mountain-induced Rossby wave train, baroclinic instability of a jet stream and the Held and Suarez simplified general circulation model. The results confirm good strong parallel scaling and demonstrate that wavetrisk can achieve grid compression ratios of several hundred times compared with an equivalent static grid model.

A Supercell Storm Simulation Using a Nonhydrostatic Cloud-Resolving Model Based on a Hybrid Isentropic-Sigma Vertical Coordinate

2013 ◽  
Vol 141 (4) ◽  
pp. 1204-1215 ◽  
Author(s):  
Michael D. Toy
Keyword(s):  
Mass Flux ◽  
Potential Temperature ◽  
Three Dimensional ◽  
Vertical Motion ◽  
Implicit Time ◽  
Vertical Coordinate ◽  
Model Based ◽  
Supercell Storm ◽  
Coordinate Model ◽  
Isentropic Coordinates

Abstract A three-dimensional simulation of a supercell storm is performed with a nonhydrostatic model based on a hybrid isentropic-sigma vertical coordinate. The coordinate is a terrain-following, height-based coordinate near the surface that smoothly transitions to potential temperature with height. Using isentropic coordinates provides the advantage of having zero cross-coordinate vertical mass flux for adiabatic flow, which virtually eliminates the numerical error in the vertical transport. The model uses an adaptive grid algorithm by which the coordinate surfaces may deviate from their target isentropes to maintain a sufficiently smooth mesh, while allowing the turbulence and vertical motion associated with convection to develop. The storm simulated by the hybrid-coordinate model compares well with simulations by Eulerian-coordinate models, but with the key difference being that the cross-coordinate mass flux is significantly smaller in much of the domain with the hybrid-coordinate model. A semi-implicit time-differencing scheme for numerically stabilizing vertically propagating acoustic modes in isentropic coordinates is also presented in the paper.

scholarly journals A New Dynamical Core of the Global Environmental Multiscale (GEM) Model with a Height-Based Terrain-Following Vertical Coordinate

2019 ◽  
Vol 147 (7) ◽  
pp. 2555-2578 ◽  
Author(s):  
Syed Zahid Husain ◽  
Claude Girard ◽  
Abdessamad Qaddouri ◽  
André Plante
Keyword(s):  
Numerical Experiments ◽  
Three Dimensional ◽  
Two Dimensional ◽  
Dynamical Core ◽  
Vertical Coordinate ◽  
Solution Approach ◽  
Global Environmental ◽  
Wide Range ◽  
Terrain Following ◽  
Global Environmental Multiscale

Abstract A new dynamical core of Environment and Climate Change Canada’s Global Environmental Multiscale (GEM) atmospheric model is presented. Unlike the existing log-hydrostatic-pressure-type terrain-following vertical coordinate, the proposed core adopts a height-based approach. The move to a height-based vertical coordinate is motivated by its potential for improving model stability over steep terrain, which is expected to become more prevalent with the increasing demand for very high-resolution forecasting systems. A dynamical core with height-based vertical coordinate generally requires an iterative solution approach. In addition to a three-dimensional iterative solver, a simplified approach has been devised allowing the use of a direct solver for the new dynamical core that separates a three-dimensional elliptic boundary value problem into a set of two-dimensional independent Helmholtz problems. The issue of dynamics–physics coupling has also been studied, and incorporating the physics tendencies within the discretized dynamical equations is found to be the most acceptable approach for the height-based vertical coordinate. The new dynamical core is evaluated using numerical experiments that include two-dimensional nonhydrostatic theoretical cases as well as 25-km resolution global forecasts. For a wide range of horizontal grid resolutions—from a few meters to up to 25 km—the results from the direct solution approach are found to be equivalent to the iterative approach for the new dynamical core. Furthermore, results from the different numerical experiments confirm that the new height-based dynamical core is equivalent to the existing pressure-based core in terms of solution accuracy.

scholarly journals Design of a Dynamical Core Based on the Nonhydrostatic “Unified System” of Equations*

2014 ◽  
Vol 142 (1) ◽  
pp. 364-385 ◽  
Author(s):  
Celal S. Konor
Keyword(s):  
Acoustic Waves ◽  
Potential Temperature ◽  
Time Derivative ◽  
System Of Equations ◽  
Dynamical Core ◽  
Vertical Coordinate ◽  
Discrete Equations ◽  
Wide Range ◽  
Supplementary Material ◽  
Vertically Propagating

Abstract This paper presents the design of a dry dynamical core based on the nonhydrostatic “unified system” of equations. The unified system filters vertically propagating acoustic waves. The dynamical core predicts the potential temperature and horizontal momentum. It uses the predicted potential temperature to determine the quasi-hydrostatic components of the Exner pressure and density. The continuity equation is diagnostic (and used to determine the vertical mass flux) because the time derivative of the quasi-hydrostatic density is obtained from the predicted potential temperature. The nonhydrostatic component of the Exner pressure is obtained from an elliptic equation. The main focus of this paper is on the integration procedure used with this unique dynamical core. In the implementation described in this paper, height is used as the vertical coordinate, and the equations are vertically discretized on a Lorenz-type grid. Cartesian horizontal coordinates are used along with an Arakawa C grid. A detailed description of the discrete equations is presented in the supplementary material, along with the rationale behind the decisions made during the discretization process. Only a short description is given here. To demonstrate that the model is capable of simulating a wide range of dynamical scales, the results from cyclone- and cloud-scale simulations are presented. The solutions obtained for the selected cloud-scale simulations are compared to those from a fully compressible, anelastic, and pseudo-incompressible models that (as far as possible) uses the same schemes used in the unified dynamical core. The results show that the unified dynamical core performs reasonably well in all these experiments.

scholarly journals Implementation of a semi-Lagrangian fully-implicit time integration of the unified soundproof system of equations for numerical weather prediction

2021 ◽  
Author(s):  
Abdessamad Qaddouri ◽  
Claude Girard ◽  
Syed Zahid Husain ◽  
Rabah Aider
Keyword(s):  
Numerical Weather Prediction ◽  
Weather Prediction ◽  
Time Integration ◽  
Implicit Time ◽  
System Of Equations ◽  
Dynamical Core ◽  
Vertical Coordinate ◽  
Numerical Weather ◽  
Fully Implicit ◽  
Computational Performance

AbstractAn alternate dynamical core that employs the unified equations of A. Arakawa and C.S. Konor (2009) has been developed within Environment and Climate change Canada’s GEM (Global Environmental Multiscale) atmospheric model. As in the operational GEM dynamical core, the novel core utilizes the same fully-implicit two-time-level semi-Lagrangian scheme for time discretization while the log-pressure-based terrain-following vertical coordinate has been slightly adapted. Overall, the new dynamical core implementation required only minor changes to the existing informatics code of the GEM model and from a computational performance perspective, the new core does not incur any significant additional cost. A broad range of tests – that include both two-dimensional idealized theoretical cases and three-dimensional deterministic forecasts covering both hydrostatic and non-hydrostatic scales–have been carried out to evaluate the performance of the new dynamical core. For all the tested cases, when compared to the operational GEM model, the new dynamical core based on the unified equations has been found to produce statistically equivalent results. These results imply that the unified equations can be adopted for operational numerical weather prediction that would employ a single soundproof system of equations to produce reliable forecasts for all meteorological scales of interest with negligible changes for the computational overhead.

scholarly journals A Fully Implicit Finite Volume Scheme for a Seawater Intrusion Problem in Coastal Aquifers

Water ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 1639
Author(s):  
Abdelkrim Aharmouch ◽  
Brahim Amaziane ◽  
Mustapha El Ossmani ◽  
Khadija Talali
Keyword(s):  
Finite Volume ◽  
Coastal Aquifers ◽  
Nonlinear Partial Differential Equations ◽  
Time Integration ◽  
Mass Conservation ◽  
Coupled System ◽  
Finite Volume Scheme ◽  
Time Step ◽  
Volume Scheme ◽  
Fully Implicit

We present a numerical framework for efficiently simulating seawater flow in coastal aquifers using a finite volume method. The mathematical model consists of coupled and nonlinear partial differential equations. Difficulties arise from the nonlinear structure of the system and the complexity of natural fields, which results in complex aquifer geometries and heterogeneity in the hydraulic parameters. When numerically solving such a model, due to the mentioned feature, attempts to explicitly perform the time integration result in an excessively restricted stability condition on time step. An implicit method, which calculates the flow dynamics at each time step, is needed to overcome the stability problem of the time integration and mass conservation. A fully implicit finite volume scheme is developed to discretize the coupled system that allows the use of much longer time steps than explicit schemes. We have developed and implemented this scheme in a new module in the context of the open source platform DuMu X . The accuracy and effectiveness of this new module are demonstrated through numerical investigation for simulating the displacement of the sharp interface between saltwater and freshwater in groundwater flow. Lastly, numerical results of a realistic test case are presented to prove the efficiency and the performance of the method.

Numerical study of multi-dimensional hyperbolic telegraph equations arising in nuclear material science via an efficient local meshless method

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Imtiaz Ahmad ◽  
Aly R. Seadawy ◽  
Hijaz Ahmad ◽  
Phatiphat Thounthong ◽  
Fuzhang Wang
Keyword(s):  
Meshless Method ◽  
Material Science ◽  
Numerical Study ◽  
Research Work ◽  
Time Integration ◽  
Three Dimensional ◽  
Nuclear Material ◽  
Telegraph Equations ◽  
Model Equations ◽  
Derivatives Of

Abstract This research work is to study the numerical solution of three-dimensional second-order hyperbolic telegraph equations using an efficient local meshless method based on radial basis function (RBF). The model equations are used in nuclear material science and in the modeling of vibrations of structures. The explicit time integration technique is utilized to semi-discretize the model in the time direction whereas the space derivatives of the model are discretized by the proposed local meshless procedure based on multiquadric RBF. Numerical experiments are performed with the proposed numerical scheme for rectangular and non-rectangular computational domains. The proposed method solutions are converging quickly in comparison with the different existing numerical methods in the recent literature.

Sign in / Sign up

Export Citation Format

Share Document