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 91 records · Page 5

Approach and Model Used to Represent a Timeline Analysis for Security Design Enhancements

Next-generation reactors will be able to use risk to inform and performance base the licensing of many aspects of the reactor, facility, and site design, including attributes of physical security. There are several factors related to security, including site topography, reactor design, and physical protection system components, to consider when designing the physical protection system into the overall facility design and plan of operation. With the versatility of advanced reactors, especially micro reactors, methods are needed to simplify and quickly evaluate potential timelines for designing a site configuration. This report describes an approach to generate qualitative and quantitative insights using a risk-informed simulation. The modeling process is described in detail, focusing on three aspects: (1) the facility mission time (the time required to control the plant until safe), (2) the attacker timeline (the time to potential sabotage), and (3) the response timeline (the time to counter the facility attack). The modeling capabilities also are extended to include facility physical phenomena such as thermal-hydraulics and heat transfer to capture realistic representation of dynamic changes to a facility. While the plant models and examples are hypothetical and do not represent a real facility, these modeling approaches could be used for future security-by-design engineering in advanced reactors. The outputs and insights from the modeling approach may be used to modify and optimize the security posture of a facility by efficiently making modifications to the model and seeing the overall impact from the modification. Lastly, use of the approach described in this report can also provide the technical basis for a physical protection program, describing how the facility and security strategy will cope with off-normal events.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Integrated system failure analysis software toolchain (IS-FAST)

Systems and methods are provided for generating faults and analyzing fault propagation and its effects. Starting from the ontologies of components, functions, flows, and faults, systems and methods are provided that describe, generate and track faults in a computer system across multiple domains throughout design and/or development. In order to construct the system and fault models, a series of concepts is introduced in the form of ontologies and their dependencies. An investigation is performed into the faults, including their type, cause, life-cycle aspects, and effect. Principles and rules are created to generate various faults based on system configurations. After the modeling process, a simulation engine is described to execute actions and simulate the process of fault generation and propagation. As a result, fault paths that impact components and functions can be obtained.

Diao, Xiaoxu↗

Lignin Utilization

Valorization of lignin has massive potential economic and sustainability benefits for the lignocellulosic biorefinery. However, challenges remain to realize lignin conversion to coproducts, especially related to selective, high yield depolymerization to monomers and quantitative analytics on lignin, the latter of which is critical for accurate process modeling. Towards these goals, the Lignin Utilization project focuses on catalytic lignin deconstruction chemistries for C-O and C-C bond cleavage, development of analytical chemistry techniques to quantitatively characterize lignin, and syntheses of requisite model compounds for understanding lignin transformations. This work is done with and supports multiple BETO projects, including the Biological Lignin Valorization (BLV) project, the Separations Consortium, and others. Outcomes of Lignin Utilization for analytics and synthesis include 1) development of new mass spectrometry (MS) methods to characterize lignin dimers/oligomers in process streams, 2) deployment of a computational-experimental tool (with Biochemical Process Modeling and Simulation) to identify lignin-derived compounds with high fidelity from MS, and 3) delivery of >40 unique compounds. From a catalysis perspective, we have developed new oxidative approaches to cleave C-O and C-C bonds and produce >50% bio-available aromatic monomers for biological funneling for the BLV project and have developed recoverable bases for base-catalyzed deconstruction of lignin.

BIOMASS FUELS↗

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗

Mean Climate and Tropical Rainfall Variability in Aquaplanet Simulations Using the Model for Prediction Across Scales‐Atmosphere

Abstract Aquaplanet experiments are important tools for understanding and improving physical processes simulated by global models; yet, previous aquaplanet experiments largely differ in their representation of subseasonal tropical rainfall variability. This study presents results from aquaplanet experiments produced with the Model for Prediction Across Scales‐Atmosphere (MPAS‐A)—a community model specifically designed to study weather and climate in a common framework. The mean climate and tropical rainfall variability simulated by MPAS‐A with varying horizontal resolution were compared against results from a recent suite of aquaplanet experiments. This comparison shows that, regardless of horizontal resolution, MPAS‐A produces the expected mean climate of an aquaplanet framework with zonally symmetric but meridionally varying sea‐surface temperature. MPAS‐A, however, has a stronger signal of tropical rainfall variability driven by convectively coupled equatorial waves. Sensitivity experiments with different cumulus parameterizations, physics packages, and vertical grids consistently show the presence of those waves, especially equatorial Kelvin waves, in phase with lower‐tropospheric convergence. Other models do not capture such rainfall‐kinematics phasing. These results suggest that simulated tropical rainfall variability depends not only on the cumulus parameterization (as suggested by previous studies) but also on the coupling between physics and dynamics of climate and weather prediction models.

Rios‐Berrios, R.↗

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗

Bayesian stability and force modeling for uncertain machining processes

Accurately simulating machining operations requires knowledge of the cutting force model and system frequency response. However, this data is collected using specialized instruments in an ex-situ manner. Bayesian statistical methods instead learn the system parameters using cutting test data, but to date, these approaches have only considered milling stability. This paper presents a physics-based Bayesian framework which incorporates both spindle power and milling stability. Initial probabilistic descriptions of the system parameters are propagated through a set of physics functions to form probabilistic predictions about the milling process. The system parameters are then updated using automatically selected cutting tests to reduce parameter uncertainty and identify more productive cutting conditions, where spindle power measurements are used to learn the cutting force model. The framework is demonstrated through both numerical and experimental case studies. Results show that the approach accurately identifies both the system natural frequency and cutting force model.

42 ENGINEERING↗

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗

Prediction of softening kinetics and recrystallization texture in non-isothermally annealed bulged tubes using CPFEM and CA models

A hierarchically coupled cellular automata (CA) model, crystal plasticity finite element method (CPFEM), and thermal finite element (FE) model is developed to predict the softening kinetics of the bulged steel tube during non-isothermal annealing. Through the developed model, the kinetics of softening mechanisms including static recovery (SRV) and static recrystallization (SRX), as well as the recrystallization texture are predicted. Later, the Johnson-Mehl-Avrami-Kohnogorov (JMAK) model based on the predicted SRX data is developed to interpret the recrystallization behavior of the material. To perform this study, diverse experimental tests including tube hydroforming (THF), annealing, uniaxial tensile test, hardness test, as well as microstructure observations through optical microscopy and Electron Backscatter Diffraction (EBSD) tests on steel tube are performed. The obtained experimental data are utilized to calibrate and verify the implemented CPFEM model for simulation of THF process, thermal FE model for prediction of the local temperature over annealing time, and CA algorithm for modeling of the softening kinetics and texture evolution throughout the annealing process. The study shows that the predicted deformation characteristics, softening kinetics, recrystallization texture and temperature profile during non-isothermal annealing are in good agreement with experimental data. During the annealing process, a total of four stages for the kinetics of softening mechanisms is observed: No softening; SRV only; SRV dominant; and SRX dominant. During the progress of SRX, the behavior of recrystallization is abruptly changed, confirming that two different mechanisms are controlling the kinetics of transformation.

36 MATERIALS SCIENCE↗

High-Resolution Model Intercomparison Project phase 2 (HighResMIP2) towards CMIP7

Abstract. Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. Consequently, the role of enhanced horizontal resolution in improved process representation in all components of the climate system continues to be of great interest. Recent simulations suggest the possibility of significant changes in both large-scale aspects of the ocean and atmospheric circulations and in the regional responses to climate change, as well as improvements in representations of small-scale processes and extremes, when resolution is enhanced. The first phase of the High-Resolution Model Intercomparison Project (HighResMIP1) was successful at producing a baseline multi-model assessment of global simulations with model grid spacings of 25–50 km in the atmosphere and 10–25 km in the ocean, a significant increase when compared to models with standard resolutions on the order of 1° that are typically used as part of the Coupled Model Intercomparison Project (CMIP) experiments. In addition to over 250 peer-reviewed manuscripts using the published HighResMIP1 datasets, the results were widely cited in the Intergovernmental Panel on Climate Change report and were the basis of a variety of derived datasets, including tracked cyclones (both tropical and extratropical), river discharge, storm surge, and impact studies. There were also suggestions from the few ocean eddy-rich coupled simulations that aspects of climate variability and change might be significantly influenced by improved process representation in such models. The compromises that HighResMIP1 made should now be revisited, given the recent major advances in modelling and computing resources. Aspects that will be reconsidered include experimental design and simulation length, complexity, and resolution. In addition, larger ensemble sizes and a wider range of future scenarios would enhance the applicability of HighResMIP. Therefore, we propose the High-Resolution Model Intercomparison Project phase 2 (HighResMIP2) to improve and extend the previous work, to address new science questions, and to further advance our understanding of the role of horizontal resolution (and hence process representation) in state-of-the-art climate simulations. With further increases in high-performance computing resources and modelling advances, along with the ability to take full advantage of these computational resources, an enhanced investigation of the drivers and consequences of variability and change in both large- and synoptic-scale weather and climate is now possible. With the arrival of global cloud-resolving models (currently run for relatively short timescales), there is also an opportunity to improve links between such models and more traditional CMIP models, with HighResMIP providing a bridge to link understanding between these domains. HighResMIP also aims to link to other CMIP projects and international efforts such as the World Climate Research Program lighthouse activities and various digital twin initiatives. It also has the potential to be used as training and validation data for the fast-evolving machine learning climate models.

54 ENVIRONMENTAL SCIENCES↗

Earth-system-model evaluation of cloud and precipitation occurrence for supercooled and warm clouds over the Southern Ocean's Macquarie Island

Abstract. Over the remote Southern Ocean (SO), cloud feedbacks contribute substantially to Earth system model (ESM) radiative biases. The evolution of low Southern Ocean clouds (cloud-top heights < ∼ 3 km) is strongly modulated by precipitation and/or evaporation, which act as the primary sink of cloud condensate. Constraining precipitation processes in ESMs requires robust observations suitable for process-level evaluations. A year-long subset (April 2016–March 2017) of ground-based profiling instrumentation deployed during the Macquarie Island Cloud and Radiation Experiment (MICRE) field campaign (54.5∘ S, 158.9∘ E) combines a 95 GHz (W-band) Doppler cloud radar, two lidar ceilometers, and balloon-borne soundings to quantify the occurrence frequency of precipitation from the liquid-phase cloud base. Liquid-based clouds at Macquarie Island precipitate ∼ 70 % of the time, with deeper and colder clouds precipitating more frequently and at a higher intensity compared to thinner and warmer clouds. Supercooled cloud layers precipitate more readily than layers with cloud-top temperatures > 0 ∘C, regardless of the geometric thickness of the layer, and also evaporate more frequently. We further demonstrate an approach to employ these observational constraints for evaluation of a 9-year GISS-ModelE3 ESM simulation. Model output is processed through the Earth Model Column Collaboratory (EMC2) radar and lidar instrument simulator with the same instrument specifications as those deployed during MICRE, therefore accounting for instrument sensitivities and ensuring a coherent comparison. Relative to MICRE observations, the ESM produces a smaller cloud occurrence frequency, smaller precipitation occurrence frequency, and greater sub-cloud evaporation. The lower precipitation occurrence frequency by the ESM relative to MICRE contrasts with numerous studies that suggest a ubiquitous bias by ESMs to precipitate too frequently over the SO when compared with satellite-based observations, likely owing to sensitivity limitations of spaceborne instrumentation and different sampling methodologies for ground- versus space-based observations. Despite these deficiencies, the ESM reproduces the observed tendency for deeper and colder clouds to precipitate more frequently and at a higher intensity. The ESM also reproduces specific cloud regimes, including near-surface clouds that account for ∼ 25 % of liquid-based clouds during MICRE and optically thin, non-precipitating clouds that account for ∼ 27 % of clouds with bases higher than 250 m. We suggest that the demonstrated framework, which merges observations with appropriately constrained model output, is a valuable approach to evaluate processes responsible for cloud radiative feedbacks in ESMs.

54 ENVIRONMENTAL SCIENCES↗

Simulation process and data flow for a large system dynamics model

This paper documents the workflow and supporting technologies that a large system dynamics model, the biomass scenario model, employs to streamline the data preparation, simulation, quality control, and analysis process at the National Renewable Energy Laboratory. The workflow centers on automation of routine aspects of the flow of data between data stores, simulations, and visualizations. It enforces quality checks on data, reproducibility of computations, and traceability of results, while maintaining complete archives of modeling and analysis artifacts. The resulting frictionless simulation/analysis environment supports large-scale sensitivity analysis, interactive creation of ensembles of simulations, and rapid visualization-based exploration of simulation results.

09 BIOMASS FUELS↗

Multi-scale Coupled Processes Modeling of Fractures as Porous, Interfacial and Granular Systems from Rock Images with the Numerical Manifold Method

Abstract The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

58 GEOSCIENCES↗

Parallel simulation via SPPARKS of on-lattice kinetic and Metropolis Monte Carlo models for materials processing

Abstract SPPARKS is an open-source parallel simulation code for developing and running various kinds of on-lattice Monte Carlo models at the atomic or meso scales. It can be used to study the properties of solid-state materials as well as model their dynamic evolution during processing. The modular nature of the code allows new models and diagnostic computations to be added without modification to its core functionality, including its parallel algorithms. A variety of models for microstructural evolution (grain growth), solid-state diffusion, thin film deposition, and additive manufacturing (AM) processes are included in the code. SPPARKS can also be used to implement grid-based algorithms such as phase field or cellular automata models, to run either in tandem with a Monte Carlo method or independently. For very large systems such as AM applications, the Stitch I/O library is included, which enables only a small portion of a huge system to be resident in memory. In this paper we describe SPPARKS and its parallel algorithms and performance, explain how new Monte Carlo models can be added, and highlight a variety of applications which have been developed within the code.

36 MATERIALS SCIENCE↗

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↗

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↗