Search NASA⌕ Search

SEARCH · Search NASA

Results for “source term”

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 505 records · Page 28

A Unified SHOC-MF Boundary Layer Parameterization in the NASA Goddard Earth Observing System Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC; Bogenschutz and Krueger, 2013). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies. Results indicate that many higher order moments are underestimated when the mass flux contribution is neglected.

Nathan Arnold↗

Coupling a Computational Fluid Dynamics Model to a Spacecraft Thermal System Model for the DraMS Instrument Thermal Analysis

The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop® (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.

Heat transfer↗

Towards Aerodynamic Shape Optimization Using an Immersed Boundary Overset Grid Method

Traditional Reynolds-averaged Navier-Stokes grid methods applied to aerodynamic shapeoptimization can struggle with the deformation of surface and volume grids at componentintersections, such as at wing-fuselage junctions. To overcome this, we propose an approachwhich utilizes curvilinear overset grids for the discretization of the domain, with the presenceof the body modeled using an immersed boundary method. This approach handles complexgeometries without the need for their explicit integration into the grid. The goal of this approachis to reduce grid generation time and allow for greater geometric freedom for component-basedaerodynamic shape optimization. Two different methods are presented: a source-term-basedand a ghost-node-based immersed boundary method. Flow analyses and adjoint solutionsobtained using the proposed methods show promising comparisons with standard body-fittedgrid methods. Preliminary aerodynamic shape optimization results obtained using one of theimmersed boundary methods are also presented.

TTT↗

Coupling a Computational Fluid Dynamics (CFD) Model to a Spacecraft Thermal System Model for the DraMS Instrument Thermal Analysis

The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.

heat transfer↗

Heated Supersonic Axisymmetric Jet Cases for the NASA Turbulence Modeling Resource

Three new supersonic jet cases were built for the NASA Turbulence Modeling Resource (TMR). These cases were taken from the 6th AIAA Propulsion Aerodynamics Workshop (PAW) nozzle test problem, which utilized data taken at NASA GRC consisting of non-intrusive measurements of jet plume velocities and temperatures, including both mean values and turbulence statistics. The three jet cases all used a Mach 1.63 axisymmetric nozzle: a heated jet for each at on-design conditions, a temperature-matched on-design condition at Mach 1.63, and an off-design heated condition for the Mach 1.63 nozzle. Computational fluid dynamics (CFD) solutions are obtained using three established Reynolds-averaged Navier-Stokes (RANS) codes: Wind-US, FUN3D, and VULCAN-CFD, all using the Menter Shear Stress Transport k- turbulence model with vorticity source term (SST-V). The NASA TMR requires at least two CFD codes yielding essentially identical results to certify that the approaches are verified, meaning that they solve the posed CFD equations as intended. The three codes generally provided very close agreement with each other for jet plume quantities, with the only exception being static temperature. This discrepancy was determined to be the result of whether or not turbulent kinetic energy was considered in the definition of total energy. Validation of the SST-V turbulence model for this case (its level of agreement/disagreement with experimental data) is also addressed.

turbulance↗

Numerical Investigation of Combustor-Turbine Interactions with a Two-Stage High-Pressure Turbine Including Cooling Airflows

Combustor–turbine interactions (CTIs) are investigated by performing three-dimensional unsteady simulations using a realistic combustor and high-pressure turbine (HPT) configurations from the Energy Efficient Engine (E3) program. To understand CTIs, we compare the numerically predicted flow fields from single-component simulations (Step 1: the combustor + the first-stage stator of turbine; Step 2: the two-stage HPT imposing the time-averaged flow solution from Step 1 as the inflow boundary condition) and a fully coupled combustor–turbine simulation (Step 3) at the sea-level takeoff (SLTO) condition. In addition to three previous simulations where the cooling airflows inside the HPT had been neglected for all Step 1, Step 2, and Step 3, two new simulations of Step 2 and Step 3 take into account the cooling airflows using the source team approach. In this approach, to mimic the cooling airflows, we impose the source term at a specific area of a cooling airflow hole at the solid surfaces without making a mesh of each hole. The objective of this study is twofold. One is to investigate the effect of the presence of the cooling airflows on the aerodynamics of the combustor and HPT as well as the HPT performance. The second is to perform a detailed comparison among the calculated flowfields by two different numerical schemes, the central-difference with the standard Jameson–Schmidt–Turkel (CD-JST) scheme and the AUSM scheme. There is a noticeable difference in the hot-streak distributions at the first-stage stator. In addition, depending on the choice of the numerical scheme and the presence of the cooling airflows, an occurrence of shock waves appearing at the first-stage stator is greatly influenced. Thus, this has a noticeable impact on the HPT performance. It is shown that the estimated turbine efficiencies from Step 3 are about 7 % less than the ones from Step 2.

Combustor-turbine interaction↗

Wall-Modeled LES of a Swept Wing with Leading-Edge Ice Using LAVA Curvilinear, Unstructured, and Cartesian Solvers

Wall-modeled large-eddy simulation (LES) of a swept wing with leading-edge ice build up is performed using three mesh paradigms and associated flow solvers and compared to experimental results. The study focuses on an 8.9% scale model of the CRM65 swept wing featuring both high-fidelity and smooth leading-edge ice shapes. The assessment is conducted using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework using the three actively supported meshing paradigms in LAVA: structured curvilinear overset, unstructured, andCartesian. For the iced configurations, these paradigms use hybrid body-fitted/source-term immersed-boundary, body-fitted, and ghost-cell immersed-boundary strategies, respectively.The unstructured and Cartesian mesh paradigms are particularly attractive for complex ice shapes as they avoid the manual mesh generation effort associated with the curvilinear approach. For the high-fidelity ice shape, good agreement with the experiment can be obtained with all three strategies; however, the curvilinear method is particularly sensitive to the source-term immersed-boundary timescale and span wise mesh resolution, and the Cartesian implementation is sensitive to the choice of numerical flux. For the smooth ice, larger discrepancies are observed across all methods. However, a mesh refinement study guided by flow visualizations leads to an improved comparison with the experiment that is particularly pronounced for the unstructured mesh paradigm.

TTT↗

The Contaminant Footprint of Landed Spacecraft: Toward an Inventory and Modelling Framework

All spacecraft generate and carry contaminants, i.e., unwanted and potentially harmful material. When a spacecraft lands and operates in vacuum, as onto Earth’s Moon, it introduces contaminants into its environment that may compromise mission science objectives and engineering performance. Contamination may degrade sites of unique value to planetary science or in situ resource utilization. This presentation will identify and compare source terms and transport vectors for contaminants – in particular, organic material – generated by landed spacecraft. An integrated modeling framework for the organic contamination footprint of spacecraft missions will be described and presented.

Gas Dynamics↗

Organic Inventory - Planetary Protection on the Moon

All spacecraft generate and carry contaminants, i.e., unwanted and potentially harmful material. When a spacecraft lands and operates in vacuum, as onto Earth’s Moon, it introduces contaminants into its environment that may compromise mission science objectives and engineering performance. Contamination may degrade sites of unique value to planetary science or in situ resource utilization. This presentation will identify and compare source terms and transport vectors for contaminants – in particular, organic material – generated by landed spacecraft. An integrated modeling framework for the organic contamination footprint of spacecraft missions will be described and presented.

Planetary Protection↗

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

Fuel Performance Modeling Status Update and Potential Model Improvements

Fuel performance modeling status update and potential model improvements overview of TRISO fuel performance modeling codes PARFUME/BISON, AGR experiment support, potential modeling improvements, BISON smeared cracking model, BISON fission product source term, and AGR-3/4 reirradiation heating test.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore

Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.

hydrogen plasma↗

Fuel Property Effects on Stochastic Preignition Events During Engine Load Transitions

Stochastic preignition (SPI) is an abnormal combustion phenomenon that can cause catastrophic engine damage. There have been several proposed mechanisms of SPI, where a uniform source is still not certain, however, SPI tendencies have been shown to be influenced by engine operating conditions, oil composition, engine age, and fuel chemical and physical properties. Laboratory research and testing for SPI propensity is challenging given the stochastic nature of events, as well as the potential for significant degradation of the engine platform and measuring equipment over time. Thus, SPI specific experiments are generally conducted under either sustained or cyclic patterning of steady-state operating conditions to avoid the influence of transient engine boundary conditions on test parameters of interest (e.g. oil additive package, fuel properties, engine speed/load, etc.). In this work a cyclically varying SPI test sequence involves a 5 min engine warmup period at a low engine load of around 4 bar gross indicated mean effective pressure (IMEPg), followed by a transition to high load (~20 bar IMEPg) at a constant 2000 rev/min engine speed for a total of 25 min. This individual test sequence load schedule is then sequentially repeated 10 times to generate significant statistical data for analysis. This work examines the influence of fuel chemical and physical properties on SPI tendency during the unsteady portion of the 10-cycle sequence (the first 5 min of the high load operation in each sequence of the loading cycle) which has been discarded from previous analyses due to the uncertainty in engine operating and thermal boundary conditions. Results from this analysis suggest an increasing trend in the ratio of SPI events during the unsteady test period relative to the steady test period with increasing fuel Reid Vapor Pressure (RVP), implying differences in uncontrolled ignition source terms, possibly from, fuel wall interactions and retention during the load transition phase of the test.

Splitter, Derek [ORNL] (ORCID:0000000174044047)↗

Meso-scale modeling of UO2 nuclear fuel to high burnup

To improve the economics of light water reactors for commercial nuclear energy generation, utility operators are seeking to obtain regulatory approval to run UO2 fuel to higher levels of burnup. One potential impediment to obtaining this approval is the phenomenon of fuel fragmentation, relocation, and dispersal (FFRD). FFRD can result when fuel experiences a rapid temperature transient, such as that occurring during a Loss Of Coolant Accident (LOCA). FFRD has historically been most associated with the rim region in UO2 fuel pellets, where the phenomenon of fragmentation is also referred to as pulverization due to the small size of the fragments. More recent evidence suggests that the so-called “dark zone” (due to its appearance in micrographs) that can be observed in the mid-radial regions of high burnup fuel is also susceptible to FFRD. Although empirical fuel performance models have been developed that can adequately predict pulverization in the rim region under typical LWR conditions, a scientific understanding of what underlies fuel restructuring and subsequent FFRD is lacking even in the rim region, and no models are currently available for the behavior the dark zone. To address these challenges, the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has employed a multi-scale modeling approach to improve scientific understanding and develop new fuel performance models. In this talk, I will focus on meso-scale efforts, which form a crucial link between atomic-scale and engineering-scale models. Phase-field modeling combined with cluster dynamics is used to predict the restructuring process in the rim region. Phase-field fracture modeling, informed by atomistic simulations, is used to predict the onset of pulverization in the rim region. Combining these techniques together allows the extent of rim pulverization to be predicted. The formation and evolution of the dark zone has also been simulated with the phase-field method, using an improved approach to vacancy source term parameterization. The work shows the important impact of microstructure on fuel performance.

fracture↗

MARIE: A Python-Based Framework for Comprehensive Fuel Recycling Modeling

One of the most pressing challenges to the continued deployment of nuclear energy systems is in the ultimate management and disposition of discharged fuel assemblies. While reprocessing and recovery of valuable materials from UNF assemblies has been considered as part of an overall strategy for minimization of the volume of reactor-based wastes to be managed, the deployment of commercial-scale reprocessing facilities presents an enormous economic challenge. The MARIE software package has been developed as a means of confronting this challenge. Representing components of a generic fuel reprocessing operation as individual physical processes, MARIE is designed as a modular framework intended to allow for analysis and cost-optimization for a hypothetical reprocessing facility while realistically accounting for the physical characteristics of the used fuel source term, such as decay heat, activity, and radiation dose (informing corresponding shielding requirements). Capabilities supported by MARIE include head-end operations such as fuel shearing, voloxidation, and dissolution; generic solvent extraction operations informed by available open-literature data; a suite of unit operations intended to represent electrochemical processing of used fuel assemblies (i.e., oxide reduction, electrorefining, and electrowinning); and finally, accounting for both costs and physical features of discharged waste streams, which can be used to inform follow-on analyses such as the feasibility of deep-borehole disposal of HLW. This paper presents an overview of the MARIE software capabilities, including how individual unit operations are implemented to enable a larger-scale optimization of a hypothetical reprocessing operation on aspects such as cost and recovery of valuable materials.

Skutnik, Steve [ORNL] (ORCID:000000016441135X)↗

Reirradiation and Post-reirradiation Heating Tests of Loose AGR-2 TRISO Particles and AGR-3/4 Fuel Compacts

To address the lack of data on the behavior of short-lived fission products I-131 (t1/2 = 8.02 d) and Xe-133 (t1/2 = 5.24 d) in the TRISO fuel system, loose AGR-2 particles and kernels were reirradiated in NRAD and subjected to post-reirradiation heating in the FACS furnace. This provided data on the release of fission products from exposed kernels, and the results were used to establish FACS condensation plate collection efficiencies for temperatures of 1000, 1200, and 1400°C where previously, the only collection efficiencies were for 1600°C tests. These new collection efficiencies were then applied to the NRAD-FACS tests of entire AGR-3/4 fuel compacts. Understanding the relative behavior of short-lived I-131, for example, is important because it is a significant contributor to offsite dose during postulated accidents, and Xe-133 has previously been used as a conservative indicator of I-131 behavior. These data will be used to support refinement of fission product transport models and HTGR source-term analyses.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessing High Burnup U-19Pu-10Zr Fuel Performance against Historical and Modeled Behavior

Advancing the deployment of sodium-cooled fast reactors (SFRs) requires thorough testing of metallic fuel pins under accident conditions to establish safe operational limits of high burnup fuel. To conduct transient testing, a comprehensive understanding of steady-state fuel behavior obtained through both experimental characterization and accurate predictive capabilities is needed. This study comparatively assesses the steady-state irradiation performance of two high burnup U-19Pu-10Zr fuel pins, DP-36 and DP-40, irradiated under prototypic fast reactor conditions in preparation for planned safety testing at the Transient Reactor Test Facility. Since DP-40 was designated for use in the test and DP-36 serves as its sibling pin, non-destructive, engineering-scale post-irradiation examinations (PIE) were conducted on both pins while destructive examinations were performed exclusively on DP-36. The results were then assessed against historical performance data from similar fuel pins irradiated in the Experimental Breeder Reactor-II. Additionally, the steady-state irradiation of each pin was modeled using the BISON fuel performance code to assess the accuracy of current modeling capabilities in predicting the baseline irradiation behavior. Non-destructive examinations included neutron radiography to measure fuel column elongation, gamma scanning to verify pin integrity and fission product migration, and profilometry to assess dimensional changes. Benchmarking against existing PIE data revealed consistent patterns in axial fuel column growth and cladding diametral strain, though both pins exhibited longer low-density “fluff” structures, which can have implications for core reactivity and source term calculations. Destructive examinations on DP-36 included fission gas release analysis and sectioning for optical microscopy, which showed more complex constituent redistribution patterns than the traditionally accepted 3-ring model. The axial evolution of fractional areas and porosities of each of the redistributed zones were quantified and presented. Modeling comparisons showed agreement in fractional fission gas release but consistently overestimated axial and radial swelling and disagreed with measured axial porosity patterns. These conservative overpredictions suggested that the pins would appear closer to failure or operational limits at the start of transient tests, potentially leading to higher strain accumulation during the transient. While conservative estimates provide safety margins, they can negatively impact fuel economics. A review of the swelling models identified areas for improvement in the gaseous swelling, solid swelling, and fuel hot-pressing models when applied to ternary fuel. The results of this study highlight the critical importance of conducting pre-test characterization on both test and sibling pins to accurately capture steady-state fuel behavior, providing a precise baseline for post-test evaluations and essential inputs for transient modeling of the planned experiments. The analysis also revealed significant data gaps that require further investigation to enhance the understanding and prediction of fuel swelling and pore dynamics. Collecting comprehensive data across different irradiation conditions, burnup levels, and fuel compositions are essential for refining existing models and developing mechanistic models for both binary and ternary metallic fuels, ultimately improving the integration of modeling and experimental approaches in accident testing.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Performance Evaluation of Different Parallel Programming Models in SCALE-Shift Sequences for Criticality and Shielding Applications [Abstract]

The SCALE code system has been widely used for nuclear criticality safety, reactor physics, radiation shielding, source term generation, and inventory analyses by researchers, industry, and regulatory bodies. Although limited support for shared- and distributed-memory parallel processing was introduced via C++ threading, OpenMP, and MPI, a hybrid parallel programming model with both distributed- and shared-memory parallelism has not been fully supported in the SCALE code system.

Nuclear Criticality Safety Program (NCSP)↗