Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulator Study”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hanford Composite Analysis Special Studies FY2025 - Groundwater Flow Simulation for use in a Comparison of Simulated Concentrations in the Saturated Zone Estimated by the Plateau-to-River Model Versions 8.3 and 9.1

In order to meet the requirements of U.S. Department of Energy (DOE) in DOE O 435.1, Radioactive Waste Management, a composite analysis (CA) must be completed for the Central Plateau at the Hanford Site. The CA requires estimates of fate and transport of radionuclides in the groundwater from multiple sources within the modeling domain. This Environmental Calculation File (ECF) details the application of the Plateau-to-River (P2R) Model (CP-57037, Rev. 3, Model Package Report for the P2R Model: Version 9.1) to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support a special analysis associated with the Hanford Site CA. The simulated flow field will support the simulation of fate and transport of contaminants for comparison to the results obtained as part of the CA to evaluate impact of the updated version of the P2R Model.

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

Particle Acceleration in Kinetic Simulations of Nonrelativistic Magnetic Reconnection with Different Ion–Electron Mass Ratios

By means of fully kinetic particle-in-cell simulations, we study whether the proton-to-electron mass ratio m i /m e influences the energy spectrum and underlying acceleration mechanism during magnetic reconnection. While kinetic simulations are critical to studying particle acceleration during magnetic reconnection, a reduced m i /m e is often used to alleviate the demanding computing resources, which leads to artificial scale separation between electron and proton scales. Recent kinetic simulations with high mass ratios have indicated new regimes of reconnection, as electron pressure anisotropy develops in the exhaust region and supports extended current layers. Here, we study whether different m i /m e changes the particle acceleration processes by performing a series of simulations with different mass ratio (m i /m e = 25–400) and guide field strength in a low-β plasma. We find that mass ratio does not strongly influence reconnection rate, magnetic energy conversion, ion internal energy gain, plasma energization processes, ion energy spectra, and the acceleration mechanisms for high-energy ions. Simulations with different mass ratios are different in electron acceleration processes, including electron internal energy gain, electron energy spectrum, and the acceleration efficiencies for high-energy electrons. We find that high-energy electron acceleration becomes less efficient when the mass ratio gets larger because the Fermi -like mechanism associated with particle curvature drift becomes less efficient. These results indicate that when particle curvature drift dominates high-energy particle acceleration, the further the particle kinetic scales are from the magnetic field curvature scales (~d i ), the weaker the acceleration will be.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Simulating strongly correlated molecules with a superconducting quantum processor

Many of the biggest challenges in expanding the nation’s access to clean and low-cost energy resources are fundamentally chemistry or materials challenges. An important case is the development of new catalysts for the up-conversion of cheap and readily available materials such as methane or water into materials suitable for use as a fuel such as methanol or oxygen. To understand and exploit such processes, computer simulations of chemical reactions provide a natural complement to experimental studies. Unfortunately, most catalytic reactions involve so-called “strongly correlated” molecules which are notoriously difficult to study with simulation algorithms that can be executed on existing (classical) computers. The recent growth in quantum information science offers an alternative potential route for simulating these difficult systems. As a result, an increasing number of computational chemists are becoming interested in quantum computing. At the same time, quantum information scientists have identified chemistry simulation as a possible first demonstration of a quantum computer providing an improvement over a classical computer. The objective of this project is to accurately simulate strongly correlated molecules on a quantum processor. To meet the high challenges of this objective, new hybrid quantum/classical algorithms will be co-designed with advanced quantum gate developments and computed on customized quantum hardware. Some of the developed techniques will be transferable to study other molecular systems, while the project as a whole will help define better strategies for advancing the quantum simulation of matter more generally.

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

The {10$\overline{1}$2} non-cozone twin-twin interactions in Mg: A stability and mobility study using 3-D atomistic simulations

Due to easy activation of tensile twinning in Mg, multiple {10$\overline{1}$2} twin variants can interact with each other and form twin-twin boundaries over the course of plastic deformation. Previous studies using 2-D settings provide only a partial understanding of these interactions, especially the non-cozone ones. Here, atomistic simulations are used to study the 3-D structural characteristic and evolution of the non-cozone {10$\overline{1}$2} twin-twin junctions. The study reveals the existence of new twin-twin boundaries (TTBs) such as TTB BP and TTB K2 , formed after the interaction between the basal prismatic and conjugate twin interfaces with the coherent twin boundary. Here, for both non-cozone twin-twin interactions, the {10$\overline{1}$2} TTB and its associated twin-twin junctions are found to play a major role in the {10$\overline{1}$2} twin's stability and mobility. Specifically, they promote the growth of the 3-D twin in both the normal and forward directions during the interaction and hinder the detwinning process upon unloading.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Probing the exact form and doping preference of magnesium in ordinary Portland cement clinker phases: A study from experiments and DFT simulations

Highlights: • The doping behaviors of Mg in OPC clinker phases were systematically studied. • About 1.3 wt% MgO enter OPC clinker phases by forming substitutional solid solutions. • In the 4 composed minerals, C{sub 3}S, C{sub 3}A and C{sub 4}AF can host more Mg ions than C{sub 2}S. • Mg ions prefer to replace Ca in C{sub 3}S, C{sub 2}S and C{sub 3}A while replace Fe and Ca in C{sub 4}AF. • Both the large structure distortions and formation energies governs the low solubility. A systematic study was performed to learn the exact existence form and doping behaviors of magnesium in ordinary Portland cement (OPC) clinker. Our results show that about 1.3 wt% MgO incorporates into OPC clinker phases by forming substitutional solid solutions while excess Mg accumulate and exist as periclase. Results from Rietveld refinement present a substitution preference in C{sub 3}S and C{sub 3}A over C{sub 2}S. More in-depth analyses from density function theoretical simulations show that in C{sub 3}S, C{sub 2}S and C{sub 3}A, the substitution occurs by replacing Ca while in C{sub 4}AF, it probably occurs by substituting Fe and Ca ions. The large structure distortions and sharp increase in formation energies with increasing MgO determines its low solubility. This work provides a clear understanding of the existence states and the intrinsic mechanism governing doping behaviors of Mg, thus should be very important in guiding the synthesis of OPC clinker by utilizing high-Mg limestone.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulation and generalized Langevin equation study of lipid subdiffusion in biomembrane phases

In biomembranes lipid mobility exhibits deviations from the classical diffusive behavior of Brownian particles, i.e. ‘anomalous’ diffusion. The question arises as to how this anomalous diffusive behavior varies in gel, ripple, and fluid biomembrane phases. Here, in this study, we perform all-atom molecular dynamics (MD) simulations of dimyristoylphosphatidylcholine (DMPC) bilayers in the three different phases and analyze the results using the framework of the Generalized Langevin Equation (GLE). This analysis emphasizes sub-diffusive behavior on the relatively short, picosecond-nanosecond timescales, capturing local molecular constraints and transient caging effects during the crossover of atomic dynamics from vibrational to incipient anharmonic motion. The ripple and gel phases are found to exhibit strong transient caging and prolonged memory effects resulting in distinct sub-diffusive behavior. The role of hydrogen bonding in lipid confinement is also examined, demonstrating its influence on phase-dependent molecular ordering and on short-time diffusional constraints. These findings demonstrate the GLE framework’s utility in characterizing molecular transport and lipid dynamics, with implications for longer timescale membrane dynamics.

Malik, Sheeba [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2024

Consequences of plane-strain and plane-stress assumptions in fully coupled chemo-mechanical Li-ion battery models

In Li-ion battery research, it is common to simulate chemo-mechanical phenomena in reduced dimensions (e.g., 2-D) as opposed to fully resolve these complex physics in 3-D. It is common to assume either (1) the out-of-plane strain is negligible (commonly referred to as plane-strain), or (2) the out-of-plane stress is negligible (commonly referred to as plane-stress). However, there is typically little consideration as to the quantitative consequences of these approximations. Furthermore, the influence of these out-of-plane assumptions can be compounded and convoluted when chemo-mechanics models implement so-called “fully coupled” formulations, where the local species concentrations influence the stress-state and the stress-state influences the local species fluxes. Here, the present manuscript explores the implications of using plane-stress and plane-strain assumptions in 2-D as compared to simulating a full 3-D electrode particle. This comparative study includes simulating both isotropic and anisotropic particle intercalation where the particles can be surrounded by either a liquid or solid electrolyte. Additionally, common Li-ion battery-model metrics such as the state-of-stress, intercalation fraction distribution, and specific capacity are compared, while also considering the effects of particle size and C-rate. As alternatives to the pure plane-strain and plane stress approximations, two modified plane-strain assumptions are found to better approximate the fully coupled chemo-mechanical 3-D behavior.

25 ENERGY STORAGE↗
DOE OSTI2022

Connecting galaxy evolution in clusters with their radial profiles and phase space distribution: results from the IllustrisTNG hydrodynamical simulations

ABSTRACT We study the population of galaxies around galaxy clusters in the hydrodynamic simulation suite IllustrisTNG 300-1 to study the signatures of their evolutionary history on observable properties. We measure the radial number density profile, phase space distribution, and splashback radius for galaxies of different masses and colours over the redshift range z = 0−1. The three primary physical effects which shape the galaxy distribution within clusters are the galaxy quenching, angular momentum distribution, and dynamical friction. We find three distinct populations of galaxies by applying a Gaussian mixture model to their distribution in colour and mass. They have distinct evolutionary histories and leave distinct signatures on their distribution around cluster haloes. We find that low-mass red galaxies show the most concentrated distribution in clusters and the largest splashback radius, while high-mass red galaxies show a less concentrated distribution and a smaller splashback radius. Blue galaxies, which mostly quench into the low-mass red population, have the shallowest distribution within the clusters, with those on radial orbits quenched rapidly before reaching pericentre. Comparison with the distribution of galaxies from the Dark Energy Survey survey around Sunyaev–Zeldovich clusters from the Atacama Cosmology Telescope and South Pole Telescope surveys shows evidence for differences in galaxy evolution between simulations and data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Interface capturing simulations of bubble population effects in PWR subchannels

As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

EXTRACTING HUMAN RELIABILITY FINDINGS FROM HUMAN FACTORS STUDIES IN THE HUMAN SYSTEMS SIMULATION LABORATORY

Modernization of U.S. nuclear power plants (NPPs) is widespread, with most plants currently replacing and transitioning equipment, control systems, and human system interfaces (HSI)s from analog to digital displays. This conversion remedies the obsolescence of analog parts along with needs for increased intuitiveness of design, safety, and capabilities. The Human Factors and Reliability team at Idaho National Laboratory (INL) carried out twelve control room modernization studies in the newly designed Human Systems Simulation Laboratory (HSSL) over nine years. The HSSL was constructed as a testbed for evaluating human factors techniques and performance measures, HSI frameworks, and cutting-edge operational concepts in NPPs. Installing a full-scope training simulator enabled direct design and evaluation work on the same instrumentation and control (I&C) and HSIs located at U.S. plants. The subsequent addition of glass top bays afforded crews opportunities to implement operations via the simulator using full-scale representations of their home NPP. Additionally, functional HSI prototypes were created, providing an environment for operator-in-the-loop benchmark studies. The HSSL has assisted in upgrades of six commercial NPP control rooms and served as an invaluable proving ground for new NPP operations technology. Human reliability analysis (HRA) was not originally the focus of the studies; however, data relating to HRA such as type and frequency of human errors can be extracted from the studies. INL is currently extracting data from the HSSL study reports to apprise how information gathered from simulation, HSI, and other related studies can create a broad look across different data sources to help inform HRA methods.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Comparison of Error Rate Depending on Operator Expertise and Simulator Complexity

This paper analyzes operator's error rate from experiments, depending on the expertise and simulator complexity. This study uses the Rancor Microworld, a simplified simulator developed by INL, and the Compact Nuclear Simulator (CNS), a less simplified simulator developed by the Korea Atomic Energy Research Institute (KAERI). The error rates were measured from simulation data of a total of 72 participants, and the collected error rate data were analyzed using analysis of variance (ANOVA) test.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Can Isoprene Oxidation Explain High Concentrations of Atmospheric Formic and Acetic Acid over Forests?

Formic and acetic acid concentrations are particularly high over forested areas of the world. However, the gas-phase mechanisms for producing these acids are poorly understood even for isoprene, the globally dominant biogenic hydrocarbon. We quantified formic and acetic acid production from reactions of hydroxyl radical (OH) (between high and low ranges of nitric oxide (NO) levels) with isoprene, methacrolein (MACR), isoprene epoxydiol (IEPOX), isoprene hydroxy hydroperoxide (ISOPOOH), and α-pinene from the focused isoprene experiments at California Institute of Technology (FIXCIT) laboratory chamber study. We find that (i) OH oxidation of MACR, IEPOX, and ISOPOOH are sources of formic acid, (ii) isoprene peroxy radical isomerization and associated photolysis oxidation products are potentially important sources of organic acids, and (iii) high levels of NO generally suppress organic acid formation from OH oxidation of isoprene. We modified existing chemical mechanisms for isoprene oxidation to account for organic acid production pathways observed in the FIXCIT study. We simulated organic acid production during the Southeastern Oxidant and Aerosol Study using the updated chemical mechanisms and represented acetic acid within a factor of 2 but still underpredicted formic acid by a factor of 6. While we cannot explain ambient formic acid with explicit chemical mechanisms, the FIXCIT results suggest that the oxidation of isoprene could account for as much as 70% of the global annual production of formic acid from gas-phase reactions.

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

Optical modeling of a radioisotope atmospheric chamber

Abstract Volatile radioisotopes represent a substantial health risk when released into the environment. To better understand the environmental fate of radioisotopes, the authors constructed a cylindrical steel-walled chamber to simulate the atmospheric processing of volatile radioactive gases. Optical modeling was performed for the properties of simulated sunlight in the chamber to better characterize atmospheric reaction studies. Optical simulations were performed using two wall materials (steel and thin-film silica) and validated against experimental measurements. This optical analysis methodology can be used to improve the fidelity of atmospheric models by accounting for optical inhomogeneities enabling a firmer grasp of radioisotopes’ environmental fate. Graphical abstract

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Hands-on with OWL: the Oak–Ridge Wang–Landau Monte Carlo software suite

The Oak–Ridge Wang–Landau (OWL) package is an open-source scientific software specialized for large-scale, Monte Carlo simulations for the study of materials properties at finite temperature. In this paper, we discuss the main features and capabilities of OWL, followed by detailed descriptions of building and running the code. The readers will be guided through the usage and functionality of the code with a few hands-on examples. This paper is based on a tutorial on OWL given at the 32nd Center for Simulational Physics Workshop on Recent Developments in Computer Simulation Studies in Condensed Matter Physics.

Li, Ying Wai↗
DOE OSTI2021

Complementary Study of Radiative Heat Transfer and Flow Physics from Moderate-scale Hydrocarbon Pool Fire Simulations

As part of the Advanced Simulation and Computing Verification and Validation (ASCVV) program, a 0.3-m diameter hydrocarbon pool fire with multiple fuels was modeled and simulated. In the study described in this report, systematic examination was performed on the radiation model used in a series of coupled Fuego/Nalu simulations. A calibration study was done with a medium-scale methanol pool fire and the effect of calibration traced throughout the radiation model. This analysis provided a more detailed understanding of the effect of radiation model parameters on each other and on other quantities in the simulations. Heptane simulation results were also examined using this approach and possible areas for further improvement of the models were identified. The effect of soot on radiative losses was examined by comparing heptane and methanol results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗