SEARCH · Engineering Papers

Results for “SOURCE TERMS”

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

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Implementation of two-phase gas transport into VERA for molten salt reactor analysis

Molten salt reactors (MSRs) are a class of next-generation nuclear reactors that have received recent industrial and research interest. A generalized species transport solver was implemented in the Virtual Environment for Reactor Applications (VERA) computing suite to extend this tool to analyze liquid-fueled MSRs. This core simulator has been extended to model the transport of fission product gases into a collection of circulating gas bubbles with the purpose of removing the gases. This paper presents the governing species transport equation, along with various nuclear source terms. Development of the source term for phase migration is discussed, along with a simplified interfacial area tracking method. Finally, a case study on a simplified MSR loop is presented in which modeling parameters were varied to assess their impact on gas removal. The steady state results show that parameters such as bubble diameter, gas injection rate and mass transfer coefficient have a low to moderate effect on the fraction of xenon in the core region. Removal efficiency has the greatest effect on the fraction in the core region. After the pump bowl, bubble diameter has a minor effect on the fraction of xenon in the gas void. Finally, these results point out that increasing parameters such as mass transfer coefficient, gas injection rate, and removal efficiency drives the xenon into the circulating gas void, while decreasing bubble diameter also drives xenon into the gas void by increasing interfacial area.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Predicting large-scale pool fire dynamics using an unsteady flamelet- and large-eddy simulation-based model suite

A low-Mach, unstructured, large-eddy-simulation-based, unsteady flamelet approach with a generalized heat loss combustion methodology (including soot generation and consumption mechanisms) is deployed to support a large-scale, quiescent, 5-m JP-8 pool fire validation study. The quiescent pool fire validation study deploys solution sensitivity procedures, i.e., the effect of mesh and time step refinement on capturing key fire dynamics such as fingering and puffing, as mesh resolutions approach O(1) cm. A novel design-order, discrete-ordinate-method discretization methodology is established by use of an analytical thermal/participating media radiation solution on both low-order hexahedral and tetrahedral mesh topologies in addition to quadratic hexahedral elements. Here, the coupling between heat losses and the flamelet thermochemical state is achieved by augmenting the unsteady flamelet equation set with a heat loss source term. Soot and radiation source terms are determined using flamelet approaches for the full range of heat losses experienced in fire applications including radiative extinction. The proposed modeling and simulation paradigm are validated using pool surface radiative heat flux, maximum centerline temperature location, and puffing frequency data, all of which are predicted within 10% accuracy. Simulations demonstrate that under-resolved meshes predict an overly conservative radiative heat flux magnitude with improved comparisons as compared to a previously deployed hybrid Reynolds-averaged Navier–Stokes/eddy dissipation concept-based methodology.

42 ENGINEERING↗

DOE OSTI2021

As-Run Physics Analysis for the EPRI Zirconium Growth C Capsule

This engineering calculation and analysis report documents the as-run nuclear-heating term, source-term, and radiation damage (measured in displacements per atom [dpa]) accumulation for the Electric Power Research Institute zirconium growth experiment's C capsule (EPRI-ZG-C). These data have been calculated using MCNP 6.1 and SCALE 6.2.3. These data can be used for: the shipping of the EPRI-ZG-C capsule, the acceptance of the capsule at the Materials and Fuels Complex, and programmatic scientific needs.

36 MATERIALS SCIENCE↗

NASA NTRS1993

A controlled variation scheme for convection treatment in pressure-based algorithm

Convection effect and source terms are two primary sources of difficulties in computing turbulent reacting flows typically encountered in propulsion devices. The present work intends to elucidate the individual as well as the collective roles of convection and source terms in the fluid flow equations, and to devise appropriate treatments and implementations to improve our current capability of predicting such flows. A controlled variation scheme (CVS) has been under development in the context of a pressure-based algorithm, which has the characteristics of adaptively regulating the amount of numerical diffusivity, relative to central difference scheme, according to the variation in local flow field. Both the basic concepts and a pragmatic assessment will be presented to highlight the status of this work.

Shyy, Wei↗

NASA NTRS2005

A New Unsteady Model for Dense Cloud Cavitation in Cryogenic Fluids

Contents include the following: Background on thermal effects in cavitation. Physical properties of hydrogen. Multi-phase cavitation with thermal effect. Solution procedure. Cavitation model overview. Cavitation source terms. New cavitation model. Source term for bubble growth. One equation les model. Unsteady ogive simulations: liquid nitrogen. Unsteady incompressible flow in a pipe. Time averaged cavity length for NACA15 flowfield.

Hosangadi, Ashvin↗

DOE OSTI2023

Experiment Design and Preparation for a Shielding Benchmark Utilizing Godiva-IV

An experiment is currently being designed to provide a high-quality shielding benchmark for Criticality Accident Alarm System (CAAS) modeling. Previous benchmarks have suffered from uncertainty primarily due to two factors. The first factor is radiation that is not directly coming from the source, also known as room return. Room return is notoriously difficult to account for in experimental data and to model in benchmarks. The second factor is uncertainty associated with the source term itself. In order to reduce these two sources of uncertainty, this experiment will utilize a room return shield [1] that will reduce the effect of room return and previous work performed to ensure the reproducibility of the source term. The experiment will use the Godiva IV [2] assembly located at the National Criticality Experiments Research Center (NCERC) to provide a neutron source representative of a criticality accident. Previous experiments [3]-[4] have shown that Godiva IV is capable of producing both bursts and steady-state emissions of equal magnitude which will allow for reduced uncertainty in the source term. The room return shield will house and isolate the shielding samples, neutron activation foils, and other detectors. Data collected will provide a comparison of the neutron and gamma dose for shielded and unshielded configurations. This data will be used to validate shielding data used for modeling criticality accident alarm systems. Shielding samples included in this benchmark are C, Fe, NaCl, Pb, high density polyethylene, and SiO 2 . This paper details the design of the room return shield, the experimental plan, challenges, and mitigation strategies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS2003

Efficient Development of High Fidelity Structured Volume Grids for Hypersonic Flow Simulations

A new technique for the control of grid line spacing and intersection angles of a structured volume grid, using elliptic partial differential equations (PDEs) is presented. Existing structured grid generation algorithms make use of source term hybridization to provide control of grid lines, imposing orthogonality implicitly at the boundary and explicitly on the interior of the domain. A bridging function between the two types of grid line control is typically used to blend the different orthogonality formulations. It is shown that utilizing such a bridging function with source term hybridization can result in the excessive use of computational resources and diminishes robustness. A new approach, Anisotropic Lagrange Based Trans-Finite Interpolation (ALBTFI), is offered as a replacement to source term hybridization. The ALBTFI technique captures the essence of the desired grid controls while improving the convergence rate of the elliptic PDEs when compared with source term hybridization. Grid generation on a blunt cone and a Shuttle Orbiter is used to demonstrate and assess the ALBTFI technique, which is shown to be as much as 50% faster, more robust, and produces higher quality grids than source term hybridization.

Alter, Stephen J.↗

DOE OSTI2022

Automated and efficient local adaptive regression for principal component-based reduced-order modeling of turbulent reacting flows

Principal Component Analysis can be used to reduce the cost of Computational Fluid Dynamics simulations of turbulent reacting flows by reducing the dimensionality of the transported variables through projection of the thermochemical state onto a lower-dimensional manifold. However, because of the nonlinearity of the principal component source terms, nonlinear regression techniques must be utilized for the source terms in terms of the principal components. Unfortunately, widely available and utilized nonlinear regression techniques can have prohibitive computational requirements and/or accuracy that is highly dependent on user experience in ad hoc tuning of model architecture and hyperparameters. Here, in this work, a new nonlinear regression approach is proposed that is both computationally efficient and automated so does not require any user input. The approach is evaluated through a priori prediction of principal component source terms using data from a Direct Numerical Simulation of a turbulent nonpremixed n-heptane/air jet flame. In particular, the proposed framework consists of local regressions whose complexity is adapted according to the local nonlinearity of the data: local linear regression when accurate enough and local Artificial Neural Networks when nonlinear regression is required. The number of local clusters for local regression is determined automatically using the Davies-Bouldin index. In addition, Bayesian optimization is utilized for model training (i.e., to select the best architectures and hyperparameters of the nonlinear regressions in an unsupervised fashion), eliminating ad hoc hand-tuning and/or expensive grid searches. Overall, compared to a single, global neural network, the new local adaptive regression approach is shown to have comparable accuracy but 69% less training time due to the utilization of local linear regression and faster training of local neural networks.

42 ENGINEERING↗

DOE OSTI2022

Panel discussion: Challenges and improvements in accident dose analysis; Regulatory and industry perspective

Full text of publication follows. Design basis accident (DBA) analyses are one of several ways nuclear power plant licensees demonstrate, and the U.S. Nuclear Regulatory Commission (NRC) staff verifies, that the design and location of a nuclear power plant are appropriate and that the protection of health and safety of the public is maintained. In July 2000, the NRC staff issued Regulatory Guide 1.183, 'Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors' (RG 1.183) [1]. RG 1.183 provides a method acceptable to the NRC staff for complying with the regulatory requirements contained in Title 10 of the Code of Federal Regulation (10 CFR) Section 50.67 (50.67), 'Accident source term' [2]. RG 1.183 also provides assumptions and parameters used to model postulated DBA radiological analyses for light-water reactors (LWR). Since the initial issuance of RG 1.183, the NRC staff and the nuclear power plant licensees have both gained substantial experience with the implementation of 10 CFR 50.67 and RG 1.183. Based on this experience, comments from licensees, the anticipation of licensing advanced LWRs, and new research, the NRC is developing a proposed revision to RG 1.183. This panel discussion will provide a discussion of that proposed revision, its potential impact on existing and future evaluations, and the industry perspective on the revised guidance. References: 1) USNRC, Regulatory Guide 1.183, 'Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors,' Revision 0, U.S. Nuclear Regulatory Commission, Washington, D.C. (2000), Electronic copies are also available in ADAMS (http://www.nrc.gov/reading-rm/adams.html), under Accession No. ML003716792. 2) U.S. Code of Federal Regulations, Title 10 (Energy), Part 50 (Domestic Licensing of Production and Utilization Facilities), Section 67, (Accident source term), Office of the Federal Register, Washington, DC (1999). (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1988

The exact calculation of quadrupole sources for some incompressible flows

This paper is concerned with the application of the acoustic analogy of Lighthill to the acoustic and aerodynamic problems associated with moving bodies. The Ffowcs Williams-Hawkings equation, which is an interpretation of the acoustic analogy for sound generation by moving bodies, manipulates the source terms into surface and volume sources. Quite often in practice the volume sources, or quadrupoles, are neglected for various reasons. Recently, Farassat, Long and others have attempted to use the FW-H equation with the quadrupole source and neglected to solve for the surface pressure on the body. The purpose of this paper is to examine the contribution of the quadrupole source to the acoustic pressure and body surface pressure for some problems for which the exact solution is known. The inviscid, incompressible, 2-D flow, calculated using the velocity potential, is used to calculate the individual contributions of the various surface and volume source terms in the FW-H equation. The relative importance of each of the sources is then assessed.

Brentner, Kenneth S.↗

DOE OSTI2022

Initial development of a generic fluoride salt-cooled reactor model

Fluoride high-temperature reactors (FHRs) are high-temperature, low-pressure reactor concepts that use tri-structural isotropic (TRISO) fuel and molten fluoride salt coolant. These reactors have the potential to provide both electrical power and high-temperature process heat. We used generic FHR parameters for a pebble-bed FHR to develop an initial model with fresh fuel for a generic FHR (gFHR) in MELCOR and SCALE (NEWT and KENO). In this paper, we present the development of our gFHR models, which will serve as the baseline for a sensitivity and uncertainty analysis to quantify the range of possible source terms for FHRs in severe accidents. We present MELCOR results for fuel and coolant temperatures through the core, a nodalization study for the steady-state thermal hydraulic model, and development of reactor physics models in SCALE. As this work progresses, these models will be used to calculate source terms for a loss-of-forced-flow accident and to conduct a sensitivity study on this accident to establish a range of possible source terms. SCALE will provide reactor physics parameters like isotopic inventory, decay heat generation, and temperature coefficients of reactivity. Using the uncertainty quantification tools within SCALE, we will generate distributions for those parameters and will use the uncertainty quantification code RAVEN or DAKOTA to sample those distributions in MELCOR to quantify the impact of reactor physics and thermal hydraulic uncertainties on FHR source terms. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

The Release of the Simplified Radionuclide Transport (SRT) Code (V.2.1)

Reactor licensing centers on the protection of the public and environment from the inadvertent release of radioactive material. Therefore, mechanistic source term analysis, or the realistic evaluation of radionuclide transport from the source to the environment for specific transient scenarios, is vital to reactor licensing efforts. Developed to resolve a gap in mechanistic source term modeling capabilities for sodium fast reactors (SFRs), the Simplified Radionuclide Transport (SRT) code created by Argonne National Laboratory (Argonne) is now utilized by advanced reactor vendors, universities, and research institutions to support a multitude of design, licensing, and research efforts. Recently, SRT version 2.1 was released, which includes improvements to code models and verification and validation (V&V) suite to support the SRT user community. The following work provides an overview of the improvements made as part of the release of SRT version 2.1. This effort is supported by the U.S. Department of Energy Office of Nuclear Energy (DOE:NE) Advanced Reactor Technologies (ART) Fast Reactor Program (FRP), as part of the program’s support of national laboratory design and safety analysis computer codes utilized by the fast reactor industry. The expansion of SRT code capabilities and improvements to code V&V associated with version 2.1 are in response to user requests and lessons learned from recent source term analyses performed by Argonne and advanced reactor vendors. They also align with the evolving role of SRT, from research and development tool to software utilized for reactor licensing calculations. The report is structured in alignment with the code improvements, as outlined in Figure 1-1. Section 2 provides background information on SRT, including its history, capabilities, and utilization. Section 3 details new code capabilities as part of version 2.1, while Section 4 focuses on the expansion of the code’s V&V suite. Lastly, Section 5 provides a summary and discussion of next steps.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

DOE OSTI2023

Review and Assessment of Available Data Regarding the Behavior of Sodium Aerosols

Recently, there has been a resurgence of interest in advanced (non-light water) reactor designs, including sodium fast reactors (SFRs). In parallel, multiple efforts are underway to develop risk informed, performance-based licensing pathways for advanced reactors. However, such licensing approaches depend on the accuracy of the associated safety assessments, including mechanistic source term (MST) analyses. Unlike historical source term assessments, which utilizes bounding estimates of radionuclide release, MSTs attempt to realistically estimate the transport and retention of radionuclides for specific transient scenarios. For SFR MST analyses, the behavior of radionuclide aerosols is a key factor, as noted in recent studies. Due to the use of sodium as the primary coolant, many potential transient scenarios involve the release of radionuclide aerosols in conjunction with sodium aerosols. For such scenarios, sodium may be the dominant aerosol species, compared to that of released radionuclides. Therefore, understanding the behavior of sodium aerosols is particularly important to an accurate assessment of aerosol transport. Although many approaches to aerosol modeling exist, it is vital to establish their capabilities in evaluating sodium aerosol behavior for SFR analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRSDate not supplied

An Improved Approach to the Predictability & Reliability of the Onset of Turbulence With Shocks

The construction of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and for (b) obtaining correct propagation speed of discontinuities in the presence of stiff source terms share one important ingredient – minimization of numerical dissipation while maintaining numerical stability. The dual requirements to achieve both numerical stability and minimal numerical dissipation are often conflicting since existing shock capturing schemes were designed mainly to be robust for rapidly developed turbulence-free flows and for shock waves without stiff source term. For the past two decades, Yee and collaborators have focused on an improved understanding of the nonlinear behavior of different high order shock-capturing methods. It was found that even very high order methods without proper nonlinear stability and numerical dissipation control can either numerically smear the onset of turbulence due to excess numerical dissipation, or induce (onset) numerical turbulence that is not physical turbulence due to lack of proper numerical dissipation to improve nonlinear stability for long time integration. Our approach is to combine (I) and (II) below for obtaining the physically correct onset of turbulence with shocks, including problems with stiff source terms: (I) Nonlinear dynamics is utilized to complement the traditional linearized stability theory (Yee & Sweby, Yee et al., Griffiths et al., Lafon & Yee, Yee, Wang et al., Kotov et al. 1990- 2015) in order to (i) Minimize numerically induced false transition to turbulence, (ii) Minimize numerical instability due to long time integration of turbulent flows, (iii) Minimize numerically induced standing wave solutions, and (iv) Minimize wrong propagation of speed of discontinuities due to the presence of stiff source terms. (II) Our recently developed physical preserving (structural preserving) high order methods with improved nonlinear stability & accuracy that are essential in minimizing spurious numerics are used.

HECC↗

DOE OSTI2020

A new form of the Saint-Venant equations for variable topography

The solution stability of river models using the one-dimensional (1D) Saint-Venant equations can be easily undermined when source terms in the discrete equations do not satisfy the Lipschitz smoothness condition for partial differential equations. Although instability issues have been previously noted, they are typically treated as model implementation issues rather than as underlying problems associated with the form of the governing equations. This study proposes a new reference slope form of the Saint-Venant equations to ensure smooth slope source terms and eliminate one source of potential numerical oscillations. It is shown that a simple algebraic transformation of channel geometry provides a smooth reference slope while preserving the correct cross-section flow area and the total Piezometric pressure gradient that drives the flow. The reference slope method ensures the slope source term in the governing equations is Lipschitz continuous while maintaining all the underlying complexity of the real-world geometry. The validity of the mathematical concept is demonstrated with the open-source Simulation Program for River Networks (SPRNT) model in a series of artificial test cases and a simulation of a small urban creek. Validation comparisons are made with analytical solutions and the Hydrologic Engineering Center's River Analysis System (HEC-RAS) model. The new method reduces numerical oscillations and instabilities without requiring ad hoc smoothing algorithms.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2018

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.

Sharma, Maitreyee↗

NASA NTRS2018

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

Sharma, Maitreyee↗