Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

MJO simulation in CMIP5 climate models: MJO skill metrics and process-oriented diagnosis

The Madden-Julian Oscillation (MJO) simulation diagnostics developed by MJO Working Group and the process-oriented MJO simulation diagnostics developed by MJO Task Force are applied to 37 Coupled Model Intercomparison Project phase 5 (CMIP5) models in order to assess model skill in representing amplitude, period, and coherent eastward propagation of the MJO, and to establish a link between MJO simulation skill and parameterized physical processes. Process-oriented diagnostics include the Relative Humidity Composite based on Precipitation (RHCP), Normalized Gross Moist Stability (NGMS), and the Greenhouse Enhancement Factor (GEF). Numerous scalar metrics are developed to quantify the results. Most CMIP5 models underestimate MJO amplitude, especially when outgoing longwave radiation (OLR) is used in the evaluation, and exhibit too fast phase speed while lacking coherence between eastward propagation of precipitation/convection and the wind field. The RHCP-metric, indicative of the sensitivity of simulated convection to low-level environmental moisture, and the NGMS-metric, indicative of the efficiency of a convective atmosphere for exporting moist static energy out of the column, show robust correlations with a large number of MJO skill metrics. The GEF-metric, indicative of the strength of the column-integrated longwave radiative heating due to cloud-radiation interaction, is also correlated with the MJO skill metrics, but shows relatively lower correlations compared to the RHCP- and NGMS-metrics. Our results suggest that modifications to processes associated with moisture-convection coupling and the gross moist stability might be the most fruitful for improving simulations of the MJO. Though the GEF-metric exhibits lower correlations with the MJO skill metrics, the longwave radiation feedback is highly relevant for simulating the weak precipitation anomaly regime that may be important for the establishment of shallow convection and the transition to deep convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Viscosity Measurements of the HLW APPS Glasses

The U.S. Department of Energy is responsible for building the Hanford Waste Treatment and Immobilization Plant (WTP) at the Hanford site in Washington to remediate 56 million gallons of radioactive waste historically stored in 177 underground tanks. The Office of River Protection has requested that the Savannah River National Laboratory (SRNL) contribute in areas of recognized capabilities and expertise for glass waste form development to support successful startup of the WTP. SRNL support of this work is defined in the Task Technical and Quality Assurance Plan.1 This report provides results from viscosity measurements on a series of simulated nuclear waste glasses designed and fabricated at Pacific Northwest National Laboratory (PNNL). The glasses were designated the high-level waste – Aspen Process Performance Simulation (HLW-APPS) study glasses. The data provided in this report is to be used in the development, validation, and implementation of enhanced property/composition models for nuclear waste glasses.

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

Fast Gaussian Process Estimation for Large-Scale In Situ Inference using Convolutional Neural Networks

Exascale computing will bring with it significant I/O limitations. One foreseeable consequence of such restrictions is that the user can save only a small fraction of complex simulation data to disk for subsequent analysis. An alternative is to fit statistical models to data in situ, that is, inside the simulation as it runs. This option requires extremely fast statistical estimation to avoid slowing down the simulation. Gaussian processes (GPs) have state-of-the-art predictive performance for modeling spatial data. However, standard estimation methods for GPs scale quite poorly to large data sets as parameter estimation requires inverting a covariance matrix to the size of the data set. In the presented work, we use a convolutional neural network (CNN) to predict the GP parameters for a spatial data set, from a simulation or otherwise, rather than optimize the parameters directly. Here, our presented case study models spatial data from E3SM, the Department of Energy’s Exascale climate model. The CNN is trained on synthetic data simulated from GP models with known parameters and then applied to data from the climate simulation. In the presented examples, the neural network scheme produces parameter estimates that compare well with standard methods such as maximum likelihood estimation in predictive performance but is obtained four orders of magnitude faster.

big data↗
DOE OSTI2022

Comptonization by reconnection plasmoids in black hole coronae II: Electron-ion plasma

ABSTRACT We perform 2D particle-in-cell simulations of magnetic reconnection in electron-ion plasmas subject to strong Compton cooling and calculate the X-ray spectra produced by this process. The simulations are performed for trans-relativistic reconnection with magnetization 1 ≤ σ ≤ 3 (defined as the ratio of magnetic tension to plasma rest-mass energy density), which is expected in the coronae of accretion discs around black holes. We find that magnetic dissipation proceeds with inefficient energy exchange between the heated ions and the Compton-cooled electrons. As a result, most electrons are kept at a low temperature in Compton equilibrium with radiation, and so thermal Comptonization cannot reach photon energies $\sim 100\,$ keV observed from accreting black holes. Nevertheless, magnetic reconnection efficiently generates $\sim 100\,$ keV photons because of mildly relativistic bulk motions of the plasmoid chain formed in the reconnection layer. Comptonization by the plasmoid motions dominates the radiative output and controls the peak of the radiation spectrum Epk. We find Epk ∼ 40 keV for σ = 1 and Epk ∼ 100 keV for σ = 3. In addition to the X-ray peak around 100 keV, the simulations show a non-thermal MeV tail emitted by a non-thermal electron population generated near X-points of the reconnection layer. The results are consistent with the typical hard state of accreting black holes. In particular, we find that the spectrum of Cygnus X-1 is well explained by electron-ion reconnection with σ ∼ 3.

Astronomy & Astrophysics↗
DOE OSTI2021

Impact Fracture and Fragmentation of Glass via the 3D Combined Finite-Discrete Element Method

A driving technical concern for the automobile industry is their assurance that developed windshield products meet Federal safety standards. Besides conducting innumerable glass breakage experiments, product developers also have the option of utilizing numerical approaches that can provide further insight into glass impact breakage, fracture, and fragmentation. The combined finite-discrete element method (FDEM) is one such tool and was used in this study to investigate 3D impact glass fracture processes. To enable this analysis, a generalized traction-separation model, which defines the constitutive relationship between the traction and separation in FDEM cohesive zone models, was introduced. The mechanical responses of a laminated glass and a glass plate under impact were then analyzed. For laminated glass, an impact fracture process was investigated and results were compared against corresponding experiments. Correspondingly, two glass plate impact fracture patterns, i.e., concentric fractures and radial fractures, were simulated. The results show that for both cases, FDEM simulated fracture processes and fracture patterns are in good agreement with the experimental observations. The work demonstrates that FDEM is an effective tool for modeling of fracture and fragmentation in glass.

42 ENGINEERING↗
DOE OSTI2024

High-Fidelity Accelerated Design of High-performance Electrochemical Systems

Large-scale electrification is vital to addressing the climate crisis, but several scientific and technological challenges remain to fully electrify both the chemical industry and transportation. In both of these areas, new electrochemical materials will be critical, but their development currently relies heavily on human-time-intensive experimental trial and error and computationally expensive first-principles, meso-scale and continuum simulations. To accelerate this process, our team has developed the AutoMat platform. AutoMat can accelerate development of new electrochemical materials along two avenues: first, automated input generation and management of simulations at multiple lengthscales as well as “handoff” of outputs from one lengthscale as inputs to the next; and second, replacement of the most computationally intensive simulation processes with machine-learned surrogate models. The crux of our team’s effort was not “reinventing the wheel” by developing entirely new techniques, but rather building a “superhighway” that allows existing state-of-the-art techniques to run faster and more smoothly than before. AutoMat can utilize tools spanning from first-principles quantum chemistry computations to automated robotic experimentation, and is driven by design space search techniques to reduce the number of iterations through the full simulation loop by rapidly targeting promising regions of design spaces such as single-atom alloy catalysts or blends of liquid electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2020

Countercurrent Actinide Lanthanide Separation Process (ALSEP) Demonstration Test with a Simulated PUREX Raffinate in Centrifugal Contactors on the Laboratory Scale

An Actinide Lanthanide Separation Process (ALSEP) for the separation of trivalent actinides (An(III)) from simulated raffinate solution was successfully demonstrated using a 32-stage 1 cm annular centrifugal contactor setup. The ALSEP solvent was composed of a mixture of 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEH[EHP]) and N,N,N',N'-tetra-(2-ethylhexyl)-diglycolamide (T2EHDGA) in n-dodecane. Flowsheet calculations and evaluation of the results were done using the Argonne’s Model for Universal Solvent Extraction (AMUSE) code using single-stage distribution data. The co-extraction of Zr(IV) and Pd(II) was prevented using CDTA (trans-1,2-diaminocyclohexane-N,N,N',N'-tetraacetic acid) as a masking agent in the feed. For the scrubbing of co-extracted Mo; citrate-buffered acetohydroxamic acid was used. The separation of An(III) from the trivalent lanthanides (Ln(III)) was achieved using citrate-buffered diethylene-triamine-N,N,N',N",N"-pentaacetic acid (DTPA), and Ln(III) were efficiently back extracted using N,N,N',N'-tetraethyl-diglycolamide (TEDGA). A clean An(III) product was obtained with a recovery of 95% americium and curium. The Ln(III) were efficiently stripped; but the Ln(III) product contained 5% of the co-stripped An(III). The carryover of Am and Cm into the Ln(III) product is attributed to too few actinide stripping stages, which was constrained by the number of centrifugal contactors available. Improved separation would be achieved by increasing the number of An strip stages. The heavier lanthanides (Pr, Nd, Sm, Eu, and Gd) and yttrium were mainly routed to the Ln product, whereas the lighter lanthanides (La and Ce) were mostly routed to the raffinate.

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

Simulation of 12-bed vacuum pressure-swing adsorption for hydrogen separation from methanol-steam reforming off-gas

This study focuses on analysis of a 12-bed vacuum pressure-swing adsorption (VPSA) process capable of purifying hydrogen from a ternary mixture (H 2 /CO 2 /CO 75/24/1 mol%) derived from methanol-steam reforming. The process produces 9 kmol H 2 /h with less than 2 ppm and 0.2 ppm of CO2 and CO, respectively, to supply a polymer electrolyte membrane fuel cell. The process model is developed in Aspen Adsorption® using the “uni-bed” approach. A parametric study of H 2 purity and recovery with respect to adsorption pressure, adsorbent height, activated carbon:zeolite ratio, feed composition, and number of beds is performed. Results show 12-bed VPSA can meet the H 2 purity goals, with H 2 recovery as high as 75.75%. Adsorption occurs at 7 bar, the column height is 1.2 m, and the adsorbent ratio is 70%:30%. Furthermore, a 4-bed VPSA can achieve the same purity goals as the 12-bed process, but H 2 recovery decreases to 61.34%.

03 NATURAL GAS↗
DOE OSTI2019

Dissociation Mechanism of Processive Cellulases Explored Through Molecular Simulation

Cellulose is the most abundant biopolymer on the planet, with significant potential as a feedstock for bioderived products. In nature, fungi and bacteria secrete cellulases to hydrolyze cellulose into smaller carbohydrates for incorporation into their own metabolic processes. For effective biomass utilization at industrial scale, these cellulases are coupled together to produce a clean sugar stream. One such cellulase is Cel7A from Trichoderma reesei, which is used extensively in biotechnological applications. However, dissociation of Cel7A from a bound cellulose strand is rate-limiting, reducing catalytic efficiency and yields within industrial applications. To explore this unknown dissociation mechanism and devise potential enzyme engineering strategies to mitigate slow dissociation, we conduct Hamiltonian replica exchange molecular dynamics simulations. Two postulated dissociation mechanisms were tested ('dethreading' of the bound cellulose strand from the binding tunnel and 'clamshell'-like opening of the binding tunnel and recrystallization of the bound cellulose). For each mechanism, free energy profiles were determined together with estimated diffusion coefficients. The kinetics for both mechanisms were then estimated and compared against experimentally measured dissociation rates, which suggest that the dethreading mechanism ought to be preferred. These findings are further supported by dissociation rates measured on engineered mutants incapable of undergoing the clamshell mechanism.

09 BIOMASS FUELS↗
DOE OSTI2025

Electronic effects in radiation-induced collision cascades in nickel

The accurate treatment of electronic effects in multi-million-atom simulations of radiation-induced collision cascades is crucial for reliable predictions of primary radiation damage. In this work, we explore the fidelity of a recently developed two-temperature molecular dynamics model implementing an electron density-dependent coupling of electronic and atomic subsystems for cascade simulations in nickel. We show that the parameter-free model realistically captures the instantaneous energy losses during all stages of the highly nonequilibrium cascade process. Our simulations predict two distinct coupling regimes, corresponding to the rapid energy losses through electronic stopping in the early stages of the cascade and to the slow equilibration through the electron-phonon coupling mechanism in the later stages, without the use of separate models or coupling terms. The intermediate stage of the cascade dynamics displays a complex energy transfer between the subsystems, which cannot be interpreted by comparison to either electronic stopping or electron-phonon coupling theories. We therefore compare the predicted atomic mixing, which is sensitive to the energy losses during the intermediate cascade stage, with experimental ion beam mixing measurements. We find good agreement with the experiments, validating the coupling model for the intermediate stage of the cascade. Predictions of final defect numbers and cluster sizes are found in line with the results from conventional electronic stopping-based methods, while significantly reducing the theoretical uncertainty in the outcomes of conventional models stemming from arbitrary choices of thresholds for different coupling terms. Our results represent a notable improvement in cascade damage predictions in nickel, providing validation of the electron density-dependent coupling model for radiation damage simulations in general. The results lead us to propose an interpretation of the electronic energy losses in the intermediate regime of velocities, where we find an effectively nonlinear dissipation.

Crystal defects↗
DOE OSTI2023

E3SMv1 simulation output from sensitivity study on aerosol process coupling

This dataset contains annual mean and 6-hourly output data from a set of E3SMv1 simulations that demonstrates the impact of numerical process coupling on the simulated dust life cycle. The simulations are documented and analyzed in the attached manuscript (pdf file). Publishing the manuscript in a peer-reviewed journal will require making the model output publicly available. We plan to use DataHub for that. The data files to be uploaded are currently saved on Compy under /compyfs/wanh895/scidac_share/cflx_2023/

Wan, Hui↗
DOE OSTI2016

Validating the simulation of large-scale parallel applications using statistical characteristics

Simulation is a widely adopted method to analyze and predict the performance of large-scale parallel applications. Validating the hardware model is highly important for complex simulations with a large number of parameters. Common practice involves calculating the percent error between the projected and the real execution time of a benchmark program. However, in a high-dimensional parameter space, this coarse-grained approach often suffers from parameter insensitivity, which may not be known a priori. Moreover, the traditional approach cannot be applied to the validation of software models, such as application skeletons used in online simulations. In this work, we present a methodology and a toolset for validating both hardware and software models by quantitatively comparing fine-grained statistical characteristics obtained from execution traces. Although statistical information has been used in tasks like performance optimization, this is the first attempt to apply it to simulation validation. Lastly, our experimental results show that the proposed evaluation approach offers significant improvement in fidelity when compared to evaluation using total execution time, and the proposed metrics serve as reliable criteria that progress toward automating the simulation tuning process.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in situ and remote sensing data sets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

Cheng, Yanyan↗
DOE OSTI2017

Mixing and trapping of dissolved CO 2 in deep geologic formations with shale layers

For dissolution trapping, the spatial variability of the geologic properties of naturally complex storage formations can significantly impact flow patterns and storage mechanisms of dissolved CO 2 . The significance of diffusive mixing that occurs in low permeability layers embedded between relatively higher permeability materials was highlighted by Agartan et al. (2015) using a highly controlled laboratory experimental study on trapping of dissolved CO 2 in multilayered systems. In this paper, we present a numerical modeling study on the impacts of low permeability layers on flow and storage of dissolved CO 2 in realistic field-scale settings. The simulator of variable-density flow used in this study was first verified using the experimental data in Agartan et al. (2015) to capture the observed processes. The simulator was then applied to a synthetic, field-scale multilayered system, with 19 sensitivity cases having variable permeability and thickness of the shale layers as well as the source strength and geometry of the source zone of dissolved CO 2 . Simulation results showed that the presence of continuous shale layers in the storage system disrupts the convective mixing by enhancing lateral spreading of dissolved CO 2 in sandstone layers and retarding the vertical mixing of dissolved CO 2 . The effectiveness of trapping of dissolved CO 2 depends on the physical properties of the shale layers and configurations of the source zone. The comparison to homogeneous cases with effective vertical permeability shows that it is important to capture these continuous thin shale layers in a storage formation and include them in the models to enhance dissolution trapping.

58 GEOSCIENCES↗
DOE OSTI2025

Numerical simulations of three-dimensional ion crystal dynamics in a Penning trap using the fast multipole method

We simulate the dynamics, including laser cooling, of three-dimensional (3-D) ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions’ equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions’ absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3-D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3-D ion crystals could be well suited as platforms for future quantum science experiments.

Zaris, John (ORCID:0009000196476323)↗
DOE OSTI2021

Intelligent Agents for the Optimization of Atomic Layer Deposition

Atomic layer deposition (ALD) is a highly controllable thin film synthesis approach with applications in computing, energy, and separations. The flexibility of ALD means that it can access a massive chemical catalogue; however, this chemical and process diversity results in significant challenges in determining processing parameters that result in stable and uniform film growth with minimal precursor consumption. In situ measurements of the ALD growth per cycle (GPC) can accelerate process development but it still requires expert intuition and time-consuming trial and error to identify acceptable processing parameters. This procedure is made more difficult by the presence of experimental noise in the GPC values and the complexity of ALD surface chemistries. Here, a need exists for efficient optimization approaches capable of autonomously determining processing conditions resulting in optimal ALD film growth. In this work, we present the development of three optimization strategies and compare their performance in optimizing four simulated ALD processes. Furthermore, the effect of noise in the GPC measurements on optimization convergence is studied.

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

Analysis of a Disturbance Event with Inverter-Based Resources Using EMT Simulations

Increasing penetration of inverter-based resources (IBRs) necessitates newer methods of planning and analysis of disturbances. The existing phasor-domain transient stability (TS) analysis may not capture the dynamics of IBRs during fault events. Here, in this paper, electromagnetic transient (EMT) simulations using high-fidelity detailed model of power grid and one of the affected photovoltaic (PV) plants during the Angeles Forest disturbance in 2018 are performed. In these simulations, the processes to develop EMT models of power grid from traditional phasor-domain TS data and PV plant from collected data are described. Thereafter, using these simulations, the response of the PV plant during the fault event in 2018 is replicated and a sensitivity analysis is performed. The sensitivity analysis consists of making changes to the components within the PV plant and in the power grid to evaluate the impact they have on the response observed by the PV plant during the fault event. This analysis provides an understanding of the components that impact the operation of a PV plant during fault events and provide guidance to system planners on the studies that need to be performed to maintain a reliable power grid as new IBR plants are integrated.

42 ENGINEERING↗
DOE OSTI2025

2025 Workshop on Envisioning Frontiers in AI and Computing for Biological Research: Position Papers

This workshop aims to identify key research directions for transforming biology using artificial intelligence (AI), machine learning (ML) and computational methods to facilitate the discovery of new behaviors, mechanisms, and designs of biological processes relevant to DOE missions, underpinning a broader U.S. bioeconomy. By developing novel AI/ML technologies to analyze and interpret complex biological data, researchers can organize and simulate biological processes at various scales as well as advance predictive understanding and manipulation of biological systems. This integration of computation, experimentation, and next-generation experimental technologies can lead to discoveries in new biological behaviors and mechanisms relevant to DOE missions. The focus is on how advanced computational and mathematical methods can impact this mission by exploring digital twins, foundation models, automated laboratory experiments, modeling of complex living systems, and data-driven approaches for the biodesign of plants and microbial systems. While data management is important, it is not the primary focus of this workshop, which will assess the current state, trends, and AI/ML challenges at the interface between biology and computational science to identify opportunities for high-impact research at their intersection. The goal is to define research needs and opportunities that align with biological sciences, computational sciences, and applied mathematics research.

59 BASIC BIOLOGICAL SCIENCES↗