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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Crack propagation in silica from reactive classical molecular dynamics simulations

Mechanistic insight into the process of crack growth can be obtained through molecular dynamics (MD) simulations. In this investigation of fracture propagation, a slit crack was introduced into an atomistic amorphous silica model and mode I stress was applied through far–field loading until the crack propagates. Atomic displacements and forces and an Irving–Kirkwood method with a Lagrangian kernel estimator were used to calculate the J –integral of classical fracture mechanics around the crack tip. The resulting fracture toughness (K IC ), 0.76 ± 0.16 MPa√m, agrees with experimental values. In addition, the stress fields and dissipation energies around the slit crack indicate the development of an inelastic region ~30Å in diameter. This is one of the first reports of KIC values obtained from up–scaled atomic–level energies and stresses through the J –integral. In conclusion, the application of the ReaxFF classical MD force field in this study provides the basis for future research into crack growth in multicomponent oxides in a variety of environmental conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A Review of Lattice-Boltzmann Models Coupled with Geochemical Modeling Applied for Simulation of Advanced Waterflooding and Enhanced Oil Recovery Processes

To maintain economic profit and improve the oil production efficiency after the primary and secondary production phase, advanced waterflooding techniques such as low salinity waterflooding in carbonate reservoirs have been investigated in numerical simulations, laboratory experiments, and field pilot tests. Multiple underlying mechanisms have been proposed based on these studies, and they are still under debate. Various numerical modeling approaches are introduced, but there exists a lack of a pore-scale comprehensive modeling scheme to fully understand the processes. Lattice-Boltzmann method (LBM) is a type of numerical fluid flow modeling technique that shows capabilities and flexibilities in modeling pore-scale fluid flow to integrate physical–chemical processes within complex structures. The intrinsic feature of LBM makes it a promising framework for simulating advanced waterflooding due to its flexibility, accuracy, and parallel efficiency. LBM works either by itself for solving reactive transport problems or by coupling with a third-party reaction solver. This review mainly introduces the LBM fluid flow and reactive transport capabilities and the concept and modeling approaches to simulate advanced waterflooding techniques. Meanwhile, an evaluation of the coupled LBM models for enhanced oil recovery (EOR) simulations is discussed with future research challenges and directions concluded.

02 PETROLEUM↗

DOE OSTI2025

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Combining Experiments and Simulations of Extraction Kinetics and Thermodynamics in Advanced Separation Processes for Used Nuclear Fuel

This 3-year project was a collaboration between University of California Irvine (UC Irvine), Pacific Northwest National Laboratory (PNNL), Idaho National Laboratory (INL), Argonne National Laboratory (ANL) and with an international collaborator at ForschungZentrum Jülich (FZJ). The project was led from UC Irvine under the direction of Profs. Mikael Nilsson and Hung Nguyen. The leads at PNNL, INL, ANL and FZJ were Dr. Liem Dang, Dr. Peter Zalupski, Dr. Nathaniel Hoyt and Dr. Giuseppe Modolo, respectively. Involved in this project at UC Irvine were three full time PhD graduate students, Tro Babikian, Ted Yoo, and Quynh Vo, and one MS student, Alba Font Bosch. The overall objective of this project was to study how the kinetics and thermodynamics of metal ion extraction can be described by molecular dynamic (MD) simulations and how the simulations can be validated by experimental data. Furthermore, the project includes the applied separation by testing the extraction systems in a single stage annular centrifugal contactor and coupling the experimental data with computational fluid dynamic (CFD) simulations. Specific objectives of the proposed research were: Study and establish a rigorous connection between MD simulations based on polarizable force fields and extraction thermodynamic and kinetic data. Compare and validate CFD simulations of extraction processes for An/Ln separation using different sizes (and types) of annular centrifugal contactors. Provide a theoretical/simulation and experimental base for scale-up of batch-wise extraction to continuous contactors. We approached objective 1 and 2 in parallel. For objective 1 we started by studying a well established extraction system with a relatively simple extraction mechanism, namely tributyl phosphate. What we found was that well optimized simulations can inform experiments and new information on TBP behavior was presented in this project, as well be discussed below. The second objective proved a larger challenge and most of the efforts were devoted to experimental studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2017

On the kinetics of transgranular particle embrittlement during simulated carburizing in steel containing grain-refining additions of aluminum and niobium plus aluminum

Here, the kinetics of toughness degradation resulting from transgranular particle embrittlement are evaluated as a function of composition and processing history for simulated carburizing operations in air-melt steel containing grain-refining additions of aluminum and aluminum plus niobium. The kinetics of particle embrittlement are inherently linked to the ripening of AlN precipitates after extended austenitization in steel containing carbon contents representative of both the case and core of a carburized component. Embrittlement in steel containing AlN occurs with an activation energy similar to the value for aluminum diffusion in austenite, although an AlN volume fraction effect on the embrittlement kinetics is manifested as decreases in activation energy with decreases in the [Al]/[N] ratio of steel. In contrast, the presence of niobium substantially retards the kinetics of particle embrittlement in steel containing 120–200 ppm N. Observations of AlN precipitates coated with Nb(C,N) indicate that the decreases in embrittlement kinetics are related to a reduction in the potential for AlN ripening during austenitization.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Formation and dissociation of shear-induced high-energy dislocations: insight from molecular dynamics simulations

Solid-phase processing (SPP) allows one to create complex microstructures, not achievable via thermal processing alone. The resulting structures exhibit a rich palette of defects, both thermal and non-thermal, including defect substructures, such as dislocation networks. It is essential to understand the mechanisms of deformation and defect structure formation to guide SPP towards achieving desired microstructures and material properties. In this study, large-scale molecular dynamics simulations are used to investigate the effects of inhomogeneous strain distribution, that mimics deformation conditions of tribological tests, on the evolution of defects under severe shear deformation in polycrystalline Al. Analysis of defect nucleation and reaction pathways reveals that strong geometric constraints suppress the nucleation and slide of low energy dislocation 1/2$\langle110\rangle${111} but promote the nucleation and slide of high energy dislocations, such as [$1\bar1$0] (001) and 1/2 [$1\bar1$$\bar2$]($1\bar1$1). Further, a rough contact surface, characteristic to tribological tests, imposes an inhomogeneous stress field leading to inhomogeneous defect substructures due to location-dependent activation of slip systems. The results suggest that high-energy dislocations can dominate the evolution of grain structures in highly constrained environments, which should be considered in modeling plastic deformation and grain refinement during SPP.

42 ENGINEERING↗

DOE OSTI2015

Integration of microseismic monitoring data into coupled flow and geomechanical models with ensemble Kalman filter

Abstract Hydraulic stimulation of low‐permeability rocks in enhanced geothermal systems, shale resources, and CO 2 storage aquifers can trigger microseismic events, also known as microearthquakes (MEQs). The distribution of microseismic source locations in the reservoir may reveal important information about the distribution of hydraulic and geomechanical rock properties. In this paper, we present a framework for conditioning heterogeneous rock permeability and geomechanical property distributions on microseismic data. To simulate the multiphysics processes in these systems, we combine a fully coupled flow and geomechanical model with the Mohr‐Coulomb type rock failure criterion. The resulting multiphysics simulation constitutes the forecast model that relates microseismic source locations to reservoir rock properties. We adopt this forward model in an ensemble Kalman filter (EnKF) data assimilation framework to jointly estimate reservoir permeability and geomechanical property distributions from injection‐induced microseismic response measurements. We show that integration of a large number of spatially correlated microseismic data with practical ensemble sizes can lead to severe underestimation of ensemble spread, and eventually ensemble collapse. To mitigate the variance underestimation issue, two low‐rank data representation schemes are presented and discussed. In the first approach, microseismic data are projected onto a low‐dimensional subspace defined by the left singular vectors of the perturbed observation matrix. The second method uses a coarser grid for representing the microseismic data. A series of numerical experiments is presented to evaluate the performance of the proposed methods and to illustrate their applicability for assimilating microseismic data into coupled flow and geomechanical forward models to estimate multiphysics rock properties.

Tarrahi, Mohammadali↗

DOE OSTI2016

Molecular dynamics simulation of the structural, elastic, and thermal properties of pyrochlores

We present a comprehensive simulation study of the effect of composition on the structural, elastic and thermal properties of 25 different compounds from the pyrochlore family. We joined a repulsive potential to an existing interatomic potential to enable molecular dynamics simulations of conditions away from equilibrium. We systematically varied the chemistry of the pyrochlore by substituting different cations in the A and B sites of the A2B2O7 formula unit. The A cations varied from Lu3+ to La3+, and the B cations from Ti4+ to Ce4+. The lattice parameter increased steadily with increasing the radius of A or B cations, but the bulk modulus showed a decreasing trend with increasing cation radius. However, the specific heat capacity and thermal expansion coefficient remained almost unchanged with increasing the radii of A and B cations. It is of interest to note that Ce on the B site significantly reduces the specific heat capacity and thermal expansion coefficient, which could have implications for annealing of radiation damage in cerate pyrochlores. The present results are consistent with the experimental measurements, which validates these potentials for simulation of dynamical processes, such as radiation damage, in pyrochlores.

computer simulation↗

DOE OSTI2024

Jet quenching: From theory to simulation

With the explosion of data on jet-based observables in relativistic heavy-ion collisions at the Large Hadron Collider and the Relativistic Heavy-Ion Collider, perturbative Quantum Chromodynamics (pQCD)-based simulations of these processes, often interacting with an expanding viscous fluid dynamical background, have taken center stage. This review is meant to bridge the gap between theory, simulation and phenomenology of jet modification in a dense medium. We will demonstrate how the existence of such end-to-end event generators with semi-realistic or even fully realistic final states allows for the most rigorous comparisons between pQCD-based jet modification theory and experiment. State-of-the-art calculations of several jet-based observables are presented. Extensions of this theory to jets in the small systems of p–A and e–A collisions are discussed.

Physics↗

DOE OSTI2019

Developing and Testing a New Gas Distribution System for the Fluidized Bed of the Integrated Waste Treatment Unit (IWTU) Steam Reformer - 19219

The Integrated Waste Treatment Unit (IWTU) has been tasked with using steam reforming technology to treat 3.4 million liters of radioactive sodium-bearing waste (SBW) at the INL. Progress in developing steam reforming for SBW treatment has been hampered by numerous processing issues associated with the existing process equipment during simulant testing. The primary issue was high differential temperatures that developed in the fluidized bed of the denitration mineralization reformer (DMR). Based on technical evaluations and computer modeling, it was determined that the primary cause of poor fluidization in the DMR was the existing fluidized gas distributor. This distributor consisted of a combination of a ring header and four removable fluidizing gas rails. The combined ring header/rail system was unable to effectively fluidize the bed at the bottom of the DMR (around and below the ring header). As a result, a decision was made to design and test a new fluidized gas distributor. Conceptual designs focused on three different distributors: a modified design of the current system; a flat plate with bubble caps; and a conical bottom with dual plena. After testing the conceptual prototypes at Hazen Research Inc., the conical bottom distributor was chosen. A 30-day test with simulant SBW feed was then performed with the new distributor. During testing, a temperature differential of less than 2 deg. C was maintained in the DMR fluidized bed. The results supported the conclusion that poor fluidized gas distribution was the primary cause of high temperature differentials in the DMR fluidized bed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2019

Simulations of post-recrystallization grain growth in monolithic U–10Mo fuel processing

In this study, the Potts Model was used to simulate the post-recrystallization grain growth behavior of U-10 wt.% Mo (U-10Mo). The fabrication process for monolithic U-10Mo fuels involves multiple stages of rolling, during which spatially heterogeneous distributions of second phase particles form in the alloy microstructure. These particles can strongly affect the post-recrystallization grain growth process. Absent particle dissolution effects, U-10Mo grain growth has been observed stagnating at increasing grain size with increasing annealing temperature. In this work, we simulate this effect in the Potts Model by making adjustments to the non-physical Monte Carlo temperature. This modification was used to account for the weakening pinning effect of the second phases with elevated annealing temperature. Simulated average grain sizes and grain size distributions were compared with experimentally obtained results for 700 °C, 800 °C, and 900 °C anneals of duration up to 240 hours. These simulations achieved reasonable agreement with experimentally observed grain growth behavior, but showed narrower grain size distributions. It is believed that this disagreement originates from the formation of bands of second-phase particles during rolling. Finally, these results demonstrate the capability to evaluate the post-recrystallization grain growth behavior for U-10Mo as a function of previous homogenization, rolling, and annealing treatments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2017

Boundary Conditions for Diffusion-Mediated Processes within Linear Nanopores: Exact Treatment of Coupling to an Equilibrated External Fluid

In this paper, we consider a variety of diffusion-mediated processes occurring within linear nanopores, but which involve coupling to an equilibrated external fluid through adsorption and desorption. By determining adsorption and desorption rates through a set of tailored simulations, and by exploiting a spatial Markov property of the models, we develop a formulation for performing efficient pore-only simulations of these processes. Coupling to the external fluid is described exactly through appropriate nontrivial boundary conditions at the pore openings. This formalism is applied to analyze the following: (i) tracer counter permeation (TCP) where different labeled particles adsorb into opposite ends of the pore and establish a nonequilibrium steady state; (ii) tracer exchange (TE) with exchange of differently labeled particles within and outside the pore; (iii) catalytic conversion reactions where a reactant in the external fluid adsorbs into the pore and converts to a product which may desorb. The TCP analysis also generates a position-dependent generalized tracer diffusion coefficient, the form of which controls behavior in the TE and catalytic conversion processes. We focus on the regime of single-file diffusion within the pore which produces the strongest correlations and largest deviations from mean-field type behavior. Finally, behavior is quantified precisely via kinetic Monte Carlo simulations but is also captured with appropriate analytic treatments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Automated Parameterized CFD Simulations of Phase-Change Material Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers is typically computationally intensive due to the relative complex time-dependent underlying physics. This paper presents a study investigating a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 37oC nominal melting temperature, leveraging automated parametric CFD simulations. The computational domain is reduced to 2D axisymmetric for a segment of the HX comprising a region in between two fins with constant wall temperature as boundary condition. This approach allows for faster evaluations of the transient thermal characteristics for a given design space of interest. Correlations and reduced order models (ROM’s) are developed. Results indicate that the ROM’s reproduce the CFD simulations with high accuracy for designs from the training data set, however, designs from cross-validation and test sets do not always have good agreement. The limited number of designs in the training set is a possible cause for the lack of robustness of the ROM’s and thus additional training designs will be added in future work. The ROM’s solve the simulations on average 500,000 times faster than CFD.

25 ENERGY STORAGE↗

DOE OSTI2025

An open-access simulated earthquake ground-motion database for an M7 Hayward Fault earthquake in the San Francisco Bay Region

Comprehensive understanding of earthquake ground motions, particularly in the near-fault region of large-magnitude events, is limited by gaps in strong-motion data. This challenge is prominent in areas with high seismic hazard but infrequent large earthquakes where data is sparse and difficult to interpret. These data limitations lead to uncertainties in the development of site-specific ground motions, which are crucial for engineering risk assessments. To address these challenges, physics-based regional-scale ground-motion simulations have been developed. With the emergence of exaflop-scale computing ecosystems, it is now possible to simulate regional earthquake processes at unprecedented fidelity and generate the large number of fault rupture realizations necessary to characterize both intra- and inter-event ground-motion variability. This article introduces a new database of simulated earthquake ground motions, created for applications in earthquake engineering, earthquake planning, and emergency response. The inaugural version of the database features simulated ground motions for a magnitude 7 Hayward Fault earthquake in the San Francisco Bay Region (SFBR), using the EarthQuake SIMulation (EQSIM) simulation framework and the Graves–Pitarka kinematic rupture model. The aim is to provide high-fidelity, spatially dense, three-component motions generated on the Department of Energy’s (DOE) newest generation of graphics processing unit (GPU)-accelerated supercomputers. These motions are being made openly available to the engineering, scientific, and disaster planning communities. In addition, this work develops protocols for the efficient dissemination of these large data sets and emphasizes community engagement to build confidence in their application. This article discusses the methodology behind the data, underlying software verification and validation, scalable data management, and a user interface for data access. The goal is to facilitate widespread use and elicit expert feedback to maximize the utility and exploitation of simulated motions. While the initial focus is on the San Francisco Region, simulations for additional regions will be added as the DOE program progresses.

Simulated ground-motion database↗

DOE OSTI2018

Error Analysis of the Thornado Supernova Simulation Code

The iron core collapse of a supernova throws 1.5 solar masses at speeds up to 35% the speed of light into a space with radius 50km, which creates a high density environment capable of generating neutrinos, producing gravitational waves, and driving chemical evolution of matter past iron. By designing a simulation of the process, we can better understand the origins and properties of these rare phenomena which are not seen anywhere else in the universe. Typically, supernovas are simulated with a finite volume method, but we want to incorporate a Nodal Discontinuous Galerkin Method to better utilize high performance computing because the DG method allows us to increase precision without increasing the number of computational elements. This project focuses on the fluid dynamics during collapse and the resulting shockwaves. To understand this, we improved on an existing simulation by comparing benchmark tests against published and theoretical results to identify and correct sources of error. The DG Method creates discontinuous functions across computational elements by solving the Euler fluid equations computed over Runge-Kutta time steps. The DG method allows for discontinuous functions to be built from polynomials defined in individual elements; however, it can lead to unphysical oscillation near discontinuities. We apply the characteristic slope limiter to reduce these affects, but when applied globally it leads to a 1% increase in absolute error. By introducing a Trouble Cell Indicator, which marks certain fluid elements based on the rate of change of observables, we can constrain the limiter to elements near shocks and manage the resulting absolute error. By constraining the limiter to discontinuous regions, we can remove what we know will be unrealistic without impacting naturally occurring oscillations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Thermal Neutron Scattering Research at ORNL – ENDF File Validation [Slides]

ENDF File Validation is a process developed to simulate VISION experiments. This method showed significant differences in inelastic spectra compared to experimental measurements. Differential data should be incorporated into thermal neutron scattering evaluation process to ensure accurate TSL files. The Polyethylene Discrepancy Issue has been found and corrected. The resulting ENDF manual change has been submitted and is currently under review.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Differences in total stem value when merchandizing with a tracked processor versus a knuckle–boom loader in Pinus taeda

Using tracked processors over knuckle–boom loaders to increase total value per tree when merchandizing timber on the landing has become a topic of interest in the southeastern region of the United States. This study compared merchantability values, product classes, and product weights of loblolly pine (Pinus taeda) for both machines to determine if there was a significant difference between machines when processing the same tree. In order to process the same tree twice, the chains from the tracked processor had to be removed from the bottom saw bar. This allowed the processor to simulate the merchandising process without actually marking or cutting the tree for a more realistic comparison. Data were analyzed using paired t–tests and two–way ANOVA models. Results depicted that when diameter and total lengths are visually estimated by the knuckle–boom loader. A significant difference in value occurs, however, once diameter and total length are modified to match the tracked processors for more accurate measurements, no difference in value was seen. Here, these results demonstrate that until mill specifications become more stringent, there is little incentive for loggers to purchase a tracked processor if their only motivation is to increase merchantability values.

42 ENGINEERING↗

DOE OSTI2021

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗