An official website of the United States government
Abstract Bin microphysics schemes are useful tools for cloud simulations and are often considered to provide a benchmark for model intercomparison. However, they may experience issues with numerical diffusion, which are not well quantified, and the transport of hydrometeors depends on the choice of advection scheme, which can also change cloud simulation results. Here, an atmospheric large‐eddy simulation model is adapted to simulate a statistically steady‐state cloud in a convection cloud chamber under well‐constrained conditions. Two bin microphysics schemes, a spectral bin method and the method of moments, as well as several advection methods for the transport of the microphysical variables are employed for model intercomparison. Results show that different combinations of microphysics and advection schemes can lead to considerable differences in simulated cloud properties, such as cloud droplet number concentration. We find that simulations using the advection scheme that suffers more from numerical diffusion tends to have a smaller droplet number concentration and liquid water content, while simulation with the microphysics scheme that suffers more from numerical diffusion tends to have a broader size distribution and thus larger mean droplet sizes. Sensitivities of simulations to bin resolution, spatial resolution, and temporal resolution are also tested. We find that refining the microphysical bin resolution leads to a broader cloud droplet size distribution due to the advection of hydrometeors. Our results provide insight for using different advection and microphysics schemes in cloud chamber simulations, which might also help understand the uncertainties of the schemes used in atmospheric cloud simulations.
more »
« less
A Non‐Column Based, Fully Unstructured Implementation of Kessler's Microphysics With Warm Rain Using Continuous and Discontinuous Spectral Elements
Abstract Numerical weather prediction is pushing the envelope of grid resolution at local and global scales alike. Aiming to model topography with higher precision, a handful of articles introduced unstructured vertical grids and tested them for dry atmospheres. The next step toward effective high‐resolution unstructured grids for atmospheric modeling requires that also microphysics is independent of any vertical columns, in contrast to what is ubiquitous across operational and research models. In this paper, we present a non‐column based continuous and discontinuous spectral element implementation of Kessler's microphysics with warm rain. We test the proposed algorithm against standard three‐dimensional benchmarks for precipitating clouds and show that the results are comparable with those presented in the literature across all of the tested effective resolutions. While presented for both continuous and discontinuous spectral elements in this paper, the method that we propose can be adapted to any numerical method used in other codes, as long as the code can already handle vertically unstructured grids.
more »
« less
- Award ID(s):
- 1953535
- PAR ID:
- 10403628
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 15
- Issue:
- 3
- ISSN:
- 1942-2466
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract A local discontinuous Galerkin (LDG) method for approximating large deformations of prestrained plates is introduced and tested on several insightful numerical examples in Bonito et al. (2022, LDG approximation of large deformations of prestrained plates. J. Comput. Phys., 448, 110719). This paper presents a numerical analysis of this LDG method, focusing on the free boundary case. The problem consists of minimizing a fourth-order bending energy subject to a nonlinear and nonconvex metric constraint. The energy is discretized using LDG and a discrete gradient flow is used for computing discrete minimizers. We first show $$\varGamma $$-convergence of the discrete energy to the continuous one. Then we prove that the discrete gradient flow decreases the energy at each step and computes discrete minimizers with control of the metric constraint defect. We also present a numerical scheme for initialization of the gradient flow and discuss the conditional stability of it.more » « less
-
Summary We construct an algebraic multigrid (AMG) based preconditioner for the reduced Hessian of a linear‐quadratic optimization problem constrained by an elliptic partial differential equation. While the preconditioner generalizes a geometric multigrid preconditioner introduced in earlier works, its construction relies entirely on a standard AMG infrastructure built for solving the forward elliptic equation, thus allowing for it to be implemented using a variety of AMG methods and standard packages. Our analysis establishes a clear connection between the quality of the preconditioner and the AMG method used. The proposed strategy has a broad and robust applicability to problems with unstructured grids, complex geometry, and varying coefficients. The method is implemented using the Hypre package and several numerical examples are presented.more » « less
-
Abstract We propose a uniform block-diagonal preconditioner for condensed $$H$$(div)-conforming hybridizable discontinuous Galerkin schemes for parameter-dependent saddle point problems, including the generalized Stokes equations and the linear elasticity equations. An optimal preconditioner is obtained for the stiffness matrix on the global velocity/displacement space via the auxiliary space preconditioning technique (Xu (1994) The Auxiliary Space Method and Optimal Multigrid Preconditioning Techniques for Unstructured Grids, vol. 56. International GAMM-Workshop on Multi-level Methods (Meisdorf), pp. 215–235). A spectrally equivalent approximation to the Schur complement on the element-wise constant pressure space is also constructed, and an explicit computable exact inverse is obtained via the Woodbury matrix identity. Finally, the numerical results verify the robustness of our proposed preconditioner with respect to model parameters and mesh size.more » « less
-
Abstract This paper reports a recent development of the high-order spectral difference method with divergence cleaning (SDDC) for accurate simulations of both ideal and resistive magnetohydrodynamics (MHD) on curved unstructured grids consisting of high-order isoparametric quadrilateral elements. The divergence cleaning approach is based on the improved generalized Lagrange multiplier, which is thermodynamically consistent. The SDDC method can achieve an arbitrarily high order of accuracy in spatial discretization, as demonstrated in the test problems with smooth solutions. The high-order SDDC method combined with the artificial dissipation method can sharply capture shock interfaces with the oscillation-free property and resolve small-scale vortex structures and density fluctuations on relatively sparse grids. The robustness of the codes is demonstrated through long time simulations of ideal MHD problems with progressively interacting shock structures, resistive MHD problems with high Lundquist numbers, and viscous resistive MHD problems on complex curved domains.more » « less
