Search NASA⌕ Search

SEARCH · Search NASA

Results for “Performance Rating Method”

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,315 records · Page 73

Catalytic direct conversion of ethane to value-added chemicals under microwave irradiation

The reaction mechanism of ethane conversion over 4 % Mo-0.5 % Fe/ZSM-5 catalyst was investigated using two different reactors with different heating modes. A microwave fixed-bed reactor (MWFB) was used and its performance was compared to a conventional thermal fixed-bed reactor (CTFB). Over 80 % ethane conversion was achieved under microwave irradiation mode at 375 °C. Here, to achieve the same ethane conversion conventionally, the reaction temperature had to be elevated to 660 °C in the CTFB. The formation rates and selectivity of the products were compared under these two heating methods. The CTFB products had a higher selectivity towards aromatic compounds; whereas, the MWFB products had a higher selectivity towards ethylene; 26.5 % higher than the CTFB. Carbon analysis over spent catalysts of CTFB showed three distinct types of coke, whereas the microwave resulted in a single higher temperature coke type.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multifunctional transitional metal-based phosphide nanoparticles towards improved polysulfide confinement and redox kinetics for highly stable lithium-sulfur batteries

The shuttle effect and the sluggish redox kinetics of lithium polysulfides (LiPSs) are the major issues impeding the practical applications of lithium-sulfur batteries (LSBs). Herein, a highly-efficient Ni 2 P electrocatalyst supported on N, P co-doped graphene (Ni 2 P@NPG) is developed via a simple “recrystallization-self-assembly” method to address the above issues. Here, the ultrafine Ni 2 P nanoparticles ensure abundant adsorption-diffusion-conversion interfaces for accelerating LiPSs transformation and Li 2 S deposition, which extremely decreases the accumulation of LiPSs in the electrolyte and therefore prevents the migration of LiPSs. Their superior catalytic performance is demonstrated by reduced Gibbs free energy changes of rate-limiting step based on the systematic theoretical calculations and the reduced shuttle effect is tested by the three-dimensional reconstructions of Raman depth profiles. Benefiting from these synergistic effects, the LSBs with Ni 2 P@NPG modified separators present a superior cycling performance with an average capacity decay rate of 0.048 % per cycle at 1C around the 400 cycles and a high-rate capacity of 731 mAh/g at 2C. Even with a high-sulfur loading of 3.53 mg cm –2 , the cell can still contain a reversible capacity of 809 mAh/g at 0.2C with a remarkable columbic efficiency of 98.4 %.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A signal compensation based cascaded PI control for an industrial heat exchange system

Since the characteristics of the industrial heat exchange system (HES) are drastically affected by variations of the steam pressure, ambient temperature and water quality, it is difficult for the classical cascade PI control system to meet the desirable control objectives. Therefore, in this paper, a deterministic linear model with unmodeled dynamics is used to describe the dynamic characteristics of the HES. It should be noted that the unmodeled dynamics at the previously sampled instant can be calculated using actual measurements and that the tracking error, which includes effects of the rate of changes of unmodeled dynamics on the controlled model, can be accurately determined. An additional compensating signal is obtained with the feedforward control together with a one-step optimal control law. Consequently, the cascade PI control strategy with the signal compensation is established. Furthermore, the simulation and industrial application show that the proposed method can effectively improve the performance of the cascade PI control system, especially when the unmodeled dynamics is subjected to unknown variations.

42 ENGINEERING↗

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗

Mechanical and electrochemical response of a LiCoO 2 cathode using reconstructed microstructures

As LiCoO 2 cathodes are charged, delithiation of the LiCoO 2 active material leads to an increase in the lattice spacing, causing swelling of the particles. When these particles are packed into a bicontinuous, percolated network, as is the case in a battery electrode, this swelling leads to the generation of significant mechanical stress. In this study we performed coupled electrochemical-mechanical simulations of the charging of a LiCoO 2 cathode in order to elucidate the mechanisms of stress generation and the effect of charge rate and microstructure on these stresses. Energy dispersive spectroscopy combined with scanning electron microscopy imaging was used to create 3D reconstructions of a LiCoO 2 cathode, and the Conformal Decomposition Finite Element Method is used to automatically generate computational meshes on this reconstructed microstructure. Replacement of the ideal solution Fickian diffusion model, typically used in battery simulations, with a more general non-ideal solution model shows substantially smaller gradients of lithium within particles than is typically observed in the literature. Using this more general model, lithium gradients only appear at states of charge where the open-circuit voltage is relatively constant. While lithium gradients do affect the mechanical stress state in the particles, the maximum stresses are always found in the fully-charged state and are strongly affected by the local details of the microstructure and particle-to-particle contacts. These coupled electrochemical-mechanical simulations begin to yield insight into the partitioning of volume change between reducing pore space and macroscopically swelling the electrode. Lastly, preliminary studies that include the presence of the polymeric binder suggest that it can greatly impact stress generation and that it is an important area for future research.

25 ENERGY STORAGE↗

Rapid assessment of the creep rupture life of metals: A model enabling experimental design

Prediction of the creep rupture life of engineering metals is critical for qualification and design of new materials. The use of long-term creep tests and the need to quantify the performance variability in a priori similar systems hinder the rapid creep assessment of a given material. Therefore, it is essential to develop methods that can extrapolate the long-term performance of alloys and the associated variability from short-term experiments. To this end, this study introduces a new model which enables the estimation of the rupture life of a material for a given stress and temperature. This model relies on two components. First, a new relation for the minimum creep rate (MCR) of materials is introduced. It includes a stress dependent stress exponent allowing the model to capture the variation of MCR across a wide range of temperatures and stresses. Second, employing the Monkman-Grant (MG) law, we establish a relation between stress, temperature and creep rupture life. Together, these two elements yield a new closed-form mathematical expression for the Larson Miller parameter as a function of stress and temperature. This expression captures the creep rupture time for many metals (Gr91, Copper, Gr122 and 347H) and compares favorably with alternate empirical approaches. The model is then used to assess the minimum duration of creep rates necessary to qualify the material up to 100000h. Furthermore, it is found that depending on the material system, creep tests as few as five limited to 5000 h for steels (Gr91, Gr122, 347H) and 100 h for copper are sufficient to model creep lifetimes. Finally, using a Bayesian inference-based approach to calibrate the model, we demonstrate that variability in rupture life can be captured via the quantification of the uncertainty in the model parameters and extrapolated from a limited number of short to moderately short creep tests; thereby paving the way for accelerated creep testing.

36 MATERIALS SCIENCE↗

Strength mechanisms and tunability in Al-Ce-Mg ternary alloys enabled by additive manufacturing

Al-Ce-based alloys are promising candidates for additive manufacturing (AM) due to their hot-cracking resistance and because they do not require heat treatment to obtain precipitation strengthening. Rapid solidification rates enabled by AM methods can lead to enhanced mechanical properties; however, the strengthening mechanisms over large composition ranges were unclear. Here, combinatorial synthesis by directed-energy deposition (DED) and hardness measurements were used to rapidly map the composition-dependent strength of the ternary Al-Ce-Mg system. Tensile testing and microstructure characterization of selected compositions were performed to elucidate the compositional dependence of the strengthening mechanisms. Al 11 Ce 3 precipitates were present in all cases, and the maximum hardness (1.25 GPa) was measured for the Al-8Ce-10Mg composition. A combination of (i) Hall-Petch strengthening, based on the FCC-matrix-phase cell size; (ii) particle strengthening, based on Al 11 Ce 3 volume fraction and size; and (iii) solid-solution strengthening, based on Mg composition of the matrix phase, were used to account for the measured strengths. Vickers hardness is shown to correlate well with ultimate tensile strength in these alloys, highlighting the value of surface-based techniques for rapid screening.

36 MATERIALS SCIENCE↗

SiO 2 Is Wasted Space in Single-Nanometer-Scale Silicon Nanoparticle-Based Composite Anodes for Li-Ion Electrochemical Energy Storage

The electrode processing conditions of silicon-based composite anodes play a pivotal role in the resulting interfacial chemical speciation and, thus, the electrochemical cycling behavior of the electrode. Systematically investigating how small chemical changes to the surface of the silicon nanoparticle (NP) affect larger, electrode-level properties is a strategy that will inform design principles to maximize electrode energy density and extend electrode lifetime. Here, we incorporate silicon nanoparticles (NPs) with an average diameter of 5.5 nm synthesized from the gas phase through a nonthermal plasma method into composite anode half-cell coin cells. We perform chemistry to functionalize the native hydride-terminated silicon NP surface with N-methylpyrrolidone (NMP) to improve slurry properties and reduce the silicon reactivity prior to composite electrode fabrication. We construct composite electrodes from these silicon NPs and observe an 86% capacity retention over 100 cycles at a rate of C/5, with an initial silicon specific capacity of 2600 mAh/g. We expose the same Si NPs to water that completely oxidizes these small-diameter NPs to SiO 2 and find that the SiO 2 NPs in the same electrode configuration exhibit no obvious lithium alloying capacity in the electrochemical potential range of lithium silicide alloy formation. As this result stands in contrast to existing literature, we provide a discussion on the origin of the discrepancies.

25 ENERGY STORAGE↗

Nucleation-controlled growth of superior lead-free perovskite Cs 3 Bi 2 I 9 single-crystals for high-performance X-ray detection

The organic-inorganic hybrid lead halide perovskites have emerged as a series of star materials for solar cells, lasers and detectors. However, the issues raised by the toxic lead element and marginal stability due to the volatile organic components have severely limited their potential applications. In this work, we develop a nucleation-controlled solution method to grow large size high-quality Cs 3 Bi 2 I 9 perovskite single crystals (PSCs). Using the technique, we harvest some centimeter-sized single crystals and achieved high device performance. We find that X-ray detectors based on PSCs exhibit high sensitivity of 1652.3 μC Gy air -1 cm -2 and very low detectable dose rate of 130 nGy air s -1 , both desired in medical diagnostics. In addition, its outstanding thermal stability inspires us to develop a high temperature X-ray detector with stable response at up to 100 °C. Furthermore, PSCs exhibit high X-ray imaging capability thanks to its negligible signal drifting and extremely high stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In situ compression artifact removal in scientific data using deep transfer learning and experience replay

The massive amount of data produced during simulation on high-performance computers has grown exponentially over the past decade, exacerbating the need for streaming compression and decompression methods for efficient storage and transfer of this data---key to realizing the full potential of large-scale computational science. Lossy compression approaches such as JPEG when applied to scientific simulation data realized as a stream of images can achieve good compression rates but at the cost of introducing compression artifacts and loss of information. This paper develops a unified framework for in situ compression artifact removal in which the fully convolutional neural network architectures are combined with scalable training, transfer learning, and experience replay to achieve superior accuracy and efficiency while significantly decreasing the storage footprint as compared with the traditional optimization-based approaches. We demonstrate the proposed approach and compare it with compressed sensing postprocessing and other baseline deep learning models using climate simulations and nuclear reactor simulations, both of which are driven by hyperbolic partial differential equations. Our approach when applied to remove the compression artifacts on the JPEG-compressed nuclear reactor simulation data (using a transfer-trained model that was pretrained on the climate simulation data and updated incrementally as the nuclear reactor simulation progressed), achieved a significant improvement---mean peak signal-to-noise ratio of 42.438 as compared with 27.725 obtained with the compressed sensing approach.

97 MATHEMATICS AND COMPUTING↗

The Influence of PV Module Materials and Design on Solder Joint Thermal Fatigue Durability

Finite element model (FEM) simulations have been performed to elucidate the effect of flat plate photovoltaic (PV) module materials and design on PbSn eutectic solder joint thermal fatigue durability. The statistical method of Latin Hypercube sampling was employed to investigate the sensitivity of simulated damage to each input variable. Variables of laminate material properties and their thicknesses were investigated. Using analysis of variance, we determined that the rate of solder fatigue was most sensitive to solder layer thickness, with copper ribbon and silicon thickness being the next two most sensitive variables. By simulating both accelerated thermal cycles (ATCs) and PV cell temperature histories through two characteristic days of service, we determined that the acceleration factor between the ATC and outdoor service was independent of the variables sampled in this study. This result implies that an ATC test will represent a similar time of outdoor exposure for a wide range of module designs. This is an encouraging result for the standard ATC that must be universally applied across all modules.

14 SOLAR ENERGY↗

Multiscale Modeling of Deformation, Fracture and Failure of Fusion Materials and Structures Final Report

The main objective of the research performed at UCLA under grant number DE-FG02-03ER54708 was to develop a physical understanding of the underlying mechanisms that control high-temperature deformation and fracture of fusion structural materials. The key aspect is the synergistic interaction between thermo-mechanical damage phenomena (e.g. creep, fatigue, and fracture), and irradiation damage phenomena (helium and hydrogen bubbles, and irradiation-induced defects). Another thrust area of the research was on the development of constitutive equations for the mechanical behavior of irradiated structures in fusion components (mainly the First Wall and Plasma-Facing components), and in utilizing them in descriptions in large-scale Finite Element Modeling of fusion structures. A parallel experimental effort on thermo-mechanical and plasma-induced damage is designed to validate and confront the developed theoretical and computational predictions. The models have been experimentally validated through experimental tests in simulated plasma transients on specially-designed W samples. This grant started at UCLA on November 14, 2003 and has had successive continuations till November 14, 2017. Research performed under the grant was to develop a range of hierarchical models for the post-elastic deformation, fracture and failure of fusion structural materials. The developed multiscale modeling approach has been based on rigorous mathematical, physical, and computational methods at the forefront of computational materials science. At the fundamental and smallest length scale (nm-μm), we developed advanced rate theory and Monte Carlo approaches to model microstructure evolution, non-equilibrium phase transformations, and dislocation-defect interactions. Microscopic and mesoscopic models of radiation hardening, ductile-to-brittle-transition, post-elastic deformation, plastic instabilities, and fracture processes are based on Dislocation Dynamics (DD) and Grain Growth Dynamics (GGD). At the continuum level, we focused on modeling fracture and failure mechanisms for Virtual Integrated Testing (VISTA), and for understanding the limits of both materials-by-design versus structural component design through advancing the new approach of multiscale-multiphysics. The overall research has been strongly coupled with the U.S. national experimental program to ascertain and verify the range of investigated phenomena. In the following, we give a summary of research progress supported by this grant since its inception till November 14, 2017. Detailed research results are available in our publications, listed in the publications section of this report.

36 MATERIALS SCIENCE↗

Fundamental Corrosion Studies in High-Temperature Molten Salt Systems for Next Generation Concentrated Solar Power Systems

The objective of this project is to increase Concentrated Solar Power (CSP) system lifetime through characterization of corrosion mechanisms in heat transfer systems and identification of materials and corrosion resistance methods applicable at temperatures above 850°C. Coupling of experimental methods with thermodynamic and kinetic modeling has given greater insight into the corrosion phenomena and enabled modeling of the corrosion in molten chloride heat transfer systems. The project has been identified the selective oxidation of Cr along the grain boundaries of high Cr content alloys as the principle corrosion mechanisms occurring in CSP systems. The corrosion rate was also found to be mass transfer limited by diffusion of chromium chlorides along the grain boundaries. It was shown that the corrosion could be increased by adding naturally convective flow to mechanism using a thermosiphon. A reaction mechanism has been proposed for selective oxidation of Cr that is driven by the electrochemical potential difference between the Ni crucible for immersion testing and the alloy coupon. This driving force will be similar to the differences caused by having varying Ni content or Cr content throughout a heat transfer system. It could also mimic the differences in electrochemical potential between various parts of the system that are at different temperatures. The reaction mechanism consists of an initiation step due to system impurities and a set of propagation steps that results in formation of a NiCr alloy at the crucible that drives corrosion and selective oxidation of Cr from the alloy. The presences of chromium chlorides in the salt are also essential to the corrosion mechanism to transport the Cr. The formation of a NiCr alloy on the crucible surface with no contact between the alloy coupons and the crucible provides support to this mechanism of having corrosion reactions occurring at both the alloy and the crucible. The thermodynamics of this reaction mechanism have been modeled and were shown to be galvanic. A corrosion model has been developed based on this mechanism and has been implemented in CFD so that it can be used to predict corrosion in a wide variety of system geometries. The model was validated using the thermosiphon data and dimensionless analysis. The corrosion rates of the samples were modeled to within 10% percent of the experimental values. The model has been extended to include the effect of Mg corrosion inhibitor on the corrosion rate in the system. The use of Mg as a metallic corrosion inhibitor has been demonstrated as a method to reduce the corrosion rate of alloys in molten chlorides to below 15 microns per year. Mg is an effective corrosion inhibitor at concentrations as low as 0.3 mol%. Corrosion potential monitoring experiments have shown that the Mg lowers the corrosion potential for the alloy below the equilibrium potential for Cr oxidation and thus prevents corrosion. Experiments have been performed on the ability of Mg to inhibit corrosion in realistic environments that include lower purity salts and welded and stressed alloys. The results have shown that for all these scenarios the corrosion can be reduced to meet the DOE targets. The electrochemical monitoring of Mg for process control has also been demonstrated.

14 SOLAR ENERGY↗

Demonstration of Microwave Resonant Cavity Transducer Performance in Fluid Flow Sensing: Development of Sensor Performance Model of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave cavity-based transducer for in-core high-temperature fluid flow sensing in molten salt cooled reactors (MSCR) and sodium fast reactors (SFR). This sensor is a hollow metallic cylindrical cavity. The principle of sensing consists of making one wall of the cylindrical cavity flexible enough so that dynamic pressure, which is proportional to fluid velocity, will cause membrane deflection. Membrane deflection causes cavity volume change, which leads to a shift in the resonant frequency. To validate sensor physics, we have performed proof-of-principle test of flow sensing in water. For this test, we have developed a cylindrical resonator for K-band, which was machined from brass. The cavity was excited through WR-42 waveguide through a subwavelength hole on the side of the wall of the cylinder. To increase the spectral signal visibility, we developed a signal processing method for baseline subtraction. A flow loop for proof-of-principle test of transducer performance in water was assembled. A commercial flow meter was installed in the loop for reference measurements. Cylindrical cavity was excited in the TEM 011 mode with resonant frequency f ≈ 17.8GHz. Frequency shift of cavity spectral response was obtained by gradually increasing water flow rate from 0 to 60gpm. Corresponding monotonic increase of resonant frequency shift by several MHz was observed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Second Annual Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with the flat wall of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the resonant frequency, which can be detected with a spectrum analyzer. Because the transducer is intended for immersion in a high temperature corrosive fluid, understanding of material degradation is crucial for estimation of transducer performance lifetime, and development of measurement interpretation algorithms. We conducted a preliminary computational investigation of relevant damage mechanisms of a stainless steel 316 cylindrical resonator in FLiBe salt. The two main damage mechanisms, creep and corrosion, were modeled using multiphysics COMSOL software. Degradation was modeled for a temperature range 500°C to 700°C. Coupling of the damage mechanisms was not considered. These models predict significant inelastic deformation at most temperatures due to creep, and qualitatively predict chromium depletion both along the liquid/solid interface and along the grain boundaries. An algorithmic approach for compensation of these degradation effects during fluid flow measurements will be developed in the future work. To validate sensor physics, we have performed proof-of-principle test of flow sensing in water. For this test, we have developed a cylindrical resonator for K-band, which was machined from brass. The cavity was excited through WR-42 waveguide through a subwavelength hole on the side of the wall of the cylinder. To increase the spectral signal visibility, we developed a signal processing method for baseline subtraction. A flow loop for proof-of-principle test of transducer performance in water was assembled. A commercial flow meter was installed in the loop for reference measurements. Cylindrical cavity was excited in the TE 011 mode with resonant frequency f ≈ 17.8GHz. Frequency shift of cavity spectral response was obtained by gradually increasing water flow rate from 0 to 60gpm. Corresponding monotonic increase of resonant frequency shift by several MHz was observed. Approximate figure of merit of sensitivity to flow rate is 100KHz/GPM. In addition, we have identified an existing liquid sodium experimental setup for demonstration of flow sensing in environment similar to that of an advanced reactor. The setup consists of a cylindrical vessel and center feed line, where transducer inserted through the lid will measure velocity of the impinging liquid jet. As a calibration experiment, we have assembled a water vessel with center feed with the same dimensions as those of the liquid sodium setup. We have also developed and insertion probe consisting of a 50cm brass waveguide enclosed in protective SS316 tube. Using the water loop, we have demonstrated feasibility of sensing the impinging liquid jet in the vessel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

First Demonstration of a Liquid Argon Time Projection Chamber Based Search for Intranuclear Neutron-Antineutron Transitions and Annihilation in $^{40}$Ar using the MicroBooNE Detector

We present a novel methodology to search for intranuclear neutron-antineutron transition ($n\rightarrow\bar{n}$) followed by $\bar{n}$-nucleon annihilation within an $^{40}$Ar nucleus, using the MicroBooNE liquid argon time projection chamber (LArTPC) detector. A discovery of $n\rightarrow\bar{n}$ transition or a new best limit on the lifetime of this process would either constitute physics beyond the Standard Model or greatly constrain theories of baryogenesis, respectively. The approach presented in this paper makes use of deep learning methods to select $n\rightarrow\bar{n}$ events based on their unique features and differentiate them from cosmogenic backgrounds. The achieved signal and background efficiencies are (70.22$\pm$6.04)\% and (0.0020$\pm$0.0003)\%, respectively. A demonstration of a search is performed with a data set corresponding to an exposure of $3.32 \times10^{26}\,$neutron-years, and where the background rate is constrained through direct measurement, assuming the presence of a negligible signal. With this approach, no excess of events over the background prediction is observed, setting a demonstrative lower bound on the $n\rightarrow\bar{n}$ lifetime in $^{40}$Ar of $\tau_{\textrm{m}} \gtrsim 1.

43 PARTICLE ACCELERATORS↗

ASCR Workshop Position Paper: Challenges and Opportunities in High Energy Physics

High energy particle physics and cosmology concern themselves with estimating fundamental parameters of nature, such as the masses and interactions of fundamental particles like the Higgs boson and the rate of expansion of the universe. In doing so, they analyze exabyte-scale datasets, some of the largest in all of science, and face many challenges in subsequent data analysis. These challenges are shared between the two disciplines, but we focus on particle physics to highlight one specific domain. In particle physics, the standard method for estimating parameters involves performing Monte Carlo (MC) integration as a function of both parameters of interest and nuisance parameters using an expensive simulator, counting the number of observed collision events (i.i.d. samples) from an experiment in the corresponding integration domains, and forming a Poisson likelihood function. This likelihood function is then used in a Frequentist manner to construct a maximum likelihood point estimate (MLE) and confidence set for the parameters. To sufficiently populate the high-dimensional integration domains, simulators consume billions of CPU-hours annually and produce hundreds of petabytes of intermediate output data. Several techniques have been developed to: optimize definitions of the integration domains so as to be maximally sensitive to a particular subset of parameters, efficiently estimate the integrals, and build robust surrogate models by interpolating between integral evaluations at different parameter points. One can view this whole endeavor as classical Simulation-Based Inference (SBI).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Towards Robust Calibration of the AWSD Reactive Burn Model

Calibration of a reactive burn model for detonation of high explosive is an important step towards predictive hy drodynamic simulations of detonation. A typical calibration consists of varying model parameters (e.g., rate constants, activation energies) until results of hydrodynamic simulations match the experimental data for a certain set of ex periments. Hydrodynamic simulations of the dependence of steady detonation velocity on the radius of a cylindrical high-explosive charge - often used in such calibrations - can be computationally expensive. In this work, we propose a method where such expensive simulations are performed infrequently, and only to parameterize and refine a surrogate model for the dependence of the detonation velocity on calibrated parameters. The method is developed, implemented and applied to an example problem - calibration of the AWSD reactive burn model for important high explosive PBX 9502. Two different flavors of the surrogate model are investigated, and the calibration is performed successfully.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗