Search NASA⌕ Search

SEARCH · Search NASA

Results for “volume fraction”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗
NASA NTRS1990

Free volume in contact lens polymers

We have measured the free volume and gas permeability in several contact lens specimens as part of a Space Commercialization Program. Free volume was measured using positron lifetime spectroscopy, while permeability for O2, N2, and CO2 gases was measured using mass spectrometry and polarography. It has been noted that permeability for all gases increases with the mean free volume cell size in the test samples. As might be expected, the specimens with the highest free volume fraction also exhibit the lowest Rockwell Hardness Number. An interesting corollary of this study is the finding that the presence of fluorine atoms in the lens chemical structure inhibits filling up of their free volume cells. This is expected to allow the lenses to breathe freely while in actual use.

Singh, Jag J.↗
DOE OSTI2025

Porosity, swelling, and composition evolution in high-burnup monolithic U-Mo fuel

The microstructural progression of very high-burnup (>8 × 10 21 fissions/cm 3 ) monolithic uranium-molybdenum (U-10wt.%Mo) was analyzed, providing crucial insights into the behavior of post-recrystallized nuclear fuel, where scant data exists. Three focused ion beam cuboids sourced from a fuel plate with varying local burnups of 8.86 × 10 21 , 9.05 × 10 21 , and 9.36 × 10 21 fissions/cm 3 were characterized. The porosity and composition of the samples were evaluated to characterize the evolution of the microstructure as a function of fission density and different locations on the fuel plate, while simultaneously isolating plate-specific parameters such as Zr diffusion barrier thickness, hot-isostatic-pressing conditions, enrichment, and reactor conditions. The porosity was segmented, and the three-dimensional distribution of the porosity was extrapolated from the two-dimensional segmentation. The composition was assessed and quantified using energy-dispersive X-ray spectroscopy areal mapping. The porosity fraction increased as a function of the burnup from 27.77±0.51, 35.12±1.54, and 37.71±0.44 % for 8.86 × 10 21 , 9.05 × 10 21 , and 9.36 × 10 21 fissions/cm 3 , respectively. When compared to literature, the porosity volume fraction plateaus at burnups greater than 6 × 10 21 fissions/cm 3 , while the pore size grows linearly as a function of fission density. The number of large pores increased in number density as a function of burnup, while the smallest pores (<0.3 µm) increased up to 9.05 × 10 21 fissions/cm 3 , followed by a decrease at 9.36 × 10 21 fissions/cm 3 . The delamination and cracking in the fuel plate propagated through an interconnected porosity sublayer identified ∼5 µm from the diffusion barrier. The local swelling of the specimens was within or near the prediction bounds of the Robinson-Williams model for local swelling. The fission products, strontium, barium, cerium, and cesium, precipitated into the pores, while neodymium accumulated adjacent to the pores. Furthermore, these findings have direct implications for the development of fuel performance codes and the accurate documentation of the microstructure evolution in high burnup U-Mo, thus enhancing the safety and efficiency of nuclear fuel usage.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1991

1200 and 1300 K slow plastic compression properties of Ni-50Al composites

XD synthesis, powder blending, and hot pressing techniques have been utilized to produce NiAl composites containing 4, 7.5, 15, and 25 vol pct alumina whiskers and hybrid composite materials with 15 vol pct Al2O3 + 10 or 20 vol pct, nominally 1 micron TiB2 particles. The resistance to slow plastic flow was determined at 1200 and 1300 K via compression testing in air under constant velocity conditions. The stress-strain behavior of the intermetallic composites depended on the fraction of second phases where the 4 and 7.5 percent Al2O3 materials flowed at a nominally constant stress after about 2 percent deformation, while all the other composites exhibited diffuse yielding followed by strain softening. The flow stress-strain rate properties increased with volume fraction of Al2O3 whiskers except for the 4 and 7.5 percent materials, which had similar strengths. The hybrid composite NiAl + 15Al2O3 + 10TiB2 was substantially stronger than the materials simply containing alumina. Deformation in these composites can be described by the Kelly and Street model of creep in perfectly bonded, rigid, discontinuous fiber materials.

Whittenberger, J. D.↗
NASA NTRS2019

Advances in Low-Cost Manufacturing and Folding of Solar Sail Membranes

Solar sail membranes must have a high area-to-mass ratio and high solid volume fraction when stowed. In order to meet mission requirements, current solar sail projects, such as NASA’s Near Earth Asteroid Scout, require metallized sail membranes with thicknesses on the order of 2-3 μm. These very thin membranes do not retain creases like thicker membranes, solar panels, or paper models. For Cubesat-class spacecraft, volume, rather than mass, is often the driving requirement for deployable structural elements. These two factors make it both difficult and highly desirable to characterize the practical differences between solar sail membrane packaging methods with laboratory demonstrations. This paper presents lessons gathered from lab work with solar sail membranes at a 10-meter scale.

Stohlman, Olive R.↗
NASA NTRS2001

Computational Thermodynamic Study to Predict Complex Phase Equilibria in the Nickel-Base Superalloy Rene N6

A previous study by Ritzert et al. on the formation and prediction of topologically closed packed (TCP) phases in the nickel-base superalloy Rene' N6 is re-examined with computational thermodynamics. The experimental data on phase distribution in forty-four alloys with a composition within the patent limits of the nickel-base superalloy Rene' N6 provide a good basis for comparison to and validation of a commercial nickel superalloy database used with ThermoCalc. Volume fraction of the phases and partitioning of the elements are determined for the forty-four alloys in this dataset. The baseline heat treatment of 400 h at 1366 K was used. This composition set is particularly interesting since small composition differences lead to dramatic changes in phase composition. In general the calculated values follow the experimental trends. However, the calculations indicated no TCP phase formation when the experimental measurements gave a volume percent of TCP phase less than 2 percent. When TCP phases were predicted, the calculations under-predict the volume percent of TCP phases by a factor of 2 to 8. The calculated compositions of the gamma and gamma' phases show fair agreement with the measurements. However, the calculated compositions of the P Phase do not agree with those measured. This may be due to inaccuracies in the model parameters for P phase and/or issues with the microprobe analyses of these phases. In addition, phase fraction diagrams and sigma and P phase solvus temperatures are calculated for each of the alloys. These calculations indicate that P phase is the primary TCP phase formed for the alloys considered here at 1366 K. Finally, a series of isopleths are calculated for each of the seven alloying elements. These show the effect of each alloying element on creating TCP phases.

Copland, Evan H.↗
NASA NTRSDate not supplied

Advances in Low-Cost Manufacturing and Folding of Solar Sail Membrane

Solar sail membranes must have a high area-to-mass ratio and high solid volume fraction when stowed. In order to meet mission requirements, current solar sail projects, such as the NASA Near Earth Asteroid Scout and Advanced Composite Solar Sail System, require metallized sail membranes with thicknesses on the order of 2-3 μm. These very thin membranes do not retain creases like thicker membranes, solar panels, or paper models. For CubeSat-class spacecraft, volume, rather than mass, is often the driving requirement for deployable structural elements. These two factors make it both difficult and highly desirable to characterize the practical differences between various membrane designs and packaging methods with laboratory demonstrations. This paper presents lessons gathered from lab work on solar sail membranes at a 9-meter scale. Among other results, laboratory deployment of a full-scale sail quadrant revealed important differences between seaming methods.

Olive R Stohlman↗
DOE OSTI2026

Influence of Rigidity–Hydration Coupling on Size-Dependent Diffusion in Hydrated Polymer Membranes

Selective ion transport in polymer membranes depends critically on how penetrant motion couples to polymer dynamics and hydration. Yet, the mechanistic interplay between polymer rigidity, water content, and penetrant size remains poorly understood, especially in the regime where the penetrant diameter, polymer Kuhn length, and correlation length are comparable. Here, we employ coarse-grained molecular dynamics simulations to systematically investigate penetrant diffusion in hydrated polymer networks across a broad range of water volume fractions, chain rigidities, and penetrant sizes. The results reveal a transition from a decoupled regime, where small penetrants diffuse nearly independently of polymer relaxation, to a coupled regime in which large penetrants require cooperative polymer motion for transport. Increasing polymer rigidity amplifies the sensitivity of diffusivity to hydration, particularly at low water content, leading to pronounced deviations from Stokes−Einstein scaling. Comparison with scaling theories and free-volume models shows that classical nanoparticle-based frameworks fail to capture this intermediate regime. To address this gap, we extend the Yasuda model to incorporate polymer rigidity through a single parameter that quantifies the dynamic contribution of chain stiffness to free-volume fluctuations. The resulting model collapses diffusivity data across all sizes, water contents, and rigidities, providing a unified description of penetrant transport in hydrated polymer matrices. Furthermore, these findings establish polymer rigidity as a key, tunable determinant of diffusion and offer a framework for interpreting size-dependent transport in ion-selective membranes.

diffusion↗
NASA NTRS2007

In-Plane Cracking Behavior and Ultimate Strength for 2D Woven and Braided Melt-Infiltrated SiC/SiC Composites Tensile Loaded in Off-Axis Fiber Directions

The tensile mechanical properties of ceramic matrix composites (CMC) in directions off the primary axes of the reinforcing fibers are important for architectural design of CMC components that are subjected to multi-axial stress states. In this study, 2D-woven melt-infiltrated (MI) SiC/SiC composite panels with balanced fiber content in the 0 degree and 90 degree directions were tensile loaded in-plane in the 0 degree direction and at 45 degree to this direction. In addition, a 2D triaxially-braided MI composite panel with balanced fiber content in the plus or minus 67 degree bias directions and reduced fiber content in the axial direction was tensile loaded perpendicular to the axial direction tows (i.e., 23 degrees from the bias fibers). Stress-strain behavior, acoustic emission, and optical microscopy were used to quantify stress-dependent matrix cracking and ultimate strength in the panels. It was observed that both off-axis loaded panels displayed higher composite onset stresses for through-thickness matrix cracking than the 2D-woven 0/90 panels loaded in the primary 0 degree direction. These improvements for off-axis cracking strength can in part be attributed to higher effective fiber fractions in the loading direction, which in turn reduces internal stresses on critical matrix flaws for a given composite stress. Also for the 0/90 panel loaded in the 45 degree direction, an improved distribution of matrix flaws existed due to the absence of fiber tows perpendicular to the loading direction. In addition, for the +67/0/-67 braided panel, the axial tows perpendicular to the loading direction were not only low in volume fraction, but were also were well separated from one another. Both off-axis oriented panels also showed relatively good ultimate tensile strength when compared to other off-axis oriented composites in the literature, both on an absolute strength basis as well as when normalized by the average fiber strength within the composites. Initial implications are discussed for constituent and architecture design to improve the directional cracking of SiC/SiC CMC components with MI matrices.

Morscher, Gregory N.↗
NASA NTRS2018

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

The Porous Microstructure Analysis (PuMA) software is a suite of tools for the analysis of porous materials and generation of material microstructures. From microstructural data, often obtained through X-ray microtomography, PuMA can determine a number of effective material properties and perform material response simulations. Version 2.2 includes capabilities for computing volume fractions, porosity, specific surface area, effective thermal and electrical conductivities, and continuum and rarefied diffusive tortuosity. PuMA can also simulate competitive diffusion/reaction processes at the micro-scale, such as surface oxidation. In this poster, recent advancements to the PuMA software are detailed, including the full refactoring of PuMA into v3.0, a new module to compute heat conduction in anisotropic materials, a particle method for simulating molecular beam experiments, a new finite-volume Laplace solver, complex fibrous material generation, woven material generation, and a coupling of PuMA with the DAKOTA software for advanced statistics.

PuMA↗
NASA NTRS2006

Micromechanics Modeling of Functionally Graded Interphase Regions in Carbon Nanotube-Polymer Composites

The effective elastic properties of a unidirectional carbon fiber/epoxy lamina in which the carbon fibers are coated with single-walled carbon nanotubes are modeled herein through the use of a multi-scale method involving the molecular dynamics/equivalent continuum and micromechanics methods. The specific lamina representative volume element studied consists of a carbon fiber surrounded by a region of epoxy containing a radially varying concentration of carbon nanotubes which is then embedded in the pure epoxy matrix. The variable concentration of carbon nanotubes surrounding the carbon fiber results in a functionally graded interphase region as the properties of the interphase region vary according to the carbon nanotube volume fraction. Molecular dynamics and equivalent continuum methods are used to assess the local effective properties of the carbon nanotube/epoxy comprising the interphase region. Micromechanics in the form of the Mori-Tanaka method are then applied to obtain the global effective properties of the graded interphase region wherein the carbon nanotubes are randomly oriented. Finally, the multi-layer composite cylinders micromechanics approach is used to obtain the effective lamina properties from the lamina representative volume element. It was found that even very small quantities of carbon nanotubes (0.36% of lamina by volume) coating the surface of the carbon fibers in the lamina can have a significant effect (8% increase) on the transverse properties of the lamina (E22, k23, G23 and G12) with almost no affect on the lamina properties in the fiber direction (E11 and v12).

Seidel, Gary D.↗
DOE OSTI2025

Soot formation and precursor chemistry in Counterflow flames of aviation fuel surrogates

To meet market demands, the aviation sector is interested in utilizing drop-in Synthetic Aviation Turbine Fuels (SATF), either as neat fuels or in blends with conventional Jet A. SATF currently approved in standard specifications may have lower aromatic content with significant fractions of normal, branched, and cyclo-alkanes. Fundamental studies on soot formation from aviation fuels (Jet A, SATF) and their surrogate components are essential to understand how fuel composition influences soot and soot precursor formation. Here, this study reports new measurements of polycyclic aromatic hydrocarbons (PAH) and soot in counterflow diffusion flames (CDFs) of aviation fuel surrogates. Both intrusive and non-intrusive diagnostics are employed to determine the profiles of temperature, gas phase species, PAHs (up to C16), and soot volume fraction (SVF) in CDFs of iso-octane and surrogate mixtures. These measurements shed light on the transition of soot precursors to primary soot particles. In addition to serving as a common surrogate component in Jet A surrogate mixtures, iso-octane is a template species for larger, less volatile branched alkanes found in SATF mixtures. The newly developed Lawrence Livermore National Laboratory (LLNL) PAH and soot model successfully captures temperature, precursor species, and SVF profiles for the mixtures and conditions discussed in this work. Finally, a high-fidelity surrogate for Jet A is proposed that matches targeted physical and chemical properties well, while leveraging the wide range of candidate fuel molecules available in the LLNL detailed chemical model. The proposed surrogate formulation is validated against newly acquired measurements of the surrogate and literature measurements of Jet A. These new experiments and simulations provide critical insights into the PAH and soot formation from aviation fuels. Reaction pathways which require further investigation are highlighted, such that future work may bridge the remaining quantitative gaps in predicting soot formation from aviation fuel surrogates and surrogate components.

Aviation fuels↗
NASA NTRS2018

Conductivity of a Copper-Carbon Covetic Composite

NASA is currently engaged in the development of highly electrically conductive wire for space and aviation power and propulsion applications. Higher levels of electrical conductivity and lower densities compared to copper and aluminum are needed to improve the efficiency of commercial aviation and space flight. Carbon, in various forms, is being considered as a component in composites to increase conductivity and decrease density in Cu and Al. One form of carbon that is currently receiving much interest is known as covectic. The technology currently known as covetics was first introduced by Shugart and Scherer whose patents claimed "Carbon can be incorporated into any of these metals by melting the metal, mixing the carbon into the molten metal and, while mixing, applying a current of sufficient amperage to the molten mixture such that the carbon becomes incorporated into the metal, thereby forming a single phase metal-carbon material." It is reported that in this process a conversion reaction occurs with the carbon forming covalent bonds with the host metal, and that covetic processing enables higher carbon solubility and improvements in properties. Although there is a growing body of work supporting covetic materials, it is still a very new and somewhat controversial field. Composition, density, and electrical conductivity measurements of a covetically processed alloy ingot of Cu and C, with a nominal composition of Cu-0.4wt.% C, were made. A pure Cu ingot was similarly processed for comparison. Sections of these ingots were worked into wire and tested. Only small amounts of C could be detected in the Cu-0.4wt.% C ingot, and trace amounts (0.006 wt.%) in the wire. All wires had densities very close to the full density of pure Cu, indicating a high level of Cu purity and low void fraction. Finally, the electrical conductivity of the Cu-0.4wt.% C wires were 9% lower than the samples similarly processed but without C additions. The wires made without C additions had conductivities approximately equal to pure Cu. The low conductivity of the Cu-0.4wt.% C material is believed to be due to: (1) a structure of interfaces in the wire created by collapsed voids lined with elemental C, (2) impurities, such as Fe in the Cu, which came from the added C, and (3) a small volume fraction of voids.

Degroh, Henry C.↗
NASA NTRS1984

SPAR X Technical Report for Experiment 76-22 Directional Solidification of Magnetic Composites

The effects of gravity on Bridgman-Stockbarger directional solidification of off-eutectic Bi/MnBi were studied in reduced gravity aboard the SPAR X flight and compared to normal-gravity investigations and previous eutectic Bi/MnBi SPAR flight experiments. The directional solidification of off-eutectic Bi/MnBi results in either a dendritic structure connected with local cooperative growth or a coupled low volume fraction faceted/non faceted aligned rod eutectic whose Mn macrosegregation, MnBi rod size, interrod spacing, and thermal and magnetic properties are sensitive functions of the solidification processing conditions. Two hypoeutectic and two hypereutectic samples were solidified during 605 sec of furnace travel, with an initial 265 sec low-gravity interval. Comparison Earth-gravity samples were solidified in the same furance assembly under identical processing conditions. Macrosegregation in the low-g samples was consistent with a metastable increase in Mn solubility in the Bi matrix, in partial agreement with previous Bi/MnBi SPAR findings of MnBi volume reduction.

Bethin, J.↗
DOE OSTI2025

The effect on soot and its gas precursors of doping ethylene with 2,2,4,6,6-pentamethyl-heptane in the nitrogen-fuel stream of a laminar non-premixed Planar Mixing Layer Flame (PMLF)

Synthetic Aviation Turbine Fuels (SATFs) are promising for reducing soot emissions from the aviation sector and diversifying Jet Fuel (JF) sources. Accurately predicting the combustion and emissions behavior of SATFs (and other JFs) necessitates robust experimental databases to elucidate the chemistry of long-chain iso-paraffins, which can compose up to two-thirds of SATF blends and whose behavior is considered to be well-represented by that of iso-dodecane isomers. Here, this study characterizes two laminar non-premixed Planar Mixing Layer Flames (PMLFs) with mild soot loads fueled by nitrogen-diluted ethylene, pure and doped with 2,2,4,6,6-pentamethyl-heptane, respectively. The two PMLFs have the same stoichiometric mixture fraction and total hydrocarbon mole fraction in the fuel stream (X F,F =X C2H4,F +X C12H26,F = 0.260), resulting in nearly the same maximum temperature (T max ≈1800 K) and simple identification of the effects of doping. Importantly, any horizontal PMLF cross-section has a self-similar structure that can be modeled as an equivalent One-Dimensional Counterflow Flame (1D-CF) with vanishingly small strain rate (a). The cross-section at a Height Above the Burner (HAB) of 50 mm is characterized in terms of C 0 -C 18 gas species using capillary sampling followed by GC-MS analyses. Laser-Induced Emission Spectroscopy (LIES) quantifies the soot volume fraction (ƒ v ) profiles at HAB=25 and 50 mm where Elastic Laser Light Scattering (E-LLS) is performed to determine the a of the equivalent 1D-CFs and the profile of the E-LLS equivalent diameter ( d 6,3 ) of soot. The substitution of 1500 ppm of ethylene with 2,2,4,6,6-pentamethyl-heptane causes an increase of ≈1.5 in the concentrations of several polycyclic aromatic hydrocarbons and fv. Concurrently, the measured d 6,3 doubles in the oxidizer stream, yet remains the same in the fuel stream, at HAB=50 mm. Instead, at HAB= 25 mm, the iso-dodecane doping does not affect the d 6,3 profile in either stream. The experimental results partially validate the chemical reactions and soot formation kinetic model developed at Lawrence Livermore National Laboratory and provide directions to further improve its predictions.

Iso-dodecane (2,2,4,6,6-Pentamethyl-Heptane)↗
DOE OSTI2024

Crystallization of Bottlebrush Statistical Copolymers of Polystyrene and Poly(ethylene oxide)

Here, bottlebrush statistical copolymers (BSCPs) with poly(ethylene oxide) (PEO) and polystyrene (PS) side chains grafted to a polynorbornene (PNB) backbone were synthesized by ring-opening metathesis polymerization (ROMP). The impact of the glassy PS side chains on the crystallization of the PEO side chains as a function of the backbone length, grafting densities, and fraction of the PS and PEO side chains is described. The bottlebrush architecture, where the side chains are anchored to the backbone, inherently constrains the mobility of PEO. Compared with the bulk crystallization temperature of PEO, the higher glass transition temperature of PS places further constraints on PEO crystallization. The limited mobility of PEO leads to cold crystallization behavior during heating. Reduced grafting densities, in turn, reduce side-chain crowding, leading to less extended pendant structures. The degree of crystallinity of PEO was found to decrease at lower grafting densities to a point where crystallization was not observed. The average distance between bottlebrush backbones increased linearly with the backbone length, suggesting that the backbone forms a distinct mesodomain. For BSCPs with asymmetric volume fractions of PS and PEO side chains, the degree of crystallinity of PEO increases with a change from cold crystallization to normal nucleation and growth during cooling, as a result of reduced constraints on PEO mobility by the glassy PS side chains.

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

A Multiscale Inelastic Internal State Variable Corrosion Model

We present a corrosion internal state variable (ISV) damage model based upon the integrated computational materials engineering (ICME) hierarchical multiscale paradigm. Structure–property experiments for magnesium alloys were used where the only inputs were the volume fractions of each element of the periodic table. This macroscale ISV corrosion model finds its basis in Horstemeyer’s mechanical damage model, which includes three separate ISVs for damage nucleation, growth, and coalescence, as well as Walton’s inclusion of corrosion, which introduces five new ISVs for pit nucleation, growth, and coalescence, along with general corrosion and intergranular corrosion. While Walton’s corrosion ISVs are phenomenological in nature, herein we develop a multiscale physical basis for the corrosion ISVs. The parameters for the macroscale corrosion ISVs were garnered from the mesoscale Butler–Volmer equations. Pure magnesium with differing amounts of aluminum were used in corrosion tests to exemplify the different pitting, general corrosion, and intergranular corrosion rates, and the macroscale ISV model was calibrated with said data, in which the only inputs to the model are the volume percentages of the elements magnesium and aluminum. Although magnesium alloys were used to motivate and calibrate the model, the model is abstract enough to possibly capture other material systems as well.

Chemistry↗