Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large system 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 379 records · Page 21

Effect of cross-beam energy transfer on target-offset asymmetry in direct-drive inertial confinement fusion implosions

The unintentional mispositioning of inertial confinement fusion (ICF) capsules from the center of laser beam convergence has long been shown in simulations to generate large ℓ = 1 asymmetry and significantly degrade implosion symmetry and fusion yields. Experimental yields on the OMEGA Laser System, however, have shown much less sensitivity to this initial target offset. This paper presents simulations of offset ICF implosions improved by including a physics model of cross-beam energy transfer (CBET), a mechanism of laser energy scattering from one beam to another. Room-temperature OMEGA implosion experiments with prescribed target offsets are simulated with and without CBET, illustrating that CBET mitigates the ℓ = 1 implosion asymmetry from target offset. Here, comparison of simulations to multiple complementary experimental observables indicates the addition of CBET physics in offset simulations is necessary to match experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of cross-beam energy transfer on target-offset asymmetry in direct-drive inertial confinement fusion implosions

The unintentional mispositioning of inertial confinement fusion (ICF) capsules from the center of laser beam convergence has long been shown in simulations to generate large ℓ = 1 asymmetry and significantly degrade implosion symmetry and fusion yields. Experimental yields on the OMEGA Laser System, however, have shown much less sensitivity to this initial target offset. This paper presents simulations of offset ICF implosions improved by including a physics model of cross-beam energy transfer (CBET), a mechanism of laser energy scattering from one beam to another. Room-temperature OMEGA implosion experiments with prescribed target offsets are simulated with and without CBET, illustrating that CBET mitigates the ℓ = 1 implosion asymmetry from target offset. Comparison of simulations to multiple complementary experimental observables indicates the addition of CBET physics in offset simulations is necessary to match experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accelerating amorphous polymer electrolyte screening by learning to reduce errors in molecular dynamics simulated properties

Polymer electrolytes are promising candidates for the next generation lithium-ion battery technology. Large scale screening of polymer electrolytes is hindered by the significant cost of molecular dynamics (MD) simulation in amorphous systems: the amorphous structure of polymers requires multiple, repeated sampling to reduce noise and the slow relaxation requires long simulation time for convergence. Here, we accelerate the screening with a multi-task graph neural network that learns from a large amount of noisy, unconverged, short MD data and a small number of converged, long MD data. We achieve accurate predictions of 4 different converged properties and screen a space of 6247 polymers that is orders of magnitude larger than previous computational studies. Further, we extract several design principles for polymer electrolytes and provide an open dataset for the community. Our approach could be applicable to a broad class of material discovery problems that involve the simulation of complex, amorphous materials.

36 MATERIALS SCIENCE↗

Scaling of Turbulent Viscosity and Resistivity: Extracting a Scale-dependent Turbulent Magnetic Prandtl Number

Turbulent viscosity νt and resistivity ηt are perhaps the simplest models for turbulent transport of angular momentum and magnetic fields, respectively. The associated turbulent magnetic Prandtl number Pr t ≡ ν t /η t has been well recognized to determine the final magnetic configuration of accretion disks. Here, we present an approach to determining these "effective transport" coefficients acting at different length scales using coarse-graining and recent results on decoupled kinetic and magnetic energy cascades. By analyzing the kinetic and magnetic energy cascades from a suite of high-resolution simulations, we show that our definitions of ν t , η t , and Pr t have power-law scalings in the "decoupled range." We observe that Pr t ≈ 1–2 at the smallest inertial-inductive scales, increasing to ≈5 at the largest scales. Furthermore, based on physical considerations, our analysis suggests that Pr t has to become scale independent and of order unity in the decoupled range at sufficiently high Reynolds numbers (or grid resolution) and that the power-law scaling exponents of velocity and magnetic spectra become equal. In addition to implications for astrophysical systems, the scale-dependent turbulent transport coefficients offer a guide for large-eddy simulation modeling.

79 ASTRONOMY AND ASTROPHYSICS↗

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗

Power System Network Reduction for Power Hardware-in-the-Loop Simulation

This paper proposes single-port equivalent and two-port equivalent network reduction methods to respectively reduce single-port and two-port areas in a large power network. Parameters of the reduced systems are rigorously derived, which guarantees that the electrical quantities at the port(s) remain unchanged over the reduction, including voltage magnitude and phase and active and reactive power injections into the area to be reduced. The proposed techniques are applied to reduce a practical Maui grid, where the total numbers of buses, lines and transformers are respectively reduced from 212, 106 and 108 to 45, 30 and 13. Dynamic behaviors between the full model and the reduced model are compared in detail to illustrate the efficacy and accuracy of the proposed network reduction.

model reduction↗

Power System Network Reduction for Power Hardware-in-the-Loop Simulation: Preprint

This paper proposes a single-port equivalent and a two-port equivalent to respectively reduce single-port and two-port areas in a large power network. Parameters of the reduced systems are rigorously derived, which guarantees that the electrical quantities at the port remain unchanged over the reduction, including voltage magnitude and phase, active and reactive power injections into the area to be reduced. The proposed techniques are applied to reduce a practical Maui grid, where the total numbers of buses, lines and transformers are respectively reduced from 212, 106 and 108 to 45, 30 and 13. Dynamic behaviors between the full model and the reduced model are compared in detail to illustrate the efficacy and accuracy of the proposed network reduction.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Enhanced dynamic equivalent identification method of large-scale power systems using multiple events

The increasing complexity of the interconnected power system makes high-fidelity dynamic simulation models computationally more intensive. To improve computation efficiency, model reduction techniques have been investigated to only preserve the dynamics in a limited area of interest (study area), while deriving equivalent representation for the external area. For this purpose, a measurement-based reduction approach using system identification techniques has been previously proposed. Accordingly, external areas are represented by dynamic equivalent loads using transfer function estimation. In this paper, in order to enhance the accuracy of the reduced model in preserving dynamics of the study area, multiple “grid” events of different types and at different locations are considered for identifying the parameters of the equivalent loads. Additionally, case studies are carried out in the NPCC 140-bus system. Comparisons are made between multiple events training and single event training, highlighting the advantages of the proposed method in providing a better representation of the grid dynamics under different operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Enhanced Dynamic Equivalent Identification Method of Large-Scale Power Systems Using Multiple Event

The increasing complexity of the interconnected power system makes high-fidelity dynamic simulation models computationally more intensive. To improve computation efficiency, model reduction techniques have been investigated to only preserve the dynamics in a limited area of interest (study area), while deriving equivalent representation for the external area. For this purpose, a measurement-based reduction approach using system identification techniques has been previously proposed. In that case, external areas are represented by dynamic equivalent loads using transfer function estimation. In this paper, in order to enhance the accuracy of the reduced model in preserving dynamics of the study area, multiple “grid” events of different types and at different locations are considered for identifying the parameters of the equivalent loads. Case studies are carried out in the NPCC 140-bus system. Comparison are made between multiple events training and single event training, highlighting the advantages of the proposed method in providing a better representation of the grid dynamics under different operating conditions.

Xue, Yaosuo↗

The Dynamic of Complex Two-Phase Mixtures During Coarsening: From Dendritic to Bicontinuous Mixtures

Dendrites are tree-like crystals that frequently form during the solidification of castings. In systems with nearly isotropic interfacial energy, these dendrites possess secondary and sometimes even tertiary arms, and lead to morphologically complex two-phase mixtures. Dealloying, in which one component is selectively leached, leads to nanoporous structures that have proven useful for many technological applications. These bicontinuous two-phase mixtures are topologically and morphologically complex with spatially varying mean and Gaussian curvatures. The evolution of these structures during coarsening is of considerable technological importance since the morphology of the dendrites and nanoporous structure sets the properties of the resulting material. During coarsening, the total surface area per volume of the mixture decreases, thus resulting in an increase in the size-scale of the two-phase mixture. The evolution of interfacial morphology of these interconnected two-phase mixtures is being studied using large-scale computer simulations. The simulations are used to develop a theory of coarsening in these systems that will elucidate the impact of the complicated morphology found in dendritic and nanoporous systems on the coarsening process. The experiments employ four-dimensional X-ray microtomography, which enables the morphology of the two-phase mixture to be followed in three dimensions and time (an additional dimension). The experiments also serve as a test of simulations and theory of the coarsening process, as well as to provide initial conditions for the simulations. Through these experiments and theory, we aim to develop a comprehensive description of coarsening in these complex and technologically important two-phase mixtures.

36 MATERIALS SCIENCE↗

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗

Attribution of the role of climate change in the forest fires in Sweden 2018

In this study, we analyse the role of climate change in the forestfires that raged through large parts of Sweden in the summer of 2018 from ameteorological perspective. This is done by studying the Canadian FireWeather Index (FWI) based on sub-daily data, both in reanalysis data sets(ERA-Interim, ERA5, the Japanese 55 year Reanalysis, JRA-55, and Modern-Era Retrospective analysis for Research and Applications version 2, MERRA-2) and three large-ensemble climatemodels (EC-Earth, weather@home, W@H, and Community Earth System Model, CESM) simulations. The FWI, based on reanalysis, correlates well with the observed burnt area in summer (r=0.6 to 0.8). We find that the maximum FWI in July 2018 had return times of ~24 years (90 % CI, confidence interval, > 10 years) for southern and northern Sweden. Furthermore, we find a negative trend of the FWI for southern Sweden over the 1979 to 2017 time period in the reanalyses, yielding a non-significant reduced probability of such an event. However, the short observational record, large uncertainty between the reanalysis products and large natural variability of the FWI give a large confidence interval around this number that easily includes no change, so we cannot draw robust conclusions from reanalysis data. The three large-ensembles with climate models point to a roughly 1.1 (0.9 to 1.4) times increased probability (non-significant) for such events in the current climate relative to preindustrial climate. For a future climate (2 °C warming), we find a roughly 2 (1.5 to 3) times increased probability for such events relative to the preindustrial climate. The increased fire weather risk is mainly attributed to the increase intemperature. The other main factor, i.e. precipitation during summer months, isprojected to increase for northern Sweden and decrease for southern Sweden. We, however, do not find a clear change in prolonged dry periods in summer months that could explain the increased fire weather risk in the climate models. In summary, we find a (non-significant) reduced probability of such eventsbased on reanalyses, a small (non-significant) increased probability due to global warming up to now and a more robust (significant) increase in the risk for such events in the future based on the climate models.

54 ENVIRONMENTAL SCIENCES↗

Scaling from Flux Towers to Ecosystem Models: Regional Constraints on Carbon Cycle Processes from Atmospheric Carbonyl Sulfide (Final Report)

DOE supported research suggests that gross primary productivity (GPP) is largely underestimated by global earth system models [Welp et al., 2011], reflecting the persistent challenge in extrapolating from local-scale GPP observations to global-scale earth system models. This poor understanding of GPP at large spatial scales is of particular concern in tropical forests. In tropical forests, some earth systems models forecast a powerful feedback between a warming climate and a decline in GPP resulting in forest dieback. While this simulated feedback is intensely debated, we lack robust large-scale constraints on GPP that are needed to resolve this debate. In particular, carbon dioxide measurements provide valuable information on net carbon flux, but not on the gross flux associated with GPP. Here we conducted a study of regional-to-global scale GPP using atmospheric carbonyl sulfide to provide a new constraint on GPP mechanisms in earth system models. Our project activities integrated modeling, in situ measurement, and remote sensing techniques to resolve GPP for the Amazon as well as global scale trends. The results of this work included initiating airborne carbonyl sulfide monitoring in the Amazon, training for postdocs and graduate students at a Hispanic Serving Institution, fundamental advances in carbonyl sulfide budgets [e.g. Hilton et al., Nature Climate Change, 2017], and high-profile publications that focused on GPP trends for the Amazon [Stinecipher et al., GRL, 2022] and global historical GPP trends [Campbell, et al., Nature, 2017]. Based on the suggestion of our DOE program manager, we published a state-of-the-science commentary to the scientific community on GPP monitoring with COS [Campbell et al., EOS, 2017] which was selected as the cover story. DOE support was acknowledged in all reports. The importance of this research to understanding climate change was communicated to the general public through community seminars (Rotary, Public Libraries, State Parks), an op-ed (SF Chronicle), and interviews in the mass media including two stories in the New York Times (4/5/17; 7/30/18), one of which was especially widely read after it was featured in the New York Time’s Quote of the Day.

54 ENVIRONMENTAL SCIENCES↗

Structural Simluation Toolkit (SST) v.12.0

The Structural Simulation Toolkit (SST) was developed to explore innovations in highly concurrent computing systems where the instruction set architecture (ISA), micro-architecture, and memory interact with the programming model and communications system. The package provides a fully modular design for extensive exploration of an individual system parameter as well as a parallel simulation environment based on message passing interface (MPI) which enable a high level of performance as well as the ability to look at large systems. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Rodrigues, ArunF.↗

High Fidelity CFD Simulations Supporting the KP-FHR

Kairos Power, LLC, is developing its version of the Fluoride-cooled High-temperature Reactor, the KP-FHR. The design uses a pebble bed core with fluoride salt as a coolant. The pebbles used in the KP-FHR have a diameter of 4 cm, with a shell fuel region where TRISO particles are embedded. A Pebble bed core design is adopted by several Gen IV reactors, They boast many benefits, such as fuel integrity, highly efficient heat transfer, and passive safety. However, it is challenging to accurately predict temperature and flow inside a pebble bed. Traditional approaches use the porous media model, which regards the pebble bed as a continuous medium, but with different temperature fields representing different levels, such as the fluid temperature, pebble surface temperature, and pebble center temperature. Empirical heat transfer correlations are adopted to calculate the heat transfer coefficient between different phases. However, empirical correlations are usually validated with experimental data, which usually lacks detail inside the pebble bed. The available experimental data is also generally at a high Reynolds number, which falls outside of the conditions of KP-FHR. Explicit computational fluid dynamics (CFD) simulations of randomly packed pebble beds have only become feasible recently. This is thanks to the rapid development of computational power and scalable algorithms. In this work, we used the Spectral Element Method (SEM) CFD code NekRS to simulate the randomly packed pebble bed in a cylindrical container. NekRS, which is the GPU variant of Nek5000, but refactored to utilize the computational power of GPUs using the OCCA library to run on hybrid architecture high performance computing systems. It was initially developed with the libParamunal library, but truncated and tuned for large-scale turbulence simulation. As a result, the SEM reaches higher precision with the same degrees of freedom by using a high-order Lagrange polynomial basis distributed on Gauss-Lobatto-Legendre quadrature inside each element, compared to lower-order methods, such the Finite Volume Method and Finite Element Method. The report is divided into five parts. We start with a general discussion of the pebble bed reactor, along with a specific investigation into the KP-FHR. The second part presents the numerical methodology. In the third part, we study a modular pebble bed with 1741 pebbles in a container of 7 pebble-diameter radius. Beyond LES simulations done by NekRS, we also leveraged the thermal radiation model in OpenFOAM to study heat transfer under no-forced-flow scenarios. Then, in the fourth part we simulated a pebble bed similar to the size of the Hermes Test Reactor. The total number of pebbles is in these simulations is 34,374. The container radius is 14 pebble-diameters. Finally, the report concludes in part five, with a discussion of future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Scalable Transmission and Distribution Co-simulation Platform for IBR-heavy Power Systems

The integration of Inverter-Based Resources (IBRs) into power systems, including both transmission and distribution networks, poses challenges for studying grid dynamic behaviours under various operating conditions and system events. This paper addresses the urgent need to investigate the influence of IBRs on grid dynamics, and their potential to enhance power system reliability and resilience. We propose a flexible, scalable transmission and distribution co-simulation platform, using opensource tools only for assessing the impact of grid-following (GFL) and grid-forming (GFM) IBRs on dynamic stability at various renewable penetration levels (up to 100%). This platform enables researchers to explore different contingencies at transmission, distribution, or both, providing a comprehensive evaluation of grid status. A series of case studies, including both small and large T&D systems with varying GFL/GFM configurations and contingencies, have been conducted. The results not only robustly validate our co-simulation framework, but also provide invaluable insights for effective IBR management in the power grid.

grid management↗

Atmospheric teleconnection associated with the Atlantic multidecadal variability in summer: assessment of the CESM1 model

Observational analysis shows that the Atlantic multidecadal variability (AMV) is associated with climate variability in the Northern Hemisphere through a zonal atmospheric teleconnection extending from the North Atlantic Ocean and propagating eastward around the Northern Hemisphere. Here, we studied the fidelity of model simulations in reproducing the observed summer AMV and the associated impacts on the mid-latitude climate by analysing simulations using the National Centre for Atmospheric Research Community Earth System Model Version 1 (CESM1), including CESM1 North Atlantic idealized and pacemaker simulations, CESM1 large ensemble twentieth century uninitialized simulations and large ensemble initialized CESM1 decadal predictions. To further compare the fidelity of CESM1, we also analysed large ensemble simulations from three other models. Our results suggest that the uninitialized large ensemble simulations from all models can produce an AMV time evolution and its regional climate impacts similar to the observations to certain degree. By initializing the observed oceanic condition in decadal prediction simulations, the simulated AMV and its regional impacts are closer to the observed ones than those in uninitialized ensemble simulations. In addition, the pacemaker simulations that nudged the time-evolving observed North Atlantic sea surface temperature anomalies produce spatiotemporal characteristics of the AMV and AMV climate impacts closer to the observed ones than the uninitialized simulations. We conclude that although coupled models can produce AMV and its regional impacts similar to observed, proper initialization and bias correction of the sea surface temperature spatial and temporal structure can improve this capability.

Initialization↗

Distributed optimization and control for autonomous distributed energy resource power dispatch and frequency regulation considering communication failures

Distributed energy resources (DERs) can significantly enhance grid flexibility by supplying power and secondary frequency regulation (SFR) services. Yet, their large number and dispersed nature limit the effectiveness of traditional centralized power dispatch and frequency control schemes, and expose the system to single-point communication failures. This paper presents a unified distributed optimization and control framework for DER power dispatch and SFR provision. A transmission-loss-aware distributed direct current optimal power flow is first formulated to optimize DER power dispatch and reserve scheduling. On top of this, a fully distributed control scheme for DER autonomous SFR provision is developed to allocate automatic generation control signals cost-effectively while ensuring robustness to communication link disruptions. The framework is validated on IEEE 14-bus and 39-bus transmission systems with large-scale distribution feeders using transmission-distribution dynamic co-simulation, demonstrating reliable and autonomous DER power dispatch and SFR provision under diverse operating conditions and communication contingencies.

24 POWER TRANSMISSION AND DISTRIBUTION↗