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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A phase field model for dislocations in hexagonal close packed crystals

In this work, an application of a phase field formulation suitable for modeling the motion of individual partial and full dislocations in hexagonal close packed (HCP) crystals is presented. The formulation incorporates periodic potentials for glide on the distinct HCP slip systems, which are informed here by density functional theory (DFT). The model is applied to simulate the dissociation process starting from an unstable perfect dislocation and ending at its final equilibrium structure for different slip planes and in different HCP metals. The structural characteristics that are predicted for these dislocations include the partial Burgers vectors, dissociation distances, core widths of the partials, and any asymmetries in these quantities. Mg is selected as one of the model materials since its dislocations are the most well studied and it is nearly elastically isotropic. For Mg, it is shown that the predictions for dissociation distances agree with those reported previously by atomic-scale calculations, including density functional theory, for dislocations on the basal < a >, prismatic < a >, and pyramidal type II slip systems. Furthermore, phase field model results are also presented for dislocations in Ti and Zr, which we find develop distinctively different equilibrium structures than Mg.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Using CASE to Exploit Process Modeling in Technology Transfer

A successful business will be one that has processes in place to run that business. Creating processes, reengineering processes, and continually improving processes can be accomplished through extensive modeling. Casewise(R) Corporate Modeler(TM) CASE is a computer aided software engineering tool that will enable the Technology Transfer Department (TT) at NASA Marshall Space Flight Center (MSFC) to capture these abilities. After successful implementation of CASE, it could then go on to be applied in other departments at MSFC and other centers at NASA. The success of a business process is dependent upon the players working as a team and continuously improving the process. A good process fosters customer satisfaction as well as internal satisfaction in the organizational infrastructure. CASE provides a method for business process success through functions consisting of systems and processes business models; specialized diagrams; matrix management; simulation; report generation and publishing; and, linking, importing, and exporting documents and files. The software has an underlying repository or database to support these functions. The Casewise. manual informs us that dynamics modeling is a technique used in business design and analysis. Feedback is used as a tool for the end users and generates different ways of dealing with the process. Feedback on this project resulted from collection of issues through a systems analyst interface approach of interviews with process coordinators and Technical Points of Contact (TPOCs).

Renz-Olar, Cheryl↗
DOE OSTI2016

Fragmentation Prediction for a Brittle Ceramic: a Two-Scale Model Approach

Glass and oxide-ceramic radioactive material forms (e.g., PuO{sub 2}) are susceptible to brittle fracture during handling and transportation. Fracture of these materials may result in the generation of fragments small enough to become airborne, and thus easily dispersed and inhaled. Entrainment and dispersal of particles with dimensions less than 10 μm is generally anticipated under normal aerodynamic conditions. The assessment of the risk to the public and the environment from these materials due to events resulting in their fracture is based on the determination of the source term for environmental release. This source term is typically quantified by an airborne release fraction (ARF) and a respirable fraction (RF) of material generated by a given event. While the recommended empirical bounding values prove satisfactory to comply with the regulatory process for handling and transportation of these materials, they are based on test data for only a few materials subjected to a limited range of loading conditions. Because the final state of fragmentation (in terms of the size and distribution of fragments) is strongly dependent not only on the properties of the material and the energetic loading conditions that lead to breakup, but also on the many intrinsic length scales associated with fragmentation, it is difficult to precisely predict fragmentation characteristics for other brittle solid wastes forms subjected to a wider range of loading conditions based on the test data alone. The purpose of this manuscript is to present a two-scale modeling approach that may be used to simulate the dynamic fragmentation of brittle solid radioactive material forms subjected to general loading and boundary conditions. The method provides the resulting fragment characteristics (average size and size distribution) across the entire fragment range of interest. At the macro-scale, the approach uses a gradient damage elastic continuum mechanics material model within a finite element framework to simulate dynamic fracture under general loading and boundary conditions. Results from the macroscopic calculations are passed to a one-dimensional (1-D) lower length scale model that incorporates the relevant fine scale effects to determine the micro-scale fragmentation. This 1-D micro-scale fragmentation model simulates the internal crack nucleation and opening processes to predict the average fragment size and the fragment size distribution as a function of the material properties, loading conditions, and microscopic characteristic length scale of the material. Results from the lower length scale model are then combined with the macro-scale model fragmentation information to determine the resulting fragment size distribution spanning both length scales. This two-scale model has been exercised on a representative problem drawn from experimental studies, consisting of an impact of a falling weight onto a brittle ceramic. (authors)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Leaf-scale carbon assimilation simulated by the multi-assumption architecture & testbed for process sensitivity analysis

Raw and processed output from the Multi-Assumption Architecture & Testbed (MAAT, https://github.com/walkeranthonyp/MAAT, tag: v1.2.1_walker2020; Walker et al., 2018). The dataset contains leaf-scale carbon assimilation data from three MAAT simulations and their post-processing and plotting scripts. Simulation 1, labelled "FvC_4proc_proc," is a four-process, process sensitivity analysis (Dai et al., 2017) of carbon assimilation based on the Farquhar et al. (1980) photosynthesis model, and its variant proposed by Collatz et al. (1991). Simulation 2, labelled "FvC_4proc_salt," is a 14 parameter sensitivity analysis (Saltelli et al., 2010) of carbon assimilation based on the Farquhar and Collatz photosynthesis models. Simulation 3, labelled "FvC," is a factorial combination of the alternative process representations and parameter values presented in those original papers across a factorial combination of incident photosynthetically active radiation (50-2000 µmol m-2 s-1, in increments of 50) and atmospheric CO2 concentration (50-1500 µmol m-2 s-1, in increments of 50). Additional scripts/functions for post-processing and plotting are part of the MAAT repository. These data were generated to evaluate and compare the Farquhar and Collatz models of leaf-scale photosynthesis across a range of CO2 and incident light conditions. Analysis of results can be found in Walker et al. (2020). This dataset is comprised of R scripts to process data and produce figures; .RDS files, a binary format specific to R and read with the 'readRDS' function, to store raw MAAT output; R and bash scripts (.bs) that help to run a MAAT simulation, .XML files that are output from MAAT describing a run, and .csv files which contain the calculated sensitivity indexes published in Walker et al. (2020).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Trapped-ion quantum simulation of electron transfer models with tunable dissipation

Electron transfer is at the heart of many fundamental physical, chemical, and biochemical processes essential for life. The exact simulation of these reactions is often hindered by the large number of degrees of freedom and by the essential role of quantum effects. Here, we experimentally simulate a paradigmatic model of molecular electron transfer using a multispecies trapped-ion crystal, where the donor-acceptor gap, the electronic and vibronic couplings, and the bath relaxation dynamics can all be controlled independently. By manipulating both the ground-state and optical qubits, we observe the real-time dynamics of the spin excitation, measuring the transfer rate in several regimes of adiabaticity and relaxation dynamics. Our results provide a testing ground for increasingly rich models of molecular excitation transfer processes that are relevant for molecular electronics and light-harvesting systems.

Science & Technology - Other Topics↗
NASA NTRS2023

Importance of Internal Variability for Climate Model Assessment

Benchmarking climate model simulations against observations of the climate is core to the process of building realistic climate models and developing accurate future projections. However, in many cases, models do not match historical observations, particularly on regional scales. If there is a mismatch between modeled and observed climate features, should we necessarily conclude that our models are deficient? Using several illustrative examples, we emphasize that internal variability can easily lead to marked differences between the basic features of the model and observed climate, even when decades of model and observed data are available. This can appear as an apparent failure of models to capture regional trends or changes in global teleconnections, or simulation of extreme events. Despite a large body of literature on the impact of internal variability on climate, this acknowledgment has not yet penetrated many model evaluation activities, particularly for regional climate. We emphasize that using a single or small ensemble of simulations to conclude that a climate model is in error can lead to premature conclusions on model fidelity. A large ensemble of multidecadal simulations is therefore needed to properly sample internal climate variability in order to robustly identify model deficiencies and convincingly demonstrate progress between generations of climate models.

Climate and Earth system modelling↗
DOE OSTI2023

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of 5 subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the 5 pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

Dim, Odera↗
DOE OSTI2023

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of five subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the five pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2019

The Necessity for Iteration in the Application of Numerical Simulation to EGS: Examples from the EGS Collab Test Bed 1

The United States Department of Energy, Geothermal Technologies Office (GTO) is funding a collaborative investigation of enhanced geothermal systems (EGS) processes at the meso-scale. This study, referred to as the EGS Collab project, is a unique opportunity for scientists and engineers to investigate the creation of fracture networks and circulation of fluids across those networks under in-situ stress conditions. The EGS Collab project is envisioned to comprise three experiments and the site for the first experiment is on the 4850 Level in phyllite of the Precambrian Poorman formation, at the Sanford Underground Research Facility, located at the former Homestake Gold Mine, in Lead, South Dakota. Principal objectives of the project are to develop a number of intermediate-scale field sites and to conduct well-controlled in situ experiments focused on rock fracture behavior and permeability enhancement. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO's flagship EGS research effort. As noted, outcomes from the EGS Collab project experiments will serve as benchmarks for computer code verification, but numerical simulation additionally plays an essential role in designing these meso-scale experiments. This paper reviews specific numerical simulations supporting the design of experiments within Test Bed 1, a volume of phyllite rock under in-situ stress conditions off the western side of the West Access Drift on the 4850 Level, near Governor's Corner. Numerical simulations were executed prior to the start of hydraulic stimulation activities within Test Bed 1 following standard practices, using best estimates of principal stress conditions, thermal conditions, and the rock petrophysical properties, including geomechanical properties. These simulations indicated notching of the borehole would promote the initiation of transverse fractures, seismic magnitudes during the hydraulic stimulation would be below 0.1 magnitude on the Richter scale, a rock thermal conductivity of 5.0 W/m K yielded agreement with kISMET borehole temperature logs, back pressure on the production borehole would increase circulation across the test bed, and the production borehole would serve to halt propagation of the hydraulic fracture to the drift. Pre-existing natural fractures, heterogeneities in the rock properties, monitoring boreholes, and overlooked mine elements have prompted a second look at numerically modeling stimulation, fluid circulation, tracer migration, and thermal breakthrough. Numerical simulation is an invaluable tool for providing insight and understanding to complex physical processes. The success of simulations, however, often depends on including all of the salient features of the system in the founding conceptual model. This paper takes a retrospective look at examples where the conceptual model and simulation was sufficient to provide accurate forecasts and those where elements were missing, necessitating rethinking of the simulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

The Necessity for Iteration in the Application of Numerical Simulation to EGS: Examples from the EGS Collab Test Bed 1

The United States Department of Energy, Geothermal Technologies Office (GTO) is funding a collaborative investigation of enhanced geothermal systems (EGS) processes at the meso-scale. This study, referred to as the EGS Collab project, is a unique opportunity for scientists and engineers to investigate the creation of fracture networks and circulation of fluids across those networks under in-situ stress conditions. The EGS Collab project is envisioned to comprise three experiments and the site for the first experiment is on the 4850 Level in phyllite of the Precambrian Poorman formation, at the Sanford Underground Research Facility, located at the former Homestake Gold Mine, in Lead, South Dakota. Principal objectives of the project are to develop a number of intermediate-scale field sites and to conduct well-controlled in situ experiments focused on rock fracture behavior and permeability enhancement. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. As noted, outcomes from the EGS Collab project experiments will serve as benchmarks for computer code verification, but numerical simulation additionally plays an essential role in designing these meso-scale experiments. This paper reviews specific numerical simulations supporting the design of experiments within Test Bed 1, a volume of phyllite rock under in-situ stress conditions off the western side of the West Access Drift on the 4850 Level, near Governor’s Corner. Numerical simulations were executed prior to the start of hydraulic stimulation activities within Test Bed 1 following standard practices, using best estimates of principal stress conditions, thermal conditions, and the rock petrophysical properties, including geomechanical properties. These simulations indicated notching of the borehole would promote the initiation of transverse fractures, seismic magnitudes during the hydraulic stimulation would be below 0.1 magnitude on the Richter scale, a rock thermal conductivity of 5.0 W/m K yielded agreement with kISMET borehole temperature logs, back pressure on the production borehole would increase circulation across the test bed, and the production borehole would serve to halt propagation of the hydraulic fracture to the drift. Pre-existing natural fractures, heterogeneities in the rock properties, monitoring boreholes, and overlooked mine elements have prompted a second look at numerically modeling stimulation, fluid circulation, tracer migration, and thermal breakthrough. Numerical simulation is an invaluable tool for providing insight and understanding to complex physical processes. The success of simulations, however, often depends on including all of the salient features of the system in the founding conceptual model. This paper takes a retrospective look at examples where the conceptual model and simulation was sufficient to provide accurate forecasts and those where elements were missing, necessitating rethinking of the simulation.

numerical simulation, EGS Collab, enhanced geother↗
NASA NTRS2000

Modeling and Analysis of the Reverse Water Gas Shift Process for In-Situ Propellant Production

This report focuses on the development of mathematical models and simulation tools developed for the Reverse Water Gas Shift (RWGS) process. This process is a candidate technology for oxygen production on Mars under the In-Situ Propellant Production (ISPP) project. An analysis of the RWGS process was performed using a material balance for the system. The material balance is very complex due to the downstream separations and subsequent recycle inherent with the process. A numerical simulation was developed for the RWGS process to provide a tool for analysis and optimization of experimental hardware, which will be constructed later this year at Kennedy Space Center (KSC). Attempts to solve the material balance for the system, which can be defined by 27 nonlinear equations, initially failed. A convergence scheme was developed which led to successful solution of the material balance, however the simplified equations used for the gas separation membrane were found insufficient. Additional more rigorous models were successfully developed and solved for the membrane separation. Sample results from these models are included in this report, with recommendations for experimental work needed for model validation.

Whitlow, Jonathan E.↗
DOE OSTI2020

Assessment of turbulent combustion models for simulating pre-chamber ignition in a natural gas engine

Lean combustion in an internal combustion engine is a promising strategy to increase thermal efficiency by leveraging a more favorable specific heat ratio of the fresh mixture and simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition events and slow flame propagation at fuel-lean condition lead to high cycle-to-cycle variability and hence limit the high-efficiency engine operating range. Pre-chamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling perspective, to date, still the science base and available simulation tools are inadequate for understanding and predicting the combustion processes in pre-chamber ignited engines. In this paper, conceptually different RANS combustion models widely adopted in the engine modeling community were used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a pre-chamber spark-ignition system. A flamelet-based turbulent combustion model, i.e., G-equation, and a multi-zone well-stirred reactor model were employed for the multi-dimensional study. Simulation results were compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two formulations respectively.

combustion model evaluation↗
NASA NTRSDate not supplied

Developing Urban Air Mobility Vehicle Models to Support Air Traffic Management Concept Development

To support Urban Air Mobility (UAM) research efforts at NASA, the Airspace Target Generator (ATG) software used in the FutureFlight Central (FFC) air traffic control tower simulator is undergoing updates to support physics-based UAM vehicle models. A process was developed to integrate UAM vertical takeoff and landing (VTOL) aircraft into the fixed-wing ATG modeling environment without significant change to the underlying equations of motion and vehicle model database. The VTOL aircraft models were converted from a six degrees-of-freedom (6-DOF) representation into a four degrees-of-freedom (4-DOF) representation for integration within ATG. Three vehicle designs from the NASA Revolutionary Vertical-Lift Technologies (RVLT) project were selected: a lift-plus-cruise (LPC) aircraft model and quadrotor, electric-powered (QEP) 1-seater and 6-seater models. With the LPC model comprised of a nonlinear force and moment build-up, and the QEP models comprised of linearized stability derivatives, two separate processes were developed to convert the lift, drag, and propulsion characteristics of each model into the ATG model database. Key aircraft performance characteristics including climb, cruise, and descent performance were preserved during the conversion process. Because ATG simulates fixed-wing aircraft through ground taxi and takeoff to approach and landing, acceleration command algorithms were developed to model the vertical takeoff and vertical landing phase of UAM operations. A strategy was then developed to transition the aircraft model to- and from- the new control mode.

urban air mobility↗
DOE OSTI2021

Multifidelity computing for coupling full and reduced order models

Hybrid physics-machine learning models are increasingly being used in simulations of transport processes. Many complex multiphysics systems relevant to scientific and engineering applications include multiple spatiotemporal scales and comprise a multifidelity problem sharing an interface between various formulations or heterogeneous computational entities. To this end, we present a robust hybrid analysis and modeling approach combining a physics-based full order model (FOM) and a data-driven reduced order model (ROM) to form the building blocks of an integrated approach among mixed fidelity descriptions toward predictive digital twin technologies. At the interface, we introduce a long short-term memory network to bridge these high and low-fidelity models in various forms of interfacial error correction or prolongation. The proposed interface learning approaches are tested as a new way to address ROM-FOM coupling problems solving nonlinear advection-diffusion flow situations with a bifidelity setup that captures the essence of a broad class of transport processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Final Technical Report for Grant No. DE-FC02-01ER54648, SciDAC Center for Simulation of Wave-Plasma Interactions

Harnessing the energy that is released in fusion reactions would provide a safe and abundant source of power to meet the growing energy needs of the world population. The next step toward the development of fusion as a practical energy source is the construction of ITER, a device capable of producing and controlling the high performance plasma required for self-sustaining fusion reactions, or “burning” plasma. The input power required to drive the ITER plasma into the burning regime will be supplied primarily with a combination of external power from radio frequency waves in the ion cyclotron range of frequencies and energetic ions from neutral beam injection sources, in addition to internally-generated Ohmic heating from the induced plasma current that also serves to create the magnetic equilibrium for the discharge. Since the success of the ITER project depends critically on the ability to create and maintain burning plasma conditions, it is absolutely necessary to have physics-based models that can accurately simulate the RF processes that affect the dynamical evolution of the ITER discharge. The Center for Simulation of Wave-Plasma Interactions (CSWPI) is a multi-institutional collaboration that has conducted ongoing research aimed at developing: (1) Coupled core-to-edge simulations that will lead to an increased understanding of parasitic losses of the applied RF power in the boundary plasma between the RF antenna and the core plasma; (2) Development of models for core interactions of RF waves with energetic electrons and ions (including fusion alpha particles and fast neutral beam ions) that include a more accurate representation of the particle dynamics in the combined equilibrium and wave fields; and (3) Development of improved algorithms that will take advantage of massively parallel computing platforms at the petascale level and beyond to achieve the needed physics, resolution, and/or statistics to address these issues. This research reported in this final technical report is focused on further verification and validation of simulation models for interactions of RF waves with energetic particles in the core with a long term goal of assessing their readiness for inclusion in integrated modeling frameworks and utilizing the world’s largest supercomputers to quantitatively predict the detailed electric field structure in the vicinity of the launching structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Dark fluxes from electromagnetic cascades

We study dark sector production in electromagnetic (EM) cascades. This problem requires accurate simulations of Standard Model (SM) and dark sector processes, both of which impact angular and energy distributions of emitted particles that ultimately determine flux predictions in a downstream detector. We describe the minimal set of QED processes which must be included to faithfully reproduce a SM cascade, and identify a universal algorithm to generate a dark sector flux given a Monte-Carlo simulation of a SM shower. We provide a new tool, "Image missing", which simulates EM cascades with associated dark vector production, and compare it against existing literature and “off the shelf” tools. The signal predictions at downstream detectors can strongly depend on the nontrivial interplay (and modelling) of SM and dark sector processes, in particular multiple Coulomb scattering and positron annihilation. We comment on potential impacts of these effects for realistic experimental setups.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗