SEARCH · Search NASA

Results for “fractured domains”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Quantitative microstructural investigation of 3D-printed and cast cement pastes using micro-computed tomography and image analysis

Highlights: • Processing alters the microstructure in 3D-printed lamellar cubic element. • A unique patterned and connected pore microstructure network exists. • Micro-channels are connected through micro-pores present at interfaces regions. • Three distinct pore volume domains exists in hardened microstructure. • Volume and frequency of pore, hydrated, and unhydrated clusters are dissimilar. Microstructural phases and mechanical properties of lamellar 3D-printed and cast hardened cement paste (hcp) elements were investigated using a lab-based X-ray microscope at two levels of magnification (0.4× and 4×). K-means clustering was used for quantitative image analysis. The entire volume of intact 3-days-old 3D-printed and cast hcp elements was characterized at 0.4× magnification. Three microstructural features (macro-pores, micro-channels, and interfacial micro-pores) were found to reside in three distinct pore size domains. The largest pores of the 3D-printed element were larger than the largest pores of the reference cast hcp element. Moreover, the smallest pore sizes of the 3D-printed element were found to be smaller than those present in the cast counterparts. Micro-channels were found to be connected to one another through the micro-pores present at interfacial regions, indicating the presence of a uniquely patterned and interconnected pore network. The role of locally weak and porous interfaces on mechanical response and fracture properties is discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Distributed Acoustic Sensing to Estimate the Permeability

Optical fiber in a borehole can be interrogated with distributed acoustic sensors (DAS) to capture fracture displacements with the potential to map surrounding fracture networks. We designed a laboratory experiment to test the capability of DAS to determine borehole flow characteristics, and we show that for the first time DAS can be used to remotely estimate permeability. Optical fiber was wrapped around a bead filled pipe and the pressure drop and flow velocity were measured to directly calculate permeability. A machine learning model using statistical features from continuous DAS estimated the bulk permeability. Fluid interactions with the permeable material demonstrate insufficient resolution using DAS amplitude-based measurements for estimating pressure drop to infer permeability. Variations in the spectral domain relate DAS measurements to the pressure drop and provide consistent permeability estimates. Resolution with DAS is sufficient to estimate permeability and provides a reliable method to monitor at depth in borehole conditions.

58 GEOSCIENCES↗

DOE OSTI2025

Shear Response of Ionizable Polymer Melts at the Crossover from Ionomers to Polyelectrolytes

Nonlinear shear response of polymers is affected by inherent barriers for diffusive motion, including entanglements and topology. In ionizable polymers, ionic clusters further constrain the intrinsic dynamics of the polymers, significantly enhancing their viscosity. Here, using fully atomistic molecular dynamics simulations, the nonlinear shear response of ionizable polymers is presented, across the transition from the ionomer regime where distinctive clusters dominate the structure to the polyelectrolyte regime where clusters percolate, in polystyrene randomly sulfonated with fractions of SO 3 − groups of f = 0.20 and 0.35, in pristine and tetrahydrofuran (THF) swollen polyelectrolyte melts. For f = 0.20, the ionic clusters first fracture into smaller clusters followed by splitting into individual ionic groups and eventually reform. At higher f, the clusters morph in shape but do not break under high shear. At very high shear rates, all of the chains stretch and recoil rapidly. As the shear rate is reduced, some chains stretch and recoil, while others remain largely unaffected by the shear. Macroscopically, for all systems, the shear viscosity displays initially an elastic response, followed by nonlinear shear stress overshoot and, eventually, a steady state. The evolution of viscosity with time and shear reflects that of the ionizable domains.

cluster chemistry↗

DOE OSTI2021

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects

Researchers at New Mexico Tech, Los Alamos National Laboratory (LANL), Sandia National Laboratory (SNL), and other collaborators are developing a methodology to improve characterization of stress in the subsurface by means of a model-based inversion of reservoir engineering data, time lapse seismic measurements, and microseismicity. Historically, various geophysical techniques have been used in efforts to understand the state of stress in the subsurface through direct and indirect imaging of stress sensitive features (faults and fractures) and observations of transient stress related observations (time variant elastic moduli and microseismicity). Standard direct and indirect techniques for seismic fault and fracture imaging suffer from detectability limits due to reliance on high data quality and multiplicity. Variations in effective elastic properties from inversion of high quality seismic (Vertical Seismic Profile (VSP)) data may be used to infer stress near the wellbore through integration with independent experimental characterization of the stress-velocity relationship. Given a sufficiently robust observation network, microseismic emissions may be inverted to characterize focal mechanisms which, together with supporting assumptions and constraints, inform estimates of in-situ stress. While each of these techniques contributes in part to the characterization of stress within limited spatial and temporal domains, no one method provides an unambiguous stress measurement or a predictive capability over a site scale spatial extent. Challenges associated with solution non-uniqueness, measurement ambiguity, and irregular sampling may be greatly minimized through combination of one or more of independent measurements within a common framework in which realistic geological, hydrodynamic, geomechanical, and seismological constituent process models may act as constraints.

02 PETROLEUM↗

DOE OSTI2023

Enhancing Lithium-Ion Battery Aging Simulations: Coupling a High-Resolution, 3D, Grain-Scale Electromechanical Model to a Single-Particle Model

One of the main goals in modeling lithium-ion batteries is to improve/predict longevity and resilience of new chemistries. Unfortunately, this requires simulation of thousands of charge/discharge cycles, which can be rather time consuming depending on the fidelity of the simulation. The purpose of this talk is to discuss a new modeling framework that couples a high-resolution, continuous damage model (CDM) to a single particle model (SPM) resulting in a good combination of speed and fidelity. In previous work, a 3D, continuum-level chemo-mechanical model was developed to investigate cracking within a single cathode particle comprised of hundreds to thousands of randomly oriented grains. The CDM predicted that particle fracture is primarily due to non-ideal grain interactions with slight dependence on high-rate charge demands. Essentially, when neighboring grains were misaligned, they expanded different rates relative to one another leading to high stresses and ultimately the formation of intraparticle cracks. The model predicted that small particles with large grains develop significantly less damage than larger particles with small grains. Finally, the model predicted most of the chemo-mechanical damage accumulates in the first charge after formation. This chemo-mechanical "damage saturation" effect indicated that initial particle fracture occurs within the first few cycles, while long-term cathode degradation is not solely chemo-mechanically induced. This led to a need for simulating fatigue-like mechanism that degrade the battery over longer time scales. In order to reach the time scales, need to resolve fatigue-like degradation, the CDM needs to be complemented by a faster model. Therefore, recent efforts have been focused on using results from the CDM to inform parameters within the SPM. These parameters are homogenization factors that are associated with diffusion, particle radius, and/or exchange current density. By coupling the CDM to the SPM, the aging simulation is broken up into two domains: short-term and long-term degradation. The short-term degradation occurs over a single cycle and is handled by the CDM due of its high fidelity, but relatively expensive computational cost. Such mechanism include break-in crack caused by mismatches in grain orientation. The long-term degradation occurs over tens of cycles and is handled by the SPM due it's computational efficiency. Most of the fatigue-like mechanism fall into this category. The eventual goal of this modeling framework is to upscale to a psudo-2D model allowing for full electrode simulation, which are informed by high-fidelity grain-scale simulations.

electrochemistry↗

NASA NTRS1991

Simplified computational methods for elastic and elastic-plastic fracture problems

This paper presents an overview of some of the recent (1984-1991) developments in computational/analytical methods in the mechanics of fracture: (i) analytical solutions for elliptical or circular cracks embedded in isotropic or transversely isotropic solids (the crack-plane being at an arbitrary angle to the axis of transverse isotropy), with crack-faces being subjected to arbitrary tractions; (ii) finite-element or boundary-element alternating methods for two- and three-dimensional crack problems; (iii) a "direct-stiffness" method for stiffened panels with flexible fasteners and with multiple cracks, using the alternating method; (iv) multiplesite-damage near a row of fastener holes; (v) analysis of cracks with bonded repair patches; (vi) methods for generation of weight-functions for 2- and 3-D crack problems, and (vii) domain-integral methods for elastic-plastic or inelastic crack mechanics.

Aging Aircraft↗

NASA NTRS1990

Macromolecular differentiation of Golgi stacks in root tips of Arabidopsis and Nicotiana seedlings as visualized in high pressure frozen and freeze-substituted samples

The plant root tip represents a fascinating model system for studying changes in Golgi stack architecture associated with the developmental progression of meristematic cells to gravity sensing columella cells, and finally to "young" and "old", polysaccharide-slime secreting peripheral cells. To this end we have used high pressure freezing in conjunction with freeze-substitution techniques to follow developmental changes in the macromolecular organization of Golgi stacks in root tips of Arabidopsis and Nicotiana. Due to the much improved structural preservation of all cells under investigation, our electron micrographs reveal both several novel structural features common to all Golgi stacks, as well as characteristic differences in morphology between Golgi stacks of different cell types. Common to all Golgi stacks are clear and discrete differences in staining patterns and width of cis, medial and trans cisternae. Cis cisternae have the widest lumina (approximately 30 nm) and are the least stained. Medial cisternae are narrower (approximately 20 nm) and filled with more darkly staining products. Most trans cisternae possess a completely collapsed lumen in their central domain, giving rise to a 4-6 nm wide dark line in cross-sectional views. Numerous vesicles associated with the cisternal margins carry a non-clathrin type of coat. A trans Golgi network with clathrin coated vesicles is associated with all Golgi stacks except those of old peripheral cells. It is easily distinguished from trans cisternae by its blebbing morphology and staining pattern. The zone of ribosome exclusion includes both the Golgi stack and the trans Golgi network. Intercisternal elements are located exclusively between trans cisternae of columella and peripheral cells, but not meristematic cells. In older peripheral cells only trans cisternae exhibit slime-related staining. Golgi stacks possessing intercisternal elements also contain parallel rows of freeze-fracture particles in their trans cisternal membranes. We propose that intercisternal elements serve as anchors of enzyme complexes involved in the synthesis of polysaccharide slime molecules to prevent the complexes from being dragged into the forming secretory vesicles by the very large slime molecules. In addition, we draw attention to the similarities in composition and apparent site of synthesis of xyloglucans and slime molecules.

NASA Program Space Biology↗

NASA NTRS1994

Analysis of the Effects of Surface Pitting and Wear on the Vibrations of a Gear Transmission System

A comprehensive procedure to simulate and analyze the vibrations in a gear transmission system with surface pitting, 'wear' and partial tooth fracture of the gear teeth is presented. An analytical model was developed where the effects of surface pitting and wear of the gear tooth were simulated by phase and magnitude changes in the gear mesh stiffness. Changes in the gear mesh stiffness were incorporated into each gear-shaft model during the global dynamic simulation of the system. The overall dynamics of the system were evaluated by solving for the transient dynamics of each shaft system simultaneously with the vibration of the gearbox structure. In order to reduce the number of degrees-of-freedom in the system, a modal synthesis procedure was used in the global transient dynamic analysis of the overall transmission system. An FFT procedure was used to transform the averaged time signal into the frequency domain for signature analysis. In addition, the Wigner-Ville distribution was also introduced to examine the gear vibration in the joint time frequency domain for vibration pattern recognition. Experimental results obtained from a gear fatigue test rig at NASA Lewis Research Center were used to evaluate the analytical model.

Choy, F. K.↗

DOE OSTI2021

Comparison of the Mechanical Properties of a Conjugated Polymer Deposited Using Spin Coating, Interfacial Spreading, Solution Shearing, and Spray Coating

The mechanical properties of π-conjugated (semiconducting) polymers are a key determinant of the stability and manufacturability of devices envisioned for applications in energy and healthcare. These properties—including modulus, extensibility, toughness, and strength—are influenced by the morphology of the solid film, which depends on the method of processing. To date, the majority of work done on the mechanical properties of semiconducting polymers has been performed on films deposited by spin coating, a process not amenable to the manufacturing of large-area films. Here, we compare the mechanical properties of thin films of regioregular poly(3-heptylthiophene) (P3HpT) produced by three scalable deposition processes—interfacial spreading, solution shearing, and spray coating—and spin coating (as a reference). Our results lead to four principal conclusions. (1) Spray-coated films have poor mechanical robustness due to defects and inhomogeneous thickness. (2) Sheared films show the highest modulus, strength, and toughness, likely resulting from a decrease in free volume. (3) Interfacially spread films show a lower modulus but greater fracture strain than spin-coated films. (4) The trends observed in the tensile behavior of films cast using different deposition processes held true for both P3HpT and poly(3-butylthiophene) (P3BT), an analogue with a higher glass transition temperature. Grazing incidence X-ray diffraction and ultraviolet–visible spectroscopy reveal many notable differences in the solid structures of P3HpT films generated by all four processes. Herein, while these morphological differences provide possible explanations for differences in the electronic properties (hole mobility), we find that the mechanical properties of the film are dominated by the free volume and surface topography. In field-effect transistors, spread films had mobilities more than 1 magnitude greater than any other films, likely due to a relatively high proportion of edge-on texturing and long coherence length in the crystalline domains. Overall, spread films offer the best combination of deformability and charge-transport properties.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Fundamentals of Scalable Nanoparticle Assembly in Engineering Polymres

Systematically improving the properties of polymer nanocomposites is contingent on effectively controlling nanoparticle (NP) dispersion; further, it is necessary to attain this at scale if these methods are to impact practical applications. While most past work has focused on amorphous polymers, ≈ 70% of all polymers sold annually are semicrystalline. Inspired by this apparent disconnect, our work in the current DOE funded grant and this renewal proposal employs a novel handle, polymer crystallization, to organize NPs into the various amorphous zones of the lamellar semicrystalline morphology (interlamellar, interfibrillar and interspherulitic domains) of stereoregular polymers. We hypothesize that the isothermal crystallization rate of a polymer host, G, and NP diffusivity, are the two relevant control parameters in this context. In the case where the NPs are miscible in the melt, they remain spatially well-dispersed when the polymer is crystallized at a rate that is fast relative to NP diffusion. In contrast, monolayer-thick NP sheets with intersheet spacing of ~10-100nm form in slowly crystallized samples. This sheet stacking is remarkably similar to that seen in natural materials, e.g., Nacre. There is also ordering of the NPs at the larger interfibrillar scale. Importantly, the relative fraction of engulfed, interlamellar and interfibrillar NPs is determined by G. Further, the modulus of these materials can increase by a factor of 2-3 due to the organization state of the NPs, while leaving fracture toughness unaffected. We have established the broad applicability of this approach in two polymers, polyethylene (PE) and polyethyleneoxide (PEO). To achieve this goal we developed a method for grafting long PE chains from NPs with controlled graft density and molecular weight, giving us the ability to tailor their dispersion in the PE melt prior to crystallization. In this renewal proposal, we will develop the fundamental, quantitative understanding that will allow us to independently tailor polymer crystallization rate, NP diffusivity and NP agglomeration state in semicrystalline polymers so as to facilely control NP ordering and hence properties. A program that tightly integrates experiment and theory is proposed, and a key outcome of this research will be a three-dimensional design map for different regimes of NP organization (and hence properties) that can be achieved by polymer crystallization. A second independent outcome is the role of directional solidification and how it competes with the rates of polymer crystallization and NP mobility to affect the organization of the NP.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Multimodal study of the impact of stimulation pH on shale pore structure, with an emphasis on organics behavior in alkaline environments

The tight nature of shale formations calls for hydraulic fracturing techniques able of altering the pore architecture to facilitate the production of hydrocarbon resources. The composition, pH, salinity, density and chemistry of Hydraulic Fracturing Fluids (HFF) vary substantially depending on reservoir characteristics. pH is among the most important stimulation fluid properties, ranging from very acidic, for carbonate reservoirs, to basic, in the case of clay-rich formations. pH regulates what geochemical reactions take place between the stimulation fluid and shale, as well as pore architecture alteration. The dissolution of carbonates and pyrite, and the precipitation of common minerals such as barite, gypsum and iron oxides in acidic environments have been extensively documented. In contrast, alkali stimulation environments and the role of organic components have received less attention. This research provides insight into the role of both minerals and organic components during alkaline stimulation, and the resulting pore architecture alterations. Here, a set of reactive experiments are performed using three different shales, with varying organic matter (OM) content, at different pH (3, 6, 8, 10, 12). The analysis of pore architecture alteration induced by the reaction was performed via nitrogen (N2) gas adsorption and Time-Domain Nuclear Magnetic Resonance (TD-NMR). Inductively Coupled Plasma - Optical Emission Spectrometry (ICP-OES) was used to measure the composition of the stimulation fluid post-reaction. The system with a larger impact on the pore architecture (pH 12) underwent a more extensive analysis through the application of Focused Ion Beam Scanning Electron Microscopy (FIB-SEM) and x-ray Diffraction (XRD). Results show that mineral dissolution contributed to the creation of secondary porosity, mainly through dissolution of clays and silicates. We show that the abundance and distribution of organic matter both play a significant role in changes in pore architecture. Most samples showed an increase in inter-organic porosity, likely due to the organic acids within kerogen reacting with the stimulation fluid. Samples with a larger amount and more widespread distribution of organic matter show the most significant alterations. Finally, we also highlight the dissolution process of iron-rich minerals, including framboidal pyrite, as another important source of secondary porosity.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2022

Utah FORGE Phase 3 Native State Model: 2022 Update

This is the Phase 3 native state model update. The Phase 3 numerical model represents a significant subsurface volume below the FORGE site footprint. The model domain of 4.0 km x 4.0 km x 4.2 km is located approximately between depths of 4000 to 4200 meters below land surface. This data archive consists of 10 files, 4 of which are simulation input files and the remaining 6 are simulation output files. There is an included readme.txt file that contains details on each of the data files. The input files include meshes, FALCON code inputs, tabulated data of water properties, temperature values, and model boundaries. The output files include simulation outfiles and point data of modeled material properties.

15 GEOTHERMAL ENERGY↗

NASA NTRS2009

A Continuum-Atomistic Analysis of Transgranular Crack Propagation in Aluminum

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain is used to study plastic processes at a crack tip in a single crystal of aluminum. The case of mode I loading is studied. A transition from deformation twinning to full dislocation emission from the crack tip is found when the crack plane is rotated around the [111] crystallographic axis. When the crack plane normal coincides with the [112] twinning direction, the crack propagates through a twinning mechanism. When the crack plane normal coincides with the [011] slip direction, the crack propagates through the emission of full dislocations. In intermediate orientations, a transition from full dislocation emission to twinning is found to occur with an increase in the stress intensity at the crack tip. This finding confirms the suggestion that the very high strain rates, inherently present in MD simulations, which produce higher stress intensities at the crack tip, over-predict the tendency for deformation twinning compared to experiments. The present study, therefore, aims to develop a more realistic and accurate predictive modeling of fracture processes.

Yamakov, V.↗

DOE OSTI2024

Inertia-induced mixing and reaction maximization in laminar porous media flows

Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.

Chen, Michael A. (ORCID:0000000294614710)↗

DOE OSTI2025

Computer Vision Pipeline for Image Analysis for Freeze‐Fracture Electron Microscopy: Rosette Cellulose Synthase Complexes Case

In materials science, plant biology, agriculture, and environmental research, the automated analysis of high-magnification, complex microscopy images, such as those generated by freeze-fracture electron microscopy (FF-TEM), remains a critical challenge that limits the scalability of data interpretation. We present a deep learning computer vision pipeline for high-throughput detection and morphological characterization analysis of cellulose synthase complexes (CSCs, or rosettes) in FF-TEM images. The pipeline integrates preprocessing, detection, human-in-the-loop verification, and semantic segmentation to quantify features such as rosette diameter and inter-lobe spacing. The approach was trained and tested on a curated dataset of high-resolution FF-TEM micrographs of Physcomitrium patens, expanded via strategic tiling and augmentation to over 650 images. We compare YOLOv8 and YOLOv9 architectures and demonstrate that YOLOv9 achieves superior performance in both localization accuracy (mAP50-95 = 0.854) and inference speed. The resulting distributions revealed biological variability consistent with prior manual studies, validating the approach for high-throughput applications. Our results show that the pipeline achieves human-expert level accuracy while dramatically reducing analysis time, enabling scalable, reproducible structural characterization of intramembrane protein complexes. The pipeline is broadly applicable to other domains requiring precise interpretation of complex microscopy data and establishes a foundation for future artificial intelligence (AI)-assisted workflows in biological imaging.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Unveiling stimulation fluid-driven alterations in shale pore architecture through combined interpretation of TD-NMR and multi-component gas adsorption

The hydraulic fracturing stimulation process can alter petrophysical and flow properties of the shale matrix in unconventional reservoirs. These alterations can result from interactions with native and exogenous fluids, and can be detrimental to reservoir productivity. Different experimental procedures have been proposed to study reactive processes and thereby understand the aforementioned interactions. Most published works do not consider the rock initial saturation with formation fluid, or the effect of organic matter (OM). To gain insights into these effects, two reactive experimental sets, differing in experimental conditions, were conducted. In the first set, at 125°C and 45 MPa, a stimulation fluid surrogate was put in contact with Baxter Shale samples, previously saturated and equilibrated with synthetic formation fluid. A baseline consisted of exposure to formation fluid exclusively. In the second set, at 0.08 MPa and 25°C, contact with stimulation fluid took place as in the first set. Two additional assays in this second set either skipped the initial formation fluid saturation or removed the OM, respectively. Petrophysical properties analyses, focused on the pore architecture, relied on the combined interpretation of gas adsorption and Time-Domain Nuclear Magnetic Resonance (TD-NMR). Results show that the presence of the formation fluid significantly restricts pore space accessibility to the stimulation fluid and mitigates the pore-architecture alteration to the macro-pore region. Mineral dissolution and precipitation dominate over surface or wettability alteration. Next, results suggest that the OM restricts accessibility to some pores. Ultimately, the equilibration process with synthetic formation fluid drives mainly surface alteration, unhindered by OM.

58 GEOSCIENCES↗

DOE OSTI2021

Predictive Data-driven Platform for Subsurface Energy Production

Subsurface energy activities such as unconventional resource recovery, enhanced geothermal energy systems, and geologic carbon storage require fast and reliable methods to account for complex, multiphysical processes in heterogeneous fractured and porous media. Although reservoir simulation is considered the industry standard for simulating these subsurface systems with injection and/or extraction operations, reservoir simulation requires spatio-temporal “Big Data” into the simulation model, which is typically a major challenge during model development and computational phase. In this work, we developed and applied various deep neural network-based approaches to (1) process multiscale image segmentation, (2) generate ensemble members of drainage networks, flow channels, and porous media using deep convolutional generative adversarial network, (3) construct multiple hybrid neural networks such as convolutional LSTM and convolutional neural network-LSTM to develop fast and accurate reduced order models for shale gas extraction, and (4) physics-informed neural network and deep Q-learning for flow and energy production. We hypothesized that physicsbased machine learning/deep learning can overcome the shortcomings of traditional machine learning methods where data-driven models have faltered beyond the data and physical conditions used for training and validation. We improved and developed novel approaches to demonstrate that physics-based ML can allow us to incorporate physical constraints (e.g., scientific domain knowledge) into ML framework. Outcomes of this project will be readily applicable for many energy and national security problems that are particularly defined by multiscale features and network systems.

58 GEOSCIENCES↗

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↗