Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum modeling”

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 811 records · Page 45

Evidence for elemental variation in the ejecta of the Tycho supernova remnant

We present an X-ray study of the Tycho supernova remnant utilizing archival data from the high-resolution imagers (HRIs) on Einstein and ROSAT, the low-energy imaging telescopes (LEITs) on EXOSAT, and spectral data from the Broad Band X-ray Telescope (BBXRT). We have made use of the differing HRI bandpasses to construct images of Tycho in two spectral bands, o.7-1.8 keV and 1.8-4.5 keV. We find that the two images differ, with the harder image showing enhanced emission along much of the south, west, and north periphery. There appears to be enhanced soft emission in the interior and in one particular knot of emission in the southeast. Besides continuum (which we model here as thermal bremsstrahlung emission), we believe the hard image shows primarily the distribution of high-ionization Si and S K-shell lines which lie in the 1.8-2.6 keV band, while the softer image has contributions from Si as well as Fe XVII to Fe XXIV L-shell lines in the 0.7-1.4 keV band. Guided by the results of nonequilibrium ionization modeling of the BBXRT spectral data, we interpret the observed contrast in hard and soft X-ray emission in terms of variations in abundance, ionization timescale, and temperature. The most likely explanation for the spectral differences are spatial variations of the relative abundances of Si, S, and Fe.

Vancura, O.↗

XTE Proposal #20102--"SS 433's High Energy Spectrum"

We observed the jet-producing compact binary system SS 433 with RXTE during three multiwavelength campaigns, the first in conjunction with ASCA observations, the second simultaneous with a VLA-VLBA-MERLIN campaign, and the third associated with a Nobeyama millimeter-band campaign. All these campaigns included optical observations. Occurring at different jet precession and binary phases, the observations also monitored the system during a radio flare. The data provide SS 433's X-ray spectrum over more than an energy decade, and track the spectral variations as the X-ray source was partially eclipsed. The continuum can be modeled as a power law with an exponential cutoff, which can be detected to approximately 50 keV. Strong line emission is evident in the 5-10 keV range which can be modeled as a broad line whose energy is precession independent and a narrow line whose energy does vary with jet precession phase; this line model is clearly an over simplification since the PCA does not have sufficient energy resolution to detect the lines ASCA observed. The eclipses are deeper at high energy and at jet precession phases when the jets are more inclined towards and away from us. A large radio flare occurred between two sets of X-ray monitoring observations; an X-ray observation at the peak of the flare found a softer spectrum with a flux approximately 1/3 that of the quiescent level.

Band, David L.↗

An Intercomparison of Radiation Codes for Retrieving Upper Tropospheric Humidity in the 6.3-micron Band: A Report from the 1st GVaP Workshop

An intercomparison of radiation codes used in retrieving upper tropospheric humidity (UTH) from observations in the v2 (6.3 microns) water vapor absorption band was performed. This intercomparison is one part of a coordinated effort within the GEWEX Water Vapor Project (GVaP) to assess our ability to monitor the distribution and variations of upper tropospheric moisture from space-borne sensors. A total of 23 different codes, ranging from detailed line-by-line (LBL) models, to coarser resolution narrow-band (NB) models, to highly-parameterized single-band (SB) models participated in the study. Forward calculations were performed using a carefully selected set of temperature and moisture profiles chosen to be representative of a wide range of atmospheric conditions. The LBL model calculations exhibited the greatest consistency with each other, typically agreeing to within 0.5 K in terms of the equivalent blackbody brightness temperature (T(sub b)). The majority of NB and SB models agreed to within +/- 1 K of the LBL models, although a few older models exhibited systematic T(sub b) biases in excess of 2 K. A discussion of the discrepancies between various models, their association with differences in model physics (e.g. continuum absorption), and their implications for UTH retrieval and radiance assimilation is presented.

Soden, B.↗

Concept for Determining the Life of Ceramic Matrix Composites Using Nondestructive Characterization Techniques

Damping measurements have been taken on ceramic matrix composite (CMC) turbopump blisks in the as fabricated, post proof testing, and post turbopump testing conditions. These results indicate that damping is able to quantify fatigue of the CMC blisk. This gives hope for the potential of determining the actual and residual life of CMC materials using a combination of nondestructive techniques. If successful, then this new paradigm for life prediction of CMCs could revolutionize the approach for designing and servicing CMC components, thereby significantly reducing costs for design, development, health monitoring, and maintenance of CMC components and systems. The Nondestructive Characterization (NDC) life prediction approach would complement life prediction using micromechanics and continuum finite element models. This paper reports on the initial concept of NDC life prediction and how changes in damping and ultrasonic elastic modulus data have established the concept as a possibility.

Effinger, Michael↗

Nondestructive Characterization Techniques Used for Ceramic Matrix Composite Life Determination

Recent results indicate that the specific damping capacity and resonant frequency measurements taken periodically during a component's lifetime is able to quantify the mechanical fatigue of CMCS. This gives hope for the potential of determining the actual and residual life of CMC materials using a combination of nondestructive techniques. If successful, then this new paradigm for life prediction of CMCs could revolutionize the approach for designing and servicing CMC components, thereby significantly reducing costs for design, development, health monitoring, and maintenance of CMC components and systems. The Nondestructive Characterization (NDC) life prediction approach would complement life prediction using micromechanics and continuum finite element models. This paper reports on the initial concept of NDC life prediction, a review of the C/SiC blisk damping data, and how changes in the specific damping capacity & ultrasonic elastic modulus data have established the concept as a possibility.

Effinger, Michael↗

Quiescent Giant Molecular Cloud Cores in the Galactic Center

We have used the Long Wavelength Spectrometer (LWS) aboard the Infrared Space Observatory (ISO) to map the far-infrared continuum emission (45-175 micrometer) toward several massive Giant Molecular Cloud (GMC) cores located near the Galactic center. The observed far-infrared and submillimeter spectral energy distributions imply low temperatures (approx. 15 - 22 K) for the bulk of the dust in all the sources, consistent with external heating by the diffuse ISRF and suggest that these GMCs do not harbor high- mass star-formation sites, in spite of their large molecular mass. Observations of FIR atomic fine structure lines of C(sub II) and O(sub I) indicate an ISRF enhancement of approx. 10(exp 3) in the region. Through continuum radiative transfer modeling we show that this radiation field strength is in agreement with the observed FIR and submillimeter spectral energy distributions, assuming primarily external heating of the dust with only limited internal luminosity (approx. 2 x 10(exp 5) solar luminosity). Spectroscopic observations of millimeter-wave transitions of H2CO, CS, and C-34S carried out with the Caltech Submillimeter Observatory (CSO) and the Institut de Radio Astronomie Millimetrique (IRAM) 30-meter telescope indicate a gas temperature of approx. 80 K, significantly higher than the dust temperatures, and density of approx. 1 x 10(exp 5)/cc in GCM0.25 + 0.01, the brightest submillimeter source in the region. We suggest that shocks caused by cloud collisions in the turbulent interstellar medium in the Galactic center region are responsible for heating the molecular gas. This conclusion is supported by the presence of wide-spread emission from molecules such as SiO, SO, and CH3OH, which are considered good shock tracers. We also suggest that the GMCs studied here are representative of the "typical", pre-starforming cloud population in the Galactic center.

Lis, D. C.↗

Modeling Fiber Kinking at the Microscale and Mesoscale

A computational micromechanics (CMM) model is employed to interrogate the assumptions of a recently developed mesoscale continuum damage mechanics (CDM) model for fiber kinking. The CMM model considers an individually discretized three dimensional fiber and surrounding matrix accounting for nonlinearity in the fiber, matrix plasticity, fiber/matrix interface debonding, and geometric nonlinearity. Key parameters of the CMM model were measured through experiments. In particular, a novel experimental technique to characterize the in situ longitudinal compressive strength of carbon fibers through indentation of micropillars is presented. The CDM model is formulated on the basis of Budiansky's fiber kinking theory (FKT) with a constitutive deformation-decomposition approach to alleviate mesh size sensitivity. In contrast to conventional mesoscale CDM models that prescribe a constitutive response directly, the response of the proposed model is an outcome of material nonlinearity and large rotations of the fiber direction following FKT. Comparison of the predictions from the CMM and CDM models shows remarkable correlation in strength, post-peak residual stress, and fiber rotation, with less than 10% difference between the two models in most cases. Additional comparisons are made with several fiber kinking models proposed in the literature to highlight the efficacy of the two models. Finally, the CMM model is exercised in parametric studies to explore opportunities to improve the longitudinal compression strength of a ply through the use of nonconventional microstructures.

Herraez, Miguel↗

NuSTAR Observation of GROJ1744−28 at Low Mass Accretion Rate

Context.Neutron stars in low-mass X-ray binaries (LMXBs) are important systems to study the physics of accretion onto compact ob-jects. The system GRO J1744−28 is particularly interesting, as it usually shows clear pulsations as well as X-ray bursts. Additionally,there are claims for a magnetic field of 5×1011G through the detection of a Cyclotron Resonant Scattering Feature (CRSF).Aims.We present the spectral analysis of GRO J1744−28 using∼29 ks of NuSTAR data taken in 2017 February at a low luminosityof 3.2×1036erg s−1(3–50 keV). Our goal is to study the variability of the source spectrum with pulse phase and to search for the claimed CRSF. Methods.The continuum spectrum is modeled with an absorbed power-law with exponential cut-off, and an additional iron line component. We find no obvious indications for a CRSF and therefore perform a detailed cyclotron line search using statistical methods.We perform this search on pulse phase-averaged as well as phase-resolved spectra. Results.GRO J1744−28 was observed at one of the lowest luminosities studied so far. The previously detected Type II X-ray bursts re absent. Clear pulsations at a period of 2.141124(9) Hz are detected. The pulse profile shows an indication of a secondary peak,which was not seen at higher flux. The 4σupper limit for the strength of a CRSF in the 3–20 keV band is 0.07 keV, lower than the strength of the line found at higher luminosity. Conclusions.The detection of pulsations shows that the source did not enter the “propeller” regime, even though the source flux of4.15×10−10erg cm−2s−1was almost one order of magnitude below the threshold for the propeller regime claimed in previous studies on this source. The transition into the propeller regime in GRO J1744−28 must therefore be below a luminosity of 3.2×1036erg s−1(3–50 keV), which implies a surface magnetic field.2.9×1011G and mass accretion rate.1.7×1016g s−1. A change of the CRSF depth as function of luminosity is not unexpected and has been observed in other sources. This result possibly implies a change in emission geometry as function of mass accretion rate to reduce the depth of the line below our detection limit.

Ole König↗

Implementation of Surface Tension on a Reacting Flow Solver, PeleLM: Preprint

In liquid rocket engines, the fuel is supplied to the combustion chamber in the liquid state though injectors. Such fuel undergoes atomization, vaporization, and combustion processes. To design reliable and efficient injectors, it is required to understand the full processes. This research is part of an effort to develop a full atomization-vaporization-combustion solver from first principles. As an initial step to tackle the atomization process, a multiphase flow solver is under development. For the development, a library of the volume of fluid scheme for multiphase, IRL is coupled with a reacting Navier-Stokes equation solver, PeleLM. Furthermore, as the surface tension has considerable effects on spray breakup. surface tension is implemented in the momentum equation using the continuum surface force model and the improved height function technique.

height function↗

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

Chris Ostrom↗

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

Chris Ostrom↗

Twistronics of Janus transition metal dichalcogenide bilayers

Twisted multilayers of two-dimensional (2D) materials are an increasingly important platform for investigating quantum phases of matter, and in particular, strongly correlated electrons. The moiré pattern introduced by the relative twist between layers creates effective potentials of long wavelength, leading to electron localization. However, in contrast to the abundance of 2D materials, few twisted heterostructures have been studied until now. Here, in this study, we develop a first-principles continuum theory to study the electronic bands introduced by moire patterns of twisted Janus transition metal dichalcogenide (TMD) homo- and heterobilayers. The model includes lattice relaxation, stacking-dependent effective mass, and Rashba spin-orbit coupling. We then perform a high-throughput generation and characterization of DFT-extracted continuum models for more than a hundred possible combinations of materials and stackings. Our model predicts that the moiré physics and emergent symmetries depend on chemical composition, vertical layer orientation, and twist angle, so that the miniband wave functions can form triangular, honeycomb, and kagome networks. Rashba spin-orbit effects, peculiar of these systems, can dominate the moiré bandwidth at small angles. Our work enables detailed investigation of Janus twisted heterostructures, allowing the discovery and control of novel electronic phenomena.

36 MATERIALS SCIENCE↗

Continuum variability of the optically-active radio-quiet quasar GQ Comae

The use of thermal accretion disk models to explain the continuum shape and variability of GQ Comae is outlined. The mass is estimated at 7 x 10 to the 8th solar masses and its accretion rate at 0.1-0.2 solar masses per year. The effects of assuming different values of the inclination angle of the disk, and of the specific angular momentum of the central black hole, are explored. Additional uncertainties in the emergent fluxes in these models are discussed.

Sitko, Michael L.↗

Development of an Open-source Alloy Selection and Lifetime Assessment Tool for Structural Components in CSP

Lack of sufficient data on high temperature mechanical and corrosion behavior of structural materials is a huge barrier in the technological maturity of current and future Concentrating Solar Power (CSP) technologies. Rapid development and selection of materials cannot be achieved by expensive and time-consuming acquisition of experimental data. The goal of the proposed work is development of an open-source alloy selection and lifetime prediction tool that will integrate validated physics-based models to describe influence of temperature, alloy composition, environment and component geometry (thickness) on mechanical and corrosion behavior of Ni and Fe-based alloys employed in molten salts/sCO 2 heat exchangers. This one-year project leveraged the extensive dataset on the creep\corrosion behavior of candidate materials generated at ORNL through past projects and input from current collaborations with industrial partners. Based on previous experience and the feedback provided by industry (Brayton Energy and Echogen), three candidate materials of interest, Ni-based alloys 740H, 282 and 625 and application-specific operating conditions (max. temperature of 730 °C and stress of 150 MPa) were identified for the heat exchanger. An extensive corrosion and creep dataset was assimilated for the relevant operating conditions and was supported by detailed characterization of about 100 metallographic cross-sections. The corrosion dataset consisted of scanning electron microscopy images (secondary electron and backscatter electron), measured concentration profiles of alloying elements using energy dispersive X-ray spectroscopy (EDS), widths of denuded zones (dissolution of strengthening phases) and depths of attack in molten KCl-MgCl 2 mixtures using image analyses. The creep dataset comprised of creep rupture data and creep strain curves (for 740H and 282). Coupled thermodynamic-kinetic microstructure-based models were employed to predict the stress-corrosion induced compositional and phase evolutions in the alloy during operation under the identified operating conditions. Reduced order models were developed from advanced physics-based models and were integrated in a user-friendly alloy selection tool. The corrosion model was able to predict the time to a critical Cr concentration at the oxide/alloy interface (chemical lifetime) within ±10% (1 standard deviation) of typical statistical variation in corrosion tests and EDS measurement errors (±0.5 wt%). The initial scope of the project was limited to predict creep rupture times (Larson-Miller parameter). Based on the input provided by industry, the mechanical lifetime of the heat exchanger is governed by accumulated creep strains (2%) rather than creep rupture. To be able to predict the times to specific creep strains, a more extensive creep model development was undertaken largely beyond the initial scope of the project. The continuum damage mechanics creep model was able to predict times to 2% creep strain, t 2% with an accuracy of ±500h. Ultimately, a screening protocol for SiC was generated to demonstrate the pathway for integration of one of the currently immature materials from a commercial adoption standpoint in the current material evaluation tool. The modeling tool developed here is accessible to the science community and stakeholders and lays the foundation for methods that will enable a rapid evaluation of optimum materials for CSP applications and reliable prediction of material degradation thereby considerably reducing operational costs, improving reliability and increasing overhaul intervals. However, the complete potential of such a tool to include a wider range of materials and test conditions can only be realized with a more concentrated combined experimental-characterization-computation effort.

14 SOLAR ENERGY↗

Inclusion of discrete-to-continuum coupling in multiphoton excitation and dissociation calculations

A nonperturbative approach to multiphoton excitation and dissociation of molecules is presented in which coupling to the continuum is treated explicitly. Transitions among continuum levels are not modeled directly, but something of their effect is represented by assuming that the continuum population density is so low as to be effectively zero at all times. Two trial applications are briefly discussed.

Zimmerman, I. H.↗

Uncovering the Spectral Energy Distribution in Active Galaxies Using High Ionization Mid-Infrared Emission Lines

The shape of the spectral energy distribution of active galaxies in the EUV soft X-ray band (13.6 eV to 1 keV) is uncertain because obscuration by dust and gas can hamper our view of the continuum. To investigate the shape of the spectral energy distribution in this energy band, we have generated a set of photoionization models which reproduce the small dispersion found in correlations between high-ionization mid-infrared emission lines in a sample of hard X-ray selected AGN. Our calculations show that a broken power-law continuum model is sufficient to reproduce the [Ne V]14.32 microns/[Ne III], [Ne V]24.32 microns/[O IV]25.89 micron and [O IV] 25.89 microns/[Ne III] ratios, and does not require the addition of a "big bump" EUV model component. We constrain the EUV-soft X-ray slope, alpha(sub i), to be between 1.5 - 2.0 and derive a best fit of alpha(sub i) approx. 1.9 for Seyfert 1 galaxies, consistent with previous studies of intermediate redshift quasars. If we assume a blue bump model, most sources in our sample have derived temperatures between T(sub BB) = 10(exp 5.18) K to 10(exp 5.7) K, suggesting that the peak of this component spans a large range of energies extending from approx. (Lambda)600 A to > (Lambda)1900 A. In this case, the best fitting peak energy that matches the mid-infrared line ratios of Seyfert 1 galaxies occurs between approx. (Lambda)700-(Lambda)1000 A. Despite the fact that our results do not rule out the presence of an EUV bump, we conclude that our power-law model produces enough photons with energies > 4 Ry to generate the observed amount of mid-infrared emission in our sample of BAT AGN.

Melendez, M.↗

Impact of the continuum on the γ decay of the lowest 2 + states in 14 C

The half-life, T 1/2 =14.6(33) fs, of the 7012-keV 2$^{+}_{1}$ state in 14 C was measured in an experiment employing the 9 Be( 6 Li,pγ) fusion-evaporation reaction and the GODDESS setup. The γ decay from the 2$^{+}_{2}$ near-threshold resonance, located 142 keV above the 8176-keV neutron-emission threshold, was also investigated. An upper limit of 4.0 × 10 –5 was established for the γ-decay branching ratio with respect to the neutron-decay channel. The B(E2) transition probabilities for these 2$^{+}_{1}$ and 2$^{+}_{2}$ states were compared to predictions from the Shell Model Embedded in the Continuum (SMEC). Significant modifications of these B(E2) probabilities, with respect to the standard shell model, are expected due to the coupling to the continuum. For calculations using the YSOX interaction, agreement was found for large negative values of V 0 , the coupling constant to the continuum. The central value V 0 = −645 MeV fm 3 results in a retardation by a factor ∼ 3.5 for the 2$^{+}_{1}$ → 0$^{+}_{1}$ transition, and an enhanced probability by a factor ∼ 2.5 for the γ-decay out of the 2$^{+}_{2}$ state. The latter factor reflects the effect of collectivization of the 2$^{+}_{2}$ excitation when the coupling to the continuum is taken into account.

14C↗

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↗