Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Direct Feed High-Level Waste APPS Model Glass Testing (DFHLW APPS) Matrix, Phase 2

This report summarizes the data collected during the batching and melting of a second matrix of Direct Feed High-Level Waste (DFHLW) glasses generated using the preliminary enhanced waste glass models (EWG2.5) and the Britton and Anderson (2024) preliminary DFHLW feed vector. The purpose of these glasses is two-fold: 1. Validate EWG2.5 glass calculations being used in the Aspen Process Performance Simulation (APPS) model. 2. Evaluate and ultimately improve the glass property models and formulation methods used for design of DFHLW glasses as part of an iterative process of data collection and model refinement. Some of the 16 APPS2 glasses tested did not satisfy all target property constraints due to the limited data on DFHLW glass supporting the EWG2.5 models. • One glass, APPS2-10, formed nepheline on canister centerline cooling (CCC) heat-treatment and failed the product consistency test (PCT) response limits. This glass also had high B and Cr release rates for the toxicity characteristic leaching procedure (TCLP). All other glasses were found to satisfy the PCT and TCLP constraints for both quenched and CCC samples. • One glass, APPS2-08, had higher than acceptable viscosity due to magnetite crystallization. • One glass, APPS2-09, formed greater than 2 vol% crystals at 950 °C. As the glass design criterion was that the temperature at 2 vol% crystal (T 2% ) be less than 950 °C, only one glass failed the criteria. However, this criterion is being reevaluated. Four additional glasses formed crystal fractions between 1 and 2 vol% at 950 °C (APPS2-03, -08, -12, and -14). • Four glasses – APPS2-01, -02, -04, and -16 – failed the Monofrax K-3 refractory neck corrosion (k neck ) design limit of 0.04 in. at 1208 °C for 6 d. This is another criterion being reevaluated. Four additional glasses (APPS2-05, -06, -11, and -13) exhibited 0.025 = k neck = 0.04 in. • All 16 glasses passed the sulfur solubility and TCLP constraints. The measured property values were compared to predicted values using EWG2.5 and a selection of other existing models. A few models (e.g., electrical conductivity, TCLP) were found to be adequate for designing DFHLW glasses in the near future, while others require refits or offsets. It is recommended that new property models be developed for EWG3.0, as a large amount of DFHLW glass property data (> 14 × existing data) is expected to be collected in the compositional spaces where no data was previously available. To enable near-term calculations and formulations for designing DFHLW glasses and processing rate estimations, a formulation algorithm with minor modifications will be developed, EWG2.6.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Physical solvents and techno-economic analysis for pre-combustion CO 2 capture: A review

In this manuscript we review the use, development, and economic performance of physical solvents for pre-combustion CO 2 capture from high pressure H 2 rich syngas streams. Commercially available physical solvents are presented, followed by an assessment of the ideal properties that are important for development of novel solvents for CO 2 capture from high-pressure syngas streams. To compare the technical and economic performance of traditional and novel physical solvents, a review of the methods, assumptions and models used in techno-economic analysis (TEA) studies was conducted. It was found that, although some novel solvents show promising technical performance in the laboratory (e.g., high CO 2 absorption capacity and low vapor pressure), other issues (e.g., solvent viscosity and cost) may limit their industrial applications. Process simulations were useful tools for modeling the technical performance of processes using traditional and novel solvents. However, model predictions are most reliable when the methods and correlations used to develop the process simulation are validated with representative experimental data, in particular highly accurate baseline models are required for fair comparison among physical solvents. The key inputs and assumptions in pre-combustion CO 2 capture TEAs have also been summarized. Some studies showed that the promising technical performance of novel physical solvents can be offset by the high and often unknown costs of these solvents. Future development of novel physical solvents for pre-combustion CO 2 capture will benefit from more studies that conduct in-depth techno-economic analysis, specifically with validated process simulations and transparent economic models.

42 ENGINEERING↗
DOE OSTI2024

Tutorial: The equations of electron emission and their evaluation

Electron emission and transport through and over potential barriers is an essential process requiring modeling and simulation to meet the design needs and characterization of an exceedingly broad range of technologically important devices and processes. The simulation and description of thermal, field, and photoemission, and the related concerns of space–charge affected electron flow, often make use of specialized formulations developed in the early days of quantum mechanics. Advancements in the utilization of electron sources and particularly the simulation of devices and applications using advanced particle-in-cell and trajectory methods for beam optics codes create a strong need for a pedagogical account of the emission models to ensure correct numerical evaluation of their equations. This Tutorial starts from simple phenomenological accounts and progressively builds to comprehensive models emphasizing straightforward and often rapid calculation. It recommends formulations to supplant the canonical Richardson–Laue–Dushman (thermal), Fowler–Nordheim (field), Fowler–DuBridge (photo), and Baroody (secondary) equations and provides a useful formulation of space–charge affected flow commonly described by the Child–Langmuir relation that takes into account cathode dependence on surface field.

Physics↗
DOE OSTI2026

Phase-field model of freeze casting

Directional solidification of water-based solutions has emerged as a versatile technique for templating hierarchical porous materials. However, the underlying mechanisms of pattern formation remain incompletely understood. In this work, we present a detailed derivation and analysis of a quantitative phase-field model for simulating this nonequilibrium process. The phase-field model extends the thin-interface formulation of dilute binary alloy solidification with antitrapping to incorporate the highly anisotropic energetic and kinetic properties of the partially faceted ice-water interface. This interface is faceted in the basal plane normal to the ⟨0001⟩ directions and atomically rough in other directions within the basal plane. On the basal plane, the model reproduces a linear or nonlinear relationship between the interface growth rate and the kinetic undercooling that can be linked to experimental measurements. In both cases, spontaneous parity breaking of the solidification front is observed when the preferred growth direction is aligned with the temperature gradient. This phenomenon leads to the formation of partially faceted ice lamellae that drift laterally in one of the ⟨0001⟩ directions. Here, we demonstrate that the drifting velocity of the ice lamellae is controlled by the kinetics on the basal plane and converges as the thickness of the diffuse solid-liquid interface decreases. Furthermore, we examine the effect of the form of the kinetic anisotropy, which is chosen here such that the inverse of the kinetic coefficient varies linearly from a finite value in the ⟨0001⟩ directions to zero in all other directions within the basal plane, consistent with the assumption that the interface grows in local thermodynamic equilibrium in this plane. Our results indicate that the drifting velocity of ice lamellae is not affected by the slope of this linear relation, and the radius and undercooling at the tip of an ice lamella converge at relatively small slope values. Consequently, the phase-field simulations remain quantitative with computationally tractable choices of both the interface thickness and the slope assumed in the form of the kinetic anisotropy.

Materials science↗
DOE OSTI2019

Arctic Cloud Microphysical Processes. Final report

The primary goal of our team’s effort under this grant was to quantify Arctic lower-troposphere processes by tying data from the large suite of NSA instruments to modeling efforts indirectly through data assimilation and directly through remote sensing and cloud modeling efforts. We maintained a distinct focus on how the ARM measurements may be used evaluate model simulations and/or physical processes. We will separate this report into two distinct sections: our data assimilation effort and the use of NSA instruments to improve our understanding and modeling of model physics to better characterize the critical surface energy budget in polar regimes.

58 GEOSCIENCES↗
DOE OSTI2022

A gradient-based deep neural network model for simulating multiphase flow in porous media

We report simulation of multiphase flow in porous media is crucial for the effective management of subsurface energy and environment-related activities. The numerical simulators used for modeling such processes rely on spatial and temporal discretization of the governing mass and energy balance partial-differential equations (PDEs) into algebraic systems via finite-difference/volume/element methods. These simulators usually require dedicated software development and maintenance, and suffer low efficiency from a runtime and memory standpoint for problems with multi-scale heterogeneity, coupled-physics processes or fluids with complex phase behavior. Therefore, developing cost-effective, data-driven models can become a practical choice, and in this work, we choose deep learning approaches as they can handle high dimensional data and accurately predict state variables with strong nonlinearity. In this paper, we describe a gradient-based deep neural network (GDNN) constrained by the physics related to multiphase flow in porous media. We tackle the nonlinearity of flow in porous media induced by rock heterogeneity, fluid properties, and fluid-rock interactions by decomposing the nonlinear PDEs into a dictionary of elementary differential operators. We use a combination of operators to handle rock spatial heterogeneity and fluid flow by advection. Since the augmented differential operators are inherently related to the physics of fluid flow, we treat them as first principles prior knowledge to regularize the GDNN training. We use the example of pressure management at geologic CO 2 storage sites, where CO 2 is injected in saline aquifers and brine is produced, and apply GDNN to construct a predictive model that is trained with physics-based simulation data and emulates the physics process. We demonstrate that GDNN can effectively predict the nonlinear patterns of subsurface responses, including the temporal and spatial evolution of the pressure and saturation plumes. We also successfully extend the GDNN to convolutional neural network (CNN), namely gradient-based CNN (GCNN), and validate its capability to improve the prediction accuracy. GDNN has great potential to tackle challenging problems that are governed by highly nonlinear physics and enable the development of data-driven models with higher fidelity.

42 ENGINEERING↗
DOE OSTI2025

Effect of toolpath in large-format additive manufacturing with bio-derived composites

There has been growing interest in integrating bio-derived composites into Large-format additive manufacturing (LFAM) feedstocks to reduce the use of petroleum-derived materials and reduce the overall carbon footprint of LFAM. However, these materials present unique challenges during manufacturing due to their variability, which can lead to unintended deformations and failures attributable to suboptimal process conditions. While numerical modelling has been extensively employed to simulate numerous manufacturing processes, its application in LFAM with bio-derived composites remains limited. This study addresses this gap by systematically developing a numerical model to simulate the LFAM process using bio-based materials, specifically wood fiber-reinforced polylactic acid (PLA/WF). Experimental investigations were conducted to characterise the thermal and mechanical properties of additively manufactured PLA/WF specimens. Numerical simulations were performed to predict temperature profiles and deformations during LFAM. The effect of varying infill patterns, internal structures, and tool paths on the temperature distribution and deformation of printed parts was explored using the developed model. This article aims to advance the utilisation of bio-derived composites in LFAM systems and provide a comprehensive understanding of the LFAM process. The findings offer valuable insights for optimising process parameters and enhancing the performance of LFAM with bio-based composites.

Large-format additive manufacturing↗
DOE OSTI2023

Understanding Dynamics and Thermodynamics of ENSO and Its Complexity Simulated by E3SM and Other Climate Models

Despite the seeming success of most state-of-the-art climate models in simulating the El Niño-Southern Oscillation (ENSO), there is strong evidence that models achieve realistic levels of ENSO activity do so often for wrong reasons. This is owing to an often occurred near cancelation of large errors in terms of contributions to ENSO growth rate from coupled dynamic and thermodynamic feedback processes. Climate models remain deficient in simulating the observed ENSO’s spatial and temporal complexity that involves interplays of coupled dynamic and thermodynamic feedbacks, interactions across multiple scales, nonlinear processes in the tropical atmosphere and ocean system, biases in mean sate and physical processes, and influences external to equatorial Pacific coupled ENSO dynamics. Our proposed research aims at advancing predictive and process-level understandings of ENSO simulated in E3SM and other climate models under current and future climate conditions with two main objectives: (i) better understanding the aforementioned broad range interactive processes and sources that control fundamental properties of ENSO in E3SM and CIMP6 outputs as well as in observational (reanalysis) data sets, using a hierarchical of coupled dynamical frameworks consisting of theoretical analysis, intermediate complexity modeling, and coupled dynamic diagnostics; (ii) to use this understanding to explore pathways towards improving E3SM’s capability of simulating ENSO and its complexity. More specifically, we will focus on four main thrusts of research: (1) ENSO’s dynamic and thermodynamic feedbacks; (2) the across-scale interactions of ENSO with annual cycle and MJO/WWB/TIW (Madden Julian Oscillation/Westerly Wind Burst/Tropical Instability Wave) activity; (3) key nonlinear processes of ENSO involving atmospheric convective thresholds, nonlinear ocean dynamic heating, and thermocline outcropping; and (4) the impacts of climate mean-state biases/changes and perturbed physical processes on simulated ENSO and its complexity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Reducing impacts of artificial ponding in modeling salt marshes using a conductivity-formulated subgrid model

The landforms of salt marsh systems are rather complex, with meandering channel networks cutting through low-lying, extensive marsh platforms. Some small-scale topographic details are critical to flooding and draining processes but cannot be reflected in the Digital Elevation Model (DEM) due to limited measurement resolution and vegetation bias. Numerical models need a grid with high enough resolution to capture the small-scale connectivity to correctly model the small-scale flow processes, making the model computationally unaffordable. Here, in this study, we derived a subgrid model, which couples two flow components, one is the regular subgrid model for simulating the subgrid-resolved flows, and the other is a conductivity-formulated porous flow model for modeling small-scale gut flows unresolved at the subgrid level. The model validation was performed using an idealized case with a narrow slot that was unresolved in the coupled flow model. Model results were compared with a slot-resolved regular subgrid model where only surface flows are solved and suggested that, with a proper model parameterization, the coupled flow model is capable of reproducing flooding and draining processes similar to the slot-resolved subgrid model of surface flow. The model was applied to simulating flooding and draining processes in a meso-tidal salt marsh. The model results demonstrate that the coupled subgrid model can be used to reduce the impacts of artificial ponding that usually occur in numerical simulations of salt marshes with numerous narrow gullies and creeks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A numerical study of process complexity in permafrost dominated regions

Numerical modeling of permafrost dynamics requires adequate representation of atmospheric and surface processes, a reasonable parameter estimation strategy, and site-specific model development. The three main research objectives of the study are: (i) to propose a novel methodology that determines the required level of surface process complexity of permafrost models by conducting parameter sensitivity and calibration, (ii) to design and compare three numerical models of increasing surface process complexity, and (iii) to calibrate and validate the numerical models at the Yakou catchment on the Qinghai-Tibet Plateau as an exemplary study site. The calibration was carried out by coupling the Advanced Terrestrial Simulator (numerical model) and PEST (calibration tool). Simulation results showed that (i) A simple numerical model that considers only subsurface processes can simulate active layer development with the same accuracy as other more complex models that include surface processes. (ii) Peat and mineral soil layer permeability, Van Genuchten alpha, and porosity are highly sensitive. (iii) Liquid precipitation aids in increasing the rate of permafrost degradation. (iv) Deposition of snow insulated the subsurface during the thaw initiation period. We have developed and released an integrated code that couples the numerical software ATS to the calibration software PEST. The numerical model can be further used to determine the impacts of climate change on permafrost degradation.

Calibration↗
DOE OSTI2023

Comparing quantile regression forest and mixture density long short-term memory models for probabilistic post-processing of satellite precipitation-driven streamflow simulations

Abstract. Deep learning (DL) and machine learning (ML) are widely used in hydrological modelling, which plays a critical role in improving the accuracy of hydrological predictions. However, the trade-off between model performance and computational cost has always been a challenge for hydrologists when selecting a suitable model, particularly for probabilistic post-processing with large ensemble members. This study aims to systematically compare the quantile regression forest (QRF) model and countable mixtures of asymmetric Laplacians long short-term memory (CMAL-LSTM) model as hydrological probabilistic post-processors. Specifically, we evaluate their ability in dealing with biased streamflow simulations driven by three satellite precipitation products across 522 nested sub-basins of the Yalong River basin in China. Model performance is comprehensively assessed using a series of scoring metrics from both probabilistic and deterministic perspectives. Our results show that the QRF model and the CMAL-LSTM model are comparable in terms of probabilistic prediction, and their performances are closely related to the flow accumulation area (FAA) of the sub-basin. The QRF model outperforms the CMAL-LSTM model in most sub-basins with smaller FAA, while the CMAL-LSTM model has an undebatable advantage in sub-basins with FAA larger than 60 000 km2 in the Yalong River basin. In terms of deterministic predictions, the CMAL-LSTM model is preferred, especially when the raw streamflow is poorly simulated and used as input. However, setting aside the differences in model performance, the QRF model with 100-member quantiles demonstrates a noteworthy advantage by exhibiting a 50 % reduction in computation time compared to the CMAL-LSTM model with the same ensemble members in all experiments. As a result, this study provides insights into model selection in hydrological post-processing and the trade-offs between model performance and computational efficiency. The findings highlight the importance of considering the specific application scenario, such as the catchment size and the required accuracy level, when selecting a suitable model for hydrological post-processing.

Geology↗
DOE OSTI2019

Model-driven Deception for Control System Environments

Deception techniques are a useful defensive mechanism that could be very valuable in control systems where updates and patches can be difficult. However, to deceive a motivated and targeted attacker it is necessary to enhance deception platforms with capabilities to model and simulate physical processes to achieve necessary fidelity. In this paper, attributes of cyberphysical decoys are discussed and design of a system with these attributes is presented. A boiler model use case is provided to demonstrate how this novel deception capability can be integrated into and used to defend a system.

Hofer, William J.↗
DOE OSTI2025

Direct Air Reactive Capture and Conversion for Utility-Scale Energy Storage (Final Report)

This final report for FEW0277 summarizes the work performed over the project performance period of October 2021 – March 2025. This project was funded under the “Reactive Capture and Conversion R&D” lab call released in FY2021. The goal of the project was to develop dual-function materials and process for capturing CO 2 from the atmosphere and converting it into CH 4 . The work was organized into four parallel tracks in 1) direct air capture materials synthesis and characterization, 2) catalysts for CO 2 conversion, 3) mechanistic investigations via ab initio simulations, and 4) process modeling, technoeconomic analysis, and lifecycle assessment. The project was split into two budget periods. The first budget period focused on development of amine-based materials, due to their known performance for CO 2 direct air capture and their potential to act synergistically with metal catalysts to enable a low-temperature methanation pathway. The second budget period focused on development of alkali-based materials and a simulated-moving-bed process for high conversion catalytic reduction of captured CO 2 to CH 4 . All project milestones were completed during the project performance period and are summarized in this report. Our work resulted in publication of eight peer-reviewed manuscripts, one patent application, and numerous presentations given at domestic and international conferences and invited academic department seminars.

03 NATURAL GAS↗
DOE OSTI2023

A Three-Tier Incremental Approach for Development of Smart Corridor Digital Twins

Development of an arterial smart corridor digital twin requires the integration of real time data streams, data storage, and simulation model execution. This process may be complex, time consuming, and susceptible to errors. To aid in smart corridor digital twin development this paper seeks to provide a framework, best practices, and development guidance. As such, a three-tier incremental approach to smart corridor digital twin development is presented. The paper highlights practical issues in digital twin construction along with key data challenges based on experiences from the development of digital twins for two large Smart Corridor deployments, one in Chattanooga, Tennessee, and the other in Atlanta, Georgia. The presented three-tier incremental approach includes: 1) development of a prepopulated offline simulation, 2) development of a pseudo digital twin that is driven by archived data, and 3) integration of real time data streams to create the online digital twin model. The three-tiered approach facilitates conducting multiple trials and scenarios with increasing complexity, allowing for incremental error processing and updating of the digital twin.

Saroj, Abhilasha↗
DOE OSTI2025

Biogenic carbon capture at pulp mills via sodium spiking and oxy-fuel calcination

Over 13 million metric tons of biogenic carbon dioxide (CO 2 ) are mineralized yearly in United States (US) pulp mill recovery boilers as molten sodium carbonate. The mineralized CO 2 is released downstream in the rotary lime kiln where it can be captured at a relatively low cost for permanent sequestration. The use of biocarbon and bioenergy at pulp mills enables the possibility for atmospheric carbon removal when biogenic CO2 is captured and sequestered. We demonstrate the feasibility of capturing biogenic CO 2 via sodium spiking coupled with oxy-fuel calcination in the rotary lime kiln at existing kraft pulp mills. Sodium (Na) spiking elevates the Na ion loading in the kraft chemical looping process by replacing chlorine-based bleaching with a highly alkaline bleaching sequence. Chemical pulping processes are modeled and simulated to understand the technical limits of implementing sodium spiking in existing pulp mills. Each 1 % of sodium added to the bleaching operations increases the rate of CO 2 mineralization by 4 %, and the maximum increase in sodium content in the kraft process is 9 %. Without sodium spiking, estimated CO 2 capture costs are $131 and $107/mt CO 2 for air and oxy-fuel combustion, respectively. With the implementation of sodium spiking, the cost of CO 2 captured decreases by 27 % and 31 %, showing costs of $96 and $74/mt CO 2 for air and oxy-fuel combustion, respectively. Oxy-fuel combustion reduces the costs of CO 2 capture compared to air combustion, but merged deployment with sodium spiking seems to be the most cost-effective pathway to integrate carbon capture.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Anthropogenic extremely low volatility organics (ELVOCs) Govern the Growth of Molecular Clusters over the Southern Great Plains during the Springtime

New particle formation (NPF) and growth govern cloud condensation nuclei (CCN) concentrations in many regions. The mechanisms governing the nucleation of molecular clusters vary substantially in different regions of the atmosphere. Additionally, the growth of these clusters from ~2 to 20 nm sizes is often governed by the availability of extremely low volatility organic vapours (ELVOCs). While the pathways to ELVOC formation from the oxidation of biogenic monoterpenes with ozone is better understood, the chemical and mechanistic pathways for ELVOC formation from oxidation of anthropogenic organics are not well understood. We integrate measurements and three-dimensional regional model simulations with the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to understand the processes governing new particle formation and growth and secondary organic aerosol (SOA) formation during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems (HI-SCALE) field campaign at the Southern Great Plains (SGP) observatory in Oklahoma, and contrast it with a site within the Bankhead National Forest (BNF), Alabama in Southeast USA, where 5-year long measurements will begin in 2024. Simulations show that nucleation rates are at least an order of magnitude higher at SGP compared to BNF during the springtime days (April 28 and May 14, 2016), largely due to lower H2SO4 concentrations at BNF, which are needed for nucleation. In addition, the larger CS at BNF (compared to SGP) increase the loss of molecular clusters by coagulation to pre-existing particles. Among the 8 different nucleation mechanisms in WRF-Chem, we find that the amine+H2SO4 nucleation mechanism dominates at the SGP site, while the pure organic ion induced nucleation mechanism dominates over BNF. Through various WRF-Chem sensitivity simulations, we find that anthropogenic ELVOCs are critical for explaining the growth of newly formed particles and the resulting number size distribution observed near the surface at the SGP site during the daytime. In addition, we show that treating organic particles as semisolid, with strong diffusion-limited uptake of organic vapours, brings model predictions into closer agreement with the observed evolution of particle size distribution. Simulations also predict that anthropogenic SOA, formed by the oxidation of aromatic volatile organic compounds (VOCs), is the dominant organic aerosol component at SGP, while biogenic SOA dominates particle composition at the BNF site in Southeast USA on these days.

Shrivastava, ManishKumar B.↗
DOE OSTI2022

Quantum signal processing for simulating cold plasma waves

Numerical modeling of radio-frequency waves in plasma with sufficiently high spatial and temporal resolution remains challenging even with modern computers. However, such simulations can be sped up using quantum computers in the future. In this work, we propose how to do such modeling for cold plasma waves, in particular, for an X wave propagating in an inhomogeneous one-dimensional plasma. The wave system is represented in the form of a vector Schrödinger equation with a Hermitian Hamiltonian. Block encoding is used to represent the Hamiltonian through unitary operations that can be implemented on a quantum computer. To perform the modeling, we apply the so-called quantum signal processing algorithm and construct the corresponding circuit. Quantum simulations with this circuit are emulated on a classical computer, and the results show agreement with traditional classical calculations. We also discuss how our quantum circuit scales with the resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗