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 559 records · Page 31

CAFQA: A Classical Simulation Bootstrap for Variational Quantum Algorithms

Classical computing plays a critical role in the advancement of quantum frontiers in the NISQ era. In this spirit, this work uses classical simulation to bootstrap Variational Quantum Algorithms (VQAs). VQAs rely upon the iterative optimization of a parameterized unitary circuit (ansatz) with respect to an objective function. Since quantum machines are noisy and expensive resources, it is imperative to classically choose the VQA ansatz initial parameters to be as close to optimal as possible to improve VQA accuracy and accelerate their convergence on today’s devices. This work tackles the problem of finding a good ansatz initialization, by proposing CAFQA, a Clifford Ansatz For Quantum Accuracy. The CAFQA ansatz is a hardware-efficient circuit built with only Clifford gates. In this ansatz, the parameters for the tunable gates are chosen by searching efficiently through the Clifford parameter space via classical simulation. The resulting initial states always equal or outperform traditional classical initialization (e.g., Hartree-Fock), and enable high-accuracy VQA estimations. CAFQA is well-suited to classical computation because: a) Clifford-only quantum circuits can be exactly simulated classically in polynomial time, and b) the discrete Clifford space is searched efficiently via Bayesian Optimization. For the Variational Quantum Eigensolver (VQE) task of molecular ground state energy estimation (up to 18 qubits), CAFQA’s Clifford Ansatz achieves a mean accuracy of nearly 99% and recovers as much as 99.99% of the molecular correlation energy that is lost in Hartree-Fock initialization. CAFQA achieves mean accuracy improvements of 6.4x and 56.8x, over the state-of-the-art, on different metrics. Here, the scalability of the approach allows for preliminary ground state energy estimation of the challenging chromium dimer (Cr2) molecule. With CAFQA’s high-accuracy initialization, the convergence of VQAs is shown to accelerate by 2.5x, even for small molecules. Furthermore, preliminary exploration of allowing a limited number of non-Clifford (T) gates in the CAFQA framework, shows that as much as 99.9% of the correlation energy can be recovered at bond lengths for which Clifford-only CAFQA accuracy is relatively limited, while remaining classically simulable.

bayesian optimization↗

Simulation and Analyses of Multi-Body Separation in Launch Vehicle Staging Environment

The development of methodologies, techniques, and tools for analysis and simulation of multi-body separation is critically needed for successful design and operation of next generation launch vehicles. As a part of this activity, ConSep simulation tool is being developed. ConSep is a generic MATLAB-based front-and-back-end to the commercially available ADAMS. solver, an industry standard package for solving multi-body dynamic problems. This paper discusses the 3-body separation capability in ConSep and its application to the separation of the Shuttle Solid Rocket Boosters (SRBs) from the External Tank (ET) and the Orbiter. The results are compared with STS-1 flight data.

Pamadi, Bandu N.↗

Residual Strength Characterization of a Curved Integrally-Stiffened Panel

Over the years, Finite-element fracture simulation methodology has been very well established at NASA Langley to predict the residual strength of damaged aircraft structures. This methodology has been experimentally verified at NASA Langley for structures ranging from laboratory coupons up to full-scale built-up structural components with single and multiple-site damage cracking. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic-plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intension of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels. In this regard, a series of fracture tests were conducted on curved aluminum-alloy integrally-stiffened panels. These curved panels were subjected to uniaxial tension and pressure loading. During the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (psi(sub c) using three-dimensional code (ZIP3D) and the plane-strain core height (h(sub c)) using two-dimensional code (STAGS).

Seshadri, B. R.↗

Practical guide to replica exchange transition interface sampling and forward flux sampling

Path sampling approaches have become invaluable tools to explore the mechanisms and dynamics of the so-called rare events that are characterized by transitions between metastable states separated by sizable free energy barriers. Their practical application, in particular to ever more complex molecular systems, is, however, not entirely trivial. Here, focusing on replica exchange transition interface sampling (RETIS) and forward flux sampling (FFS), we discuss a range of analysis tools that can be used to assess the quality and convergence of such simulations, which is crucial to obtain reliable results. The basic ideas of a step-wise evaluation are exemplified for the study of nucleation in several systems with different complexities, providing a general guide for the critical assessment of RETIS and FFS simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

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↗

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↗

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↗

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↗

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↗

Self-Consistent Simulation of the Magnetotail

Kinetic effects may be critically important at various stages of magnetotail evolution prior to and after the break up of substorms. The quasi-adiabaticity of ion dynamics results in the appearance of pronounced off diagonal terms in the pressure tensor (an effect that cannot be reproduced using the standard MHD (magnetohydrodynamic) approach). Only the appearance of these terms could support the formation of thin current sheets that recent experimental data have repeatedly indicated. We present in this paper the results of a self-consistent large-scale kinetic simulation of the Earth's magnetotail. We launch ions from the plasma mantle and self-consistently adjust the background magnetic field by computing the perturbations magnetic field from the particle current. We find that a thin current sheet with a half-thickness of approximately 0.3 R(sub E) quickly develops in the simulation. Also, the magnetotail evolves into a periodic state where the location of the X-line in the tail oscillates between x approximately 40 R(sub E) and x approximately 70 R(sub E) on a 5-minute time scale. Evidence of this periodic motion is also seen in the auroral precipitation profile of ions.

Peroomian, Vahe↗

A Machine Learning Examination of Hydroxyl Radical Differences Among Model Simulations for CCMI-1

Hydroxyl radical (OH) plays critical roles within the troposphere, such as determining the lifetime of methane (CH4), yet is challenging to model due to its fast cycling and dependence on a multitude of sources and sinks. As a result, the reasons for variations in OH and the resulting CH4 lifetime (TCH4), both between models and in time, are difficult to diagnose. We apply a neural network (NN) approach to address this issue within a group of models that participated in the Chemistry-Climate Model Initiative (CCMI). Analysis of the historical specified dynamics simulations performed for CCMI indicates that the primary drivers of TCH4 differences among ten models are the flux of UV light to the troposphere (indicated by the photolysis frequency JO1D) due mostly to clouds, mixing ratio of tropospheric ozone (O3), the abundance of nitrogen oxides (NOx≡NO+NO2), and details of the various chemical mechanisms that drive OH. Water vapor, carbon monoxide (CO), the ratio of NO:NOx, and formaldehyde (HCHO) explain moderate differences in TCH4, while isoprene, CH4, the photolysis frequency of NO2 by visible light (JNO2), overhead O3 column, and temperature account for little-to-no model variation in CH4. We also apply the NNs to analysis of temporal trends in OH from 1980 to 2015. All models that participated in the specified dynamics historical simulation for CCMI demonstrate a decline in CH4 during the analysed timeframe. The significant contributors to this trend, in order of importance, are tropospheric O3, JO1D, NOx, and H2O, with CO also causing substantial interannual variability in OH burden. Finally, the identified trends in TCH4 are compared to calculated trends in the tropospheric mean OH concentration from previous work, based on analysis of observations. The comparison reveals a robust result for the effect of rising water vapor on OH and CH4, imparting an increasing and decreasing trend of about 0.5% decade(exp -1), respectively. The responses due to NOx, O3 column, and temperature are also in reasonably good agreement between the two studies, though a discrepancy in the CH4 response highlights a need for further examination of the CH4 feedback on the abundance of OH.

Julie M Nicely↗

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↗

Numerical fracture simulation of compact and bend specimens

An elastic-plastic finite-element analysis with a critical crack-tip-opening displacement criterion was used to simulate fracture of various size compact and bend specimens made of HY-130 steel. From the calculated load-crack-extension and load-displacement curves, J-resistance (J-R) curves were determined by several methods. The simulated 3-R curves were insensitive to specimen size up to maximum load but were sensitive to specimen configuration for crack extensions greater than 10 percent of the initial uncracked ligament length.

Shivakumar, K. N.↗

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↗

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↗