An official website of the United States government
Abstract We present a heuristic model to explain the suppression of deep convection in convection‐resolving models (CRMs) with a small number of grid columns, such as those used in super‐parameterized or multi‐scale modeling framework (MMF) general circulation models (GCM) of the atmosphere. Domains with few grid columns require greater instability to sustain convection because they force a large convective fraction, driving strong compensating subsidence warming. Updraft dilution, which is stronger for reduced horizontal grid spacing, enhances this effect. Thus, suppression of deep convection in CRMs with few grid columns can be reduced by increasing grid spacing. Radiative‐convective equilibrium simulations using standalone CRM simulations with the System for Atmospheric Modeling (SAM) and using GCM‐coupled CRM simulations with the Energy Exascale Earth System Model (E3SM)‐MMF confirm the heuristic model results.
more »
« less
Global System for Atmospheric Modeling: Model Description and Preliminary Results
Abstract The extension of a cloud‐resolving model, the System for Atmospheric Modeling (SAM), to global domains is described. The resulting global model, gSAM, is formulated on a latitude‐longitude grid. It uses an anelastic dynamical core with a single reference profile (as in SAM), but its governing equations differ somewhat from other anelastic models. For quasihydrostatic flows, they are isomorphic to the primitive equations (PE) in pressure coordinates but with the globally uniform reference pressure playing the role of actual pressure. As a result, gSAM can exactly maintain steady zonally symmetric baroclinic flows that have been specified in pressure coordinates, produces accurate simulations when initialized or nudged with global reanalyses, and has a natural energy conservation equation despite the drawbacks of using the anelastic system to model global scales. gSAM employs a novel treatment of topography using a type of immersed boundary method, the Quasi‐Solid Body Method, where the instantaneous flow velocity is forced to stagnate in grid cells inside a prescribed terrain. The results of several standard tests designed to evaluate the accuracy of global models with and without topography as well as results from real Earth simulations are presented.
more »
« less
- Award ID(s):
- 1743753
- PAR ID:
- 10445937
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 14
- Issue:
- 6
- ISSN:
- 1942-2466
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
A finite difference scheme is used to develop a numerical method to solve the flow of an unbounded viscoelastic fluid with zero to moderate inertia around a prolate spheroidal particle. The equations are written in prolate spheroidal coordinates, and the shape of the particle is exactly resolved as one of the coordinate surfaces representing the inner boundary of the computational domain. As the prolate spheroidal grid is naturally clustered near the particle surface, good resolution is obtained in the regions where the gradients of relevant flow variables are most significant. This coordinate system also allows large domain sizes with a reasonable number of mesh points to simulate unbounded fluid around a particle. Changing the aspect ratio of the inner computational boundary enables simulations of different particle shapes ranging from a sphere to a slender fiber. Numerical studies of the latter particle shape allow testing of slender body theories. The mass and momentum equations are solved with a Schur complement approach allowing us to solve the zero inertia case necessary to isolate the viscoelastic effects. The singularities associated with the coordinate system are overcome using L’Hopital’s rule. A straightforward imposition of conditions representing a time-varying combination of linear flows on the outer boundary allows us to study various flows with the same computational domain geometry. For the special but important case of zero fluid and particle inertia we obtain a novel formulation that satisfies the force- and torque-free constraint in an iteration-free manner. The numerical method is demonstrated for various flows of Newtonian and viscoelastic fluids around spheres and spheroids (including those with large aspect ratio). Good agreement is demonstrated with existing theoretical and numerical results.more » « less
-
Abstract This paper proposes and investigates the two-grid stabilized lowest equal-order finite element method for the time-independent dual-permeability-Stokes model with the Beavers-Joseph-Saffman-Jones interface conditions. This method is mainly based on the idea of combining the two-grid and the two local Gauss integrals for the dual-permeability-Stokes system. In this technique, we use a difference between a consistent mass matrix and an under-integrated mass matrix for the pressure variable of the dual-permeability-Stokes model using the lowest equal-order finite element quadruples. In the two-grid scheme, the global problem is solved using the standard$$ P_1-P_1-P_1-P_1 $$ finite element approximations only on a coarse grid with grid sizeH. Then, a coarse grid solution is applied on a fine grid of sizehto decouple the interface terms and the mass exchange terms for solving the three independent subproblems such as the Stokes equations, microfracture equations, and the matrix equations on the fine grid. On the other hand, microfracture and matrix equations are decoupled through the mass exchange terms. The weak formulation is reported, and the optimal error estimate is derived for the two-grid schemes. Furthermore, the numerical results validate that the two-grid stabilized lowest equal-order finite element method is effective and has the same accuracy as the coupling scheme when we choose$$ h=H^2 $$ .more » « less
-
SUMMARY Full-waveform inversion (FWI) methods rely on accurate numerical simulation of wave propagation in the analysed medium. Acoustic or elastic wave equations are often used to model seismic wave propagation. These types of simulations do not account for intrinsic attenuation effects due to material anelasticity, and thus correction techniques have been utilized in practice to partially compensate the anelasticity. These techniques often only consider the waveform amplitude correction based on averaging of overall amplitude response over the entire data set, and ignore the phase correction. Viscoelastic wave equations account for the anelastic response in both waveform amplitude and phase, and are therefore a more suitable alternative. In this study, we present a novel 3-D Gauss–Newton viscoelastic FWI (3-D GN-VFWI) method. To address the main challenge of the Gauss–Newton optimization, we develop formulas to compute the Jacobian efficiently by the convolution of virtual sources and backward wavefields. The virtual sources are obtained by directly differentiating the viscoelastic wave equations with respect to model parameters. In order to resolve complex 3-D structures with reasonable computational effort, a homogeneous attenuation (Q factor) is used throughout the analysis to model the anelastic effects. Synthetic and field experiments are performed to demonstrate the utility of the method. The synthetic results clearly demonstrate the ability of the method in characterizing a challenging velocity profile, including voids and reverse velocity layers. The field experimental results show that method successfully characterizes the complex substructure with two voids and undulating limestone bedrock, which are confirmed by invasive tests. Compared to 3-D elastic FWI results, the presented viscoelastic method produces more accurate results regarding depths of the voids and bedrock. This study suggests that the improvement of imaging accuracy would warrant the widespread use of viscoelastic wave equations in FWI problems. To our best knowledge, this is the first reported study on 3-D GN-VFWI at any scale. This study provides the new theory and formulation for the use of Gauss–Newton optimization on the 3-D viscoelastic problem.more » « less
-
Summary In this paper, a three‐dimensional numerical solver is developed for suspensions of rigid and soft particles and droplets in viscoelastic and elastoviscoplastic (EVP) fluids. The presented algorithm is designed to allow for the first time three‐dimensional simulations of inertial and turbulent EVP fluids with a large number particles and droplets. This is achieved by combining fast and highly scalable methods such as an FFT‐based pressure solver, with the evolution equation for non‐Newtonian (including EVP) stresses. In this flexible computational framework, the fluid can be modeled by either Oldroyd‐B, neo‐Hookean, FENE‐P, or Saramito EVP models, and the additional equations for the non‐Newtonian stresses are fully coupled with the flow. The rigid particles are discretized on a moving Lagrangian grid, whereas the flow equations are solved on a fixed Eulerian grid. The solid particles are represented by an immersed boundary method with a computationally efficient direct forcing method, allowing simulations of a large numbers of particles. The immersed boundary force is computed at the particle surface and then included in the momentum equations as a body force. The droplets and soft particles on the other hand are simulated in a fully Eulerian framework, the former with a level‐set method to capture the moving interface and the latter with an indicator function. The solver is first validated for various benchmark single‐phase and two‐phase EVP flow problems through comparison with data from the literature. Finally, we present new results on the dynamics of a buoyancy‐driven drop in an EVP fluid.more » « less
