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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS2008

Simulation of Ground Winds Time Series for the NASA Crew Launch Vehicle (CLV)

Simulation of wind time series based on power spectrum density (PSD) and spectral coherence models for ground wind turbulence is described. The wind models, originally developed for the Shuttle program, are based on wind measurements at the NASA 150-m meteorological tower at Cape Canaveral, FL. The current application is for the design and/or protection of the CLV from wind effects during on-pad exposure during periods from as long as days prior to launch, to seconds or minutes just prior to launch and seconds after launch. The evaluation of vehicle response to wind will influence the design and operation of constraint systems for support of the on-pad vehicle. Longitudinal and lateral wind component time series are simulated at critical vehicle locations. The PSD model for wind turbulence is a function of mean wind speed, elevation and temporal frequency. Integration of the PSD equation over a selected frequency range yields the variance of the time series to be simulated. The square root of the PSD defines a low-pass filter that is applied to adjust the components of the Fast Fourier Transform (FFT) of Gaussian white noise. The first simulated time series near the top of the launch vehicle is the inverse transform of the adjusted FFT. Simulation of the wind component time series at the nearest adjacent location (and all other succeeding next nearest locations) is based on a model for the coherence between winds at two locations as a function of frequency and separation distance, where the adjacent locations are separated vertically and/or horizontally. The coherence function is used to calculate a coherence weighted FFT of the wind at the next nearest location, given the FFT of the simulated time series at the previous location and the essentially incoherent FFT of the wind at the selected location derived a priori from the PSD model. The simulated time series at each adjacent location is the inverse Fourier transform of the coherence weighted FFT. For a selected design case, the equations, the process and the simulated time series at multiple vehicle stations are presented.

Adelfang, Stanley I.↗

DOE OSTI2023

Mechanistic Insights into Peptide Binding and Deactivation of an Adhesion G Protein-Coupled Receptor

Adhesion G protein-coupled receptors (ADGRGs) play critical roles in the reproductive, neurological, cardiovascular, and endocrine systems. In particular, ADGRG2 plays a significant role in Ewing sarcoma cell proliferation, parathyroid cell function, and male fertility. In 2022, a cryo-EM structure was reported for the active ADGRG2 bound by an optimized peptide agonist IP15 and the Gs protein. The IP15 peptide agonist was also modified to antagonists 4PH-E and 4PH-D with mutations of the 4PH residue to Glu and Asp, respectively. However, experimental structures of inactive antagonist-bound ADGRs remain to be resolved, and the activation mechanism of ADGRs such as ADGRG2 is poorly understood. Here, we applied Gaussian accelerated molecular dynamics (GaMD) simulations to probe conformational dynamics of the agonist- and antagonist-bound ADGRG2. By performing GaMD simulations, we were able to identify important low-energy conformations of ADGRG2 in the active, intermediate, and inactive states, as well as explore the binding conformations of each peptide. Moreover, our simulations revealed critical peptide-receptor residue interactions during the deactivation of ADGRG2. In conclusion, through GaMD simulations, we uncovered mechanistic insights into peptide (agonist and antagonist) binding and deactivation of the ADGRG2. These findings will potentially facilitate rational design of new peptide modulators of ADGRG2 and other ADGRs.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Equivalent Fundamental-Mode Source Simulations for Spherical Uranium and Plutonium Systems

It is important in neutron noise measurements and reactor startup to understand how the source distribution in a system compares to the equivalent fundamental-mode (EFM) source. Deviation from the fundamental mode source distribution increases as the system $\kappa$ eff increases. A correction factor, g *, is a measure of this deviation, and is the focus of this work. Two methods were utilized to simulate g *. These methods are compared to each other and (where applicable) previous results. The effects associated with nuclear material inner radius, materials of reflectors, and reflector thickness are also investigated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Phase-Field Modeling of Thermally-Grown Oxide and the Induced Damage Evolution in Environmental Barrier Coatings

The advent of next-generation hydrogen-based engines necessitates materials capable of withstanding temperatures beyond the reach of current superalloys. SiC-based ceramic matrix composites, augmented with environmental barrier coatings (EBCs), present a promising materials solution. Given the active search for effective and durable EBCs, there is a pressing need for modeling tools to understand and predict damage evolution in these materials to help accelerate their development. This study introduces a phase-field model (PFM) designed to simulate the critical role of thermally grown oxides (TGO) in the degradation and failure of EBCs. The model accounts for the severe volume expansion due to oxidation, alongside phase transformations and microstructural evolution during thermal cycling, offering a comprehensive view of the damage processes. Simulation results are validated against experimental findings reported in the literature, establishing the model's potential as a significant tool for understanding and improving the resilience of EBCs in cyclic oxidative environments.

Cheng, Tianle↗

NASA NTRS1991

Bench-level characterization of a CMOS standard-cell D-latch using alpha-particle sensitive test circuits

A methodology is described for predicting the SEU susceptibility of a standard-cell D-latch using an alpha-particle sensitive SRAM, SPICE critical charge simulation results, and alpha-particle interaction physics. Measurements were made on a 1.6-micron n-well CMOS 4-kb test SRAM irradiated with an Am-241 alpha-particle source. A collection depth of 6.09 micron was determined using these results and TRIM computer code. Using this collection depth and SPICE derived critical charge results on the latch design, an LET threshold of 34 MeV sq cm/mg was predicted. Heavy ion tests were then performed on the latch and an LET threshold of 41 MeV sq cm/mg was determined.

Blaes, B. R.↗

NASA NTRSDate not supplied

Digital Lunar Exploration Sites (DLES)

It has been almost 50 years since humans last set foot on the Moon. With NASA’s Artemis program, the United States and its international and commercial partners are embarking on a new endeavor to explore the lunar surface. Before we return, we will have simulated every aspect of these future missions. Many of these models and simulations (M&S) will rely on well-known and commonly-used technologies, some of which trace their origins back to the Apollo program. However, M&S has advanced significantly, as have the underlying computational capabilities. As a result, we are able to model many more aspects of the Artemis vehicles and support systems with significantly improved detail and confidence. Accurately and efficiently modeling the lunar environment will be critical to simulating the Artemis elements and mission activities. This includes characterizing and modeling lunar topography, smaller craters, exposed surface rocks, lunar regolith, surface lighting, and ambient thermal environment. These are all necessary for understanding fundamental behaviors and performance of vehicles and support systems in the lunar environment and are often determining factors in the selection of exploration sites and defining mission profiles. The Astromaterials Research and Exploration Sciences (ARES) team in the Exploration Integration and Science Directorate at NASA’s Johnson Space Center (JSC) and the NASA Exploration Systems Simulations (NExSyS) team in the Simulation and Graphics Branch at JSC are developing and maintaining the Digital Lunar Exploration Sites (DLES) data, documentation, and software packages. This data is being fed directly into a diverse collection of graphics and simulation environments, where it is used to construct the closest known truth for numerous potential Lunar South Pole landing sites. DLES is available to projects across NASA and particularly the Artemis program to support coordinated digital representations of the lunar environment. This paper describes the fundamental need for DLES, the science data sets that are going into DLES, some of the processes used to integrate this data into DLES products, the basic products that constitute DLES, and some examples of DLES in use.

Lunar↗

DOE OSTI2026

Design of Controller Hardware-In-the-Loop Model of Microgrid with Modular Building Blocks and Automated Design Script

The scalability of controller hardware-in-the-loop (CHIL) simulation is critical for validating control coordination and energy management in microgrids with distributed energy resources, especially as these modern systems become more complex and decentralized. This paper presents a CHIL modeling methodology that combines modular building blocks with an automated design script to streamline the development of high-fidelity microgrid models. Standardized subsystem templates for resources, converters, and buses are integrated with a Python-based script that compiles structured JSON configuration files into simulation-ready initialization code. The proposed approach reduces development time, improves model consistency, and enhances simulation fidelity. The methodology is validated on a Typhoon HIL604 platform and is broadly applicable to real-time simulation of complex, networked microgrid systems. This framework establishes a foundation for automated, scalable CHIL validation and accelerates the design of next-generation distributed energy systems.

Kim, Namwon [ORNL] (ORCID:0000000200438489)↗

DOE OSTI2025

A round bobbin critical current measurement and thermal runaway simulation of REBCO coated conductors

A round bobbin test has been widely used to measure critical current Ic of practical superconductors including Nb–Ti, Nb3Sn, and Bi-2212 round strands at 4.2 K and in high magnetic fields but rarely used for rare earth-barium-copper-oxide (REBCO) coated conductor tapes. Here, we applies this method to REBCO tapes and test their Ic at 77 K, self-field, and 4.2 K with a magnetic field of 14 T applied parallel to the tape. The sample carries very high current densities; it is 1.3 m long and has a dense array of voltage taps that allows probing localized thermal runaways and correlation with Ic variations along the length. At 77 K, the Ic values of all sections were measured and found to be rather uniform. However, at 4.2 K, we found that a 3 cm section has 90% of the overall voltage drop and it went into thermal runaways upon further increasing current. Thus, at 4.2 K, we were prevented by localized thermal runaways and unable to determine Ic of most of the REBCO conductor sections; this behavior is different from Bi-2212 and Bi-2223 for which Ic of all turns of a round barrel sample can be determined. We simulate REBCO tapes during current ramping to a thermal runaway and analyze impacts of various factors including cooling, conductor n-values, and thickness of stabilizer on the thermal runaway current Iq. The combined experiments and simulations identify n-value as a crucial factor that leads to the observed differences between Bi-2212 and REBCO. Together with cooling and the standard deviation of Ic variations along the length, the n-value determines whether Iq is above or below the average Ic of a long piece of conductor and their ratio Iq/Ic,ave, greatly affecting quench detection and thermal runaway prevention strategies of Bi-2212 and REBCO magnets.

Mosat, Marek↗

DOE OSTI2025

Viscosity Measurements in Extreme Conditions

Radiation-hydrodynamics simulations are critical to the NNSA complex for understanding high energy density physics and for addressing problems in national security science. However, current rad-hydro codes typically do not account for viscosity, which generally leads to large uncertainties in problems involving mixing of materials under shock loading. This report details progress made towards experimental measurements of fluid viscosity at high pressures and temperatures achieved by shock wave compression techniques, as well as complementary material model development and hydrodynamic simulations.

36 MATERIALS SCIENCE↗

NASA NTRS2017

Direct Estimation of the Global Distribution of Vertical Velocity Within Cirrus Clouds

Cirrus clouds determine the radiative balance of the upper troposphere and the transport of water vapor across the tropopause.The representation of vertical wind velocity, W, in atmospheric models constitutes the largest source of uncertainty in the calculation of the cirrus formation rate. Using global atmospheric simulations with a spatial resolution of 7 kilometers we obtain for the first time a direct estimate of the distribution of W at the scale relevant for cirrus formation, validated against long-term observations at two different ground sites. The standard deviation in W, sigma (sub W), varies widely over the globe with the highest values resulting from orographic uplift and convection, and the lowest occurring in the Arctic. Globally about 90 of the simulated sigma (sub W) values are below 0.1 meters per second and about one in 10 (sup 4) cloud formation events occur in environments with sigma (sub W) greater than 0.8 meters per second. Combining our estimate with reanalysis products and an advanced cloud formation scheme results in lower homogeneous ice nucleation frequency than previously reported, and a decreasing average ice crystal concentration with decreasing temperature. These features are in agreement with observations and suggest that the correct parameterization of sigma (sub W) is critical to simulate realistic cirrus properties.

Radiative Balance↗

NASA NTRS1984

Chip level modeling of LSI devices

The advent of Very Large Scale Integration (VLSI) technology has rendered the gate level model impractical for many simulation activities critical to the design automation process. As an alternative, an approach to the modeling of VLSI devices at the chip level is described, including the specification of modeling language constructs important to the modeling process. A model structure is presented in which models of the LSI devices are constructed as single entities. The modeling structure is two layered. The functional layer in this structure is used to model the input/output response of the LSI chip. A second layer, the fault mapping layer, is added, if fault simulations are required, in order to map the effects of hardware faults onto the functional layer. Modeling examples for each layer are presented. Fault modeling at the chip level is described. Approaches to realistic functional fault selection and defining fault coverage for functional faults are given. Application of the modeling techniques to single chip and bit slice microprocessors is discussed.

Armstrong, J. R.↗

DOE OSTI2021

Preparing an Incompressible-Flow Fluid Dynamics Code for Exascale-Class Wind Energy Simulations

The U.S. Department of Energy has identified exascale-class wind farm simulation as critical to wind energy scientific discovery. A primary objective of the ExaWind project is to build high-performance, predictive computational fluid dynamics (CFD) tools that satisfy these modeling needs. GPU accelerators will serve as the computational thoroughbreds of next-generation, exascale-class supercomputers. Here, we report on our efforts in preparing the ExaWind unstructured mesh solver, Nalu-Wind, for exascale-class machines. For computing at this scale, a simple port of the incompressible-flow algorithms to GPUs is insufficient. To achieve high performance, one needs novel algorithms that are application aware, memory efficient, and optimized for the latest-generation GPU devices. The result of our efforts are unstructured-mesh simulations of wind turbines that can effectively leverage thousands of GPUs. In particular, we demonstrate a first-of-its-kind, incompressible-flow simulation using Algebraic Multigrid solvers that strong scales to more than 4000 GPUs on the Summit supercomputer.

algebraic multigrid↗

DOE OSTI2022

Metastable doubly threaded [3]rotaxanes with a large macrocycle

Ring size is a critically important parameter in many interlocked molecules as it directly impacts many of the unique molecular motions that they exhibit. Reported herein are studies using one of the largest macrocycles reported to date to synthesize doubly threaded [3]rotaxanes. A large ditopic 46 atom macrocycle containing two 2,6-bis(N-alkyl-benzimidazolyl)pyridine ligands has been used to synthesize several metastable doubly threaded [3]rotaxanes in high yield (65–75% isolated) via metal templating. Macrocycle and linear thread components were synthesized and self-assembled upon addition of iron(II) ions to form the doubly threaded pseudo[3]rotaxanes that could be subsequently stoppered using azide–alkyne cycloaddition chemistry. Following demetallation with base, these doubly threaded [3]rotaxanes were fully characterized utilizing a variety of NMR spectroscopy, mass spectrometry, size-exclusion chromatography, and all-atom simulation techniques. Critical to the success of accessing a metastable [3]rotaxane with such a large macrocycle was the nature of the stopper group employed. By varying the size of the stopper group it was possible to access metastable [3]rotaxanes with stabilities in deuterated chloroform ranging from a half-life of <1 minute to ca. 6 months at room temperature potentially opening the door to interlocked materials with controllable degradation rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Non-equilibrium critical scaling and universality in a quantum simulator

Universality and scaling laws are hallmarks of equilibrium phase transitions and critical phenomena. However, extending these concepts to non-equilibrium systems is an outstanding challenge. Despite recent progress in the study of dynamical phases, the universality classes and scaling laws for non-equilibrium phenomena are far less understood than those in equilibrium. In this work, using a trapped-ion quantum simulator with single-spin resolution, we investigate the non-equilibrium nature of critical fluctuations following a quantum quench to the critical point. We probe the scaling of spin fluctuations after a series of quenches to the critical Hamiltonian of a long-range Ising model. With systems of up to 50 spins, we show that the amplitude and timescale of the post-quench fluctuations scale with system size with distinct universal critical exponents, depending on the quench protocol. While a generic quench can lead to thermal critical behavior, we find that a second quench from one critical state to another (i.e. a double quench) results in a new universal non-equilibrium behavior, identified by a set of critical exponents distinct from their equilibrium counterparts. Our results demonstrate the ability of quantum simulators to explore universal scaling beyond equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Temperature effects on adsorption and capillarity water retention mechanisms in constrained unsaturated soils

Abstract This paper focuses on the impact of elevated temperatures on the adsorptive and capillarity water retention mechanisms of unsaturated soils under constrained (constant volume) conditions. This topic is critical for simulating the thermo-hydraulic behavior of soils in hydrogeological or geotechnical applications, including climate change effects on near surface soils, energy piles or soil borehole thermal energy storage systems in unsaturated soil layers, and buffers for geological nuclear waste repositories. A nonisothermal soil water retention curve (SWRC) that separately considers the temperature-dependency of the key parameters governing adsorptive and capillarity water retention mechanisms and soil physical parameters (e.g., surface tension, contact angle, adsorption capacity, cation exchange capacity, mean cavitation suction, air entry value and equilibrium film thickness) was developed to provide insights into the impact of temperature on water retention over the full suction range. The nonisothermal SWRC was validated using experimental data on high plasticity clays, with a good prediction of temperature effects on adsorption and capillarity water retention mechanisms in constrained unsaturated soils.

Lu, Yu↗

DOE OSTI2022

Multiscale petrophysical characterization and flow unit classification of the Minnelusa eolian sandstones

Integration of petrophysical and geological information is critical to simulation of subsurface carbon storage (GCS). In this sense, two depositional facies were identified from the core description and well-log interpretation, namely massive (MS) and cross-bedded (CB) facies groups. Additionally, pore-scale characteristics were studied by a combination of techniques, e.g. Nuclear Magnetic Resonance (NMR) and mercury intrusion capillary pressure (MICP). Scanning electron microscope (SEM) and petrographic analyses show that the pore structure is dominantly controlled by the depositional environment and dolomite cementation. NMR-T 2 distributions of MS and CB facies show triple and quadruple modes, respectively. In addition, MICP of high- and low-permeability MS facies samples, and their CB facies group mixtures were collected. The MS sample pore-throat size distribution is uni-modal, while the triple-modal characteristic of the mixtures indicates heterogeneous pore structures at the sub-core scale for CB facies. The reliably estimates of porosity and permeability for both facies groups via NMR techniques and the MLR (Multiple Linear Regression) approach demonstrate the applicability of these techniques to eolian sandstone. Moreover, irreducible water saturation via the T 2 -cutoff method correlates strongly with T 2LM instead of porosity. Finally, the rock quality index and flow zone indicator were calculated based on Combinable Magnetic Resonance (CMR) log interpretations. This provides direct connection to properties measured in the well. Four flow units were classified for both facies groups. Finally, results show that better reservoir quality with significant heterogeneities is observed in the CB facies. This study highlights the importance integrating multiscale petrophysical properties including facies, pore architecture and diagenesis analysis with core- to log-scale property characterization. The results herein validate our reservoir characterization and flow unit classification in eolian reservoirs.

58 GEOSCIENCES↗

DOE OSTI2025

Soot formation and precursor chemistry in Counterflow flames of aviation fuel surrogates

To meet market demands, the aviation sector is interested in utilizing drop-in Synthetic Aviation Turbine Fuels (SATF), either as neat fuels or in blends with conventional Jet A. SATF currently approved in standard specifications may have lower aromatic content with significant fractions of normal, branched, and cyclo-alkanes. Fundamental studies on soot formation from aviation fuels (Jet A, SATF) and their surrogate components are essential to understand how fuel composition influences soot and soot precursor formation. Here, this study reports new measurements of polycyclic aromatic hydrocarbons (PAH) and soot in counterflow diffusion flames (CDFs) of aviation fuel surrogates. Both intrusive and non-intrusive diagnostics are employed to determine the profiles of temperature, gas phase species, PAHs (up to C16), and soot volume fraction (SVF) in CDFs of iso-octane and surrogate mixtures. These measurements shed light on the transition of soot precursors to primary soot particles. In addition to serving as a common surrogate component in Jet A surrogate mixtures, iso-octane is a template species for larger, less volatile branched alkanes found in SATF mixtures. The newly developed Lawrence Livermore National Laboratory (LLNL) PAH and soot model successfully captures temperature, precursor species, and SVF profiles for the mixtures and conditions discussed in this work. Finally, a high-fidelity surrogate for Jet A is proposed that matches targeted physical and chemical properties well, while leveraging the wide range of candidate fuel molecules available in the LLNL detailed chemical model. The proposed surrogate formulation is validated against newly acquired measurements of the surrogate and literature measurements of Jet A. These new experiments and simulations provide critical insights into the PAH and soot formation from aviation fuels. Reaction pathways which require further investigation are highlighted, such that future work may bridge the remaining quantitative gaps in predicting soot formation from aviation fuel surrogates and surrogate components.

Aviation fuels↗

DOE OSTI2022

Insights into Sorption and Molecular Transport of Aqueous Glucose into Zeolite Nanopores

Liquid-phase heterogeneous catalysis using zeolites is important for biomass conversion to fuels and chemicals. There is a substantial body of work on gas-phase sorption in zeolites with different topologies; however, studies investigating the diffusion of complex molecules in liquid medium into zeolitic nanopores are scarce. Here, we present a molecular dynamics study to understand the sorption and diffusion of aqueous β-d-glucose into β-zeolite silicate at T = 395 K and P = 1 bar. Through 2-μs-long molecular dynamics trajectories, we reveal the role of the solvent, the kinetics of the pore filling, and the effect of the water model on these properties. We find that the glucose and water loading is a function of the initial glucose concentration. Although the glucose concentration increases monotonically with the initial glucose concentration, the water loading exhibits a nonmonotonic behavior. At the highest initial concentration (~20 wt %), we find that the equilibrium loading of glucose is approximately five molecules per unit cell and displays a weak dependence on the water model. Glucose molecules follow a single-file diffusion in the nanopores due to confinement. Here, the dynamics of glucose and water molecules slows significantly at the interface. The average residence time for glucose molecules is an order of magnitude larger than that in the bulk solution, while it is about twice as large for the water molecules. Our simulations reveal critical molecular details of the glucose molecule’s local environment inside the zeolite pore relevant to catalytic conversion of biomass to valuable chemicals.

09 BIOMASS FUELS↗