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 973 records · Page 54

Relativistic modeling capabilities in PERSEUS extended MHD simulation code for HED plasmas

We discuss the incorporation of relativistic modeling capabilities into the PERSEUS extended MHD simulation code for high-energy-density (HED) plasmas, and present the latest hybrid X-pinch simulation results. The use of fully relativistic equations enables the model to remain self-consistent in simulations of such relativistic phenomena as X-pinches and laser-plasma interactions. By suitable formulation of the relativistic generalized Ohm’s law as an evolution equation, we have reduced the recovery of primitive variables, a major technical challenge in relativistic codes, to a straightforward algebraic computation. Our code recovers expected results in the non-relativistic limit, and reveals new physics in the modeling of electron beam acceleration following an X-pinch. Through the use of a relaxation scheme, relativistic PERSEUS is able to handle nine orders of magnitude in density variation, making it the first fluid code, to our knowledge, that can simulate relativistic HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

phq: A Fortran code to compute phonon quasiparticle properties and dispersions

We here introduce a Fortran code that computes anharmonic free energy of solids from first-principles based on our phonon quasiparticle approach. In this code, phonon quasiparticle properties, i.e., renormalized phonon frequencies and lifetimes, are extracted from mode-projected velocity auto-correlation functions (VAF) of modes sampled by molecular dynamics (MD) trajectories. Using renormalized frequencies as input, the code next constructs an effective harmonic force constant matrix to calculate anharmonic phonon dispersions over the whole Brillouin zone and thus the anharmonic free energy in the thermodynamic limit ( N → ∞ ). A detailed description of all the input parameters and the subroutines is provided as well. We illustrate the use of this code to compute ab initio temperature-dependent anharmonic phonons of Si in the diamond structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Program optimizations: The interplay between power, performance, and energy

Practical considerations for future supercomputer designs will impose limits on both instantaneous power consumption and total energy consumption. Working within these constraints while providing the maximum possible performance, application developers will need to optimize their code for speed alongside power and energy concerns. This paper analyzes the effectiveness of several code optimizations including loop fusion, data structure transformations, and global allocations. A per component measurement and analysis of different architectures is performed, enabling the examination of code optimizations on different compute subsystems. Using an explicit hydrodynamics proxy application from the U.S. Department of Energy, LULESH, we show how code optimizations impact different computational phases of the simulation. This provides insight for simulation developers into the best optimizations to use during particular simulation compute phases when optimizing code for future supercomputing platforms. Here, we examine and contrast both x86 and Blue Gene architectures with respect to these optimizations.

97 MATHEMATICS AND COMPUTING↗

The magnetotail boundary and energy transfer processes

A particle code is used to simulate the magnetopause region in the high latitude geomagnetic tail in which the magnetic field undergoes a significant increase in going from the magnetosheath to the magnetotail lobe. The simulation indicates that plasma can flow from the magnetosheath to the lobe, which is accompanied by a drop in pressure and density. In the earth's inertial frame, the particles do work against the convection electric field. Hence the magnetopause region serves as a dynamo. The simulation also shows that the width of the transition region increases with time. In the earth's inertial frame this is seen as an expansion of the magnetopause thickness in the antisunward direction.

Swift, D. W.↗

Recent Progress in the Development of a Multi-Layer Green's Function Code for Ion Beam Transport

To meet the challenge of future deep space programs, an accurate and efficient engineering code for analyzing the shielding requirements against high-energy galactic heavy radiation is needed. To address this need, a new Green's function code capable of simulating high charge and energy ions with either laboratory or space boundary conditions is currently under development. The computational model consists of combinations of physical perturbation expansions based on the scales of atomic interaction, multiple scattering, and nuclear reactive processes with use of the Neumann-asymptotic expansions with non-perturbative corrections. The code contains energy loss due to straggling, nuclear attenuation, nuclear fragmentation with energy dispersion and downshifts. Previous reports show that the new code accurately models the transport of ion beams through a single slab of material. Current research efforts are focused on enabling the code to handle multiple layers of material and the present paper reports on progress made towards that end.

Tweed, John↗

Replacement of Legacy Analytical Codes at the Advanced Test Reactor

For each operating cycle of the Advanced Test Reactor (ATR) at Idaho National Laboratory, a Core Safety Assurance Package (CSAP) is necessary to demonstrate compliance with the safety basis approved by the United State Department of Energy (DOE). Certain computer codes are used in CSAP development, most of them developed in-house. This work describes replacement of a large set of these codes and updates previous work. Replacement of legacy codes is necessary due to computer hardware failure but also has generally improved user-friendliness.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Replacement of Legacy Analytical Codes at the Advanced Test Reactor

For each operating cycle of the Advanced Test Reactor (ATR) at Idaho National Laboratory, a Core Safety Assurance Package (CSAP) is necessary to demonstrate compliance with the safety basis approved by the United State Department of Energy (DOE). Certain computer codes are used in CSAP development, most of them developed in-house. This work describes replacement of a large set of these codes and updates previous work. Replacement of legacy codes is necessary due to computer hardware failure but also has generally improved user-friendliness.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

ARRA Energy Storage Demonstration Projects: Lessons Learned and Recommendations

The American Recovery and Reinvestment Act (ARRA) of 2009 (Recovery Act) provided funding for 16 energy storage demonstration projects. The projects ranged in scope from feasibility studies and technology demonstrations to full-scale, operational energy storage plants. This investment had a significant positive impact on the grid-connected energy storage industry. The goal of this report is to summarize the lessons learned from the ARRA projects, and to make recommendations for future Department of Energy (DOE) investments. Information for this report primarily came from three sources: a questionnaire and interview with each project team; DOE energy storage program peer review presentations; and DOE reports required as part of the ARRA project. Some lessons learned were common to many projects. Development of standards, codes and protocols specific to energy storage systems will mitigate uncertainty over code compliance, streamline permitting, and should be a priority (especially related to safety). Removal of regulatory barriers that preclude optimal operation of an energy storage system with multiple applications would immediately enable further deployment (e.g., FERC standards preclude the marketing department from reliability activities). Maturity of the approach to monetization varies substantially between applications with frequency regulation as an ancillary service leading other applications. Finally, developers focused on ramp mitigation and time shifting envision a reference plant. This reference plant would scale up from the current demonstration systems and would lead to the deployment of 50 MW-scale peaker plants.

25 ENERGY STORAGE↗

Uncertainty and Sensitivity Analyses of a Pebble Bed HTGR Loss of Cooling Event

The Very High Temperature Reactor Methods Development group at the Idaho National Laboratory identified the need for a defensible and systematic uncertainty and sensitivity approach in 2009. This paper summarizes the results of an uncertainty and sensitivity quantification investigation performed with the SUSA code, utilizing the International Atomic Energy Agency CRP 5 Pebble Bed Modular Reactor benchmark and the INL code suite PEBBED-THERMIX. Eight model input parameters were selected for inclusion in this study, and after the input parameters variations and probability density functions were specified, a total of 800 steady state and depressurized loss of forced cooling (DLOFC) transient PEBBED-THERMIX calculations were performed. The six data sets were statistically analyzed to determine the 5% and 95% DLOFC peak fuel temperature tolerance intervals with 95% confidence levels. It was found that the uncertainties in the decay heat and graphite thermal conductivities were the most significant contributors to the propagated DLOFC peak fuel temperature uncertainty. No significant differences were observed between the results of Simple Random Sampling (SRS) or Latin Hypercube Sampling (LHS) data sets, and use of uniform or normal input parameter distributions also did not lead to any significant differences between these data sets.

Strydom, Gerhard↗

A convective and radiative heat transfer analysis for the FIRE II forebody

A Navier-Stokes flowfield solution method (LAURA code) using finite-rate chemistry and two-temperature thermal nonequilibrium was used in combination with two nonequilibrium radiative heat transfer codes to calculate heating for the FIRE II vehicle. An axisymmetric model of the actual body shape was used. One radiative heating code (NEQAIR) was used in uncoupled fashion with the flowfield solver's energy equations, while the other code (LORAN) was used in both coupled and uncoupled variations. Several trajectory points ranging from highly nonequilibrium flow to near-equilibrium flow were used for a study of both convective and radiative heating over the vehicle. Considerable variation in radiative heating was seen at the extremes, while agreement was good in the intermediate trajectory points. Total heat transfer calculations gave good comparison until the peak heating trajectory points were encountered, and returned to good agreement for the last two equilibrium points.

Greendyke, Robert B.↗

Review of Hybrid (Deterministic/Monte Carlo) Radiation Transport Methods, Codes, and Applications at Oak Ridge National Laboratory

This paper provides a review of the hybrid (Monte Carlo/deterministic) radiation transport methods and codes used at the Oak Ridge National Laboratory and examples of their application for increasing the efficiency of real-world, fixed-source Monte Carlo analyses. The two principal hybrid methods are (1) Consistent Adjoint Driven Importance Sampling (CADIS) for optimization of a localized detector (tally) region (e.g., flux, dose, or reaction rate at a particular location) and (2) Forward Weighted CADIS (FW-CADIS) for optimizing distributions (e.g., mesh tallies over all or part of the problem space) or multiple localized detector regions (e.g., simultaneous optimization of two or more localized tally regions). The two methods have been implemented and automated in both the MAVRIC sequence of SCALE 6 and ADVANTG, a code that works with the MCNP code. As implemented, the methods utilize the results of approximate, fast-running 3-D discrete ordinates transport calculations (with the Denovo code) to generate consistent space- and energy-dependent source and transport (weight windows) biasing parameters. These methods and codes have been applied to many relevant and challenging problems, including calculations of PWR ex-core thermal detector response, dose rates throughout an entire PWR facility, site boundary dose from arrays of commercial spent fuel storage casks, radiation fields for criticality accident alarm system placement, and detector response for special nuclear material detection scenarios and nuclear well-logging tools. Substantial computational speed-ups, generally O(102-4), have been realized for all applications to date. This paper provides a brief review of the methods, their implementation, results of their application, and current development activities, as well as a considerable list of references for readers seeking more information about the methods and/or their applications.

97 MATHEMATICS AND COMPUTING↗

The Effect of Second-Order Hydrodynamics on a Floating Offshore Wind Turbine

The design of offshore floating wind turbines uses design codes that can simulate the entire coupled system behavior. At the present, most codes include only first-order hydrodynamics, which induce forces and motions varying with the same frequency as the incident waves. Effects due to second- and higher-order hydrodynamics are often ignored in the offshore industry, because the forces induced typically are smaller than the first-order forces. In this report, first- and second-order hydrodynamic analysis used in the offshore oil and gas industry is applied to two different wind turbine concepts--a spar and a tension leg platform.

17 WIND ENERGY↗

SU-E-T-796: Variation of Surface Photon Energy Spectra On Bone Heterogeneity and Beam Obliquity Between Flattened and Unflattened Beam

Purpose: This study investigates the spectra of surface photon energy and energy fluence in the bone heterogeneity and beam obliquity using flattened and unflattened photon beams. The spectra were calculated in a bone and water phantom using Monte Carlo simulation (the EGSnrc code). Methods: Spectra of energy, energy fluence and mean energy of the 6 MV flattened and unflattened photon beams (field size = 10 × 10 cm{sup 2}) produced by a Varian TrueBEAM linear accelerator were calculated at the surfaces of a bone and water phantom using Monte Carlo simulations. The spectral calculations were repeated with the beam angles turned from 0° to 15°, 30° and 45° in the phantoms. Results: It is found that the unflattened photon beams contained more photons in the low-energy range of 0 – 2 MeV than the flattened beams with a flattening filter. Compared to the water phantom, both the flattened and unflattened beams had slightly less photons in the energy range < 0.4 MeV when a bone layer of 1 cm is present under the phantom surface. This shows that the presence of the bone decreased the low-energy photons backscattered to the phantom surface. When the photon beams were rotated from 0° to 45°, the number of photon and mean photon energy increased with the beam angle. This is because both the flattened and unflattened beams became more hardened when the beam angle increased. With the bone heterogeneity, the mean energies of both photon beams increased correspondingly. This is due to the absorption of low-energy photons by the bone, resulting in more significant beam hardening. Conclusion: The photon spectral information is important in studies on the patient’s surface dose enhancement when using unflattened photon beams in radiotherapy.

60 APPLIED LIFE SCIENCES↗

MT71x: Multi-Temperature Library Based on ENDF/B-VII.1

The Nuclear Data Team has released a multitemperature transport library, MT71x, based upon ENDF/B-VII.1 with a few modifications as well as additional evaluations for a total of 427 isotope tables. The library was processed using NJOY2012.39 into 23 temperatures. MT71x consists of two sub-libraries; MT71xMG for multigroup energy representation data and MT71xCE for continuous energy representation data. These sub-libraries are suitable for deterministic transport and Monte Carlo transport applications, respectively. The SZAs used are the same for the two sub-libraries; that is, the same SZA can be used for both libraries. This makes comparisons between the two libraries and between deterministic and Monte Carlo codes straightforward. Both the multigroup energy and continuous energy libraries were verified and validated with our checking codes checkmg and checkace (multigroup and continuous energy, respectively) Then an expanded suite of tests was used for additional verification and, finally, verified using an extensive suite of critical benchmark models. We feel that this library is suitable for all calculations and is particularly useful for calculations sensitive to temperature effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Leveraging MARVEL and SPHERE to Demonstrate NEAMS Thermal Hydraulics Codes

The Single Primary Heat Extraction and Removal Emulator (SPHERE) and Microreactor Applications Research Validation and Evaluation (MARVEL) were utilized to demonstrate thermal hydraulics codes of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. A recent long-duration test of a sodium heat pipe conducted at the SPHERE facility was utilized for the validation of the heat pipe code Sockeye. The first 10 hours of operation featured a startup of the heat pipe from room temperature; this portion was modeled with Sockeye and compared to the experimental data. Measurements included several thermocouple measurements across the length of the heat pipe surface. Two of Sockeye’s heat pipe models were applied to experiment, as well as several models of the cooling at the condenser end of the heat pipe. Sockeye results gave a decent match to the experimental data, with notable discrepancy in the condenser section, where experimental data indicated a significant inactive length of the heat pipe. The MARVEL microreactor was modeled through a coupled simulation of System Analysis Module (SAM) and SubChannel Module (SCM). This multiscale coupled simulation was performed as an exercise to demonstrate the compatibility and applicability of these codes for thermal-hydraulic analysis of microreactors. The coupling between the SCM and SAM codes for multiscale modeling of MARVEL, was achieved through a domain overlapping approach. In this example a transient SAM simulation provided boundary conditions to SCM, and SCM provided pressure drop to SAM.

97 - MATHEMATICS AND COMPUTING↗

Charge generation by heavy ions in power MOSFETs, burnout space predictions, and dynamic SEB sensitivity

The transport, energy loss, and charge production of heavy ions in the sensitive regions of IRF 150 power MOSFETs are described. The dependence and variation of transport parameters with ion type and energy relative to the requirements for single event burnout in this part type are discussed. Test data taken with this power MOSFET are used together with analyses by means of a computer code of the ion energy loss and charge production in the device to establish criteria for burnout and parameters for space predictions. These parameters are then used in an application to predict burnout rates in a geostationary orbit for power converters operating in a dynamic mode. Comparisons of rates for different geometries in simulating SEU (single event upset) sensitive volumes are presented.

Stassinopoulos, E. G.↗