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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

CompF2: Theoretical Calculations and Simulation Topical Group Report

This report summarizes the work of the Computational Frontier topical group on theoretical calculations and simulation for Snowmass 2021. We discuss the challenges, potential solutions, and needs facing six diverse but related topical areas that span the subject of theoretical calculations and simulation in high energy physics (HEP): cosmic calculations, particle accelerator modeling, detector simulation, event generators, perturbative calculations, and lattice QCD (quantum chromodynamics). The challenges arise from the next generations of HEP experiments, which will include more complex instruments, provide larger data volumes, and perform more precise measurements. Calculations and simulations will need to keep up with these increased requirements. The other aspect of the challenge is the evolution of computing landscape away from general-purpose computing on CPUs and toward special-purpose accelerators and coprocessors such as GPUs and FPGAs. These newer devices can provide substantial improvements for certain categories of algorithms, at the expense of more specialized programming and memory and data access patterns.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Energy Conversion and Electron Acceleration and Transport in 3D Simulations of Solar Flares

Recent observations and simulations indicate that solar flares undergo extremely complex 3D evolution, making 3D particle transport models essential for understanding electron acceleration and interpreting flare emissions. In this study, we investigate this problem by solving Parker’s transport equation with 3D MHD simulations of solar flares. By examining energy conversion in the 3D system, we evaluate the roles of different acceleration mechanisms, including reconnection current sheet (CS), termination shock (TS), and supra-arcade downflows (SADs). We find that large-amplitude turbulent fluctuations are generated and sustained in the 3D system. The model results demonstrate that a significant number of electrons are accelerated to hundreds of keV and even a few MeV, forming power-law energy spectra. These energetic particles are widely distributed, with concentrations at the TS and in the flare looptop region, consistent with results derived from recent hard X-ray (HXR) and microwave (MW) observations. By selectively turning particle acceleration on or off in specific regions, we find that the CS and SADs effectively accelerate electrons to several hundred keV, while the TS enables further acceleration to MeV. However, no single mechanism can independently account for the significant number of energetic electrons observed. Instead, the mechanisms work synergistically to produce a large population of accelerated electrons. Our model provides spatially and temporally resolved electron distributions in the whole flare region and at the flare footpoints, enabling synthetic HXR and MW emission modeling for comparison with observations. These results offer important insights into electron acceleration and transport in 3D solar flare regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

How Accurate Are Simulations and Experiments for the Lattice Energies of Molecular Crystals?

Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a delicate interplay of intermolecular interactions such as hydrogen bonding and Van der Waals dispersion forces. Here, by exploiting recent algorithmic developments, we report the first set of diffusion Monte Carlo lattice energies for all 23 molecular crystals in the popular and widely used X23 dataset. Comparisons with previous state-of-the-art lattice energy predictions (on a subset of the dataset) and a careful analysis of experimental sublimation enthalpies reveals that high-accuracy computational methods are now at least as reliable as (computationally derived) experiments for the lattice energies of molecular crystals. Overall, this work demonstrates the feasibility of high-level explicitly correlated electronic structure methods for broad benchmarking studies in complex condensed phase systems, and signposts a route towards closer agreement between experiment and simulation. Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

A Framework for Optimal Placement of Rooftop Photovoltaic: Maximizing Solar Production and Operational Cost Savings in Residential Communities

Optimizing the placement of photovoltaic (PV) panels on residential buildings has the potential to significantly increase energy efficiency benefits to both homeowners and communities. Strategic PV placement can lower electricity costs by reducing the electricity fed from the grid during on-peak hours, while maintaining PV panel efficiency in terms of the amount of solar radiation received. In this article, we present a framework that identifies the ideal location of PV panels on residential rooftops. Our framework combines energy and environmental simulation, parametric modeling, and optimization to inform PV placement as it relates to and affects the entire community (in terms of both energy use and financial cost), as well as individual buildings. Ensuring that our framework accounts for shading from nearby buildings, different utility rate structures, and different buildings’ energy demand profiles means that existing communities and future housing developments can be optimized for energy savings and PV efficiency. The framework comprises two workflows, each contributing to optimal PV placement with a unique target: (a) maximizing PV panel efficiency (i.e., solar generation) and (b) minimizing operational energy cost considering utility rate structures for operational energy. We apply our framework to a residential community in Fort Collins, Colorado, to demonstrate the optimal PV placement, considering the two workflow targets. Here, we present our results and illustrate the effect of PV location and orientation on solar energy production efficiency and operational energy cost.

14 SOLAR ENERGY↗
DOE OSTI2023

Validation of the single-event method for low-energy electron transport via stopping power calculations with $\mathrm{MCNP}$

Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Forward and backward stimulated Raman scattering in multi-speckled beams: Density dependence and effects on cross-beam energy transfer

Particle-in-cell simulations of cross-beam energy transfer (CBET) between multi-speckled beams reveal a potentially prominent energy loss from forward stimulated Raman scattering (FSRS) in the amplified seed beam [Yin et al., Phys. Plasmas 26, 082708 (2019)]. In this work, an intensity scan is used to examine this energy loss in the CBET-amplified seed beam at intensities of a few times 1014 W/cm2. Electron tracers demonstrate both the trapping in the electron plasma waves that is sustained over multiple speckles and the subsequent acceleration to velocities well beyond the phase velocity. Motivated by expanding our understanding of this energy loss mechanism, we isolate the FSRS density dependence by using a 2×1014 W/cm2 pulse to scan densities between 0.06 ncr and 0.12 ncr; a spectral analysis characterizes the energy losses from scattering with the changing densities. At saturation, both greater beam attenuation is measured at the laser frequency and larger scattering FSRS angles are observed as density increases, along with a rapid transition to primarily backscatter SRS between 0.10 ncr and 0.12 ncr, earlier than predicted by linear theory. Density gradients are shown to effectively mitigate oblique FSRS and beam energy losses in the two lower density cases, whereas 0° forward scattering is primarily suppressed in the higher density cases. Finally returning to the two-beam CBET system of interest, the interplay and mitigation of the secondary instabilities at larger domain and beam spatial scales are discussed, here also examining backward stimulated Brillouin scattering.

Stark, D. J. (ORCID:0000000203642733)↗
DOE OSTI2021

Probabilistic Resilience of DER Systems - A Simulation Assisted Optimization Approach

Energy systems resilience is becoming increasingly important as the frequency of major grid outages increases. In this work, we present a methodology to optimize a behind-the-meter distributed energy resource system to sustain a site's critical loads during a pre-defined outage period. With the fixed system design, we then propose an outage simulation approach to estimate the resilience potential of the DER system to sustain loads beyond the fixed outage period - a yearlong resilience performance analysis. We apply statistical analysis to assess the system's resilience performance over a broader parametric problem space on an hourly, monthly, and yearly basis. We demonstrate the impact of the pre-defined outage period on the resilience performance through a case study. Results show that the probability of surviving a random outage of a given duration changes from 20% to 95% when the outage is modeled for a weekday instead of a weekend for the given load-profile.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Formation and dissociation of shear-induced high-energy dislocations: insight from molecular dynamics simulations

Solid-phase processing (SPP) allows one to create complex microstructures, not achievable via thermal processing alone. The resulting structures exhibit a rich palette of defects, both thermal and non-thermal, including defect substructures, such as dislocation networks. It is essential to understand the mechanisms of deformation and defect structure formation to guide SPP towards achieving desired microstructures and material properties. In this study, large-scale molecular dynamics simulations are used to investigate the effects of inhomogeneous strain distribution, that mimics deformation conditions of tribological tests, on the evolution of defects under severe shear deformation in polycrystalline Al. Analysis of defect nucleation and reaction pathways reveals that strong geometric constraints suppress the nucleation and slide of low energy dislocation 1/2$\langle110\rangle${111} but promote the nucleation and slide of high energy dislocations, such as [$1\bar1$0] (001) and 1/2 [$1\bar1$$\bar2$]($1\bar1$1). Further, a rough contact surface, characteristic to tribological tests, imposes an inhomogeneous stress field leading to inhomogeneous defect substructures due to location-dependent activation of slip systems. The results suggest that high-energy dislocations can dominate the evolution of grain structures in highly constrained environments, which should be considered in modeling plastic deformation and grain refinement during SPP.

42 ENGINEERING↗
DOE OSTI2023

Coupled Decay Heat and Thermal Hydraulic Capability for Loss-of-Coolant Accident Simulations

As the nuclear energy industry considers ways to achieve improved economics in the current fleet of light-water reactors (LWRs), one possible approach is to operate each cycle for longer durations. This causes a greater portion of the fuel to be burned and reduces the frequency of outages, which ultimately reduces the cost to operate the reactor. However, this also leads to higher burnup fuels than has traditionally been allowed in these reactors. Thus, there are concerns about integrity of high-burnup (HBu) fuel, especially during accident conditions such as loss-of-coolant accidents (LOCAs), as shown by Capps et al.. To investigate these concerns, advanced modeling and simulation capabilities are being leveraged to determine the susceptibility of HBu fuel to fuel fragmentation, relocation, and dispersion (FFRD). Improvements have previously been made to fuel performance capabilities to more accurately model these phenomena; multiphysics simulations have also been conducted to determine the power and burnup histories of the HBu fuel, which are needed as inputs for the fuel performance calculations. Most recently, new statistical approaches have been developed to identify a subset of fuel rods that have greater FFRD susceptibility, reducing the total number of fuel performance simulations required. Prior LOCA simulations have relied on the TRACE systems code, which can model the core and primary loop during accident conditions. TRACE includes many models for various aspects of the primary loop, but two sets of models are important for this report. First, TRACE uses a lumped-fuel approach for modeling the core. This approximates the ~50,000 fuel rods in the core with a much smaller number of rods. The rods can be lumped in various ways as determined by the user. For example, one lumped rod may be used to represent all rods in an assembly, sometimes with an additional rod representing the hottest fuel rod. However, due to runtime constraints and complexity of modeling, a more common approach is to group several assemblies or larger regions of the core into single lumped rods. These lumping schemes apply not only to fuel rods but to flow channels as well. Second, TRACE has several different models for treating decay heat, ranging from pregenerated decay heat curves based on an ANSI/ANS-5.1 standard (hereinafter abbreviated simply as ANSI) to explicit time-dependent heat inputs from the user. None of these models account for differences in isotopics between different rods, which is an approximation the work in this report seeks to eliminate. This report focuses on the implementation of coupled decay heat capabilities in the Virtual Environment for Reactor Applications (VERA) code suite to address a gap identified in previous LOCA simulations. This constitutes an improvement for both the lumped-fuel and decay heat models in TRACE. VERA has been developed to perform high-fidelity, whole-core multiphysics simulations for LWRs. Previously, during the Consortium for Advanced Simulation of LWRs (CASL) program, the emphasis was on providing accurate steady-state analysis—with a secondary focus on reactivity insertion accident (RIA) analysis—to address operational challenges in the nuclear energy industry. Under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, these capabilities are being extended to a broader range of transient analyses with the goal of quantifying the risk of fuel failures such as FFRD. To properly model such conditions with VERA, decay heat calculations have been integrated with the multiphysics to enable rod-by-rod thermal hydraulic (TH) conditions to be driven by the decay heat in long-running accidents such as LOCAs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

The IDAES process modeling framework and model library—Flexibility for process simulation and optimization

Abstract Energy systems and manufacturing processes of the 21st century are becoming increasingly dynamic and interconnected, which require new capabilities to effectively model and optimize their design and operations. Such next generation computational tools must leverage state‐of‐the‐art techniques in optimization and be able to rapidly incorporate new advances. To address these requirements, we have developed the Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform, which builds on the strengths of both process simulators (model libraries) and algebraic modeling languages (advanced solvers). This paper specifically presents the IDAES Core Modeling Framework (IDAES‐CMF), along with a case study demonstrating the application of the framework to solve process optimization problems. Capabilities provided by this framework include a flexible, modifiable, open‐source platform for optimization of process flowsheets utilizing state‐of‐the‐art solvers and solution techniques, fully open and extensible libraries of dynamic unit operations models and thermophysical property models, and integrated support for superstructure‐based conceptual design and optimization under uncertainty.

42 ENGINEERING↗
DOE OSTI2020

Feasibility Analysis of Converter-Interfaced Combined Heat and Power System

As a promising new design concept, the converter-interfaced combined heat and power (CHP) system is coupled to the bulk grid through a rectifier and a grid-ready tied inverter. Compared to the traditional directly-coupled CHP system, it removes the requirement for oversizing the CHP generator, limits the short-circuit contribution of the generator and simplifies the grid integration process of CHP system. This paper evaluates the economic benefits of this concept by calculating the annualized Return-on-Investment (ROI) and comparing it to the directly-coupled system. The economic analysis includes timeseries simulations to compute energy transactions with the bulk grid as well as sizing the equipment to calculate the capital and operational costs. Obtained results indicate that in majority of user cases evaluated, the converter-interfaced CHP systems can provide better ROI than directly-coupled systems. Given the additional technical benefits provided by inverter-based distributed energy resources (DERs), the proposed concept is proved to be technically viable and economically feasible.

03 NATURAL GAS↗
DOE OSTI2020

Causal Interactions between Southern Ocean Polynyas and High-Latitude Atmosphere–Ocean Variability

Abstract Weddell Sea open-ocean polynyas have been observed to occasionally release heat from the deep ocean to the atmosphere, indicating that their sporadic appearances may be an important feature of high-latitude atmosphere–ocean variability. Yet, observations of the phenomenon are sparse and many standard-resolution models represent these features poorly, if at all. We use a fully coupled, synoptic-scale preindustrial control simulation of the Energy Exascale Earth System Model (E3SMv0-HR) to effectively simulate open-ocean polynyas and investigate their role in the climate system. Our approach employs statistical tests of Granger causality to diagnose local and remote drivers of, and responses to, polynya heat loss on interannual to decadal time scales. First, we find that polynya heat loss Granger causes a persistent increase in surface air temperature over the Weddell Sea, strengthening the local cyclonic wind circulation. Along with responding to polynyas, atmospheric conditions also facilitate their development. When the Southern Ocean experiences a rapid poleward shift in the circumpolar westerlies following a prolonged negative phase of the southern annular mode (SAM), Weddell Sea salinity increases, promoting density destratification and convection in the water column. Finally, we find that the reduction of surface heat fluxes during periods of full ice cover is not fully compensated by ocean heat transport into the high latitudes. This imbalance leads to a buildup of ocean heat content that supplies polynya heat loss. These results disentangle the complex, coupled climate processes that both enable the polynya’s existence and respond to it, providing insights to improve the representation of these highly episodic sea ice features in climate models.

Kaufman, Zachary S.↗
DOE OSTI2023

Thermal energy storage in concrete: Review, testing, and simulation of thermal properties at relevant ranges of elevated temperature

This study examines the thermal performance of concrete used for thermal energy storage (TES) applications. The influence of concrete constituents (aggregates, cementitious materials, and fibers) on the thermal conductivity and specific heat are summarized based on literature and via experimentation at elevated temperatures. It is indicated that concrete with siliceous aggregate, low water/cement ratios, and steel fibers produces better overall thermal properties. The results of literature review and laboratory testing are used for numerical simulation to evaluate the impact of concrete thermal property variability on TES capabilities. The TES module consists of the heat exchanger and surrounding concrete as the storage media. Both idealized cases without heat loss and realistic cases with insulation and heat losses are considered. Furthermore, the result shows that conventional concrete mixes can be tailored to decrease the module charge/discharge durations and increase the amount of thermal energy storage.

25 ENERGY STORAGE↗
DOE OSTI2025

Scaling relations in the phase-space structure of dark matter haloes

We present new scaling relations for the isotropic phase-space distribution functions (DFs) and energy distributions of simulated dark matter haloes. These relations are inspired by those for the singular isothermal sphere with density profile $\rho (r)\propto r^{-2}$, for which the DF satisfies $f(E) \propto r_{\mathrm{max}}^{-2}(E)$ and the energy distribution satisfies $\mathrm{ d}M/\mathrm{ d}E \propto r_{\mathrm{max}}(E)$, with $r_{\mathrm{max}}(E)$ being the radius where the gravitational potential equals energy $E$. For the simulated haloes, we find $f(E)\propto r_{\mathrm{max}}^{-2.08}(E)$ and $\mathrm{ d}M/\mathrm{ d}E \propto r_{\mathrm{max}}(E)$ across broad energy ranges. In addition, the proportionality coefficients depend on the gravitational constant and the parameters of the best-fitting Navarro–Frenk–White density profile. These scaling relations are satisfied by haloes over a wide mass range and provide an efficient method to approximate their DFs and energy distributions. Understanding the origin of these relations may shed more light on halo formation.

galaxies: haloes↗
DOE OSTI2026

Multiple Coulomb scattering for muon energy reconstruction in SBND Monte Carlo simulations

The Short Baseline Near Detector (SBND) is a liquid argon time projection chamber (LArTPC) located 110 meters downstream of the Booster Neutrino Beam target at Fermilab \cite{SBN}. The LArTPC is a detector technology with unprecedented position and calorimetric resolution for precise neutrino studies. The most dominant neutrino interaction topology in SBND is charged-current quasi-elastic muon neutrino-argon interactions, where the final state consists of a muon and a proton. Here, muons often travel long distances, sometimes exiting the detector before stopping. This loss of information creates a challenge for accurate energy estimation. There have been efforts in the MicroBooNE and ICARUS collaborations to improve energy estimation using a phenomenon known as Multiple Coulomb Scattering (MCS), which describes the jagged path taken by muons through a detector due to continuous Coulomb interactions with argon atoms. Building on these efforts, my project aims to improve an MCS-based approach to estimating the energies of exiting muons in SBND.

Puig, Xilonem Angeles [Mt. Holyoke Coll.; Fermilab↗
DOE OSTI2021

Effect of resistive wall on thermal quench in JET disruptions

In this paper, experimental data, simulations and theory are presented of a JET tokamak thermal quench. The emphasis is on the timescale of the bulk plasma thermal energy loss. The simulations suggest that the thermal energy loss is caused by a resistive wall tearing mode, and experimental data is consistent with this conclusion. The time scale of the thermal quench is the inverse of the mode growth rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗