Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Energy Simulation”

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,693 records · Page 94

Nuclear recoil scintillation linearity of a high pressure 4 He gas detector

We investigate scintillation linearity of a commercial high pressure 4 He gas detector using monoenergetic 2.8 MeV neutrons from a deuterium-deuterium fusion neutron generator. The scintillation response of the detector was measured for a range of recoil energies between 83 keV and 626 keV by tagging neutrons scattering into fixed angles with a far-side organic scintillator detector. Detailed Monte Carlo simulations were compared to experimental data to determine the linearity of the detector response by comparing the scaling of the energy deposits in the simulations to the detector output. In this analysis, a linear scintillation response corresponds to a consistent value for the scaling factor between simulated energy deposits and experimental data for several different scattering angles. Here we demonstrate that the detector can be used to detect fast neutron interactions down to 83 keV recoil energies and can be used to characterize low-energy neutron sources, one of its potential applications.

47 OTHER INSTRUMENTATION↗

Parallel Algorithms and Software for Nuclear, Energy, and Environmental Applications Part I: Multiphysics Algorithms

There is a growing trend within energy and environmental simulation to consider tightly coupled solutions to multiphysics problems. This can be seen in nuclear reactor analysis where analysts are interested in coupled flow, heat transfer and neutronics, and in nuclear fuel performance simulation where analysts are interested in thermomechanics with contact coupled to species transport and chemistry. In energy and environmental applications, energy extraction involves geomechanics, flow through porous media and fractured formations, adding heat transport for enhanced oil recovery and geothermal applications, and adding reactive transport in the case of applications modeling the underground flow of contaminants. These more ambitious simulations usually motivate some level of parallel computing. Many of the physics coupling efforts to date utilize simple code coupling or first-order operator splitting, often referred to as loose coupling. While these approaches can produce answers, they usually leave questions of accuracy and stability unanswered. Additionally, the different physics often reside on distinct meshes and data are coupled via simple interpolation, again leaving open questions of stability and accuracy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Rotational-translational energy transfer in rarefied nonequilibrium flows

A new model for simulating the transfer of energy between the translational and rotational modes is derived for a homogeneous gas of diatomic molecules. The model has been developed specifically for use in discrete particle simulation methods where molecular motion and intermolecular collisions are treated at the molecular level. A temperature dependence is introduced which has been predicted by theory and observed in experiment. The new model is applied to the relaxation of rotational temperature, and is found to produce significant differences in comparison with the model normally employed at both high and low temperatures. Calculations have also been performed for a Mach 7 normal shock wave.

Boyd, Iain D.↗

Voltage-controlled low-energy switching of nanomagnets through Ruderman-Kittel-Kasuya-Yosida interactions for magnetoelectric device applications

In this article, we consider through simulation low-energy switching of nanomagnets via electrostatically gated inter-magnet Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions on the surface of three-dimensional topological insulators, for possible memory and nonvolatile logic applications. We model the possibility and dynamics of RKKY-based switching of one nanomagnet by coupling to one or more nanomagnets of set orientation. Potential applications to both memory and nonvolatile logic are illustrated. Sub-attojoule switching energies, far below conventional spin transfer torque (STT)-based memories and even below CMOS logic appear possible. Switching times on the order of a few nanoseconds, comparable to times for STT switching, are estimated for ferromagnetic nanomagnets, but the approach also appears compatible with the use of antiferromagnets which may allow for faster switching.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Chicago microclimate and building energy use data

These data comprise three elements: - High resolution, 90 m simulated weather data for 1 year at 15 min. intervals (with known gaps toward the end of each month). These files are in .csv format. - A mapping of individual buildings with individual IDs, their latitude/longitude location, and height. (Excel file) - Energy simulation output of these individual buildings, at 15 min. intervals for a whole year. (.json and other files)

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Chicago microclimate and building energy use data

These data comprise three elements: - High resolution, 90 m simulated weather data for 1 year at 15 min. intervals (with known gaps toward the end of each month). These files are in .csv format. - A mapping of individual buildings with individual IDs, their latitude/longitude location, and height. (Excel file) - Energy simulation output of these individual buildings, at 15 min. intervals for a whole year. (.json and other files)

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Shielding implications for secondary neutrons and photons produced within the patient during IMPT

Purpose: Intensity modulated proton therapy (IMPT) uses a combination of computer controlled spot scanning and spot-weight optimized planning to irradiate the tumor volume uniformly. In contrast to passive scattering systems, secondary neutrons and photons produced from inelastic proton interactions within the patient represent the major source of emitted radiation during IMPT delivery. Various published studies evaluated the shielding considerations for passive scattering systems but did not directly address secondary neutron production from IMPT and the ambient dose equivalent on surrounding occupational and nonoccupational work areas. Thus, the purpose of this study was to utilize Monte Carlo simulations to evaluate the energy and angular distributions of secondary neutrons and photons following inelastic proton interactions within a tissue-equivalent phantom for incident proton spot energies between 70 and 250 MeV.Methods: Monte Carlo simulation methods were used to calculate the ambient dose equivalent of secondary neutrons and photons produced from inelastic proton interactions in a tissue-equivalent phantom. The angular distribution of emitted neutrons and photons were scored as a function of incident proton energy throughout a spherical annulus at 1, 2, 3, 4, and 5 m from the phantom center. Appropriate dose equivalent conversion factors were applied to estimate the total ambient dose equivalent from secondary neutrons and photons.Results: A reference distance of 1 m from the center of the patient was used to evaluate the mean energy distribution of secondary neutrons and photons and the resulting ambient dose equivalent. For an incident proton spot energy of 250 MeV, the total ambient dose equivalent (3.6 Multiplication-Sign 10{sup -3} mSv per proton Gy) was greatest along the direction of the incident proton spot (0 Degree-Sign -10 Degree-Sign ) with a mean secondary neutron energy of 71.3 MeV. The dose equivalent decreased by a factor of 5 in the backward direction (170 Degree-Sign -180 Degree-Sign ) with a mean energy of 4.4 MeV. An 8 Multiplication-Sign 8 Multiplication-Sign 8 cm{sup 3} volumetric spot distribution (5 mm FWHM spot size, 4 mm spot spacing) optimized to produce a uniform dose distribution results in an ambient dose equivalent of 4.5 Multiplication-Sign 10{sup -2} mSv per proton Gy in the forward direction.Conclusions: This work evaluated the secondary neutron and photon emission due to monoenergetic proton spots between 70 and 250 MeV, incident on a tissue equivalent phantom. Example calculations were performed to estimate concrete shield thickness based upon appropriate workload and shielding design assumptions. Although lower than traditional passive scattered proton therapy systems, the ambient dose equivalent from secondary neutrons produced by the patient during IMPT can be significant relative to occupational and nonoccupational workers in the vicinity of the treatment vault. This work demonstrates that Monte Carlo simulations are useful as an initial planning tool for studying the impact of the treatment room and maze design on surrounding occupational and nonoccupational work areas.

61 RADIATION PROTECTION AND DOSIMETRY↗

Using Co-Optimized Machine Learned Manifolds for Modeling Chemically Reacting Flows

Chemically reacting flows play a key role in a wide range of engineered systems, from chemical and polymer processing to combustion-based energy conversion technologies. Simulations of these flows involve solving a coupled set of partial differential equations for mass, momentum, energy, and all relevant chemical species in the system. Chemical reaction pathways may be extremely complex and involve hundreds or more intermediate species, with reactions that occur over timescales varying by several orders of magnitude - presenting a significant numerical stiffness challenge. The combination of these factors makes simulation of chemically reacting flows vastly more expensive than nonreactive simulations, and often makes direct solution of the governing equations intractable. It is necessary to apply lower-fidelity models in place of the detailed governing equations in order to reduce computational cost to enable reacting flow simulation tools to be used in the engineering design process. Many of the models employed for this purpose are based on reducing the dimension of the thermochemical state, motivated by the observation that the observed thermochemical states in a system lie on a low-dimensional manifold in thermochemical state space. This behavior occurs due to the fast equilibration of certain reactive and transport processes, and physics-based manifold models rely on idealized assumptions about the balance of timescales and the way in which chemistry and transport are coupled. In this work, we apply a novel method for data-driven manifold-based modeling that can leverage data from high-fidelity reacting flow simulations to improve model accuracy in cases where the physics-based modeling assumptions break down. The approach is designed to be broadly applicable across chemically reacting flow systems but is applied here to turbulent combustion modeling.

machine learning↗

On The Origin Of High Energy Correlations in Gamma-ray Bursts

I investigate the origin of the observed correlation between a gamma-ray burst's {nu}F{sub {nu}} spectral peak E{sub pk} and its isotropic equivalent energy E{sub iso} through the use of a population synthesis code to model the prompt gamma-ray emission from GRBs. By using prescriptions for the distribution of prompt spectral parameters as well as the population's luminosity function and co-moving rate density, I generate a simulated population of GRBs and examine how bursts of varying spectral properties and redshift would appear to a gamma-ray detector here on Earth. I find that a strong observed correlation can be produced between the source frame Epk and Eiso for the detected population despite the existence of only a weak and broad correlation in the original simulated population. The energy dependance of a gamma-ray detector's flux-limited detection threshold acts to produce a correlation between the source frame E{sub pk} and E{sub iso} for low luminosity GRBs, producing the left boundary of the observed correlation. Conversely, very luminous GRBs are found at higher redshifts than their low luminosity counterparts due to the standard Malquest bias, causing bursts in the low E{sub pk}, high E{sub iso} regime to go undetected because their E{sub pk} values would be redshifted to energies at which most gamma-ray detectors become less sensitive. I argue that it is this previously unexamined effect which produces the right boundary of the observed correlation. Therefore, the origin of the observed correlation is a complex combination of the instrument's detection threshold, the intrinsic cutoff in the GRB luminosity function, and the broad range of redshifts over which GRBs are detected. Although the GRB model presented here is a very simplified representation of the complex nature of GRBs, these simulations serve to demonstrate how selection effects caused by a combination of instrumental sensitivity and the cosmological nature of an astrophysical population can act to produce an artificially strong correlation between observed properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The performance of ensemble-based free energy protocols in computing binding affinities to ROS1 kinase

Optimization of binding affinities for compounds to their target protein is a primary objective in drug discovery. Herein we report on a collaborative study that evaluates a set of compounds binding to ROS1 kinase. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to rank the binding free energies. The predicted binding free energies from ESMACS simulations show good correlations with experimental data for subsets of the compounds. Consistent binding free energy differences are generated for TIES and ESMACS. Although an unexplained overestimation exists, we obtain excellent statistical rankings across the set of compounds from the TIES protocol, with a Pearson correlation coefficient of 0.90 between calculated and experimental activities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic scale etching of diamond: insights from molecular dynamics simulations

Diamond is a promising material for multiple applications in quantum information processing and sensing as well as applications in microelectronics. However, diamond devices can be limited by surface defects that compromise charge stability and spin coherence, among others. Improved strategies in plasma etching of diamond could play an important role in minimizing or eliminating these defects. In this work, we explore plasma-assisted atomic scale etching of diamond using argon ions (Ar + ), hydrogen ions (H + ) and hydrogen atoms (H). We employ classical molecular dynamics (MD) simulations and test several interatomic potentials based on the Reactive Empirical Bond Order (REBO) form with comparisons to a variety of published experimental results. We performed MD simulations of low-energy hydrogen ($\leqslant$50 eV) and argon ( $\leqslant$200 eV) ion bombardment of diamond surfaces. Ar + bombardment can be used to locally smooth initially rough diamond surfaces via the formation of an amorphous C layer, the thickness of which increases with argon ion energy. Subsequent exposure with hydrogen ions (or fast neutrals) will selectively etch this amorphous C layer, leaving the underlying diamond layer mostly intact if the H energy is maintained below about 10 eV. The simulations suggest that combining Ar + smoothing with selective, near threshold energy H removal of amorphous C can be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices.

74 ATOMIC AND MOLECULAR PHYSICS↗

3D Simulations Capture the Persistent Low-Mode Asymmetries Evident in Laser-Direct-Drive Implosions on OMEGA

Spherical implosions in Inertial Confinement Fusion (ICF) are inherently sensitive to perturbations that may arise from experimental constraints and errors. Control and mitigation of low-mode (long wavelengths) perturbations is a key milestone to improving implosion performances. Here, we present the first 3-D radiation-hydrodynamic simulations of directly driven ICF implosions with an inline package for polarized Crossed-Beam Energy Transfer (CBET). Simulations match bang times, yields (separately accounting for laser-induced high modes and fuel age), hot spot flow velocities and direction, for which polarized CBET contributes to the systematic flow orientation evident in the OMEGA implosion database. Current levels of beam mispointing, imbalance, target offset and asymmetry from polarized CBET degrade yields by more than 40%. The effectiveness of two mitigation strategies for low-modes is explored.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design of a cavity heat pipe receiver experiment

A cavity heat pipe experiment has been designed to test the critical issues involved with incorporating thermal energy storage canisters into a heat pipe. The experiment is a replication of the operation of a heat receiver for a Brayton solar dynamic power cycle. The heat receiver is composed of a cylindrical receptor surface and an annular heat pipe with thermal energy storage canisters and gaseous working fluid heat exchanger tubes surrounding it. Hardware for the cavity heat pipe experiment will consist of a sector of the heat pipe, complete with gas tube and thermal energy storage canisters. Thermal cycling tests will be performed on the heat pipe sector to simulate the normal energy charge/discharge cycle of the receiver in a spacecraft application.

Schneider, Michael G.↗

First Steps Toward Incorporating Image Based Diagnostics Into Particle Accelerator Control Systems Using Convolutional Neural Networks

At present, a variety of image-based diagnostics are used in particle accelerator systems. Often times, these are viewed by a human operator who then makes appropriate adjustments to the machine. Given recent advances in using convolutional neural networks (CNNs) for image processing, it should be possible to use image diagnostics directly in control routines (NN-based or otherwise). This is especially appealing for non-intercepting diagnostics that could run continuously during beam operation. Here, we show results of a first step toward implementing such a controller: our trained CNN can predict multiple simulated downstream beam parameters at the Fermilab Accelerator Science and Technology (FAST) facility's low energy beamline using simulated virtual cathode laser images, gun phases, and solenoid strengths.

43 PARTICLE ACCELERATORS↗

NEAMS Update: Quarterly Report for October - December 2011

The Advanced Modeling and Simulation Office within the DOE Office of Nuclear Energy (NE) has been charged with revolutionizing the design tools used to build nuclear power plants during the next 10 years. To accomplish this, the DOE has brought together the national laboratories, U.S. universities, and the nuclear energy industry to establish the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program. The mission of NEAMS is to modernize computer modeling of nuclear energy systems and improve the fidelity and validity of modeling results using contemporary software environments and high-performance computers. NEAMS will create a set of engineering-level codes aimed at designing and analyzing the performance and safety of nuclear power plants and reactor fuels. The truly predictive nature of these codes will be achieved by modeling the governing phenomena at the spatial and temporal scales that dominate the behavior. These codes will be executed within a simulation environment that orchestrates code integration with respect to spatial meshing, computational resources, and execution to give the user a common 'look and feel' for setting up problems and displaying results. NEAMS is building upon a suite of existing simulation tools, including those developed by the federal Scientific Discovery through Advanced Computing and Advanced Simulation and Computing programs. NEAMS also draws upon existing simulation tools for materials and nuclear systems, although many of these are limited in terms of scale, applicability, and portability (their ability to be integrated into contemporary software and hardware architectures). NEAMS investments have directly and indirectly supported additional NE research and development programs, including those devoted to waste repositories, safeguarded separations systems, and long-term storage of used nuclear fuel. NEAMS is organized into two broad efforts, each comprising four elements. The quarterly highlights October-December 2011 are: (1) Version 1.0 of AMP, the fuel assembly performance code, was tested on the JAGUAR supercomputer and released on November 1, 2011, a detailed discussion of this new simulation tool is given; (2) A coolant sub-channel model and a preliminary UO{sub 2} smeared-cracking model were implemented in BISON, the single-pin fuel code, more information on how these models were developed and benchmarked is given; (3) The Object Kinetic Monte Carlo model was implemented to account for nucleation events in meso-scale simulations and a discussion of the significance of this advance is given; (4) The SHARP neutronics module, PROTEUS, was expanded to be applicable to all types of reactors, and a discussion of the importance of PROTEUS is given; (5) A plan has been finalized for integrating the high-fidelity, three-dimensional reactor code SHARP with both the systems-level code RELAP7 and the fuel assembly code AMP. This is a new initiative; (6) Work began to evaluate the applicability of AMP to the problem of dry storage of used fuel and to define a relevant problem to test the applicability; (7) A code to obtain phonon spectra from the force-constant matrix for a crystalline lattice has been completed. This important bridge between subcontinuum and continuum phenomena is discussed; (8) Benchmarking was begun on the meso-scale, finite-element fuels code MARMOT to validate its new variable splitting algorithm; (9) A very computationally demanding simulation of diffusion-driven nucleation of new microstructural features has been completed. An explanation of the difficulty of this simulation is given; (10) Experiments were conducted with deformed steel to validate a crystal plasticity finite-element code for bodycentered cubic iron; (11) The Capability Transfer Roadmap was completed and published as an internal laboratory technical report; (12) The AMP fuel assembly code input generator was integrated into the NEAMS Integrated Computational Environment (NiCE). More details on the planned NEAMS computing environment is given; and (13) The NEAMS program website (neams.energy.gov) is nearly ready to launch.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design of specially adapted reactive coordinates to economically compute potential and kinetic energy operators including geometry relaxation

Quantum dynamics simulations require prior knowledge of the potential energy surface as well as the kinetic energy operator. Typically, they are evaluated in a low-dimensional subspace of the full configuration space of the molecule as its dimensionality increases proportional to the number of atoms. This entails the challenge to find the most suitable subspace. We present an approach to design specially adapted reactive coordinates spanning this subspace. In addition to the essential geometric changes, these coordinates take into account the relaxation of the non-reactive coordinates without the necessity of performing geometry optimizations at each grid point. The method is demonstrated for an ultrafast photoinduced bond cleavage in a commonly used organic precursor for the generation of electrophiles. The potential energy surfaces for the reaction as well as the Wilson G-matrix as part of the kinetic energy operator are shown for a complex chemical reaction, both including the relaxation of the non-reactive coordinates on equal footing. A microscopic interpretation of the shape of the G-matrix elements allows to analyze the impact of the non-reactive coordinates on the kinetic energy operator. Additionally, we compare quantum dynamics simulations with and without the relaxation of the non-reactive coordinates included in the kinetic energy operator to demonstrate its influence.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Wall effects observed in tissue-equivalent proportional counters from 1.05 GeV/nucleon iron-56 particles

Tissue-equivalent proportional counters (TEPCs) have been used to measure energy deposition in simulated volumes of tissue ranging in diameter from 0.1 to 10 microm. There has been some concern that the wall used to define the volume of interest could influence energy deposition within the sensitive volume because it has a density significantly greater than that of the cavity gas. These effects become important for high-velocity heavy ions. Measurements of energy deposition were made for 1 GeV/nucleon iron particles in a TEPC simulating a 1-microm-diameter sphere of tissue. The TEPC was nested within a particle spectrometer that provided identification and flight path of individual particles. Energy deposition was studied as a function of pathlength through the TEPC. Approximately 30% of the energy transfer along trajectories through the center of the detector escapes the sensitive volume. The response of the TEPC, for trajectories through the detector, is always larger than calculations for energy loss in a homogeneous medium. This enhancement is greatest for trajectories near the cavity/wall interface. An integration of the response indicates that charged-particle equilibrium is essentially achieved for a wall thickness of 2.54 mm. However, estimates of the linear energy transfer for the incident particles are influenced by these wall effects.

Non-NASA Center↗

Developing Equivalent Surface Heat Transfer Input Parameters for Updated Standard 140/BESTEST Thermal Fabric Test Cases

ANSI/ASHRAE Standard 140, Standard Method of Test for the Evaluation of Building Energy Analysis Computer Programs, is widely cited for evaluating software for use with performance-path energy efficiency analysis in conjunction with well-known energy-efficiency standards, including: ASHRAE Standard 90.1, the International Energy Conservation Code, and other international standards. This paper summarizes recent work to develop improved equivalent surface heat transfer input parameters for the building thermal fabric comparative test cases that comprise a portion of Standard 140. The software-to-software comparative test method allows comparison of simulation results from a wide variety of building energy simulation programs to each other. In this type of test method, it is of great importance to specify inputs that are as well-equivalenced as possible to accommodate a wide variety of input schemes for the tested programs, and thereby reduce the potential for disagreements among them in the tests. Results include new alternative constant surface coefficients that are appropriate to the test cases and demonstrate their ability to approximate the more detailed surface heat transfer models from which they are derived.

BESTEST↗