Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code of Accounts”

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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Tailoring composition and deformation modes at the microstructural level for next generation low-cost high-strength austenitic stainless steels

The objective of this project is to enable deliberate development of cost-effective, hydrogen resistant alloys by establishing detailed relationships specific to the effects of alloy composition, short-range order (SRO), and microsegregation in the presence of hydrogen on the transition between homogeneous deformation and localized plasticity in shear bands. In collaboration with the International Institute for Carbon-Neutral Energy Research, I2CNER, at Kyushu University in Japan, we conceptualized, designed, and manufactured four austenitic alloys that maintain corrosion resistance and ensure lower cost relative to baseline commercial alloys. The mechanical properties and deformation modes of the novel alloys (KU alloys) were assessed in the presence of hydrogen (H). Correlations between composition and performance revealed that two of the KU alloys are suitable replacements for 316 steel, while another is a viable replacement for 304 steel at room temperature. We found that, in the presence of other austenite stabilizing elements namely Mn and N, replacing Ni with Cu does not lead to martensite formation as has been previously reported.1–3 Furthermore, we found that the addition of Cu leads to an earlier onset of multiple slip resulting in an relative earlier onset of a higher work hardening rate (WHR). Greater understanding of the relationships between alloy composition and SRO required the development of a novel advanced electron diffraction methodology to characterize SRO in complex FCC alloys. This innovative approach, which combines fluctuation and correlation analyses of diffuse-scattering signals, successfully differentiated between SRO and long-range ordering (LRO). Further investigations into annealed austenitic stainless steels could provide insights into manipulating SRO and its effects on material properties. Atomistic simulations provided understanding of SRO behavior that was difficult to capture experimentally. This project created the first spin cluster expansion model that is able to capture and describe SRO effects in Fe-Ni-Cr FCC alloys, accounting for the non-negligible effects of magnetism. An automated computational workflow was established to provide reliable predictions of SRO in Fe-Ni-Cr austenitic alloys, both with and without the presence of H atoms. Analysis of the propensity for SRO in Fe-Ni-Cr alloys revealed that H tends to cluster with specific, well-defined SRO domains. The computational framework is general purpose and can be extended to realistic stainless steels across diverse composition ranges. With confidence that SRO is possible in austenitic stainless steels, we developed a discrete dislocation finite element code to understand the interaction of dislocations with SRO in the presence of H. By incorporating H effects on the dislocation emission and SRO stress field we show that the critical stress for the dislocation pileup to breakthrough the SRO domain decreases in the presence of H, which directly contributes localized deformation at the macroscale. Through the simulation of a uniaxial tension test, we demonstrated that H-induced weakening of SRO stress field and H-enhanced dislocation emission can lead to the onset of shear localization at lower macroscopic strains. As a whole, this project identified three novel alloys that show improvements in performance and cost efficiency for H-facing applications by studying correlations between alloy chemistry and deformation behavior. We also made significant advancements to experimental and computational methodologies necessary to study the chemistry and distribution of SRO across a range of alloys, which in turn allowed us to demonstrate how deformation mechanisms change due to the contributions of SRO in austenitic alloys in the presence of H. The combined advancements in fundamental understanding with novel alloy development in this project has increased the viability of next generation H-technologies for the broader public through accessible low-cost alloys and accelerated development towards future H-infrastructure.

08 HYDROGEN↗
NASA NTRS1994

A Theoretical Evaluation of Secondary Atomization Effects on Engine Performance for Aluminum Gel Propellants

A one-dimensional model of a gel-fueled rocket combustion chamber has been developed. This model includes the processes of liquid hydrocarbon burnout, secondary atomization. aluminum ignition, and aluminum combustion. Also included is a model of radiative heat transfer from the solid combustion products to the chamber walls. Calculations indicate that only modest secondary atomization is required to significantly reduce propellant burnout distances, aluminum oxide residual size and radiation heat wall losses. Radiation losses equal to approximately 2-13 percent of the energy released during combustion were estimated. A two-dimensional, two-phase nozzle code was employed to estimate radiation and nozzle two-phase flow effects on overall engine performance. Radiation losses yielded a 1 percent decrease in engine I(sub sp). Results also indicate that secondary atomization may have less effect on two-phase losses than it does on propellant burnout distance and no effect if oxide particle coagulation and shear induced droplet breakup govern oxide particle size. Engine I(sub sp) was found to decrease from 337.4 to 293.7 seconds as gel aluminum mass loading was varied from 0-70 wt percent. Engine I(sub sp) efficiencies, accounting for radiation and two-phase flow effects, on the order of 0.946 were calculated for a 60 wt percent gel, assuming a fragmentation ratio of 5.

Mueller, D. C.↗
NASA NTRS2006

PAB3D: Its History in the Use of Turbulence Models in the Simulation of Jet and Nozzle Flows

This is a review paper for PAB3D s history in the implementation of turbulence models for simulating jet and nozzle flows. We describe different turbulence models used in the simulation of subsonic and supersonic jet and nozzle flows. The time-averaged simulations use modified linear or nonlinear two-equation models to account for supersonic flow as well as high temperature mixing. Two multiscale-type turbulence models are used for unsteady flow simulations. These models require modifications to the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k-epsilon) model with a RANS/LES transition function, dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) formulation. All of these models are implemented in the three-dimensional Navier-Stokes code PAB3D. This paper discusses computational methods, code implementation, computed results for a wide range of nozzle configurations at various operating conditions, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions.

Abdol-Hamid, Khaled S.↗
NASA NTRS2016

Hall2De Simulations with a First-Principles Electron Transport Model Based on the Electron Cyclotron Drift Instability

Several years of work have combined plasma measurements in a laboratory Hall thruster and r-z numerical simulations with Hall2De to isolate the spatial variation of the anomalous collision frequency needed in Ohm’s law to produce the observed thruster behavior. This numerical solution is used here to test the validity of a first-principles model of the anomalous transport in these devices before such model is implemented self-consistently in r-z fluid codes like Hall2De. The first-principles model employs quasi-linear theory and is based on the hypothesis that the Electron Cyclotron Drift Instability (ECDI) excites ion acoustic turbulence that, in turn, enhances the effective collision frequency in these devices. We find that a model of the ECDI that assumes Maxwellian velocity distributions for electrons and singly-charged, main-beam, cold ions (Ti=0.07 eV) is insufficient to explain the expected variation of the anomalous collision frequency both in the interior and exterior of the acceleration channel. When warm ions (~0.5-3 eV) are accounted for, the ECDI model in the channel interior appears more promising but fails by orders of magnitude in the near plume region due to the much higher Landau damping of the ion acoustic waves there. This implies that either (a) some process allows the ECDI instability to remain uninhibited by classical Landau damping or, (b) that a different instability (or instabilities) altogether, also insusceptible to Landau damping, is/are active in this region. A previous hypothesis, that convection of wave energy generated by the ECDI in the channel plays a significant role in the near plume, is not supported by the results of the simulations.

Mikellides, Ioannis G.↗
DOE OSTI2016

SU-G-IeP2-04: Dosimetric Accuracy of a Monte Carlo-Based Tool for Cone-Beam CT Organ Dose Calculation: Validation Against OSL and XRQA2 Film Measurements

Purpose: The intensive use of Cone-Beam Computed Tomography (CBCT) during radiotherapy treatments raise some questions about the dose to healthy tissues delivered during image acquisitions. We hence developed a Monte Carlo (MC)-based tool to predict doses to organs delivered by the Elekta XVI kV-CBCT. This work aims at assessing the dosimetric accuracy of the MC tool, in all tissue types. Methods: The kV-CBCT MC model was developed using the PENELOPE code. The beam properties were validated against measured lateral and depth dose profiles in water, and energy spectra measured with a CdTe detector. The CBCT simulator accuracy then required verification in clinical conditions. For this, we compared calculated and experimental dose values obtained with OSL nanoDots and XRQA2 films inserted in CIRS anthropomorphic phantoms (male, female, and 5-year old child). Measurements were performed at different locations, including bone and lung structures, and for several acquisition protocols: lung, head-and-neck, and pelvis. OSLs and film measurements were corrected when possible for energy dependence, by taking into account for spectral variations between calibration and measurement conditions. Results: Comparisons between measured and MC dose values are summarized in table 1. A mean difference of 8.6% was achieved for OSLs when the energy correction was applied, and 89.3% of the 84 dose points were within uncertainty intervals, including those in bones and lungs. Results with XRQA2 are not as good, because incomplete information about electronic equilibrium in film layers hampered the application of a simple energy correction procedure. Furthermore, measured and calculated doses (Fig.1) are in agreement with the literature. Conclusion: The MC-based tool developed was validated with an extensive set of measurements, and enables the organ dose calculation with accuracy. It can now be used to compute and report doses to organs for clinical cases, and also to drive strategies to optimize imaging protocols.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

MBMS1.0: An Open-Source Code for Modeling and Simulation of Membrane-Based Dehumidification and Energy Recovery

Air dehumidification processes are important and significantly impact both the performance of building ventilation and air conditioning systems, and the health of occupants. In a membrane-based dehumidifier, a permselective membrane is adopted to separate the two flow streams at feed and permeate sides. The feed side is continuously supplied with humid air at atmospheric pressure, while the permeate side is connected to a vacuum pump to produce vacuum or a low permeate vapor pressure as a driving force for water vapor permeation. The software MBMS 1.0 enables modeling and simulation of the design and configuration optimization of a membrane-based air-to-air (A2A) or air-to-vapor (A2V) device. A first-order physical-constant model was developed to reasonably account for the fundamental heat and mass transfer processes involved.

Gao, Zhiming↗
DOE OSTI2019

Consensus Control (Wind Farm Collective Yaw Control) [SWR-19-11]

Wind turbines in a wind farm typically operate individually to maximize their own performance and do not take into account information from nearby turbines. In an autonomous wind farm, enabling cooperation to achieve farm-level objectives, turbines will need to use information from nearby turbines to optimize performance, ensure resiliency when other sensors fail, and adapt to changing local conditions. A key element of achieving an autonomous wind farm is to develop algorithms that provide necessary information to ensure reliable, robust, and efficient operation of wind turbines in a wind plant using local sensor information that is already being collected, such as supervisory control and data acquisition (SCADA) data, local meteorological stations, and nearby radars/sodars/lidars. This code implements a framework for implementing an autonomous wind farm that incorporates information from local sensors in real time to better align turbines in a wind farm. Oftentimes, measurements made at an individual turbine are noisy and unreliable. By incorporating measurements from multiple nearby turbines, a more robust estimate of the wind direction can be obtained at an individual turbine. Results indicate that this estimate of the wind direction can be used to improve the turbine's knowledge of the wind direction and could decrease dynamic yaw misalignment, decrease the amount of time a turbine spends yawing due to a more robust input to the yaw controller, and increase resiliency to faulty wind-vane measurements.

King, Jennifer↗
DOE OSTI2014

RAVEN. Dynamic Event Tree Approach Level III Milestone

Conventional Event-Tree (ET) based methodologies are extensively used as tools to perform reliability and safety assessment of complex and critical engineering systems. One of the disadvantages of these methods is that timing/sequencing of events and system dynamics are not explicitly accounted for in the analysis. In order to overcome these limitations several techniques, also know as Dynamic Probabilistic Risk Assessment (DPRA), have been developed. Monte-Carlo (MC) and Dynamic Event Tree (DET) are two of the most widely used D-PRA methodologies to perform safety assessment of Nuclear Power Plants (NPP). In the past two years, the Idaho National Laboratory (INL) has developed its own tool to perform Dynamic PRA: RAVEN (Reactor Analysis and Virtual control ENvironment). RAVEN has been designed to perform two main tasks: 1) control logic driver for the new Thermo-Hydraulic code RELAP-7 and 2) post-processing tool. In the first task, RAVEN acts as a deterministic controller in which the set of control logic laws (user defined) monitors the RELAP-7 simulation and controls the activation of specific systems. Moreover, the control logic infrastructure is used to model stochastic events, such as components failures, and perform uncertainty propagation. Such stochastic modeling is deployed using both MC and DET algorithms. In the second task, RAVEN processes the large amount of data generated by RELAP-7 using data-mining based algorithms. This report focuses on the analysis of dynamic stochastic systems using the newly developed RAVEN DET capability. As an example, a DPRA analysis, using DET, of a simplified pressurized water reactor for a Station Black-Out (SBO) scenario is presented.

42 ENGINEERING↗
DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1980

Destruction of N/2D/ by O2 - A major source of 6300 A dayglow emission

The paper discusses the N(2D) + O2 yields NO + O(1D), for k1 reaction, and examines under different geophysical conditions, the effect of this source of 6300 A surface brightness on the dayglow redline. In calculating the 6300 A volume emission rate, a computer code is used which provides the thermal electron density, photoelectron flux, ion concentrations, electron and ion temperatures and odd nitrogen densities required to evaluate the O(1D) sources and sinks. It is found that under most conditions this source is the dominant mechanism responsible for the production of the 6300 A dayglow above 150 km. The volume emission rate of this source in summer is 2.5 times larger than the sum of all other sources at 180 km, and 1.5 times as large in winter. In addition, when the N(2D) source of O(1D) is taken into account, the rate coefficient for quenching of O(1D) by N2 is estimated to be 4 plus or minus 1 x 10 to the -11th cu cm per sec.

Torr, D. G.↗
NASA NTRS2015

Constraints on the Gas Content of the Fomalhaut Debris Belt: Can Gas-Dust Interactions Explain the Belt's Morphology?

Context. The 440 Myr old main-sequence A-star Fomalhaut is surrounded by an eccentric debris belt with sharp edges. This sort of a morphology is usually attributed to planetary perturbations, but the orbit of the only planetary candidate detected so far, Fomalhaut b, is too eccentric to efficiently shape the belt. Alternative models that could account for the morphology without invoking a planet are stellar encounters and gas-dust interactions. Aims. We aim to test the possibility of gas-dust interactions as the origin of the observed morphology by putting upper limits on the total gas content of the Fomalhaut belt. Methods. We derive upper limits on the CII158 microns and OI 63 microns emission by using non-detections from the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel Space Observatory. Line fluxes are converted into total gas mass using the non-local thermodynamic equilibrium (non-LTE) code radex. We consider two different cases for the elemental abundances of the gas: solar abundances and abundances similar to those observed for the gas in the Beta Pictoris debris disc. Results. The gas mass is shown to be below the millimetre dust mass by a factor of at least approx. 3 (for solar abundances) respectively approx. 300 (for Beta Pic-like abundances). Conclusions. The lack of gas co-spatial with the dust implies that gas-dust interactions cannot efficiently shape the Fomalhaut debris belt. The morphology is therefore more likely due to a yet unseen planet (Fomalhaut c) or stellar encounters.

Fomalhaut b↗
NASA NTRS2018

On the Natural Structure of Amino Acid Patterns in Families of Protein Sequences

All known terrestrial proteins are coded as continuous strings of ≈20 amino acids. The patterns formed by the repetitions of elements in groups of finite sequences describes the natural architectures of protein families. We present a method to search for patterns and groupings of patterns in protein sequences using a mathematically precise definition for “repetition”, an efficient algorithmic implementation and a robust scoring system with no adjustable parameters. We show that the sequence patterns can be well-separated into disjoint classes according to their recurrence in nested structures. The statistics of the occurrences of patterns indicate that short repetitions are sufficient to account for the differences between natural families and randomized groups of sequences by more than 10 standard deviations, while contiguous sequence patterns shorter than 5 residues are effectively random in their occurrences. A small subset of patterns is sufficient to account for a robust ”familiarity” definition between arbitrary sets of sequences.

Pablo Turjanski↗
DOE OSTI2016

Dynamic Simulation with on the fly Doppler Broadening treatment in RMC

In conventional transient analysis, thermal hydraulics calculations are generally coupled to deterministic neutronics solvers. Due to the limited accuracy of deterministic methods and the rapid development of computer power, Monte Carlo neutronics calculations have received increasing attention nowadays. One of challenges for Monte Carlo methods coupled to thermal hydraulics calculations is the temperature dependent cross sections. Taking the detailed temperature distribution feedback into account, each nuclide may be in a large range of temperatures, pre-generated libraries at every temperature point are no longer applicable owing to unaffordable memory footprint. Different solutions have been proposed to deal with this problem. TMS, accounting for Doppler broadening effect by the explicit treatment of thermal motion, was used to model fast transients in Serpent 2. Stochastic mixing, using two sets of libraries for each nuclide and mixing them to get the correct temperature, was used in transient analysis by Sjenitzer and Hoogenboom. Besides, there are other on the fly Doppler broadening methods being developed, including series expansion of cross sections, multipole representation, different integration algorithms based on Doppler broadening equation. In this paper, on the fly Doppler broadening treatment based on integration algorithms was applied to transient calculation in Reactor Monte Carlo code RMC . In the following section, the integration algorithms implemented and dynamic simulation method in RMC is briefly introduced. Finally, a pure physical transient fuel assembly model was chosen for numerical validation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Effects of electron recirculation on a hard x-ray source observed during the interaction of a high intensity laser pulse with thin Au targets

The interaction of a high intensity laser pulse on the preplasma of a high-Z solid target produced by the pulse's pedestal generates high-energy electrons. These electrons subsequently penetrate inside the solid target and produce bremsstrahlung photons, generating an x-ray source which can be used for photonuclear studies or to radiograph high area density objects. The source characteristics are compared for targets with thin (20 μm) and thick (100 μm) Au foils on the Omega EP laser at Laboratory for Laser Energetics. Simulations using the particle-in-cell code CALDER show that for a 20 μm thickness Au target, electrons perform multiple round-trips in the target under the effect of the laser ponderomotive potential and the target electrostatic potential. These relativistic electrons have random transverse displacements, with respect to the target normal, attributed to electrostatic fluctuation fields. As a result, the x-ray spot size is increased by a factor 2 for thin target compared to thick targets, in agreement with experimental results. In addition, the computed doses agree with the measured ones provided that electron recirculation in the thin target is taken into account. A dose increase by a factor 1.7 is then computed by allowing for recirculation. In the 100 μm target case, on the other hand, this effect is found to be negligible.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

Electron-impact ionization of P-like ions forming Si-like ions

We have calculated electron-impact ionization (EII) for P-like systems from P to Zn{sup 15+} forming Si-like ions. The work was performed using the flexible atomic code (FAC) which is based on a distorted-wave approximation. All 3ℓ → nℓ' (n = 3-35) excitation-autoionization (EA) channels near the 3p direct ionization threshold and 2ℓ → nℓ' (n = 3-10) EA channels at the higher energies are included. Close attention has been paid to the detailed branching ratios. Our calculated total EII cross sections are compared both with previous FAC calculations, which omitted many of these EA channels, and with the available experimental results. Moreover, for Fe{sup 11+}, we find that part of the remaining discrepancies between our calculations and recent measurements can be accounted for by the inclusion of the resonant excitation double autoionization process. Lastly, at the temperatures where each ion is predicted to peak in abundances in collisional ionization equilibrium, the Maxwellian rate coefficients derived from our calculations differ by 50%-7% from the previous FAC rate coefficients, with the difference decreasing with increasing charge.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

SU-E-T-586: Field Size Dependence of Output Factor for Uniform Scanning Proton Beams: A Comparison of TPS Calculation, Measurement and Monte Carlo Simulation

Purpose: Output dependence on field size for uniform scanning beams, and the accuracy of treatment planning system (TPS) calculation are not well studied. The purpose of this work is to investigate the dependence of output on field size for uniform scanning beams and compare it among TPS calculation, measurements and Monte Carlo simulations. Methods: Field size dependence was studied using various field sizes between 2.5 cm diameter to 10 cm diameter. The field size factor was studied for a number of proton range and modulation combinations based on output at the center of spread out Bragg peak normalized to a 10 cm diameter field. Three methods were used and compared in this study: 1) TPS calculation, 2) ionization chamber measurement, and 3) Monte Carlos simulation. The XiO TPS (Electa, St. Louis) was used to calculate the output factor using a pencil beam algorithm; a pinpoint ionization chamber was used for measurements; and the Fluka code was used for Monte Carlo simulations. Results: The field size factor varied with proton beam parameters, such as range, modulation, and calibration depth, and could decrease over 10% from a 10 cm to 3 cm diameter field for a large range proton beam. The XiO TPS predicted the field size factor relatively well at large field size, but could differ from measurements by 5% or more for small field and large range beams. Monte Carlo simulations predicted the field size factor within 1.5% of measurements. Conclusion: Output factor can vary largely with field size, and needs to be accounted for accurate proton beam delivery. This is especially important for small field beams such as in stereotactic proton therapy, where the field size dependence is large and TPS calculation is inaccurate. Measurements or Monte Carlo simulations are recommended for output determination for such cases.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

Anisotropic Elastic Resonance Scattering model for the Neutron Transport equation

The resonance scattering transfer cross-section has been reformulated to account for anisotropic scattering in the center-of-mass of the neutron-nucleus system. The main innovation over previous implementations is the relaxation of the ubiquitous assumption of isotropic scattering in the center-of-mass and the actual effective use of scattering angle distributions from evaluated nuclear data files in the computation of the angular moments of the resonant scattering kernels. The formulas for the high order anisotropic moments in the laboratory system are also derived. A multi-group numerical formulation is derived and implemented into a module incorporated within the NJOY nuclear data processing code. An ultra-fine energy mesh cross section library was generated using these new theoretical models and then was used for fuel assembly calculations with the PARAGON lattice physics code. The results obtained indicate a strong effect of this new model on reactivity, multi-group fluxes and isotopic inventory during depletion.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Demonstration of MELCOR and MACCS Capabilities for Molten Salt Reactor Decay Heat Removal During both Normal Operations and Salt Spill Scenarios

This report provides a demonstration of MELCOR and MELCOR Accident Consequence Code System (MACCS) capabilities to perform a dose assessment for a Molten Salt Reactor (MSR) off-gas system. A primary containment system salt spill is used as the off-normal scenario, along with a normal operation dose assessment for comparison. This report discusses the tools, methods, and information used in this assessment so that it may be utilized as a starting point for future advanced reactor consequence analyses. This report also highlights several gaps, to include the need for reactor inventory information specific to advanced reactors, and the need for specific atmospheric transport models that take into account the unique deposition behaviors of tritium and carbon-14, and makes recommendations for closing these gaps. This report satisfies the DOE NE Milestone M4RD-21SN0601062.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗