Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical core”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effect of the Grell-Freitas Deep Convection Scheme in Quasi-Uniform and Variable-Resolution Aquaplanet CAM Simulations

The Grell-Freitas deep convection scheme (GF) has been implemented in the Community Atmosphere Model (CAM) with the nonhydrostatic Model for Prediction Across Scales (MPAS) dynamical core. Aquaplanet simulations have been performed with a global quasi-uniform resolution mesh with ~120 km grid spacing and a variable resolution mesh employing a circular refined region with ~30 km grid spacing centering at the equator and ~120 km elsewhere. GF produces overall patterns of the general circulations comparable to those of other aquaplanet experiments, including those produced using the default Zhang-McFarlane convection scheme (ZM) in CAM Version 5.4. GF alleviates the overestimated occurrence of weak intensity precipitation characteristic of ZM due to GF’s different closure methods and deep-convection trigger conditions. Compared to ZM, the GF characteristics of the equatorial Kelvin waves are closer to observations. GF also better simulates precipitation at variable resolution, with no unrealistic wave distortion by the refinement. Overall, GF shows less sensitivity of convective precipitation and heating/moistening tendencies to changes in horizontal resolution compared to ZM. However, vertically-overlapped total cloud fraction simulated with GF in the tropics is more sensitive to the resolution than ZM due to the larger resolution-sensitivity of high-level cloud fraction in the tropics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Global Atmosphere‐aerosol Model ICON‐A‐HAM2.3–Initial Model Evaluation and Effects of Radiation Balance Tuning on Aerosol Optical Thickness

Abstract The Hamburg Aerosol Module version 2.3 (HAM2.3) from the ECHAM6.3‐HAM2.3 global atmosphere‐aerosol model is coupled to the recently developed icosahedral nonhydrostatic ICON‐A (icon‐aes‐1.3.00) global atmosphere model to yield the new ICON‐A‐HAM2.3 atmosphere‐aerosol model. The ICON‐A and ECHAM6.3 host models use different dynamical cores, parameterizations of vertical mixing due to sub‐grid scale turbulence, and parameter settings for radiation balance tuning. Here, we study the role of the different host models for simulated aerosol optical thickness (AOT) and evaluate impacts of using HAM2.3 and the ECHAM6‐HAM2.3 two‐moment cloud microphysics scheme on several meteorological variables. Sensitivity runs show that a positive AOT bias over the subtropical oceans is remedied in ICON‐A‐HAM2.3 because of a different default setting of a parameter in the moist convection parameterization of the host models. The global mean AOT is biased low compared to MODIS satellite instrument retrievals in ICON‐A‐HAM2.3 and ECHAM6.3‐HAM2.3, but the bias is larger in ICON‐A‐HAM2.3 because negative AOT biases over the Amazon, the African rain forest, and the northern Indian Ocean are no longer compensated by high biases over the sub‐tropical oceans. ICON‐A‐HAM2.3 shows a moderate improvement with respect to AOT observations at AERONET sites. A multivariable bias score combining biases of several meteorological variables into a single number is larger in ICON‐A‐HAM2.3 compared to standard ICON‐A and standard ECHAM6.3. In the tropics, this multivariable bias is of similar magnitude in ICON‐A‐HAM2.3 and in ECHAM6.3‐HAM2.3. In the extra‐tropics, a smaller multivariable bias is found for ICON‐A‐HAM2.3 than for ECHAM6.3‐HAM2.3.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗
DOE OSTI2014

Co-evolution of galactic nuclei and globular cluster systems

We revisit the hypothesis that dense galactic nuclei are formed from inspiraling globular clusters. Recent advances in the understanding of the continuous formation of globular clusters over cosmic time and the concurrent evolution of the galaxy stellar distribution allow us to construct a simple model that matches the observed spatial and mass distributions of clusters in the Galaxy and the giant elliptical galaxy M87. In order to compare with observations, we model the effects of dynamical friction and dynamical evolution, including stellar mass loss, tidal stripping of stars, and tidal disruption of clusters by the growing galactic nucleus. We find that inspiraling globular clusters form a dense central structure, with mass and radius comparable to the typical values in observed nuclear star clusters (NSCs) in late-type and low-mass early-type galaxies. The density contrast associated with the NSC is less pronounced in giant elliptical galaxies. Our results indicate that the NSC mass as a fraction of mass of the galaxy stellar spheroid scales as M{sub NSC}/M{sub ∗}≈0.0025 M{sub ∗,11}{sup −0.5}. Thus disrupted globular clusters could contribute most of the mass of NSCs in galaxies with stellar mass below 10{sup 11} M {sub ☉}. The inner part of the accumulated cluster may seed the growth of a central black hole via stellar dynamical core collapse, thereby relieving the problem of how to form luminous quasars at high redshift. The seed black hole may reach ∼10{sup 5} M {sub ☉} within ≲ 1 Gyr of the beginning of globular cluster formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A new multi-dimensional general relativistic neutrino hydrodynamics code for core-collapse supernovae. IV. The neutrino signal

Considering six general relativistic, two-dimensional (2D) supernova (SN) explosion models of progenitor stars between 8.1 and 27 M {sub ☉}, we systematically analyze the properties of the neutrino emission from core collapse and bounce to the post-explosion phase. The models were computed with the VERTEX-COCONUT code, using three-flavor, energy-dependent neutrino transport in the ray-by-ray-plus approximation. Our results confirm the close similarity of the mean energies, (E), of ν-bar {sub e} and heavy-lepton neutrinos and even their crossing during the accretion phase for stars with M ≳ 10 M {sub ☉} as observed in previous 1D and 2D simulations with state-of-the-art neutrino transport. We establish a roughly linear scaling of 〈E{sub ν-bar{sub e}}〉 with the proto-neutron star (PNS) mass, which holds in time as well as for different progenitors. Convection inside the PNS affects the neutrino emission on the 10%-20% level, and accretion continuing beyond the onset of the explosion prevents the abrupt drop of the neutrino luminosities seen in artificially exploded 1D models. We demonstrate that a wavelet-based time-frequency analysis of SN neutrino signals in IceCube will offer sensitive diagnostics for the SN core dynamics up to at least ∼10 kpc distance. Strong, narrow-band signal modulations indicate quasi-periodic shock sloshing motions due to the standing accretion shock instability (SASI), and the frequency evolution of such 'SASI neutrino chirps' reveals shock expansion or contraction. The onset of the explosion is accompanied by a shift of the modulation frequency below 40-50 Hz, and post-explosion, episodic accretion downflows will be signaled by activity intervals stretching over an extended frequency range in the wavelet spectrogram.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Simple Cloud-Resolving E3SM Atmosphere Model Running on the Frontier Exascale System

We present an efficient and performance portable implementation of the Simple Cloud Resolving E3SM Atmosphere Model (SCREAM). SCREAM is a full featured atmospheric global circulation model with a nonhydrostatic dynamical core and state-of-the-art parameterizations for microphysics, moist turbulence and radiation. It has been written from scratch in C++ with the Kokkos library used to abstract the on-node execution model for both CPUs and GPUs. SCREAM is one of only a few global atmosphere models to be ported to GPUs. As far as we know, SCREAM is the first such model to run on both AMD GPUs and NVIDIA GPUs, as well as the first to run on nearly an entire Exascale system (Frontier). On Frontier, we obtained a record setting performance of 1.26 simulated years per day for a realistic cloud resolving simulation.

Taylor, Mark↗
DOE OSTI2025

RDycore/RDycore

A River dynamical core for E3SM

Johnson, Jeffrey N. [Cohere Consulting LLC]↗
DOE OSTI2023

Effects from Time Dependence of Ice Nucleus Activity for Contrasting Cloud Types

Abstract The role of time-dependent freezing of ice nucleating particles (INPs) is evaluated with the “Aerosol–Cloud” (AC) model in 1) deep convection observed over Oklahoma during the Midlatitude Continental Convective Cloud Experiment (MC3E), 2) orographic clouds observed over North California during the Atmospheric Radiation Measurement (ARM) Cloud Aerosol Precipitation Experiment (ACAPEX), and 3) supercooled, stratiform clouds over the United Kingdom, observed during the Aerosol Properties, Processes And Influences on the Earth’s climate (APPRAISE) campaign. AC uses the dynamical core of the WRF Model and has hybrid bin–bulk microphysics and a 3D mesoscale domain. AC is validated against coincident aircraft, ground-based, and satellite observations for all three cases. Filtered concentrations of ice (>0.1–0.2 mm) agree with those observed at all sampled levels. AC predicts the INP activity of various types of aerosol particles with an empirical parameterization (EP), which follows a singular approach (no time dependence). Here, the EP is modified to represent time-dependent INP activity by a purely empirical approach, using our published laboratory observations of time-dependent INP activity. In all simulated clouds, the inclusion of time dependence increases the predicted INP activity of mineral dust particles by 0.5–1 order of magnitude. However, there is little impact on the cloud glaciation because the total ice is mostly (80%–90%) from secondary ice production (SIP) at levels warmer than about −36°C. The Hallett–Mossop process and fragmentation in ice–ice collisions together initiate about 70% of the total ice, whereas fragmentation during both raindrop freezing and sublimation contributes <10%. Overall, total ice concentrations and SIP are unaffected by time-dependent INP activity. In the simulated APPRAISE case, the main causes of persistence of long-lived clouds and precipitation are predicted to be SIP in weak embedded convection and reactivation following recirculation of dust particles in supercooled layer cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Nonlinear Response of Atmospheric Blocking to Early Winter Barents–Kara Seas Warming: An Idealized Model Study

Wintertime Ural blocking (UB) has been shown to play an important role in cold extremes over Eurasia, and thus it is useful to investigate the impact of warming over the Barents–Kara Seas (BKS) on the behavior of Ural blocking. In our study the response of UB to stepwise tropospheric warming over the BKS is examined using a dry dynamic core model. Nonlinear responses are found in the frequency and local persistence of UB. The frequency and local persistence of the UB increase with the strength of BKS warming in a less strong range and decrease with the further increase of BKS warming, which is linked to the UB propagation influenced by upstream background atmospheric circulation. For a weak BKS warming, the UB becomes more persistent due to its less westward movement associated with intensified upstream zonal wind and meridional potential vorticity gradient (PVy) in the North Atlantic mid-high latitudes, which corresponds to a negative height response over the North Atlantic high latitudes. When BKS warming is strong, a positive height response appears in the early winter stratosphere, and its subsequent downward propagation leads to a negative NAO response or increased Greenland blocking events, which reduces zonal wind and PVy in the high latitudes from North Atlantic to Europe, thus enhancing the westward propagation of UB and reducing its local persistence. The transition to the negative NAO phase and the retrogression of UB are not found when numerically suppressing the downward influence of weakened stratospheric polar vortex, suggesting a crucial role of the stratospheric pathway in nonlinear responses of UB to the early winter BKS warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Physics-Dynamics Coupling in Weather, Climate, and Earth System Models: Challenges and Recent Progress

Geophysical models of the atmosphere and ocean invariably involve parameterizations. These represent two distinct areas: a) Subgrid processes which the model cannot (yet) resolve, due to its discrete resolution, and b) sources in the equation, due to radiation for example. Hence coupling between these physics parameterizations and the resolved fluid dynamics and also between the dynamics of the different fluids in the system (air and water) is necessary. This coupling is an important aspect of geophysical models. However, often model development is strictly segregated into either physics or dynamics. Hence, this area has many more unanswered questions than in-depth understanding. Furthermore, recent developments in the design of dynamical cores (significant increase of resolution, move to non-hydrostatic equation sets etc), extended process physics (prognostic micro physics, 3D turbulence, non-vertical radiation etc) and predicted future changes of the computational infrastructure (Exascale with its need for task parallelism, data locality and asynchronous time stepping for example) is adding even more complexity and new questions.

Gross, Markus↗
DOE OSTI2019

Physics–Dynamics Coupling with Element-Based High-Order Galerkin Methods: Quasi-Equal-Area Physics Grid

Atmospheric modeling with element-based high-order Galerkin methods presents a unique challenge to the conventional physics–dynamics coupling paradigm, due to the highly irregular distribution of nodes within an element and the distinct numerical characteristics of the Galerkin method. The conventional coupling procedure is to evaluate the physical parameterizations ( physics) on the dynamical core grid. Evaluating the physics at the nodal points exacerbates numerical noise from the Galerkin method, enabling and amplifying local extrema at element boundaries. Grid imprinting may be substantially reduced through the introduction of an entirely separate, approximately isotropic finite-volume grid for evaluating the physics forcing. Integration of the spectral basis over the control volumes provides an area-average state to the physics, which is more representative of the state in the vicinity of the nodal points rather than the nodal point itself and is more consistent with the notion of a “large-scale state” required by conventional physics packages. This study documents the implementation of a quasi-equal-area physics grid into NCAR’s Community Atmosphere Model Spectral Element and is shown to be effective at mitigating grid imprinting in the solution. The physics grid is also appropriate for coupling to other components within the Community Earth System Model, since the coupler requires component fluxes to be defined on a finite-volume grid, and one can be certain that the fluxes on the physics grid are, indeed, volume averaged.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Maximization of permanent trapping of CO{sub 2} and co-contaminants in the highest-porosity formations of the Rock Springs Uplift (Southwest Wyoming): experimentation and multi-scale modeling

Under this project, a multidisciplinary team of researchers at the University of Wyoming combined state-of-the-art experimental studies, numerical pore- and reservoir-scale modeling, and high performance computing to investigate trapping mechanisms relevant to geologic storage of mixed scCO{sub 2} in deep saline aquifers. The research included investigations in three fundamental areas: (i) the experimental determination of two-phase flow relative permeability functions, relative permeability hysteresis, and residual trapping under reservoir conditions for mixed scCO{sub 2}-­brine systems; (ii) improved understanding of permanent trapping mechanisms; (iii) scientifically correct, fine grid numerical simulations of CO{sub 2} storage in deep saline aquifers taking into account the underlying rock heterogeneity. The specific activities included: (1) Measurement of reservoir-­conditions drainage and imbibition relative permeabilities, irreducible brine and residual mixed scCO{sub 2} saturations, and relative permeability scanning curves (hysteresis) in rock samples from RSU; (2) Characterization of wettability through measurements of contact angles and interfacial tensions under reservoir conditions; (3) Development of physically-­based dynamic core-­scale pore network model; (4) Development of new, improved high-­performance modules for the UW-­team simulator to provide new capabilities to the existing model to include hysteresis in the relative permeability functions, geomechanical deformation and an equilibrium calculation (Both pore-­ and core-­scale models were rigorously validated against well-­characterized core-­ flooding experiments); and (5) An analysis of long term permanent trapping of mixed scCO{sub 2} through high-­resolution numerical experiments and analytical solutions. The analysis takes into account formation heterogeneity, capillary trapping, and relative permeability hysteresis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

A Flexible Atmospheric Modeling Framework for the CESM

We have created two global dynamical cores based on the unified system of equations and Z-grid staggering on an icosahedral grid, which are collectively called UZIM (Unified Z-grid Icosahedral Model). The z-coordinate version (UZIM-height) can be run in hydrostatic and nonhydrostatic modes. The sigma-coordinate version (UZIM-sigma) runs in only hydrostatic mode. The super-parameterization has been included as a physics option in both models. The UZIM versions with the super-parameterization are called SUZI. With SUZI-height, we have completed aquaplanet runs. With SUZI-sigma, we are making aquaplanet runs and realistic climate simulations. SUZI-sigma includes realistic topography and a SiB3 model to parameterize the land-surface processes.

icosahedral grid↗
DOE OSTI2015

Collaborative project. Ocean-atmosphere interaction from meso-to planetary-scale. Mechanisms, parameterization, and variability

This project aims to improve long term global climate simulations by resolving ocean mesoscale activity and the corresponding response in the atmosphere. The main computational objectives are; i) to perform and assess Community Earth System Model (CESM) simulations with the new Community Atmospheric Model (CAM) spectral element dynamical core; ii) use static mesh refinement to focus on oceanic fronts; iii) develop a new Earth System Modeling tool to investigate the atmospheric response to fronts by selectively filtering surface flux fields in the CESM coupler. The climate research objectives are 1) to improve the coupling of ocean fronts and the atmospheric boundary layer via investigations of dependency on model resolution and stability functions: 2) to understand and simulate the ensuing tropospheric response that has recently been documented in observations: and 3) to investigate the relationship of ocean frontal variability to low frequency climate variability and the accompanying storm tracks and extremes in high resolution simulations. This is a collaborative multi-institution project consisting of computational scientists, climate scientists and climate model developers. It specifically aims at DOE objectives of advancing simulation and predictive capability of climate models through improvements in resolution and physical process representation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Final Project Report

Most climate models are currently run with grid spacings of around 100km, which, with today’s computing power, allows for long (up to 1000 year) simulations, or ensembles of simulations to explore climate change and variability. However this grid spacing does not resolve important components of the weather/climate system such as atmospheric fronts and mesoscale systems, and ocean boundary currents and eddies. The overall aim of this project has been to look at the effect of these small-scale features on the weather/climate system using a suite of high and low resolution climate models, idealized models and observations. This project was only possible due to the highly scalable aspect of the CAM Spectral Element dynamical core, and the significant resources allocated at Yellowstone and NERSC for which we are grateful.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

MPAS Atmospheric Boundary Layer Simulation under Selected Stability Conditions: Evaluation Using the SWIFT Datasen

Modeling the transition from mesoscale to microscale is necessary in order to model different processes that affect a wind farm and to develop forecasting tools that operate at the farm scale. The mesoscale-to-microscale coupling (MMC) project is an A2e (Atmosphere-toelectrons) coordinated activity for developing modeling capabilities at the wind farm scale. By moving the focus of the research from a single wind turbine to the collection of turbines that comprise a wind farm, A2e extends the range of spatial and timescales that need representation in a model from tens of meters to hundreds of kilometers and timescales from a few seconds to days (Bokharaie et al. 2016). In the atmosphere, these scales are represented by mesoscale-tomicroscale models. The modeling available at these scales has differed in its representation of various physical processes. The MMC group is responsible for evaluating existing models at these scales and recommending a set of options for coupling the mesoscale and microscale with the best-performing models. The group was organized in 2015 and will explore options for coupling strategies with real-world test problems in fiscal year (FY) 2017. The model of choice for this exercise is WRF (Weather Research Forecasting) for mesoscale and WRF-LES (Large Eddy Simulation) for microscale simulations. The MPAS (Model Prediction Across Scales) variable mesh model that can be continuously refined; it has dynamic core and physics options adopted from WRF, which offer an alternative platform for modeling the mesoscale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Predicting Ice Sheet and Climate Evolution on Extreme Scales (PISCEES) Final Report

Predicting Ice Sheet and Climate Evolution at Extreme Scales (PISCEES) PISCEES is a SciDAC Earth System Modeling project with the following goals: To develop and apply robust, accurate, and scalable dynamical cores for ice sheet modeling on structured and unstructured meshes with adaptive refinements To evaluate ice sheet models using new tools and data sets for verification and validation (V&V) and uncertainty quantification (UQ) To integrate these models and tools into DOE's Accelerated Climate Model for Energy (ACME) using improved estimates of ice sheet initial conditions, we will simulate decade-to-century-scale evolution of the Greenland and Antarctic ice sheets, running PISCEES ice sheet models both in standalone mode and coupled to ACME. We aim to provide useful, credible predictions, including uncertainty ranges, of future ice-sheet mass loss and resulting changes in climate and sea level. PISCEES is jointly funded by the Office of Biological and Environmental Research (BER) and the Office of Advanced Scientific Computing Research (ASCR) of the DOE Office of Science. Participating institutions: Los Alamos National Laboratory (LANL) Lawrence Berkeley National Laboratory (LBNL) Oak Ridge National Laboratory (ORNL) Sandia National Laboratories (SNL) Florida State University (FSU) Massachusetts Institute of Technology (MIT) University of South Carolina (USC) University of Texas at Austin (UT) National Center for Atmospheric Research (NCAR) (https://www.scidac.gov/PISCEES/).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

A Non-hydrostatic Variable Resolution Atmospheric Model in ACME

The Accelerated Climate Modeling for Energy (ACME) Project is designed to perform coupled climate simulations at high resolution on current and next generation DOE leadership class supercomputers. Prior to this project, the ACME atmospheric model employed the primitive equations of motion, which constrained its maximum resolution and limited the physical phenomena that it could resolve. Our team removed this constraint by implementing a non-hydrostatic (NH) atmospheric model in the ACME Framework. Collaboration with SciDAC FASTMath institutes facilitated improvements to key technologies needed for the non-hydrostatic model. This model was tested in a variable resolution framework, which enables high resolution simulation regions at reasonable computational costs. The correct performance of the NH model was verified with standard dynamical-core tests, and by demonstrating that it reproduced expected climate statistics in the hydrostatic regime. The completed model was applied to evaluate the limitations of existing turbulent mountain stress and gravity wave drag parameterizations, in order to demonstrate the power of the new nonhydrostatic model to model new phenomena and enhance the quality of present day 1/4 degree simulations.

Hall, David Matthew↗