Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeling 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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Decomposing a Compound Flood Event in an Urban Pacific Northwest Estuary: Primary Drivers and Projections for the Future

The Duwamish River Estuary (DRE) of Washington is prone to compound flooding during atmospheric river (AR) events. The processes contributing to such flooding (coastal and fluvial) have remained opaque to municipalities that are increasingly impacted. Here, we conduct a suite of coupled atmosphere-hydrology-ocean model simulations with varying forcing combinations (tide, surge, and/or river discharge) to identify the primary drivers of compound flooding during a recent AR event. We also test year 2,100 climate forcing to project how flooding and drivers may change in the future for the same event. We identify a clear distinction between dynamics in the downstream, engineered portion of the DRE compared to the upstream, “natural” river. Downstream, tides dominate water levels but contributions from storm surge and nonlinear tide-surge interaction elevate tide-only high waters from no flooding to major flooding. Upstream, total water levels during the event are ∼6 cm higher than downstream due to an increasing influence of river discharge over surge and tides. Notably, nonlinear surge-river and tide-river interactions act to reduce upstream water levels up to 50% compared to estimates which linearly sum tides, surge, and river, likely reducing flood vulnerability. Under two future climate scenarios: one with only sea level rise (SLR) and another with SLR plus atmospheric warming, we find little change in mechanism contributions to water levels. Expanded flooding in both cases is largely due to SLR, as a ∼50% increase to river discharge under the warming scenario has no impact downstream and marginally increases (∼3 cm) water level upstream.

atmospheric river↗
DOE OSTI2026

Recent Tropical Cyclone Outer Size Increases in the Western North Atlantic

The outer size of a tropical cyclone (TC) plays a pivotal role in its impact through a modulation of wind duration, rainfall extent and storm surge. While global and basin scale trends in TC size have been examined previously, whether there have been regional changes over time remains unclear. Here, using a suite of observations and reanalyses, we detect an increasing TC outer size over the western North Atlantic during 1979–2022. Over the western subtropical North Atlantic where the largest changes are found, the area-averaged increase in the radius of 17 ms -1 winds is about 7.5% over this period. Analysis of the large-scale environment reveals that the air-sea thermodynamic state has become more favorable for storm intensification, and thereby expansion, over that region. Further, numerical model simulations link changes in the environment with those in storm size, and suggest a combination of natural variability and external forcing in driving them.

Balaguru, Karthik [Pacific Northwest National Labo↗
DOE OSTI2025

Convergence of Cloud Droplet Spectral Relative Dispersion During Entrainment‐Mixing Based on Particle‐Resolved Direct Numerical Simulations

Entrainment-mixing processes critically impact cloud microphysical properties, but their effects on the relative dispersion (d) of cloud droplet size distributions (CDSDs) remain elusive. A direct numerical simulation model is initialized with different CDSDs to fill the gap. These results show that d decreases for broad CDSDs and increases for narrow ones, ultimately converging to approximately 0.5 regardless of initial CDSDs during the evaporation-dominated mixing stage. The supersaturation fluctuation and the shape of CDSDs jointly influence the convergence behavior of d. Further sensitivity tests show that the initial microphysical/dynamical/thermodynamical conditions exert negligible effects on the final converged value of d but affect the convergence rate (k). The k generally increases with increasing droplet number concentration and dissipation rate, and increases with decreasing liquid water content, relative humidity of entrained air, and mixing fraction of cloudy air. A conceptual model with two timescales is proposed; k and the timescales are negatively correlated, meaning that slow mixing and/or evaporation process results in slow convergence of d. In conclusion, this finding provides an important reference for improving understanding and parameterization of d during the entrainment-mixing processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Atmospheric Feedbacks Reverse the Sensitivity of Modeled Photosynthesis to Stomatal Function

Stomata mediate fluxes of carbon and water between terrestrial plants and the atmosphere. These fluxes are governed by stomatal function and can be modulated in many Earth system models by an empirical parameter within the calculation of stomatal conductance, the stomatal slope (𝑔 1⁢𝑀 ). Intuitively, 𝑔 1⁢𝑀 represents the marginal water cost of carbon, relating it to the emergent plant property of water use efficiency. Observations show that 𝑔 1⁢𝑀 can range widely across and within plant types in varying environments, and this distribution of 𝑔 1⁢𝑀 is not captured within Earth system models which represent each plant type with a single 𝑔 1⁢𝑀 value. Here we examine how 𝑔 1⁢𝑀 influences photosynthesis using coupled Earth system model simulations by perturbing 𝑔 1⁢𝑀 to observed 5⁢t⁢h and 95⁢t⁢h percentiles for each plant type. We find that high 𝑔 1⁢𝑀 reduces photosynthesis nearly everywhere, while low 𝑔 1⁢𝑀 has regionally dependent responses. Under fixed atmospheric conditions, low 𝑔 1⁢𝑀 increases photosynthesis in the Amazon and central North America but decreases photosynthesis in boreal Canada. These responses reverse when the atmosphere responds interactively due to spatially differing sensitivity to increases in temperature and vapor pressure deficit. Choice of 𝑔 1⁢𝑀 also influences photosynthetic response to changes in atmospheric carbon dioxide (CO 2 ), with lower and higher 𝑔1⁢𝑀 modifying total global response to elevated 2x preindustrial CO 2 by +6.4% and −9.6%, respectively. Our work demonstrates that atmospheric feedbacks are critical for determining the photosynthetic response to 𝑔 1⁢𝑀 assumptions and some regions are particularly sensitive to choice of 𝑔 1⁢𝑀 .

Liu, Amy X. [University of Washington, Seattle, WA↗
DOE OSTI2024

Learning protocols for the fast and efficient control of active matter

Exact analytic calculation shows that optimal control protocols for passive molecular systems often involve rapid variations and discontinuities. However, similar analytic baselines are not generally available for active-matter systems, because it is more difficult to treat active systems exactly. Here we use machine learning to derive efficient control protocols for active-matter systems, and find that they are characterized by sharp features similar to those seen in passive systems. We show that it is possible to learn protocols that effect fast and efficient state-to-state transformations in simulation models of active particles by encoding the protocol in the form of a neural network. We use evolutionary methods to identify protocols that take active particles from one steady state to another, as quickly as possible or with as little energy expended as possible. Our results show that protocols identified by a flexible neural-network ansatz, which allows the optimization of multiple control parameters and the emergence of sharp features, are more efficient than protocols derived recently by constrained analytical methods. Our learning scheme is straightforward to use in experiment, suggesting a way of designing protocols for the efficient manipulation of active matter in the laboratory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Future increase in compound soil drought-heat extremes exacerbated by vegetation greening

Compound soil drought and heat extremes are expected to occur more frequently with global warming, causing wide-ranging socio-ecological repercussions. Vegetation modulates air temperature and soil moisture through biophysical processes, thereby influencing the occurrence of such extremes. Global vegetation cover is broadly expected to increase under climate change, but it remains unclear whether vegetation greening will alleviate or aggravate future increases in compound soil drought-heat events. Here, using a suite of state-of-the-art model simulations, we show that the projected vegetation greening will increase the frequency of global compound soil drought-heat events, equivalent to 12-21% of the total increment at the end of 21st century. This increase is predominantly driven by reduced albedo and enhanced transpiration associated with increased leaf area. Although greening-induced transpiration enhancement has counteracting cooling and drying effects, the excessive water loss in the early growing season can lead to later soil moisture deficits, amplifying compound soil drought-heat extremes during the subsequent warm season. These changes are most pronounced in northern high latitudes and are dominated by the warming effect of CO 2 . Our study highlights the necessity of integrating vegetation biophysical effects into mitigation and adaptation strategies for addressing compound climate risks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Robust increase in South Asian monsoon rainfall under warming driven by extratropical clouds and ocean

The responses of South Asian Monsoon (SAM) circulation under global warming are known to be highly uncertain, leading to the wide spread of SAM rainfall projections among models. Here, we show that the uncertain SAM circulation in Coupled Model Intercomparison Project Phase 6 models consists of two robust components that partly offset each other: a weakening component linked to a global thermodynamic constraint and a northward shift component understood through a regional 2D energetic perspective. We further attribute the robust northward shift of SAM circulation to positive cloud feedback over the Eurasia Continent and heat uptake in the Southern Ocean. A set of climate model simulations supports the finding that SAM rainfall increase is primarily due to the northward shift of circulation driven by extratropical processes. This energetic perspective opens new avenues for predicting monsoon rainfall by connecting circulation changes to radiative forcing, feedbacks, and ocean heat uptake.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Projected changes in cross-equatorial northerly surges and their hydrological impacts in the near future

Recent decades have witnessed unprecedented extreme precipitation and catastrophic flooding in some of the most populated regions in Southeast Asia during boreal winter. These extreme events are often influenced by the cross-equatorial northerly surge (CENS), characterized by a strengthening of northerly moist monsoon winds south of the equator in the western Maritime Continent. However, the potential future changes in CENS and its hydrological impacts remain underexplored. Here, using an ensemble of high-resolution climate model simulations (CMIP6 HighResMIP), we show that the regional impacts of CENS on the intensity and frequency of extreme precipitation events are projected to significantly increase in the near future (2030–2050), particularly over the adjacent coastal regions of southern Indonesia and northwestern Australia. More specifically, the risk of CENS-related extreme precipitation is projected to increase by up to 39 ± 1.2% relative to the seasonal probability, despite no apparent changes in the CENS characteristics in the near future. Such stronger impact is attributed to the enhanced moistening efficiency and moist static instability due to a more humid and warmer environment, which leads to more intense CENS convection. Our results suggest the need for effective monitoring and disaster managements to respond to the increasing severity of such events in the near future.

Climate sciences↗
DOE OSTI2025

Cooling performance of direct expansion system using multiple stages and adjustable heat exchange areas for fresh air handling

Providing fresh air to rooms can ensure indoor air quality. The direct expansion fresh air handling units are widely used, and most of them are single-stage treatment units, which are less energy-efficient under variable operating conditions. To overcome the shortcomings of existing systems, a direct expansion system using multiple stages and adjustable heat exchange areas for fresh air handling is proposed in this study. The fresh air is handled by coils stage by stage to achieve high efficiency, and the dampers of the idle coils are opened when some of the coils are working, which results in the performance improvement under partial load conditions. The simulation model is established, and operating conditions at load ratio of 95.70%, 45.38%, 31.76% and 0 are selected as typical operating conditions. The proposed system performances under different conditions are compared with traditional system. The results show that: (1) For four conditions, as the load rate decreases, the proposed system can save energy by 6.2%, 26.1%, 46.5%, and 25.0% compared to the traditional system respectively; (2) The proposed system can bypass the filter when the outdoor air is good, which can further save energy by 1.3%, 1.8%, 3.4%, and 25.9% under four conditions, respectively.

Zha, Fuhai↗
DOE OSTI2025

Data Quality Assessment Process for Real-Time Data-Driven Traffic Microsimulation of Smart Corridor

Smart corridor digital twins are often created for the development and evaluation of emerging intelligent transportation systems and Connected and Autonomous Vehicle (CAV) technologies. However, limited guidance exists for data quality assessment for digital twin development. To address this, this paper discusses the data quality assessment utilized to develop data-driven real-time microscopic simulation models, i.e., digital twins, for two separate smart corridors: the North Avenue Smart Corridor in Atlanta, GA, and the Martin Luther King Smart Corridor in Chattanooga, Tennessee. This paper provides a summary of the author’s investigations of data requirements and data characteristics for the given smart corridor digital twin development efforts. With a focus on data, this summary includes a description of the data investigation process, key data issues observed, and strategies to address observed issues. Discussion is provided to help expand the lessons from these studies to other digital twin development efforts.

Saroj, Abhilasha [ORNL] (ORCID:0000000191178063)↗
DOE OSTI2026

Radiation-induced segregation in dilute Fe–Cr: A rate-theory framework for the Cr enrichment–depletion transition at the grain boundary

Radiation-induced segregation (RIS) poses a significant challenge for ferritic Fe–Cr alloys under irradiation, as it can compromise mechanical integrity and increase susceptibility to intergranular corrosion. Yet, the mechanisms governing Cr segregation remain incompletely understood. Here, in this study, we present a physics-based rate-theory model parameterized using self-consistent mean field theory-based Onsager transport coefficients to investigate RIS at the grain boundary in dilute Fe-(0.1 at. %) Cr. Under equal production rates of vacancies and self-interstitial atoms (SIA), and their equal absorption rates by bulk dislocations, the model simulates the experimentally observed transition from Cr enrichment at low temperatures to depletion at higher temperatures. Under these unbiased conditions, systematic investigation reveals that while temperature-dependent transport properties dictate the segregation direction, dose rate, grain size, and dislocation density only influence the magnitude and spatial extent of Cr segregation. However, under more realistic conditions of preferential vacancy production within damage cascade and/or preferential SIA absorption by bulk dislocations, the enrichment-to-depletion transition shifts to lower temperatures. Our findings demonstrate that RIS predictions based solely on transport coefficients are valid only under symmetric point defect flux conditions, and that biases in defect production and absorption must be considered for accurate predictions. This work provides a mechanistic framework for understanding RIS in ferritic alloys and informs alloy design for advanced nuclear systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Competition response of cloud supersaturation explains diminished Twomey effect for smoky aerosol in the tropical Atlantic

The Twomey effect brightens clouds by increasing aerosol concentrations, which activates more droplets and decreases cloud supersaturation in response to more competition for water vapor. To quantify this competition response, we used marine low cloud observations in clean and smoky conditions at Ascension Island in the tropical South Atlantic during the Layered Aerosol Smoke Interactions with Cloud (LASIC) campaign. These observations show similar increases in droplet number for increased accumulation-mode particles from surface-based and satellite cloud retrievals, demonstrating the importance of below-cloud aerosol measurements for retrieving aerosol–cloud interactions (ACI) in clean and smoky aerosol conditions. Four methods for estimating cloud supersaturation from aerosol–cloud measurements were compared, with cloud scene-based and parcel-based methods showing sufficient variability for a strong dependence on both aerosol accumulation number concentration and cloud-base updraft velocities. Decomposing aerosol-related changes in cloud albedo and optical depth shows the calculated competition response accounts for dampening the activation response by 12 to 35%, explaining the diminished Twomey effect at high aerosol concentrations observed for smoky conditions at LASIC and previously around the world. This result was consistent for independent supersaturation retrievals by cloud scene-based droplet number and cloud condensation nuclei and parcel-based multimode size-resolving Lagrangian methods. Translating aerosol effects to local radiative forcing with clean conditions as a proxy for preindustrial and smoky conditions for present-day showed that the competition response reduces cooling from the Twomey radiative forcing by 12 to 35%, providing an essential process-specific constraint for improving the representation of aerosol competition in climate model simulation of indirect aerosol forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Measurement of post-disruption runaway electron kinetic energy and pitch angle during final loss instability in DIII-D

Post-disruption runaway electron (RE) kinetic energy K and pitch angle sin$\vartheta$ are critical parameters for determining resulting first wall material damage during wall strikes, but are very challenging to measure experimentally. During the final loss instability, confined RE K and sin$\vartheta$ are reconstructed during center-post wall strikes for both high impurity (high-Z) and low impurity (low-Z) plasmas by combining soft x-ray, hard x-ray, synchrotron emission, and total radiated power measurements. Deconfined (wall impacting) RE sin$\vartheta$ is then reconstructed for these shots by using time-decay analysis of infra-red imaging. Additionally, deconfined RE K and sin$\vartheta$ are reconstructed for a low-Z downward loss shot by analyzing resulting damage to a sacrificial graphite dome limiter. The damage analysis uses multi-step modeling simulating plasma instability, RE loss orbits, energy deposition, and finally material expansion (MARS-F, KORC, GEANT-4, and finally COMSOL). Overall, mean kinetic energies are found to be in the range $\langle$$K$$\rangle$ $≈ 3-4$ MeV for confined REs. KORC simulations indicate that the final loss instability process does not change individual RE kinetic energy K. In conclusion, confined RE pitch angles are found to be fairly low initially pre-instability, $\langle$sin$\vartheta$$\rangle$ $≈ 0.1 – 0.2$, but appear to increase roughly 2x, to $\langle$sin$\vartheta$$\rangle$ $≈ 0.3 – 0.4$ for both confined and deconfined REs during instability onset in the low-Z case; this increase is not observed in the high-Z case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The effect of African Easterly Wave suppression by periodicity on Atlantic tropical cyclones

Research has shown that suppressing African Easterly Waves (AEWs) does not reduce basin-wide North Atlantic tropical cyclone (TC) frequency but can enhance TC environmental favorability. We investigated the AEW-TC relationship further by examining the effects of suppressing the two AEW periodicities individually on TC activity. Using regional model simulations, AEWs were prescribed or suppressed in the 2–6 d or 6–10 d ranges through the lateral boundary conditions. Seasonal TC frequency increased significantly when either AEW periodicity was suppressed, with a larger increase when the 2–6 d waves were suppressed. We also found that suppressing the 2–6 d waves increased mid-tropospheric moisture by up to 8%, as well as overall atmospheric instability, near the western coast of northern Africa. Furthermore, the convective disturbances that developed into TCs exhibited stronger rotation, increased ascending motion, and higher rainfall. Our results suggest that reduced 2–6 d AEW activity may trigger a more active TC season.

Danso, Derrick K↗
DOE OSTI2024

Measurements of the z > 5 Lyman-α forest flux autocorrelation functions from the extended XQR-30 data set

We present the first observational measurements of the Lyman-α (Ly α) forest flux autocorrelation functions in ten redshift bins from 5.1 ≤ z ≤ 6.0. We use a sample of 35 quasar sightlines at z > 5.7 from the extended XQR-30 data set; these data have signal-to-noise ratios of >20 per spectral pixel. We carefully account for systematic errors in continuum reconstruction, instrumentation, and contamination by damped Ly α systems. With these measurements, we introduce software tools to generate autocorrelation function measurements from any simulation. Our measurements of the smallest bin of the autocorrelation function increase with redshift when normalizing by the mean flux, $\langle{F}\rangle$. This increase may come from decreasing $\langle{F}\rangle$ or increasing mean free path of hydrogen-ionizing photons, λmfp. Recent work has shown that the autocorrelation function from simulations at z > 5 is sensitive to λmfp, a quantity that contains vital information on the ending of reionization. For an initial comparison, we show our autocorrelation measurements with simulation models for recently measured λmfp values and find good agreements. Further work in modelling and understanding the covariance matrices of the data is necessary to get robust measurements of λmfp from this data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Investigating event-shape methods in the search for the chiral magnetic effect in relativistic heavy ion collisions

The chiral magnetic effect (CME) is a phenomenon in which electric charge is separated by a strong magnetic field from local domains of chirality imbalance and parity violation in quantum chromodynamics. The CME-sensitive observable, the charge-dependent three-point azimuthal correlator Δ⁢𝛾 , is contaminated by a major physics background proportional to the particle's elliptic flow anisotropy 𝑣 2 . Event-shape engineering (ESE) binning events in dynamical fluctuations of 𝑣 2 and event-shape selection (ESS) binning events in statistical fluctuations of 𝑣 2 are two methods to search for the CME by projecting Δ⁢𝛾 to the measured anisotropy 𝑣 2 = 0 intercept. Here, we conduct a systematic study of these two methods using physics models as well as toy model simulations. It is observed that the ESE method fulfills the general premise of measuring the CME but is statistically hungry. It is found that the intercept from the ESS method depends on the details of the event content, such as the mixtures of background-contributing sources, because of statistical fluctuations of intertwining variables used in the method, and is thus not practically useful to measure the CME.

Relativistic heavy-ion collisions↗
DOE OSTI2024

Neural network emulation of flow in heavy-ion collisions at intermediate energies

Applications of new techniques in machine learning are speeding up progress in research in various fields. In this work, we construct and evaluate a deep neural network (DNN) to be used within a Bayesian statistical framework as a faster and more reliable alternative to the Gaussian process (GP) emulator of an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions at intermediate beam energies. We found strong evidence of the DNN being able to emulate the IBUU simulator's prediction on the strengths of protons' directed and elliptical flow very efficiently even with small training datasets and with accuracy about ten times higher than the GP. Here, limitations of our present work and future improvements are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗