Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic modeling”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Estimating the Importance of Viral Contributions to Soil Carbon Dynamics

ABSTRACT Biogeochemical models for predicting carbon dynamics increasingly include microbial processes, reflecting the importance of microorganisms in regulating the movement of carbon between soils and the atmosphere. Soil viruses can redirect carbon among various chemical pools, indicating a need for quantification and development soil carbon models that explicitly represent viral dynamics. In this opinion, we derive a global estimate of carbon potentially released from microbial biomass by viral infections in soils and synthesize a quantitative soil carbon budget from existing literature that explicitly includes viral impacts. We then adapt known mechanisms by which viruses influence carbon cycles in marine ecosystems into a soil‐explicit framework. Finally, we explore the diversity of virus–host interactions during infection and conceptualize how infection mode may impact soil carbon fate. Our synthesis highlights key knowledge gaps hindering the incorporation of viruses into soil carbon cycling research and generates specific hypotheses to test in the pursuit of better quantifying microbial dynamics that explain ecosystem‐scale carbon fluxes. The importance of identifying critical drivers behind soil carbon dynamics, including these elusive but likely pervasive viral mechanisms of carbon redistribution, becomes more pressing with climate change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A deep generative model for deciphering cellular dynamics and in silico drug discovery in complex diseases

Human diseases are characterized by intricate cellular dynamics. Single-cell transcriptomics provides critical insights, yet a persistent gap remains in computational tools for detailed disease progression analysis and targeted in silico drug interventions. Here we introduce UNAGI, a deep generative neural network tailored to analyse time-series single-cell transcriptomic data. This tool captures the complex cellular dynamics underlying disease progression, enhancing drug perturbation modelling and screening. When applied to a dataset from patients with idiopathic pulmonary fibrosis, UNAGI learns disease-informed cell embeddings that sharpen our understanding of disease progression, leading to the identification of potential therapeutic drug candidates. Validation using proteomics reveals the accuracy of UNAGI’s cellular dynamics analysis, and the use of the fibrotic cocktail-treated human precision-cut lung slices confirms UNAGI’s predictions that nifedipine, an antihypertensive drug, may have anti-fibrotic effects on human tissues. UNAGI’s versatility extends to other diseases, including COVID, demonstrating adaptability and confirming its broader applicability in decoding complex cellular dynamics beyond idiopathic pulmonary fibrosis, amplifying its use in the quest for therapeutic solutions across diverse pathological landscapes.

Neural Network↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2024

A variational phase-field framework for thermal softening and dynamic ductile fracture

A variational phase field model for dynamic ductile fracture is presented. Here, the model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.

42 ENGINEERING↗
DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗
DOE OSTI2025

Modeling and verification of dynamic field ionization for laser-target interactions

Integrating field ionization models into kinetic plasma simulations is required for a variety of applications, especially when field strengths vary from low to high regimes, such as in laser-target interactions. The introduction of new physics models into kinetic codes requires a rigorous verification of their accuracy through well-defined verification problems. In this work, the field ionization model that has been included in the kinetic plasma code VPIC is presented, along with the detailed approach adopted for its integration. This model includes a comprehensive range of field ionization mechanisms: multiphoton ionization, tunneling ionization, and barrier suppression ionization. New verification problems employed to evaluate the ionization model's fidelity are outlined, and the simulation parameters that affect the accuracy of simulation results are explored. Additionally, this work addresses the impact of field ionization on computational performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Subject-specific multi-scale modeling of the fate of inhaled aerosols

Determining the fate of inhaled aerosols in the respiratory system is essential in assessing the potential toxicity of inhaled airborne materials, responses to airborne pathogens, or in improving inhaled drug delivery. The availability of high-resolution clinical lung imaging and advances in the reconstruction of lung airways from CT images have led to the development of subject-specific in-silico 3D models of aerosol dosimetry, often referred to as computational fluid-particle-dynamics (CFPD) models. As CFPD models require extensive computing resources, they are typically confined to the upper and large airways. These models can be combined with lower-dimensional models to form multiscale models that predict the transport and deposition of inhaled aerosols in the entire respiratory tract. Understanding where aerosols deposit is only the first of potentially several key events necessary to predict an outcome, being a detrimental health effect or a therapeutic response. To that end, multiscale approaches that combine CFPD with physiologically-based pharmacokinetics (PBPK) models have been developed to evaluate the absorption, distribution, metabolism, and excretion (ADME) of toxic or medicinal chemicals in one or more compartments of the human body. CFPD models can also be combined with host cell dynamics (HCD) models to assess regional immune system responses. Here, this paper reviews the state of the art of these different multiscale approaches and discusses the potential role of personalized or subject-specific modeling in respiratory health.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Prioritizing urban heat adaptation infrastructure based on multiple outcomes: Comfort, health, and energy

Globally, cities face increasing extreme heat, impacting comfort, health, and energy consumption. Infrastructure-based heat adaptation strategies can improve these outcomes, but each strategy has a unique mix of benefits and drawbacks. Here, we apply an urbanized meteorological model (WRF) with the newly integrated multilayer BEP-Tree street tree model to dynamically downscale Earth System Model projections and a 3-D microclimate model (TUF-Pedestrian) to simulate the street-scale radiation environment impacting pedestrians. We evaluate the performance of five heat adaptation strategies (street trees, cool roofs, green roofs, rooftop photovoltaics (PV), and reflective pavements) during extreme heat events in three cities with contrasting background climates (Toronto, Phoenix, and Miami), under contemporary and end-of-century projected climates, based on three metrics: outdoor heat stress, air conditioning (AC) energy use, and ventilation of vehicular air pollution. No single adaptation strategy improves all three outcomes. While street trees inhibit ventilation, they reduce outdoor heat stress four times more effectively than the next best strategy via shade provision, fully offsetting heat stress increases under a high-emissions end-of-century climate scenario in all cities studied. Cool roofs and green roofs moderately reduce heat stress and energy use. Alternatively, rooftop PV with energy storage can generate sufficient power for space cooling but have marginal effects on heat stress. Reflective pavements are the least effective across metrics. Where the ventilation of street-level emissions is of less concern, our results clearly support the combination of street trees and rooftop PV as a highly complementary and effective means of adaptive mitigation across different climates and neighborhood densities.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Forward variable selection enables fast and accurate dynamic system identification with Karhunen-Loève decomposed Gaussian processes

A promising approach for scalable Gaussian processes (GPs) is the Karhunen-Loève (KL) decomposition, in which the GP kernel is represented by a set of basis functions which are the eigenfunctions of the kernel operator. Such decomposed kernels have the potential to be very fast, and do not depend on the selection of a reduced set of inducing points. However KL decompositions lead to high dimensionality, and variable selection thus becomes paramount. This paper reports a new method of forward variable selection, enabled by the ordered nature of the basis functions in the KL expansion of the Bayesian Smoothing Spline ANOVA kernel (BSS-ANOVA), coupled with fast Gibbs sampling in a fully Bayesian approach. It quickly and effectively limits the number of terms, yielding a method with competitive accuracies, training and inference times for tabular datasets of low feature set dimensionality. Theoretical computational complexities are O ( N P 2 ) in training and O ( P ) per point in inference, where N is the number of instances and P the number of expansion terms. The inference speed and accuracy makes the method especially useful for dynamic systems identification, by modeling the dynamics in the tangent space as a static problem, then integrating the learned dynamics using a high-order scheme. The methods are demonstrated on two dynamic datasets: a ‘Susceptible, Infected, Recovered’ (SIR) toy problem, along with the experimental ‘Cascaded Tanks’ benchmark dataset. Comparisons on the static prediction of time derivatives are made with a random forest (RF), a residual neural network (ResNet), and the Orthogonal Additive Kernel (OAK) inducing points scalable GP, while for the timeseries prediction comparisons are made with LSTM and GRU recurrent neural networks (RNNs) along with the SINDy package.

Hayes, Kyle↗
DOE OSTI2024

Topographical Controls on Hillslope-Scale Hydrology Drive Shrub Distributions on the Seward Peninsula, Alaska: Modeling Archive

The dataset is a model output used in a study that examined the role topography plays in determining shrub expansion by applying a coupled transect version of a mechanistic ecosystem model (ecosys) in a tundra hillslope site in the Seward Peninsula, Alaska. The study highlights the importance of mechanistically representing hillslope-scale coupled surface and sub-surface hydrology to accurately predict current plant distributions and future trajectories in Arctic ecosystems. The dataset contain ecosys model outputs used in a study: Article Title - Topographical Controls on Hillslope-Scale Hydrology Drive Shrub Distributions on the Seward Peninsula, Alaska. It includes the data files (ten *.csv files) used in creating figures in the publication and three *.pdf files (user guide and two Supplementary Information documents).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Data-scarce surrogate modeling of shock-induced pore collapse process

Understanding the mechanisms of shock-induced pore collapse is of great interest in various disciplines in sciences and engineering, including materials science, biological sciences, and geophysics. However, numerical modeling of the complex pore collapse processes can be costly. To this end, a strong need exists to develop surrogate models for generating economic predictions of pore collapse processes. Here, in this work, we study the use of a data-driven reduced-order model, namely dynamic mode decomposition, and a deep generative model, namely conditional generative adversarial networks, to resemble the numerical simulations of the pore collapse process at representative training shock pressures. Since the simulations are expensive, the training data are scarce, which makes training an accurate surrogate model challenging. To overcome the difficulties posed by the complex physics phenomena, we make several crucial treatments to the plain original form of the methods to increase the capability of approximating and predicting the dynamics. In particular, physics information is used as indicators or conditional inputs to guide the prediction. In realizing these methods, the training of each dynamic mode composition model takes only around 30 s on CPU. In contrast, training a generative adversarial network model takes 8 h on GPU. Moreover, using dynamic mode decomposition, the final-time relative error is around 0.3% in the reproductive cases. We also demonstrate the predictive power of the methods at unseen testing shock pressures, where the error ranges from 1.3 to 5% in the interpolatory cases and 8 to 9% in extrapolatory cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Autonomous alloy composition optimization using molecular dynamics guided by a large language model

Here, we present an autonomous materials discovery framework that couples a large language model (LLM) with molecular dynamics (MD) simulations to optimize Fe–Cr–Mn alloy compositions for tensile strength. Starting from six distinct compositions, the LLM operated as an intelligent agent, iteratively proposing changes based on prior simulation results and constraints. Over 50 iterations per case, the LLM adaptively explored the composition space, identifying high-strength regions, not easily accessible by conventional methods. The highest strength, 18.7 GPa, was achieved with Fe 71 Cr 25 Mn 4 composition, identified from a Fe 75 Cr 20 Mn 5 starting point. The LLM autonomously adjusted its strategy in real time, demonstrating closed-loop decision-making using commodity hardware. This approach showcases the potential of LLMs as scientific co-pilots, capable of accelerating materials discovery and generalizable to other domains like biology and drug design.

Autonomy↗
DOE OSTI2026

Divergent carbon use efficiency-growth rate tradeoff in popular biological growth models

Carbon use efficiency (CUE) is an important trait emerging from processes regulating biological growth. CUE can be computed either based on the growth of structural biomass or total biomass divided by substrate uptake rate. Nonequilibrium thermodynamics and observations suggest that, for an exponentially growing population of cells, structural biomass CUE should first increase, then peak, and finally decrease with specific growth rate; meanwhile, total biomass CUE increases asymptotically with specific growth rate. We compared predictions from six popular models that are often used for plant and microbial growth in existing ecosystem models. We found that, for an exponentially growing population of biological cells, (1) the source-driven Pirt and Compromise models predict that structural biomass CUE increase asymptotically with growth rate; (2) the apparent sink-driven modified Droop model predicts that structural biomass CUE decreases with growth rate; and (3) the sink-driven variable internal storage model and two dynamic energy budget models predict that structural biomass CUE first increases, then peaks, and finally decreases with growth rate. Moreover, the modified Droop model predicts that total biomass CUE is constant with growth rate, while all other five models predict that total biomass CUE increases with growth rate asymptotically. For non-exponential biological growth, we show that there is no static relationship between total biomass CUE or structural biomass CUE with respect to either growth rate or temperature. Therefore, we contend that biological growth models should explicitly represent interactions between substrate acquisition, substate transformation, and maintenance respiration to better capture observed CUE dynamics, and the sink-driven model should be preferred for general ecosystem biogeochemistry modeling.

Tang, Jinyun [Lawrence Berkeley National Laborator↗
DOE OSTI2026

Microwave-Assisted Plastic Upcycling: Dynamic Data Reconciliation, Parameter Estimation, and Kinetic Modeling

Microwave (MW)-assisted catalytic pyrolysis offers a promising pathway for efficient plastic upcycling. This work develops an integrated modeling framework combining dynamic data reconciliation, a temperature-dependent rate model, and a yield model to represent the time-varying production rate of components in MW-assisted LDPE pyrolysis conducted in a batch reactor. An Arrhenius-type rate model with a temperature-dependent reaction order is developed. A biexponential correlation is proposed for the yield of gaseous products that enables to capture the evolving product formation behavior during conversion. In the yield correlation, one term is used to represent the initial increase in yield, reflecting the rapid formation of intermediate or primary products at the early stages of the reaction when a larger fraction of the reactant remains available. As conversion progresses, the influence of this term gradually diminishes. The other term accounts for the subsequent decrease in the predicted yield, representing secondary reactions such as further cracking or coke formation that reduce the concentration of certain products at higher conversion. The model is found to accurately represent reconciled experimental flow rate profiles from an in-house MW-assisted catalytic batch reactor for major products, including ethylene, ethane, 1-butene, and benzene, across 250−350 °C. Ethylene remains the dominant product but decreases from about 41.95% at 250 °C to 30.14% at 350 °C, while heavier products increase significantly, with 1-butene rising to nearly 8.37% and benzene reaching 2.17% at intermediate temperatures. The model shows that the ethylene production rate can be maximized at around 270 °C. The models developed in this work can be utilized for process optimization, reactor design and scale-up of microwave-assisted plastic conversion technologies, and economic analysis.

Damahe, Harish [West Virginia Univ., Morgantown, W↗
DOE OSTI2026

An Integrated Modeling Framework for Sediment Dynamics During Urban Flooding: Application to Hurricane Harvey in Houston

Floodwater can mobilize and redistribute large volumes of sediment from upland to downstream urban areas, threatening infrastructure, water quality, and ecosystem health. However, existing modeling approaches often fail to capture sediment dynamics in urban floodplains due to the lack of integration between upland hydrological processes and riverine sediment transport. This study presents the first integrated modeling framework that couples the Energy Exascale Earth System Model (E3SM) land component, which simulates runoff and hillslope erosion, with TELEMAC-GAIA, a two-dimensional hydrodynamic and sediment transport model. This framework enables the fully distributed, process-based simulation of high-resolution (as fine as 30 m) sediment dynamics from hillslopes to floodplains. Applied to a highly urbanized watershed in Houston during Hurricane Harvey, this framework reproduced observed water levels at 16 USGS gauges (median R 2 = 0.83 and KGE = 0.78), key sediment dynamics such as sediment transport and deposition processes, and reproduced spatial deposition patterns consistent with LiDAR-derived data. Based on the simulation, we estimate 8.0 million m 3 of event-scale sediment deposition, including 5.7 million m 3 trapped in the flood-control reservoirs and 2.3 million m 3 deposited along major channels and floodplains. Using a representative unit removal cost, this corresponds to an estimated dredging cost of $581 million for total deposition. These results provide a first-order, physically based quantification of Harvey-scale sediment impacts. This study provides a valuable tool for the holistic analysis of sediment dynamics triggered by extreme urban flooding, supporting flood-resilience planning. More broadly, it highlights the importance of integrating physically based hydrological processes for urban flooding and sediment research.

Hurricane Harvey↗
DOE OSTI2024

Operation of Grid Forming Converters as Self Excited Induction Generators Under Non-Ideal Loading Conditions

Self-excited induction generators offer a robust solution for power production for standalone as well as grid-connected systems. In general, self-excited induction generators require excitation capacitors which make use of the machine magnetization characteristics for voltage build up process as well as operation at a specific frequency. In this paper, a self-excited induction machine is modeled with both the electrical and mechanical dynamics. This modeled virtual machine's dynamics are utilized for voltage build up process for a standalone photovoltaic converter connected to a local load for a microgrid application. The modeled machine's parameters are used from the name plate rating from the manufacturer. However, in a microgrid the accommodation of unbalanced and/or nonlinear harmonic rich load is a necessity, therefore, in this work the virtual self-excitation capacitors of the modeled machine are varied based on the machine characteristics. With the objective of ensuring harmonic free point of common coupling voltage, the modeled virtual self-excitation capacitors are varied to accomplish change in terminal frequency and the virtual load torque is varied to obtain voltage magnitude change. To verify the efficacy, the overall system is modeled in MATLAB/Simulink and PLECS domain and most important case studies are presented.

grid forming converters (GFM)↗
DOE OSTI2025

Ionization dynamics of intense laser-produced argon plasmas revealed by NLTE modeling

The ionization dynamics and transient behavior of under-dense plasma irradiated by an intense laser are investigated. We report two significant effects in the ionization behavior: (1) a surprisingly large delay in ionization response and (2) a stepwise ionization process which involves collisional and laser-driven photoionization (LDP) processes. Ionization induced by intense lasers can exhibit delayed responses due to rapid changes in conditions, particularly when atomic transition processes occur more slowly than the relevant time scales. Furthermore, modeling reveals that the two-step ionization process—collisional excitation followed by LDP—plays an important role in this ionization delay, with collisional excitation acting as the bottleneck. Even low-energy photons (∼3.5 eV) can predominantly ionize plasmas, challenging the conventional belief that such energies are insufficient to overcome the binding energy of bound electrons. These findings underscore the necessity of including such processes into plasma simulations for various laser-plasma experiments.

Collisional excitation↗
DOE OSTI2024

Open Source Modeling for Power Systems Dynamics

Time-domain simulations for power systems (both phasor and electromagnetic transients) are mostly dominated by commercial tools, such as PSS/E or PSCAD. However, with the increasing penetration of inverter-based resources and changes in model paradigm, new open-source tools are looking into become a possible alternative for a range of stakeholders. In this presentation we will focus on how to improve communication with industry and other stakeholders to properly showcase the value of open source simulators. Participants will identify current barriers facing the adoption of open-source time-domain tools, and will discuss how to improve the bench-marking processes across different tools, while maintaining software modularity to include new features and models as new challenges continue to happen in energy systems.

ENERGY PLANNING, POLICY, AND ECONOMY,MATHEMATICS A↗