Search NASA⌕ Search

SEARCH · Search NASA

Results for “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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cloud droplet diffusional growth in homogeneous isotropic turbulence: bin microphysics versus Lagrangian super-droplet simulations

The increase in the spectral width of an initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated by applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or a Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase in the spectral width. We apply a large computational domain (64 3 m 3 ), one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence integral timescale is much larger than the droplet phase relaxation timescale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 µm is only sufficient for the lowest droplet concentration (26 cm -3 ). For the highest droplet concentration (650 cm -3 ), an order of magnitude smaller bin size is barely sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high-TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and a relatively small number of Lagrangian particles that are used in the particle-based microphysics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Evaluation of interactive and prescribed agricultural ammonia emissions for simulating atmospheric composition in CAM-chem

Abstract. Ammonia (NH3) plays a central role in the chemistry of inorganic secondary aerosols in the atmosphere. The largest emission sector for NH3 is agriculture, where NH3 is volatilized from livestock wastes and fertilized soils. Although the NH3 volatilization from soils is driven by the soil temperature and moisture, many atmospheric chemistry models prescribe the emission using yearly emission inventories and climatological seasonal variations. Here we evaluate an alternative approach where the NH3 emissions from agriculture are simulated interactively using the process model FANv2 (Flow of Agricultural Nitrogen, version 2) coupled to the Community Atmospheric Model with Chemistry (CAM-chem). We run a set of 6-year global simulations using the NH3 emission from FANv2 and three global emission inventories (EDGAR, CEDS and HTAP) and evaluate the model performance using a global set of multi-component (atmospheric NH3 and NH4+, and NH4+ wet deposition) in situ observations. Over East Asia, Europe and North America, the simulations with different emissions perform similarly when compared with the observed geographical patterns. The seasonal distributions of NH3 emissions differ between the inventories, and the comparison to observations suggests that both FANv2 and the inventories would benefit from more realistic timing of fertilizer applications. The largest differences between the simulations occur over data-scarce regions. In Africa, the emissions simulated by FANv2 are 200 %–300 % higher than in the inventories, and the available in situ observations from western and central Africa, as well as NH3 retrievals from the Infrared Atmospheric Sounding Interferometer (IASI) instrument, are consistent with the higher NH3 emissions as simulated by FANv2. Overall, in simulating ammonia and ammonium concentrations over regions with detailed regional emission inventories, the inventories based on these details (HTAP, CEDS) capture the atmospheric concentrations and their seasonal variability the best. However these inventories cannot capture the impact of meteorological variability on the emissions, nor can these inventories couple the emissions to the biogeochemical cycles and their changes with climate drivers. Finally, we show with sensitivity experiments that the simulated time-averaged nitrate concentration in air is sensitive to the temporal resolution of the NH3 emissions. Over the CASTNET monitoring network covering the US, resolving the NH3 emissions hourly instead monthly reduced the positive model bias from approximately 80 % to 60 % of the observed yearly mean nitrate concentration. This suggests that some of the commonly reported overestimation of aerosol nitrate over the US may be related to unresolved temporal variability in the NH3 emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Long-term variability in immersion-mode marine ice-nucleating particles from climate model simulations and observations

Abstract. Ice-nucleating particles (INPs) in the Southern Ocean (SO) atmosphere have significant impacts on cloud radiative and microphysical properties. Yet, INP prediction skill in climate models remains poorly understood, in part because of the lack of long-term measurements. Here we show, for the first time, how model-simulated INP concentrations compare with year-round INP measurements during the Macquarie Island Cloud Radiation Experiment (MICRE) campaign from 2017–2018. We simulate immersion-mode INP concentrations using the Energy Exascale Earth System Model version 1 (E3SMv1) by combining simulated aerosols with recently developed deterministic INP parameterizations and the native classical nucleation theory (CNT) for mineral dust in E3SMv1. Because MICRE did not collect aerosol measurements of super-micron particles, which are more effective ice nucleators, we evaluate the model's aerosol fields at other high-latitude sites using long-term in situ observations of dust and sea spray aerosol. We find that the model underestimates dust and overestimates sea spray aerosol concentrations by 1 to 2 orders of magnitude for most of the high-latitude sites in the Southern Hemisphere. We next compare predicted INP concentrations with concentrations of INPs collected on filter samples (typically for 2 or 3 d) and processed offline using the Colorado State University ice spectrometer (IS) in immersion freezing mode. We find that when deterministic parameterizations for both dust and sea spray INPs are used, simulated INPs are within a factor of 10 of observed INPs more than 60 % of the time during summer. Our results also indicate that the E3SM's current treatment of mineral dust immersion freezing in the SO is impacted by compensating biases – an underprediction of dust amount was compensated by an overprediction of its effectiveness as INPs. We also perform idealized droplet freezing experiments to quantify the implications of the time-dependent behavior assumed by the E3SM's CNT-parameterization and compare with the ice spectrometer observations. We find that the E3SM CNT 10 s diagnostic used in this study is a reasonable approximation of the exact formulation of CNT, when applied to ice spectrometer measurements in low-INP conditions similar to Macquarie Island. However, the linearized 10 s diagnostic underestimates the exact formula by an order of magnitude or more in places with high-INP conditions like the Sahara. Overall, our findings suggest that it is important to correct the biases in E3SM's simulated dust life cycle and update E3SM's INP parameterizations. INP prediction errors of 2 to 3 orders of magnitude can have considerable impacts on the simulated cloud and radiative properties in global climate models. On comparing INP concentrations during MICRE against ship-based campaigns, Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) and Antarctic Circumnavigation Expedition (ACE), we find that INPs from the latter are significantly higher only in regions closer to Macquarie Island. This alludes to the fact that physical, chemical and biological processes affecting INP concentrations as stimulated by the island could be partly responsible for the high INP concentrations observed at Macquarie Island during the MICRE campaign. Therefore, improvements to both aerosol simulation and INP parameterizations are required to adequately simulate INPs and their cloud impacts in E3SM. It will be helpful to include a parallel measurement of the size-resolved aerosol composition and explore opportunities for long-term measurement platforms in future field campaigns studying INP sources in remote marine regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Incorporation of aerosol into the COSPv2 satellite lidar simulator for climate model evaluation

Abstract. Atmospheric aerosol has substantial impacts on climate, air quality and biogeochemical cycles, and its concentrations are highly variable in space and time. A key variability to evaluate within models that simulate aerosol is the vertical distribution, which influences atmospheric heating profiles and aerosol–cloud interactions, to help constrain aerosol residence time and to better represent the magnitude of simulated impacts. To ensure a consistent comparison between modeled and observed vertical distribution of aerosol, we implemented an aerosol lidar simulator within the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package version 2 (COSPv2). We assessed the attenuated total backscattered (ATB) signal and the backscatter ratios (SRs) at 532 nm in the U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1). The simulator performs the computations at the same vertical resolution as the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), making use of aerosol optics from the E3SMv1 model as inputs and assuming that aerosol is uniformly distributed horizontally within each model grid box. The simulator applies a cloud masking and an aerosol detection threshold to obtain the ATB and SR profiles that would be observed above clouds by CALIOP with its aerosol detection capability. Our analysis shows that the aerosol distribution simulated at a seasonal timescale is generally in good agreement with observations. Over the Southern Ocean, however, the model does not produce the SR maximum as observed in the real world. Comparison between clear-sky and all-sky SRs shows little differences, indicating that the cloud screening by potentially incorrect model clouds does not affect the mean aerosol signal averaged over a season. This indicates that the differences between observed and simulated SR values are due not to sampling errors, but to deficiencies in the representation of aerosol in models. Finally, we highlight the need for future applications of lidar observations at multiple wavelengths to provide insights into aerosol properties and distribution and their representation in Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

How well are hazards associated with derechos reproduced in regional climate simulations?

Abstract. A 15-member ensemble of convection-permitting regional simulations of the fast-moving and destructive derecho of 29–30 June 2012 that impacted the northeastern urban corridor of the USA is presented. This event generated 1100 reports of damaging winds, generated significant wind gusts over an extensive area of up to 500 000 km2, caused several fatalities, and resulted in widespread loss of electrical power. Extreme events such as this are increasingly being used within pseudo-global-warming experiments to examine the sensitivity of historical, societally important events to global climate non-stationarity and how they may evolve as a result of changing thermodynamic and dynamic contexts. As such it is important to examine the fidelity with which such events are described in hindcast experiments. The regional simulations presented herein are performed using the Weather Research and Forecasting (WRF) model. The resulting ensemble is used to explore simulation fidelity relative to observations for wind gust magnitudes, spatial scales of convection (as is manifest in high composite reflectivity, cREF), and both rainfall and hail production as a function of model configuration (microphysics parameterization, lateral boundary conditions (LBCs), start date, use of nudging, compiler choice, damping, and number of vertical levels). We also examine the degree to which each ensemble member differs with respect to key mesoscale drivers of convective systems (e.g., convective available potential energy and vertical wind shear) and critical manifestations of deep convection, e.g., vertical velocities, cold-pool generation, and how those properties relate to the correct characterization of the associated atmospheric hazards (wind gusts and hail). Use of a double-moment, seven-class scheme with number concentrations for all species (including hail and graupel) results in the greatest fidelity of model-simulated wind gusts and convective structure to the observations of this event. All ensemble members, however, fail to capture the intensity of the event in terms of the spatial extent of convection and the production of high near-surface wind gusts. We further show very high sensitivity to the LBCs employed and specifically that simulation fidelity is higher for simulations nested within ERA-Interim compared to ERA5. Excess convective available potential energy (CAPE) in all ensemble members after the derecho passage leads to excess production of convective cells, wind gusts, cREF > 40 dBZ, and precipitation during a frontal passage on the subsequent day. This event proved very challenging to forecast in real time and to reproduce in the 15-member hindcast simulation ensemble presented here. Future work could examine if simulations with other initial and lateral boundary conditions can achieve greater fidelity.

Shepherd, Tristan (ORCID:0000000186276419)↗
DOE OSTI2022

IM3/HyperFACETS Thermodynamic Global Warming (TGW) Simulation Datasets

Publication For a thorough description of the methods, see the peer-reviewed paper: Jones, A.D., Rastogi, D., Vahmani, P. et al. Continental United States climate projections based on thermodynamic modification of historical weather. Sci Data 10, 664 (2023). https://doi.org/10.1038/s41597-023-02485-5 Overview The IM3 / HyperFACETS climate simulations provide 40-year historical (1980-2019) as well as four 80-year future simulations (2020-2099) over the U.S. The future simulations are split into near (2020-2059) and far future (2060-2099) segments. The future scenarios span a range of plausible changes in future climate (both Global Circulation Model (GCM) and Representative Concentration Pathways/Shared Socioeconomic Pathway (RCP/SSP) dimensions). The simulations provide climate variables with high spatiotemporal resolution (25 hourly variables and 207 3-hourly variables at 12 km2). The datasets are generated using dynamical downscaling with the WRF (Weather Research and Forecasting) model (version 4.2.1) and therefore preserve physical consistency across variables. WRF is a state-of-the-art, fully compressible, non-hydrostatic, mesoscale numerical weather prediction model. WRF is coupled with an urban canopy model (UCM), which resolves urban surfaces. The future scenarios were developed using a thermodynamic global warming approach where past events are replayed under a range of future warming conditions. These scenarios therefore provide a perspective on potential increases in extreme event intensity, geographic scope, and duration, with previously non-extreme conditions potentially crossing new thresholds to be considered extreme by today's standards. This approach is not intended to estimate future changes in extreme event frequency that might result from changes in large-scale atmospheric dynamics. This dataset has NOT been bias corrected. A bias corrected version of selected variables is under development and will be released here when available. Scenarios Files Data for each scenario is provided in weekly NetCDF files. 25 variables are available at hourly resolution, and 207 variables are available at three-hourly resolution. Spatial resolution is 12km and spans the conterminous United States (CONUS), including some areas of Canada and Mexico, resulting in a grid of 424 by 299 cells. The spatial projection is a Lambert Conformal Conic with the following proj-string: "+proj=lcc +lat_0=40.0000076293945 +lon_0=-97 +lat_1=30 +lat_2=45 +x_0=0 +y_0=0 +R=6370000 +units=m +no_defs". The available scenarios and simulation periods are listed below: historical | 1980 - 2019 rcp45cooler | 2020 - 2059 rcp45cooler | 2060 - 2099 rcp45hotter | 2020 - 2059 rcp45hotter | 2060 - 2099 rcp85cooler | 2020 - 2059 rcp85cooler | 2060 - 2099 rcp85hotter | 2020 - 2059 rcp85hotter | 2060 - 2099 * The first year (1979, 2019, and 2059) of data within each scenario represents a model warmup period and should not be used. These are located in the `spinup_files` directory. Historical year 2020 is considered an extra year of data beyond the simulation period and can be found in the `additional_files` directory. For information on specific variables and a more in-depth discussion of methodology, please refer to the data landing page at https://tgw-data.msdlive.org. Delta Warming Files The global and CONUS warming deltas for each scenario are provided in degrees Celsius annually and monthly. Restart Files Yearly restart files are provided for each scenario which can be used to restart the WRF model at a particular point in time. Spinup Files The first year of data within each simulation period represents a model warmup period and should not be used. The files are provided here for the sake of reproducibility. Additional Files Additional years of data are provided as an extension of the historic simulation.

Jones, Andrew D.↗
DOE OSTI2017

Real-Time Digital Simulation of Inertial Response with Hardware-in-the-Loop Implementation on the CART3 Wind Turbine at the National Wind Technology Center

With increasing penetrations of wind power on electric grids, the stability and reliability of interconnected power systems may be impacted. In some countries that have developed renewable energy sources and systems, grid codes have been revised to require wind power plants (WPPs) to provide ancillary services to support the power system frequency in case of severe grid events. To do this, wind turbine generators (WTGs) should be deloaded to reserve a certain amount of active power for primary frequency response; however, deloading curtails annual energy production, and the market for this type of service needs to be further developed. In this report, we focus on the temporary frequency support provided by WTGs through inertial response. WTGs have potential to provide inertial response, but appropriate control methods should be implemented. With the implemented inertial control methods, wind turbines are capable of increasing their active power output by releasing some of their stored kinetic energy when a frequency excursion occurs. Active power can be temporarily boosted above the maximum power points, after which the rotor speed decelerates, and subsequently an active power output reduction restores the kinetic energy. In this report, we develop two types of models for wind power systems: the first is common, based on the wind power aerodynamic equation, and the power coefficient can be regressed using nonlinear functions; the second is much more complicated, wherein the wind turbine system is modeled using the Fatigue, Aerodynamics, Structures, and Turbulence Modeling (FAST) tool with several degrees of freedoms. A nine-bus test power system is built in Simulink and the Real-Time Digital Simulator, respectively, which are used to evaluate the frequency support performance of the WPPs. We implement two distinct types of inertial control methods in the modeled wind turbines: frequency-based inertial control (FBIC) and stepwise inertial control (SIC). We compare the performances of the two methods in terms of their frequency nadirs, rates of change of frequency, and recovery times. We conclude the results under various wind speeds and penetration cases, which provide insight into designing the inertial response of WTGs. Further, we discuss the impact of the parameters on the performance of the inertial control methods. We evaluate both the scaling factors for the FBIC method and the slope values for the TLIC methods. The simulation work shows the characteristics of different inertial responses compared to conventional synchronous generators. Based on the simulation results, we modify, improve, and test the inertial control methods under a more realistic wind turbine model based on FAST. We then validate the inertial responses under highly turbulent wind conditions generated by TurbSim, and we examine their influences on the turbine mechanical components. The extensive simulation proves the effectiveness of the proposed inertial control methods as well as the nine-bus test power system. We then reconsider the parameters. We rebuild the same test power system using Real time Simulator Computer Aided Design (RSCAD), and we implement the inertial control methods in the real Controls Advanced Research Turbine (CART3), which is prepared for the hardware-in-the-loop field-test simulation. After the setups for the hardware and software hybrid simulation platform are complete, the inertial response is further tested on a real wind turbine for the first time, in which CART3 release the controlled inertial response against the emulated frequency excursion, provided by the real-time simulated power system test bed in RTDS.

17 WIND ENERGY↗
DOE OSTI2017

Simulating, visualizing, and searching traffic in a software defined network

The present disclosure pertains to systems and methods for simulating data packet routing within a software defined network (“SDN”), visualizing the results of the simulation, and permitting a user to search the resulting simulation. In one specified embodiment, a system may receive from a user a simulation parameter associated with a packet to be simulated in the SDN. A packet based on the at least one simulation parameter may be generated. A response of the SDN to the packet may be simulated by identifying applicable traffic routing rules and identifying a subsequent destination based on the applicable traffic routing rules. A record of the subsequent destination may be added to the simulation result, and the process may continue until a terminating condition is satisfied.

42 ENGINEERING↗
DOE OSTI2020

Hybrid EMT-TS Simulation Strategies to Study High Bandwidth MMC-Based HVdc Systems

Modular multilevel converters (MMCs) are widely used in the design of modern high-voltage direct current (HVdc) transmission system. High-fidelity dynamic models of MMCs-based HVdc system require small simulation time step and can be accurately modeled in electro-magnetic transient (EMT) simulation programs. The EMT program exhibits slow simulation speed and limitation on the size of the model and brings certain challenges to test the high-fidelity HVdc model in system-level simulations. This paper presents the design and implementation of a hybrid simulation framework, which enables the co-simulation of the EMT model of Atlanta-Orlando HVdc line and the transient stability (TS) model of the entire Eastern Interconnection system. This paper also introduces the implementation of two high-fidelity HVdc line models simulated at different time steps and discusses a dedicated method for sizing the buffer areas on both sides of the HVdc line. The simulation results of the two HVdc models with different sizes of buffer areas are presented and compared.

Liu, Yuan↗
DOE OSTI2020

Real-time simulator and controller of power system using distributed data streaming server

Systems and methods for simulating and controlling a power system in real time, using a controller, are provided. The controller includes a simulation layer to simulate an operation of the power system, a disturbance generation layer to provide data to the simulation layer to disturb the simulated operation of the power system, an application layer to display the simulated operation of the power system and generate a control signal to control, based on the simulated operation of the power system, at least one element of the power system, and a distributed data streaming server (D2S2) to allow interoperability among the simulation layer, the disturbance layer, and the application layer.

Li, Fangxing↗
DOE OSTI2021

Comparison of Electromagnetic Transient and Phasor-based Simulation for the Stability of Grid-Forming-Inverter-based Microgrids

The increasing interest in studying feeder-level, inverter-based microgrids requires appropriate simulation methods that can effectively simulate the dynamics of distribution feeders that typically have thousands of nodes. To examine the accuracy of phasor-based simulation for the stability of feeder-level, inverter-based microgrids, a comparative study between electromagnetic transient (EMT) and phasor-based simulation was conducted. Case studies of a two-inverter small microgrid and a modified version of the IEEE 34-Node Test System that has seven grid-forming inverters show that phasor-based simulation agrees well with EMT simulation when the P-f droop gains of grid-forming inverter controllers are within normal ranges. However, under extremely high P-f droop gains, phasor-based simulation can result in over-optimistic stability results. This study provides insight into the applicability of phasor-based simulation for the stability of feeder-level, inverter-based microgrids.

Du, Wei↗
DOE OSTI2021

Incorporation of Thermal Hydraulic Models for Thermal Power Dispatch into a PWR Power Plant Simulator

This report describes the development, modeling, and results of a generic pressurized water reactor power plant simulator that incorporates coupled electrical and thermal power dispatch to an industrial process located approximately one kilometer from the nuclear power plant. The simulator is a commercial PWR simulator that has been modified to include thermal power dispatch as described in past milestone reports [ , ]. The commercial PWR simulator is a generic simulator available from GSE SYSTEMS® (Sykesville, MD, USA) that is built using RELAP5-HDTM Real-Time Solution and in-house software developed by GSE Systems. This generic PWR (GPWR) simulator performs real-time simulation of the complete power plant from the reactor neutronics to the electricity generation and distribution. All primary, secondary, and auxiliary systems are modeled including all control logic in order to provide the most accurate representation of actual nuclear power plant (NPP) operation, and the simulator results have been rigorously verified by an actual NPP operating at approximately 1 GWe. This report is a continuation of worked performed in previous years, and supplemental information from previous reports is included in the appendix for reference.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Snowmass2021 Computational Frontier White Paper: Cosmological Simulations and Modeling

Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the "Cosmic Frontier", which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will be leveraged to present new classes of cosmic probes as well as to minimize systematic biases present in individual surveys. Success of this observational program requires actively pairing it with a well-matched state-of-the-art simulation and modeling effort. Next-generation cosmological modeling will increasingly focus on physically rich simulations able to model outputs of sky surveys spanning multiple wavebands. These simulations will have unprecedented resolution, volume coverage, and must deliver guaranteed high-fidelity results for individual surveys as well as for the cross-correlations across different surveys. The needed advances are as follows: (1) Development of scientifically rich and broadly-scoped simulations, which capture the relevant physics and correlations between probes (2) Accurate translation of simulation results into realistic image or spectral data to be directly compared with observations (3) Improved emulators and/or data-driven methods serving as surrogates for expensive simulations, constructed from a finite set of full-physics simulations (4) Detailed and transparent verification and validation programs for both simulations and analysis tools. (Abridged)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Efficient Subset Simulation using Hamiltonian Neural Network enhanced Markov Chain Monte Carlo Methods

The Monte Carlo method delivers an unbiased estimate of the probability of failure. However, the variance of the estimate depends on the number of evaluated samples. This number must be very large for estimations of a low probability of failure. If the evaluation of each sample is computationally expensive, the crude Monte Carlo simulation strategy is impracticable. Therefore, subset simulations are used to reduce the required number of evaluations. Subset simulations require a Markov Chain Monte Carlo sampler, such as the random walk Metropolis-Hastings algorithm. The algorithm, however, struggles with sampling in low-probability regions, especially if they are narrow. As a consequence, advanced Markov Chain Monte Carlo simulations have been developed. In particular, the Hamiltonian Monte Carlo method explores the target distribution rapidly. Driven by the idea of Hamiltonian dynamics, this sampler provides a non-random walk through the target distribution. The incorporation of subset simulation and Hamiltonian Monte Carlo methods has shown promising results for reliability analysis. One downside of the Hamiltonian Monte Carlo method is that gradient evaluations are computationally expensive, especially when dealing with high-dimensional problems and evaluating long trajectories. We show that integrating Hamiltonian neural networks in Hamiltonian Monte Carlo simulations significantly speeds up the sampling task. Furthermore, the enhancement of adaptive trajectory length within the Hamiltonian Monte Carlo results in the efficient proposal of the following states. Based on this recent enhancement, we provide a fast sampling strategy for subset simulations using Hamiltonian neural networks to replace the evaluation of the gradient and significantly speed up the Hamiltonian Monte Carlo simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Buffer-IPyC separation process in TRISO fuel particles simulated with Bison code

During High Temperature Gas-cooled Reactor (HTGR) operation, tristructural isotropic (TRISO) coated-particle fuel undergoes irradiation-induced changes in morphology and thermomechanical properties. Experimental results from the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program show, among other things, the mechanism of gap formation between the buffer and inner pyrolytic carbon (IPyC) layers, which could be explored further via computational simulations using the Bison code. Two simulation models were developed, the debonding restricted model, where no gap formation between buffer and IPyC layers is permitted, and the debonding enabled model, where the gap between those layers is created. The inputs of the simulated models are based on the irradiation conditions from the AGR-1 experiment. The research included simulations on spherical and aspherical fuel types. Under the specific temperatures and fluences of the AGR-1 irradiation experiment, and the Bison simulations, it was concluded that the most common scenario is a gap formation along the buffer-IPyC interface, while the least possible scenario is the situation where there is no gap formation at the buffer-IPyC junction. The computational results confirmed that the sphericity of the fuel influences the thickness of the gap that occurs at the buffer-IPyC junction, in a way that with increasing aspect ratio the gap thickness increases. The results obtained for spherical and aspherical fuel are nearly identical. Finally, performed simulations match conclusions observed from the AGR-1 experiment, which as such shows that the Bison code is a good computational method for simulating the TRISO fuel. Future simulations will include the validation of performed research and comparison of the results between Bison and PARFUME codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

The Importance of Being Adaptable: An Exploration of the Power and Limitations of Domain Adaptation for Simulation-Based Inference with Galaxy Clusters

The application of deep machine learning methods in astronomy has exploded in the last decade, with new models showing remarkably improved performance on benchmark tasks. Not nearly enough attention is given to understanding the models' robustness, especially when the test data are systematically different from the training data, or "out of domain." Domain shift poses a significant challenge for simulation-based inference, where models are trained on simulated data but applied to real observational data. In this paper, we explore domain shift and test domain adaptation methods for a specific scientific case: simulation-based inference for estimating galaxy cluster masses from X-ray profiles. We build datasets to mimic simulation-based inference: a training set from the Magneticum simulation, a scatter-augmented training set to capture uncertainties in scaling relations, and a test set derived from the IllustrisTNG simulation. We demonstrate that the Test Set is out of domain in subtle ways that would be difficult to detect without careful analysis. We apply three deep learning methods: a standard neural network (NN), a neural network trained on the scatter-augmented input catalogs, and a Deep Reconstruction-Regression Network (DRRN), a semi-supervised deep model engineered to address domain shift. Although the NN improves results by 17% in the Training Data, it performs 40% worse on the out-of-domain Test Set. Surprisingly, the Scatter-Augmented Neural Network (SANN) performs similarly. While the DRRN is successful in mapping the training and Test Data onto the same latent space, it consistently underperforms compared to a straightforward Yx scaling relation. These results serve as a warning that simulation-based inference must be handled with extreme care, as subtle differences between training simulations and observational data can lead to unforeseen biases creeping into the results.

Ntampaka, Michelle [Baltimore, Space Telescope Sci↗
DOE OSTI2016

Large-scale large eddy simulation of nuclear reactor flows: Issues and perspectives

Numerical simulation has been an intrinsic part of nuclear engineering research since its inception. In recent years a transition is occurring toward predictive, first-principle-based tools such as computational fluid dynamics. Even with the advent of petascale computing, however, such tools still have significant limitations. In the present work some of these issues, and in particular the presence of massive multiscale separation, are discussed, as well as some of the research conducted to mitigate them. Petascale simulations at high fidelity (large eddy simulation/direct numerical simulation) were conducted with the massively parallel spectral element code Nek5000 on a series of representative problems. These simulations shed light on the requirements of several types of simulation: (1) axial flow around fuel rods, with particular attention to wall effects; (2) natural convection in the primary vessel; and (3) flow in a rod bundle in the presence of spacing devices. Finally, the focus of the work presented here is on the lessons learned and the requirements to perform these simulations at exascale. Additional physical insight gained from these simulations is also emphasized.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Python Library for Monte Carlo Simulations with Ab Initio and Machine-Learned Interatomic Potentials

There is a growing need in the simulation community for software that provides a transparent, reproducible, usable, and extensible (TRUE) Monte Carlo (MC) simulation framework employing energies from ab initio methods and machine-learning interatomic potentials (MLIPs). We introduce a Python library (ASE-MC) that adds Monte Carlo functionality to the Atomic Simulation Environment (ASE) package. Now, we can combine the powerful tools used to build systems and perform ab initio and MLIP in ASE with MC simulation algorithms to sample the configurational space with a concise Python script. After presenting the design philosophy, we demonstrate the flexibility of our approach using selected examples. These example simulations include liquid water described with a message-passing MLIP in the canonical and isothermal–isobaric ensembles, sampling the characteristic dihedral angle of biphenyl and comparing an MLIP to first-principles calculations, and a grand canonical Monte Carlo simulation of ammonia adsorption on Pt(111). These examples showcase the main features of the software, which include flexibility in the choice of ab initio or MLIP engine, ab initio or MLIP grand canonical MC with cavity bias insertions and deletions, the ability to add custom MC moves to the move set, and how users can condense complex MC workflows into a single Python script. Finally, this library serves as a framework for reproducible Monte Carlo simulations, facilitating easy reproduction of the work and application to new systems.

97 MATHEMATICS AND COMPUTING↗