Search NASA⌕ Search

SEARCH · Search NASA

Results for “two-phase models”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Self-equilibrated backstresses induce compensation between hardening and softening: Micromechanical and microstructural features

For crystal plasticity formulation neglecting the strain gradients at the dislocation scale, backstresses are necessary to satisfy the mesoscopic stress equilibrium and strain compatibility requirements. Backstresses in a material point and those between neighboring points are treated separately. Physically backstresses can dynamically induce either hardening or softening effects to adjust the dislocation hardening of slip systems, whereas the sum of the stresses resulting in backstresses in the representative volume element (RVE) is zero. Backstress hardening suppresses certain slip activities, while backstress softening triggers excess dislocations. A two-phase brass is studied. Backstresses make the hard phase deform smaller and the soft phase deform more than those without employing backstresses in a model. In the former, more homogeneous lattice strain and higher yield stress are obtained. Backstress has limited effects on the evolution of the critical resolved shear stress and crystallographic texture. Furthermore, it is a major cause of the banded structures featuring sharp crystallographic contrasts which are not parallel to the slip planes. Tests indicate that backstress induces a serrated behavior at the early stages of deformation and induces the geometrically necessary bands (GNBs), which are parallel to the macroscopic axes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Impact of particle loading and phase coupling on gas–solid flow dynamics: A case study of a two-phase, gas–solid flow in an annular pipe

The present study is devoted to a two-phase, gas–solid flow in an annular pipe (hollow cylinder) at an elevated pressure of 15 bars and moderate Reynolds number of circa 6000. The influence of the particle loading, the interaction between the phases, and turbulence dispersion on the flow dynamics is systematically studied by means of computational fluid dynamics simulations, employing the Ansys FLUENT commercial package. The cases with a particle volumetric fraction of 1.2% are referred to as “high particle loading,” and those with 0.13% are denoted as “low particle loading.” The following cases are investigated: (1) pure gas flow; (2) low particle loading two-phase flow with one-way coupling and with turbulence dispersion; (3) low particle loading two-phase flow with two-way coupling but without turbulence dispersion; (4) low particle loading two-phase flow with two-way coupling and with turbulence dispersion; (5) high particle loading two-phase flow with one-way coupling and with turbulence dispersion; (6) high particle loading two-phase flow with two-way coupling but without turbulence dispersion; and (7) high particle loading two-phase flow with two-way coupling and with turbulence dispersion. Additionally, the boundary layer is found to grow without fluctuations of the turbulent kinetic energy (TKE) for cases 1, 2, and 5. For case 4, the TKE fluctuations have been identified, although they appear to be less substantial than those in cases 6 and 7. The authors attribute the semi-chaotic nature of the TKE fluctuations to the particle loading and two-way coupling. In addition, the onset and development of the flow instability have been observed at a random axial distance in cases 4, 6, and 7. Such instability is also attributed to the two-way coupling with turbulence dispersion in the flow. It is concluded that the particle loading, one-way, or two-way coupling between the phases, and the turbulence dispersion models significantly influence the development of the flow dynamics with the same inlet and boundary conditions. Consequently, it is not a trivial question, which result a user should trust. The present computational results inspire to perform verification as well as experimental validation of the simulations, so the simulation results can subsequently be used with confidence for design analysis.

42 ENGINEERING↗
DOE OSTI2022

Constraining the release of Sn to the ambient melting point following shock loading using time-resolved x-ray diffraction

The formation of liquid following release from a shocked state governs the transition from spall to cavitation and the formation of ejecta in metals. In order to build physics-based models of these processes, it is necessary to critically evaluate the relative importance of kinetics and entropy generation during the release along with the accuracy of multiphase equations of state. Tin (Sn) has served as a testbed for a variety of experiments examining strength and ejecta due to its accessible melt boundary and solid–solid phase transitions. This work presents experiments examining the phase evolution of high purity Sn following the shock and release to ambient stress near the melting point. Sn is found to release to states between its ambient solidus and liquidus from approximately 19 to 33 GPa under uniaxial loading, with the two-phase region being characterized by a reduction in the intensity of the (220), (211) β -Sn doublet. Jetting experiments performed at 27–28 GPa exhibit comparable diffraction patterns with what is observed following the uniaxial release. The solid fractions of β-Sn in the ambient mixed phase region are found to decrease linearly with increasing shock stress as increasing liquid Sn is formed. Here, the results provide much needed information for interpreting measurements of dynamic strength at a high strain rate and experiments examining cavitation and shallow bubble collapse in Sn.

, spall↗
DOE OSTI2021

Pore-scale simulation of drying in porous media using a hybrid lattice Boltzmann: pore network model

In this work, a hybrid method coupling a pseudo-potential lattice Boltzmann model (LBM) and a pore network model (PNM) to simulate drying in porous media is proposed. Based on the watershed method, the porous medium is firstly decomposed into pore regions. According to the liquid–vapor phase distribution at a given time, the pore regions are further divided into four pore types, namely two-phase pores where a liquid–vapor interface exists, buffer pores next to the two-phase pores, single-liquid and single-vapor phase pores. The pseudo-potential LBM is used in the two-phase and buffer pores to simulate liquid drying and track the movement of the interfaces, while the single-phase PNM simulations are conducted in the buffer and single-phase pores to simulate vapor or liquid flow. LBM and PNM are coupled in the buffer pores through exchange of boundary information. The hybrid method is applied to simulate liquid drying in a porous medium. The whole-domain LBM simulation is considered as the reference solution to validate the hybrid method. Liquid saturation variation during the drying process and detailed phase and pressure distributions obtained by the two methods match quite well, demonstrating the accuracy of the hybrid method. For the specific case studied, the hybrid method saves more than 60% computational time compared to the whole-domain LBM simulation. In addition, the speedup of the hybrid method becomes more significant for a larger computational domain. In summary, the hybrid method developed in this work combines the accuracy of LBM and the efficiency of PNM to simulate drying in porous media at pore scale and can lead to significant reduction of computation time, thus allowing the pore-scale consideration of drying in larger porous systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The electric field dependence of single electron emission in the PIXeY two-phase xenon detector

Dual phase xenon detectors are widely used in experimental searches for galactic dark matter particles. The origin of single electron backgrounds following prompt scintillation and proportional scintillation signals in these detectors is not fully understood, although there has been progress in recent years. Here, in this paper, we describe single electron backgrounds in 83m Kr calibration events and their correlation with drift and extraction fields, using the Particle Identification in Xenon at Yale (PIXeY) dual-phase xenon time projection chamber. The single electron background induced by the Fowler-Nordheim (FN) effect is measured, and its electric field dependence is quantified. The photoionization of grids and impurities by prompt scintillation and proportional scintillation also contributes to the single electron background.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

SILCC – VII. Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities

ABSTRACT We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σgas = 10, 30, 50, and 100 $\mathrm{M_\odot \, pc^{-2}}$. The silcc project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density ΣSFR. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $\sigma _\mathrm{z} \propto \Sigma _\mathrm{SFR}^a$, with a ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼2.2 higher than that of the CNM, in agreement with local observations. For $\Sigma _\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\, \mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$ the WIM motions become supersonic.

Astronomy & Astrophysics↗
DOE OSTI2022

Early-type dwarf galaxies in the local universe. Evidence of ex situ growth

We report the discovery of a rare early-type dwarf galaxy (dE), SDSS J125651.47+163024.2 (hereafter, dE1256), possessing a tidal feature that was likely built up by accretion of an even smaller dwarf galaxy. dE1256 is located in a nearly isolated environment, at the outskirt of the Virgo cluster. A detailed morphological examination reveals that the accreted stellar population is mainly deposited in the outer part of dE1256, where the tidal tail is most prominent. The inner part of dE1256 is perfectly modelled with a simple Sérsic function of index n = 0.63 and half-light radius R h = 0.6 kpc, but in contrast, the entire galaxy has a size of R h = 1.2 kpc. The mass ratio between the host and the putative accreted dwarf galaxy is calculated to be 5:1, assuming that the observed two components, inner Sérsic and outer tidal tail residual, represent the host’s and accreted galaxy’s stellar populations, respectively. We suggest that while the accretion contributes only 20 per cent of the overall stellar population, the size of dE1256 grew by a factor of two via the accretion event. Our results provide, for the first time, strong observational evidence that a dE is undergoing a two-phase growth, a common phenomenon for massive galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Three and Two Phase Rotating Field Inductive Couplers for Wireless Power Transfer with One Phase per Layer Windings

Multiphase inductive wireless charging coils have been proposed recently to improve coupler surface power density, reduce component stress and size, and provide near-constant power delivery to charge mobile electric systems. Several aspects for the fundamental characterization of multiphase coils are explored up to six phases including approximate mutual inductance with size and turn variation, induced voltage, and output power estimation. The relative component stress and size of passive components for resonant operation are compared between the multiphase variants. A combination of an experimentally validated 3D electromagnetic finite element analysis (FEA) and power electronic co-simulations are used to validate the estimated quantities approximated with a mixture of analytical equations and an artificial neural network model for mutual inductance. A novel three-phase transmitter, two-phase receiver coil pair is also proposed for electric vehicle charging to reduce the number of connections and compensation complexity on the vehicle-side with improved power output compared to a two-phase configuration.

Lewis, Donovin D. [University of Kentucky]↗
DOE OSTI2021

Flow Assisted Evaporative Cooling for Electric Motor

This paper examines a novel concept of flow assisted latent heat driven two-phase evaporative cooling (EC) confined in-between slot liner and active-winding of electric motor. Wicking micro-structure enhanced PDMS liner axially sucks coolant in the form of thin film between the PDMS liner and active-winding and eventually enables thin film evaporation on the outer surface of the active-winding. Therefore, EC based thermal management eliminates contact resistance between the winding and slot-liner and, enhances the heat extraction from the winding without compromising the electro-magnetic performance. Two-way coupled electro-magnetic (EM) – computational fluid dynamics/heat transfer (CFD/HT) and EM - lumped parameter thermal network (LPTN) models have been developed to assess the electro-thermal performance of the EC under steady and transient conditions. Taking a case study of a 125 kW jacket cooled BMW i3 motor and dielectric coolant FC-84, EC is shown to be capable of handling a maximum steady state rms current density of 26 A/mm2 at a evaporative heat transfer coefficient of 5,000 W/m2.K, which is about 78.7% higher compared to the traditional jacket cooling (JC). In case of EC, a maximum steady state and transient rms current density of 30 A/mm2 (106.2% higher compared to the JC), and 40.8 A/mm2 have been realized by using high thermal conductivity epoxy (1.9 W/m.K) impregnation material. Thermal performance of the EC is also assessed and compared with JC over a dynamic drive cycle. Lastly, a motorette testing has been performed to demonstrate the applicability of the proposed EC method and to validate the developed modeling framework.

42 ENGINEERING↗
DOE OSTI2025

Switchgrass ( Panicum Virgatum ) and Miscanthus ( Miscanthus × Giganteus ) Long-Term Yield Patterns Reveal Consistent Productivity Declines

Perennial grasses like switchgrass ( Panicum virgatum ) and miscanthus ( Miscanthus × giganteus ) are expected to supply a substantial amount of the United States bioeconomy's feedstock demand. However, uncertainties around their long-term yields challenge the viability of their potential and limit their wider adoption. To resolve their long-term yield patterns, we analyzed over 200 plantings of switchgrass and miscanthus across Michigan and Wisconsin, USA, measured over 5–15 years. We found a consistent two-phase long-term yield dynamic; during a yield-building phase , peak yields occurred within 4–5 years after planting, followed by a yield-decline phase in which switchgrass and miscanthus lost 30%–47% and 14%–40% of peak yields, respectively. Among the potential drivers of this dynamic and the yield decline, we found that weather conditions had little impact, as the variation across years was not large enough to drive the observed yield differences. Added nitrogen increased peak yields by 10%–20% and attenuated the yield decline by 20%–50%. However, since fertilized stands still showed a yield decline, other factors became limiting as stands aged. This conserved long-term yield dynamic has direct implications on management. A farm-to-gate economic analysis suggests replanting switchgrass and miscanthus 5 and 9 years following their peak yields maximizes profit over a 30-year time horizon. Results call for further management and breeding strategies to mitigate the yield-decline phase, and for reparameterization of global bioenergy models with carbon capture and storage, which may overestimate yields and the economic and environmental benefits of crops grown for bioenergy feedstocks.

bioenergy↗
DOE OSTI2021

Modeling Electrokinetics of Oxygen Electrodes in Solid Oxide Electrolyzer Cells

A microscale model is presented in this study to simulate electrode kinetics of the oxygen electrode in a solid oxide electrolyzer cell (SOEC). Two mixed ionic/electronic conducting structures are examined for the oxygen producing electrode in this work: single layer porous lanthanum strontium cobalt ferrite (LSCF), and bilayer LSCF/SCT (strontium cobalt tantalum oxide) structures. A yttrium-stabilized zirconia (YSZ) electrolyte separates the hydrogen and oxygen electrodes, as well as a gadolinium doped-ceria (GDC) buffer layer on the oxygen electrode side. Electrochemical reactions occurring at the two-phase boundaries (2PBs) and three-phase boundaries (3PBs) of single-layer LSCF and bilayer LSCF/SCT oxygen electrodes are modeled under various SOEC voltages with lattice oxygen stoichiometry as the key output. The results reveal that there exists a competition in electrode kinetics between 2PBs and 3PBs, but 3PBs are the primary reactive sites for single-layer LSCF oxygen electrode under high voltages. These locations experience the greatest oxygen stoichiometry variations and are therefore the most likely locations for dimensional changes. By applying an active SCT layer over LSCF, the 2PBs become activated to compete with the 3PBs, thus alleviating oxygen stoichiometry variations and reducing the likelihood of dimensional change. This strategy could reduce lattice structural expansion, proving to be valuable for electrode-electrolyte delamination prevention and will be the focus of future work.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Understanding the Effect of Electrochemical Properties and Microstructure on the Microgalvanic Corrosion of Mg Alloys via Phase-Field Simulations

In this work, we apply a phase-field model to gain insights into the corrosion behavior of Mg alloys. Specifically, we study the effect of electrochemical properties of the materials and their environment and the spatial distribution of second phases in a magnesium alloy. To this end, we perform sensitivity analyses in which we separately vary the exchange current density and corrosion potential of the second phase, as well as the electrolyte conductivity. We describe the effects of these parameters on the corrosion current and the resulting corrosion-front morphology. We also confirm that the ratio between the Wagner length, calculated for the anodic phase, and the initial width of the anodic domain is a good predictor of the morphological characteristic of the corroding surface. In addition, we examine the effect of the microstructure on the corrosion behavior in both 2D and 3D systems by varying the shape, size, and distribution of the second-phase regions. Our findings provide insights into the influence of the two-phase microstructure on the corrosion behavior in a magnesium alloy.

Electrochemistry↗
DOE OSTI2025

Combining Distributed Electrochemical Measurements with a Semi-Empirical Model to Identify Local Variations in Overpotentials for PEM Water Electrolyzers

Abstract A segmented cell and PCB approach was used to demonstrate an advanced technique for measuring and understanding the electrochemical behavior of PEM water electrolyzers. By employing a semi-empirical model, discrete modes of overpotential and their spatial distributions were identified in an operating cell. Distributed area-specific resistance (DASR measurements indicate that local dehydration of the polymer electrolyte membrane is the primary driver of increased polarization, behaving similarly to mass transport-limited currents at low flow rates due to ineffective gas-phase bubble removal. To improve quantification, the overall ohmic overpotential is divided into background (η Ohm,b ) and mass transport-coupled (η Ohm,mt ) components. Two types of anode diffusion media, a baseline PTL and a thin foil liquid/gas diffusion layer (LGDL), were examined. The results show significant variations in η Ohm,b and highlight that ineffective bubble removal leads to pronounced local dehydration. The effect of cell temperature on various overpotentials was also examined. Negligible differences in diffusion overpotential were observed when comparing PEWE performance using the PTL and LGDL, despite major differences in two-phase fluid transport within each diffusion media. This underscores the critical role of the catalyst layer's micro-porous structure in diffusion processes. The study highlights the complex nature and intricacies of mass transport limitations within PEWEs.

Roy, Anirban (ORCID:0000000163803386)↗
DOE OSTI2021

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗
DOE OSTI2021

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Report on Year-4 of Water NSTF Matrix Testing: Facility Maintenance and Accident Testing

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report serves as a summary of maintenance and experimental activities during the program’s fourth year of water-based operation. A planned six-month maintenance period began in August 2021, during which major inspections, repairs, cleaning, and installation of new instrumentation and data acquisition hardware were conducted. Most significantly, two heaters that faulted during Year-3 were repaired, allowing the facility to resume use of the full heated section area and full range of available electric power. The remainder of the year consisted of eight months of test operations, during which the facility logged 211 hours of active heating across one bake-out (following the maintenance period) and seven matrix tests; five classifieds as Accepted per NQA-1, one as Trending, and one as Failed. Testing began by performing two repeat cases to confirm expected facility response and behavior during both single- and two-phase flow conditions, ensuring no changes were introduced during the maintenance period that might have altered the thermal-hydraulic characteristics of the facility. In continuation of the power parametric series initiated in previous years, a high-power test case was then performed examining heat removal performance at a decay heat load equivalent to 2.4 MWt, full-scale, a level exceeding maximum design targets. Additional testing then introduced various blockage scenarios along the network piping, examining the effects of partial and complete blockages of the flow paths on the system behavior and heat removal performance. A study of static boiling tests directed at understanding the geysering two-phase instability was also conducted. The loop was filled only to the bottom of the tank outlet, creating an open loop configuration that prevents any natural circulation flow from occurring, and the heaters were powered on until the facility reached saturation conditions. Following, a series of quasi-steady-state conditions were introduced by adjusting the inventory level in the adiabatic chimney piping at decreasingly lower elevations above the heated region. A strong correlation of geysering characteristics and loop level was observed, with flow and temperature excursions decreasing in intensity, but increasing in frequency, as the fill were reduced to lower elevations along the chimney piping. Once the level fell very low in the chimney, at points near the top of the heated section, the system reached a stable state of continuous boiling without any occurrence of geysering eruptions. A final significant testing accomplishment this year was successful completion of an “accident scenario” test, whose operating conditions were based on a prototypic decay heat curve provided by Framatome and scaled for the NSTF. This test began by establishing steady-state, single-phase “normal operation” conditions, before simulating an accident trip where the availability of active cooling systems was lost. Loop temperatures gradually increased until reaching saturation and subsequent two-phase boiling flow. Over the course of an extended operational period along the defined decay heat curve, steam boil-off caused gradual but continued depletion of liquid inventory until reaching a critically low level causing flow stagnation and cessation of natural circulation heat removal. At this point, after nearly 72 hours of continuous operation, a cold refill was performed to replenish the system inventory and allow the facility to re-establish closed loop natural circulation flow and return to a safe operational state.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Deflagration to Detonation Transition Update: XDDT Code Modularization

A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

A Modeling and Neutron Diffraction Study of the High Temperature Properties of Sub-Stoichiometric Yttrium Hydride for Novel Moderator Applications

Low-enriched-uranium (LEU) reactor systems utilize moderators to improve neutron economy. Solid yttrium hydride is one of the primary moderator candidates for high-temperature (>700 °C) nuclear reactor applications. This is due to its ability to retain hydrogen at elevated temperatures compared to other metal hydrides. For reactor modeling purposes, both neutronic and thermos-mechanical modeling, several high-temperature properties for sub-stoichiometric yttrium hydride (YH2−x) are needed. In this paper, we present an atomistics and a neutron diffraction study of the high-temperature properties of Y and YH2−x. Specifically, we focus on the thermal lattice expansion effects in yttrium metal and yttrium hydride, which also govern bulk thermal expansion. Previously reported physical and mechanical properties for sub-stoichiometric yttrium hydride at ambient conditions are expanded using lattice dynamics to take into account high-temperature effects. Accordingly, an array of newly generated properties is presented that enables high-fidelity neutronics, and thermomechanical modeling. These properties include various elastic moduli, thermal expansion parameters for yttrium and yttrium hydride, and single-phase (YH2−x) and two-phase (Y + YH2−x) density as a function of stoichiometry and density.

36 MATERIALS SCIENCE↗