Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source term calculation”

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 271 records · Page 15

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗

An assumed pdf approach for the calculation of supersonic mixing layers

In an effort to predict the effect that turbulent mixing has on the extent of combustion, a one-equation turbulence model is added to an existing Navier-Stokes solver with finite-rate chemistry. To average the chemical-source terms appearing in the species-continuity equations, an assumed pdf approach is also used. This code was used to analyze the mixing and combustion caused by the mixing layer formed by supersonic coaxial H2-air streams. The chemistry model employed allows for the formation of H2O2 and HO2. Comparisons are made with recent measurements using laser Raman diagnostics. Comparisons include temperature and its rms, and concentrations of H2, O2, N2, H2O, and OH. In general, good agreement with experiment was noted.

Baurle, R. A.↗

Thermal Radiation Analyzer

Thermal Radiation Analyzer System TRASYS II solves radiation related aspects of thermal analysis problems. TRASYS II provides for calculation of internode radiation interchange data and for calculation of incident and absorbed heat rate data originating from environmental radiant heat sources. When TRASYS II used in conjunction with generalized thermal analysis program such as Systems Improved Numerical Differencing Analyzer (SINDA), any thermal problem expressed in terms of lumped parameter RC thermal network solved.

Vogt, R. A.↗

Structural, Criticality, and Radiation Dose Calculations to support SNF Loading into a DOE Standard Canister

The DOE Standard Canister Demonstration Project includes the development of an internal support structure (ISS) for the 4.6-m long and 45.7-cm diameter canister. This study presents structural evaluations, criticality safety assessments, and dose rate calculations conducted within the scope of the ISS design process to support smooth canister loading operations and to ensure safe storage, transportation, and disposal of Peach Bottom 1 Core II (PB2) and Fort St. Vrain (FSV) SNF, currently stored at the CPP-603 facility of the INL. The ISS includes a 316L stainless steel basket that holds twelve PB2 SNF rods. Further, six individual, 78.7-cm long, 316L stainless steel columns are equally spaced and welded to the inner canister wall at the lower end of the shell. These columns represent the FSV basket and can hold one FSV SNF element. An A92014 T6 aluminum spacer disc is bolted to the bottom plate of the PB2 basket to vertically restrain the FSV SNF element after the PB2 basket is placed inside the canister on top of the FSV basket. The structural evaluations of the ISS followed applicable ASME BPVC.III.3 guidelines and included finite element (FE) analyses of the PB2 basket structure; analyses of welds and bolds; buckling analyses of selected components, and acceptability assessments of the expected basket deformations under loading operations. The criticality safety assessments used the Monte Carlo N-Particle (MCNP) software architecture version 6.2 including ENDF/B-V continuous energy cross-section libraries, considering intact SNF in a single storage overpack or two multi-storage overpack configurations, and intact or failed SNF configured for disposal. The dose rate computations are based on source terms taken from the DOE Spent Fuel Database. The isotopic composition was decay corrected for the year 2022 using the ORIGEN module in the SCALE suite. A 19-group photon spectrum and 27-group neutron-source spectra were generated and used in MCNP to calculate estimated dose-equivalent rates, both on DOE Standard Canister contact and at a radial distance of from the canister surface. The results of this study indicate a structurally sound system that can uphold its criticality safety functions throughout its intended operational phases. Further, they increase confidence that sufficient radiological protection is technically achievable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Structural, Criticality, and Radiation Dose Calculations to Support SNF Loading into a DOE Standard Canister

The DOE Standard Canister Demonstration Project includes the development of an internal support structure (ISS) for the 4.6-m-long and 45.7-cm-diameter canister. This study presents structural evaluations, criticality safety assessments, and dose rate calculations conducted within the scope of the ISS design process to support smooth canister loading operations and to ensure safe storage, transportation, and disposal of Peach Bottom 1 Core II (PB2) and Fort St. Vrain (FSV) SNF, currently stored at INL’s CPP-603 facility. The ISS includes a 316L stainless-steel basket that holds 12 PB2 SNF rods. The PB2 basket rests on top of six individual, 78.7-cm-long, 316L stainless-steel columns that are equally spaced and welded to the inner canister wall at the lower end of the shell. These columns represent the FSV basket and can hold one FSV SNF element. An A92014 T6 aluminum spacer disc is bolted to the bottom plate of the PB2 basket to vertically restrain the FSV SNF element after the PB2 basket is placed inside the canister above the FSV basket. The structural evaluations of the ISS followed applicable ASME BPVC.III.3 guidelines and included finite element (FE) analyses of the PB2 basket structure; analyses of welds and bolds; buckling analyses of selected components, and acceptability assessments of the expected basket deformations under loading operations. The criticality safety assessments used the Monte Carlo N-Particle (MCNP) software architecture Version 6.2, including ENDF/B-V continuous-energy cross-section libraries, considering intact SNF in a single storage overpack or two multistorage overpack configurations and intact or failed SNF configured for disposal. The dose rate computations are based on source terms taken from the DOE Spent Fuel Database. The isotopic composition was decay corrected for the year 2022 using the ORIGEN module in the SCALE suite. A 19-group photon spectrum and a 27-group neutron-source spectrum were generated and used in MCNP to calculate estimated dose-equivalent rates, both on DOE Standard Canister contact and at a radial distance of 1 m from the canister surface. The results of this study indicate a structurally sound system that can uphold its criticality safety functions throughout its intended operational phases. Furthermore, this study provides confidence that sufficient radiological protection is technically achievable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The incidence of X-ray selected AGN in nearby galaxies

ABSTRACT We present the identification and analysis of an unbiased sample of active galactic nuclei (AGN) that lie within the local galaxy population. Using the MPA-JHU catalogue (based on SDSS DR8) and 3XMM DR7 we define a parent sample of 25 949 local galaxies (z ≤ 0.33). After confirming that there was strictly no AGN light contaminating stellar mass and star-formation rate calculations, we identified 917 galaxies with central, excess X-ray emission likely originating from an AGN. We analysed their optical emission lines using the BPT diagnostic and confirmed that such techniques are more effective at reliably identifying sources as AGN in higher mass galaxies: rising from 30 per cent agreement in the lowest mass bin to 93 per cent in the highest. We then calculated the growth rates of the black holes powering these AGN in terms of their specific accretion rates (∝LX/M*). Our sample exhibits a wide range of accretion rates, with the majority accreting at rates $\le 0.5\ \mathrm{ per \, cent}$ of their Eddington luminosity. Finally, we used our sample to calculate the incidence of AGN as a function of stellar mass and redshift. After correcting for the varying sensitivity of 3XMM, we split the galaxy sample by stellar mass and redshift and investigated the AGN fraction as a function of X-ray luminosity and specific black hole accretion rate. From this we found the fraction of galaxies hosting AGN above a fixed specific accretion rate limit of 10−3.5 is constant (at $\approx 1\ \mathrm{ per \, cent}$) over stellar masses of 8 < log M*/M⊙ < 12 and increases (from $\approx 1\ \mathrm{ per \, cent}$ to 10 per cent) with redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

Transformational challenge reactor design characteristics

The Transformational Challenge Reactor (TCR) program was conceived with the goal to reduce costs and time frames associated with advanced reactor deployment by leveraging developments in advanced manufacturing, advanced materials, data science, and rapid prototyping and testing. The final deliverable of the TCR program was to be an operational test of a novel reactor design. The TCR core design incorporates a dense tri-structural-isotropic/SiC fuel form and volumetrically efficient yttrium hydride moderator, both of which were manufactured and characterized under the TCR program. The TCR is a 3 MW{sub th} He-cooled experimental nuclear reactor designed to reach a total integrated burnup of less than 24 effective full-power hours to keep the radioactive source term to a very low level. TCR design process revealed a positive moderator coefficient; however, the negative doppler coefficients for the fuel and thermal expansion of fuel, moderator, and core support plate yield an overall negative reactivity coefficient. Calculated fuel element temperatures and stresses are well within safety margins. The maximum hypothetical accident (i.e., de-pressurized loss of forced cooling) yields only a modest increase in reactor temperatures that are all within safety margins. This paper summarizes the high-level TCR design characteristics, which were derived from neutronics, thermohydraulics, thermomechanics, and safety analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Engineering description of the ascent/descent bet product

The Ascent/Descent output product is produced in the OPIP routine from three files which constitute its input. One of these, OPIP.IN, contains mission specific parameters. Meteorological data, such as atmospheric wind velocities, temperatures, and density, are obtained from the second file, the Corrected Meteorological Data File (METDATA). The third file is the TRJATTDATA file which contains the time-tagged state vectors that combine trajectory information from the Best Estimate of Trajectory (BET) filter, LBRET5, and Best Estimate of Attitude (BEA) derived from IMU telemetry. Each term in the two output data files (BETDATA and the Navigation Block, or NAVBLK) are defined. The description of the BETDATA file includes an outline of the algorithm used to calculate each term. To facilitate describing the algorithms, a nomenclature is defined. The description of the nomenclature includes a definition of the coordinate systems used. The NAVBLK file contains navigation input parameters. Each term in NAVBLK is defined and its source is listed. The production of NAVBLK requires only two computational algorithms. These two algorithms, which compute the terms DELTA and RSUBO, are described. Finally, the distribution of data in the NAVBLK records is listed.

Seacord, A. W., II↗

Aeroacoustic computation of cylinder wake flow

Sound generation by uniform flow over a cylinder at the Reynolds number of 200 is calculated. The incompressible, time dependent flow field is first computed by standard computational fluid dynamics techniques using a stream function/vorticity formulation. The sound radiation is then obtained by integrating over the flow field with the Coriolis acceleration as the source term and a low frequency Green's function technique. Time histories, spectra, and directivity of the radiated sound are obtained. Comparisons with experimental data, where possible, indicate the feasibility of this computational aeroacoustics approach.

Hardin, J. C.↗

Comparison of Single-Event Transients in an Epitaxial Silicon Diode Resulting from Heavy Ion-, Focused X-Ray-, and Pulsed Laser-Induced Charge Generation

Heavy ion, focused x-ray, and pulsed laser single event transient experiments are performed on a silicon epitaxial diode. Collected charge, transient rise times, and transient fall times are calculated and compared between the different sources. It is observed that these transient shape characteristics depend on the source (ion, x-ray, or laser), even when similar amounts of charge are generated. Finally, the observed differences are examined and explained in terms of basic charge collection mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gravitational waves from rotating and precessing rigid bodies. 2: General solutions and computationally useful formulae

The classical mechanics results for free precession which are needed in order to calculate the weak field, slow-motion, quadrupole-moment gravitational waves are reviewed. Within that formalism, algorithms are given for computing the exact gravitational power radiated and waveforms produced by arbitrary rigid-body freely-precessing sources. The dominant terms are presented in series expansions of the waveforms for the case of an almost spherical object precessing with a small wobble angle. These series expansions, which retain the precise frequency dependence of the waves, may be useful for gravitational astronomers when freely-precessing sources begin to be observed.

Zimmerman, M.↗

The origin and evolution of atmospheric oxygen

This paper discusses the chemical processes involved in the evolution of the earth's atmospheric oxygen and ozone, as well as the sources, sinks, and transfer rates of oxygen in the present atmosphere. Special attention is given the evolution of atmospheric O3 as a function of the buildup of O2, with the results of calculations presented as the vertical profiles of O3, in terms of the present atmospheric level (PAL) oxygen values. Calculations show that the total O3 column density that is approximately half of the present level was reached when atmospheric oxygen level reached 0.1 PAL. At this level of ozone, the biological shielding of the earth's surface from the UV radiation is believed to have been achieved.

Levine, Joel S.↗

Infrared Limb Sounding with Cassini CIRS: Optimal Viewing Strategy Using Horizon Nodes

In this paper we investigate a question of science optimization during Cassini flybys of Titan. The Composite Infrared Spectrometer (CIRS) makes limb during the closest approach period when the visible horizon circle is moving swiftly across the planet. We have sought to discover if any points on the horizon are preferred for limb sounding due to having minimum movement relative to the surface. By numerical calculation, backed by geometric analysis, we find that two limited regions on the horizon are continuously visible during the entire encounter. We term these "limb nodes" and show how they may be employed by CIRS to optimize science by minimizing the source of systematic error due to spatial smear. These conclusions are applicable to many similar scenarios of spacecraft limb sounding during hyperbolic flyby encounters.

Mixon, Conor A.↗

The generation of gravitational waves. I - Weak-field sources

This paper derives and summarizes a 'plug-in-and-grind' formalism for calculating the gravitational waves emitted by any system with weak internal gravitational fields. If the internal fields have negligible influence on the system's motions, the formalism reduces to standard 'linearized theory'. Independent of the effects of gravity on the motions, the formalism reduces to the standard 'quadrupole-moment formalism' if the motions are slow and internal stresses are weak. In the general case, the formalism expresses the radiation in terms of a retarded Green's function for slightly curved spacetime and breaks the Green's function integral into five easily understood pieces: direct radiation, produced directly by the motions of the source; whump radiation, produced by the 'gravitational stresses' of the source; transition radiation, produced by a time-changing time delay ('Shapiro effect') in the propagation of the nonradiative 1/r field of the source; focusing radiation, produced when one portion of the source focuses, in a time-dependent way, the nonradiative field of another portion of the source; and tail radiation, produced by 'back-scatter' of the nonradiative field in regions of focusing.

Thorne, K. S.↗

CPU and memory efficient coherent mode decomposition for the partially coherent x-ray simulations

In this work, the method of the Coherent Mode Decomposition (CMD) is applied to numerical wave propagation calculations for partially-coherent X-rays, using the Fourier optics and compatible methods. Its CPU and memory efficiency is discussed in various cases of the wavefront at the source and the beam waist. With the absence of the quadratic phase terms, the required sampling density of the electric fields is effectively reduced. The problem size is thus moderate and the method is feasible to be implemented on a single-node CPU server. In other cases, the same argument holds with proper treatments of the quadratic phase terms. Tests on CMD and the modes propagation are done for the case of the Coherent Hard X-ray beamline of the National Synchrotron Light Source II, using the Synchrotron Radiation Workshop software. We observe a few hundred or less dominant decomposed modes that resemble the electric fields converge to the wavefront intensity at a high accuracy of over 99%.

36 MATERIALS SCIENCE↗

A morphological study of silicon carbide prepared by chemical vapor deposition

Surface morphologies of SiC deposits obtained via chemical vapor deposition (CVD) were studied as functions of substrate surface temperature and the concentration of the silicon and carbon source material, methyltrichlorosilane. Substrates of graphite and alpha-SiC crystals were used. Explanations of the observed morphologies on graphite substrates and their marked changes with temperature are given in terms of chemical kinetics and mass transport arguments. The results of thermodynamic calculations were used to help explain the observed morphologies of the deposits on alpha-SiC substrates.

Tsui, P.↗

The electron distribution function upstream from the earth's bow shock

A general analytic theory for the distribution function of particles backstreaming from an arbitrary shock in a magnetized plasma is presented. Particle motions are shown to be restricted to two-dimensional planes. A general form for the source term describing upstreaming particles accelerated at the shock is presented along with explicit source terms for planar and parabolic shocks. The origin and form of the escape velocity cutoff for a shock in a magnetized plasma are discussed. The distribution function upstream of a finite planar shock and a parabolic shock in two dimensions is calculated, and an analytic approximation of the cutoff velocity at points in the upstream region is derived. The theory is then applied to the earth's bow shock, calculating distribution functions and the nature and spatial variation of the cutoff velocity. The theory for the particle distribution upstream of the bow shock is compared with the model of Filbert and Kellogg (1979).

Cairns, Iver H.↗

Active Control of Fan Noise: Feasibility Study: Theoretical Analysis for Coupling of Active Noise Control Actuator Ring Sources to an Annular Duct with Flow - Volume 6

The objective of this effort is to develop an analytical model for the coupling of active noise control (ANC) piston-type actuators that are mounted flush to the inner and outer walls of an annular duct to the modes in the duct generated by the actuator motion. The analysis will be used to couple the ANC actuators to the modal analysis propagation computer program for the annular duct, to predict the effects of active suppression of fan-generated engine noise sources. This combined program will then be available to assist in the design or evaluation of ANC systems in fan engine annular exhaust ducts. An analysis has been developed to predict the modes generated in an annular duct due to the coupling of flush-mounted ring actuators on the inner and outer walls of the duct. The analysis has been combined with a previous analysis for the coupling of modes to a cylindrical duct in a FORTRAN computer program to perform the computations. The method includes the effects of uniform mean flow in the duct. The program can be used for design or evaluation purposes for active noise control hardware for turbofan engines. Predictions for some sample cases modeled after the geometry of the NASA Lewis ANC Fan indicate very efficient coupling in both the inlet and exhaust ducts for the m = 6 spinning mode at frequencies where only a single radial mode is cut-on. Radial mode content in higher order cut-off modes at the source plane and the required actuator displacement amplitude to achieve 110 dB SPL levels in the desired mode were predicted. Equivalent cases with and without flow were examined for the cylindrical and annular geometry, and little difference was found for a duct flow Mach number of 0.1. The actuator ring coupling program will be adapted as a subroutine to the cylindrical duct modal analysis and the exhaust duct modal analysis. This will allow the fan source to be defined in terms of characteristic modes at the fan source plane and predict the propagation to the arbitrarily-located ANC source plane. The actuator velocities can then be determined to generate the anti-phase mode. The resulting combined fan source/ANC pressure can then be calculated at any desired wall sensor position. The actuator velocities can be determined manually or using a simulation of a control system feedback loop. This will provide a very useful ANC system design and evaluation tool.

Kraft, R. E.↗