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,549 records · Page 86

Calculation and evaluation of energy levels and radiative rates for 2s{sup 2}2p{sup 6}-2s2p{sup 6}np and 2s{sup 2}2p{sup 6}-2s{sup 2}2p{sup 5}nd(4{<=}n{<=}20) transitions in Cu XX

Energy levels, oscillator strengths, and radiative rates have been calculated for 2s{sup 2}2p{sup 6}-2s2p{sup 6}np and 2s{sup 2}2p{sup 6}-2s{sup 2}2p{sup 5}nd transitions with 4{<=}n{<=}20 in Cu XX ions. In the calculations, the Flexible Atomic Code has been used, and systematic analysis and evaluation of the accuracy of the calculated results for the two types of transitions are also given with different indicators. It is shown that the results for the transitions with a 2s hole are greatly improved by using the present method comparing to those from RELAC code. Furthermore, our results show that the energy levels of the 2s2p{sup 6}(n+1)l and 2s{sup 2}2p{sup 5}(n+2)l configurations closely interact with those of the 2s2p{sup 6}np series and it is essential to include the 2s{sup 2}2p{sup 5}(n+1)l configurations in the calculations for the 2s{sup 2}2p{sup 5}nd series.

74 ATOMIC AND MOLECULAR PHYSICS↗

Modeling, Performance Assessment, and Nodal Data Analysis of TRISO-Fueled Systems with Shift

This technical report documents several enhancements to the Shift Monte Carlo (MC) code under the US Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in fiscal year (FY) 2022. Performance enhancements were added to Shift specifically for tristructural isotropic (TRISO)–fueled reactor systems and guided based on performance analysis in FY 2021. For the pebble performance model developed in previous studies, the runtime improved by ~ 91× compared to the original model and ~ 2× compared to the user-optimized model. Compared to Serpent, Shift is ~ 3× slower if Serpent delta-tracking is enabled but ~ 2× faster when delta-tracking is disabled. The multigroup cross section generation was improved through simplifying tally input definitions, porting several post-processing tally operations from Python scripts into the Shift code base, and accounting for production reactions in the scattering multiplicity. Progress was also made on two emerging capabilities: (1) the development of Titan (a Shift reactor physics user interface) and (2) initial investigation into path-length tallies for computing multigroup scattering matrices.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Calculation of electron and isotopes dose point kernels with fluka Monte Carlo code for dosimetry in nuclear medicine therapy

Purpose: The calculation of patient-specific dose distribution can be achieved by Monte Carlo simulations or by analytical methods. In this study, fluka Monte Carlo code has been considered for use in nuclear medicine dosimetry. Up to now, fluka has mainly been dedicated to other fields, namely high energy physics, radiation protection, and hadrontherapy. When first employing a Monte Carlo code for nuclear medicine dosimetry, its results concerning electron transport at energies typical of nuclear medicine applications need to be verified. This is commonly achieved by means of calculation of a representative parameter and comparison with reference data. Dose point kernel (DPK), quantifying the energy deposition all around a point isotropic source, is often the one. Methods: fluka DPKs have been calculated in both water and compact bone for monoenergetic electrons (10{sup -3} MeV) and for beta emitting isotopes commonly used for therapy ({sup 89}Sr, {sup 90}Y, {sup 131}I, {sup 153}Sm, {sup 177}Lu, {sup 186}Re, and {sup 188}Re). Point isotropic sources have been simulated at the center of a water (bone) sphere, and deposed energy has been tallied in concentric shells. fluka outcomes have been compared to penelope v.2008 results, calculated in this study as well. Moreover, in case of monoenergetic electrons in water, comparison with the data from the literature (etran, geant4, mcnpx) has been done. Maximum percentage differences within 0.8{center_dot}R{sub CSDA} and 0.9{center_dot}R{sub CSDA} for monoenergetic electrons (R{sub CSDA} being the continuous slowing down approximation range) and within 0.8{center_dot}X{sub 90} and 0.9{center_dot}X{sub 90} for isotopes (X{sub 90} being the radius of the sphere in which 90% of the emitted energy is absorbed) have been computed, together with the average percentage difference within 0.9{center_dot}R{sub CSDA} and 0.9{center_dot}X{sub 90} for electrons and isotopes, respectively. Results: Concerning monoenergetic electrons, within 0.8{center_dot}R{sub CSDA} (where 90%-97% of the particle energy is deposed), fluka and penelope agree mostly within 7%, except for 10 and 20 keV electrons (12% in water, 8.3% in bone). The discrepancies between fluka and the other codes are of the same order of magnitude than those observed when comparing the other codes among them, which can be referred to the different simulation algorithms. When considering the beta spectra, discrepancies notably reduce: within 0.9{center_dot}X{sub 90}, fluka and penelope differ for less than 1% in water and less than 2% in bone with any of the isotopes here considered. Complete data of fluka DPKs are given as Supplementary Material as a tool to perform dosimetry by analytical point kernel convolution. Conclusions: fluka provides reliable results when transporting electrons in the low energy range, proving to be an adequate tool for nuclear medicine dosimetry.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Delayed-Gamma Energy Biasing with Exact Energy Sampling in MCNP 6.2.0

Characterization of delayed-particle signatures is of the utmost importance for detecting, identifying, and quantifying special nuclear material (SNM). Various non-destructive assay (NDA) techniques used in nuclear material safeguards rely on characteristic gamma-ray or neutron emissions to infer properties of a material. The all-particle, all-energy Monte-Carlo radiation transport code MCNP6 is capable of simulating delayed-particle signatures from SNM in addition to modeling the transport of emitted particles to a detector or other instrument. In some cases, NDA techniques require the measurement of low probability emissions from a given material. In physical scenarios, this problem can be bypassed by longer count-times or may be simply irrelevant due to the mass, and therefore the total activity, of the material being measured. However, when applying the Monte-Carlo method to such scenarios, sampling of low probability emissions requires a large number of histories to converge to a statistically acceptable solution. To address this issue, the delayed-neutron and gamma energy biasing feature was introduced in MCNP6 by H. Armstrong et al. in 2013 which allows the user to set up energy windows in which delayed-particle emission sampling can be adjusted. For delayed-gamma sampling, this biasing method was based on the bin-wise structure of gamma emission data for both the multi-bin gamma and line-emission sampling modes. Recent improvements have been made to the line-emission sampling algorithm which upgrade the previous mini-bin sampling treatment to an exact energy-line treatment. This upgrade was not compatible with the biasing method implemented by H. Armstrong et al. due to the sampled data no longer having a bin-wise structure in line-emission mode. For this reason the delayed-gamma energy biasing feature has been rewritten to support the exact line-sampling treatment. This paper provides a description of the delayed-gamma energy biasing (DGEB) feature in MCNP6 as well as results and discussion for a gamma-ray spectroscopy example used to demonstrate the biasing utility. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Multi-variable, Multi-value Hydrodynamic Modeling Approach to Support Marine Energy and Coastal Resilience Applications

Funded by the Marine Energy Seedling Laboratory Program, the study reported on here aimed to develop a multi-variable, multi-value hydrodynamic modeling approach based on the Finite-Volume Community Ocean Model (FVCOM) to support a wide range of marine renewable energy and coastal resilience applications. We migrated the marine hydrokinetic energy module to simulate tidal energy extraction by tidal turbines under wave-current interactions. We explored the particle tracking module in FVCOM and applied it to Sequim Bay, Washington, to illustrate the surface current field and mimic the pollutant transport near the surface. A tracer-water age module was added to FVCOM to examine water exchange through the entrance. Lastly, a user-specified vertical coordinate option was added to FVCOM to improve simulation of the effect of floating turbines on surface currents and scalar transport. The above-described model development has been tested using a simplified model domain and demonstrated in real-world coastal systems. The outcome of this project will provide a useful modeling code and approach for marine renewable energy and coastal resilience research and applications.

16 TIDAL AND WAVE POWER↗

Dissipation measures in weakly collisional plasmas

ABSTRACT The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell vpic, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov–Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

A new time dependent density functional algorithm for large systems and plasmons in metal clusters

A new algorithm to solve the Time Dependent Density Functional Theory (TDDFT) equations in the space of the density fitting auxiliary basis set has been developed and implemented. The method extracts the spectrum from the imaginary part of the polarizability at any given photon energy, avoiding the bottleneck of Davidson diagonalization. The original idea which made the present scheme very efficient consists in the simplification of the double sum over occupied-virtual pairs in the definition of the dielectric susceptibility, allowing an easy calculation of such matrix as a linear combination of constant matrices with photon energy dependent coefficients. The method has been applied to very different systems in nature and size (from H{sub 2} to [Au{sub 147}]{sup −}). In all cases, the maximum deviations found for the excitation energies with respect to the Amsterdam density functional code are below 0.2 eV. The new algorithm has the merit not only to calculate the spectrum at whichever photon energy but also to allow a deep analysis of the results, in terms of transition contribution maps, Jacob plasmon scaling factor, and induced density analysis, which have been all implemented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distributed Energy Resource Optimization Using a Software as Service (SaaS) Approach at the University of California, Davis Campus

Together with OSIsoft LLC as its private sector partner and matching sponsor, the Lawrence Berkeley National Laboratory (Berkeley Lab) won an FY09 Technology Commercialization Fund (TCF) grant from the U.S. Department of Energy. The goal of the project is to commercialize Berkeley Lab's optimizing program, the Distributed Energy Resources Customer Adoption Model (DER-CAM) using a software as a service (SaaS) model with OSIsoft as its first non-scientific user. OSIsoft could in turn provide optimization capability to its software clients. In this way, energy efficiency and/or carbon minimizing strategies could be made readily available to commercial and industrial facilities. Specialized versions of DER-CAM dedicated to solving OSIsoft's customer problems have been set up on a server at Berkeley Lab. The objective of DER-CAM is to minimize the cost of technology adoption and operation or carbon emissions, or combinations thereof. DER-CAM determines which technologies should be installed and operated based on specific site load, price information, and performance data for available equipment options. An established user of OSIsoft's PI software suite, the University of California, Davis (UCD), was selected as a demonstration site for this project. UCD's participation in the project is driven by its motivation to reduce its carbon emissions. The campus currently buys electricity economically through the Western Area Power Administration (WAPA). The campus does not therefore face compelling cost incentives to improve the efficiency of its operations, but is nonetheless motivated to lower the carbon footprint of its buildings. Berkeley Lab attempted to demonstrate a scenario wherein UCD is forced to purchase electricity on a standard time-of-use tariff from Pacific Gas and Electric (PG&E), which is a concern to Facilities staff. Additionally, DER-CAM has been set up to consider the variability of carbon emissions throughout the day and seasons. Two distinct analyses of value to UCD are possible using this approach. First, optimal investment choices for buildings under the two alternative objectives can be derived. Second, a week-ahead building operations forecaster has been written that executes DER-CAM to find an optimal operating schedule for buildings given their expected building energy services requirements, electricity prices, and local weather. As part of its matching contribution, OSIsoft provided a full implementation of PI and a server to install it on at Berkeley Lab. Using the PItoPI protocol, this gives Berkeley Lab researchers direct access to UCD's PI data base. However, this arrangement is in itself inadequate for performing optimizations. Additional data not included in UCD's PI database would be needed and the campus was not able to provide this information. This report details the process, results, and lessons learned of this commercialization project.

25 ENERGY STORAGE↗

Automated Comparison of Building Energy Simulation Engines [Slides]

This presentation describes the BEopt comparative test suite, which is a tool that facilitates the automated comparison of building energy simulation engines. It also demonstrates how the test suite is improving the accuracy of building energy simulation programs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Review of Life Cycle Cost Analysis Tools

Life cycle costing (LCC) is a vital aspect of decision-making in building retrofitting. It allows stakeholders to thoroughly assess all expenses involved throughout the lifespan of these retrofits, covering initial costs; ongoing operation, maintenance, and repair costs; and disposal costs of the system. The Federal Energy Management Program of the US Department of Energy has established clear guidelines for conducting LCC analysis, particularly for investments aimed at energy and water conservation, as well as renewable energy projects. These guidelines are detailed in the Code of Federal Regulations, 10 CFR 436, Subpart A, which outlines the methodology and procedure for LCC analysis (Department of Energy 1990). The National Institute of Standards and Technology (NIST) has developed the NIST Handbook 135 (Joshua Kneifel 2022) to standardize the process of LCC analysis within the building industry. Whereas the handbook is tailored to the needs of federal agencies, its principles and methodologies can be valuable for other organizations and industries interested in conducting life cycle costing analyses for their building projects.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A NUMERICAL TREATMENT OF ANISOTROPIC RADIATION FIELDS COUPLED WITH RELATIVISTIC RESISTIVE MAGNETOFLUIDS

We develop a numerical scheme for solving fully special relativistic, resistive radiation magnetohydrodynamics. Our code guarantees conservation of total mass, momentum, and energy. The radiation energy density and the radiation flux are consistently updated using the M-1 closure method, which can resolve an anisotropic radiation field, in contrast to the Eddington approximation, as well as the flux-limited diffusion approximation. For the resistive part, we adopt a simple form of Ohm's law. The advection terms are explicitly solved with an approximate Riemann solver, mainly the Harten-Lax-van Leer scheme; the HLLC and HLLD schemes are also solved for some tests. The source terms, which describe the gas-radiation interaction and the magnetic energy dissipation, are implicitly integrated, relaxing the Courant-Friedrichs-Lewy condition even in an optically thick regime or a large magnetic Reynolds number regime. Although we need to invert 4 Multiplication-Sign 4 matrices (for the gas-radiation interaction) and 3 Multiplication-Sign 3 matrices (for the magnetic energy dissipation) at each grid point for implicit integration, they are obtained analytically without preventing massive parallel computing. We show that our code gives reasonable outcomes in numerical tests for ideal magnetohydrodynamics, propagating radiation, and radiation hydrodynamics. We also applied our resistive code to the relativistic Petschek-type magnetic reconnection, revealing the reduction of the reconnection rate via radiation drag.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

THE AGORA HIGH-RESOLUTION GALAXY SIMULATIONS COMPARISON PROJECT. II. ISOLATED DISK TEST

Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt–Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of ∼3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

79 ASTRONOMY AND ASTROPHYSICS↗

Spatially distributed scintillator arrays for diagnosing runaway electron transport and energy behavior in tokamaks

We present details of a new bismuth germanate [Bi{sub 4}Ge{sub 3}O{sub 12} (BGO)] scintillator array used to diagnose the transport and energy behavior of runaway electrons (REs) in DIII-D. BGO exhibits important properties for these compact detectors including high light yield which sufficiently excites photodiode detectors (8500 photons/MeV), high density and atomic numbers of constituent materials which maximizes sensitivity, and relative neutron blindness which minimizes complications in data interpretation. The detectors observe primarily hard x-ray radiation emitted in a forward beamed pattern by RE when they strike first wall materials or bulk ions and neutrals in the plasma, although we also address photoneutron signals. The arrangement of the array enables time resolved location of x-ray emission and associated asymmetries which help identify instabilities and confinement properties of RE. By shielding a subset of detectors with different thicknesses of lead, and with interpretative support of the code EGSNRC, we also measure RE energy, although due to the often distributed nature of RE strike points and the forward beamed character of emitted hard x-rays, we restrict interpretation as a lower bound for RE energy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Energy Production Demonstrator and Material Testing Station optimization for Megawatt proton beams

A simulation study of the Energy Production Demonstrator (EPD) concept – a solid heavy metal target irradiated by GeV-range intense proton beams and producing more energy than consuming – is carried out. Neutron production, fission, energy deposition, energy gain, testing volume and helium production are simulated using the MARS15 code for tungsten, thorium, and natural uranium targets in the proton energy range 0.5–120 GeV. This study shows that the proton energy range of 2–4 GeV is optimal for both a nat U EPD and the tungsten-based testing station for proton accelerator facilities. Simulation-based conservative estimates not including breeding and fission of plutonium suggest the proton beam current sufficient to produce 1 GW of thermal output power with the nat U EPD while supplying a relatively small fraction of that power to operate the accelerator. The thermal analysis has been performed and has shown that the EPD with a parallel proton beam feeding the core has a potential problem due to a possible core meltdown. In conclusion, a scheme has been proposed with a beam steering on the outer surface of the target; its thermal analysis indicates that slicing the target in shorter parts and maintaining a temperature of 25 °C on their surfaces would to avoid the target meltdown.

43 PARTICLE ACCELERATORS↗

Saturation of Alfvén modes in tokamaks

Here, the growth of Alfvén modes driven unstable by a distribution of high energy particles up to saturation is investigated with a guiding center code, using numerical eigenfunctions produced by linear theory and a numerical high energy particle distribution, in order to make detailed comparison with experiment and with models for saturation amplitudes and the modification of beam profiles. Two innovations are introduced. First, a very noise free means of obtaining the mode-particle energy and momentum transfer is introduced, and secondly, a spline representation of the actual beam particle distribution is used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Saturation of Alfven modes in tokamaks

Growth of Alfven modes driven unstable by a distribution of high energy particles up to saturation is investigated with a guiding center code, using numerical eigenfunctions produced by linear theory and a numerical high energy particle distribution, in order to make detailed comparison with experiment and with models for saturation amplitudes and the modification of beam profiles. Two innovations are introduced. First, a very noise free means of obtaining the mode-particle energy and momentum transfer is introduced, and secondly, a spline representation of the actual beam particle distribution is used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AFIT and LANL Collaboration Research Proposals [Slides]

MCNP is a general-purpose, continuous-energy, generalized-geometry, time-dependent, Monte Carlo radiation-transport code designed to track many particle types over broad ranges of energies. MCNP has the capability to track particles on unstructured polyhedral meshes that are embedded into the legacy constructive solid geometry (CSG).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling and Simulation of HVAC Faulty Operations and Performance Degradation due to Maintenance Issues

Almost half of the total energy used in the U.S. buildings is consumed by heating, ventilation and air conditionings (HVAC) according to EIA statistics. Among various driving factors to energy performance of building, operations and maintenance play a significant role. Many researches have been done to look at design efficiencies and operational controls for improving energy performance of buildings, but very few study the impacts of HVAC systems maintenance. Different practices of HVAC system maintenance can result in substantial differences in building energy use. If a piece of HVAC equipment is not well maintained, its performance will degrade. If sensors used for control purpose are not calibrated, not only building energy usage could be dramatically increased, but also mechanical systems may not be able to satisfy indoor thermal comfort. Properly maintained HVAC systems can operate efficiently, improve occupant comfort, and prolong equipment service life. In the paper, maintenance practices for HVAC systems are presented based on literature reviews and discussions with HVAC engineers, building operators, facility managers, and commissioning agents. We categorize the maintenance practices into three levels depending on the maintenance effort and coverage: 1) proactive, performance-monitored maintenance; 2) preventive, scheduled maintenance; and 3) reactive, unplanned or no maintenance. A sampled list of maintenance issues, including cooling tower fouling, boiler/chiller fouling, refrigerant over or under charge, temperature sensor offset, outdoor air damper leakage, outdoor air screen blockage, outdoor air damper stuck at fully open position, and dirty filters are investigated in this study using field survey data and detailed simulation models. The energy impacts of both individual maintenance issue and combined scenarios for an office building with central VAV systems and central plant were evaluated by EnergyPlus simulations using three approaches: 1) direct modeling with EnergyPlus, 2) using the energy management system feature of EnergyPlus, and 3) modifying EnergyPlus source code. The results demonstrated the importance of maintenance for HVAC systems on energy performance of buildings. The research is intended to provide a guideline to help practitioners and building operators to gain the knowledge of maintaining HVAC systems in efficient operations, and prioritize HVAC maintenance work plan. The paper also discusses challenges of modeling building maintenance issues using energy simulation programs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗