Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “planetary dynamics”

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

No evidence for a strong decrease of planetesimal accretion in old white dwarfs

ABSTRACT A large fraction of white dwarfs are accreting or have recently accreted rocky material from their planetary systems, thereby ‘polluting’ their atmospheres with elements heavier than helium. In recent years, the quest for mechanisms that can deliver planetesimals to the immediate vicinity of their central white dwarfs has stimulated a flurry of modelling efforts. The observed time evolution of the accretion rates of white dwarfs through their multi-Gyr lifetime is a crucial test for dynamical models of evolved planetary systems. Recent studies of cool white dwarf samples have identified a significant decrease of the mass accretion rates of cool, old white dwarfs over Gyr time-scales. Here, we revisit those results using updated white dwarf models and larger samples of old polluted H- and He-atmosphere white dwarfs. We find no compelling evidence for a strong decrease of their time-averaged mass accretion rates for cooling times between 1 and 8 Gyr. Over this period, the mass accretion rates decrease by no more than a factor of the order of 10, which is one order of magnitude smaller than the decay rate found in recent works. Our results require mechanisms that can efficiently and consistently deliver planetesimals inside the Roche radius of white dwarfs over at least 8 Gyr.

Blouin, Simon (ORCID:0000000296321436)↗
DOE OSTI2026

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗
DOE OSTI2022

Nature of the bonded-to-atomic transition in liquid silica to TPa pressures

First-principles calculations and analysis of the thermodynamic, structural, and electronic properties of liquid SiO 2 characterize the bonded-to-atomic transition at 0.1–1.6 TPa and 10 4 –10 5 K (1–9 eV), the high-energy-density regime relevant to understanding planetary interiors. We find strong ionic bonds that become short-lived due to high kinetics during the transition, with sensitivity of the transition temperature to pressure, and our calculated Hugoniots agree with past experimental data. Furthermore, these results reconcile previous experimental and theoretical findings by clarifying the nature of the bond dissociation process in early Earth and “rocky” (oxide) constituents of large planets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Exploring the Meteorological Impacts of Surface and Rooftop Heat Mitigation Strategies Over a Tropical City

Different heat mitigation technologies have been developed to improve the thermal environment in cities. However, the regional impacts of such technologies, especially in the context of a tropical city, remain unclear. The deployment of heat mitigation technologies at city-scale can change the radiation balance, advective flow, and energy balance between urban areas and the overlying atmosphere. We used the mesoscale Weather Research and Forecasting model coupled with a physically based single-layer urban canopy model to assess the impacts of five different heat mitigation technologies on surface energy balance, standard surface meteorological fields, and planetary boundary layer (PBL) dynamics for premonsoon typical hot summer days over a tropical coastal city in the month of April in 2018, 2019, and 2020. Results indicate that the regional impacts of cool materials (CMs), super-cool broadband radiative coolers, green roofs (GRs), vegetation fraction change, and a combination of CMs and GRs (i.e., “Cool city (CC)”) on the lower atmosphere are different at diurnal scale. Results showed that super-cool materials have the maximum potential of ambient temperature reduction of 1.6°C during peak hour (14:00 LT) compared to other technologies in the study. During the daytime hours, the PBL height was considerably lower than the reference scenario with no implementation of strategies by 700 m for super-cool materials and 500 m for both CMs and CC cases; however, the green roofing system underwent nominal changes over the urban area. During the nighttime hours, the PBL height increased by CMs and the CC strategies compared to the reference scenario, but minimal changes were evident for super-cool materials. The changes of temperature on the vertical profile of the heat mitigation implemented city reveal a stable PBL over the urban domain and a reduction of the vertical mixing associated with a pollution dome. This would lead to crossover phenomena above the PBL due to the decrease in vertical wind speed. Therefore, assessing the coupled regional impact of urban heat mitigation over the lower atmosphere at city-scale is urgent for sustainable urban planning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Observation of a highly conductive warm dense state of water with ultrafast pump–probe free-electron-laser measurements

The electrical conductivity of water under extreme temperatures and densities plays a central role in modeling planetary magnetic fields. Experimental data are vital to test theories of high-energy-density water and assess the possible development and presence of extraterrestrial life. These states are also important in biology and chemistry studies when specimens in water are confined and excited using ultrafast optical or free-electron lasers (FELs). Here we utilize femtosecond optical lasers to measure the transient reflection and transmission of ultrathin water sheet samples uniformly heated by a 13.6 nm FEL approaching a highly conducting state at electron temperatures exceeding 20 000 K. The experiment probes the trajectory of water through the high-energy-density phase space and provides insights into changes in the index of refraction, charge carrier densities, and AC electrical conductivity at optical frequencies. At excitation energy densities exceeding 10 MJ/kg, the index of refraction falls to n = 0.7, and the thermally excited free-carrier density reaches ne = 5 × 1027 m−3, which is over an order of magnitude higher than that of the electron carriers produced by direct photoionization. Significant specular reflection is observed owing to critical electron density shielding of electromagnetic waves. The measured optical conductivity reaches 2 × 104 S/m, a value that is one to two orders of magnitude lower than those of simple metals in a liquid state. At electron temperatures below 15 000 K, the experimental results agree well with the theoretical calculations using density-functional theory/molecular-dynamics simulations. With increasing temperature, the electron density increases and the system approaches a Fermi distribution. In this regime, the conductivities agree better with predictions from the Ziman theory of liquid metals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Global peak water limit of future groundwater withdrawals

Over the past 50 years, humans have extracted the Earth’s groundwater stocks at a steep rate, largely to fuel global agro-economic development. Given society’s growing reliance on groundwater, we explore ‘peak water limits’ to investigate whether, when and where humanity might reach peak groundwater extraction. Using an integrated global model of the coupled human–Earth system, we simulate groundwater withdrawals across 235 water basins under 900 future scenarios of global change over the twenty-first century. Here we find that global non-renewable groundwater withdrawals exhibit a distinct peak-and-decline signature, comparable to historical observations of other depletable resources (for example, minerals), in nearly all (98%) scenarios, peaking on average at 625 km 3 yr -1 around mid-century, followed by a decline through 2100. The peak and decline occur in about one-third (82) of basins, including 21 that may have already peaked, exposing about half (44%) of the global population to groundwater stress. Most of these basins are in countries with the highest current extraction rates, including the United States, Mexico, Pakistan, India, China, Saudi Arabia and Iran. Furthermore, these groundwater-dependent basins will probably face increasing costs of groundwater and food production, suggesting important implications for global agricultural trade and a diminished role for groundwater in meeting global water demands during the twenty-first century.

58 GEOSCIENCES↗
DOE OSTI2021

Investigation into the transport properties of planetary interiors through inelastic X-ray scattering experiments and quantum molecular dynamics (Final Technical Report)

This was a two-year research project for the period 08/15/2018 - 08/14/2020 conducted at the University of Nevada, Reno using Pronghorn, a newly built high-performance cluster located at the University. The short-term project goal was to provide computational support for the experimental campaigns at LCLS through data analysis of previous and upcoming scheduled experiments, along with performing state-of-the-art atomistic simulations of dense plasmas. The longer-term goal was to create a new computational high energy density physics group located at the University of Nevada, Reno. Guided by results from past and future high-resolution scattering experiments, we aimed to research and develop non-equilibrium and non-adiabatic atomistic simulations of dense plasmas that went beyond the Born-Oppenheimer approximation

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Scale‐Dependent Analysis of the Barotropic Vorticity Budget in a Global Ocean Simulation

Abstract The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak.

Khatri, Hemant↗
DOE OSTI2021

The Principal Hugoniot of Iron-Bearing Olivine to 1465 GPa

Shock compression experiments on natural compositions are imperative to accurately model planetary accretion and the interior dynamics of planets. Combining shock compression experiments from the Sandia Z Machine and the OMEGA EP laser facility with density functional theory-based molecular dynamics calculations, we report the first pressure-density-temperature (P-ρ-T) relationship of natural iron (Fe)-bearing olivine ((Mg 0.91 Fe 0.09 ) 2 SiO 4 ) on the principal Hugoniot between 166 and 1,465 GPa. Additionally, we report the first reflectivities of natural olivine liquid in this pressure range. Compared to the magnesium-endmember forsterite (Mg 2 SiO 4 ), the presence of Fe in typical mantle abundance (~9 wt% FeO) alters the U S -u P relation of olivine. Furthermore, the shock temperature and reflectivity of olivine are indistinguishable from forsterite where experimental conditions overlap. Both forsterite and olivine increase in reflectivity (and hence optical conductivity) with increasing temperature, with a maximum reflectivity of ~31% at shock velocities greater than 22 km/s (~800 GPa).

58 GEOSCIENCES↗
DOE OSTI2024

Thermoelastic Anomaly of Iron Carbonitride Across the Spin Transition and Implications for Planetary Cores

Carbon and nitrogen are considered as candidate light elements present in planetary cores. However, there is limited understanding regarding the structure and physical properties of Fe-C-N alloys under extreme conditions. Here diamond anvil cell experiments were conducted, revealing the stability of hexagonal-structured Fe 7 (N 0.75 C 0.25 ) 3 up to 120 GPa and 2100 K, without undergoing any structural transformation or dissociation. Notably, the thermal expansion coefficient and Grüneisen parameter of the alloy exhibit a collapse at 55–70 GPa. First-principles calculations suggest that such anomaly is associated with the spin transition of iron within Fe 7 (N 0.75 C 0.25 ) 3 . Our modeling indicates that the presence of ∼1.0 wt% carbon and nitrogen in liquid iron contributes to 9–12% of the density deficit of the Earth's outer core. The thermoelastic anomaly of the Fe-C-N alloy across the spin transition is likely to affect the density and seismic velocity profiles of (C,N)-rich planetary cores, thereby influencing the dynamics of such cores.

Geology↗
DOE OSTI2021

Ultrafast X-ray Diffraction Study of a Shock-Compressed Iron Meteorite above 100 GPa

Natural kamacite samples (Fe92.5Ni7.5) from a fragment of the Gibeon meteorite were studied as a proxy material for terrestrial cores to examine phase transition kinetics under shock compression for a range of different pressures up to 140 GPa. In situ time-resolved X-ray diffraction (XRD) data were collected of a body-centered cubic (bcc) kamacite section that transforms to the high-pressure hexagonal close-packed (hcp) phase with sub-nanosecond temporal resolution. The coarse-grained crystal of kamacite rapidly transformed to highly oriented crystallites of the hcp phase at maximum compression. The hcp phase persisted for as long as 9.5 ns following shock release. Comparing the c/a ratio with previous static and dynamic work on Fe and Fe-rich Fe-Ni alloys, it was found that some shots exhibit a larger than ideal c/a ratio, up to nearly 1.65. This work represents the first time-resolved laser shock compression structural study of a natural iron meteorite, relevant for understanding the dynamic material properties of metallic planetary bodies during impact events and Earth’s core elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The TESS-Keck Survey. II. An Ultra-short-period Rocky Planet and Its Siblings Transiting the Galactic Thick-disk Star TOI-561

We report the discovery of TOI-561, a multiplanet system in the galactic thick disk that contains a rocky, ultra-short-period planet. This bright (V = 10.2) star hosts three small transiting planets identified in photometry from the NASA TESS mission: TOI-561 b (TOI-561.02, P = 0.44 days, R {sub p} = 1.45 ± 0.11 R {sub ⊕}), c (TOI-561.01, P = 10.8 days, R {sub p} = 2.90 ± 0.13 R {sub ⊕}), and d (TOI-561.03, P = 16.3 days, R {sub p} = 2.32 ± 0.16 R {sub ⊕}). The star is chemically ([Fe/H] = −0.41 ± 0.05, [α/Fe] = +0.23 ± 0.05) and kinematically consistent with the galactic thick-disk population, making TOI-561 one of the oldest (10 ± 3 Gyr) and most metal-poor planetary systems discovered yet. We dynamically confirm planets b and c with radial velocities from the W. M. Keck Observatory High Resolution Echelle Spectrometer. Planet b has a mass and density of 3.2 ± 0.8 M {sub ⊕} and 5.5{sub −1.6}{sup +2.0} g cm{sup −3}, consistent with a rocky composition. Its lower-than-average density is consistent with an iron-poor composition, although an Earth-like iron-to-silicates ratio is not ruled out. Planet c is 7.0 ± 2.3 M {sub ⊕} and 1.6 ± 0.6 g cm{sup −3}, consistent with an interior rocky core overlaid with a low-mass volatile envelope. Several attributes of the photometry for planet d (which we did not detect dynamically) complicate the analysis, but we vet the planet with high-contrast imaging, ground-based photometric follow-up, and radial velocities. TOI-561 b is the first rocky world around a galactic thick-disk star confirmed with radial velocities and one of the best rocky planets for thermal emission studies.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Inelastic scattering in isotopologues of O 2 –Ar: the effects of mass, symmetry, and density of states

The two species considered here, O 2 (oxygen molecule) and Ar (argon-atom), are both abundant components of Earth's atmosphere and hence familiar collision partners in this medium. O 2 is quite reactive and extensively involved in atmospheric chemistry, including Chapman's cycle of the formation and destruction of ozone; while Ar, like N 2 , typically plays the nevertheless crucial role of inert collider. Inert species can provide stabilization to metastable encounter-complexes through the energy transfer associated with inelastic collisions. The interplay of collision frequency and energy transfer efficiency, with state lifetimes and species concentrations, contributes to the rich and varied chemistry and dynamics found in diverse environments ranging from planetary atmospheres to the interstellar and circumstellar media. The nature and density of bound and resonance states, coupled electronic states, symmetry, and nuclear spin-statistics can all play a role. Here, we systematically investigate some of those factors by looking at the O 2 –Ar system, comparing rigorous quantum-scattering calculations for the 16 O 16 O– 40 Ar, 18 O 16 O– 40 Ar, and 18 O 18 O– 40 Ar isotope combinations. A new accurate potential energy surface was constructed for this purpose holding the O 2 bond distance at its vibrationally averaged distance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Measuring the melting curve of iron at super-Earth core conditions

The discovery of more than 4500 extrasolar planets has created a need for modeling their interior structure and dynamics. Given the prominence of iron in planetary interiors, we require accurate and precise physical properties at extreme pressure and temperature. A first-order property of iron is its melting point, which is still debated for the conditions of Earth’s interior. Here, we used high-energy lasers at the National Ignition Facility and in situ x-ray diffraction to determine the melting point of iron up to 1000 gigapascals, three times the pressure of Earth’s inner core. We used this melting curve to determine the length of dynamo action during core solidification to the hexagonal close-packed (hcp) structure. We find that terrestrial exoplanets with four to six times Earth’s mass have the longest dynamos, which provide important shielding against cosmic radiation.

58 GEOSCIENCES↗
DOE OSTI2024

Role of atmospheric resonance and land–atmosphere feedbacks as a precursor to the June 2021 Pacific Northwest Heat Dome event

We demonstrate an indirect, rather than direct, role of quasi-resonant amplification of planetary waves in a summer weather extreme. We find that there was an interplay between a persistent, amplified large-scale atmospheric circulation state and soil moisture feedbacks as a precursor for the June 2021 Pacific Northwest “Heat Dome” event. An extended resonant planetary wave configuration prior to the event created an antecedent soil moisture deficit that amplified lower atmospheric warming through strong nonlinear soil moisture feedbacks, favoring this unprecedented heat event.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hydrodynamic instability at impact interfaces and planetary implications

Abstract Impact-induced mixing between bolide and target is fundamental to the geochemical evolution of a growing planet, yet aside from local mixing due to jetting – associated with large angles of incidence between impacting surfaces – mixing during planetary impacts is poorly understood. Here we describe a dynamic instability of the surface between impacting materials, showing that a region of mixing grows between two media having even minimal initial topography. This additional cause of impact-induced mixing is related to Richtmyer-Meshkov instability (RMI), and results from pressure perturbations amplified by shock-wave refraction through the corrugated interface between impactor and target. However, unlike RMI, this new impact-induced instability appears even if the bodies are made of the same material. Hydrocode simulations illustrate the growth of this mixing zone for planetary impacts, and predict results suitable for experimental validation in the laboratory. This form of impact mixing may be relevant to the formation of stony-iron and other meteorites.

58 GEOSCIENCES↗
DOE OSTI2021

Evolving CO 2 rather than SST leads to a factor of ten decrease in GCM convergence time

The high computational cost of Global Climate Models (GCMs) is a problem that limits their use in many areas. Recently an inverse climate modeling (InvCM) method, which fixes the global mean sea surface temperature (SST) and evolves the CO 2 mixing ratio to equilibrate climate, has been implemented in a cloud resolving model. In this paper, we apply InvCM to ExoCAM GCM aquaplanet simulations, allowing the SST pattern to evolve while maintaining a fixed global-mean SST. We find that InvCM produces the same climate as normal slab-ocean simulations but converges an order of magnitude faster. We then use InvCM to calculate the equilibrium CO 2 for SSTs ranging from 290 K to 340 K at 1 K intervals and reproduce the large increase in climate sensitivity at an SST of about 315 K at much higher temperature resolution. The speedup provided by InvCM could be used to equilibrate GCMs at higher spatial resolution or to perform broader parameter space exploration in order to gain new insight into the climate system. Additionally, InvCM could be used to find unstable and hidden climate states, and to find climate states close to bifurcations such as the runaway greenhouse transition.

54 ENVIRONMENTAL SCIENCES↗