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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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 OSTI2020

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
DOE OSTI2021

In Situ Machine Learning for Intelligent Data Capture on Exascale Platforms. Final Report

In many dynamic systems, interesting events occur locally in time and space. Examples of such systems include ignition events in combustion simulations, material fractures in mechanics simulations, and extreme weather events in climate simulations. Due to memory constraints and data I/O costs, current simulation workflows save data at regularly spaced time-steps, at a fixed rate determined before the start of the simulation. Often this mode of operation results in missed events of interest, necessitating a simulation restart from before an event occurred with more frequent data saves. This data saving workflow is grossly inefficient and is already a bottleneck in the computing process. We propose to develop machine learning algorithms that can detect when interesting dynamical events are occurring, triggering data saves. These machine learning algorithms will perform in situ anomaly detection to flag regions with different dynamical properties than those previously recorded. The adaptive data saves would be local in time and space to match the event of interest, thereby enabling a much more efficient workflow that will reduce data I/O costs and data storage memory requirements. The algorithms will be tested on two applications: auto-ignition simulations and climate simulations. A critical component of this project will be developing machine learning algorithms that can be deployed efficiently in situ on HPC platforms with out-of-the-box functionality. The development of in situ machine learning methods to detect anomalous events would enable a more efficient and effective workflow, in which all the relevant data are saved in a single simulation run, without re-starts or scientist intervention.

42 ENGINEERING↗
DOE OSTI2021

Representing methane emissions from wet tropical forest soils using microbial functional groups constrained by soil diffusivity

Tropical ecosystems contribute significantly to global emissions of methane (CH 4 ), and landscape topography influences the rate of CH 4 emissions from wet tropical forest soils. However, extreme events such as drought can alter normal topographic patterns of emissions. Here we explain the dynamics of CH 4 emissions during normal and drought conditions across a catena in the Luquillo Experimental Forest, Puerto Rico. Valley soils served as the major source of CH 4 emissions in a normal precipitation year (2016), but drought recovery in 2015 resulted in dramatic pulses in CH 4 emissions from all topographic positions. Geochemical parameters including (i) dissolved organic carbon (C), acetate, and soil pH and (ii) hydrological parameters like soil moisture and oxygen (O 2 ) concentrations varied across the catena. During the drought, soil moisture decreased in the slope and ridge, and O 2 concentrations increased in the valley. We simulated the dynamics of CH 4 emissions with the Microbial Model for Methane Dynamics-Dual Arrhenius and Michaelis–Menten (M3D-DAMM), which couples a microbial functional group CH 4 model with a diffusivity module for solute and gas transport within soil microsites. Contrasting patterns of soil moisture, O 2 , acetate, and associated changes in soil pH with topography regulated simulated CH 4 emissions, but emissions were also altered by rate-limited diffusion in soil microsites. Changes in simulated available substrate for CH 4 production (acetate, CO 2 , and H 2 ) and oxidation (O 2 and CH 4 ) increased the predicted biomass of methanotrophs during the drought event and methanogens during drought recovery, which in turn affected net emissions of CH 4 . A variance-based sensitivity analysis suggested that parameters related to aceticlastic methanogenesis and methanotrophy were most critical to simulate net CH 4 emissions. This study enhanced the predictive capability for CH 4 emissions associated with complex topography and drought in wet tropical forest soils.

54 ENVIRONMENTAL SCIENCES↗
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↗
DOE OSTI2020

Optimization of a Cyclone Using Multiphase Flow Computational Fluid Dynamics

The U.S. Department of Energy National Energy Technology Laboratory's (NETL) 50 kW th chemical looping reactor (CLR) has an underperforming cyclone, which was designed using empirical correlations. To improve the performance of this cyclone using computational fluid dynamics (CFD)-based modeling simulations, four critical design parameters including the vortex tube radius and length, barrel radius, and the inlet width and height were optimized. Here, NETL's open source multiphase flow with interphase exchange (MFiX) CFD code has been used to model a series of cyclones by systematically varying the geometric design parameters. To perform the optimization process, the surrogate modeling and sensitivity analysis followed by the optimization capability in nodeworks was used. The basic methodology for the process is to employ a statistical design of experiments (DOE) method to generate sampling simulations that fill the design space. Corresponding CFD models are then created, executed, and postprocessed. A response surface is created to characterize the relationship between input parameters and the quantities of interest (QoI). Finally, the CFD-surrogate is used by an optimization method to find the optimal design condition based on the objective and constraints prescribed. The resulting optimal cyclone has a larger diameter and longer vortex tube, a larger diameter barrel, and a taller and narrower solids inlet. The improved design has a predicted pressure drop 11 times lower than the original design while reducing the mass loss by a factor of 2.3.

42 ENGINEERING↗
DOE OSTI2020

Disorder effects in nitride semiconductors: impact on fundamental and device properties

Abstract Semiconductor structures used for fundamental or device applications most often incorporate alloy materials. In “usual” or “common” III–V alloys, based on the InGaAsP or InGaAlAs material systems, the effects of compositional disorder on the electronic properties can be treated in a perturbative approach. This is not the case in the more recent nitride-based GaInAlN alloys, where the potential changes associated with the various atoms induce strong localization effects, which cannot be described perturbatively. Since the early studies of these materials and devices, disorder effects have indeed been identified to play a major role in their properties. Although many studies have been performed on the structural characterization of materials, on intrinsic electronic localization properties, and on the impact of disorder on device operation, there are still many open questions on all these topics. Taking disorder into account also leads to unmanageable problems in simulations. As a prerequisite to address material and device simulations, a critical examination of experiments must be considered to ensure that one measures intrinsic parameters as these materials are difficult to grow with low defect densities. A specific property of nitride semiconductors that can obscure intrinsic properties is the strong spontaneous and piezoelectric fields. We outline in this review the remaining challenges faced when attempting to fully describe nitride-based material systems, taking the examples of LEDs. The objectives of a better understanding of disorder phenomena are to explain the hidden phenomena often forcing one to use ad hoc parameters, or additional poorly defined concepts, to make simulations agree with experiments. Finally, we describe a novel simulation tool based on a mathematical breakthrough to solve the Schrödinger equation in disordered potentials that facilitates 3D simulations that include alloy disorder.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2020

Data Federation Challenges in Remote Near-Real-Time Fusion Experiment Data Processing

Fusion energy experiments and simulations provide critical information needed to plan future fusion reactors. As next-generation devices like ITER move toward long-pulse experiments, analyses, including AI and ML, should be performed in a wide range of time and computing constraints, from near-real-time constraints, between-shot analysis, and to campaign-wide long-term analysis. However, the data volume, velocity, and variety make it extremely challenging for analyses using only local computational resources. Researchers need the ability to compose and execute workflows spanning edge resources to large-scale high-performance computing facilities.We present Delta, a system to address data analysis challenges, including AI/ML, in fusion science, by leveraging the ADIOS I/O library and middleware, to support executing science workflows over the wide area network for near-real-time streaming. We discuss the data federation challenges in performing remote workflows, focusing on on-going research work in (1) managing, reducing, and streaming data to minimize I/O and data movement overheads, (2) decompressing and reorganizing data for analysis, and (3) executing workflows for automated data analysis. We introduce examples for deep-learning based data analysis for the fusion domain and demonstrate how we use Delta to construct end-to-end workflows for a fusion device in Korea, connecting a remote DOE facility in the USA. The capability demonstrated by this project is the basis for improving the state of the art for near-real-time data federation amongst remote facilities.

Choi, Jong Youl↗
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 OSTI2020

Should Sea-Ice Modeling Tools Designed for Climate Research Be Used for Short-Term Forecasting?

In theory, the same sea-ice models could be used for both research and operations, but in practice, differences in scientific and software requirements and computational and human resources complicate the matter. Although sea-ice modeling tools developed for climate studies and other research applications produce output of interest to operational forecast users, such as ice motion, convergence, and internal ice pressure, the relevant spatial and temporal scales may not be sufficiently resolved. For instance, sea-ice research codes are typically run with horizontal resolution of more than 3 km, while mariners need information on scales less than 300 m. Certain sea-ice processes and coupled feedbacks that are critical to simulating the Earth system may not be relevant on these scales; and therefore, the most important model upgrades for improving sea-ice predictions might be made in the atmosphere and ocean components of coupled models or in their coupling mechanisms, rather than in the sea-ice model itself. This paper discusses some of the challenges in applying sea-ice modeling tools developed for research purposes for operational forecasting on short time scales, and highlights promising new directions in sea-ice modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The effect of local chemical ordering on Frank-Read source activation in a refractory multi-principal element alloy

In this work, we investigate the operation of Frank-Read (FR) sources in a refractory multi-principal element alloy (MPEA). Simulations of discrete dislocation motion in MPEAs is enabled by the development of a phase field dislocation dynamics model that treats the atomic-scale fluctuations in lattice energies across the glide plane, present due to local ordering in the chemical composition within the nominally random MPEA atomic structure. We consider, through simulation, a range of length scales over which ordering occurs, varying from short-range lengths, a few times dislocation core width, to long-range lengths, an order of magnitude longer than the core. Characteristic of this body-centered cubic MPEA, the simulations also include screw/edge character dependence in glide resistance, as informed by atomic scale simulation. The critical stresses to activate the source for the same source size are found to be statistically distributed, as a direct consequence of the underlying variation in lattice energy. Furthermore, analysis of critical state for activating edge and screw FR sources in the MPEA reveals that FR source operation occurs via a two-step mechanism, involving athermal kink-pair formation, unlike the conventional FR source operation in a material with no composition fluctuations. This mechanism lowers the average critical stress required to activate the FR source and causes the statistical dispersion in critical stress to depend on the range of composition ordering. More importantly, it leads to a more severe dependence of source strength on FR source length than predicted by line tension alone.

36 MATERIALS SCIENCE↗
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↗