Incorporating Condensational Heating into a Nonhydrostatic Atmospheric Model Based on a Hybrid Isentropic–Sigma Vertical Coordinate

2011 ◽  
Vol 139 (9) ◽  
pp. 2940-2954 ◽  
Author(s):  
Michael D. Toy

Using isentropic coordinates in atmospheric models has the advantage of eliminating the cross-coordinate vertical mass flux for adiabatic flow, and virtually eliminating the associated numerical error in the vertical transport. This is a significant benefit since much of the flow in the atmosphere is approximately adiabatic. Nonadiabatic processes, such as condensational heating, result in a nonzero vertical velocity [Formula: see text] in isentropic coordinates. A method for incorporating condensational heating into a nonhydrostatic atmospheric model based on a hybrid isentropic–sigma vertical coordinate is presented. The model is tested with various 2D moist simulations and the results are compared with those using a traditional terrain-following, height-based sigma coordinate. With the hybrid coordinate, there are improvements in the representation of the developing cloud field in a mountain wave experiment. In a simulation of deep convection, the adaptive hybrid coordinate successfully simulates the turbulent nature of the convection, while maintaining the quasi-Lagrangian nature of the isentropic coordinate in the surrounding dry air. The vertical cross-coordinate mass flux is almost zero in the environmental air, as well as in the stratosphere above the convective tower.

2013 ◽  
Vol 141 (4) ◽  
pp. 1204-1215 ◽  
Author(s):  
Michael D. Toy

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.


2009 ◽  
Vol 137 (7) ◽  
pp. 2305-2330 ◽  
Author(s):  
Michael D. Toy ◽  
David A. Randall

The isentropic system of equations has particular advantages in the numerical modeling of weather and climate. These include the elimination of the vertical velocity in adiabatic flow, which simplifies the motion to a two-dimensional problem and greatly reduces the numerical errors associated with vertical advection. The mechanism for the vertical transfer of horizontal momentum is simply the pressure drag acting on isentropic coordinate surfaces under frictionless, adiabatic conditions. In addition, vertical resolution is enhanced in regions of high static stability, which leads to better resolution of features such as the tropopause. Negative static stability and isentropic overturning frequently occur in finescale atmospheric motion. This presents a challenge to nonhydrostatic modeling with the isentropic vertical coordinate. This paper presents a new nonhydrostatic atmospheric model based on a generalized vertical coordinate. The coordinate is specified in a manner similar to that of Konor and Arakawa, but “arbitrary Eulerian–Lagrangian” (ALE) methods are used to maintain coordinate monotonicity in regions of negative static stability and to return the coordinate surfaces to their isentropic “targets” in statically stable regions. The model is mass conserving and implements a vertical differencing scheme that satisfies two additional integral constraints for the limiting case of z coordinates. The hybrid vertical coordinate model is tested with mountain-wave experiments including a downslope windstorm with breaking gravity waves. The results show that the advantages of the isentropic coordinate are realized in the model with regard to vertical tracer and momentum transport. Also, the isentropic overturning associated with the wave breaking is successfully handled by the coordinate formulation.


2010 ◽  
Vol 138 (9) ◽  
pp. 3683-3689 ◽  
Author(s):  
Daniel Leuenberger ◽  
Marcel Koller ◽  
Oliver Fuhrer ◽  
Christoph Schär

Abstract Most atmospheric models use terrain-following coordinates, and it is well known that the associated deformation of the computational mesh leads to numerical inaccuracies. In a previous study, the authors proposed a new terrain-following coordinate formulation [the smooth level vertical (SLEVE) coordinate], which yields smooth vertical coordinate levels at mid and upper levels and thereby considerably reduces numerical errors in the simulation of flow past complex topography. In the current paper, a generalization of the SLEVE coordinate is presented by using a modified vertical decay of the topographic signature with height. The new formulation enables an almost uniform thickness of the lowermost computational layers, while preserving the fast transition to smooth levels in the mid and upper atmosphere. This allows for a more consistent and more stable coupling with planetary boundary layer schemes, while retaining the advantages over classic sigma coordinates at upper levels. The generalized SLEVE coordinate is implemented and successfully tested in real-case simulations using an operational nonhydrostatic atmospheric model.


Icarus ◽  
2006 ◽  
Vol 182 (1) ◽  
pp. 259-273 ◽  
Author(s):  
Timothy E. Dowling ◽  
Mary E. Bradley ◽  
Edward Colón ◽  
John Kramer ◽  
Raymond P. LeBeau ◽  
...  

2021 ◽  
Vol 4 (1) ◽  
Author(s):  
Wen-yih Sun

In the terrain following coordinate, Gal-Chen and Somerville (1975) and other proposed a vertical coordinate  z*=(z-zb)/(zt-zb) and constant spatial intervals of dx* and  dy*along the other directions.  Because the variation of  and  was ignored, their coordinate does not really follow the terrain.  It fails to reproduce the divergence and curl over a complex terrain.  Aligning the coordinate with real terrain, the divergence and curl we obtained from the curvilinear coordinate are consistent with the Cartesian coordinate.  With a modification, the simulated total mass, energy, and momentum from the Navier-Stokes equations are conserved and in agreement with those calculated from Cartesian coordinate.


2021 ◽  
Vol 30 (1) ◽  
pp. 24-35
Author(s):  
Wenhui Cui ◽  
Wei Qu ◽  
Min Jiang ◽  
Gang Yao

Abstract Traditional atmospheric models are based on the analysis and fitting of various factors influencing the space atmosphere density. Neural network models do not specifically analyze the polynomials of each influencing factor in the atmospheric model, but use large data sets for network construction. Two traditional atmospheric model algorithms are analyzed, the main factors affecting the atmospheric model are identified, and an atmospheric model based on neural networks containing various influencing factors is proposed. According to the simulation error, the Levenberg-Marquardt algorithm is used to iteratively realize the rapid network weight correction, and the optimal neural network atmospheric model is obtained. The space atmosphere is simulated and calculated with an atmospheric model based on neural networks, and its average error rate is lower than that of traditional atmospheric models such as the DTM2013 model and the MSIS00 model. At the same time, the calculation complexity of the atmospheric model based on the neural networks is significantly simplified than that of the traditional atmospheric model.


Sign in / Sign up

Export Citation Format

Share Document