SEARCH · Search NASA

Results for “lateral mixing”

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

DOE OSTI2025

Electronic effects in radiation-induced collision cascades in nickel

The accurate treatment of electronic effects in multi-million-atom simulations of radiation-induced collision cascades is crucial for reliable predictions of primary radiation damage. In this work, we explore the fidelity of a recently developed two-temperature molecular dynamics model implementing an electron density-dependent coupling of electronic and atomic subsystems for cascade simulations in nickel. We show that the parameter-free model realistically captures the instantaneous energy losses during all stages of the highly nonequilibrium cascade process. Our simulations predict two distinct coupling regimes, corresponding to the rapid energy losses through electronic stopping in the early stages of the cascade and to the slow equilibration through the electron-phonon coupling mechanism in the later stages, without the use of separate models or coupling terms. The intermediate stage of the cascade dynamics displays a complex energy transfer between the subsystems, which cannot be interpreted by comparison to either electronic stopping or electron-phonon coupling theories. We therefore compare the predicted atomic mixing, which is sensitive to the energy losses during the intermediate cascade stage, with experimental ion beam mixing measurements. We find good agreement with the experiments, validating the coupling model for the intermediate stage of the cascade. Predictions of final defect numbers and cluster sizes are found in line with the results from conventional electronic stopping-based methods, while significantly reducing the theoretical uncertainty in the outcomes of conventional models stemming from arbitrary choices of thresholds for different coupling terms. Our results represent a notable improvement in cascade damage predictions in nickel, providing validation of the electron density-dependent coupling model for radiation damage simulations in general. The results lead us to propose an interpretation of the electronic energy losses in the intermediate regime of velocities, where we find an effectively nonlinear dissipation.

Crystal defects↗

DOE OSTI2022

M-Star® Software Test and Verification for Impeller Mixing in a Tank

Before a sample of the Slurry Mix Evaporator (SME) can be taken, the SME product sampling procedure requires that the agitator power be stable between 20 and 30 kW for at least one hour. The SME transfer to Melter Feed Tank (MFT) procedure also requires the SME agitator power be stabilized between 20 and 30 kW prior to transfer. These requirements are specified to ensure samples are homogeneous, as discussed in the Waste Form Qualification Report. During SME Batch 804, the agitator power dropped to 19 kW and struggled to achieve and maintain 20 kW. It was later discovered that the cause of the power drop was due to one of the bottom blades of the agitator breaking off. Both the sample and the transfer occurred without the power stabilizing between 20 and 30 kW. Therefore, the quality of SME Batch 804 is indeterminate, and homogeneity was questionable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

Basin-scale analysis of Mokelumne River Formation for multi-well CO 2 injection

Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.

Basin-scale↗

DOE OSTI2024

Cumulonimbus Clouds Convert a Smaller Fraction of CAPE into Kinetic Energy in a Warmer Atmosphere

Abstract This study investigates how entrainment’s diluting effect on cumulonimbus updraft buoyancy is affected by the temperature of the troposphere, which is expected to increase by the end of the century. A parcel model framework is constructed that allows for independent variations in the temperature ( T ), the entrainment rate ε , the free-tropospheric relative humidity (RH), and the convective available potential energy (CAPE). Using this framework, dilution of buoyancy is evaluated with T and RH independently varied and with CAPE either held constant or increased with temperature. When CAPE is held constant, buoyancy decreases as T increases, with parcels in warmer environments realizing substantially smaller fractions of their CAPE as kinetic energy (KE). This occurs because the increased moisture difference between an updraft and its surroundings at warmer temperatures drives greater updraft dilution. Similar results are found in midlatitude and tropical conditions when CAPE is increased with temperature. With the expected 6%–7% increase in CAPE per kelvin of warming, KE only increases at 2%–4% K −1 in narrow updrafts but tracks more closely with CAPE at 4%–6% in wider updrafts. Interestingly, the rate of increase in the KE with T becomes larger than that of CAPE when the later quantity increases at more than 10% K −1 . These findings emphasize the importance of considering entrainment in studies of moist convection’s response to climate change, as the entrainment-driven dilution of buoyancy may partially counteract the influence of increases in CAPE on updraft intensity. Significance Statement Cumulonimbus clouds mix air with their surrounding environment through a process called entrainment, which controls how efficiently environmental energy is converted into upward speed in thunderstorm updrafts. Our research shows that warmer temperatures will exacerbate the moisture difference between cumulonimbus updrafts and their surroundings, leading to greater mixing and less efficient conversion of environmental energy into updraft speeds. This effect should be considered in future research that investigates how climate change will affect cumulonimbus clouds.

Peters, John M.↗

DOE OSTI2020

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗

DOE OSTI2025

Can increasing the size and flexibility of a molecule reduce decoherence and prolong charge migration?

Coherent superposition of electronic states, created by ionizing a molecule, can initiate ultrafast dynamics of the electron density. Correlation between nuclear and electron motions, however, typically dissipates the electronic coherence in only a few femtoseconds, especially in larger and more flexible molecules. We, therefore, use ab initio semiclassical dynamics to study decoherence in a sequence of analogous organic molecules of increasing size and find, surprisingly, that extending the carbon skeleton in propynal analogs slows down decoherence and prolongs charge migration. To elucidate this observation, we decompose the overall decoherence into contributions from individual vibrational modes and show that: 1) The initial decay of electronic coherence is caused by high- and intermediate-frequency vibrations via momentum separation of nuclear wavepackets evolving on different electronic surfaces. 2) At later times, the coherence disappears completely due to the increasing position separation in the low-frequency modes. 3) In agreement with another study, we observe that only normal modes that preserve the symmetry of the molecule induce decoherence. All together, we justify the enhanced charge migration by a combination of increased hole-mixing and the disappearance of decoherence contributions from specific vibrational modes—CO stretching in butynal and various H rockings in pentynal.

Science & Technology - Other Topics↗

DOE OSTI2020

A New Estimator to Correct for Bias from Tag Rate Expansion on Natural-Origin Fish Attributes in Mixed-Stock Analysis Using Parentage-Based Tagging

Abstract In fisheries where hatchery- and natural-origin conspecifics occur as mixed stocks, it is often important to estimate both the natural-origin proportion of the mixture and the composition of attributes within the natural-origin portion (e.g., genetic stock, sex, and age-class). These estimates are facilitated by parentage-based tagging, which allows large numbers of hatchery fish to be efficiently tagged and later identified. When tag rates are less than 100%, the untagged fish in the mixture include both untagged hatchery fish and natural-origin fish. Unbiased estimation of the abundance and attribute composition of the natural-origin portion requires an estimator that accounts for tag rates. Two estimators are described: one “accounting-style” estimator similar to previously described approaches and one maximum likelihood method. These estimators were evaluated and compared using simulations mimicking estimation of the composition of fish migrating past a dam. The two estimators performed similarly at reducing bias when tag rates were high, but the maximum likelihood method had smaller mean square error when tag rates were low. We provide an R package to allow usage of this estimator in a wide variety of fisheries applications.

Delomas, Thomas A. (ORCID:000000015154759X)↗

DOE OSTI2026

An Optimized Parameterization of Sub‐Grid Scale Advection for Convection Permitting Models

Convection‐permitting models (CPMs) explicitly resolve deep convection yet under‐resolve the organized lateral exchanges among drafts and their environment that control entrainment/detrainment, precipitation efficiency, and mesoscale structure. In this work, we introduce the Optimized Advection Scheme (OAS), which introduces a small rotation of the Cartesian frame of reference for the horizontal winds relative to other variables used in advection that induces cross‐gradient transport to mimic under‐resolved convective mixing. The rotation angle is selected to minimize the Kullback–Leibler divergence between the simulated and satellite observed precipitation intensity distributions, yielding a physically consistent perturbation that is computationally inexpensive and portable. Optimized Advection Scheme is implemented in WRF and evaluated over Amazon (April 2014). It shifts precipitation–precipitable‐water joint distributions toward lighter rain, reduces overly intense rates, and improves mesoscale convective system (MCS) lifetime and propagation. Mechanistically, the added cross‐gradient transport promotes convective detrainment and environmental mixing, which cools and moistens the mid‐troposphere, weakens downward momentum transport, alleviates excessive downwelling shortwave biases, and warms the surface temperature. The optimized rotation angle yields comparable improvements at 4‐km and 1‐km grid spacing, demonstrating resolution‐independent benefits across the CPM gray zone. By targeting the dynamical root of under‐mixed convective circulations, rather than tuning model microphysics or closures, OAS delivers robust, scale‐aware improvements in precipitation statistics, cloud vertical structure, and characteristics of MCS (MCSs), offering a practical pathway to more reliable CPM simulations for weather and climate applications.

CPM↗

DOE OSTI2022

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Deterministic Fabrication of Large-Area, High-Crystallinity Oxide Moiré Superlattices

Oxide twistronics extends moiré engineering beyond van der Waals materials, offering a promising platform for accessing emergent interfacial phenomena arising from the strong coupling of lattice, charge, and orbital degrees of freedom in complex oxides. However, deterministic fabrication of high-crystallinity oxide moiré superlattices over large lateral dimensions remains challenging due to the three-dimensional bonding network of oxides. Here, we demonstrate a scalable, generalized fabrication strategy that enables the formation of high-crystallinity oxide moiré superlattices with clean, chemically bonded interfaces and precisely controlled twist angles down to nominal values of 0.1°, achieving subdegree twist-angle accuracy across large contiguous lateral dimensions approaching the millimeter scale. Using NaNbO3 as a model system, we show that the resulting interlayer coupling drives pronounced structural reconstruction that modifies both the phase structure and ferroelectric domain configuration. Synchrotron-based X-ray 3D reciprocal space mapping reveals the emergence of a single-phase state in twisted bilayers, in contrast to the mixed-phase structure observed in single-layer membranes prior to twist assembly. The structural signatures are further consistent with gradual lattice rotation distributed along the thickness direction that may accommodate interfacial shear strain, distinct from reconstruction observed in van der Waals moiré systems which primarily occurs through in-plane stacking rearrangement. This collective lattice response is correlated with twist-dependent nanoscale electromechanical modulations observed by piezoresponse force microscopy. These results establish a scalable materials platform for oxide twistronics and support the implementation of twist-engineered functionalities in practical, macroscale device architectures.

Ghanbari, Reza [North Carolina State University (N↗

DOE OSTI2020

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗

DOE OSTI2021

Amount and reactivity of dissolved organic matter export are affected by land cover change from old‐growth to second‐growth forests in headwater ecosystems

Abstract Headwater forest ecosystems of the western USA generate a large portion of the dissolved organic matter (DOM) transported in streams across North America. Land cover changes that alter forest structure and species composition affect the quantity and composition of DOM transferred to aquatic ecosystems. Clear‐cut harvesting affects ~1% of the forest area of North America annually, leaving most western forests in varying stages of regrowth and the total area of old‐growth forest is decreasing. The consequences of this widespread management practice on watershed carbon cycling remain unknown. We investigated the role of land cover change, because of clear‐cut harvesting, from mixed‐species old‐growth to lodgepole pine‐dominated second‐growth forest on the character and reactivity of hillslope DOM exports. We evaluated inputs of DOM from litter leachates and export of DOM collected at the base of trenched hillslopes during a 3‐year period (2016–2018) at the Fraser Experimental Forest in north‐central Colorado, USA. Dissolved organic carbon and total dissolved nitrogen were higher in lateral subsurface flow draining old‐ versus second‐growth forest. Fluorescence spectroscopy showed that the DOM exported from the old‐growth forest was more heterogeneous and aromatic and that proteinaceous, microbially processed DOM components were more prevalent in the second‐growth forest. Biological oxygen demand assays revealed much lower microbial metabolism of DOM in litter leachate and subsurface exports from the old‐growth forest relative to second growth. Old‐growth and second‐growth forests are co‐mingled in managed ecosystems, and our findings demonstrate that land cover change from a mixture of conifer species to lodgepole pine dominance influences DOM inputs that can increase the reactivity of DOM transferred from terrestrial to aquatic ecosystems.

Fegel, Timothy S.↗

DOE OSTI2021

Amount and reactivity of dissolved organic matter export are affected by land cover change from old-growth to second-growth forests in headwater ecosystems

Headwater forest ecosystems of the western USA generate a large portion of the dissolved organic matter (DOM) transported in streams across North America. Land cover changes that alter forest structure and species composition affect the quantity and composition of DOM transferred to aquatic ecosystems. Clear-cut harvesting affects ~1% of the forest area of North America annually, leaving most western forests in varying stages of regrowth and the total area of old-growth forest is decreasing. The consequences of this widespread management practice on watershed carbon cycling remain unknown. We investigated the role of land cover change, because of clear-cut harvesting, from mixed-species old-growth to lodgepole pine-dominated second-growth forest on the character and reactivity of hillslope DOM exports. We evaluated inputs of DOM from litter leachates and export of DOM collected at the base of trenched hillslopes during a 3-year period (2016–2018) at the Fraser Experimental Forest in north-central Colorado, USA. Dissolved organic carbon and total dissolved nitrogen were higher in lateral subsurface flow draining old- versus second-growth forest. Fluorescence spectroscopy showed that the DOM exported from the old-growth forest was more heterogeneous and aromatic and that proteinaceous, microbially processed DOM components were more prevalent in the second-growth forest. Biological oxygen demand assays revealed much lower microbial metabolism of DOM in litter leachate and subsurface exports from the old-growth forest relative to second growth. Old-growth and second-growth forests are co-mingled in managed ecosystems, and our findings demonstrate that land cover change from a mixture of conifer species to lodgepole pine dominance influences DOM inputs that can increase the reactivity of DOM transferred from terrestrial to aquatic ecosystems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Ultrathin ferroic HfO 2 –ZrO 2 superlattice gate stack for advanced transistors

With the scaling of lateral dimensions in advanced transistors, an increased gate capacitance is desirable both to retain the control of the gate electrode over the channel and to reduce the operating voltage. This led to a fundamental change in the gate stack in 2008, the incorporation of high-dielectric-constant HfO 2 , which remains the material of choice to date. Here we report HfO 2 –ZrO 2 superlattice heterostructures as a gate stack, stabilized with mixed ferroelectric–antiferroelectric order, directly integrated onto Si transistors, and scaled down to approximately 20 ångströms, the same gate oxide thickness required for high-performance transistors. The overall equivalent oxide thickness in metal–oxide–semiconductor capacitors is equivalent to an effective SiO 2 thickness of approximately 6.5 ångströms. Such a low effective oxide thickness and the resulting large capacitance cannot be achieved in conventional HfO 2 -based high-dielectric-constant gate stacks without scavenging the interfacial SiO 2 , which has adverse effects on the electron transport and gate leakage current. Accordingly, our gate stacks, which do not require such scavenging, provide substantially lower leakage current and no mobility degradation. This work demonstrates that ultrathin ferroic HfO 2 –ZrO 2 multilayers, stabilized with competing ferroelectric–antiferroelectric order in the two-nanometre-thickness regime, provide a path towards advanced gate oxide stacks in electronic devices beyond conventional HfO 2 -based high-dielectric-constant materials.

36 MATERIALS SCIENCE↗

DOE OSTI2024

SAM Code Enhancements for Modeling of Liquid Metal-Cooled Fast Reactor Concepts

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. This report summarizes two major efforts in addressing the modeling gaps in SAM for liquid-metal-cooled fast reactor (LMFR) applications, i.e. thermal mixing and stratification phenomena in large pools and corrosion-oxidation of components in flowing lead. A new one-dimensional model for thermal mixing and stratification effects in large pools and enclosures is developed and implemented. Thermal mixing and stratification occur when fluid enters a pool at a temperature different than the bulk fluid itself, a scenario often encountered during transients in pool-type systems. These phenomena are critical for the safety of reactors, impacting phenomena like natural circulation, which is essential for passive cooling. The improved model in SAM addresses limitations of state-of-the-art approaches by combining one-dimensional (1D) channels, representing the coolant jet flow, with lumped-parameter zero-dimensional (0D) pools, representing the rest of coolant in the tank. Energy exchange between the 1D jet and the 0D pools is based on heat transfer correlations calibrated against 3D simulations. It is verified that this model can handle various flow configurations, including hot jets in colder pools, cold jets in hotter pools, and the presence of features like ceilings, free surfaces, and obstacles. Additionally, validation against experimental data demonstrates the ability of the model to capture mixing and stratification effects in a wide range of conditions. The flexibility and improved accuracy of the new model make it a valuable tool for reactor safety analysis, allowing for the simulation of different geometries encountered in advanced reactors. A system-level corrosion modeling capability is developed and implemented in SAM to support Lead Fast Reactor (LFR) development. Although the initial focus of this capability will be on LFR application, this can later be expanded to include other liquid metals such as Lead-Bismuth Eutectic (LBE) and PbLi. This report summarizes the common corrosion mitigation strategies and outlines the progress on implementing and validating a corrosion-oxidation model in SAM. Verification and validation of the corrosion-oxidation portion of the model was performed using analytical solution and measured data from samples tested in the non-isothermal pumped lead loop at IPPE Obninsk. The iron transport and corrosion/precipitation portion of the model was assessed using an analytical model and measured corrosion depths from a natural convection lead loop experiment performed at CEA. It is demonstrated that the model implemented in SAM performed well in these assessments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Convective and Turbulent Motions in Nonprecipitating Cu. Part II: LES Simulated Cloud Represented by a Starting Plume

The dynamic structure of a small trade wind cumulus (Cu) is analyzed using a novel approach. Cu developing in a shear-free environment is simulated by 10-m-resolution LES model with spectral bin microphysics. The aim is to clarify the dynamical nature of cloud updraft zone (CUZ) including entrainment and mixing in growing Cu. The validity of concept stating that a cloud at developing state can be represented by a parcel or a jet is tested. To investigate dynamical entrainment in CUZ performed by motions with scales larger than the turbulence scales, the modeled fields of air velocity were filtered by wavelet filter that separated convective motions from turbulent ones. Two types of objects in developing cloud were investigated: small volume ascending at maximal velocity (point parcel) and CUZ. It was found that the point parcel representing the upper part of cloud core is adiabatic. The motion of the air in this parcel ascending from cloud base determines cloud-top height. The top-hat (i.e., averaged) values of updraft velocity and adiabatic fraction in CUZ are substantially lower than those in the point parcel. Evaluation of the terms in the dynamical equation typically used in 1D cloud parcel models show that this equation can be applied for calculation of vertical velocities at the developing stage of small Cu, at least up to the heights of the inversion layer. Dynamically, the CUZ of developing cloud resembles the starting plume with the tail of nonstationary jet. Both the top-hat vertical velocity and buoyancy acceleration linearly increase with the height, at least up to the inversion layer. An important finding is that lateral entrainment of convective (nonturbulent) nature has a little effect on the top-hat CUZ velocity and cannot explain the vertical changes of conservative variables qt and θl. In contrast, entrained air lifting inside CUZ substantially decreases top-hat liquid water content and its adiabatic fraction. Possible reasons of these effects are discussed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Dynamics and observational signatures of core-collapse supernovae with central engines: hydrodynamics simulations with Monte Carlo post-processing

A long-lived central engine embedded in expanding supernova ejecta can alter the dynamics and observational signatures of the event, producing an unusually luminous, energetic, and/or rapidly evolving transient. We use 2D hydrodynamics simulations to study the effect of a central energy source, varying the amount, rate, and isotropy of the energy deposition. We post-process the results with a time-dependent Monte Carlo radiation transport code to extract observational signatures. The engine excavates a bubble at the centre of the ejecta, which becomes Rayleigh–Taylor unstable. Sufficiently powerful engines are able to break through the edge of the bubble and accelerate, shred, and compositionally mix the entire ejecta. The breakout of the engine-driven wind occurs at distinct rupture points, and the outflowing high-velocity gas may eventually give rise to radio emission. The dynamical impact of the engine leads to faster rising optical light curves, with photon escape facilitated by the faster expansion of the ejecta and the opening of low-density channels. For models with strong engines, the spectra are initially hot and featureless, but later evolve to resemble those of broad-line Ic supernovae. Under certain conditions, line emission from ionized, low-velocity material near the centre of the ejecta may be able to escape and produce narrow emission similar to that seen in interacting supernovae. We discuss how variability in the engine energy reservoir and injection rate could give rise to a heterogeneous set of events spanning multiple observational classes, including the fast blue optical transients, broad-line Ic supernovae, and superluminous supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗