A nonnegative and shape‐preserving global transport model on cubed sphere using high‐order conservative collocation scheme

Author(s):  
Pei Huang ◽  
Chungang Chen ◽  
Xingliang Li ◽  
Xueshun Shen ◽  
Feng Xiao
2018 ◽  
Vol 144 (714) ◽  
pp. 1611-1627 ◽  
Author(s):  
Jie Tang ◽  
Chungang Chen ◽  
Xingliang Li ◽  
Xueshun Shen ◽  
Feng Xiao

2012 ◽  
Vol 4 (11) ◽  
pp. 3528-3543 ◽  
Author(s):  
Nick Schutgens ◽  
Makiko Nakata ◽  
Teruyuki Nakajima

2010 ◽  
Vol 67 (12) ◽  
pp. 1993-2014 ◽  
Author(s):  
C. G. Chen ◽  
F. Xiao ◽  
X. L. Li ◽  
Y. Yang
Keyword(s):  

2015 ◽  
Vol 143 (7) ◽  
pp. 2937-2954 ◽  
Author(s):  
Kiran K. Katta ◽  
Ramachandran D. Nair ◽  
Vinod Kumar

Abstract This paper presents two finite-volume (FV) schemes for solving linear transport problems on the cubed-sphere grid system. The schemes are based on the central-upwind finite-volume (CUFV) method, which is a class of Godunov-type method for solving hyperbolic conservation laws, and combines the attractive features of the classical upwind and central FV methods. One of the CUFV schemes is based on a dimension-by-dimension approach and employs a fifth-order one-dimensional (1D) Weighted Essentially Nonoscillatory (WENO5) reconstruction method. The other scheme employs a fully two-dimensional (2D) fourth-order accurate reconstruction method. The cubed-sphere grid system imposes several computational challenges due to its patched-domain topology and nonorthogonal curvilinear grid structure. A high-order 1D interpolation procedure combining cubic and quadratic interpolations is developed for the FV schemes to handle the discontinuous edges of the cubed-sphere grid. The WENO5 scheme is compared against the fourth-order Kurganov–Levy (KL) scheme formulated in the CUFV framework. The performance of the schemes is compared using several benchmark problems such as the solid-body rotation and deformational-flow tests, and empirical convergence rates are reported. In addition, a bound-preserving filter combined with an optional positivity-preserving filter is tested for nonsmooth problems. The filtering techniques considered are local, inexpensive, and effective. A fourth-order strong stability preserving explicit Runge–Kutta time-stepping scheme is used for integration. The results show that schemes are competitive to other published FV schemes in the same category.


2018 ◽  
Vol 146 (7) ◽  
pp. 2047-2064 ◽  
Author(s):  
Hyun-Gyu Kang ◽  
Hyeong-Bin Cheong

Abstract A high-order filter for a cubed-sphere spectral element model was implemented in a three-dimensional spectral element dry hydrostatic dynamical core. The dynamical core incorporated hybrid sigma–pressure vertical coordinates and a third-order Runge–Kutta time-differencing method. The global high-order filter and the local-domain high-order filter, requiring numerical operation with a huge sparse global matrix and a locally assembled matrix, respectively, were applied to the prognostic variables, except for surface pressure, at every time step. Performance of the high-order filter was evaluated using the baroclinic instability test and quiescent atmosphere with underlying topography test presented by the Dynamical Core Model Intercomparison Project. It was revealed that both the global and local-domain high-order filters could better control the numerical noise in the noisy circumstances than the explicit diffusion, which is widely used for the spectral element dynamical core. Furthermore, by adopting the high-order filter, the effective resolution of the dynamical core could be increased, without weakening the stability of the dynamical core. Computational efficiency of the high-order filter was demonstrated in terms of both the time step size and the wall-clock time. Because of the nature of an implicit diffusion, the dynamical core employing this filter can take a larger time step size, compared to that using the explicit diffusion. The local-domain high-order filter was computationally more efficient than the global high-order filter, but less efficient than the explicit diffusion.


2011 ◽  
Vol 24 (15) ◽  
pp. 4037-4055 ◽  
Author(s):  
Saroj K. Mishra ◽  
Mark A. Taylor ◽  
Ramachandran D. Nair ◽  
Peter H. Lauritzen ◽  
Henry M. Tufo ◽  
...  

Abstract The NCAR Community Climate System Model, version 4 (CCSM4), includes a new dynamical core option based on NCAR’s High-Order Method Modeling Environment (HOMME). HOMME is a petascale-capable high-order element-based conservative dynamical core developed on the cubed-sphere grid. Initial simulations have been completed in an aquaplanet configuration of the Community Atmosphere Model, version 4 (CAM4), the atmospheric component of CCSM4. The authors examined the results of this simulation and assessed its fidelity in simulating rainfall, which is one of the most important components of the earth’s climate system. For this they compared the results from two other dynamical cores of CAM4: the finite volume (FV) and Eulerian (EUL). Instantaneous features of rainfall in HOMME are similar to FV and EUL. Similar to EUL and FV, HOMME simulates a single-peak intertropical convergence zone (ITCZ) over the equator. The strength of the ITCZ is found to be almost the same in HOMME and EUL but more in FV. It is observed that in HOMME and EUL, there is higher surface evaporation, which supplies more moisture to the deep tropics and gives more rainfall over the ITCZ. The altitude of maximum precipitation is found to be at almost the same level in all three dynamical cores. The eastward propagation of rainfall bands is organized and more prominent in FV and HOMME than in EUL. The phase speed of the eastward propagation in HOMME is found to be higher than in FV. The results show that, in general, the rainfall simulated by HOMME falls in a regime between that of FV and EUL. Hence, they conclude that the key aspects of rainfall simulation with HOMME falls into an acceptable range, as compared to the existing dynamical cores used in the model.


Sign in / Sign up

Export Citation Format

Share Document