Search NASA⌕ Search

SEARCH · Search NASA

Results for “Energy Code”

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,405 records · Page 78

On the Statistical Uncertainty of Monte Carlo-Calculated Scattering Sensitivities

Sensitivity coefficients calculated with Monte Carlo codes are widely used for nuclear data uncertainty quantification in the modeling and simulation of complex 3D reactor systems. This study systematically compares sensitivity coefficients and associated statistical uncertainties for the multiplication factor and fuel temperature reactivity across multiple Monte Carlo codes (SCALE/KENO, SCALE/Shift, MCNP, and Serpent) using simple models representing light-water reactors and advanced reactor concepts. For multiplication factor sensitivities, statistical uncertainties are generally acceptable, although scattering sensitivities show significantly larger statistical uncertainties than, for example, fission and capture reactions. Fuel temperature reactivity sensitivities show significantly larger statistical uncertainties across all reactions. Elastic scattering sensitivities are the most problematic: all Monte Carlo codes fail to resolve energy-dependent coefficients, and they produce dramatically different energy-collapsed values. Critically, the use of these sensitivity coefficients in nuclear data uncertainty propagation leads to reduced statistical uncertainties in individual uncertainty contributions. This can lead to the masking of unusable sensitivity coefficients and producing misleading uncertainty results. The findings of this study show that new or enhanced methods are needed to improve Monte Carlo elastic scattering sensitivity calculations. Additionally, this study shows the relevance of verifying sensitivity coefficients through direct perturbation calculations for individual nuclide reactions, instead of only for total cross sections as commonly done.

Bostelmann, Rike [ORNL] (ORCID:0000000165968088)↗

Coronal closed structures. IV - Hydrodynamical stability and response to heating perturbations

The response of magnetically confined atmospheres to perturbations in the temperature and density distribution, and the local heating rate by means of a one-dimensional time-dependent hydrodynamical code, which incorporates the full energy, momentum and mass conservation equations is studied. These studies extend the linear instability analysis of Habbal and Rosner (1979) into the finite-amplitude regime, and generalize the confined atmosphere models of Serio et al., to the time-dependent domain. The results show that closed coronal atmospheres are stable against finite-amplitude perturbations if the chromospheric response is taken into account; and observed correlated increases in coronal density and temperature can only be achieved under quiescent conditions by increasing the heat deposition rate relatively more in the chromosphere than in the corona.

Peres, G.↗

Electron cyclotron thruster new modeling results preparation for initial experiments

The following topics are discussed: a whistler-based electron cyclotron resonance heating (ECRH) thruster; cross-field coupling in the helicon approximation; wave propagation; wave structure; plasma density; wave absorption; the electron distribution function; isothermal and adiabatic plasma flow; ECRH thruster modeling; a PIC code model; electron temperature; electron energy; and initial experimental tests. The discussion is presented in vugraph form.

Hooper, E. Bickford↗

Trajectory Based Heating and Ablation Calculations for MESUR Pathfinder Aeroshell

Based on the geometry of Mars Environment Survey (MESUR) Pathfinder aeroshell and an estimated Mars entry trajectory, two-dimensional axisymmetric time dependent calculations have been obtained using GIANTS (Gauss-Siedel Implicit Aerothermodynamic Navier-Stokes code with Thermochemical Surface Conditions) code and CMA (Charring Material Thermal Response and Ablation) Program for heating analysis and heat shield material sizing. These two codes are interfaced using a loosely coupled technique. The flowfield and convective heat transfer coefficients are computed by the GIANTS code with a species balance condition for an ablating surface, and the time dependent in-depth conduction with surface blowing is simulated by the CMA code with a complete surface energy balance condition. In this study, SLA-561V has been selected as heat shield material. The solutions, including the minimum heat shield thicknesses over aeroshell forebody, pyrolysis gas blowing rates, surface heat fluxes and temperature distributions, flowfield, and in-depth temperature history of SLA-561V, are presented and discussed in detail.

Chen, Y. K.↗

Validity of the Aluminum Equivalent Approximation in Space Radiation Shielding

The origin of the aluminum equivalent shield approximation in space radiation analysis can be traced back to its roots in the early years of the NASA space programs (Mercury, Gemini and Apollo) wherein the primary radiobiological concern was the intense sources of ionizing radiation causing short term effects which was thought to jeopardize the safety of the crew and hence the mission. Herein, it is shown that the aluminum equivalent shield approximation, although reasonably well suited for that time period and to the application for which it was developed, is of questionable usefulness to the radiobiological concerns of routine space operations of the 21 st century which will include long stays onboard the International Space Station (ISS) and perhaps the moon. This is especially true for a risk based protection system, as appears imminent for deep space exploration where the long-term effects of Galactic Cosmic Ray (GCR) exposure is of primary concern. The present analysis demonstrates that sufficiently large errors in the interior particle environment of a spacecraft result from the use of the aluminum equivalent approximation, and such approximations should be avoided in future astronaut risk estimates. In this study, the aluminum equivalent approximation is evaluated as a means for estimating the particle environment within a spacecraft structure induced by the GCR radiation field. For comparison, the two extremes of the GCR environment, the 1977 solar minimum and the 2001 solar maximum, are considered. These environments are coupled to the Langley Research Center (LaRC) deterministic ionized particle transport code High charge (Z) and Energy TRaNsport (HZETRN), which propagates the GCR spectra for elements with charges (Z) in the range I <= Z <= 28 (H -- Ni) and secondary neutrons through selected target materials. The coupling of the GCR extremes to HZETRN allows for the examination of the induced environment within the interior' of an idealized spacecraft as approximated by a spherical shell shield, and the effects of the aluminum equivalent approximation for a good polymeric shield material such as genetic polyethylene (PE). The shield thickness is represented by a 25 g/cm spherical shell. Although one could imagine the progression to greater thickness, the current range will be sufficient to evaluate the qualitative usefulness of the aluminum equivalent approximation. Upon establishing the inaccuracies of the aluminum equivalent approximation through numerical simulations of the GCR radiation field attenuation for PE and aluminum equivalent PE spherical shells, we Anther present results for a limited set of commercially available, hydrogen rich, multifunctional polymeric constituents to assess the effect of the aluminum equivalent approximation on their radiation attenuation response as compared to the generic PE.

Badavi, Francis F.↗

Practical Applications of Cosmic Ray Science: Spacecraft, Aircraft, Ground Based Computation and Control Systems and Human Health and Safety

In this paper we review the discovery of cosmic ray effects on the performance and reliability of microelectronic systems as well as on human health and safety, as well as the development of the engineering and health science tools used to evaluate and mitigate cosmic ray effects in earth surface, atmospheric flight, and space flight environments. Three twentieth century technological developments, 1) high altitude commercial and military aircraft; 2) manned and unmanned spacecraft; and 3) increasingly complex and sensitive solid state micro-electronics systems, have driven an ongoing evolution of basic cosmic ray science into a set of practical engineering tools (e.g. ground based test methods as well as high energy particle transport and reaction codes) needed to design, test, and verify the safety and reliability of modern complex electronic systems as well as effects on human health and safety. The effects of primary cosmic ray particles, and secondary particle showers produced by nuclear reactions with spacecraft materials, can determine the design and verification processes (as well as the total dollar cost) for manned and unmanned spacecraft avionics systems. Similar considerations apply to commercial and military aircraft operating at high latitudes and altitudes near the atmospheric Pfotzer maximum. Even ground based computational and controls systems can be negatively affected by secondary particle showers at the Earth's surface, especially if the net target area of the sensitive electronic system components is large. Accumulation of both primary cosmic ray and secondary cosmic ray induced particle shower radiation dose is an important health and safety consideration for commercial or military air crews operating at high altitude/latitude and is also one of the most important factors presently limiting manned space flight operations beyond low-Earth orbit (LEO).

Atwell, William↗

Operational advances in ring current modeling using RAM-SCB

The Ring current Atmosphere interaction Model with Self-Consistently calculated 3D Magnetic field (RAM-SCB) combines a kinetic model of the ring current with a force-balanced model of the magnetospheric magnetic field to create an inner magnetospheric model that is magnetically self consistent. RAM-SCB produces a wealth of outputs that are valuable to space weather applications. For example, the anisotropic particle distribution of the KeV-energy population calculated by the code is key for predicting surface charging on spacecraft. Furthermore, radiation belt codes stand to benefit substantially from RAM-SCB calculated magnetic field values and plasma wave growth rates - both important for determining the evolution of relativistic electron populations. RAM-SCB is undergoing development to bring these benefits to the space weather community. Data-model validation efforts are underway to assess the performance of the system. 'Virtual Satellite' capability has been added to yield satellite-specific particle distribution and magnetic field output. The code's outer boundary is being expanded to 10 Earth Radii to encompass previously neglected geosynchronous orbits and allow the code to be driven completely by either empirical or first-principles based inputs. These advances are culminating towards a new, real-time version of the code, rtRAM-SCB, that can monitor the inner magnetosphere conditions on both a global and spacecraft-specific level. This paper summarizes these new features as well as the benefits they provide the space weather community.

47 OTHER INSTRUMENTATION↗

MC21/CTF and VERA-CS Mutiphysics Solutions to VERA Core Physics Benchmark Problems 6 and 7

The continuous energy Monte Carlo neutron transport code, MC21, was coupled to the CTF subchannel thermal-hydraulics code using a combination of CASL tools and in-house Python scripts. An MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 6 demonstrated excellent agreement with MC21/COBRA-IE and VERA solutions. The MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 7, Watts Bar Unit 1 at beginning of cycle (BOC) hot full power (HFP) equilibrium xenon conditions, is the first published coupled Monte Carlo neutronics / subchannel T-H solution for this problem. MC21/CTF predicted a critical boron concentration of 854.5 ppm, yielding a critical eigenvalue of 0.99994 ± 6.8E-6 (95% CI). Excellent agreement with a VERA solution of Problem 7 was also demonstrated for integral and local power and temperature parameters.

Aviles, Brian↗

Analysis of SCALE Criticality and Sensitivity Calculations for Reflected HEU Cylinders [Abstract]

The SCALE code package offers several nuclear data libraries to support Monte Carlo (MC) transport, as well as MC-based derivation of $\kappa$ eff sensitivity and uncertainty (S/U) data. The CSAS sequence using the KENO MC code can utilize continuous-energy (CE) cross sections, or pre-generated multigroup (MG) cross section libraries. The use of MG libraries introduces bias into calculations in exchange for faster transport solutions. The TSUNAMI-3D sequence also utilizes KENO MC calculations. TSUNAMI-3D has two CE calculational methods: the Iterated Fission Probability (IFP) method, and the Contribution-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method. Previous work has shown poor agreement between CLUTCH and confirmatory direct perturbation calculations in specific applications, e.g., fissionable and polyethylene reflectors

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of burnup and isotopic compositions of BWR 9 x 9 UO{sub 2} fuel assemblies

In order to extend isotopic composition data focusing on fission product nuclides, measurements are progressing using facilities of JAEA for five samples taken from high burnup BWR 9 x 9 UO{sub 2} fuel assemblies. Neutronics analysis with an infinite assembly model was applied to the preliminary measurement data using a continuous-energy Monte Carlo burnup calculation code MVP-BURN with nuclear libraries based on JENDL-3.3 and JENDL-4.0. The burnups of the samples were determined to be 28.0, 39.3, 56.6, 68.1, and 64.0 GWd/t by the Nd-148 method. They were compared with those calculated using node-average irradiation histories of power and in-channel void fractions which were taken from the plant data. The comparison results showed that the deviations of the calculated burnups from the measurements were -4 to 3%. It was confirmed that adopting the nuclear data library based on JENDL-4.0 reduced the deviations of the calculated isotopic compositions from the measurements for {sup 238}Pu, {sup 144}Nd, {sup 145}Nd, {sup 146}Nd, {sup 148}Nd, {sup 134}Cs, {sup 154}Eu, {sup 152}Sm, {sup 154}Gd, and {sup 157}Gd. On the other hand, the effect of the revision in the nuclear. data library on the neutronics analysis was not significant for major U and Pu isotopes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

New Features of TRIPOLI-4{sup R} Version 10

TRIPOLI-4{sup R} is a three-dimensional and continuous-energy Monte-Carlo transport code developed by the 'Service d'Etudes des Reacteurs et de Mathematiques Appliquees' (SERMA), at the CEA Paris-Saclay Center. It is devoted to shielding, reactor physics, criticality safety and nuclear instrumentation. TRIPOLI-4{sup R} has been developed from the mid 1990's in C++. A new version of the code is typically released every two years. In this paper, we present the new features of TRIPOLI- 4{sup R} version 10, to be soon released to NEA and RSICC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Calculating Control Rod Differential Worth Using Adjoint-Weighted Method

Recently, there are more and more research on using Monte Carlo method to calculate reactor important parameters. Control rod worth is one of the important parameters as it can determine the critical control rod position. There are two presentations of the control rod worth for utilization purpose, integral and differential. The former is the total reactivity change, while the latter is the reactivity change caused by per unit movement. In Monte Carlo calculation, differential control rod worth of a particular position is estimated by finite difference of reactivity between a small gap divided by gap distance. Though easy to implemented, the problem of this approach lies in that too large gap will enlarge error between 'differential' and 'difference', while too small gap will be disturbed by error caused by random fluctuation. Monte Carlo perturbation theory has been used to calculate the k{sub eff} sensitivity to certain parameters for a long time. By calculating the differential coefficient along with the k{sub eff}, this method is quite time-consuming and relatively accurate compared to finite difference method. Nevertheless, most of perturbation research focus on the perturbation of nuclide density or nuclear data. For geometrical perturbation, limited research was done before. Takahashi estimated the k{sub eff} sensitivity to reflector movement of pulsed fast reactor. MCNP6 used adjoint-weighted method for Godiva, Zeus and uranium oxyfluoride solution container benchmarks. This paper tries to apply the perturbation method to the determination of control rod differential worth. In this paper, based on RMC, a continuous-energy Reactor Monte Carlo transport code, Iterated Fission Probability method (IFP) is researched to calculate control rod differential worth on assembly level. IFP was originally developed by Hurwitz, and has been adopted by most of Monte Carlo code to compute adjoint-weighted parameters. This paper is organized with four parts. The first is introduction. Adjoint-weighted method is introduced in the second part. Problem 4 of benchmark from CASL is built to verify the method. Conclusion is given in final part.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Local Power Impact Experiment Design for a New Fuel Type for use in the Advanced Test Reactor

The Advanced Test Reactor (ATR), and complimentary zero-power ATR Critical (ATRC) reactor, located at Idaho National Labs (INL), are undergoing conversion from Highly Enriched Uranium (HEU) to Low Enriched Uranium (LEU). Both have a variety of testing locations that can receive large variations in flux due to its unique serpentine design, consisting of five lobes (see Figure 1). Initial criticality and power distribution throughout the core are controlled by core-external outer shim control cylinders (OSCCs). Distinct test loops allow for testing at specific temperatures, pressures, and irradiation conditions. The ATR is one of the key nuclear engineering research and testing facilities within the DOE National Laboratory Complex, and the ATRC supports its operation [1]. Currently, the Office of Material Management and Minimization (M3) within the National Nuclear Security Administration of the DOE is working to convert the remaining research reactors, including the ATR, from 93% HEU fuel to 19.75% LEU fuel (LEU) to support non-proliferation [2]. Extensive materials testing at INL and internationally has demonstrated that a high-density uranium molybdenum (U 10Mo) alloy can meet the performance requirements of the remaining high powered research reactors. However, there are many technical challenges to address before the conversion to LEU can be successful, including the accurate characterization of the reactor core physics with LEU fuel. Reactor physics safety evaluations currently use Monte Carlo for the 21st Century (MC21), a continuous-energy Monte Carlo radiation transport code [3]. Existing MC21 models of the ATR and ATRC cores have a validation basis for use in neutronics analyses with HEU fuel. The models are used to support safety analyses that include comparisons to the safety requirements for the reactors. However, the use of the LOWE element in the ATR and ATRC is not currently covered by the current model validation basis. To deploy the new fuel type, extensive computational reactor physics support is necessary to support the use of LOWE in the ATR and ATRC. Therefore, LOWE requires a rigorous validation basis, aligned with that of HEU fuel, that takes advantage of the existing software tools and processes currently used for the ATR and ATRC. The experiment to validate of the MC21 models for determining power, the Power Impact Validation Experiment, will consist of two flux runs in the ATRC, one with fully HEU loading and one with a single LOWE element. Both flux runs will be instrumented with 20 sets of azimuthal fission wires and 3 sets of axial fission wires, as shown in Figure 4. Standard flux run methodology will be used [4]. Power Impact Validation Experiment data will be compared against MC21 calculated data, both for absolute fission rate accuracy and to determine the relative change in fission rates between the two runs. The results of the Power Impact Validation Experiment and subsequent evaluations will provide the validation basis for MC21 for use with LOWE elements. Key features of the Power Impact Validation Experiment include: (1) Two flux runs to allow for LOWE perturbed measurements to be compared to already validated measurements taken from a full core of HEU fuel, (2) Optimization of instrumentation to balance analytical needs with practical considerations (e.g., limited time window to count beta particles from fission products), and (3) Standard ATRC core loading, including both driver positions and flux traps, to minimize cost while remaining representative of typical ATR core loading.

42 ENGINEERING↗

Portability: A Necessary Approach for Future Scientific Software

Today's world of scientific software for High Energy Physics (HEP) is powered by x86 code, while the future will be much more reliant on accelerators like GPUs and FPGAs. The portable parallelization strategies (PPS) project of the High Energy Physics Center for Computational Excellence (HEP/CCE) is investigating solutions for portability techniques that will allow the coding of an algorithm once, and the ability to execute it on a variety of hardware products from many vendors, especially including accelerators. We think without these solutions, the scientific success of our experiments and endeavors is in danger, as software development could be expert driven and costly to be able to run on available hardware infrastructure. We think the best solution for the community would be an extension to the C++ standard with a very low entry bar for users, supporting all hardware forms and vendors. We are very far from that ideal though. We argue that in the future, as a community, we need to request and work on portability solutions and strive to reach this ideal.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Experimental and Monte Carlo studies of fluence corrections for graphite calorimetry in low- and high-energy clinical proton beams

Purpose: The aim of this study was to determine fluence corrections necessary to convert absorbed dose to graphite, measured by graphite calorimetry, to absorbed dose to water. Fluence corrections were obtained from experiments and Monte Carlo simulations in low- and high-energy proton beams. Methods: Fluence corrections were calculated to account for the difference in fluence between water and graphite at equivalent depths. Measurements were performed with narrow proton beams. Plane-parallel-plate ionization chambers with a large collecting area compared to the beam diameter were used to intercept the whole beam. High- and low-energy proton beams were provided by a scanning and double scattering delivery system, respectively. A mathematical formalism was established to relate fluence corrections derived from Monte Carlo simulations, using the FLUKA code [A. Ferrari et al., “FLUKA: A multi-particle transport code,” in CERN 2005-10, INFN/TC 05/11, SLAC-R-773 (2005) and T. T. Böhlen et al., “The FLUKA Code: Developments and challenges for high energy and medical applications,” Nucl. Data Sheets 120, 211–214 (2014)], to partial fluence corrections measured experimentally. Results: A good agreement was found between the partial fluence corrections derived by Monte Carlo simulations and those determined experimentally. For a high-energy beam of 180 MeV, the fluence corrections from Monte Carlo simulations were found to increase from 0.99 to 1.04 with depth. In the case of a low-energy beam of 60 MeV, the magnitude of fluence corrections was approximately 0.99 at all depths when calculated in the sensitive area of the chamber used in the experiments. Fluence correction calculations were also performed for a larger area and found to increase from 0.99 at the surface to 1.01 at greater depths. Conclusions: Fluence corrections obtained experimentally are partial fluence corrections because they account for differences in the primary and part of the secondary particle fluence. A correction factor, F(d), has been established to relate fluence corrections defined theoretically to partial fluence corrections derived experimentally. The findings presented here are also relevant to water and tissue-equivalent-plastic materials given their carbon content.

60 APPLIED LIFE SCIENCES↗

Intense Magnetized Plasma-Wall Interaction

This research project studied wall-plasma interactions relevant to fusion science. Such interactions are a critical aspect of Magneto-Inertial Fusion (MIF) because flux compression by a pusher material, in particular the metal for the liner approach to MIF, involves strong eddy current heating on the surface of the pusher, and probably interactions and mixing of the pusher with the interior fuel during the time when fusion fuel is being burned. When the pusher material is a metal liner, high-energy-density conditions result in fascinating behavior. For example, "warm dense matter" is produced, for which material properties such as resistivity and opacity are not well known. In this project, the transformation into plasma of metal walls subjected to pulsed megagauss magnetic fields was studied with an experiment driven by the UNR 1 MA Zebra generator. The experiment was numerically simulated with using the MHRDR code. This simple, fundamental high-energy-density physics experiment, in a regime appropriate to MIF, has stimulated an important and fascinating comparison of numerical modeling codes and tables with experiment. In addition, we participated in developing the FRCHX experiment to compress a field-reversed-configuration (FRC) plasma with a liner, in collaboration with researchers from Air Force Research Laboratory and Los Alamos National Lab, and we helped develop diagnostics for the Plasma Liner Experiment (PLX) at LANL. Last, but not least, this project served to train students in high-energy-density physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation as a sustainable trading zone: Aiming at intergenerational justice

The paper, drawing on the example of simulation codes used in nuclear physics and high-energy physics, seeks to call attention to the ethical implications of discontinuing support for simulation codes and the loss of knowledge embodied in them. Predicated on the concept of trading zones and actor network models, the paper addresses the problem of extinction of simulation codes and attempts to understand their evolution and development within those frameworks. We show that simulation codes of closed type develop to the level of creoles, becoming local languages and standards of scientific centers and disappearing as their few main developers leave, whereas codes of open types become universal languages, imposing problemsolving patterns on the entire community and crowding out other codes. The paper suggests that because of simulations’ reliance on tacit knowledge, practices entrenched in codes cannot be exhaustively explicated or transmitted through writing alone; on the contrary, the life cycle of a simulation code is determined by the life cycle of its trading zone. We examine the extent to which both of these phenomena pose a risk to the preservation of knowledge. Bearing upon intergenerational ethics, we draw analogies between the pure intergenerational problem (PIP) and the problem of preserving the knowledge implemented in simulation codes and transmitting it to future generations. We argue that for the complete transfer of knowledge, it is necessary to develop and maintain inhabitability and sustainability of simulation trading zones in a controllable way, at least until the demand for these codes is warranted to cease in the future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coulomb effects in low-energy nuclear fragmentation

Early versions of the Langley nuclear fragmentation code NUCFRAG (and a publicly released version called HZEFRG1) assumed straight-line trajectories throughout the interaction. As a consequence, NUCFRAG and HZEFRG1 give unrealistic cross sections for large mass removal from the projectile and target at low energies. A correction for the distortion of the trajectory by the nuclear Coulomb fields is used to derive fragmentation cross sections. A simple energy-loss term is applied to estimate the energy downshifts that greatly alter the Coulomb trajectory at low energy. The results, which are far more realistic than prior versions of the code, should provide the data base for future transport calculations. The systematic behavior of charge-removal cross sections compares favorably with results from low-energy experiments.

Wilson, John W.↗