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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

An Improved Neutron Transport Algorithm for HZETRN

Long term human presence in space requires the inclusion of radiation constraints in mission planning and the design of shielding materials, structures, and vehicles. In this paper, the numerical error associated with energy discretization in HZETRN is addressed. An inadequate numerical integration scheme in the transport algorithm is shown to produce large errors in the low energy portion of the neutron and light ion fluence spectra. It is further shown that the errors result from the narrow energy domain of the neutron elastic cross section spectral distributions, and that an extremely fine energy grid is required to resolve the problem under the current formulation. Two numerical methods are developed to provide adequate resolution in the energy domain and more accurately resolve the neutron elastic interactions. Convergence testing is completed by running the code for various environments and shielding materials with various energy grids to ensure stability of the newly implemented method.

Slaba, Tony C.↗
DOE OSTI2024

The Impact of the Halden Reactor Project on the Development of the Bison Fuel Performance Code

The Bison fuel performance code has been under development at Idaho National Laboratory by various programs within the U.S. Department of Energy since 2009. As part of any fuel performance code development, validation is necessary. To this end, several Halden experiments have been used. Early on, the focus was on validating the thermal mechanical behaviour, including pellet cladding mechanical interaction and fission gas release models implemented in Bison. Support for design-basis accident capabilities began in approximately 2014. Also in 2014, the development of models for accident tolerant fuel (ATF) concepts also began. This includes dopedUO2 fuel for which Halden has some experiments with temperature and fission gas release measurements. The Bison team has been fortunate to have two senior developers be secondees to Halden for a year: the late Dr. Giovanni Pastore from May 2015 to April 2016 and the now retired Dr. Richard Williamson in 2018. The paper will provide examples of Bison validation using Halden data for standard UO2 fuel under normal operation and loss-of-coolant accident condition, and ATF concepts. A tribute to the legacy that both the late Giovanni and retired Richard left on the Bison code is provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Preliminary neutronic and thermal analysis of Inertial Fusion Energy systems with a focus on thick liquid wall concepts

This paper provides an assessment of the neutronic behavior of Inertial Fusion Energy systems using the Serpent 2 Monte Carlo code. The study explores different types of reactor concepts, including dry-wall and thick liquid wall systems. In particular, it characterizes the tritium breeding and energy multiplication performance, as well as the effect of long-term irradiation on structural materials. For this purpose, the displacement per atom (DPA) and the gas production rate in structural materials are calculated. A simplified geometry with varying geometrical parameters is used as a representation for the different kinds of devices. Additionally, the case of the HYLIFE-II concept is used to assess modeling inaccuracies arising from certain approximations. A comparison between the simplified geometry and a realistic one is performed, thanks to Serpent’s capability to run transport calculations on unstructured meshes. The study then addresses the impact of these discrepancies on the temperature distribution in the first wall using the GeN-Foam code for conjugate heat transfer analysis. The results highlight the necessity to consider the heterogeneous nature of the liquid wall, as simplified models tend to underestimate radiation damage metrics and heat loads.

Conjugate heat transfer↗
DOE OSTI2025

Enhanced Monte Carlo Simulations for Electron Energy Loss Mitigation in Real-Space Nanoimaging of Thick Biological Samples and Microchips

High-resolution imaging using Transmission Electron Microscopy (TEM) is essential for applications such as grain boundary analysis, microchip defect characterization, and biological imaging. However, TEM images are often compromised by electron energy spread and other factors. In TEM mode, where the objective and projector lenses are positioned downstream of the sample, electron–sample interactions cause energy loss, which adversely impacts image quality and resolution. This study introduces a simulation tool to estimate the electron energy loss spectrum (EELS) as a function of sample thickness, covering electron beam energies from 300 keV to 3 MeV. Leveraging recent advances in MeV-TEM/STEM technology, which includes a state-of-the-art electron source with 2-picometer emittance, an energy spread of 3 × 10 -5 , and optimized beam characteristics, we aim to minimize energy spread. By integrating EELS capabilities into the BNL Monte Carlo (MC) simulation code for thicker samples, we evaluate electron beam parameters to mitigate energy spread resulting from electron–sample interactions. Based on our simulations, we propose an experimental procedure for quantitively distinguishing between elastic and inelastic scattering. The findings will guide the selection of optimal beam settings, thereby enhancing resolution for nanoimaging of thick biological samples and microchips.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nuclear Computational Low Energy Initiative (NUCLEI)

The NUCLEI project, as defined by the scope of work, developed, implemented and run codes for large-scale computations of many topics in low-energy nuclear physics. Physics studied include the properties of nuclei and nuclear decays, nuclear structure and reactions, and the properties of nuclear matter. The computational techniques used include Quantum Monte Carlo, Configuration Interaction, Coupled Cluster, and Density Functional methods. The research program emphasized areas of high interest to current and possible future DOE nuclear physics facilities, including ATLAS and FRIB (nuclear structure and reactions, and nuclear astrophysics), TJNAF (neutron distributions in nuclei, few body systems, and electroweak processes), NIF (thermonuclear reactions), MAJORANA and FNPB (neutrinoless double-beta decay and physics beyond the Standard Model), and LANSCE (fission studies).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2006

NASA Radiation Protection Research for Exploration Missions

The HZETRN code was used in recent trade studies for renewed lunar exploration and currently used in engineering development of the next generation of space vehicles, habitats, and EVA equipment. A new version of the HZETRN code capable of simulating high charge and energy (HZE) ions, light-ions and neutrons with either laboratory or space boundary conditions with enhanced neutron and light-ion propagation is under development. Atomic and nuclear model requirements to support that development will be discussed. Such engineering design codes require establishing validation processes using laboratory ion beams and space flight measurements in realistic geometries. We discuss limitations of code validation due to the currently available data and recommend priorities for new data sets.

Wilson, John W.↗
DOE OSTI2025

PETSc/TAO developments for GPU-based early exascale systems

The Portable Extensible Toolkit for Scientific Computation (PETSc) library provides scalable solvers for nonlinear time-dependent differential and algebraic equations and for numerical optimization via the Toolkit for Advanced Optimization (TAO). PETSc is used in dozens of scientific fields and is an important building block for many simulation codes. During the U.S. Department of Energy’s Exascale Computing Project, the PETSc team has made substantial efforts to enable efficient utilization of the massive fine-grain parallelism present within exascale compute nodes and to enable performance portability across exascale architectures. We recap some of the challenges that designers of numerical libraries face in such an endeavor, and then discuss the many developments we have made, which include the addition of new GPU backends, features supporting efficient on-device matrix assembly, better support for asynchronicity and GPU kernel concurrency, and new communication infrastructure. In conclusion, we evaluate the performance of these developments on some pre-exascale systems as well as the early exascale systems Frontier and Aurora, using compute kernel, communication layer, solver, and mini-application benchmark studies, and then close with a few observations drawn from our experiences on the tension between portable performance and other goals of numerical libraries.

Exascale Computing Project (ECP)↗
DOE OSTI2024

Proton beam power limits for stationary water-cooled tungsten target with different cladding materials

The proton beam power limit for a solid-tungsten spallation target is largely determined by beam induced thermomechanical structural loads and decay heat power deposition, while its lifetime is limited by radiation damage and fatigue life of the target materials. In this paper, we studied the power limits of a stationary water-cooled solid tungsten target concept. Tantalum clad tungsten was considered as a reference case. Being a low activation material, zircaloy 2 cladding option was studied and its decay heat driven power limit was compared with the reference case. Zirconium alloys have proven operations records in spallation target and nuclear fission environments, supported by materials data obtained from post irradiation examinations. Recent study also demonstrated feasibility of diffusion bonding zirconium to tungsten using vanadium foil inter layer. Particle transport simulations code FLUKA was used to calculate energy deposition and decay heat power deposition in the target, based on the beam parameters technically feasible at the Second Target Station of the Spallation Neutron Source at Oak Ridge National Laboratory. The energy deposition data were used for flow, thermal, and structural analyses to determine the beam intensity limit on the target concept studied. The decay heat deposition data were used to calculate the transient temperature evolution in the tungsten volumes in a loss of coolant accident (LOCA) scenario to determine its beam power limit. For a 1.3 GeV proton beam, the power limit on a stationary target was 400 kW for a tantalum clad target model and 800 kW for a zircaloy 2 clad target model.

Lee, Yong Joong↗
DOE OSTI2024

NEAMS Reactor Physics Project Management and Coordination Activities

The following report summarizes the coordination activities and collaboration efforts by research staff at Oak Ridge National Laboratory, Idaho National Laboratory, and Argonne National Laboratory for the development of Shift and Griffin computer codes under the US Department of Energy’s Nuclear Energy Advanced Modeling and Simulation program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Thermochemical Nonequilibrium Modeling in a Continuous-Galerkin, Finite-Element Framework

The presented work discusses the implementation, verification, and validation of Park's two-temperature model in a scalable, computational fluid dynamics (CFD) code developed at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The implementation of Park's two-temperature model was verified through 0D test cases involving an adiabatic reactor and a nitrogen thermal bath. The implementation was then validated through comparisons with other validated CFD codes and experimental data on a hypersonic cylinder and double cones. These are standard benchmark test cases for thermochemical non-equilibrium (TCNE) modeling, and all data are shared publicly. The verification and validation results showed that ORNL's in-house CFD code could model complex, high-speed flow problems with and without TCNE modeling. This work is essential for future research involving 3D shock wave/boundary layer interactions (SBLIs).

Nutter, Nicole↗
NASA NTRS1974

Developing a gas rocket performance prediction technique

A simple, semi-empirical performance correlation/prediction technique applicable to gaseous and liquid propellant rocket engines is presented. Excellent correlations were attained for over 100 test firings by adjusting the computation of the gaseous mixing of an unreactive, coaxial jet using a correlation factor, F, which resulted in prediction of the experimental combustion efficiency for each firing. Static pressure, mean velocity and turbulence intensity in the developing region of non-reactive coaxial jets, typical of those of coaxial injector elements were determined. Detailed profiles were obtained at twelve axial locations (extending from the nozzle exit for a distance of five diameters) downstream from a single element of the Bell Aerospace H2/O2 19-element coaxial injector. These data are compared with analytical predictions made using both eddy viscosity and turbulence kinetic energy mixing models and available computer codes. Comparisons were disappointing, demonstrating the necessity of developing improved turbulence models and computational techniques before detailed predictions of practical coaxial free jet flows are attempted.

Morgenthaler, J. H.↗
NASA NTRS1984

Magnetically suspended flywheel system study

A program to study the application of a graphite/epoxy, magnetically suspended, pierced disk flywheel for the combined function of spacecraft attitude control and energy storage (ACES) is described. Past achievements of the program include design and analysis computer codes for the flywheel rotor, a magnetically suspended flywheel model, and graphite/epoxy rotor rings that were successfully prestressed via interference assembly. All hardware successfully demonstrated operation of the necessary subsystems which form a complete ACES design. Areas of future work include additional rotor design research, system definition and control strategies, prototype development, and design/construction of a UM/GSFC spin test facility. The results of applying design and analysis computer codes to a magnetically suspended interference assembled rotor show specific energy densities of 42 Wh/lb (92.4 Wh/kg) are obtained for a 1.6 kWh system.

Kirk, J. A.↗
NASA NTRS1985

The heating of the solar wind by the interstellar neutral gas

Solar wind heating by ionization and assimilation of the interstellar neutral wind is investigated using a numerical model of the interaction. The model is a time-dependent, one-dimensional, spherically symmetric, one-fluid code which includes mass, momentum, and energy sources due to the efficient incorporation of ionized interstellar hydrogen. Solar wind data at 1 AU from October 12, 1978 to February 25, 1980 were input to the inner boundary of the model and propagated out to 10 AU. Then, 52-day averages of proton temperature as functions of radius were produced for three values of the interstellar hydrogen density: 0.0, 0.03, and 0.1 per cu cm. It is concluded that, within the context of the model, the observations of solar wind proton temperatures at Voyager 1 and Pioneer 11 are consistent with heating of the solar wind by an interaction with inflowing interstellar neutral hydrogen. The density of this inflowing gas is near 0.03/cu cm.

Isenberg, P. A.↗
NASA NTRS1987

An implicit flux-split algorithm to calculate hypersonic flowfields in chemical equilibrium

An implicit, finite-difference, shock-capturing algorithm that calculates inviscid, hypersonic flows in chemical equilibrium is presented. The flux vectors and flux Jacobians are differenced using a first-order, flux-split technique. The equilibrium composition of the gas is determined by minimizing the Gibbs free energy at every node point. The code is validated by comparing results over an axisymmetric hemisphere against previously published results. The algorithm is also applied to more practical configurations. The accuracy, stability, and versatility of the algorithm have been promising.

Palmer, Grant↗
NASA NTRS1988

BRYNTRN: A baryon transport computer code, computation procedures and data base

The development is described of an interaction data base and a numerical solution to the transport of baryons through the arbitrary shield material based on a straight ahead approximation of the Boltzmann equation. The code is most accurate for continuous energy boundary values but gives reasonable results for discrete spectra at the boundary with even a relatively coarse energy grid (30 points) and large spatial increments (1 cm in H2O).

Wilson, John W.↗
NASA NTRS1988

The efficiency of convective energy transport in the sun

Mixing length theory (MLT) utilizes adiabatic expansion (as well as radiative transport) to diminish the energy content of rising convective elements. Thus in MLT, the rising elements lose their energy to the environment most efficiently and consequently transport heat with the least efficiency. On the other hand Malkus proposed that convection would maximize the efficiency of energy transport. A new stellar envelope code is developed to first examine this other extreme, wherein rising turbulent elements transport heat with the greatest possible efficiency. This other extreme model differs from MLT by providing a small reduction in the upper convection zone temperatures but greatly diminished turbulent velocities below the top few hundred kilometers. Using the findings of deep atmospheric models with the Navier-Stokes equation allows the calculation of an intermediate solar envelope model. Consideration is given to solar observations, including recent helioseismology, to examine the position of the solar envelope compared with the envelope models.

Schatten, Kenneth H.↗
NASA NTRS1991

Energy coupling in catastrophic collisions

The prediction of events leading to the catastrophic collisions and disruption of solar system bodies is fraught with the same difficulties as are other theories of impact events; since one simply cannot perform experiments in the regime of interest. In the catastrophic collisions of asteroids that regime involves bodies of a few tons to hundred of kilometers in diameter, and velocities of several kilometers pre second. For hundred kilometer bodies, gravitational stresses dominate material fracture strengths, but those gravitational stresses are essentially absent for laboratory experiments. Only numerical simulations using hydrocodes can in principle analyze the true problems, but they have their own major uncertainties about the correctness of the physical models and properties. The question of the measure of the impactor and its energy coupling is investigated using numerical code calculations. The material model was that of a generic silicate rock, including high pressure melt and vapor phases, and includes material nonlinearity and dissipation via a Mie-Gruniesen model. A series of calculations with various size ratios and impact velocities are reported.

Holsapple, K. A.↗
NASA NTRS1992

Three-dimensional MHD modeling of magnetotail dynamics for different polytropic indices

The paper analyzes the 3D resistive tearing instability in a magnetotail configuration by means of a resistive MHD code for various forms of the energy equation, including ohmic and compressional heating. A complete derivation of the set of MHD equations from a kinetic formalism is presented, with all simplifying assumptions noted. Emphasis is placed on the changes in magnetotail evolution caused by a variation of the assumed value of the polytropic exponent gamma, representing the ratio of specific heat in the energy equation which closes the set of MHD equations. Within the gamma range of 2/3-4, a tearing mode is found to develop in all cases. Certain properties, such as the formation and ejection of a plasmoid, associated with the formation of the reversal region of the north-south component of the magnetic field, and strong earthward and tailward flow fields prevail as qualitative features independent of the value of gamma chosen. Ohmic dissipation is found to be the major source of the heating, while adiabatic or nonadiabatic compression may contribute primarily in the closed field region earthward of the reconnection site.

Hesse, Michael↗