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 325 records · Page 18

Integration Plan for NEAMS Structural Materials and Chemistry MSR Modeling

As part of its broad objective to support advanced reactor development, the US Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has a strong interest in developing simulation capabilities to address chemistry and corrosion issues for molten salt reactors (MSRs). To that end, the Structural Materials & Chemistry (SMC) technical area within the NEAMS program is developing a set of codes and databases to describe the behavior of molten salts and their interactions with structural materials. This document reviews the current activities with respect to MSR-relevant codes and is an effort to improve their integration in simulating MSR systems. Interoperability between the codes being developed is critical for the success of the NEAMS-developed MSR simulation capability. This is greatly facilitated by having codes based on the MOOSE framework when possible, and by developing interfaces to MOOSE-based tools for those tools that are not based on that framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Non-LWR Fuel Cycle Scenarios for SCALE and MELCOR Modeling Capability Demonstration

To assess the modeling and simulation capabilities for radionuclide characterization, criticality, and shielding in the nuclear fuel cycle of non-light–water reactor (LWR) technologies, the US Nuclear Regulatory Commission (NRC) initiated a collaborative project between the NRC, Sandia National Laboratories (SNL), and Oak Ridge National Laboratory (ORNL) with the goal to demonstrate capabilities of MELCOR and SCALE to calculate accident scenarios in all stages of the nuclear fuel cycle for relevant non-LWRs. The first project task was to develop representative nuclear fuel cycles and identify potential hazards and accident scenarios in the individual fuel cycle stages based on publicly available information. Because the nuclear fuel cycle is not established for any non-LWR concept, many assumptions were made, and it is anticipated that the details of the fuel cycles will eventually look different.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Increasing Al-Pair Abundance in SSZ-13 Zeolite via Zeolite Synthesis in the Presence of Alkaline Earth Metal Hydroxide Produces Hydrothermally Stable Co-, Cu- and Pd-SSZ-13 Materials

Replacing alkaline for alkaline-earth metal hydroxide in the synthesis gel during the synthesis of siliceous SSZ-13 zeolite (Si/Al~10) yields SSZ-13 with novel, advantageous properties. Its NH 4 -form ion-exchanges higher amount of isolated divalent M(II) ions than the conventional one: this is the consequence of an increased number of Al pairs in the structure induced by the +2 charge of Sr(II) cations in the synthesis gel that force two charge-compensating AlO 4 – motives to reside closer together. We characterize the +2 state of Co(II) ions in these materials with infra-red spectroscopy and X-ray absorption spectroscopy measurements and show their utility for NOx pollutant adsorption from ambient air: the ones derived from SSZ-13 with higher Al pair content contain more isolated cobalt(II) and, thus, perform better as ambient-air NOx adsorbers. Notably, Co(II)/SSZ-13 with an increased number of Al pairs is significantly more hydrothermally stable than its NaOH-derived analogue. Loading Pd(II) into Co-SSZ-13(Sr) produces an active NOx adsorber (PNA) material that can be used for NOx adsorption from simulated diesel engine exhaust. The critical issue for these applications is hydrothermal stability of Pd-zeolites. Pd/SSZ-13 synthesized in the presence of Sr(OH) 2 does not lose its PNA capacity after extremely harsh aging at 850 and 900 °C (10 h in 10% H 2 O/air flow) and loses only ~55% capacity after hydrothermal aging at 930 °C. This can be extended to other divalent metals for catalytic applications, such as copper: we show that Cu/SSZ-13 catalyst can survive hydrothermal aging at 920 °C without losing its catalytic properties, metal dispersion and crystalline structure. Thus, we provide a new, simple, and scalable strategy for making remarkably (hydro)thermally stable metal-zeolite materials/catalysts with a number of useful applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

H I Depletion Begins Well Beyond the Virial Radius: A FAST Stacking Study of 36 Galaxy Clusters to 5 × R 200

Abstract We present a stacking study of the neutral atomic hydrogen (H i ) content in and around 36 local galaxy clusters at z < 0.07, using a combination of the FAST All Sky H i survey (FASHI) and the extensive spectroscopic catalog mainly from the Dark Energy Spectroscopic Instrument (DESI). We employ spectral stacking techniques to probe the average H i mass and HI-to-stellar mass ratio ( M HI / M * ) for member galaxies down to stellar masses of M * ∼ 10 9 M ⊙ , spanning a projected cluster-centric distance of up to 5 R 200 . Our analysis reveals a pronounced environmental effect; both M HI and M HI / M * decrease steadily toward the cluster center, dropping by ∼0.5 dex on average from the outskirts to the core. Crucially, we find that M HI / M * of galaxies remain lower than the field galaxies even at the 5 R 200 . This provides direct, statistical evidence for substantial gas stripping and preprocessing in the cluster outskirts, likely occurring in infalling groups and large-scale filaments. By further splitting the sample by g − r color, we show that the H i deficiency persists at fixed galaxy color; even the bluest cluster members exhibit ∼0.5 dex lower M HI / M * than field galaxies of similar color, reflecting environmental effects on the cold gas reservoir prior to full optical transformation. The total H i mass within clusters and their outskirts agrees broadly with predictions from cosmological simulations. Our results underscore the critical role of the extended cluster environment in quenching galaxies by depleting their cold gas reservoirs well before they enter the dense cluster core.

Cheng, Cheng [Chinese Academy of Sciences South Am↗

BISON Methods and Contact Algorithms

This paper reviews five critical algorithms for BISON’s simulation of nuclear fuel: the burnup model, the coolant channel model, fuel cracking models, mechanical contact, and thermal contact.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Insights into the Physical-Chemical Properties of a CO2-Responsive Fracturing Fluid

Here we determine the phase behavior of StimuFrac, a CO2 responsive fracturing fluid, under geothermal wellbore conditions. StimuFrac is an aqueous poly(allylamine) fluid that crosslinks in the presence of CO2. StimuFrac significantly reduces the net pressure required to induce fractures, relative to other fracture fluids, and has potential to reduce water use. However, the phase behavior and equations of state to describe StimuFrac’s phase behavior remain unavailable. Here we determine the density and molar volume of the fluid as a function of geothermal relevant temperatures, pressures, and weight fractions of StimuFrac added. In general, experiments find that StimuFrac’s density decreases as temperature and pressure increase. Using these results, equations of saturated state and phase diagrams for different polymer concentrations are reported. These results are critical inputs for planned numerical simulation efforts.

Pease, Leonard F.↗

Interpreting Primal-Dual Algorithms for Constrained Multiagent Reinforcement Learning

Constrained multiagent reinforcement learning (C-MARL) is gaining importance as MARL algorithms find new applications in real-world systems ranging from energy systems to drone swarms. Most C-MARL algorithms use a primal-dual approach to enforce constraints through a penalty function added to the reward. In this paper, we study the structural effects of this penalty term on the MARL problem. First, we show that the standard practice of using the constraint function as the penalty leads to a weak notion of safety. However, by making simple modifications to the penalty term, we can enforce meaningful probabilistic (chance and conditional value at risk) constraints. Second, we quantify the effect of the penalty term on the value function, uncovering an improved value estimation procedure. We use these insights to propose a constrained multiagent advantage actor critic (C-MAA2C) algorithm. Simulations in a simple constrained multiagent environment affirm that our reinterpretation of the primal-dual method in terms of probabilistic constraints is effective, and that our proposed value estimate accelerates convergence to a safe joint policy.

chance constraints↗

Dynamics modeling of molten salt reactor with reduced and expanded representations of delayed neutron precursors

Molten salt reactors (MSRs) present unique challenges in dynamic behavior due to the mobility of their fuel. In these reactors, delayed neutron precursors (DNPs) drift with the fuel circulation through the primary loop. As a result, a fraction of DNPs decays outside the core, effectively reducing the available delayed neutron population for reactivity control. Consequently, precise modeling of the distribution and behavior of DNPs is critical for accurate reactor dynamics simulations. In this study, the System Dynamics Analysis Tool (SDAT) was used to simulate a thermal-spectrum MSR under steady-state conditions and following transients. The effects of using reduced and expanded representations of DNPs with fewer or more groups than the conventional 6-group model were investigated. Their impact on the simulated distribution of precursors in the primary loop, reactivity loss value, and reactor response to transients was analyzed. Simulation results showed that reduced models lead to the loss of the actual DNPs distribution data, resulting in less accurate estimates of reactivity loss. Reactor power predictions using these reduced models showed significant deviations compared to those using the conventional 6-group model in transient simulations. Expanded models offered a more accurate representation of the distribution of DNPs and reactivity loss estimates. Reactor power predictions using expanded models showed minimal deviation from the conventional 6-group model during the simulated transients.

analysis↗

Simulated Precipitation Diurnal Variation With a Deep Convective Closure Subject to Shallow Convection in Community Atmosphere Model Version 5 Coupled With CLUBB

In order to improve the physical consistency between shallow and deep convection, we modify the deep convective closure in the Community Atmosphere Model version 5 (CAM5) coupled with a third-order turbulence closure parameterization (i.e., Cloud Layers Unified by Binormals; CLUBB). The revised closure reserves a portion of the total convective available potential energy for shallow convection via utilizing the heating and moistening profiles from CLUBB to distribute moisture and energy between shallow and deep convection. Simulations at two resolutions (i.e. 2° and 0.5°, respectively) are conducted to investigate the impacts of convective closure on the simulated precipitation diurnal variations. Results from low-resolution simulations show that the revised closure suppresses deep convection until the lower troposphere is sufficiently moistened by shallow convection, which improves the precipitation diurnal variation simulations compared to the default closure, with the precipitation diurnal peak over tropical lands delayed from 12LST to 19LST. The revised closure better simulates the diurnal variations for precipitation over the Asian monsoon region, such as the delayed precipitation onset, but still fails to well capture the nighttime peak for precipitation there. This deficiency is alleviated to some extent when applying the revised closure in highresolution simulations, but nighttime precipitation is still underestimated probably because key processes responsible for nighttime convection are missing. Altogether, our results indicate that establishing the consistency between shallow and deep convection is critical for the precipitation diurnal cycle simulations.

54 ENVIRONMENTAL SCIENCES↗

Quantum criticality in coupled hybrid metal-semiconductor islands

We show that the combined effects of a dynamical Coulomb blockade and integer quantum Hall effect in a coupled hybrid metal-semiconductor setup provide a pathway for realizing resonant tunneling in Luttinger liquids. This hybrid setup can be brought to the quantum critical regime by varying gate voltages and contact resistances. We explore the nature of quantum criticality, Kondo effect, charge fractionalization, and transport in such a hybrid setup, and verify their robust non-Fermi-liquid behaviors. In conclusion, our work opens a promising route for quantum simulating exotic zero-temperature quantum critical phenomena associated with Luttinger-liquid physics in a nanoengineered electronic circuit with well-defined quantum Hall channels.

Coulomb blockade↗

Convergence in simulating global soil organic carbon by structurally different models after data assimilation

Abstract Current biogeochemical models produce carbon–climate feedback projections with large uncertainties, often attributed to their structural differences when simulating soil organic carbon (SOC) dynamics worldwide. However, choices of model parameter values that quantify the strength and represent properties of different soil carbon cycle processes could also contribute to model simulation uncertainties. Here, we demonstrate the critical role of using common observational data in reducing model uncertainty in estimates of global SOC storage. Two structurally different models featuring distinctive carbon pools, decomposition kinetics, and carbon transfer pathways simulate opposite global SOC distributions with their customary parameter values yet converge to similar results after being informed by the same global SOC database using a data assimilation approach. The converged spatial SOC simulations result from similar simulations in key model components such as carbon transfer efficiency, baseline decomposition rate, and environmental effects on carbon fluxes by these two models after data assimilation. Moreover, data assimilation results suggest equally effective simulations of SOC using models following either first‐order or Michaelis–Menten kinetics at the global scale. Nevertheless, a wider range of data with high‐quality control and assurance are needed to further constrain SOC dynamics simulations and reduce unconstrained parameters. New sets of data, such as microbial genomics‐function relationships, may also suggest novel structures to account for in future model development. Overall, our results highlight the importance of observational data in informing model development and constraining model predictions.

54 ENVIRONMENTAL SCIENCES↗

Spatiotemporal quenches for efficient critical ground state preparation in the two-dimensional transverse field Ising model

Quantum simulators have the potential to shed light on the study of quantum many-body systems and materials, offering unique insights into various quantum phenomena. Although adiabatic evolution has been conventionally employed for state preparation, it faces challenges when the system evolves too quickly or the coherence time is limited. In such cases, shortcuts to adiabaticity, such as spatiotemporal quenches, provide a promising alternative. This paper numerically investigates the application of spatiotemporal quenches in the two-dimensional transverse field Ising model with ferromagnetic interactions, focusing on the emergence of the ground state and its correlation properties at criticality when the gap vanishes. We demonstrate the effectiveness of these quenches in rapidly preparing ground states in critical systems. Our simulations reveal the existence of an optimal quench front velocity at the emergent speed of light, leading to minimal excitation energy density and correlation lengths of the order of finite system sizes we can simulate. These findings emphasize the potential of spatiotemporal quenches for efficient ground state preparation in quantum systems, with implications for the exploration of strongly correlated phases and programmable quantum computing.

2-dimensional systems↗

Critical pitting temperature of selective laser melted 316L stainless steel: A mechanistic approach

In this work, the critical pitting temperature (CPT) of selective laser melted (SLM) 316 L stainless steel in 1.0 M NaCl was measured and compared with a commercial wrought alloy. Potentiostatic measurements determined a mean CPT value of 16 ± 0.7 °C, 27.5 ± 0.8 °C and 31 ± 1 °C for the wrought alloy, the SLM alloy normal to the build direction and parallel to the build direction, respectively. The lead-in pencil electrode technique was used to study the pit chemistry of the two alloys and to explain the higher CPT values observed for the SLM alloy. A lower critical current density required for passivation in a simulated pit solution was measured for the SLM alloy. Moreover, the ratio of the critical concentration to saturated concentration of dissolving metal cations was found to be higher for the SLM alloy, which was related to its different salt film properties, possibly as a result of the SLM’s distinct microstructure.

36 MATERIALS SCIENCE↗

A Methodology for Estimating the Energy and Moisture Budget of the Convective Boundary Layer Using Continuous Ground-Based Infrared Spectrometer Observations

Abstract Land–atmosphere interactions play a critical role in both the atmospheric water and energy cycles. Changes in soil moisture and vegetation alter the partitioning of surface water and energy fluxes, influencing diurnal evolution of the planetary boundary layer (PBL). The mixing-diagram framework has proven useful in understanding the evolution of the heat and moisture budget within the convective boundary layer (CBL). We demonstrate that observations from the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains site provide all of the needed inputs needed for the mixing-diagram framework, allowing us to quantify the impact from the surface fluxes, advection, radiative heating, encroachment, and entrainment on the evolution of the CBL. Profiles of temperature and humidity retrieved from the ground-based infrared spectrometer [Atmospheric Emitted Radiance Interferometer (AERI)] are a critical component in this analysis. Large-eddy simulation results demonstrate that mean mixed-layer values derived are shown to be critical to close the energy and moisture budgets. A novel approach demonstrated here is the use of network of AERIs and Doppler lidars to quantify the advective fluxes of heat and moisture. The framework enables the estimation of the entrainment fluxes as a residual, providing a way to observe the entrainment fluxes without using multiple lidar systems. The high temporal resolution of the AERI observations enables the morning, midday, and afternoon evolution of the CBL to be quantified. This work provides a new way to use observations in this framework to evaluate weather and climate models. Significance Statement The energy and moisture budget of the planetary boundary layer (PBL) is influenced by multiple sources, and accurately representing this evolution in numerical models is critical for weather forecasts and climate predictions. The mixing-diagram approach, driven by profiling observations as illustrated here, provides a powerful way to quantify the contributions from each of these sources. In particular, the energy and moisture mixed into the PBL from above the PBL can be determined accurately from ground-based remote sensors using this approach.

54 ENVIRONMENTAL SCIENCES↗

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↗

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↗