Search NASA⌕ Search

SEARCH · Search NASA

Results for “fracture geometry”

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

Benchmark study of a new simplified DFN model for shearing of intersecting fractures and faults

It is challenging to quantitatively predict shearing of intersecting fractures/faults because of dynamic frictional contacts accompanied by possible nonlinear rock deformation. To address such challenges, a new conceptual model—the simplified DFN model—was proposed and validated by Hu et al. 46 to use major paths (MPs) to represent complicated DFNs for calculation of shearing. In this work, we conducted a benchmark study for three examples that involve different levels of complexity of intersecting fractures, and correspondingly different numbers of MPs. The codes and software that were used in the benchmark cover a range of continuum, discontinuum and hybrid numerical methods: NMM (LBNL), FLAC3D (LBNL), GBDEM (KIGAM), FRACOD (DynaFrax), and CASRock (CAS). The general consistency between DFN and MP cases as predicted by all the codes/software demonstrates that major paths can be used to simplify the geometry of DFNs in a wide range of software. Disagreement in results made by some software and potential future improvements are discussed. We show that (1) shearing of one or multiple major fractures can be reduced if there are multiple smaller intersecting fractures in that area, which is a useful basis for understanding and controlling induced seismicity and merits further analysis, and (2) the agreement achieved in the benchmark examples provide confidence that the simplified DFN model is a promising conceptual model that can be used for different types of numerical approaches and software for simplifying the analysis of the shearing of intersecting fractures and faults.

58 GEOSCIENCES↗
DOE OSTI2022

Directional Cooling-Induced Fracturing: Directional Cooling Data Under True-Triaxial Stress State

Data includes Directional Cooling-Induced Fracturing (DCIF) testing data using westerly granite blocks. This submission includes data from two samples of westerly granite, lab sample 7 and 8. Files contain stress, temperature and acoustic emission data acquired during polyaxial, laboratory testing of westerly granite blocks for each sample. FILES: .tradb -- files containing acoustic emission waveforms; sqlite3 database .pridb -- files containing basic acoustic emission information (no waveforms); sqlite3 database .geom -- geometry of the AE sensor network (ASCII)

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

DEM analysis of the effect of interface transition zone on dynamic splitting tensile behavior of high-strength concrete based on multi-phase model

Highlights: • Crushable aggregates with the real geometry shape are generated by using clump-cluster method. • Flat-Joint Model is adopted as the constitutive law. • The tensile strength is almost not affected by ITZ strength when its strength exceeds 0.8 times that of mortar. • A modified formula that can reflect the functional relationship between the ITZ strength and the TDIF is proposed. The tensile strength is one of the main parameters in the design of concrete structures and has an important influence on the cracking performance and durability of concrete materials. As the weak region in concrete, the strength of the interface transition zone (ITZ) directly affects the dynamic mechanical response of concrete. Meanwhile, considering that the tensile performance of concrete is far weaker than its compressive performance, the study of the strength of ITZ is of great significance to the dynamic tensile fracture characteristics and mechanical properties of concrete. In this paper, a method combining 3D scanning technology and Clump-Cluster was proposed to generate a real geometry shape model of crushable aggregate. On this basis, the 3D discrete element model (DEM) of high-strength concrete was constructed using the Flat-Joint Model (FJM). The effects of the ITZ strength on the failure mode, number of microcracks, tensile strength and strain rate effect of high-strength concrete under strain rates from 0.01/s to 100/s were studied by splitting tensile test. The simulation results indicate that the number of microcracks generated decreases with the increase of ITZ strength, and the proportion of ITZ cracks decreases obviously. When the ITZ strength exceeds the critical value of 0.8 times that of mortar, the tensile strength is almost not affected by ITZ strength. In addition, based on the numerical simulation results, a modified formula that can reflect the functional relationship between the ITZ strength and the dynamic increasing factor of the concrete tensile strength (TDIF) is proposed, and the rationality of the TDIF modified formula is verified.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanical evaluation of WC-Co materials with varying microstructures

Tungsten carbide–cobalt materials are useful in a variety of extreme applications due to a desirable blend of properties, yet the technology has not significantly changed since their initial development in the 1920s. The mechanical properties of this class of materials is highly dependent on two variables, the size of the tungsten carbide grains, and the amount of binder phase present in the final body. In this study, the amount of binder phase is isolated across three commercial materials from the same manufacturer with three different grain sizes to investigate the effect on mechanical properties. The mechanical properties investigated are indentation hardness, flexure and tensile strength, as well as fracture toughness. In general, an increase in hardness and tensile strength with decreasing grain size was observed, while the fracture toughness showed the opposite trend with toughness increasing with increasing grain size. The flexure strength results did not show a correlation to grain size. Fractographic analysis identified the dominant strength-limiting flaw for each sample, which largely were in the form of porosity. Other flaws types, such as inclusions from the milling process, clusters of large grains, and machining cracks from the surface finishing process, were also identified. Finally, Weibull analysis was performed and deemed appropriate for analysis of these materials, but strength-size scaling was not conducted due to the variability in the strength-limiting feature between the different specimen geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Intrinsic and extrinsic effects on the fracture toughness of ductile metallic glasses

The potential application of metallic glasses (MGs) as structural materials demands accurate measurements of their toughness and the understanding of the underlying factors affecting it. Currently, it is challenging to precisely measure the toughness of MGs, especially ductile MGs. Measured toughness values are widely scattered, even for MGs with identical compositions. That is attributed to the combined effect of intrinsic and extrinsic factors including processing, sample geometries, and loading conditions. A fundamental understanding of the influences of these elements is thus of great significance. In the present study, molecular dynamics simulations are performed to investigate the influence of intrinsic and extrinsic effects on the fracture toughness of CuZr MGs. In particular, focused is placed on the effects of cooling rate, notch shape and size, strain rate, and temperature. The results indicate that the fracture toughness of a MG scales with the cooling rate used to prepare it. This is attributed to increased free volume content generated at high cooling rates, which enhances plastic deformation and amplifies the associated energy dissipation during plastic deformation events. Fracture toughness also increases with the strain rate, arguably due to strain rate-induced crack extension delay. Overall, the results demonstrate that the largest fracture toughness are achieved when MG samples are fabricated at high cooling rates and subjected to high strain rate deformation. In addition, results suggest that the fracture toughness decreases with increasing temperature, due to the significant decrease in strength. Here, the correlations revealed between these crucial intrinsic and extrinsic parameters and the calculated MG fracture toughness support the development of a framework to understand the root of the discrepancies in the measurement of the toughness of MGs and provide insights into the design of tough MGs for structural applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A reactive tracer method for predicting EGS reservoir geometry and thermal lifetime: development and field validation

The project summarized here (DOE Award No. DE-EE0006764) was intended to develop a methodology for predicting advective heat transfer in fracture-dominated crystalline rock. Our goal was to determine if a combination of inert and reactive tracers could adequately constrain the effective heat transfer surface area between an injection-production well pair. Our approach consisted of: 1. developing a novel computational framework; 2. Performing heat and tracer experiments at meso-scale; and 3. Comparing predictions of advective heat transfer to the “true” thermal breakthrough measured at the Altona site. Below is a summary of project activities/findings, a summary of project tasks, and a conclusion

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Flow Channeling in Fracture Networks: Characterizing the Effect of Density on Preferential Flow Path Formation

Flow channelization is a commonly observed phenomenon in fractured subsurface media where the flow of fluids is restricted primarily to highly transmissive fracture networks surrounded by a low-permeability rock matrix. The multiscale structural heterogeneity of these networks results in multiscale flow channelization where preferential flow paths form at length scales ranging from the entire system down to the subfracture size. We present an analysis of how one of the largest scales in fractured media, the network density, influences the degree of flow channeling that occurs using an ensemble of semigeneric three-dimensional discrete fracture network (DFN) simulations. We construct 10 DFNs, whose fracture lengths follow a power law distribution, at four densities for a total of 40 networks. We characterize their structure in terms of the network topology and geometry. Eulerian and Lagrangian observations of the steady-state flow fields obtained within the networks are used to quantify the degree of flow channelization at the network scale. We introduce a measure for the importance-ranking/hierarchy of different flow paths in the network using graph-based analysis of Lagrangian transport by which the degree of flow channeling between networks is compared. These flow observations are then linked to the structural properties of the networks. In general, network-scale flow channeling decreases as the network density increases. However, at low densities, there is more uniform flow within the entire connected network than in high-density networks. We also demonstrate how standard transport observables can be used to infer the degree of flow channelization occurring within a fracture network.

58 GEOSCIENCES↗
DOE OSTI2024

Optimal Geometry for Focused Ion Beam-Milled Samples for Direct-Pull Micro-Tensile Testing Performed In Situ in a Scanning Electron Microscope

A thorough procedure was developed to efficiently manufacture dogbone samples using focused ion beam (FIB) milling for micro-tensile testing. A Bruker PI 89 PicoIndenter, Billerica, MA, USA, was used as a case study, although the analysis and results are applicable to other micro-mechanical testing systems capable of mounting a standard, Ø12.7 mm × Ø3.2 mm pin, scanning electron microscopy (SEM) pin stub (Ted Pella, Redding, CA, USA). Nine dogbones were made from an Fe-45Cu alloy additively manufactured using powder-fed laser-directed energy deposition (DED-LB). Testing showed that fracture was confined to the gauge section for all dogbones and that the fracture mode, ductile vs. brittle, was entirely dependent on the grain orientation relative to the loading direction. The analysis showed that the measured plastic strain to failure can vary from >11% (optimal geometry) to <1% (non-optimal geometry) in micro-tensile testing of high-tensile-strength (>1 GPa) metallic materials. Subsequently, a finite element analysis (FEA) was conducted to identify the improved dogbone geometries. A total of ten thousand dogbone geometries were tested, and their dimensions were defined by a set of four adjustable parameters (corner radius, load surface angle, load surface length, and dogbone head length). The gauge width and gauge length were fixed to 4 µm and 10 µm, respectively. Three-dimensional surface plots of the stress concentration as a function of two parameters were used to identify the optimal ranges of parameter values. The addition of maximum width and length constraints, measuring 25 µm and 30 µm, respectively, allowed us to identify an optimal geometry at load surface angles of 30° and 45°. Their respective dimensions (corner radius, load surface length, and dogbone head length) are, in µm, 12, 6, and 7 and 10, 7, and 7. Testing these two optimal geometries with a range of gauge lengths from 4 to 20 µm showed that smaller gauge lengths only slightly reduced the detrimental stress concentration outside the gauge section. However, smaller gauge lengths will notably improve the FIB surface polishing step as tapering is reduced with smaller dogbone lengths.

Chemistry↗
DOE OSTI2022

Investigation of fracture in porous materials: a phase-field fracture study informed by ReaxFF

Microscopic features (e.g, pore shapes, sizes, and distribution) in porous material substantially affect the overall mechanical properties such as stiffness and strength. In turn, these material properties determine the macroscopic behaviors of fracture in the porous material. In certain cases, macroscopic properties can be derived from the porous skeleton and void ratio (i.e., porosity), but in many other cases, such derivation is a challenging task. This paper presents a numerical investigation of microporosity and micropore shapes effect on the macrofracture behavior in porous amorphous silica. For this study, we extend the recently-proposed combined molecular dynamic (MD) and phase-field (PF) fracture modeling approach by including different pore shapes in the atomistic domain. In the MD simulations, we adopt ReaxFF to evaluate the material properties, where four different micropore cases are considered. Based on the material properties derived from MD simulations, the macrofracture propagation of porous media is studied using hybrid PF simulation. In the characterization of the pore structure, the concept of pore ligament is proposed to relate the pore shape and the critical energy release rate. Two classical fracture problems were used to evaluate the effect of pore shape on the macrofracture behavior. The results of the case study show that although the micropore shapes change the macrofracture behaviors, these effects vary with the geometry and loading conditions of macroscopic boundary value problems. Furthermore, the case study also shows that the influence of micropore structure can be captured at the macroscopic level through the material properties derived from the MD simulations.

42 ENGINEERING↗
DOE OSTI2022

High strain-rate compression behavior of polymeric rod and plate Kelvin lattice structures

The compressive high strain-rate behavior of polymeric Kelvin lattice structures with rod-based or plate-based unit cells was investigated through experimental techniques and finite element simulations. Polymeric lattice structures with 5x5x5 unit cell geometries were manufactured on the millimeter scale using vat polymerization additive manufacturing and tested at low (0.001/s) and high (1000/s) strain-rates. High strain-rate experiments were performed and validated for a viscoelastic split-Hopkinson (Kolsky) pressure bar system (SHPB) coupled with high-speed imaging and digital image correlation (DIC). Experimental results at both low and high strain-rates show the formation of a localized deformation band which was more prevalent in low relative density specimens and low strain-rate experiments. Strain-rate effects of lattice specimens strongly correlate with effects of the base polymer material; both bulk polymer and lattice specimen demonstrated strain-rate hardening, strain-rate stiffening, and decreased fracture strain under dynamic loading. Results show mechanical failure properties and energy absorption depended strongly on the relative density of the lattice specimen and exhibited distinct scaling between relative density and geometry type (rod, plate) and loading rate. High relative density plate-lattices demonstrated inferior mechanical properties to rod-lattices; however, there exists a critical relative density for a given mechanical property (17%- 28%) below which plate-lattices outperform rod-lattices of similar mass. As a result, high strain-rate explicit finite element simulations were performed and showed good agreement with the mechanical failure trends and deformation modes observed in the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shoulder fillet effects in strength distributions of microelectromechanical system components

The failure forces and fracture strengths of polysilicon microelectromechanical system (MEMS) components in the form of stepped tensile bars with shoulder fillets were measured using a sequential failure chain methodology. Approximately 150 specimens for each of four fillet geometries with different stress concentration factors were tested. The resulting failure force and strength distributions of the four geometries were related by a common sidewall flaw population existing within different effective stressed lengths. The failure forces, strengths, and flaw population were well described by a weakest-link based analytical framework. Finite element analysis was used to verify body-force based expressions for the stress concentration factors and to provide insight into the variation of specimen effective length with fillet geometry. Monte Carlo simulations of flaw size and location, based on the strength measurements, were also used to provide insight into fillet shape and size effects. The successful description of the shoulder fillet specimen strengths provides further empirical support for application of the strength and flaw framework in MEMS fabrication and design optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Deep Direct-Use Feasibility Study Reservoir Productivity Uncertainty Analysis for the Tuscarora Sandstone, Morgantown, WV

This dataset contains figures that summarize the Tuscarora Sandstone core permeability data collected from the Preston 119 well in Preston County, WV, and summary results of a stochastic analysis that was used to estimate reservoir productivity for the currently unexplored Tuscarora Sandstone below Morgantown, West Virginia. Uncertainties in reservoir productivity considered the thickness of the reservoir, rock permeability, and fluid viscosity. A Monte Carlo analysis of these uncertain properties was used to predict reservoir flow productivity for the case of a matrix-dominated reservoir, and a fracture-dominated reservoir. These results are summarized in figures, spreadsheets, and maps.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

FY20 Verification of BISON Using Analytic and Manufactured Solutions

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution Modeling and Simulation (M&S) tools for modeling of Light Water Reactors (LWRs). This hub, Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including BISON, CTF, MPACT, and MAMBA. BISON is a fuel performance code which models the thermo-mechanical behavior of nuclear fuel using high performance M&S. It is capable of modeling traditional LWR fuel rods, fuel plates, and TRi-structural ISOtropic (TRISO) fuel particles. It can employ three-dimensional Cartesian, two-dimensional axisymmetric cylindrical, or one-dimensional radial spherical geometry. It includes empirical models for a large variety of fuel physics: temperature- and burnup-dependent thermal properties, fuel swelling and densification, fission gas production, cladding creep, fracture, cladding plasticity, and gap/plenum models. This document details a series of code verification test problems that are used to test BISON. These problems add confidence that the BISON code is a faithful representation of its underlying mathematical model. The suite of verification tests are mapped to the underlying conservation equations solved by the code: heat conduction, mechanics, and species conservation. Twenty-two problems are added for the heat conduction solution, two for the mechanics solution, and none for species conservation. Method of Manufactured Solutions (MMS) capability is demonstrated with three problems, and temperature drops across the fuel gap are tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Size-dependent stochastic tensile properties in additively manufactured 316L stainless steel

Recent work in metal additive manufacturing (AM) suggests that mechanical properties may vary with feature size; however, these studies do not provide a statistically robust description of this phenomenon, nor do they provide a clear causal mechanism. Because of the huge design freedom afforded by 3D printing, AM parts typically contain a range of feature sizes, with particular interest in smaller features, so the size effect must be well understood in order to make informed design decisions. This work investigates the effect of feature size on the stochastic mechanical performance of laser powder bed fusion tensile specimens. A high-throughput tensile testing method was used to characterize the effect of specimen size on strength, elastic modulus and elongation in a statistically meaningful way. The effective yield strength, ultimate tensile strength and modulus decreased strongly with decreasing specimen size: all three properties were reduced by nearly a factor of two as feature dimensions were scaled down from 6.25 mm to 0.4 mm. Hardness and microstructural observations indicate that this size dependence was not due to an intrinsic change in material properties, but instead the effects of surface roughness on the geometry of the specimens. Finite element analysis using explicit representations of surface topography shows the critical role surface features play in creating stress concentrations that trigger deformation and subsequent fracture. The experimental and finite element results provide the tools needed to make corrections in the design process to more accurately predict the performance of AM components.

316L stainless steel↗
DOE OSTI2024

Distributed Strain Sensing based Reservoir Monitoring for Carbon Sequestration Reservoirs

Hydraulic or shear fractures propagation can significantly impact the integrity of the cap rock in carbon sequestration reservoirs. This presentation will evaluate distributed strain sensing (DSS) techniques to detect and provide early warnings for fracture propagation and cap rock integrity issues. We will highlight the advantages of DSS techniques over other methods for measuring cap rock integrity. We will then examine how the geometry of monitoring wells affects detection results and suggest optimal well designs for monitoring. Finally, we will present real-life field examples and propose future research directions.

Jin, Ge↗
DOE OSTI2021

A variational phase-field model For ductile fracture with coalescence dissipation

A novel phase-field model for ductile fracture is presented. The model is developed within a consistent variational framework in the context of finite-deformation kinematics. A novel coalescence dissipation introduces a new coupling mechanism between plasticity and fracture by degrading the fracture toughness as the equivalent plastic strain increases. The proposed model is compared with a recent alternative where plasticity and fracture are strongly coupled. Several representative numerical examples motivate specific modeling choices. In particular, a linear crack geometric function provides an “unperturbed” ductile response prior to crack initiation, and Lorentz-type degradation functions ensure that the critical fracture strength remains independent of the phase-field regularization length. In addition, the response of the model is demonstrated to converge with a vanishing phase-field regularization length. Additionally, the model is then applied to calibrate and simulate a three-point bending experiment of an aluminum alloy specimen with a complex geometry. The effect of the proposed coalescence dissipation coupling on simulations of the experiment is first investigated in a two-dimensional plane strain setting. The calibrated model is then applied to a three-dimensional calculation, where the calculated load-deflection curves and the crack trajectory show excellent agreement with experimental observations. Finally, the model is applied to simulate crack nucleation and growth in a specimen from a recent Sandia Fracture Challenge.

42 ENGINEERING↗
DOE OSTI2023

Predictive machine learning approaches for the microstructural behavior of multiphase zirconium alloys

Abstract Zirconium alloys are widely used in harsh environments characterized by high temperatures, corrosivity, and radiation exposure. These alloys, which have a hexagonal closed packed (h.c.p.) structure thermo-mechanically degrade, when exposed to severe operating environments due to hydride formation. These hydrides have a different crystalline structure, than the matrix, which results in a multiphase alloy. To accurately model these materials at the relevant physical scale, it is necessary to fully characterize them based on a microstructural fingerprint, which is defined here as a combination of features that include hydride geometry, parent and hydride texture and crystalline structure of these multiphase alloys. Hence, this investigation will develop a reduced order modeling approach, where this microstructural fingerprint is used to predict critical fracture stress levels that are physically consistent with microstructural deformation and fracture modes. Machine Learning (ML) methodologies based on Gaussian Process Regression, random forests, and multilayer perceptrons (MLP) were used to predict material fracture critical stress states. MLPs, or neural networks, had the highest accuracy on held-out test sets across three predetermined strain levels of interest. Hydride orientation, grain orientation or texture, and hydride volume fraction had the greatest effect on critical fracture stress levels and had partial dependencies that were highly significant, and in comparison hydride length and hydride spacing have less effects on fracture stresses. Furthermore, these models were also used accurately predicted material response to nominal applied strains as a function of the microstructural fingerprint.

36 MATERIALS SCIENCE↗
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↗