Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Observing System Simulation”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Importance of different parameterization changes for the updated dust cycle modeling in the Community Atmosphere Model (version 6.1)

Abstract. The Community Atmosphere Model (CAM6.1), the atmospheric component of the Community Earth System Model (CESM; version 2.1), simulates the life cycle (emission, transport, and deposition) of mineral dust and its interactions with physio-chemical components to quantify the impacts of dust on climate and the Earth system. The accuracy of such quantifications relies on how well dust-related processes are represented in the model. Here we update the parameterizations for the dust module, including those on the dust emission scheme, the aerosol dry deposition scheme, the size distribution of transported dust, and the treatment of dust particle shape. Multiple simulations were undertaken to evaluate the model performance against diverse observations, and to understand how each update alters the modeled dust cycle and the simulated dust direct radiative effect. The model–observation comparisons suggest that substantially improved model representations of the dust cycle are achieved primarily through the new more physically-based dust emission scheme. In comparison, the other modifications induced small changes to the modeled dust cycle and model–observation comparisons, except the size distribution of dust in the coarse mode, which can be even more influential than that of replacing the dust emission scheme. We highlight which changes introduced here are important for which regions, shedding light on further dust model developments required for more accurately estimating interactions between dust and climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Quantum computing universal thermalization dynamics in a (2 + 1)D Lattice Gauge Theory

Simulating non-equilibrium phenomena in strongly-interacting quantum many-body systems, including thermalization, is a promising application of near-term and future quantum computation. By performing experiments on a digital quantum computer consisting of fully-connected optically-controlled trapped ions, we study the role of entanglement in the thermalization dynamics of a Z 2 lattice gauge theory in 2+1 spacetime dimensions. Using randomized-measurement protocols, we efficiently learn a classical approximation of non-equilibrium states that yields the gap-ratio distribution and the spectral form factor of the entanglement Hamiltonian. These observables exhibit universal early-time signals for quantum chaos, a prerequisite for thermalization. Our work, therefore, establishes quantum computers as robust tools for studying universal features of thermalization in complex many-body systems, including in gauge theories.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Achieving record hot spot energies with large HDC implosions on NIF in HYBRID-E

HYBRID-E is an inertial confinement fusion implosion design that increases energy coupled to the hot spot by increasing the capsule scale in cylindrical hohlraums while operating within the current experimental limits of the National Ignition Facility. HYBRID-E reduces the hohlraum scale at a fixed capsule size compared to previous HYBRID designs, thereby increasing the hohlraum efficiency and energy coupled to the capsule, and uses the cross-beam energy transfer (CBET) to control the implosion symmetry by operating the inner (23° and 30°) and outer (44° and 50°) laser beams at different wavelengths (Δλ> 0). Small case to capsule ratio designs can suffer from insufficient drive at the waist of the hohlraum. We show that only a small amount of wavelength separation between the inner and outer beams (Δλ1-2 Å) is required to control the symmetry in low-gas-filled hohlraums (0.3 mg/cm3 He) with enough drive at the waist of the hohlraum to symmetrically drive capsules 1180 μm in outer radius. This campaign is the first to use the CBET to control the symmetry in 0.3 mg/cm3 He-filled hohlraums, the lowest gas fill density yet fielded with Δλ> 0. We find a stronger sensitivity of hot spot P2 in μm per Angstrom (40–50 μm/Å wavelength separation) than observed in high-gas-filled hohlraums and previous longer pulse designs that used a hohlraum gas fill density of 0.6 mg/cm3. There is currently no indication of transfer roll-off with increasing Δλ, indicating that even longer pulses or larger capsules could be driven using the CBET in cylindrical hohlraums. We show that the radiation flux symmetry is well controlled during the foot of the pulse, and that the entire implosion can be tuned symmetrically in the presence of the CBET in this system, with low levels of laser backscatter out of the hohlraum and low levels of hot electron production from intense laser–plasma interactions. Radiation hydrodynamic simulations can accurately represent the early shock symmetry and be used as a design tool, but cannot predict the late-time radiation flux symmetry during the peak of the pulse, and semi-empirical models are used to design the experiments. Deuterium–tritium (DT)-layered tests of 1100 μm inner radius implosions showed performance close to expectations from simulations at velocities up to ~360 km/s, and record yields at this velocity, when increasing the DT fuel layer thickness to mitigate hydrodynamic mixing of the ablator into the hot spot as a result of defects in the ablator. However, when the implosion velocity was increased, mixing due to these defects impacted performance. The ratio of measured to simulated yield for these experiments was directly correlated with the level of observed mixing. These simulations suggest that reducing the mixing, e.g., by improving the capsule defects, could result in higher performance. In addition, future experiments are planned to reduce the coast time at this scale, delay between the peak compression and the end of the laser, to increase the hot spot convergence and pressure. To reduce the coast time by several hundred ps compared to the 1100 μm inner radius implosions, HYBRID-E has also fielded 1050 μm inner radius capsules, which resulted in higher hot spot pressure and a fusion energy yield of ~170 kJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Comparative Study of Cloud Properties Between Northern and Southern High Latitudes Based on ARM Observations and EAMv2 Simulations

Earth's high-latitude regions show a considerable amount of cloud coverage, which significantly affects global radiation budgets. Here, in this study, we examine ground- and ship-based remote sensing observations from two Atmospheric Radiation Measurement (ARM) field campaigns conducted in the high-latitude regions of Northern and Southern Hemispheres (NH and SH). A comparison of macrophysical properties for low, mid-level, high, and deep clouds shows that the cloud top heights and thickness values are higher in the Arctic for each cloud type compared with the Southern Ocean and Antarctic regions, possibly due to higher sea surface temperatures in the north Atlantic. The Arctic also shows lower ice phase occurrence frequencies compared with the SH high latitudes for low clouds. For high, mid-level and deep clouds, the cloud phase distributions as a function of temperature show small geographical variabilities. The observations are then compared with the Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) to evaluate the representations of these cloud properties. Here, the results show that the model captures cloud base and top heights well and reproduces the differences in cloud macrophysical properties between the two hemispheres. But the model underestimates ice phase frequencies for all cloud types. Such model cloud phase bias is strongly correlated with positive relative humidity biases and negative aerosol number concentration biases, indicating an important role of relative humidity and aerosols in modulating cloud phase partitioning in the model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Computational and experimental study of different brines in temperature swing solvent extraction desalination with amine solvents

Rapid global urbanization and high-salinity wastewater disposal from industrial activities have exerted significant pressure on water resources. Over the past few years, temperature swing solvent extraction (TSSE) has been identified as a promising technique to desalinate hypersaline brines. Despite its potential, the TSSE desalination literature has been mainly based on empirical insights, and the limited molecular simulation studies have primarily focused on NaCl brines. Herein, we use molecular dynamics (MD) simulations to study the TSSE desalination of four different brines, namely, KCl, KBr, NaCl, and NaBr using diisopropylamine as the solvent. Based on both bulk and interfacial brine-diisopropylamine MD results, here we investigate the qualitative and quantitative performance of the simulations by benchmarking these results against our experimental evaluations of these same systems. MD results provide satisfactory qualitative agreement with the experimental data of water solubilization in the organic phase and amine solubilization in the aqueous phase for the KBr, KCl, and NaBr brines. Also, the molecular mechanism of solvation of ionic species by water molecules over diisopropylamine suggested by the MD simulations is in agreement with our experimental data. However, larger qualitative and quantitative deviations were observed for the NaCl brines, and this is likely due to polarization and charge transfer effects, as quantified by our quantum chemical calculations.

42 ENGINEERING↗
DOE OSTI2023

Evaluation of distributed process-based hydrologic model performance using only a priori information to define model inputs

Fully distributed, integrated surface–subsurface hydrological models (ISSHMs) have seen renewed interest due to availability of better software, high performance computing facilities, and high-resolution, spatially extensive data products. ISSHMs are valuable as tools for advancing system understanding as they can resolve multiple processes defined on the plot scale including three-dimensional interaction of surface water and groundwater. Here, we evaluated the performance of an ISSHM, the Advanced Terrestrial Simulator (ATS), on seven diverse catchments across the continental US using widely available data products to define model inputs without calibration. We compare the ATS-simulated streamflow and evapotranspiration with gauge observations and MODIS-derived evapotranspiration, respectively. Using the Kling-Gupta Efficiency (KGE) as metric, ATS with default data products performed reasonably well at 6 of 7 catchments for streamflow. However, in one of those 6 catchments ATS had poor performance on baseflow and ATS’s overall performance was thus judged to be inadequate despite the acceptable KGE. ATS performance for evapotranspiration was good in all 7 catchments using default data products. In the two catchments where ATS streamflow performance using default data products was not acceptable, the performance was significantly improved by using local information on subsurface properties below the soil. We also compare the model-simulated streamflow and evapotranspiration with the Sacramento soil moisture accounting (SAC-SMA) model, a semi-distributed model that was calibrated on a catchment-by-catchment basis. Uncalibrated ATS performance is comparable to the calibrated SAC-SMA model in terms of streamflow while ATS performance is similar to or better (much better in certain catchments) in reproducing MODIS-derived evapotranspiration. Reasonably good performance of ATS without catchment-specific calibration provides new confidence in the ISSHM class of models and community data products as tools for advancing understanding of watershed function in a changing environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Rain-Induced Stratification of the Equatorial Indian Ocean and Its Potential Feedback to the Atmosphere

Surface freshening through precipitation can act to stably stratify the upper ocean, forming a rain layer (RL). RLs inhibit subsurface vertical mixing, isolating deeper ocean layers from the atmosphere. This process has been studied using observations and idealized simulations. The present ocean modeling study builds upon this body of work by incorporating spatially resolved and realistic atmospheric forcing. Fine-scale observations of the upper ocean collected during the Dynamics of the Madden-Julian Oscillation field campaign are used to verify the General Ocean Turbulence Model (GOTM). Spatiotemporal characteristics of equatorial Indian Ocean RLs are then investigated by forcing a 2D array of GOTM columns with realistic and well-resolved output from an existing regional atmospheric simulation. RL influence on the ocean-atmosphere system is evaluated through analysis of RL-induced modification to surface fluxes and sea surface temperature (SST). This analysis demonstrates that RLs cool the ocean surface on time scales longer than the associated precipitation event. A second simulation with identical atmospheric forcing to that in the first, but with rainfall set to zero, is performed to investigate the role of rain temperature and salinity stratification in maintaining cold SST anomalies within RLs. Approximately one third, or 0.1°C, of the SST reduction within RLs can be attributed to rain effects, while the remainder is attributed to changes in atmospheric temperature and humidity. The prolonged RL-induced SST anomalies enhance SST gradients that have been shown to favor the initiation of atmospheric convection. These findings encourage continued research of RL feedbacks to the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An Assessment of Nonhydrostatic and Hydrostatic Dynamical Cores at Seasonal Time Scales in the Energy Exascale Earth System Model (E3SM)

Abstract In global atmospheric modeling, the differences between nonhydrostatic (NH) and hydrostatic (H) dynamical cores are negligible in dry simulations when grid spacing is larger than 10 km. However, recent studies suggest that those differences can be significant at far coarser resolution when moisture is included. To better understand how NH and H differences manifest in global fields, we perform and analyze an ensemble of 28 and 13 km seasonal simulations with the NH and H dynamical cores in the Energy Exascale Earth System Model global atmosphere model, where the differences between H and NH configurations are minimized. A set of idealized rising bubble experiments is also conducted to further investigate the differences. Although NH and H differences are not significant in global statistics and zonal averages, significant differences in precipitation amount and patterns are observed in parts of the tropics. The most prominent differences emerge near India and the Western Pacific in the boreal summer, and the central‐southern Indian Ocean and Pacific in the boreal winter. Tropical differences influence surrounding regions through modification of the regional circulation and can propagate to the extratropics, leading to significant temperature and geopotential differences over the middle to high latitudes. While the dry bubble experiments show negligible deviation between H and NH dynamics until grid spacing is below 6.25 km, precipitation amount and vertical velocity are different in the moist case even at 25 km resolution.

58 GEOSCIENCES↗
DOE OSTI2023

Anthropogenic Aerosols Contribute to the Recent Decline in Precipitation Over the U.S. Southwest

Abstract The winter‐spring precipitation over the Southwestern United States (SWUS) decreased since 1980. It is frequently attributed to Pacific internal decadal variability, but recent studies found anthropogenic aerosols (AA) can also induce a transition to a negative Pacific Decadal Variability (PDV) phase. We revisit the attribution of SWUS drying by quantifying the contributions of anthropogenically forced decadal Pacific Sea Surface Temperatures (SSTs). Applying a low‐frequency component analysis to observations, Community Earth System Model version 2 (CESM2) all‐forcings and single‐forcing large ensembles, we find up to 42% of the observed precipitation trend to be related to the AA‐induced negative PDV‐like pattern, which is driven by the emission shift from the Western to the Eastern Hemisphere. In CESM2, other radiative forcings counteract the influence of AA, but it remains unclear whether the model correctly simulates this balance. This implies that the near‐future trajectories of these forcings, in particular Asian aerosols, are important for projections of SWUS precipitation.

Geology↗
DOE OSTI2025

The Role of Leaf Area Changes Within Plant CO 2 Physiological Impacts on the Global Hydrological Cycle

Rising atmospheric CO 2 concentrations enhance greenhouse warming and drive changes to plant physiology, leading to innumerable climate impacts. This study explores the impacts of plant responses on hydrological cycling at 2x preindustrial CO 2 concentrations by analyzing simulations that isolate plant physiological effects using the Community Earth System Model versions 1 and 2. We find that leaf area growth increases canopy evaporation, which offsets transpiration declines, and dampens changes in global mean evapotranspiration, precipitation, and runoff in a CESM2 experiment with dynamic leaf area. These leaf area impacts are also evident in the differences between CESM1 and CESM2, with CESM2 better capturing observed leaf area magnitudes but potentially overestimating leaf area-CO 2 sensitivity, highlighting the importance of plant CO 2 physiology on hydrological cycle changes and the need to improve its representation in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Optimizing the critical temperature and superfluid density of a metal-superconductor bilayer

A promising path to realizing higher superconducting transition temperatures 𝑇 c is the strategic engineering of artificial heterostructures. For example, quantum materials could, in principle, be coupled with other materials to produce a more robust superconducting state. Here, in this work, we add numerical support to the hypothesis that a strongly interacting superconductor weakened by phase fluctuations can boost its 𝑇 c by hybridizing the system with a metal. Using determinant quantum Monte Carlo, we simulate a two-dimensional bilayer composed of an attractive Hubbard model and a metallic layer in two regimes of the interaction strength −|𝑈|. In the strongly interacting regime, we find that increasing the interlayer hybridization 𝑡 ⊥ results in a nonmonotonic enhancement of 𝑇 c , with an optimal value comparable to the maximum 𝑇 c observed in the single-layer attractive Hubbard model, confirming trends inferred from other approaches. In the intermediate coupling regime, when −|𝑈| is close to the value associated with the maximum 𝑇 c of the single-layer model, increasing 𝑡 ⊥ tends to decrease 𝑇 c , implying that the correlated layer was already optimally tuned. Importantly, we demonstrate that the mechanism behind these trends is related to enhancement in the superfluid stiffness, as was initially proposed by Kivelson [Phys. B: Condens. Matter 318, 61 (2002)].

Zhang, Yutan [Univ. of California, Davis, CA (Unit↗
DOE OSTI2021

Energy-Resolved Information Scrambling in Energy-Space Lattices

Here, weakly interacting Fermi gases simulate spin lattices in energy space, offering a rich platform for investigating information spreading and spin coherence in a large many-body quantum system. We show that the collective spin vector can be determined as a function of energy from the measured spin density, enabling general energy-space resolved protocols. We measure an out-of-time-order correlation function in this system and observe the energy dependence of the many-body coherence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin↗
DOE OSTI2022

Tree hydrodynamic modelling of the soil–plant–atmosphere continuum using FETCH3

Modelling the water transport along the soil–plant–atmosphere continuum is fundamental to estimating and predicting transpiration fluxes. A Finite-difference Ecosystem-scale Tree Crown Hydrodynamics model (FETCH3) for the water fluxes across the soil–plant–atmosphere continuum is presented here. The model combines the water transport pathways into one vertical dimension, and assumes that the water flow through the soil, roots, and above-ground xylem can be approximated as flow in porous media. This results in a system of three partial differential equations, resembling the Richardson–Richards equation, describing the transport of water through the plant system and with additional terms representing sinks and sources for the transfer of water from the soil to the roots and from the leaves to the atmosphere. The numerical scheme, developed in Python 3, was tested against exact analytical solutions for steady state and transient conditions using simplified but realistic model parameterizations. The model was also used to simulate a previously published case study, where observed transpiration rates were available, to evaluate model performance. With the same model setup as the published case study, FETCH3 results were in agreement with observations. Through a rigorous coupling of soil, root xylem, and stem xylem, FETCH3 can account for variable water capacitance, while conserving mass and the continuity of the water potential between these three layers. FETCH3 provides a ready-to-use open access numerical model for the simulation of water fluxes across the soil–plant–atmosphere continuum.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Breakdown of superdiffusion in perturbed quantum integrable spin chains and ladders

Superdiffusive transport with dynamical exponent z=3/2 has been firmly established at finite temperature for a class of integrable systems with a non-Abelian global symmetry G. On the inclusion of integrability-breaking perturbations, diffusive transport with z=2 is generically expected to hold in the limit of late time. Recent studies of the classical Haldane-Ishimori-Skylanin model have found that perturbations that preserve the global symmetry lead to a much slower time scale for the onset of diffusion, albeit with uncertainty over the exact scaling exponent. That is, for perturbations of strength λ, the characteristic timescale for diffusion goes as t*∼λ−α for some α. Using large-scale matrix product state simulations, we investigate this behavior for perturbations to the canonical quantum model showing superdiffusion: the S=1/2 quantum Heisenberg chain. We consider a ladder configuration and look at various perturbations that either break or preserve the SU(2) symmetry, leading to scaling exponents consistent with those observed in one classical study [McCarthy , Phys. Rev. B 110, L180301 (2024)2469-995010.1103/PhysRevB.110.L180301]: α=2 for symmetry-breaking terms and α=6 for symmetry-preserving terms. We also consider perturbations from another integrable point of the ladder model with G=SU(4) and find consistent results. Finally, we consider a generalization to an SU(3) ladder and find that the α=6 scaling appears to be universal across superdiffusive systems when the perturbations preserve the non-Abelian symmetry G.

Wang, Kevin↗
DOE OSTI2023

An ab initio molecular dynamics study of varied compositions of the $\text{LiF-NaF-KF}$ molten salt

With increasing interest in molten salt reactors, there becomes a demand for investigations into thermophysical properties of salt systems. The LiF-NaF-KF (FLiNaK) salt system is a primary candidate for use in these reactors. However, the thermophysical properties of compositions outside the eutectic composition are still largely unknown. In this article, properties of ten unique compositions, including four ternary compositions, are investigated using ab initio molecular dynamics simulations across five temperatures between 900 K and 1300 K. The properties of interest are the density, thermal expansion, bulk modulus, compressibility, heat capacity, and enthalpy of mixing. In general, the results were found to be in good agreement with other literature and experimental results. The density and heat capacity had a tendency to be slightly underpredicted. No conclusions could be drawn about the bulk modulus and compressibility in terms of compositional dependence. The thermal expansion had a negative trend with respect to the LiF concentration and no trends were observed for the NaF or KF concentration. The enthalpy of mixing shows minima for the ternary compositions, with the near-equiatomic composition exhibiting the lowest values. Here this work shows the potential for compositional tailoring in the FLiNaK system to optimize thermophysical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dissociating the phononic, magnetic and electronic contributions to thermal conductivity: a computational study in alpha-iron

Computational tools to study thermodynamic properties of magnetic materials have, until recently, been limited to phenomenological modeling or to small domain sizes limiting our mechanistic understanding of thermal transport in ferromagnets. Herein, we study the interplay of phonon and magnetic spin contributions to the thermal conductivity in a-iron utilizing non-equilibrium molecular dynamics simulations. It was observed that the magnetic spin contribution to the total thermal conductivity exceeds lattice transport for temperatures up to two-thirds of the Curie temperature after which only strongly coupled magnon-phonon modes become active heat carriers. Characterizations of the phonon and magnon spectra give a detailed insight into the coupling between these heat carriers, and the temperature sensitivity of these coupled systems. Comparisons to both experiments and ab initio data support our inferred electronic thermal conductivity, supporting the coupled molecular dynamics/spin dynamics framework as a viable method to extend the predictive capability for magnetic material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular simulation of the structural and thermodynamic properties of n-alkane/brine interfacial systems with nonionic surfactants

All-atom molecular dynamics (MD) simulations are used to study the structural and thermodynamic properties of water + n-heptane interfacial systems in the presence of salinity (NaCl, NaI, KCl, and KI) and two different nonionic surfactants. Excellent qualitative and quantitative agreement with interfacial tension experimental data is obtained, but this requires applying a charge scaling factor to the ionic species, suggesting the likely role of polarizability effects. Tensoactive behavior is observed for the NaI and KI systems, coinciding with layering of the iodide species near the interface. Further, there are significant interaction differences between the surfactants and the different ionic species, but the structural and thermodynamic behavior of the surfactant molecules is relatively unaffected by the specific salt species. Our analysis suggest that the surfactant-surfactant interactions play a critical role in determining the interfacial behavior, with a much smaller impact associated with the composition of the brine phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗