Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial velocity”

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.

70 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Interfacial Nanostructure and Hydrogen Bond Networks of Choline Chloride and Glycerol Mixtures Probed with X-ray and Vibrational Spectroscopies

The molecular distribution at the liquid-vapor interface and evolution of the hydrogen bond interactions in mixtures of glycerol and choline chloride are investigated using X-ray photoelectron spectroscopy (XPS) and Raman spectroscopy. Nanoscale depth profiles of supersaturated deep eutectic solvent (DES) mixtures up to ~2 nm measured by ambient-pressure XPS show the enhancement of choline cation (Ch + ) concentration by a factor of 2 at the liquid-vapor interface compared to the bulk. In addition, Raman spectral analysis of a wide range of DES mixtures reveals the conversion of gauche-conformer Ch + into the anti-conformer in relatively lower ChCl concentrations. Finally, the depletion of Ch + from the interface (probing depth = 0.4 nm) is demonstrated by aerosol-based velocity map imaging XPS measurements of glyceline and water mixtures. The nanostructure of liquid-vapor interfaces and structural rearrangement by hydration can provide critical insight into the molecular origin of the deep eutectic behavior and gas-capturing application of DESs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Gold-Induced Chemical Perturbations in CdTe-Based Photovoltaic Cells

Back-contacting p-type CdTe has been identified as one of the major areas of loss in CdTe photovoltaic (PV) power conversion efficiency (PCE). In research settings, Au is a common contact material due to its ease of use and decent performance. This work provides a detailed investigation into using gold for back contacting As-doped, CdCl 2 -treated, polycrystalline CdTe that has been exposed to air after absorber processing, another routine practice. First, X-ray photoemission spectroscopy (XPS) is used to determine the native oxide to be 1.6 nm of CdTeO 3 by using a combination of angle-resolved XPS and the cadmium-modified Auger parameter. During gold metallization of CdTe, oxygen and oxidized tellurium are eliminated from the thin CdTeO 3 native oxide. The fates of the released oxygen and possibly cadmium and tellurium are not known, but these reaction byproducts can enter the absorber bulk or grain boundaries, stay at the interface, or dissolve in the Au. Interfacial hole barriers between CdTe and Au are measured for samples with and without the native oxide present prior to metallization. Results show that the thin CdTeO 3 alleviates the downward band bending by 40 meV from 470 to 430 meV even though it is consumed during interface formation. The implications of these chemical reactions on the device are assessed through photoluminescence (PL) spectroscopy, which shows losses in internal open circuit voltage (iVoc) from 820 to 795 meV, carrier lifetime from 123 to 45 ns, and PL quantum yield from 2.9 × 10 −5 to 1.2 × 10 −5 . Modeling timeresolved PL lifetimes demonstrates the back surface recombination velocity due to metallization reduces minority carrier lifetimes. These results identify the native oxide and show that it plays an important role in mediating downward band bending along with how the back interface reaction can negatively impact device-scale parameters and reduce PV PCE.

14 SOLAR ENERGY↗
DOE OSTI2026

ThermoPI—an Online Tool to Calculate Heat Transfer Through Foam Insulation

Thermal insulation materials with ultra-low effective thermal conductivity are crucial for a multitude of applications. Over the years, a significant amount of experimental work has been dedicated to creating new insulation materials with low thermal conductivity. Similarly, substantial efforts have been made to enhance the theoretical understanding of thermal transport mechanisms in thermal insulation materials and to push the boundaries of lower thermal conductivity. However, ample room remains for enhancing the thermal resistivity of closed-cell foam insulations. To aid in the development of ultra-low effective thermal conductivity foam insulations, Oak Ridge National Laboratory has introduced a unique online tool, ThermoPI. This tool calculates not only the effective thermal conductivity but also the thermal conductivity components for porous materials, including gases, solids, and radiative, based on the materials’ structural information (e.g., porosity, pore size, gas species, pressure, solid species, pore geometry, and temperature).This tool collected and improved the existing theoretical models for thermal insulation materials’ gas, solid, and radiative thermal conductivity. These improved models have been validated with experimental data. Effective medium theory and sound velocity softening effects are considered for solid thermal conductivity. The built-in solid materials include polystyrene, polyurethane, polyethylene, and silica. Users are allowed to input new solid materials that are not built inside the tool. For gas thermal conductivity, the subcontinuum Knudsen effect is considered. In addition to several built-in gases, users can input new gases. For solid and radiative thermal conductivity calculations, several models are available to select, and the tool can determine the best model to choose based on the materials that users input. Users are also allowed to change selections manually. In addition, the tool can also calculate the effect of interfacial resistance on the overall thermal conductivity of layered materials. This work will elaborate on the tool and discuss how it can guide the development of new insulation products.

Shrestha, Som [ORNL] (ORCID:0000000183993797)↗
DOE OSTI2020

Probing the meta-stability of oxide core/shell nanoparticle systems at atomic resolution

Hybrid nanoparticles allow exploiting the interplay of confinement, proximity between different materials and interfacial effects. However, to harness their properties an in-depth understanding of their (meta)stability and interfacial characteristics is crucial. This is especially the case of nanosystems based on functional oxides working under reducing conditions, which may severely impact their properties. Here, the in-situ electron-induced selective reduction of Mn 3 O 4 to MnO is studied in magnetic Fe 3 O 4 /Mn 3 O 4 and Mn 3 O 4 /Fe 3 O 4 core/shell nanoparticles by means of high-resolution scanning transmission electron microscopy combined with electron energy-loss spectroscopy. Such in-situ transformation allows mimicking the actual processes in operando environments. A multi-stage image analysis using geometric phase analysis combined with particle image velocity enables direct monitoring of the relationship between structure, chemical composition and strain relaxation during the Mn 3 O 4 reduction. In the case of Fe 3 O 4 /Mn 3 O 4 core/shell the transformation occurs smoothly without the formation of defects. However, for the inverse Mn 3 O 4 /Fe 3 O 4 core/shell configuration the electron beam-induced transformation occurs in different stages that include redox reactions and void formation followed by strain field relaxation via formation of defects. This study highlights the relevance of understanding the local dynamics responsible for changes in the particle composition in order to control stability and, ultimately, macroscopic functionality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Data for "Gold-Induced Chemical Perturbations in CdTe-Based Photovoltaic Cells"

Back contacting p-type CdTe has been identified as one of the major areas of loss in CdTe photovoltaic (PV) power conversion efficiency (PCE). In research settings, Au is a common contact material due to its ease of use and decent performance. This work provides a detailed investigation into using gold for back contacting As-doped, CdCl2 treated, polycrystalline CdTe that has been exposed to air after absorber processing, another routine practice. First, X-ray photoemission spectroscopy (XPS) is used to determine the native oxide to be 1.6 nm of CdTeO3 using a combination of angle-resolved XPS and the cadmium modified Auger parameter. During gold metallization of CdTe, oxygen and oxidized tellurium are eliminated from the thin CdTeO3 native oxide. The fate of the released oxygen and possibly cadmium and tellurium are not known, but these reaction byproducts can enter the absorber bulk or grain boundaries, stay at the interface, or dissolve in the Au. Interfacial hole barriers between CdTe and Au are measured for samples with and without the native oxide present prior to metallization. Results show that the thin CdTeO3 alleviates the downward band bending by 40 meV from 470 meV to 430 meV even though it is consumed during interface formation. The implications of these chemical reactions on the device are assessed through photoluminescence (PL) spectroscopy which shows losses in internal open circuit voltage (iVoc) from 820 meV to 795 meV, carrier lifetime from 123 ns to 45 ns, and PL quantum yield from 2.9x10-5 to 1.2x10-5. Modeling time-resolved PL lifetimes demonstrates the back surface recombination velocity due to metallization reduces minority carrier lifetimes. These results identify the native oxide and show that it plays an important role in mediating downward band bending along with how the back interface reaction can negatively impact device-scale parameters and reduce PV PCE.

14 SOLAR ENERGY↗
DOE OSTI2023

Characterization of local parameters of inclined upward air–water bubbly two-phase flows in a round pipe

The present work seeks to develop a detailed experimental database and characterize the local effects of inclination on horizontal and upward bubbly air–water two-phase flows. Experimental conditions are selected to investigate bubbly flows at a variety of orientations while reducing the influence of flow regime transition. Detailed time-averaged data is collected using state-of-the-art four-sensor conductivity probes for measurements across the flow area for five flow conditions at five inclination orientations in a 25.4 mm inner diameter pipe. Using the newly established database, effects of inclination on bubbly flows are analyzed, as are changes in liquid and gas flow rate at different orientations. The observations include (1) bubble migration and agglomeration at the top pipe wall as orientation changes from vertical to horizontal; (2) reducing the angle from vertical reduces the local gas velocity, and correspondingly reduces the global slip from greater than unity at 90° to less than unity for 60° and lower; (3) changes in the bubble size distribution are primarily at the interface between the bubble cluster and the liquid region of the pipe, not coalescence within the bubble cluster. In conclusion, physical mechanisms for these observations are proposed, and the present trends at different orientations support previous experiments in horizontal two-phase flows.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Colloquium : Sliding and pinning in structurally lubric 2D material interfaces

A plethora of two-dimensional (2D) materials have been introduced in physics and engineering in the past two decades. Their robust, membranelike sheets permit (mostly require) deposition, giving rise to solid-solid dry interfaces whose mobility, pinning, and general tribological properties under shear stress are currently being understood and controlled, both experimentally and theoretically. Here, in this Colloquium simulated case studies of twisted graphene systems are used as a prototype workhorse tool to demonstrate and discuss the general picture of 2D material interface sliding. First highlighted is the crucial mechanical difference, often overlooked, between small and large incommensurabilities, which corresponds to, for example, small and large twist angles in graphene interfaces. In both cases, focusing on flat, structurally lubric or “superlubric” geometries, the generally separate scalings with the area of static friction in pinned states and of kinetic friction during sliding are elucidated and reviewed, tangled as they are with the effects of velocity, temperature, load, and defects. The roles of island boundaries and elasticity are also discussed, and compared when possible to results in the literature for systems other than graphene. It is proposed that the resulting picture of pinning and sliding should be applicable to interfaces in generic 2D materials that are of importance for the physics and technology of existing and future bilayer and multilayer systems.

2-dimensional systems↗
DOE OSTI2025

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Influence of sizing concentration on strength, stiffness, and porosity in textile grade carbon fiber (TCF)-Epoxy composites: Revealing inverse trends

The effect of fiber sizing (i.e., surface treatment) concentration (0 %, 1.36 %, 1.52 %, 1.94 %, and 2.13 %) on the mechanical properties (tensile, flexural, interlaminar shear strength (ILSS), and low velocity impact) of textile grade carbon fiber (TCF)-epoxy composite is examined. An inverse relationship between the strength and stiffness of the composite is observed with increased sizing concentration. The root mean square (RMS) roughness of the fiber surface increased from 17.8 nm (unsized) to 22.7 nm with 2.13 % sizing concentration. It was found that the tensile strength increased by 131 % from 221.4 ± 18.5 MPa (unsized) to 510.8 ± 28.05 MPa (for 1.36 % sizing) and further by 155 %–563.7 ± 14.95 MPa at 2.13 % sizing. On the contrary, the stiffness is initially increased by 126 % from 33.52 ± 7.80 GPa (unsized) to 75.9 ± 3.21 GPa (for 1.36 % sizing) but reduced with further increase in the sizing concentration. A single fiber pull-out test is simulated using the finite element method to validate the reverse trend in strength and stiffness. The varying sizing concentration is simulated by introducing an interface of varying thickness between fibers and matrix. Simulation results confirm that a thicker interface, corresponding to a higher sizing concentration, decreases interfacial shear stress, enhancing material strength while decreasing stiffness. The reverse trend in strength and stiffness with the sizing concentration aligns with experimental observations. In conclusion, the present study emphasizes the importance of sizing concentration for mechanical properties and provide a design criterion for customized high-strength and high-stiffness applications.

Porosity↗
DOE OSTI2022

Computing the viscous effect in early-time drop impact dynamics

The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. ( Phys. Rev. Lett. , vol. 112, no. 13, 2014 b , p. 134501) show that the liquid–air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier–Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square-root scaling relationship. We also predict the effect of interfacial tension at the liquid–gas interface on the drop impact, showing that it mediates the lift-off behaviour.

42 ENGINEERING↗
DOE OSTI2020

Wettability and capillary effects: Dynamics of pinch-off in unconstricted straight capillary tubes

We study the interfacial evolution of immiscible two-phase flow within a capillary tube in the partial wetting regime using direct numerical simulation. We investigate the flow patterns resulting from the displacement of a more viscous fluid by a less viscous one under a wide range of wettability conditions. We find that beyond a wettability dependent critical capillary number, a uniform displacement by a less viscous fluid can transition into a growing finger that eventually breaks up into discrete blobs by a series of pinch-off events for both wetting and non-wetting contact angles. This study validates previous experimental observations of pinch-off for wetting contact angles and extends those to non-wetting contact angles. We find that the blob length increases with the capillary number. We observe that the time between consecutive pinch-off events decreases with the capillary number and is greater for more wetting conditions in the displaced phase. Here, we further show that the blob separation distance as a function of the difference between the inlet velocity and the contact line speed collapses into two monotonically decreasing curves for wetting and non-wetting contact angles. For the phase separation in the form of pinch-off, this work provides a quantitative study of the emerging length and time scales and their dependence on the wettability conditions, capillary effects, and viscous forces.

42 ENGINEERING↗
DOE OSTI2021

High-speed impact of micron-sized diesel drop trains—Splashing dynamics, secondary droplet formation, and effects of pre-existing film thickness

This work investigates surface impingement of monodispersed diesel drop trains using computational simulations. The three-dimensional simulations include high-speed impact of micrometer-sized fuel drops onto initially dry and wetted stainless steel substrates. The drop size and impact velocity are representative of fuel injection conditions in internal combustion engines (ICEs). The drop trains serve as a simplified representation of fuel spray. To accurately capture the film that naturally develops on initially dry substrates, a dynamic diesel-stainless steel contact angle model was developed using experiments of single drop impact under ICE representative temperature. Simulations are presented for a highly splashing case and the effects of a pre-existing film on splashing dynamics were investigated, including the temporal evolution of splashed mass and film thickness. It was concluded that for thin films, the effect of pre-existing film thickness is less significant than in single drop impingement. Secondary droplet characterization was performed on simulation results, leading to instantaneous and time-averaged distributions of secondary droplet size, velocity magnitude, and trajectory angle. It was found that for each drop impingement, approximately 58% of the splashed mass is from the impinging drop itself, while the remainder of the splashed mass is composed of film liquid. From a high resolution case, details of secondary droplet formation are observed and three distinct phases of secondary droplet formation are identified. The detailed analysis of drop train impingement under engine-relevant conditions serves as a first step toward a robust understanding of fuel impingement and the development of cleaner and more efficient ICEs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Correlation of Early-Stage Growth Process Conditions with Dislocation Evolution in MOCVD-Based GaP/Si Heteroepitaxy

To identify the complex relationships between early-stage growth processes and the resultant defect microstructure in GaP/Si heteroepitaxy, a holistic study of several key metal-organic chemical vapor deposition (MOCVD) parameters was conducted, focusing on Si surface preparation and GaP atomic layer epitaxy (ALE) based nucleation processes. Here, crystalline defects related to the lattice mismatch and/or interfacial heterovalency, namely misfit dislocations (MD), threading dislocations (TD), and stacking fault pyramids (SFP), were quantitatively characterized via electron channeling contrast imaging (ECCI) and correlated against the different process variations. Choice of Si surface preparation method between the two examined (dilute SiH4 annealing versus Si2H6 based homoepitaxy) had little impact on resultant GaP film morphology and defect content, whereas differing GaP ALE nucleation conditions produced much more substantial changes. In particular, the initial precursor species (tert-butylphosphine versus triethygallium) and ALE cycle purge times both yielded significant influence over threading dislocation densities (TDD) in thin (100 nm), post-critical thickness GaP/Si films, with TDD spanning two orders of magnitude, from 6.7×10 7 cm -2 to 7.1×10 5 cm -2 , depending on the specific process conditions employed. SFP densities were also found to follow a similar trend, ranging from 2.0×10 7 cm -2 to 1.8×10 5 cm -2 , but with no apparent causal relationship between SFP density and TDD. To help explain the dramatic differences observed, detailed, large-area MD network characterization was used to provide statistically-relevant quantitative analyses of the critical dislocation dynamics (introduction rates and glide velocities) associated with the different process variants. These extracted values are then correlated against the ALE process variants to provide insight into the potential mechanistic roles of the different growth processes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗