Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

A Vertically Flow-Following, Icosahedral Grid Model for Medium-Range and Seasonal Prediction. Part 1: Model Description

A hydrostatic global weather prediction model based on an icosahedral horizontal grid and a hybrid terrain following/ isentropic vertical coordinate is described. The model is an extension to three spatial dimensions of a previously developed, icosahedral, shallow-water model featuring user-selectable horizontal resolution and employing indirect addressing techniques. The vertical grid is adaptive to maximize the portion of the atmosphere mapped into the isentropic coordinate subdomain. The model, best described as a stacked shallow-water model, is being tested extensively on real-time medium-range forecasts to ready it for possible inclusion in operational multimodel ensembles for medium-range to seasonal prediction.

mathematical models↗
NASA NTRS2017

Mars Propellant Liquefaction and Storage Performance Modeling using Thermal Desktop with an Integrated Cryocooler Model

NASAs current Mars architectures are assuming the production and storage of 23 tons of liquid oxygen on the surface of Mars over a duration of 500+ days. In order to do this in a mass efficient manner, an energy efficient refrigeration system will be required. Based on previous analysis NASA has decided to do all liquefaction in the propulsion vehicle storage tanks. In order to allow for transient Martian environmental effects, a propellant liquefaction and storage system for a Mars Ascent Vehicle (MAV) was modeled using Thermal Desktop. The model consisted of a propellant tank containing a broad area cooling loop heat exchanger integrated with a reverse turbo Brayton cryocooler. Cryocooler sizing and performance modeling was conducted using MAV diurnal heat loads and radiator rejection temperatures predicted from a previous thermal model of the MAV. A system was also sized and modeled using an alternative heat rejection system that relies on a forced convection heat exchanger. Cryocooler mass, input power, and heat rejection for both systems were estimated and compared against sizing based on non-transient sizing estimates.

fluid modeling↗
NASA NTRS2024

Future CH4 Budgets as Modelled By A Fully Coupled Earth System Model Using Prescribed GHG Concentrations vs. Interactive CH4 Sources and Sinks

We have used the NASA Goddard Institute for Space Studies (GISS) Earth system model GISS-E2.1 to study the future budgets and trends of global and regional CH4 under different emission scenarios. GISS-E2.1 is one of the few ESMs that can be driven by anthropogenic CH4 emissions, as well as interactive natural sources such as wetlands, and can simulate the tropospheric CH4 chemistry. In frame of the recent short-lived climate forcers (SLCFs) assessment by the Arctic Monitoring and Assessment Programme (AMAP), we used the GISS-E2.1 model with prescribed long-lived greenhouse gas (GHG) concentrations. In the present study, we have supplemented these simulations using the interactive CH4 sources and sinks in order to quantify the model performance and the sensitivity to CH4 sources and sinks. We have used the Current Legislation (CLE) and the Maximum Feasible Reduction (MFR) emission scenarios from the Eclipse V6b emission database to simulate the future chemical composition and climate impacts from 2015 to 2050. We have also simulated 1995-2014 in order to evaluate the model performance following the AMAP-SLCF protocol. The prescribed GHG version underestimates the Global Atmospheric Watch (GAW) surface CH4 observations during the period between 1995 and 2023 by 1% [-8.4%-2.0%], with a correlation (r) of 0.71 [-0.41 0.99]. The largest underestimations are over the continental emission regions such as North America, Europe, and Asia, while biases are smallest over oceans. On the other hand, the simulation with interactive sources and sinks underestimates the GAW observations more than the prescribed simulation, by 18.5% [-25% -10.4%], with a lower r of 0.36 [-0.82 0.93]. Opposite to the prescribed simulation, the biases are largest over oceans and smaller over the continents, however they are still larger over land than the prescribed simulation. The interactive simulation, with large sources virtually over land and strong sink over oceans, has a land/ocean ratio larger than 1 while the prescribed simulation has this ratio equal to 1 as it distributes the global prescribed CH4 concentration equally in longitude over a given latitude. This clearly shows that the interactive sources and sinks should be represented in models in order to realistically simulate the chemical composition and the oxidative capacity of the atmosphere. As expected, the MFR scenario simulates lower global surface CH4 concentrations and burdens compared to the CLE scenario, however in both cases, global surface CH4 and burden continue to increase through 2050 compared to present day. In the CLE scenario, increases are largest over the equatorial belt, in particular over India and East China, while the MFR scenario shows increases over the whole Southern Hemisphere, however much smaller compared to CLE. Finally, the interactive simulation shows that the chemical CH4 sink increases in the CLE scenario, while it slightly decreases in the MFR, leading to a larger CH4 lifetime in the MFR scenario compared to in the CLE scenario.

CH4↗
NASA NTRS2024

Enabling Advanced Snow Physics within Land Surface Models Through an Interoperable Model-Physics Coupling Framework

Accurate estimation of snow accumulation and melt is a critical part of decision-making in snow-dominated watersheds. In this study, we demonstrate a flexible methodology to couple a detailed snow model, Crocus, separately to two different land surface models (LSMs), Noah-MP and Noah. The original LSMs and the coupled models (Noah-MP-Crocus and Noah-Crocus) are used to simulate snow depth, snow water equivalent, and other water and energy states and fluxes. The results of simulations are compared against a wide range of independent gridded and point scale reference datasets. Our results show that coupling the detailed snow model, Crocus, with the LSMs improves the snow depth and snow water equivalent relative to independent observations. Overall, larger improvements are obtained with coupling Crocus to the Noah LSM, with the coupled Noah-Crocus configuration reducing the RMSE and bias of snow depth from 2-12% and 57-75%, respectively, relative to Snow Data Assimilation System (SNODAS) and snow product from the University of Arizona. On the other hand, smaller improvements are obtained by coupling Crocus with Noah-MP. The Coupled Noah-MP-Crocus reduces the snow depth bias but slightly degrades the RMSE of snow depth and snow water equivalent. The corresponding impacts in other water budget terms such as evapotranspiration, soil moisture, and streamflow, however, are mixed, pointing to the significant need to improve the coupling assumptions of these processes within land models. Overall, the interoperable coupling framework demonstrated here offers the opportunity to include more detailed snow physics and processes, and to advance data assimilation systems through improved exploitation of information from snow remote sensing instruments.

Snow model↗
DOE OSTI2024

Fuel Performance Modeling Status Update and Potential Model Improvements

Fuel performance modeling status update and potential model improvements overview of TRISO fuel performance modeling codes PARFUME/BISON, AGR experiment support, potential modeling improvements, BISON smeared cracking model, BISON fission product source term, and AGR-3/4 reirradiation heating test.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Modeling atmospheric data and identifying dynamics Temporal data-driven modeling of air pollutants

Atmospheric modelling has recently experienced a surge with the advent of deep learning. Most of these models, however, predict concentrations of pollutants following a data-driven approach in which the physical laws that govern their behaviors and relationships remain hidden. Furthermore, with the aid of real-world air quality data collected hourly in different stations throughout Madrid, we present a case study using a series of data-driven techniques with the following goals: (1) Find systems of ordinary differential equations that model the concentration of pollutants and their changes over time; (2) assess the performance and limitations of our model using stability analysis; (3) reconstruct the time series of chemical pollutants not measured in certain stations using delay coordinate embedding results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Model Data Archive for Manuscript Titled "Evaluation of a Coupled Surface–Subsurface Hydrologic Model Using Dense Water‑Level Sensors in a Mixed Urban–Rural Watershed"

This archive provides scripts, input files, and datasets used for the implementation and evaluation of a fully coupled surface–subsurface hydrologic model in the Neches River Basin, southeast Texas. The study uses the Advanced Terrestrial Simulator (ATS) to simulate coupled surface–subsurface hydrologic processes over a mixed urban–rural watershed and evaluates model performance using a dense network of 136 in situ water-level sensors, nine U.S. Geological Survey (USGS) stream gauges, and SSEBop-derived evapotranspiration estimates during the period October 2014–June 2024. The workflow is implemented primarily in Python 3 using the Watershed Workflow package. The Jupyter notebooks can be executed using open-source software such as Anaconda JupyterLab or Visual Studio Code. Other data files include TXT, CSV, XML, SHP, TIF, NetCDF, HDF5, and ExodusII files, which can be processed using the provided Python scripts. ATS input files are provided in XML format and can be edited using any commonly used text editor. This archive contains: *Scripts and input files used to generate the ATS model setup, including watershed discretization, mesh generation, parameter mapping, and model configuration. *Jupyter notebooks used for preprocessing observational data, evaluating streamflow, water levels, and evapotranspiration, computing performance metrics, and generating the figures presented in the manuscript. *ATS simulation outputs and processed observational datasets, including OneRain and DD6 water-level sensors, USGS streamflow observations, GIS data, and supporting spatial datasets used throughout the study.

Dense water-level sensor network↗
DOE OSTI2020

Learning compact physics‐aware delayed photocurrent models using dynamic mode decomposition

Abstract Radiation‐induced photocurrent in semiconductor devices can be simulated using complex physics‐based models, which are accurate, but computationally expensive. This presents a challenge for implementing device characteristics in high‐level circuit simulations where it is computationally infeasible to evaluate detailed models for multiple individual circuit elements. In this work we demonstrate a procedure for learning compact delayed photocurrent models that are efficient enough to implement in large‐scale circuit simulations, but remain faithful to the underlying physics. Our approach utilizes dynamic mode decomposition (DMD), a system identification technique for learning reduced‐order discrete‐time dynamical systems from time series data based on singular value decomposition. To obtain physics‐aware device models, we simulate the excess carrier density induced by radiation pulses by solving numerically the ambipolar diffusion equation, then use the simulated internal state as training data for the DMD algorithm. Our results show that the significantly reduced‐order delayed photocurrent models obtained via this method accurately approximate the dynamics of the internal excess carrier density—which can be used to calculate the induced current at the device boundaries—while remaining compact enough to incorporate into larger circuit simulations.

Hanson, Joshua↗
DOE OSTI2021

A microstructure-based modeling approach to predict the mechanical properties of Zr alloy with hydride precipitates

In nuclear reactors, hydrides can form in fuel cladding due to hydrogen absorption in ZIrcaloy and cause embrittlement. This work presents a microstructure-based finite element model to predict the stress-strain response of Zircaloy containing hydrides. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. The constitutive material model for zircaloy in this study is based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as brittle material along with a damage model. Hydride formation inside the zircoloy matrix results in residual stress. This phenomenon is also captured in this model. A parametric study has been conducted to understand the effect of volume fraction, orientation, and lamellae thickness of the hydride phase on the mechanical properties of the overall material.

Kulkarni, Shank S.↗
DOE OSTI2022

Leaching Model of Radionuclides in Metal-Organic Framework Particles

Metal-organic frameworks (MOFs) have been used to sequester radionuclides and seal them inside of porous scaffolds using postsynthetic modification procedures. Experiments show that certain Zr-MOF with different capping linkers significantly affects the radionuclide release kinetics. In this work, we developed a leaching model of radionuclides in Zr-MOF particles. The model assumes that uranyl species occupy two energetically favored sites: the metal node and the MOF pores. For a given overall concentration of uranyl species, the partitions of uranyl species at the metal nodes and within the pores are determined by their chemical potentials. The model also considers the effect of particle surface and concentration on chemical potentials and diffusivity. The effects of spatial and structural dependent chemical potentials and diffusivity as well as particle sizes on leaching kinetics are investigated with the model. Predicted and measured uranyl leaching kinetics in Zr-MOF particles under batch experiments are compared. The results demonstrate the model’s capability for exploring the mechanisms of leaching and provide guidance for material design.

Metal-Organic Framework, Diffusion, uranyl, Leachi↗
DOE OSTI2021

A mechanistic damage model for solid oxide fuel cell ceramic materials - Part I: Constitutive modeling

A multiscale mechanistic damage model for solid oxide fuel cell (SOFC) ceramics that combines micromechanics of stiffness reduction due to material porosity change and microcracking with a continuum damage mechanics (CDM) description for the evolution of microcracks up to fracture is developed in this work. The model also accounts for volumetric swelling of the anode due to redox cycling. Porosity change and swelling during redox combined with normal operating thermomechanical loads could increase the stress, strain and damage distributions in SOFC stacks leading them to failure. The proposed model involves three governing parameters: (i) the porosity in terms of the pore or void volume fraction, (ii) material swelling magnitude prescribed at a given loading step, and (iii) the damage variable that describes microcracking caused by thermomechanical loads and swelling. Pores and microcracks in the ceramics are modeled as randomly distributed ellipsoidal inclusions with negligible stiffness by an Eshelby-Mori-Tanaka formulation combined with an inclusion orientation distribution method. Microcracking damage evolution is described by a CDM formulation. Volumetric swelling is treated in a similar way to thermal expansion in the constitutive relations. After validation and numerical checks, the damage model is used to analyze a simplified SOFC positive-electrode/electrolyte/negative-electrode (PEN) structure subjected a redox cycle in addition to thermomechanical loads experienced during normal operation of a SOFC stack.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Molecular Simulation of Lithium Carbonate Reactive Vapor–Liquid Equilibria Using a Deep Potential Model

We developed a first-principles machine learning model for the reactive vapor–liquid phase behavior of molten Li 2 CO 3 . The model was trained on ab initio electronic density functional theory data using the Deep Potential (DP) methodology, and its accuracy was evaluated by comparing model predictions of density and viscosity to experimental measurements. Direct coexistence simulations with the DP model over time scales of tens of nanoseconds were used to observe equilibrium dissociation of Li 2 CO 3 into CO 2 residing primarily in the vapor phase and Li 2 O which remains dissolved in the liquid. The simulations covered a range of temperatures, overall system sizes, and vapor-to-liquid volume ratios. Results were analyzed in terms of the observed chemical composition of the liquid and vapor phases, product structure, and CO 2 partial pressures. In addition, we calculated equilibrium constants for the dissociation reaction by assuming ideal-solution behavior for the liquid. As expected on the basis of thermodynamic arguments and prior experiments for this system, the observed partial pressure of CO 2 in the gas phase depends on both the temperature and the ratio of vapor to liquid volumes, while the calculated equilibrium constants only depend on temperature. DP model predictions for the equilibrium constant of the reaction are generally consistent with the available experimental measurements. Furthermore, the present study establishes the validity of the DP methodology for the description of reactive, multiphase equilibria from first principles, with possible applications to many other systems of scientific and technological interest even in the absence of relevant experimental measurements.

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

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to stop unphysical rearrangements during geometry relaxation. A simple fixed-atom or “coordinate-lock” approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

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

Exploring topography-based methods for downscaling subgrid precipitation for use in Earth System Models

Topography exerts major control on land surface processes. To improve representation of topographic impacts on land surface processes, a new topography-based subgrid structure has been introduced to the Energy Exascale Earth System Model to represent the subgrid heterogeneity of surface elevation. To take advantage of the new subgrid structure for improving land surface modeling, this study explores four variations of topography-based methods for downscaling grid precipitation to the corresponding subgrids. In the first three methods, the deviation of the subgrid precipitation from the grid’s is equal to the grid precipitation multiplied by the ratio of the elevation difference between the subgrid and grid mean to a specified elevation equals to the grid elevation, the difference between the maximum and minimum subgrid elevation, and the maximum subgrid elevation, respectively. The second method limits the ratio to 0.5 to avoid extreme values on mountains and the third method accounts for the slope effect. The fourth method is similar to the third method except that the Froude Number is used to limit the numerator in a blocking regime of the ambient flow. The downscaled precipitation is evaluated using the PRISM precipitation data over the U.S. using statistical metrics. Results show that by accounting for topographic slope besides elevation, the third and fourth methods show clear advantages over the first and second methods. Furthermore, introducing the Froude Number in the fourth method improves downscaling skill and shows consistent advantages over the third method in areas with larger subgrid heterogeneity across different grids sizes.

Tesfa, Teklu K.↗
DOE OSTI2022

Using an Explainable Machine Learning Approach to Characterize Earth System Model Errors: Application of SHAP Analysis to Modeling Lightning Flash Occurrence

Abstract Computational models of the Earth System are critical tools for modern scientific inquiry. Efforts toward evaluating and improving errors in representations of physical and chemical processes in these large computational systems are commonly stymied by highly nonlinear and complex error behavior. Recent work has shown that these errors can be effectively predicted using modern Artificial Intelligence (A.I.) techniques. In this work, we go beyond these previous studies to apply an explainable A.I. technique to not only predict model errors but also move toward understanding the underlying reasons for successful error prediction. We use XGBoost classification trees and SHapley Additive exPlanations analysis to explore the errors in the prediction of lightning occurrence in the NASA Goddard Earth Observing System model, a widely used Earth System Model. This explainable error prediction system can effectively predict the model error and indicates that the errors are strongly related to convective processes and the characteristics of the land surface.

54 ENVIRONMENTAL SCIENCES↗