Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeled”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A heuristic model inversion for coupled thermo-hydro-mechanical modelling of triaxial experiments

As multiphysics geomechanical models get developed, their increasing complexity and number of parameters make it particularly difficult to calibrate against experimental data. In this contribution, we present a heuristic workflow to invert for parameters of a coupled Thermo-Hydro-Mechanical (THM) model in a way that helps the theoretical modellers refine their definition of the underlying elasto-visco-plastic model itself. We apply this workflow to the calibration of deviatoric and volumetric data for two sets of triaxial experiments on mudstone and sandstone. Here, we show that beyond the calibration of well-defined parameters of the THM model, one or two scaling factors of the experimentally obtained yield surface and two more factors to capture the confinement and pore pressure dependency of the flow law lead to the satisfactory matching of series of experiments at varying confinements for two different rock types. Using this physical model, we also show that tracking the volumetric component of the mechanical power in the numerical simulations might allow reducing the number of experiments required to calibrate the model.

58 GEOSCIENCES↗
DOE OSTI2020

Non-local model for surface tension in fluid-fluid simulations

We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.

Howard, Amanda A.↗
DOE OSTI2021

Deep glassy state dynamic data challenge glass models: Configurational entropy models

This work was the result of the student's (Dongjie Chen) qualifying exam response. In it the experimental observation from our labs that the relaxation times of ultrastable glasses (a 20 million year old amber and a vapor deposited amorphous teflon) do not follow the expected Vogel-Fulcher or WLF types of divergence as the temperatures fall below the glass transition temperature. In the work we evaluated several theories by considering them in different ways. We used them as generally used and show that they do not capture the non-diverging behavior. Because the models were all entropy-based models, we also considered the possibility that the temperature dependence of the entropy used in the models is incorrect and chose two ways of estimating the entropy that do not follow the classical theoretical models. First we used the actually measured calorimetric response for the two ultrastable glasses, and we used the Milchev model, which does not go to zero above 0 Kelvin. Both of these approaches led to improved agreement between the theories and the experimental data. Finally, we also asked what the temperature dependence of the entropy would have to be so that the models agreed with the observed data. In some cases the entropy so-calculated was reasonable, while in others it was unphysical. In sum, the important point of the work is that the observation of non-diverging time-scales that suggest that there is no ideal glass transition can be reasonably captured by current theories if they take more realistic models of the entropy than is usually the case.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling information flow in a computer processor with a multi-stage queuing model

In this paper, we introduce a nonlinear stochastic model to describe the propagation of information inside a computer processor. In this model, a computational task is divided into stages, and information can flow from one stage to another. The model is formulated as a spatially-extended, continuous-time Markov chain where space represents different stages. This model is equivalent to a spatially-extended version of the M/M/s queue. The main modeling feature is the throttling function which describes the processor slowdown when the amount of information falls below a certain threshold. We derive the stationary distribution for this stochastic model and develop a closure for a deterministic ODE system that approximates the evolution of the mean and variance of the stochastic model. In conclusion, we demonstrate the validity of the closure with numerical simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Load‐Balancing Intense Physics Calculations to Embed Regionalized High‐Resolution Cloud Resolving Models in the E3SM and CESM Climate Models

Abstract We design a new strategy to load‐balance high‐intensity sub‐grid atmospheric physics calculations restricted to a small fraction of a global climate simulation's domain. We show why the current parallel load balancing infrastructure of Community Earth System Model (CESM) and Energy Exascale Earth Model (E3SM) cannot efficiently handle this scenario at large core counts. As an example, we study an unusual configuration of the E3SM Multiscale Modeling Framework (MMF) that embeds a binary mixture of two separate cloud‐resolving model grid structures that is attractive for low cloud feedback studies. Less than a third of the planet uses high‐resolution (MMF‐HR; sub‐km horizontal grid spacing) relative to standard low‐resolution (MMF‐LR) cloud superparameterization elsewhere. To enable MMF runs with Multi‐Domain cloud resolving models (CRMs), our load balancing theory predicts the most efficient computational scale as a function of the high‐intensity work's relative overhead and its fractional coverage. The scheme successfully maximizes model throughput and minimizes model cost relative to precursor infrastructure, effectively by devoting the vast majority of the processor pool to operate on the few high‐intensity (and rate‐limiting) high‐resolution (HR) grid columns. Two examples prove the concept, showing that minor artifacts can be introduced near the HR/low‐resolution CRM grid transition boundary on idealized aquaplanets, but are minimal in operationally relevant real‐geography settings. As intended, within the high (low) resolution area, our Multi‐Domain CRM simulations exhibit cloud fraction and shortwave reflection convergent to standard baseline tests that use globally homogenous MMF‐LR and MMF‐HR. We suggest this approach can open up a range of creative multi‐resolution climate experiments without requiring unduly large allocations of computational resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Improved Boundary Conditions for Coupled Geospace Models: An Application in Modeling Spacecraft Surface Charging Environment

Abstract Spacecraft surface charging in the inner magnetosphere often occurs in the pre‐midnight to the dawn sector when electron fluxes of tens of keV increase. Inner magnetosphere ring current models can be used to simulate/predict the spacecraft surface charging environment, with their outer boundary conditions specified either based on observations or provided by other models, such as MHD models. In the latter approach, using MHD quantities, the flux spectrum at the outer boundary is commonly assumed to follow a Kappa or Maxwellian distribution function. Such a method however often departs greatly from the realistic spectrum at E < tens of keV, a crucial energy range in the surface charging anomaly. In order to achieve a better representation of the surface charging environment, we propose to combine the MHD‐parameterized flux spectrum with an empirical electron flux model of E < 40 keV to set the electron flux boundary condition. Results indicate that as opposed to the case where the MHD‐parameterized flux distribution is solely used at the model boundary, simulations with the new boundary condition yields a more intense surface charging environment. The integrated electron flux between 10 < E < 50 keV, a measure of the severity of the surface charging environment, is significantly enhanced by 1‐2 orders of magnitude, leading to a much better agreement with Van Allen Probes measurements. This study hence demonstrates a reasonable solution to the setting of outer boundary conditions for inner magnetosphere models and is recommended for coupled geospace circulation models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Incorporating an Interactive Fire Plume-Rise Model in the DOE's Energy Exascale Earth System Model Version 1 (E3SMv1) and Examining Aerosol Radiative Effect

The vertical distribution of biomass burning aerosol (BBA) is important in regulating their impacts on weather and climate. The plume-rise process affects the injection height of BBA and interacts with the air parcel lifting and cloud processes. However, these processes are not represented in most global climate models. In this study, we replaced the fixed vertical profiles of monthly BBA emissions in the Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) with an interactive fire plume-rise model. The vertical distribution of BBA emissions was calculated as a function of ambient thermodynamic conditions from the host E3SMv1, with distributions of fire sizes and sensible heat fluxes derived from the observations. The maximum fire radiative power (FRP) technique was used to determine the fire size. Scaling-FRP technique is used to calculate the wildfire heat release. Daily BBA emission, superimposed with a fire diurnal cycle retrieved from the satellite observation, was included in model simulations. The model shows improved agreement with satellite retrievals and in situ measurement during the National Oceanic and Atmospheric Administration Wildfire Experiment for Cloud chemistry, Aerosol absorption, and Nitrogen campaign. The model-observation comparison demonstrates the importance of the plume-rise model and fire diurnal cycle assumption in determining the BBA fields. We also find that E3SMv1 with new features produces a larger carbonaceous aerosol burden, leading to 0.13 W m –2 warming at the top of atmosphere compared to the default E3SMv1. This highlights the importance of accurately representing the BBA injection height and indicates a no-linear nature in the BBA-induced radiative effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗
DOE OSTI2023

Towards fast, accurate predictions of RF simulations via data-driven modeling: Forward and lateral models

Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modeling and real-time control applications. More accurate simulations with fast electron diffusion are even slower, requiring multiple hours of run time with parallel processing. The machine learning models use a database of 16,000+ GEN-RAY/CQL3D simulations for training, validation, and testing. Latin hypercube sampling methods implemented in πScope ensure that the database covers the range of 9 input parameters (n e0 , T e0 , I p , B t , R 0 , n ∥︀ , Z e f f , V loop , P LHCD ) with sufficient density in all regions of parameter space. The surrogate models reduce the computation time from minutes-hours to ms with high accuracy across the input parameter space. Data-driven surrogate models also allow for solving inverse and “lateral” problems. A surrogate model for the inverse problem maps from a desired current drive or power deposition profile to a set of input parameters that would result in such a profile, while a surrogate model for the lateral problem maps from a measured experimental quantity such as hard x-ray emission to a current drive or power deposition profile. In conclusion, the πScope database creation workflow is flexible and applicable to other RF simulation codes such as TORIC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Wind plant wake losses: Disconnect between turbine actuation and control of plant wakes with engineering wake models

Wake losses from neighboring plants may become a major factor in wind plant design and control as additional plants are constructed in areas with high wind resource availability. Because plant wakes span a large range of physical scales, from turbine rotor diameter to tens of kilometers, it is unclear whether conventional wake models or turbine control strategies are effective at the plant scale. Wake steering and axial induction control are evaluated in the current work as means of reducing the impact of neighboring wind plants on power and levelized cost of electricity. FLOw Redirection and Induction in Steady State (FLORIS) simulations were performed with the Gauss–Curl Hybrid and TurbOPark wake models as well as two operation and maintenance models to investigate control setpoint sensitivity to wake representation and economic factors. Both wake models estimate losses across a range of atmospheric conditions, although the wake loss magnitude is dependent on the wake model. Annual energy production and levelized cost of electricity are driven by wind direction frequency, with frequently aligned plants experiencing the greatest losses. However, both wake steering and axial induction are unable to mitigate the impact of upstream plants. Wake steering is constrained by plant geometry, since wake displacement is much less than the plant wake width, while axial induction requires curtailing the majority of turbines in upstream plants. Individual turbine strategies are limited by their effective scale and model representation. New wake models that include plant-scale physics are needed to facilitate the design of effective plant wake control strategies.

Scott, Ryan (ORCID:0000000328107574)↗
DOE OSTI2024

Uncertainty quantification of mass models using ensemble Bayesian model averaging

Developments in the description of the masses of atomic nuclei have led to various nuclear mass models that provide predictions for masses across the whole chart of nuclides. These mass models play an important role in understanding the synthesis of heavy elements in the rapid neutron capture ( r ) process. However, it is still a challenging task to estimate the size of uncertainty associated with the predictions of each mass model. In this work, a method called ensemble Bayesian model averaging (EBMA) is introduced to quantify the uncertainty of one-neutron separation energies (S 1 n ) which are directly relevant in the calculations of r -process observables. Here, this Bayesian method provides a natural way to perform model averaging, selection, and uncertainty quantification, by combining the mass models as a mixture of normal distributions whose parameters are optimized against the experimental data, employing the Markov chain Monte Carlo method using the no-u-turn sampler. The EBMA model optimized with all the experimental S 1 n from the AME2003 nuclides are shown to provide reliable uncertainty estimates when tested with the new data in the AME2020.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Real-time Nonlinear Model Predictive Control (NMPC) Strategies using Physics-Based Models for Advanced Lithium-ion Battery Management System (BMS)

Optimal operation of lithium-ion batteries requires robust battery models for advanced battery management systems (ABMS). A nonlinear model predictive control strategy is proposed that directly employs the pseudo-two-dimensional (P2D) model for making predictions. Using robust and efficient model simulation algorithms developed previously, the computational time of the nonlinear model predictive control algorithm is quantified, and the ability to use such models for nonlinear model predictive control for ABMS is established.

25 ENERGY STORAGE↗
DOE OSTI2021

Toward more mechanistic representations of biogeochemical processes in stream and river networks: Implementation and demonstration of a multiscale model : Modeling Archive

This Modeling Archive supports the ORNL-SFA and IDEAS-Waterhsed publication submitted to Environmental Modeling and Software. A recently introduced multiscale model for representing the combined effect of hyporheic exchange flows and small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The ADELS (Advection Dispersion Equation with Lagrangian Subgrid) model uses a one-dimensional advection-dispersion-reaction equation for the channel network and couples that equation at each channel location to a one-dimensional advection-reaction subgrid model representing an ensemble of streamlines that are diverted into the hyporheic zone before returning to the channel. The subgrid model is written in a computationally advantageous Lagrangian form with hyporheic age replacing the hyporheic travel distance. In the paper, we summarized implementation in the integrated surface/subsurface hydrology modeling system Advanced Terrestrial Simulator (ATS).

East Fork Poplar Creek↗
DOE OSTI2025

ORBIT-2 Dataset for Scaling Exascale Vision Foundation Models for Weather and Climate Downscaling

This dataset release corresponds to the work conducted in ORBIT-2: Scaling Exascale Vision Foundation Models for Weather and Climate Downscaling, where large-scale AI methods were applied to improve climate and weather resolution. The collection integrates four widely used, publicly available datasets: ERA5, PRISM, DAYMET, and IMERG. To prepare the data for ORBIT-2 model training and evaluation, we applied a preprocessing pipeline that generates paired low-resolution and high-resolution samples, enabling supervised downscaling experiments. The transformation from coarse to fine scales was performed using bilinear regridding, consistent with the procedures described in WeatherBench2, a community benchmark for weather and climate AI models. This dataset supports the development and evaluation of foundation models designed for weather and climate downscaling at exascale. Additional details on methodology and applications can be found in Wang et al., ORBIT-2 (arXiv:2505.04802, 2025).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Empirical Validation of Multi-Zone HVAC System Model: Evaluation of Existing Infiltration Models used in Building Energy Simulation

Infiltration can have a significant impact on building loads. Studies have shown that infiltration can account for 15-40% of annual space conditioning needs in commercial buildings (Emmerich et al. 2019; Younes et al. 2012). The driving force of infiltration is the pressure difference across the building envelope caused by wind, the stack effect (known as buoyancy effect), and the operation of ventilation equipment. Wind pressure is governed by wind direction, speed, building shape, and other structures around the building. The stack effect is a function of the building height and air density differences of ambient air (Han, 2015). The effect of wind is dominant in low-rise residential buildings, and the stack effect is dominant in highrise buildings (ASHRAE 2017). In building energy simulation programs (e.g., EnergyPlus), various empirical infiltration models (e.g., the effective leakage area model, the flow coefficient model) are available to simulate infiltration rates. To help users in selecting a proper infiltration model for modeling of the two-story Flexible Research Platform (FRP), the team evaluates the existing infiltration models in EnergyPlus based on field measurements from the FRP. The blower door and tracer gas decay tests were performed in the FRP. The blower door test result was used to estimate input parameters required in the infiltration models. The actual infiltration rates were estimated with the tracer gas decay test results.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Modeling Molten Salt Spreading and Heat Transfer using MELTSPREAD – Model Development Updates

This report summarizes spreading and heat transfer results calculated using MELTSPREAD for a scenario of molten FLiNaK spilling onto a flat stainless steel substrate that was inspired by the “maximum credible accident” scenario of the Molten Salt Reactor Experiment (MSRE). MELTSPREAD was developed at Argonne to model the one-dimensional flowing and freezing of molten corium and was applied to model molten salts for the first time during FY21. The molten salt MELTSPREAD model that was developed during FY21 was updated as part of this work using a corrected heat of fusion value of FLiNaK and to accommodate a larger salt spill volume. The model was run with and without the inclusion of contributions from decay heat and a sensitivity analysis of initial spill conditions was performed to determine the importance of those factors on model outcome. The results provide insight into the expected spreading and heat transfer behavior of simulated and irradiated fuel salt that has been spilled and indicate aspects of the model to be revised in future versions of the MELTSPREAD model for applications to molten salt.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗