Engineering Papers⌕ Search

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.

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

At least 73 records · Page 4

Transition of SRT Source Term Software to Modern Version Control System to Enable More Efficient Verification, Validation, and User Engagement

The Simplified Radionuclide Transport (SRT) code assists in reactor design and licensing by providing a detailed analysis of radionuclide transport under specific transient scenarios. This analysis is essential for ensuring that reactors meet safety and regulatory standards before they are authorized to operate. Mechanistic source term analysis, which involves a realistic evaluation of radionuclide behavior from the source to the environment, is a key component of this process. It provides a comprehensive understanding of potential release scenarios, supporting the development of robust safety measures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Unclassified Source Term and Radionuclide Data for the Groundwater Flow and Contaminant Transport Model for Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nye County, Nevada, Revision 1

This report documents the evaluation of information and data available on the unclassified source term and radionuclide (RN) contamination for Central and Western Pahute Mesa (PM): Corrective Action Units (CAUs) 101 and 102. The Central and Western PM CAUs constitute two areas of the Nevada National Security Site (NNSS) used for underground nuclear testing (Figure 1-1). Eighty-five underground nuclear tests were conducted at the PM CAUs between 1965 and 1992 (NNSA/NFO, 2015b). Underground nuclear tests deposit radioactive elements in the subsurface, resulting in groundwater contamination in the vicinity of the underground test locations. At the PM CAUs, the testing resulted in more than 2.68 × 10 7 curies (Ci) of radioactivity (decay-corrected to September 30, 2012) remaining underground (Finnegan et al., 2016). Of the 85 underground nuclear tests conducted at the PM CAUs, three were shallow tests detonated within 200 meters (m) of the surface and referred to as the “near-surface tests”; the remaining 82 were detonated deeper than 200 m below ground surface (bgs) and are referred to as the “deep underground tests” in this document.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

LU-SGS implicit algorithm for three-dimensional incompressible Navier-Stokes equations with source term

A numerical method is developed for solving the incompressible Navier-Stokes equations using the concept of pseudocompressibility. A lower-upper symmetric-Gauss-Seidel implicit scheme is developed for three-dimensional incompressible viscous flow computations. The present algorithm offers additional advantages when solving the flow equations with source terms. Complete vectorizability of the algorithm on oblique planes of sweep in three-dimensions is accomplished in a new flow solver, INS3D-LU code. Spatial differencing is a second-order accurate semi-discrete finite-volume method augmented by a third-order accurate numerical dissipation model which is based on spectral-radii. Comparison of numerical solutions for a curved duct with experimental data shows good agreement. The method is applied to calculate the inducer flow of the Space Shuttle Main Engine turbopump.

Yoon, Seokkwan↗

Advanced Modeling and Simulation to Characterize Advanced Boiling Water Reactor Source Terms to Support a Regulatory Approval Pathway for Right-Sized Emergency Planning Zone. Advanced Boiling Water Reactor (A-BWR) Project Phase 2

The United States Department of Energy (US DOE) is currently supporting the development of various small modular reactor (SMR), microreactor and advanced reactor designs. These reactors have improved safety features as compared to conventional Large Light Water Reactors (LLWRs), which include features that improve potential reduction of radiological source terms in the event of design and beyond-design basis accidents. Specifically, some reactors feature a relatively smaller containment volume with respect to the available surface area for fission product deposition; these include Integrated Pressurized Water Reactors (iPWRs), such as the NuScale SMR and SMR-type Boiling Water Reactors (BWRs) such as General Electric-Hitachi’s BWRX-300 design.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

LES Simulation of Cooling Airflow of High-Pressure Turbine Using the Source Term Approach

In this study, the Open National Combustion code (OpenNCC) is applied to simulate the airflow inside the high-pressure turbine (HPT) of the energy efficient engine (EEE). The main objective of this study is to validate the proposed methodology of a three-dimensional unsteady calculation (LES) of the HPT with a moving mesh capability to capture a relative motion at the stator and rotor interface. In our previous study, the adiabatic wall assumption was made at all the solid surfaces of HPT (i.e., the cooling airflows were not included.) To relax this assumption, we implement the surface source team approach, in which we impose the source term at a specific area of a cooling airflow exit at the solid surfaces by specifying an injection angle, temperature, turbulent intensity, and mass flowrate of each cooling airflow. The beneficial feature of this approach is that it does not require a mesh for each hole. For the validation of the model, the data for EEE model is obtained from the General Electric (GE) test campaign, and the full scale warm-air rig test condition representing a test point close to the integrated core/low spool design condition is considered.

CFD Hot-streaks↗

Numerical Investigation of the Effect of Cooling Airflow inside High-Pressure Turbine using the Source Term Approach

In this study, the Open National Combustion code (OpenNCC) is applied to simulate the airflow inside the high-pressure turbine (HPT) of the energy efficient engine (EEE). The main objective of this study is to compare the characteristics of the HPT performance obtained by using three different approaches (Case 1: the realistic combustion products are considered at the HPT inlet, but the cooling airflows are not included, Case 2: the working fluid is assumed to be a single gas (i.e., air), and the cooling airflows are considered, Case 3: the cooling airflow and the combustion products at the HPT inlet are included). To mimic the cooling airflows at the HPT surface, we use the surface source team approach, in which we impose the source term at a specific area of a cooling airflow hole at the solid surfaces by specifying an injection angle, temperature, turbulent intensity, and mass flowrate of each cooling airflow. For the validation of the model, the data for the EEE model was obtained from the General Electric (GE) test campaign (RDG10), and the full scale warm-air rig test condition representing a test point close to the integrated core/low spool design condition was considered. This validation test was done using the stationary uniform inflow condition (i.e., the reported total temperature and total pressure). Secondly, we considered the Sea-Level Take-off (SLTO) engine condition and separately performed the combustor simulation to obtain a more realistic time-averaged, spatially-nonuniform HPT inflow condition. This work provides a better understanding the effect of the cooling airflow model with different assumptions on the HPT performance. It is found that using different working gas have a minor effect on the performance metric, however the presence of the cooling airflows significantly affects the HPT performance such as the turbine efficiency.

CFD Hot-streaks↗

Application of MELCOR for Simulating Molten Salt Reactor Accident Source Terms

Molten Salt Reactor (MSR) systems can be divided into two basic categories: liquid-fueled MSRs in which the fuel is dissolved in the salt, and solid-fueled systems such as the Fluoride-salt-cooled High-temperature Reactor (FHR). The molten salt provides an impediment to fission product release as actinides and many fission products are soluble in molten salt. Nonetheless, under accident conditions, some radionuclides may escape the salt by vaporization and aerosol formation, which may lead to release into the environment. Here, we present recent enhancements to MELCOR to represent the transport of radionuclides in the salt and releases from the salt. Some soluble but volatile radionuclides may vaporize and subsequently condense to aerosol. Insoluble fission products can deposit on structures. Thermochimica, an open-source Gibbs Energy Minimization (GEM) code, has been integrated into MELCOR. With the appropriate thermochemical database, Thermochimica provides the solubility and vapor pressure of species as a function of temperature, pressure, and composition, which are needed to characterize the vaporization rate and the state of the salt with fission products. Since thermochemical databases are still under active development for molten salt systems, thermodynamic data for fission product solubility and vapor pressure may be user specified. This enables preliminary assessments of fission product transport in molten salt systems. In this paper, we discuss modeling of soluble and insoluble fission product releases in a MSR with Thermochimica incorporated into MELCOR. Separate-effects experiments performed as part of the Molten Salt Reactor Experiment in which radioactive aerosol was released are discussed as needed for determining the source term.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Porting the WAVEWATCH III (v6.07) wave action source terms to GPU

Abstract. Surface gravity waves play a critical role in several processes, including mixing, coastal inundation, and surface fluxes. Despite the growing literature on the importance of ocean surface waves, wind–wave processes have traditionally been excluded from Earth system models (ESMs) due to the high computational costs of running spectral wave models. The development of the Next Generation Ocean Model for the DOE’s (Department of Energy) E3SM (Energy Exascale Earth System Model) Project partly focuses on the inclusion of a wave model, WAVEWATCH III (WW3), into E3SM. WW3, which was originally developed for operational wave forecasting, needs to be computationally less expensive before it can be integrated into ESMs. To accomplish this, we take advantage of heterogeneous architectures at DOE leadership computing facilities and the increasing computing power of general-purpose graphics processing units (GPUs). This paper identifies the wave action source terms, W3SRCEMD, as the most computationally intensive module in WW3 and then accelerates them via GPU. Our experiments on two computing platforms, Kodiak (P100 GPU and Intel(R) Xeon(R) central processing unit, CPU, E5-2695 v4) and Summit (V100 GPU and IBM POWER9 CPU) show respective average speedups of 2× and 4× when mapping one Message Passing Interface (MPI) per GPU. An average speedup of 1.4× was achieved using all 42 CPU cores and 6 GPUs on a Summit node (with 7 MPI ranks per GPU). However, the GPU speedup over the 42 CPU cores remains relatively unchanged (∼ 1.3×) even when using 4 MPI ranks per GPU (24 ranks in total) and 3 MPI ranks per GPU (18 ranks in total). This corresponds to a 35 %–40 % decrease in both simulation time and usage of resources. Due to too many local scalars and arrays in the W3SRCEMD subroutine and the huge WW3 memory requirement, GPU performance is currently limited by the data transfer bandwidth between the CPU and the GPU. Ideally, OpenACC routine directives could be used to further improve performance. However, W3SRCEMD would require significant code refactoring to make this possible. We also discuss how the trade-off between the occupancy, register, and latency affects the GPU performance of WW3.

58 GEOSCIENCES↗

Numerical Dissipation and Wrong Propagation Speed of Discontinuities for Stiff Source Terms

In compressible turbulent combustion/nonequilibrium flows, the constructions of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and (b) obtaining correct propagation speed of discontinuities for stiff reacting terms on coarse grids share one important ingredient - minimization of numerical dissipation while maintaining numerical stability. Here coarse grids means standard mesh density requirement for accurate simulation of typical non-reacting flows. This dual requirement to achieve both numerical stability and accuracy with zero or minimal use of numerical dissipation is most often conflicting for existing schemes that were designed for non-reacting flows. The goal of this paper is to relate numerical dissipations that are inherited in a selected set of high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities as a function of stiffness of the source term and the grid spacing.

Yee, H. C.↗

Accident Progression and Source Term Demonstration Calculations for the Advanced Burner Test Reactor (ABTR)

As advanced reactor designers approach broad deployment of new reactor technologies and the associated applications, a critical milestone on this path includes licensing or authorization of these designs in the current regulatory environment. Commercial licensing and DOE authorization both consist of numerous aspects to not only certify the designed plant will operate as intended, but also related to confirmation of the safe operation of these facilities that comply with specific standards and requirements. To that end, domestic industry and regulators require toolsets and methods which can be used to confirm that advanced reactor designs not only align with expected performance during normal operation, but also a series of transient scenarios. These codes require a satisfactory validation basis, and in the absence of such a basis, code-to-code comparisons are often utilized as a surrogate to confirm model adequacy. To support this need for fast reactor technologies, a series of analyses using SAS4A/SASSYS- 1 and the Simplified Radionuclide Transport (SRT) code have been conducted for the Advanced Burner Test Reactor (ABTR). These analyses are expected to not only provide a baseline for transient and source term analyses of ABTR, but to also provide context for the key physics that determine evolution of specific transient scenarios in a sodium-cooled fast reactor (SFR). These codes are utilized for this work given their adoption in the SFR industry and by key SFR stakeholders.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Radiological Source Term Estimation and Isotopic Identification with Parallel Log Domain Particle Filters

This paper presents a parallel log-domain particle filtering algorithm combined with gamma spectrum unfolding to perform localization, identification, and evaluation of multiple point sources of various isotopes in an environment with attenuating obstacles. The method uses sets of precomputed attenuation kernels that map the attenuation characteristics of the environment. These kernels are specific to the energy level of a photopeak of interest. The spectral measurements are deconvolved into count measurements of each photopeak. These count measurements are fed into a set of parallel particle filters using attenuation kernels computed for that photopeak’s energy level. The individual regularized particle filters perform all likelihood calculations in the logarithmic domain to mitigate the effects of particle degeneracy. The output of each particle filter is combined to estimate which isotopes are present as well as their positions and strengths. The performance of the algorithm is characterized in a lab-scale environment using a mobile robot equipped with a gamma ray spectrometer in the presence of up to three different radioactive isotopes simultaneously. The sources were localized to within 10 cm, and their strengths were estimated within 10% of their true values. Furthermore, the isotopes were all correctly identified, and no spurious sources were reported.

42 ENGINEERING↗

A Source-Term Based Boundary Layer Bleed/Effusion Model for Passive Shock Control

A modeling framework for boundary layer effusion has been developed based on the use of source (or sink) terms instead of the usual practice of specifying bleed directly as a boundary condition. This framework allows the surface boundary condition (i.e. isothermal wall, adiabatic wall, slip wall, etc.) to remain unaltered in the presence of bleed. This approach also lends itself to easily permit the addition of empirical models for second order effects that are not easily accounted for by simply defining effective transpiration values. Two effusion models formulated for supersonic flows have been implemented into this framework; the Doerffer/Bohning law and the Slater formulation. These models were applied to unit problems that contain key aspects of the flow physics applicable to bleed systems designed for hypersonic air-breathing propulsion systems. The ability of each model to predict bulk bleed properties was assessed, as well as the response of the boundary layer as it passes through and downstream of a porous bleed system. The model assessment was performed with and without the presence of shock waves. Three-dimensional CFD simulations that included the geometric details of the porous plate bleed systems were also carried out to supplement the experimental data, and provide additional insights into the bleed flow physics. Overall, both bleed formulations fared well for the tests performed in this study. However, the sample of test problems considered in this effort was not large enough to permit a comprehensive validation of the models.

Baurle, Robert A.↗

Source term estimation using noble gas and aerosol samples

Algorithms that estimate the location, time, and magnitude of a point-source atmospheric release using remotely sampled air concentrations typically use data for a single chemical or radioactive isotope. Here, a Bayesian algorithm is presented that uses data from multiple radioactive isotopes that are all released in the same short-duration event. Data from noble gas and aerosol samplers can be used simultaneously in the model. Application to a large synthetic data set using four isotopes shows the new algorithm generally gives more accurate location and time estimates than a comparable model using a single isotope.

54 ENVIRONMENTAL SCIENCES↗

Dirty bomb source term characterization and downwind dispersion: Review of experimental evidence

Dirty bombs are considered one of the easiest forms of radiological terrorism, a form of terrorism based on the deliberate use of radiological material to cause adverse effects in a target population. One U.S. Government official has even described a dirty bomb attack as “all but inevitable”. While people in the vicinity of the blast may experience acute radiation effects, people downwind may unknowingly be contaminated by the radioactive airborne particulate and face increased long-term cancer risk. The likelihood of increased cancer risk depends on the radionuclide used and its specific activity, its aerosolization potential, the particle sizes generated in the blast, and where a person is with respect to the detonation. Different studies have reported that plausible radionuclides for dirty bomb include 60 Co, 90 Sr, 137 Cs, 192 Ir, 241 Am based on their availability in commercial sources as well as safeguards, the amount needed for adverse health effects, previous mishandling of radionuclides and malicious uses. In order to have increased long-term cancer risk, the radionuclide would have to deposit inside the body by entering the respiratory tract and then possibly migrate to other organs or bones (ground shine is not considered in this paper because areas affected by the event will likely become inaccessible). This implies that the particles will have to be smaller than 10 μm to be inhaled. Experiments involving the detonation of dirty bombs have shown that particles or droplets smaller than 10 μm are generated, independently from the initial radionuclide or its state (e.g., powder, solution). Atmospheric tests have shown that in unobstructed terrain, the radionuclide laden cloud can travel kilometers downwind even for relatively small amounts of explosives. Furthermore, buildings in the path of the cloud can change the dose rate. For instance, in one experiment with a single building, the dose rate was 1–2 orders of magnitude lower behind the obstacle compared to its front face. For people walking around, the amount of particulate deposited on them and inhaled will depend on their path relative to the cloud, resulting in the counterintuitive result that the closer people may actually not be the ones more at risk because they could simply miss the bulk of the cloud in their wandering. In summary, the long-term cancer risk for people caught in a dirty bomb cloud away from the detonation requires considering where and when the people are, which radionuclide was used, and the layout of the obstacles (e.g., buildings, vegetation) in the path of the cloud.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗