Search NASA⌕ Search

SEARCH · Search NASA

Results for “non-hydrostatic”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Convection-permitting hindcasting of diurnal variation of Mei-yu rainfall over East China with a global variable-resolution model

Two sets of 48-hr hindcast experiments for Mei-yu season using non-hydrostatic global variable-resolution model (MPAS-Atmosphere) with regional refinement at 4km resolution (V4km) and uniform 60km resolution (U60km) are first-time conducted to attempt at simulating the Mei-yu rainfall in 2015 over East China. The analysis focuses on the seasonal average of diurnal variation of precipitation amount (PA), intensity (PI), and frequency (PF). Both simulations reasonably reproduce the spatial distribution of PA. V4km is more skillful in simulating the spatial distributions and magnitudes of PF and PI. The diurnal cycle of Mei-yu rainfall shows a major early morning peak and a minor afternoon peak, contributed by precipitation during two sub-periods with distinct synoptic circulations. With strong Mei-yu in the first sub-period, the diurnal variation of PA is controlled by nocturnal southwesterly jet. V4km overestimates the morning peak mainly due to its bias in simulating boundary layer inertial oscillation. Although differences in PA between the simulations are small, U60km overestimates PF and underestimates PI. With weak Mei-yu in the second sub-period, the diurnal variation of PA is controlled by both synoptic circulation and local convection. At both resolutions, deviations in the large-scale circulation modulated by a few typhoons lead to positive biases in the morning peak of PA. After removing the typhoon impacts, V4km captures the observed diurnal cycle of PA well, while U60km significantly underestimates PA and PI particularly in the afternoon. Future studies focusing on advancing modeling of southwesterly jet and typhoons may further improve convection permitting simulation of Mei-yu rainfall.

Xu, Mingyue↗
DOE OSTI2021

Convection-Permitting Simulations With the E3SM Global Atmosphere Model

This paper describes the first implementation of the Δx = 3.25 km version of the Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40-day prescribed-sea-surface-temperature simulation (January 20 through February 28, 2020). This simulation was performed as part of the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) Phase 2 model intercomparison. Effective resolution is found to be ~ 6x the horizontal dynamics grid resolution despite using a coarser grid for physical parameterizations. Despite this new model being in an immature and untuned state, moving to 3.25 km grid spacing solves several long-standing problems with the E3SM model. In particular, Amazon precipitation is much more realistic, the frequency of light and heavy precipitation is improved, agreement between the simulated and observed diurnal cycle of tropical precipitation is excellent, and the vertical structure of tropical convection and coastal stratocumulus look good. In addition, the new model is able to capture the frequency and structure of important weather events (e.g., tropical cyclones, extratropical cyclones including atmospheric rivers, and cold air outbreaks). Interestingly, this model does not get rid of the erroneous southern branch of the intertropical convergence zone nor the tendency for strongest convection to occur over the Maritime Continent rather than the West Pacific, both of which are classic climate model biases. Several other problems with the simulation are identified, underscoring the fact that this model is a work in progress.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Bimodality in simulated precipitation frequency distributions and its relationship with convective parameterizations

Abstract Bimodality in precipitation frequency distributions is often evident in atmospheric models, but rarely in observations. This study i) proposes a metric to objectively quantify the bimodality in precipitation distributions, ii) evaluates model simulations contributed to the Coupled Model Intercomparison Project (CMIP) phase 5 (CMIP5), phase 6 (CMIP6), and the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) project by comparing them to satellite-based and reanalysis precipitation products, and iii) investigates possible origins of bimodal precipitation distributions. Our results reveal that about 83% (20 out of 24) of CMIP5 and 70% (21 out of 30) of CMIP6 models used in this study exhibit bimodal distributions. The few DYAMOND models that use a deep convective parameterization also show bimodal distributions, while most DYAMOND models do not. Predictably, the bimodality originates from the separation of precipitation process between resolved grid-scale and parameterized subgrid-scale. However, in a larger number of models bimodality arises from the parameterized subgrid-scale convective precipitation alone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Structural and mechanical properties of magnesium aluminate nanoceramics under high pressure

Nanoceramics may have different structural and physical properties compared to their coarse-grained counterparts. Here, we report the high-pressure study of micro- and nano-crystalline MgAl2O4 in order to examine the effect of particle size on the structural stability. A reversible pressure-induced phase transition (cubic to tetragonal) is observed in MgAl2O4 nanocrystals under non-hydrostatic pressure at room temperature, in contrast to the previously reported structural transition of MgAl2O4 at high pressure and high temperature. It is also found that the compressed MgAl2O4 microcrystals do not fracture further below 60 nm, suggesting a plastic deformation mechanism transition. MgAl2O4 with a grain size above ∼60 nm exhibits normal cracking behaviors, but shows metal-like plastic deformation behaviors below this critical size. It is implied that combined ductility and strength can be achieved in nanoceramic MgAl2O4.

Xu, Jianing (ORCID:0000000325474344)↗
DOE OSTI2021

The ultrahigh pressure stability of silver: An experimental and theoretical study

We measured the atomic volume of Ag in a toroidal diamond anvil cell to a maximum pressure of 416 GPa and calculated the atomic volume and elastic constants of Ag up to 750 and 460 GPa, respectively. Our density functional theory calculations at 0 K utilize an all-electron fully relativistic method and agree well with our volume measurements, particularly at pressures above ~75 GPa. In this work, we corrected our experimental results for non-hydrostaticity using a line shift analysis, and the resulting Vinet equation of state (EOS) parameters are reported. We find that the uniaxial stress sustained by Ag increases linearly up to 4.5 GPa at a pressure of 416 GPa. Our experimental results indicate that the fcc structure of Ag remains stable to at least 416 GPa at room temperature. Our theoretical results show that C44 increases as pressure increases and reaches a maximum at ~100 GPa above which it begins to decrease, a sign that the fcc structure of Ag is becoming unstable, and at V/V 0 = 0.30, the bcc structure is lower in energy than fcc.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shear strength measurements and hydrostatic compression of rhenium diboride under high pressures

Shear strength measurements have been carried out on rhenium diboride, ReB 2 , to a pressure of 74 GPa using a Radial X-ray Diffraction (R-XRD) technique in a diamond anvil cell using platinum as an internal x-ray pressure standard. The R-XRD result has provided a unique insight into the deformation of hexagonal crystal lattice under non-hydrostatic compression and variation of shear strength with increasing pressure. From R-XRD data, we have estimated hydrostatic component of compression to determine an equation of state of rhenium diboride yielding a bulk modulus of K 0 = 366 ± 25 GPa with a pressure derivative $K^{'}_{0}$ = 4.3 ± 0.5 in good agreement with hydrostatic density functional theory calculations. Here, the average lower bound of shear strength (τ) from various diffraction planes was then calculated using the measured interplanar d-spacing (d m ) and hydrostatic component of d-spacing (d p ) to be shown to approach 6.7 ± 0.4 GPa at 70 GPa. Our results show that the anisotropic compression effects observed in ReB 2 under hydrostatic compression are correlated to electronic structure changes under compression as predicted by theoretical calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Raman scattering of rhenium for secondary pressure calibration

With the increasing number of 100 s GPa experiments in the diamond anvil cell (DAC), improved accuracy in secondary pressure calibrations to extreme pressures is essential. The rhenium equation of state has been proposed as a pressure calibrant via x-ray diffraction with potentially broad applications as it is commonly used as a gasket material in DAC experiments. In this work, we conducted Raman spectroscopy experiments on rhenium in the DAC and report the pressure shift of the E2g mode, a refined high-pressure C44 and mode-Grüneisen parameter above 200 GPa. We used flat, beveled, and toroidal diamond anvils under quasi-hydrostatic and non-hydrostatic conditions. By measuring the E2g mode from the culet edge to the center, we analyzed pressure distribution based on culet type and distance from the anvil center. The shift in the E2g mode can be expressed as a function of pressure, and diamond edge measurements appear reliable across all anvil types. Comparing the center and edge pressures reveals anvil cupping, offering insights into predicting or preventing anvil failure during materials properties measurements at extreme conditions.

Diamond anvil cells↗
DOE OSTI2025

High pressure Raman spectroscopy of boron-rich boron carbides up to 50 GPa

We investigated the effects of increased boron content on the resistance to non-hydrostatic stress-induced (NHSI) local amorphous zones in boron-rich boron carbide (BxC) compounds. Using high-pressure Raman spectroscopy, we subjected B 4.3 C, B 6.4 C, and B 10.4 C to pressures up to 50 GPa and monitored their responses. Our results show that higher boron content delays the onset of NHSI local amorphous zone formation, shifting it from 35 GPa in B 4.3 C to approximately 50 GPa in B 10.4 C. This enhanced resistance is attributed to a reduction in the formation of the B 12 (CCC) polytype, which is susceptible to amorphization, and the greater flexibility of B–B–B chains. Furthermore, alternative mechanisms, such as boron vacancy-driven C–C bond formation, provide additional insights into defect-mediated structural changes that may influence the amorphization process. In conclusion, these findings highlight the dual role of boron content and defect mechanisms in improving the structural stability of B x C materials under extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sensitivity of organized convective storms to model grid spacing in current and future climates

Mesoscale convective systems (MCSs) are complexes of thunderstorms that become organized and cover hundreds of kilometres over several hours. MCSs are prolific rain producers in the tropics and mid-latitudes and are the major cause of warm-season flooding. Traditionally, climate models have difficulties in simulating MCSs partly due to the misrepresentation of complex process interactions that operate across a large range of scales. Significant improvements in simulating MCSs have been found in kilometre-scale models that explicitly simulate deep convection. However, these models operate in the grey zone of turbulent motion and have known deficiencies in simulating small-scale processes (e.g. entrainment, vertical mass transport). Here, we perform mid-latitude idealized ensemble MCS simulations under current and future climate conditions in three atmospheric regimes: hydrostatic (12 km horizontal grid spacing; Δx), non-hydrostatic (Δx = 4, 2 and 1 km) and large eddy scale (Δx = 500 m and 250 m). Our results show a dramatic improvement in simulating MCS precipitation, movement, cold pools, and cloud properties when transitioning from 12 km to 4 km Δx. Decreasing Δx beyond 4 km results in modest improvements except for up- and downdraft sizes, average vertical mass fluxes, and cloud top height and temperature, which continue to change. Most important for climate modelling is that Δx = 4 km simulations reliably capture most MCS climate change signals compared to those of the Δx = 250 m runs. Significantly different climate change signals are found in Δx = 12 km runs that overestimate extreme precipitation changes by up to 100%.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A protocol and analysis of year-long simulations of global storm-resolving models and beyond

We propose a protocol to evaluate and analyze year-long simulations of global storm-resolving models (GSRMs). The proposed protocol complements an earlier 40-day simulation protocol under the DYAMOND (DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains) project to allow the analysis of the seasonal cycle and associated climatic relevant phenomena. This intercomparison aims to reveal how GSRMs, which can simulate mesoscale convective systems (MCSs) in the global domain, reproduce atmospheric large-scale structures related to convection beyond month-long simulations. The intercomparison for one-year simulations is conducted by either atmosphere-only models or atmosphere–ocean coupled models with atmospheric horizontal mesh sizes less than 5 km. We recommend the continuous four seasons from March 2020 to February 2021 as a target period for the intercomparison but with options for many groups to join more flexibly. The output variables are collected at 0.25° resolution, and archives of a small set of native grid variables are encouraged to analyze tropical cyclones and MCSs. Through the proposed global storm-resolving simulation, we will evaluate the climatological distributions of the atmospheric large-scale circulations, such as the Intertropical Convergence Zone (ITCZ), monsoon, midlatitude jets, their time evolution, and the upscale impacts on them. We present sample analyses from a one-year simulation using the 3.5 km mesh Nonhydrostatic Icosahedral Atmospheric Model (NICAM), revealing the realistic zonal contrast of tropical precipitation, no double ITCZ structure, the reasonable midlatitude jet position and intensity but a weak bias of storm track activities, and a warm bias over the Eurasia during boreal winter. We also clarify the cross-scale interaction, such as the effects of cold pools on mean precipitation over the Maritime Continent through the precipitation diurnal cycle and the effects of resolved gravity waves on midlatitude mean flows. The proposed one-year simulation protocol is referred to as the “Sendai Protocol.” This protocol is not unique or definite for evaluating GSRMs; we prospect a hierarchical set of experiments from short-term to multi-year simulations as GSRM intercomparisons.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Object-Based Evaluation of Tropical Precipitation Systems in DYAMOND Simulations over the Maritime Continent

The present study uses an object-based evaluation metric to examine the precipitation bias over the Maritime Continent in the global cloud-resolving models. We specifically focus on the difference between the models that directly resolve convection and those using convection parameterization. The 40-day hindcast experiments of the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domain (DYAMOND) intercomparison project are evaluated against the high-resolution satellite rainfall products. The hindcast of the Central Weather Bureau Global Forecast System (CWBGFS) under the DYAMOND protocol is also included. The results indicate that most models simulate insufficient numbers of large precipitation system [object-based precipitation system (OPS), > 370 km in scale], indicating weaker convection organization. The observation indicates that the maximum precipitation within the OPS intensifies with increasing object size. All of the models capture this positive relationship, but most of them overestimate the sensitivity. Most of the models overestimate both the frequency and intensity of small OPS (< 160 km), except for the models with convection parameterization [i.e., CWBGFS, European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS)-9 km]. Although most of the models can reproduce the observed peak time of diurnal precipitation over the land area in the Maritime Continent, the simulated fractional contribution of different sizes of OPS to the total precipitation varies from model to model, and their peak times do not follow the observed ones with delayed peak times as the size of OPS increases from small, mid-size, to large categories. Most of the models reasonably capture the mean diurnal cycle peak time, but only the models with convection parameterization and Model for Prediction Across Scales (MPAS) can represent the diurnal evolution of fractional contribution from different OPSs. The implications of the current results to the upscale processes of the tropical convection systems in the global models are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Global Teleconnections: Atmospheric Interactions Across Scales [Slides]

This project aims at an improved understanding of non-hydrostatic dynamics in atmospheric modeling using the DOE’s next generation climate model. One aspect that we have been focusing on is baroclinic instability which is closely related to mid-latitude cyclone genesis. The other aspect that we recently developed is E3SM simulations under Radiative-Convective Equilibrium (RCE), a prototype of tropical climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Lagrangian Study of the Transition from Shallow to Deep Convection Using ASR Observations and Les Simulations (Final Report)

With the support from this grant, we have made the following accomplishments. We have found that a formula for fractional entrainment rate, a key aspect for parameterization moist convection, that works for both shallow and deep convection, which provide a promising direction for unified treatment of entrainment for both shallow and deep convection. We also examined the formulation of vertical velocity and showed that the effect of non-hydrostatic pressure and entrainment on vertical velocity are not constants, as assumed in past studies, but can be computed using a simple model for the cloudy updrafts. These results also apply to both deep and shallow convection, which opens the door for a more physically based parameterization of the updraft vertical velocity. We have also conducted extensive analysis of the RACORO aircraft data and have published the results in Wei et al. (2020). We have made extensive effort scrutinizing the data in detail and examining the various hypotheses with the aid of the LES simulations. This was the first application of the modified Paluch diagram on modern aircraft data and the first comprehensive analysis of mixing in clouds using real data (instead of only LES output) in a very long time. These analyses enhanced the science return of the ASR aircraft measurements. This grant also supported a study of the collisions among cold pools and an effort to build parameterizations of moist convection using tools developed in control theory.

42 ENGINEERING↗
DOE OSTI2021

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗
DOE OSTI2019

Experimental Evidence for Partially Dehydrogenated ε-FeOOH

Hydrogen in hydrous minerals becomes highly mobile as it approaches the geotherm of the lower mantle. Its diffusion and transportation behaviors under high pressure are important in order to understand the crystallographic properties of hydrous minerals. However, they are difficult to characterize due to the limit of weak X-ray signals from hydrogen. In this study, we measured the volume changes of hydrous ε-FeOOH under quasi-hydrostatic and non-hydrostatic conditions. Its equation of states was set as the cap line to compare with ε-FeOOH reheated and decompression from the higher pressure pyrite-FeO 2 Hx phase with 0 < x < 1. We found the volumes of those re-crystallized ε-FeOOH were generally 2.2% to 2.7% lower than fully hydrogenated ε-FeOOH. Our observations indicated that ε-FeOOH transformed from pyrite-FeO 2 Hx may inherit the hydrogen loss that occurred at the pyrite-phase. Hydrous minerals with partial dehydrogenation like ε-FeOOHx may bring it to a shallower depth (e.g., < 1700 km) of the lower mantle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Experimental and Computational Studies on Superhard Material Rhenium Diboride under Ultrahigh Pressures

An emerging class of superhard materials for extreme environment applications are compounds formed by heavy transition metals with light elements. In this work, ultrahigh pressure experiments on transition metal rhenium diboride (ReB 2 ) were carried out in a diamond anvil cell under isothermal and non-hydrostatic compression. Two independent high-pressure experiments were carried out on ReB 2 for the first time up to a pressure of 241 GPa (volume compression V/V 0 = 0.731 ± 0.004), with platinum as an internal pressure standard in X-ray diffraction studies. The hexagonal phase of ReB 2 was stable under highest pressure, and the anisotropy between the a-axis and c-axis compression increases with pressure to 241 GPa. The measured equation of state (EOS) above the yield stress of ReB 2 is well represented by the bulk modulus K 0 = 364 GPa and its first pressure derivative K 0 ´ = 3.53. Corresponding density-functional-theory (DFT) simulations of the EOS and elastic constants agreed well with the experimental data. DFT results indicated that ReB 2 becomes more ductile with enhanced tendency towards metallic bonding under compression. The DFT results also showed strong crystal anisotropy up to the maximum pressure under study. The pressure-enhanced electron density distribution along the Re and B bond direction renders the material highly incompressible along the c-axis. Our study helps to establish the fundamental basis for anisotropic compression of ReB 2 under ultrahigh pressures

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of the tangent linear and adjoint models of the global online chemical transport model MPAS-CO 2 v7.3

We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.

54 ENVIRONMENTAL SCIENCES↗