Search NASA⌕ Search

SEARCH · Search NASA

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,765 records · Page 98

Computing the reconnection rate in turbulent kinetic layers by using electron mixing to identify topology

Three-dimensional kinetic simulations of magnetic reconnection for parameter regimes relevant to the magnetopause current layer feature the development of turbulence, driven by the magnetic and velocity shear, and dominated by coherent structures including flux ropes, current sheets, and flow vortices. Here, we propose a new approach for computing the global reconnection rate in the presence of this complexity. The mixing of electrons originating from separate sides of the magnetopause layer is used as a proxy to rapidly identify the magnetic topology and track the evolution of magnetic flux. The details of this method are illustrated for an asymmetric current layer relevant to the subsolar magnetopause and for a flow shear dominated layer relevant to the lower latitude magnetopause. While the three-dimensional reconnection rates show a number of interesting differences relative to the corresponding two-dimensional simulations, the time scale for the energy conversion remains very similar. These results suggest that the mixing of field lines between topologies is more easily influenced by kinetic turbulence than the physics responsible for the energy conversion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hardware-in-the-Loop Simulation for the Proposed Slider-Crank Wave Energy Conversion Device

In this work, a novel model of a wave energy converter that utilizes a slider crank power take-off system (PTOS) is simulated in a hardware-in-the-loop framework. The simulated parts in this framework include a buoy, slider crank linkage, and a gearbox. The slider crank linkage acting as the linear to rotary converter is represented by a DC motor coupled with a physical generator. With a regular sinusoidal wave excitation force in simulation, the hydrodynamic forces are converted into a torque, emulated by a physical motor that is directly coupled with another identical physical DC machine as the generator. By utilizing the dSPACE Control Desk software and the associated dSPACE hardware, the control algorithm in the device uses a power electronics board to drive the machines. The results show that the control system applied to the hardware-in-the-loop framework can successfully extract energy for a scaled down version of the proposed WEC device.

16 TIDAL AND WAVE POWER↗

Transforming Windows from Energy Liabilities to Zero-Energy Assets: Next-Generation Solutions for Buildings

Windows have traditionally contributed to a building's HVAC load, but they can also become a source of net energy gain or even operate as zero-energy components. For heating applications, highly insulating windows can harness more solar heat than the energy lost through them, transforming windows from energy liabilities to assets. Dynamic glazings provide further benefits by regulating solar heat gain, reducing cooling loads in summer and heating demands in winter. This simulation study focuses on developing the next generation of zero-energy windows (ZEW) for residential new construction. Through annual energy simulations across climate zones 1-8, ZEW performance benchmarks were established based on current code-level buildings, and we've identified the regions where meeting ZEW standards are most achievable. This work evaluates both static and dynamic window technologies, assessing their effects on annual energy use and cost. Key findings demonstrate that ZEW performance is achievable across diverse climate zones, with specific regional requirements. Most climate zones from 3-8 can achieve ZEW with specific configurations, while some warm climates (1-2) appear challenging for ZEW implementation. Climate zones 4-6 consistently allow for zero energy window implementation, offering multiple pathways through either static or dynamic window technologies. Colder climate zones (7-8) ZEW products allow for higher SHGC values while requiring low U-values.

Yu, Lili↗

Simulation budgeting for hybrid effective field theories

In this work, we forecast the number of, and requirements on, N-body simulations needed to train hybrid effective field theory (HEFT) emulators for a range of use cases, using a hybrid of HMcode and perturbation theory as a surrogate model. Our accuracy goals, determined with careful consideration of statistical and systematic uncertainties, are 1% accurate in the high-likelihood range of cosmological parameters, and 2% accurate over a broader parameter space volume for k < 1 h Mpc -1 and z < 3. Focusing in part on the 8-parameter w 0 w a CDM+ m ν cosmological model, we find that < 225 simulations are required to meet our error goals over our wide parameter space, including models with rapidly evolving dark energy, given our simulation and emulator recommendations. For a more restricted parameter space volume, as few as 80 simulations are sufficient. We additionally present simulation forecasts for example use cases, and make the code used in our analyses publicly available. These results offer practical guidance for efficient emulator design and simulation budgeting in future cosmological analyses.

cosmological parameters from LSS↗

Deforestation reshapes land-surface energy-flux partitioning

Land-use and land-cover change significantly modify local land-surface characteristics and water/energy exchanges, which can lead to atmospheric circulation and regional climate changes. In particular, deforestation accounts for a large portion of global land-use changes, which transforms forests into other land cover types, such as croplands and grazing lands. Many previous efforts have focused on observing and modeling land–atmosphere–water/energy fluxes to investigate land–atmosphere coupling induced by deforestation. However, interpreting land–atmosphere–water/energy-flux responses to deforestation is often complicated by the concurrent impacts from shifts in land-surface properties versus background atmospheric forcings. In this study, we used 29 paired FLUXNET sites, to improve understanding of how deforested land surfaces drive changes in surface-energy-flux partitioning. Each paired sites included an intact forested and non-forested site that had similar background climate. We employed transfer entropy, a method based on information theory, to diagnose directional controls between coupling variables, and identify nonlinear cause–effect relationships. Transfer entropy is a powerful tool to detective causal relationships in nonlinear and asynchronous systems. The paired eddy covariance flux measurements showed consistent and strong information flows from vegetation activity (gross primary productivity (GPP)) and physical climate (e.g. shortwave radiation, air temperature) to evaporative fraction (EF) over both non-forested and forested land surfaces. More importantly, the information transfers from radiation, precipitation, and GPP to EF were significantly reduced at non-forested sites, compared to forested sites. We then applied these observationally constrained metrics as benchmarks to evaluate the Energy Exascale Earth System Model (E3SM) land model (ELM). ELM predicted vegetation controls on EF relatively well, but underpredicted climate factors on EF, indicating model deficiencies in describing the relationships between atmospheric state and surface fluxes. Moreover, changes in controls on surface energy flux partitioning due to deforestation were not detected in the model. We highlight the need for benchmarking model simulated surface-energy fluxes and the corresponding causal relationships against those of observations, to improve our understanding of model predictability on how deforestation reshapes land surface energy fluxes.

54 ENVIRONMENTAL SCIENCES↗

Estimating Cosmic Ray Spectral Parameters From Simulated Detector Responses With Detector Design Implications

A simple power law model consisting of a single spectral index alpha (sub 1), is believed to be an adequate description of the galactic cosmic ray (GCR) proton flux at energies below 10(exp 13) eV, with a transition at knee energy E(sub k) to a steeper spectral index alpha(sub 2) greater than alpha(sub 1) above E(sub k). The maximum likelihood procedure is developed for estimating these three spectral parameters of the broken power law energy spectrum from simulated detector responses. These estimates and their surrounding statistical uncertainty are being used to derive the requirements in energy resolution, calorimeter size, and energy response of a proposed sampling calorimeter for the Advanced Cosmic ray Composition Experiment for the Space Station (ACCESS). This study thereby permits instrument developers to make important trade studies in design parameters as a function of the science objectives, which is particularly important for space-based detectors where physical parameters, such as dimension and weight, impose rigorous practical limits to the design envelope.

Howell, L. W.↗

Estimating Cosmic-Ray Spectral Parameters from Simulated Detector Responses with Detector Design Implications

A simple power law model consisting of a single spectral index (alpha-1) is believed to be an adequate description of the galactic cosmic-ray (GCR) proton flux at energies below 10(exp 13) eV, with a transition at knee energy (E(sub k)) to a steeper spectral index alpha-2 > alpha-1 above E(sub k). The maximum likelihood procedure is developed for estimating these three spectral parameters of the broken power law energy spectrum from simulated detector responses. These estimates and their surrounding statistical uncertainty are being used to derive the requirements in energy resolution, calorimeter size, and energy response of a proposed sampling calorimeter for the Advanced Cosmic-ray Composition Experiment for the Space Station (ACCESS). This study thereby permits instrument developers to make important trade studies in design parameters as a function of the science objectives, which is particularly important for space-based detectors where physical parameters, such as dimension and weight, impose rigorous practical limits to the design envelope.

Howell, L. W.↗

Distinct interactions of Na{sup +} and Ca{sup 2+} ions with the selectivity filter of the bacterial sodium channel Na{sub V}Ab

Highlights: ► Ca{sup 2+} translocates slowly in the filter, due to lack of “loose” knock-on mechanism. ► Identification of a high affinity binding site in Na{sub V}Ab selectivity filter. ► Changes of EEEE locus triggered by electrostatic interactions with Ca{sup 2+} ions. -- Abstract: Rapid and selective ion transport is essential for the generation and regulation of electrical signaling pathways in living organisms. In this study, we use molecular dynamics simulations and free energy calculations to investigate how the bacterial sodium channel Na{sub V}Ab (Arcobacter butzleri) differentiates between Na{sup +} and Ca{sup 2+} ions. Multiple nanosecond molecular dynamics simulations revealed distinct binding patterns for these two cations in the selectivity filter and suggested a high affinity calcium binding site formed by backbone atoms of residues Leu-176 and Thr-175 (S{sub CEN}) in the sodium channel selectivity filter.

60 APPLIED LIFE SCIENCES↗

Simulation and evaluation of latent heat thermal energy storage

The relative value of thermal energy storage (TES) for heat pump storage (heating and cooling) as a function of storage temperature, mode of storage (hotside or coldside), geographic locations, and utility time of use rate structures were derived. Computer models used to simulate the performance of a number of TES/heat pump configurations are described. The models are based on existing performance data of heat pump components, available building thermal load computational procedures, and generalized TES subsystem design. Life cycle costs computed for each site, configuration, and rate structure are discussed.

Sigmon, T. W.↗

Warthog: Coupling Nek5000 Thermal Hydraulics to BISON Fuel Performance through the Giraffe Interface

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program under the US Department of Energy (DOE) consists of three product lines. The Fuels Product Line (FPL), hosted at Idaho National Laboratory (INL), has produced a multiphysics framework with the name Multiphysics Object Oriented Simulation Environment (MOOSE) (Gaston et al. 2015). The FPL has also produced several high-fidelity physics codes including the fuel performance code Bison (Hales et al. 2016). The Reactors Product Line (RPL), hosted at Argonne National Laboratory (Argonne), has developed the high-order spectral-element Computational Fluid Dynamics (CFD) code Nek5000 (Fischer et al. 2015). The Integration Product Line (IPL), hosted at Oak Ridge National Laboratory (ORNL), is tasked with ensuring tools from the other product lines can be coupled together so that users may perform high-fidelity multiphysics simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth

The response of the cross-beam energy transfer instability (CBET) to laser bandwidth is investigated through a combination of theory and simulation. Existing linear theory is generalized to treat broadband lasers, demonstrating that CBET is most effectively suppressed when the bandwidth exceeds the ion-acoustic wave (IAW) frequency. It is shown that for such bandwidths, reverse (seed to pump) transfer becomes possible, which reduces the net energy transfer rapidly as bandwidth is increased. The CBET gain exponent in this regime scales with bandwidth (Δω) as Δω−3 for Gaussian or Lorentzian laser spectra with different scalings possible for other spectra. Comparison of our theory with linearized fluid and particle-in-cell simulations, performed with the laser-plasma simulation environment (LPSE) and vector particle in cell (VPIC) codes, respectively, finds that the model is accurate in the absence of nonlinear processes. However, linear analysis also finds that the IAW energy density scales as Δω−1, implying that nonlinear effects may be more difficult to control than the CBET scaling would suggest. Indeed, nonlinear effects are found to be present in VPIC simulations with high-intensity lasers, despite minimal apparent CBET. Nonlinear processes in the VPIC cases include particle trapping, the two-ion wave decay, and ion wave self-focusing. In some high intensity VPIC cases, these effects lead to net energy transfer from seed to pump and increases to backscatter stimulated Brillouin scattering reflectivities. Finally, for a given bandwidth, we show that improved control of nonlinear processes can be achieved via smoothing by spectral dispersion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stopping and transport of fast electrons in superdense matter

Studied is the stopping and transport of relativistic fast electrons in the vicinity of compressed dense plasma core relevant to fast ignition. Electromagnetic cascade Monte-Carlo is coupled to 2D-PIC simulation. The 2D PIC simulates input electron energy spectrum and angular dependence. The electron energy distributions after passing through the plasma core are calculated at different viewing angles, which well agree with the experiment below several MeV energy range. The implications of calculated results as to collisional damping on several MeV electrons are discussed with the theory based on the stopping power model. The spatial distribution of plasma temperature is also estimated via deposited energy by fast electrons, showing the strong heating at the core surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neutrino vertex reconstruction with in-ice radio detectors using surface reflections and implications for the neutrino energy resolution

Ultra high energy neutrinos (E ν > 10 16.5 eV) are efficiently measured via radio signals following a neutrino interaction in ice. An antenna placed1Script O(15 m) below the ice surface will measure two signals for the vast majority of events (90% at E ν =10 18 eV): a direct pulse and a second delayed pulse from a reflection off the ice surface. This allows for a unique identification of neutrinos against backgrounds arriving from above. Furthermore, the time delay between the direct and reflected signal (D'n'R) correlates with the distance to the neutrino interaction vertex, a crucial quantity to determine the neutrino energy. In a simulation study, we derive the relation between time delay and distance and study the corresponding experimental uncertainties in estimating neutrino energies. We find that the resulting contribution to the energy resolution is well below the natural limit set by the unknown inelasticity in the initial neutrino interaction. We present an in-situ measurement that proves the experimental feasibility of this technique. Continuous monitoring of the local snow accumulation in the vicinity of the transmit and receive antennas using this technique provide a precision of Script O(1 mm) in surface elevation, which is much better than that needed to apply the D'n'R technique to neutrinos.

79 ASTRONOMY AND ASTROPHYSICS↗

Energy Partition in Collisionless Counterstreaming Plasmas

Fast, counterstreaming plasma outflows drive magnetic field amplification, plasma heating, and particle acceleration in numerous astrophysical environments, from supernova remnant shocks to active galactic nuclei jets. Understanding how, in the absence of Coulomb collisions, energy is redistributed between the different plasma species remains a fundamental open question. We use three-dimensional, fully kinetic simulations to investigate energy partition in weakly magnetized counterpropagating plasmas. Our results reveal a complex interplay between different processes, where at early times the Weibel instability drives a first stage of magnetic field amplification and at late times the kinking of current filaments drives a second amplification stage via a dynamo-type mechanism. Electrons are heated primarily during the latter phase through magnetic pumping. By the time the flows thermalize, we observe that the final temperature ratio T e /T i and energy partition depend on the ion-to-electron mass ratio. For electron–proton flows, the electron thermal energy only reaches up to a few percent of the initial ion kinetic energy.

Marret, Alexis [SLAC National Accelerator Laborato↗

Dielectric Properties of Poly(ethylene oxide) from Molecular Dynamics Simulations

The order, conformations and dynamics of poly(oxyethylene) (POE) melts have been investigated through molecular dynamics simulations. The potential energy functions were determined from detailed ab initio electronic structure calculations of the conformational energies of the model molecules 1,2-dimethoxyethane (DME) and diethylether. The x-ray structure factor for POE from simulation will be compared to experiment. In terms of conformation, simulations reveal that chains are extended in the melt relative to isolated chains due to the presence of strong intermolecular O...H interactions, which occur at the expense of intramolecular O...H interactions. Conformational dynamics about the C-C bond were found to be significantly faster than in polymethylene, while conformational dynamics about the C-O bond even faster than the C-C dynamics. The faster local dynamics in POE relative to polymethylene is consistent with C-13 NMR spin-lattice relaxation experiments. Conformational transitions showed significant second-neighbor correlation, as was found for polymethylene. This correlation of transitions with C-C neighbors was found to be reduced relative to C-O neighbors. Dielectric relaxation from simulation will also be compared with experiment.

Smith, Grant D.↗

Comparison of Martian Surface Radiation Predictions to the Measurements of Mars Science Laboratory Radiation Assessment Detector (MSL/RAD)

For the analysis of radiation risks to astronauts and planning exploratory space missions, detailed knowledge of particle spectra is an important factor. Detailed measurements of the energetic particle radiation environment on the surface of Mars have been made by the Mars Science Laboratory Radiation Assessment Detector (MSL-RAD) on the Curiosity rover since August 2012, and particle fluxes for a wide range of ion species (up to several hundred MeV/u) and high energy neutrons (8 - 1000 MeV) have been available for the first 200 sols. Although the data obtained on the surface of Mars for 200 sols are limited in the narrow energy spectra, the simulation results using the Badhwar-O'Neill galactic cosmic ray (GCR) environment model and the high-charge and energy transport (HZETRN) code are compared to the data. For the nuclear interactions of primary GCR through Mars atmosphere and Curiosity rover, the quantum multiple scattering theory of nuclear fragmentation (QMSFRG) is used, which includes direct knockout, evaporation and nuclear coalescence. Daily atmospheric pressure measurements at Gale Crater by the MSL Rover Environmental Monitoring Station are implemented into transport calculations for describing the daily column depth of atmosphere. Particles impinging on top of the Martian atmosphere reach the RAD after traversing varying depths of atmosphere that depend on the slant angles, and the model accounts for shielding of the RAD by the rest of the instrument. Calculations of stopping particle spectra are in good agreement with the RAD measurements for the first 200 sols by accounting changing heliospheric conditions and atmospheric pressure. Detailed comparisons between model predictions and spectral data of various particle types provide the validation of radiation transport models, and thus increase the accuracy of the predictions of future radiation environments on Mars. These contributions lend support to the understanding of radiation health risks to astronauts for the planning of various mission scenarios.

Kim, Myung-Hee Y.↗

Quantum simulations of nuclei and nuclear pasta with the multiresolution adaptive numerical environment for scientific simulations

Neutron star and supernova matter at densities just below the nuclear matter saturation density is expected to form a lattice of exotic shapes. These so-called nuclear pasta phases are caused by Coulomb frustration. Their elastic and transport properties are believed to play an important role for thermal and magnetic field evolution, rotation, and oscillation of neutron stars. Furthermore, they can impact neutrino opacities in core-collapse supernovae. In this work, we present proof-of-principle three-dimensional (3D) Skyrme Hartree-Fock (SHF) simulations of nuclear pasta with the Multi-resolution ADaptive Numerical Environment for Scientific Simulations (MADNESS). We perform benchmark studies of 16 O, 208 Pb, and 238 U nuclear ground states and calculate binding energies via 3D SHF simulations. Results are compared with experimentally measured binding energies as well as with theoretically predicted values from an established SHF code. The nuclear pasta simulation is initialized in the so-called waffle geometry as obtained by the Indiana University Molecular Dynamics (IUMD) code. The size of the unit cell is 24 fm with an average density of about ρ = 0.05fm -3 , proton fraction of Yp = 0.3, and temperature of T = 0 MeV. Overall, our calculations reproduce the binding energies and shapes of light and heavy nuclei with different geometries. For the pasta simulation, we find that the final geometry is very similar to the initial waffle state. We compare calculations with and without spin-orbit forces. We find that while subtle differences are present, the pasta phase remains in the waffle geometry. Within the MADNESS framework, we can successfully perform calculations of inhomogeneous nuclear matter. By using pasta configurations from IUMD it is possible to explore different geometries and test the impact of self-consistent calculations on the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Structural features of ISG borosilicate nuclear waste glasses revealed from high-energy X-ray diffraction and molecular dynamics simulations

Vitrification is a widely accepted method to immobilize nuclear waste. Detailed structural information is critical to understand the physical and chemical behaviors, including the long-term chemical durability of these glasses that are to be stored in geological storage sites. High-energy X-ray diffraction studies are used to obtain accurate structural information of the International Simple Glass (ISG), a six-component borosilicate model nuclear waste glass with a composition (mol%) of 60.2SiO 2 -16.0B 2 O 3 -12.6Na 2 O-3.8Al 2 O 3 -5.7CaO-1.7ZrO 2 . Classical molecular dynamics (MD) simulations with recently developed effective potentials are utilized to generate structural models of the ISG glasses to provide interpretation of the high-energy X-ray structural data and additional short and medium range structural details. The atomic structure model generated from MD simulations shows an excellent agreement with high-energy X-ray diffraction, measured for the first time for ISG, with an R x value of 5.5%. Systematic investigations of the medium-range structural features of ISG are performed, built upon recent study of the short-range structure. In particular, the glass former cation oxygen polyhedral analysis show that [SiO 4 ]-[SiO 4 ] and [BO 4 ]-[SiO 4 ] connections are the most abundant while [ZrO 6 ] octohedra having higher probability to connect to another [SiO 4 ], [BO 4 ] and [BO 3 ]. Moreover, we have studied the ISG surface structure using MD simulations. Noticeable differences in terms of both chemical compositions and structural features between the bulk and surface are observed. There exist increased small-sized rings, under coordinated Si species, high 3-coordinated boron content (86.5%) on the surface. Very importantly, sodium ions are found to be enriched on the glass surface, consistent with recent experimental results. Furthermore, both bulk and surface structural information are discussed in terms of the dissolution and corrosion mechanisms of ISG and related nuclear waste glasses.

36 MATERIALS SCIENCE↗