Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “criticality simulations”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A comprehensive model to capture electrical discharge and spark channel evolution during spark-ignition processes

Emissions reduction through engine efficiency improvements is a priority for automakers who have turned to unconventional engine operation such as highly dilute, boosted, and stratified charge. Given its importance to flame initiation and sustained turbulent flame propagation, reliable and accurate spark ignition models are necessary to design ignition systems that reduce cyclic indicated mean effective pressure (IMEP) variability and increase engine efficiency in these operation modes. In this paper, secondary electric circuit, short-circuit, blowout, and re-strike sub-models are added to the Lagrangian-Eulerian spark ignition (LESI) model to simulate electrical discharge and spark channel elongation in an inert cross-flow combustion vessel. First, the physics of the underlying sub-models are described and the governing equations discussed with the spark channel voltage expression playing a critical role. Then, the experimental and simulations setups are presented. The results section begins with the derivation of spark channel voltage from experimental results. Then, electrical discharge and spark channel elongation simulations in inert flows are carried out using LESI and compared against experimental results. Finally, the results validate the model's ability to accurately predict spark channel elongation, as well as the occurrence of short-circuits, blowouts, re-strikes, and end of discharge.

42 ENGINEERING↗
DOE OSTI2024

Critical role of slags in pitting corrosion of additively manufactured stainless steel in simulated seawater

Abstract Pitting corrosion in seawater is one of the most difficult forms of corrosion to identify and control. A workhorse material for marine applications, 316L stainless steel (316L SS) is known to balance resistance to pitting with good mechanical properties. The advent of additive manufacturing (AM), particularly laser powder bed fusion (LPBF), has prompted numerous microstructural and mechanical investigations of LPBF 316L SS; however, the origins of pitting corrosion on as-built surfaces is unknown, despite their utmost importance for certification of LPBF 316L SS prior to fielding. Here, we show that Mn-rich silicate slags are responsible for pitting of the as-built LPBF material in sodium chloride due to their introduction of deleterious defects such as cracks or surface oxide heterogeneities. In addition, we explain how slags are formed in the liquid metal and deposited at the as-built surfaces using high-fidelity melt pool simulations. Our work uncovers how LPBF changes surface oxides due to rapid solidification and high-temperature oxidation, leading to fundamentally different pitting corrosion mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tundra SGS – Critical Importance of Site-Specific Injectivity Testing for Carbon Dioxide Injection Well Simulation Calibration

Paper presented at 16th International Conference on Greenhouse Gas Control Technologies (GHGT-16), October 23–27, 2022, Lyon, France. Brine injectivity tests were conducted on both the Broom Creek (BC) and Deadwood (DW) storage formations in one of the cored stratigraphic test wells prior to finalized simulations. The injectivity tests entailed injection of brine into the formations, followed by shutting in the wells. The shut-in duration represents the fall-off period which was used to carry out pressure transient analyses of the tests. The analysis of injectivity tests indicated that core-indicated and modeled permeability distributions should be conservatively multiplied by 2.5× and 5× in the BC and DW geologic models, respectively, prior to numerical simulations of CO 2 injection using Computer Modelling Group’s (CMG’s) GEM equation-of-state reservoir simulator. The scaler permeability multipliers allowed the fluid flow simulations to fully represent well test-measured injection well rates in the respective formations.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Hydration-driven stabilization and volume collapse of grain boundaries in Mg2SiO4 forsterite predicted by first-principles simulations

Grain boundaries in mantle minerals are of critical importance to geophysical and geochemical processes of the Earth’s interior. One of the fundamental issues is to understand how the water (H 2 O) component influences the properties of grain boundaries in silicate materials. Here, we report the results of the structure and stability of several tilt grain boundaries in Mg 2 SiO 4 forsterite over the pressure range 0 to 15 GPa using density functional theory-based first-principles simulations. Here, the results suggest greater energetic stability and hydration-driven volume collapse (negative excess volume) at zero pressure for the majority of hydrous grain boundaries relative to the anhydrous (dry) ones. All the hydrous grain boundaries become increasingly favorable at elevated pressures as the calculated hydration enthalpy systematically decreases with increasing pressure. The hydrous components at the interfacial regions are predominantly in the hydroxyl form and, to a lesser extent, in the molecular H 2 O form. Their calculated ratio ranges from 1.6 to 8.7 among the different grain boundary configurations. Our structural analysis also reveals that the hydroxyls are bound to either both Mg and Si or to Mg only. In comparison, the molecular species are bound only to Mg sites. Besides direct oxygen-hydrogen bonding, intermolecular hydrogen bonding becomes important with compression. On the basis of our results, we suggest that local atomic rearrangements caused by dissociative adsorption of water facilitate efficient compaction of the boundary interfaces, which, in turn, results in greater relative stability of hydrous grain boundaries. This means that water prefers to be incorporated within the grain boundaries over the bulk of silicate materials.

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

Accelerating transients with NekRS: GPU overlapping domain implementation and multi-rate timestepping

The simulation of nuclear transients using Computational Fluid Dynamics (CFD) presents significant computational challenges due to the inherent complexity and the wide separation in temporal scales between various flow physical phenomena. These disparities lead to high computational costs, often making the simulation of transients impractical without advanced techniques. Consequently, multiple research initiatives are being pursued by the NEAMS thermal-hydraulic area, some driven by academic institutions and some by national laboratories. Overall, they are exploring novel methods to make transient simulations more feasible and efficient. This report delves into recent advancements within the CFD code NekRS, specifically those achieved in Fiscal Year 2024 under the CONNECT effort, aimed at improving the performance and feasibility of transient simulations. The first major advancement involves the porting of NekRS to Aurora, one of the Department of Energy’s (DOE) most powerful supercomputers. Additionally, the report discusses the implementation of an overlapping domain capability within NekRS. This novel GPU-accelerated capability allows different spatial regions of the domain to be solved independently, enhancing the code’s efficiency, particularly when running large-scale simulations in complex domains. The scalability of this approach is demonstrated, highlighting its potential to transform how transients are approached in CFD simulations. Lastly, the report focuses on how this overlapping domain capability specifically accelerates transient simulations through multi-rate timestepping. By decoupling different regions and facilitating faster computations, this method offers a promising pathway to making nuclear transient simulations more computationally feasible, addressing one of the critical bottlenecks in the field. Together, these advancements represent a significant leap forward in transient simulation technology, bringing closer the possibility of handling highly complex nuclear scenarios with greater efficiency and accuracy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A Visual Analytic Platform for Interactive Validation of Human Mobility Simulations

Human mobility insights guide domain experts in an array of decisions, including critical infrastructure design, disaster response, epidemic modeling, national security, and policy making. Due to the inherent noise and privacy concerns in real-world individual-level mobility data, it is often preferred to leverage simulators that generate synthetic mobility data instead. However, it is critical to inspect and validate the output of such simulators to ensure the synthetic data is aligned with the characteristics of the population and the area of interest known to domain experts. While there exist many quantitative approaches for validating synthetic data, we argue it is also important to also validate such data qualitatively to capture aspects that are known to domain experts but difficult to quantify. In this work, we demonstrate a visual analytic platform that empowers domain experts to interact with their simulation outputs along spatial and temporal dimensions. By augmenting automated techniques and human skills, our visual analytic platform is a step towards interactive capabilities for model steering and quality control of mobility simulators.

Monadjemi, Shayan↗
DOE OSTI2025

Genetic algorithm optimization of nuclear criticality experiment for reduction of intermediate-energy 239 Pu nuclear data uncertainties

Nuclear criticality experiments are conducted to investigate specific nuclear data important for safe handling and storage of fissile materials, reactor design and operation, and the validation of radiation transport codes. Incorrect or uncertain nuclear data can prohibitively impact operational safety limits, reactor licensing, and predictive simulation capability; therefore, integral measurements from criticality experiments are necessary and should be performed frequently. To maximize the impact of the integral measurements, it is important to consider experiment geometry, material selection, and component dimensions. When taking these considerations into account, the experiment design process becomes iterative and very time intensive. This work utilizes a genetic algorithm to efficiently explore potential nuclear criticality experiment designs for the Laboratory Directed Research & Development project PARADIGM (PARallel Approach of Differential and InteGral Measurements) at Los Alamos National Laboratory. In this paper, the building blocks of the genetic algorithm are discussed in detail, the genetic algorithm methodology is verified, and the genetic algorithm is used to produce three candidate experiment models for the final PARADIGM design. The three candidate models produced by the genetic algorithm consist of copper-reflected assemblies containing 14 repeating units of alumina, graphite, boron, and plutonium plates. Furthermore, in addition to the optimization results, final design considerations are also discussed for designs with a height and/or weight very close to or slightly above assembly machine operational limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A Participation Factor-Based Approach for Defining the EMT Model Boundary for Power System Simulations with Inverter-Based Resources

The increasing penetration of inverter-based resources (IBRs) introduces new challenges to power system simulations, particularly with the emergence of fast electromagnetic transient (EMT) dynamics and sub-synchronous oscillations (SSO) that require time-consuming EMT simulations. To reduce the time cost for simulating a large-scale power grid with IBRs, this paper proposes a novel participation factor-based approach for defining a critical zone for detailed EMT modeling and simulations, which includes the IBRs, synchronous generators, and the network components participating significantly in simulated contingencies. Both model-based and response-based methods are introduced for the estimation of participation factors (PFs). The case study on the 240-bus Western Electricity Coordinating Council (WECC) system demonstrates that the EMT zone determined by the proposed approach can effectively capture power system dynamics involving IBRs.

EMT simulation↗
DOE OSTI2025

Lost Muons Analysis

The Muon g-2 experiment at FermiLab represents a significant scientific endeavor aimed at measuring the anomalous magnetic moment of the muon, $a_µ$, with a precision of 140 parts per billion (ppb). During my internship at FermiLab, I focused on the simulation and analysis of Lost Muons, a critical aspect to ensure the accuracy of the experimental results. I conducted a detailed study of the beam dynamics using the existing simulation, I generated a new simulation sample enrich of events with lost muons, I analyzed it and compared the results I obtained with those from the Run-2 dataset acquired by the experiment. In this report I will briefly describe my work, discuss the results and suggest future improvements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Interpretation of wide-angle x-ray diffraction patterns from mesophase pitch-based carbon fibers – a simulation and experimental study

Carbon fiber is a critical material in a wide range of industries, where it is highly valued for its high specific strength/stiffness, excellent wear resistance, efficient electrical and thermal transport properties, chemical resistance, and low coefficient of thermal expansion. The properties of a specific carbon fiber are closely tied to its structural characteristics at all length scales. Here, in this work, we applied wide-angle x-ray diffraction to a set of heat-treated mesophase pitch-based carbon fibers, with the goal of elucidating the crystalline structures as a function of fiber orientation. To assist with analysis and interpretation of the experimental data, we employed diffraction pattern simulations using the scalar and vector forms of the Debye scattering equation to determine the influence of basal plane orientation, crystalline ordering (turbostratic-graphitic), and basal plane asymmetry on the diffraction patterns. The results presented here suggest that growth of the transverse crystallites in mesophase pitch-based carbon fiber is fixed until graphitization temperatures are reached. The work completed here provides a framework for the analysis of carbon fiber and other oriented carbon-based materials via diffraction.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗
DOE OSTI2023

Non-linear Wire Rope Isolator Model to Enhance Transportation Simulation of Fragile Equipment

Abstract Introduction Wire rope isolators are often used as means of vibration isolation for the transportation of fragile machine equipment. When deployed as suspension systems in a transportation scenario, a simple spring–damper model is often used to predict the dynamic loads that act on the structure being transported. However, this simple model has proved to be too simplistic to be used in the development of an accurate numerical model that can predict the vibration levels experienced by the inner components of fragile equipment during its transportation. This paper describes the experimental tests conducted on a wire rope isolator used for the transportation of the prototype SSR1 cryomodule. Materials and methods A hysteretic Bouc–Wen model has been used to analytically describe the force–deformation relationship of the wire rope isolators. The developed model of the isolator has been implemented in a larger model to simulate the actual transportation of the prototype SSR1 cryomodule, a section of the new PIP-II linear accelerator under construction at Fermilab. A series of multibody dynamic simulations with rigid and flexible components was used to numerically determine the acceleration of some critical components. Results and conclusions An actual experimental transportation was simulated using two numerical models: the developed Bouc–Wen model and a conventional spring–damper model. It is shown how the Bouc–Wen formulation of the isolator characteristics drastically improves the correspondence between experimental and simulated results if compared to a spring–damper model, especially in the range of 0–30 Hz which is the most critical for transportation problems.

42 ENGINEERING↗
DOE OSTI2022

Transient Storage Model Parameter Optimization Using the Simulated Annealing Method

Hyporheic exchange in streams is critical to ecosystem functions such as nutrient cycling along river corridors, especially for slowly moving or small stream systems. The transient storage model (TSM) has been widely used for modeling of hyporheic exchange. TSM calibration, for hyporheic exchange, is typically used to estimate four parameters, including the mass exchange rate coefficient, the dispersion coefficient, stream cross-sectional area, and hyporheic zone cross-sectional area. Prior studies have raised concerns regarding the non-uniqueness of the inverse problem for the TSM, that is, the occurrence of different parameter vectors resulting in TSM solution that reproduces the observed in-stream tracer break through curve (BTC) with the same error. This leads to practical non-identifiability in determining the unknown parameter vector values even when global-optimal values exist, and the parameter optimization becomes practically non-unique. To address this problem, we applied the simulated annealing method to calibrate the TSM to BTCs, because it is less susceptible to local minima-induced non-identifiability. A hypothetical (or synthetic) tracer test data set with known parameters was developed to demonstrate the capability of the simulated annealing method to find the global minimum parameter vector, and it identified the “hypothetically-true” global minimum parameter vector even with input data that were modified with up to 10% noise without increasing the number of iterations required for convergence. The simulated annealing TSM was then calibrated using two in-stream tracer tests conducted in East Fork Poplar Creek, Tennessee. In conclusion, simulated annealing was determined to be appropriate for quantifying the TSM parameter vector because of its search capability for the global minimum parameter vector.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Identifying Hydrometeorological Factors Influencing Reservoir Releases Using Machine Learning Methods

Simulation of reservoir releases plays a critical role in social-economic functioning and our nation's security. How-ever, it is challenging to predict the reservoir release accurately because of many influential factors from natural environments and engineering controls such as the reservoir inflow and storage. Moreover, climate change and hydrological intensification causing the extreme precipitation and temperature make the accurate prediction of reservoir releases even more challenging. Machine learning (ML) methods have shown some successful applications in simulating reservoir releases. However, previous studies mainly used inflow and storage data as inputs and only considered their short-term influences (e.g, previous one or two days). In this work, we use long short-term memory (LSTM) networks for reservoir release prediction based on four input variables including inflow, storage, precipitation, and temperature and consider their long-term influences. We apply the LSTM model to 30 reservoirs in Upper Colorado River Basin, United States. We analyze the prediction performance using six statistical metrics. More importantly, we investigate the influence of the input hydrometeorological factors, as well as their temporal effects on reservoir release decisions. Results indicate that inflow and storage are the most influential factors but the inclusion of precipitation and temperature can further improve the prediction of release especially in low flows. Additionally, the inflow and storage have a relatively long-term effect on the release. These findings can help optimize the water resources management in the reservoirs.

Fan, Ming↗
DOE OSTI2022

Platinum nanoparticle compression: Combining in situ TEM and atomistic modeling

The mechanical behavior of nanoparticles governs their performance and stability in many applications. However, the small sizes of technologically relevant nanoparticles, with diameters in the range of 10 nm or less, significantly complicate experimental examination. These small nanoparticles are difficult to manipulate onto commercial test platforms and deform at loads that are below the typical noise floor of the testing instruments. Here, we synthesized small platinum nanoparticles directly onto a mechanical tester and used a modified nanomanipulator to enhance load resolution to the nanonewton scale. We demonstrated the in situ compression of an 11.5-nm platinum nanoparticle with simultaneous high-resolution measurements of load and particle morphology. Molecular dynamics simulations were performed on similarly sized particles to achieve complementary measurements of load and morphology, along with atomic resolution of dislocations. The experimental and simulation results revealed comparable values for the critical resolved shear stress for failure, 1.28 and 1.15 GPa, respectively. Altogether, this investigation demonstrated the promise of, and some initial results from, the combination of atomistic simulations and in situ experiments with an unprecedented combination of high spatial resolution and high load resolution to understand the behavior of metal nanoparticles under compression.

42 ENGINEERING↗
DOE OSTI2025

Uncertainty Visualization of Critical Points of 2D Scalar Fields for Parametric and Nonparametric Probabilistic Models

This paper presents a novel end-to-end framework for closed-form computation and visualization of critical point uncertainty in 2D uncertain scalar fields. Critical points are fundamental topological descriptors used in the visualization and analysis of scalar fields. The uncertainty inherent in data (e.g., observational and experimental data, approximations in simulations, and compression), however, creates uncertainty regarding critical point positions. Uncertainty in critical point positions, therefore, cannot be ignored, given their impact on downstream data analysis tasks. Here, in this work, we study uncertainty in critical points as a function of uncertainty in data modeled with probability distributions. Although Monte Carlo (MC) sampling techniques have been used in prior studies to quantify critical point uncertainty, they are often expensive and are infrequently used in production-quality visualization software. We, therefore, propose a new end-to-end framework to address these challenges that comprises a threefold contribution. First, we derive the critical point uncertainty in closed form, which is more accurate and efficient than the conventional MC sampling methods. Specifically, we provide the closed-form and semianalytical (a mix of closed-form and MC methods) solutions for parametric (e.g., uniform, Epanechnikov) and nonparametric models (e.g., histograms) with finite support. Second, we accelerate critical point probability computations using a parallel implementation with the VTK-m library, which is platform portable. Finally, we demonstrate the integration of our implementation with the ParaView software system to demonstrate near-real-time results for real datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Evaluated Nuclear Data Requirements in Support of Fission and Fusion Applications [Slides]

This presentation touches on thermal scattering cross section and scattering lengths including criticality benchmarks (fission). Secondly, development and completeness of charged ­particle libraries: criticality benchmarks (fission), neutron source (fission/fusion), simulation benchmarks and light elements generation (fusion). Additionally, this presentation details on the covariances (methodology and completeness) with adjustment and optimization of nuclear data libraries (fission/fusion). Finally nuclear data reproducibility and conclusions are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Multireference Approach to Electron and Electron–Nuclear Dynamics in Nanomaterials (Final Report)

Many important chemical and physical phenomena involve dynamics on large number of electronic states. Thus, there is a critical need to develop methods to simulate dynamics in dense manifolds of states. Towards this end, we have: a) developed the multiple cloning in dense manifolds of states (MCDMS) method, which is capable of accurately modeling the quantum mechanical coherence between populations on a large number of electronic states, b) implemented MCDMS into the free, open-source PySpawn software package, c) developed graphics processing unit-accelerated algorithms modeling electron dynamics in light fields via Floquet time-dependent configuration interaction (F-TDCI), and d) critically compared different orbital bases in order to achieve an accurate and efficient F-TDCI expansion. This grant ended in August 2020, when our group moved to from Michigan State University to Stony Brook University, where this project continues under grant number DE-SC0021643.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗