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 901 records · Page 50

Evaluation for 4S core nuclear design method through integration of benchmark data

The 4S is a sodium-cooled small fast reactor which is reflector-controlled for operation through core lifetime about 30 years. The nuclear design method has been selected to treat neutron leakage with high accuracy. It consists of a continuous-energy Monte Carlo code, discrete ordinate transport codes and JENDL-3.3. These two types of neutronic analysis codes are used for the design in a complementary manner. The accuracy of the codes has been evaluated by analysis of benchmark critical experiments and the experimental reactor data. The measured data used for the evaluation is critical experimental data of the FCA XXIII as a physics mockup assembly of the 4S core, FCA XVI, FCA XIX, ZPR, and data of experimental reactor JOYO MK-1. Evaluated characteristics are criticality, reflector reactivity worth, power distribution, absorber reactivity worth, and sodium void worth. A multi-component bias method was applied, especially to improve the accuracy of sodium void reactivity worth. As the result, it has been confirmed that the 4S core nuclear design method provides good accuracy, and typical bias factors and their uncertainties are determined. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Benchmarking ENDF/B-VII.1, JENDL-4.0 and JEFF-3.1

Three nuclear data libraries have been tested extensively using criticality safety benchmark calculations. The three libraries are the new release of the US library ENDF/B-VII.1 (2011), the new release of the Japanese library JENDL-4.0 (2011), and the OECD/NEA library JEFF-3.1 (2006). All calculations were performed with the continuous-energy Monte Carlo code MCNP (version 4C3, as well as version 6-beta1). Around 2000 benchmark cases from the International Handbook of Criticality Safety Benchmark Experiments (ICSBEP) were used. The results were analyzed per ICSBEP category, and per element. Overall, the three libraries show similar performance on most criticality safety benchmarks. The largest differences are probably caused by elements such as Be, C, Fe, Zr, W. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Calculation of neutron and gamma fluxes in support to the interpretation of measuring devices irradiated in the core periphery of the OSIRIS Material Testing Reactor

Technological irradiations carried out in material testing reactors (MTRs) are used to study the behavior of materials under irradiation conditions required by different types of nuclear power plants (NPPs). For MTRs, specific instrumentation is required for the experiment monitoring and for the characterization of irradiation conditions, in particular the flux of neutrons and photons. To measure neutron and photon flux in experimental locations, different sensors can be used, such as SPNDs (self-powered neutron detectors), SPGDs (self-powered gamma detectors) and ionization chambers. These sensors involve interactions producing ultimately a measurable electric current. Various sensors have been recently tested in the core periphery of the OSIRIS reactor (located at the CEA-Saclay center) in order to qualify their responses to the neutron and the photon flux. One of the key input data for this qualification is to have a relevant evaluation of neutron and gamma fluxes at the irradiation location. The objective of this work is to evaluate the neutron and the gamma flux in the core periphery of the OSIRIS reactor. With this intention, specific neutron-photonic three-dimensional calculations have been performed and are mainly based on the TRIPOLI-4{sup R} three-dimensional continuous-energy Monte Carlo code, developed by CEA (Saclay Center) and extensively validated against reactor dosimetry benchmarks. In the case of the OSIRIS reactor, TRIPOLI-4{sup R} code has been validated against experimental results based on neutron flux and nuclear heating measurements performed in ex-core and in-core experiments. In this work, simultaneous contribution of neutrons and gamma photons in the core periphery is considered using neutron-photon coupled transport calculations. Contributions of prompt and decay photons have been taken into account for the gamma flux calculation. Specific depletion codes are used upstream to provide the decay-gamma sources required by TRIPOLI-4 calculations. In this paper, the neutron-photonic calculation scheme used for the neutron-gamma flux evaluation in the OSIRIS core periphery will be first detailed. Then, a comparison of calculated and measured values of neutron and gamma fluxes will be introduced. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Coupled dynamics analysis of wind energy systems

A qualitative description of all key elements of a complete wind energy system computer analysis code is presented. The analysis system addresses the coupled dynamics characteristics of wind energy systems, including the interactions of the rotor, tower, nacelle, power train, control system, and electrical network. The coupled dynamics are analyzed in both the frequency and time domain to provide the basic motions and loads data required for design, performance verification and operations analysis activities. Elements of the coupled analysis code were used to design and analyze candidate rotor articulation concepts. Fundamental results and conclusions derived from these studies are presented.

Hoffman, J. A.↗

A new moment-based general-relativistic neutrino-radiation transport code: Methods and first applications to neutron star mergers

ABSTRACT We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in υ/c, retaining all non-linear neutrino–matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino–matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim }70\,$ ms. The latter is the longest dynamical space-time, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the post-merger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries, and equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗

An Update on NiCE Support for BISON

The Nuclear Energy Advanced Modeling and Simulation program (NEAMS) from the Department of Energy s Office of Nuclear Energy has funded the development of a modeling and simulation workflow environment to support the various codes in its nuclear energy scientific computing toolkit. This NEAMS Integrated Computational Environment (NiCE) provides extensible tools and services that enable efficient code execution, input generation, pre-processing visualizations, and post-simulation data analysis and visualization for a large portion of the NEAMS Toolkit. A strong focus for the NiCE development team throughout FY 2015 has been support for the Multiphysics Object Oriented Simulation Environment (MOOSE) and the NEAMS nuclear fuel performance modeling application built on that environment, BISON. There is a strong desire in the program to enable and facilitate the use of BISON throughout nuclear energy research and industry. A primary result of this desire is the need for strong support for BISON in NiCE. This report will detail improvements to NiCE support for BISON. We will present a new and improved interface for interacting with BISON simulations in a variety of ways: (1) improved input model generation, (2) embedded mesh and solution data visualizations, and (3) local and remote BISON simulation launch. We will also show how NiCE has been extended to provide support for BISON code development.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

(γ, 2n)-Reaction cross-section calculations of several even-even lanthanide nuclei using different level density models

There are several level density models that can be used to predict photo-neutron cross sections. Some of them are Constant Temperature + Fermi Gas Model (CTFGM), Back-Shifted Fermi Gas Model (BSFM), Generalized Superfluid Model (GSM), Hartree-Fock-Bogoliubov microscopic Model (HFBM). In this study, the theoretical photo-neutron cross sections produced by (γ, 2n) reactions for several eveneven lanthanide nuclei such as {sup 140,142}Ce, {sup 142,144,146,148,150}Nd, {sup 144,148,150,152,154}Sm, and {sup 160}Gd have been calculated on the different level density models as mentioned above by using TALYS 1.6 and EMPIRE 3.1 computer codes for incident photon energies up to 30 MeV. The obtained results have been compared with each other and available experimental data existing in the EXFOR database. Generally, at least one level density model cross-section calculations are in agreement with the experimental results for all reactions except {sup 144}Sm(γ, 2n){sup 142}Sm along the incident photon energy, TALYS 1.6 BSFM option for the level density model cross-section calculations can be chosen if the experimental data are not available or are improbable to be produced due to the experimental difficulty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

WEC-SIM Phase 1 Validation Testing -- Numerical Modeling of Experiments: Preprint

The Wave Energy Converter Simulator (WEC-Sim) is an open-source code jointly developed by Sandia National Laboratories and the National Renewable Energy Laboratory. It is used to model wave energy converters subjected to operational and extreme waves. In order for the WEC-Sim code to be beneficial to the wave energy community, code verification and physical model validation is necessary. This paper describes numerical modeling of the wave tank testing for the 1:33-scale experimental testing of the floating oscillating surge wave energy converter. The comparison between WEC-Sim and the Phase 1 experimental data set serves as code validation. This paper is a follow-up to the WEC-Sim paper on experimental testing, and describes the WEC-Sim numerical simulations for the floating oscillating surge wave energy converter.

16 TIDAL AND WAVE POWER↗

A Binary-Encounter-Bethe Approach to Simulate DNA Damage by the Direct Effect

The DNA damage is of crucial importance in the understanding of the effects of ionizing radiation. The main mechanisms of DNA damage are by the direct effect of radiation (e.g. direct ionization) and by indirect effect (e.g. damage by.OH radicals created by the radiolysis of water). Despite years of research in this area, many questions on the formation of DNA damage remains. To refine existing DNA damage models, an approach based on the Binary‐Encounter‐Bethe (BEB) model was developed[1]. This model calculates differential cross sections for ionization of the molecular orbitals of the DNA bases, sugars and phosphates using the electron binding energy, the mean kinetic energy and the occupancy number of the orbital. This cross section has an analytic form which is quite convenient to use and allows the sampling of the energy loss occurring during an ionization event. To simulate the radiation track structure, the code RITRACKS developed at the NASA Johnson Space Center is used[2]. This code calculates all the energy deposition events and the formation of the radiolytic species by the ion and the secondary electrons as well. We have also developed a technique to use the integrated BEB cross section for the bases, sugar and phosphates in the radiation transport code RITRACKS. These techniques should allow the simulation of DNA damage by ionizing radiation, and understanding of the formation of double‐strand breaks caused by clustered damage in different conditions.

Plante, Ianik↗

Los Alamos RAGE Simulations of the HAIV Mission Concept

The mitigation of potentially hazardous objects (PHOs) can be accomplished by a variety of methods including kinetic impactors, gravity tractors and several nuclear explosion options. Depending on the available lead time prior to Earth impact, non- nuclear options can be very effective at altering a PHOs orbit. However if the warning time is short nuclear options are generally deemed most effective at mitigating the hazard. The NIAC mission concept for a nuclear mission has been presented at several meetings, including the last PDC (2013).We use the adaptive mesh hydrocode RAGE to perform detailed simulations of this Hypervelocity Asteroid Intercept Vehicle (HAIV) mission concept. We use the RAGE code to simulate the crater formation by the kinetic impactor as well as the explosion and energy coupling from the follower nuclear explosive device (NED) timed to detonate below the original surface to enhance the energy coupling. The RAGE code has been well validated for a wide variety of applications. A parametric study will be shown of the energy and momentum transfer to the target 100 m diameter object: 1) the HAIV mission as planned; 2) a surface explosion and 3) a subsurface (contained) explosion; both 2) and 3) use the same source energy as 1).Preliminary RAGE simulations show that the kinetic impactor will carve out a surface crater on the object and the subsequent NED explosion at the bottom of the crater transfers energy and momentum to the target effectively moving it off its Earth crossing orbit. Figure 1 shows the initial (simplified) RAGE 2D setup geometry for this study. Figure 2 shows the crater created by the kinetic impactor and Figure 3 shows the time sequence of the energy transfer to the target by the NED.

Asteroid Mitigation↗

You can hide but you have to run: Direct detection with vector mediators

We study direct detection in simplified models of Dark Matter (DM) in which interactions with Standard Model (SM) fermions are mediated by a heavy vector boson. We consider fully general, gauge-invariant couplings between the SM, the mediator and both scalar and fermion DM. We account for the evolution of the couplings between the energy scale of the mediator mass and the nuclear energy scale. This running arises from virtual effects of SM particles and its inclusion is not optional. We compare bounds on the mediator mass from direct detection experiments with and without accounting for the running. In some cases the inclusion of these effects changes the bounds by several orders of magnitude, as a consequence of operator mixing which generates new interactions at low energy. We also highlight the importance of these effects when translating LHC limits on the mediator mass into bounds on the direct detection cross section. For an axial-vector mediator, the running can alter the derived bounds on the spin-dependent DM-nucleon cross section by a factor of two or more. Lastly, we provide tools to facilitate the inclusion of these effects in future studies: general approximate expressions for the low energy couplings and a public code runDM to evolve the couplings between arbitrary energy scales.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

WEC-Sim (Wave Energy Converter - SIMulator)

WEC-Sim (Wave Energy Converter SIMulator) is a code developed by Sandia National Laboratories and the National Renewable Energy Laboratory to model wave energy converters (WECs) when they are subject to operational waves. The code is a time-domain modeling tool developed in MATLAB/Simulink using the multi-body dynamics solver SimMechanics. In WEC-Sim, WECs are modeled by connecting rigid bodies to one another with joint or constraint blocks from the WEC-Sim library. WEC-Sim is a publicly available, open-source code to model WECs.

Source record↗

WEC-Sim (Wave Energy Converter - SIMulator)

WEC-Sim (Wave Energy Converter SIMulator) is a code developed by Sandia National Laboratories and the National Renewable Energy Laboratory to model wave energy converters (WECs) when they are subject to operational waves. The code is a time-domain modeling tool developed in MATLAB/Simulink using the multi-body dynamics solver SimMechanics. In WEC-Sim, WECs are modeled by connecting rigid bodies to one another with joint or constraint blocks from the WEC-Sim library. WEC-Sim is a publicly available, open-source code to model WECs.

Source record↗

Computational Particle Dynamic Simulations on Multicore Processors (CPDMu) Final Report Phase I

Statement of Problem - Department of Energy has many legacy codes for simulation of computational particle dynamics and computational fluid dynamics applications that are designed to run on sequential processors and are not easily parallelized. Emerging high-performance computing architectures employ massively parallel multicore architectures (e.g., graphics processing units) to increase throughput. Parallelization of legacy simulation codes is a high priority, to achieve compatibility, efficiency, accuracy, and extensibility. General Statement of Solution - A legacy simulation application designed for implementation on mainly-sequential processors has been represented as a graph G. Mathematical transformations, applied to G, produce a graph representation {und G} for a high-performance architecture. Key computational and data movement kernels of the application were analyzed/optimized for parallel execution using the mapping G {yields} {und G}, which can be performed semi-automatically. This approach is widely applicable to many types of high-performance computing systems, such as graphics processing units or clusters comprised of nodes that contain one or more such units. Phase I Accomplishments - Phase I research decomposed/profiled computational particle dynamics simulation code for rocket fuel combustion into low and high computational cost regions (respectively, mainly sequential and mainly parallel kernels), with analysis of space and time complexity. Using the research team's expertise in algorithm-to-architecture mappings, the high-cost kernels were transformed, parallelized, and implemented on Nvidia Fermi GPUs. Measured speedups (GPU with respect to single-core CPU) were approximately 20-32X for realistic model parameters, without final optimization. Error analysis showed no loss of computational accuracy. Commercial Applications and Other Benefits - The proposed research will constitute a breakthrough in solution of problems related to efficient parallel computation of particle and fluid dynamics simulations. These problems occur throughout DOE, military and commercial sectors: the potential payoff is high. We plan to license or sell the solution to contractors for military and domestic applications such as disaster simulation (aerodynamic and hydrodynamic), Government agencies (hydrological and environmental simulations), and medical applications (e.g., in tomographic image reconstruction). Keywords - High-performance Computing, Graphic Processing Unit, Fluid/Particle Simulation. Summary for Members of Congress - Department of Energy has many simulation codes that must compute faster, to be effective. The Phase I research parallelized particle/fluid simulations for rocket combustion, for high-performance computing systems.

97 MATHEMATICS AND COMPUTING↗

Model for transport and reaction of defects and carriers within displacement cascades in gallium arsenide

A model is presented for recombination of charge carriers at evolving displacement damage in gallium arsenide, which includes clustering of the defects in atomic displacement cascades produced by neutron or ion irradiation. The carrier recombination model is based on an atomistic description of capture and emission of carriers by the defects with time evolution resulting from the migration and reaction of the defects. The physics and equations on which the model is based are presented, along with the details of the numerical methods used for their solution. The model uses a continuum description of diffusion, field-drift and reaction of carriers, and defects within a representative spherically symmetric cluster of defects. The initial radial defect profiles within the cluster were determined through pair-correlation-function analysis of the spatial distribution of defects obtained from the binary-collision code MARLOWE, using recoil energies for fission neutrons. Properties of the defects are discussed and values for their parameters are given, many of which were obtained from density functional theory. The model provides a basis for predicting the transient response of III-V heterojunction bipolar transistors to displacement damage from energetic particle irradiation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Diffusive Transport of Several Hundred keV Electrons in the Earth's Slot Region

Here, we investigate the gradual diffusion of energetic electrons from the inner edge of the outer radiation belt into the slot region. The Van Allen Probes observed slow inward diffusion and decay of ~200–600 keV electrons following the intense geomagnetic storm that occurred on 17 March 2013. During the 10 day nondisturbed period following the storm, the peak of electron fluxes gradually moved from L ~ 2.7 to L ~ 2.4, and the flux levels decreased by a factor of ~2–4 depending on the electron energy. We simulated the radial intrusion and decay of electrons using a three–dimensional diffusion code, which reproduced the energy–dependent transport of electrons from ~100 keV to 1 MeV in the slot region. At energies of 100–200 keV, the electrons experience fast transport across the slot region due to the dominance of radial diffusion; at energies of 200–600 keV, the electrons gradually diffuse and decay in the slot region due to the comparable rate of radial diffusion and pitch angle scattering by plasmaspheric hiss; at energies of E > 700 keV, the electrons stopped diffusing near the inner edge of outer radiation belt due to the dominant pitch angle scattering loss. In addition to plasmaspheric hiss, magnetosonic waves and VLF transmitters can cause the loss of high pitch angle electrons, relaxing the sharp “top–hat” shaped pitch angle distributions created by plasmaspheric hiss. Our simulation indicates the importance of balance between radial diffusion and loss through pitch angle scattering in forming the diffusive intrusion of energetic electrons across the slot region.

79 ASTRONOMY AND ASTROPHYSICS↗

Collaborative Research: Improved Efficiency and Coupling of the Radiation Code in the ACME Earth System Model. Final Report

This final report details all work performed on the project by both project partners. This project provided support to properly couple RTE+RRTMGP, a high-performance broadband radiation code, within DOE’s Energy Exascale Earth System Model (E3SM). RTE+RRTMGP is a successor to the RRTMG radiation code, which has been widely accepted for its speed and accuracy by the global modeling community, and has been in use in the NCAR CESM for many years and was implemented in the initial version of E3SM. However, the computational cost of RRTMG remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, issues that were rectified by the development of RTE+RRTMGP. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. One focus of the project was to enhance the code’s optimization on the limited number of emerging computing systems on which the model is expected be used, including Many Integrated Core (MIC) architectures and Graphics Processing Unit (GPU) hardware. We also developed additional capabilities for RTE+RRTMGP that E3SM scientists identified as important for the planned applications of the model. The result of our project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.

54 ENVIRONMENTAL SCIENCES↗

Neutron environments on the Martian surface

Radiation is a primary concern in the planning of a manned mission to Mars. Recent studies using NASA Langley Research Center's HZETRN space radiation transport code show that the low energy neutron fluence on the Martian surface is larger than previously expected. The upper atmosphere of Mars is exposed to a background radiation field made up of a large number of protons during a solar particle event and mixture of light and heavy ions caused by galactic cosmic rays at other times. In either case, these charged ions interact with the carbon and oxygen atoms of the Martian atmosphere through ionization and nuclear collisions producing secondary ions and neutrons which then interact with the atmospheric atoms in a similar manner. In the past, only these downward moving particles have been counted in evaluating the neutron energy spectrum on the surface. Recent enhancements in the HZETRN code allow for the additional evaluation of those neutrons created within the Martian regolith through the same types of nuclear reactions, which rise to the surface. New calculations using this improved HZETRN code show that these upward moving neutrons contribute significantly to the overall neutron spectrum for energies less than 10 MeV.

NASA Center LaRC↗