Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surrogates”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Uncertainty Quantification Framework for Predicting Material Response with Large Number of Parameters: Application to Creep Prediction in Ferritic-Martensitic Steels Using Combined Crystal Plasticity and Grain Boundary Models

This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗
DOE OSTI2024

Bayes_Opt-SWMM: A Gaussian process-based Bayesian optimization tool for real-time flood modeling with SWMM

Real-time flood model plays a pivotal role in averting urban flood damage, particularly when there is minimal lead time for preparatory measures. However, urban flood modeling in real-time often contends with inherent uncertainties arising from input data uncertainty and parameter ambiguities. Here this study introduces a real-time calibration (RTC) tool called Bayes_Opt-SWMM, specifically tailored for real-time urban flood modeling and uncertainty optimization. This tool leverages the Gaussian process-based Bayesian optimization algorithm and interfaces seamlessly with the Stormwater Management Model (SWMM). It integrates real-time model forcing data and flood monitoring collected through sensors and gauges which are strategically placed within critical locations of urban drainage systems. Our approach hinges on the Surrogate Model based Uncertainty Optimization (SMUO) concept, providing an avenue for enhancing real-time flood modeling. Bayes_Opt-SWMM runs the optimization process using a surrogate model called Gaussian Process emulator with two inference methods: (1) the Gaussian Process (GP) model and (2) Markov Chain Monte Carlo (MCMC) algorithm in GP model (GP_MCMC). Furthermore, three acquisition functions, namely Expected Improvement (EI), Maximum Probability of Improvement (MPI), and Lower Confidence Bound (LCB), facilitate optimal parameter fitting within the surrogate models. The efficiency of GP-based surrogate models in learning SWMM model parameters, leads to an improved uncertainty quantification and accelerated real-time flood modeling in urban areas. Overall, Bayes_Opt-SWMM emerges as a cost-effective and valuable tool for real-time flood modeling and monitoring, with significant potential for managing intelligent storm water systems in urban environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Exploring the Margin Recipe for Online Adaptive Radiation Therapy for Intermediate-Risk Prostate Cancer: An Intrafractional Seminal Vesicles Motion Analysis

Purpose: To provide a benchmark for seminal vesicle (SV) margin selection to account for intrafractional motion and to investigate the effectiveness of 2 motion surrogates in predicting intrafractional SV coverage. Methods and Materials: Fifteen prostate patients were studied. Each patient had 5 pairs (1 patient had 4 pairs) of pretreatment and posttreatment cone beam CTs (CBCTs). Each pair of CBCTs was registered on the basis of prostate fiducial markers. All pretreatment SVs were expanded with 1-, 2-, 3-, 4-, 5-, and 8-mm isotropic margins to form a series of planning target volumes, and their intrafractional coverage to the posttreatment SV determined the “ground truth” for exact coverage. Two motion surrogates, the center of mass (COM) and the border of contour, were evaluated by the use of Pearson product-moment correlation coefficient and exponential fitting for predicting SV underdosage. Action threshold of each surrogate was calculated. The margin for each surrogate was calculated according to a traditional margin recipe. Results: Ninety-five percent posttreatment SV coverage was achieved in 9%, 53%, 73%, 86%, 95%, and 97% of fractions with 1-, 2-, 3-, 4-, 5-, and 8-mm margins, respectively. The 5-mm margins provided 95% intrafractional SV coverage in over 90% of fractions. The correlation between the COM and border was weak, moderate, and strong in the left-right (L-R), anterior-posterior (A-P), and superior-inferior (S-I) directions, respectively. Exponential fitting gave the underdosage threshold of 4.5 and 7.0 mm for the COM and border. The Van Herk margin recipe recommended 0-, 0.5-, and 0.8-mm margins in the L-R, A-P, and S-I directions based on the COM, and 1.2-, 3.9-, and 2.5-mm margins based on the border. Conclusions: Five-millimeter isotropic margins for the SV constitute the minimum required to mitigate the intrafractional motion. Both the COM and the border are acceptable predictors for SV underdosage with 4.5- and 7.0-mm action threshold. Traditional margin based on the COM or border underestimates the margin.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Development of multi-objective core optimization framework and application to sodium-cooled fast test reactors

The optimization of a Sodium-cooled Fast Reactor (SFR) core is a challenging process, due to the large number of design parameters, the nonlinearities among inputs and outputs, and the complicated correlation among output parameters. This study attempts to develop a generalized framework for the SFR core optimization by coupling the sensitivity analysis, advanced optimization algorithm, and optionally the surrogate modeling. The framework is built based on the fast reactor modeling capability of the Argonne Reactor Computation (ARC) suite and the sensitivity analysis and optimization modules embedded in the DAKOTA code, both have been integrated within the NEAMS Workbench. The genetic algorithm is selected as the optimization method for its robustness, while the option of surrogate modeling is also explored to alleviate the computational burden caused by employing the ARC direct core physics simulation and thus enhance the efficiency of the optimization. Finally, the normalized deviations of performance parameters of the near-optimal solution from those of the ideal core are calculated and used as criteria to down select the final core design. The developed framework is applied to the Advanced Burner Test Reactor (ABTR) core, and optimal solutions are determined by balancing various objectives simultaneously, including peak fast flux, core volume, power, reactivity swing, plutonium mass feed, while at the same time satisfying the predefined constraints due to safety or economics considerations. The optimal ABTR core design obtained using the direct physical simulation and surrogated model are compared and discussed. It is found that using the accurately constructed surrogate models can significantly reduce the required computational time while maintaining satisfactory accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Mid-century climate change impacts on tornado-producing tropical cyclones

Tornadoes are a co-occurring extreme that can be produced by landfalling tropical cyclones (TCs). These tornadoes can exacerbate the loss of life and property damage caused by the TC from which they were spawned. It is uncertain how the severe weather environments of landfalling TCs may change in a future climate and how this could impact tornado activity from TCs. In this study, we investigated four TCs that made landfall in the U.S. and produced large tornado outbreaks. We performed four-member ensembles of convective-allowing (4-km resolution) regional climate model simulations representing each TC in the historical climate and a mid-twenty-first century future climate. To identify potentially tornadic storms, or TC-tornado (TCT) surrogates, we used thresholds for three-hourly maximum updraft helicity and radar reflectivity, as tornadoes are not resolved in the model. We found that the ensemble-mean number of TCT-surrogates increased substantially (56–299%) in the future, supported by increases in most-unstable convective available potential energy, surface-to-700-hPa bulk wind shear, and 0–1-km storm-relative helicity in the tornado-producing region of the TCs. On the other hand, future changes in most-unstable convective inhibition had minimal influence on future TCT-surrogates. This provides robust evidence that tornado activity from TCs may increase in the future. Furthermore, TCT-surrogate frequency between 00Z and 09Z increased for three of the four cases, suggesting enhanced tornado activity at night, when people are asleep and more likely to miss warnings. All of these factors indicate that TC-tornadoes may become more frequent and a greater hazard in the future, compounding impacts from future increases in TC winds and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Data-driven model for divertor plasma detachment prediction

We present a fast and accurate data-driven surrogate model for divertor plasma detachment prediction leveraging the latent feature space concept in machine learning research. Our approach involves constructing and training two neural networks: an autoencoder that finds a proper latent space representation (LSR) of plasma state by compressing the multi-modal diagnostic measurements and a forward model using multi-layer perception (MLP) that projects a set of plasma control parameters to its corresponding LSR. By combining the forward model and the decoder network from autoencoder, this new data-driven surrogate model is able to predict a consistent set of diagnostic measurements based on a few plasma control parameters. In order to ensure that the crucial detachment physics is correctly captured, highly efficient 1D UEDGE model is used to generate training and validation data in this study. The benchmark between the data-driven surrogate model and UEDGE simulations shows that our surrogate model is capable of providing accurate detachment prediction (usually within a few per cent relative error margin) but with at least four orders of magnitude speed-up, indicating that performance-wise, it has the potential to facilitate integrated tokamak design and plasma control. Comparing with the widely used two-point model and/or two-point model formatting, the new data-driven model features additional detachment front prediction and can be easily extended to incorporate richer physics. This study demonstrates that the complicated divertor and scrape-off-layer plasma state has a low-dimensional representation in latent space. Understanding plasma dynamics in latent space and utilising this knowledge could open a new path for plasma control in magnetic fusion energy research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Impact of Selected High-Performance Fuel Blends on Three-Way Catalyst Light Off under Synthetic Spark-Ignition Engine-Exhaust Conditions

The U.S. Department of Energy funded Co-Optimization of Fuels and Engines initiative aims to simultaneously develop advanced engines along with high-performance fuels to reduce petroleum consumption. The engine exhaust of spark-ignited light-duty vehicles contains pollutants such as nitrogen oxides, carbon monoxide (CO), and non-methane organic gases. When operated above their “light-off” temperature, three-way catalysts (TWCs) efficiently control the emissions of these pollutants from the vehicle exhaust. However, below the catalyst light-off temperature, during cold start, the TWCs are not effective. Thus, the stringent environmental regulations necessitate cold-start compliance of advanced engines operating on novel fuels for commercialization. Exhaust composition strongly impacts the effectiveness of TWCs. Hence, ensuring that the high-performance fuels under consideration do not have detrimental effects on current emissions control technology is necessary. To mitigate cold-start emissions, a low light-off temperature of the fuel on the TWC is desirable. As real-world fuels are multicomponent blends, we conducted investigations into the light-off behavior of representative fuel mixtures on a three-mode redox-aged commercial TWC under synthetic engine-exhaust conditions. The high-performance fuels in this study included 10–30% volumetric levels of ethanol, isobutanol, diisobutylene, and an aromatic mixture. Each of these high-performance fuel components was mixed into a gasoline surrogate blendstock for oxygenate blending (BOB). Our results showed that aromatics and alkenes in the surrogate BOB inhibit low-temperature reactivity of alkanes, alcohols, and CO on the TWC and dominate the blend light-off behavior. All the high-performance fuel blends had a very similar light-off behavior to the surrogate gasoline BOB, indicating that blending up to 30% (by vol.) of high-performance blendstocks in a gasoline base fuel can potentially reduce greenhouse gas emissions through improved engine efficiency and petroleum displacement without jeopardizing the ability to meet emissions regulations. While some high-performance blendstocks demonstrated lower light-off temperatures than a surrogate gasoline blend, taking advantage of the higher catalytic reactivity of these blendstocks to reduce cold-start emissions would require reducing the aromatic content in petroleum-based market fuels.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2017

Nanoparticle Precipitation in Irradiated and Annealed Ceria Doped with Metals for Emulation of Spent Fuels

Epsilon alloy precipitates have been observed with varied compositions and sizes in spent nuclear fuels, such as UO2. Presence of the inclusions, along with other oxide precipitates, gas bubbles and irradiation-induced structural defects, can significantly degrade the physical properties of the fuel. To predict fuel performance, a fundamental study of the precipitation processes is needed. This study uses ceria (CeO2) as a surrogate for UO2. Polycrystalline CeO2 films doped with Mo, Ru, Rh, Pd and Re (surrogate for Tc) were grown at 823 K using pulsed laser deposition, irradiated at 673 K with He+ ions, and subsequently annealed at higher temperatures. A number of methods, including transmission electron microscopy and atom probe tomography, were applied to characterize the samples. The results indicate that there is a uniform distribution of the doped metals in the as-grown CeO2 film. Pd particles of ~3 nm in size appear near dislocation edges after He+ ion irradiation to ~13 dpa. Thermal annealing at 1073 K in air leads to formation of precipitates with Mo and Pd around grain boundaries. Further annealing at 1373 K produces 70 nm sized precipitates with small grains at cavities.

Jiang, Weilin↗
DOE OSTI2019

In Situ Study of Particle Precipitation in Metal-Doped CeO2 during Thermal Treatment and Ion Irradiation for Emulation of Irradiating Fuels

Metallic particles formed in oxide fuels (e.g., UO 2 ) during neutron irradiation have an adverse impact on fuel performance. A fundamental investigation of particle precipitation is needed to predict the fuel performance and potentially improve fuel designs and operations. In this study, we report on the precipitation of Mo-dominant β-phase particles in polycrystalline CeO 2 (surrogate for UO 2 ) films doped with Mo, Pd, Rh, Ru, and Re (surrogate for Tc). In situ heating scanning transmission electron microscopy indicates that particle precipitation starts at ~1073 K with a limited particle growth to ~10 nm. While particle concentration increases with increasing temperature, particle size remains largely unchanged up to 1273 K. There is a dramatic change in the microstructure following vacuum annealing at 1373 K, probably due to phase transition of reduced cerium oxide. At the high temperature, particles grow up to 75 nm or larger with distinctive facets. The particles are predominantly composed of Mo with a body-centered cubic structure (β phase). An oxide layer was observed after storage at ambient conditions. In situ heating X-ray photoelectron spectroscopy reveals an increasing reduction of Ce charge state from 4+ to 3+ in the doped CeO 2 film at temperatures from 673 to 1273 K. In situ ion irradiation transmission electron microscopy with 2 MeV Al 2+ ions up to a dose of ~20 displacements per atom at nominally room temperature does not lead to precipitation of visible particles. However, irradiation with 1.7 MeV Au 3+ ions to ~10 dpa at 973 K produces ~2 nm sized pure Pd particles; Au 3+ irradiation at 1173 K appears to result in precipitates of ~6 nm in size. Some of the defects produced by ion irradiation could be nucleation sites for precipitation, leading to generation of smaller particles with a higher concentration.

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

Deep Learning for Subsurface Flow: A Comparative Study of U‐Net, Fourier Neural Operators, and Transformers in Underground Hydrogen Storage

Subsurface flow research is essential for the sustainable management of natural resources and the environment. Deep learning (DL) has significantly advanced this field by developing efficient and accurate surrogate models to replace computationally expensive physics‐based simulations. These surrogate models are commonly used to predict the spatiotemporal evolution of state variables, such as gas saturation and reservoir pressure, in heterogeneous geological formations. Despite the various DL models applied to this task, there is a lack of studies systematically comparing their performance. This absence of comparative analysis leads to somewhat arbitrary DL model selection in subsurface flow research, resulting in suboptimal performance and potentially inaccurate predictions. To bridge this gap, we conduct a systematic comparison study of three popular DL architectures—U‐Net, Fourier Neural Operators (FNO), and Segmentation Transformer (SETR)—in surrogate modeling of underground hydrogen storage (UHS). We focus on UHS due to its promise of enhancing clean energy resilience and its cyclic operational conditions that represent common scenarios in various subsurface applications. We evaluate the models based on accuracy, training cost, and inference speed. The comparison shows that U‐Net achieves the highest accuracy, followed by SETR and FNO. Despite its lower accuracy, FNO has the highest inference speed. SETR offers competitive accuracy with the least training memory usage, demonstrating the potential of transformers in learning subsurface flow. Our results provide guidance for selecting DL models for surrogate modeling in a wide range of subsurface flow problems.

42 ENGINEERING↗
DOE OSTI2024

An active learning framework for the rapid assessment of galvanic corrosion

The current present in a galvanic couple can define its resistance or susceptibility to corrosion. However, as the current is dependent upon environmental, material, and geometrical parameters it is experimentally costly to measure. To reduce these costs, Finite Element (FE) simulations can be used to assess the cathodic current but also require experimental inputs to define boundary conditions. Due to these challenges, it is crucial to accelerate predictions and accurately predict the current output for different environments and geometries representative of in-service conditions. Machine learned surrogate models provides a means to accelerate corrosion predictions. However, a one-time cost is incurred in procuring the simulation and experimental dataset necessary to calibrate the surrogate model. Therefore, an active learning protocol is developed through calibration of a low-cost surrogate model for the cathodic current of an exemplar galvanic couple (AA7075-SS304) as a function of environmental and geometric parameters. The surrogate model is calibrated on a dataset of FE simulations, and calculates an acquisition function that identifies specific additional inputs with the maximum potential to improve the current predictions. This is accomplished through a staggered workflow that not only improves and refines prediction, but identifies the points at which the most information is gained, thus enabling expansion to a larger parameter space. The protocols developed and demonstrated in this work provide a powerful tool for screening various forms of corrosion under in-service conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Data-driven emulation of modal aerosol microphysics via neural operator-based modeling

The complexity and the small characteristic scales of aerosol microphysical processes pose a big challenge for accurate and efficient Earth system simulations at regional and global scales. In this work, we construct and evaluate a surrogate model: the aerosol deep operator network (ADON), a physics-inspired dual-net architecture for emulating the aerosol microphysics parameterization suite in the version 2 of the Energy Earth System Model (E3SMv2). The current version of the surrogate model is trained on a dataset comprising 9.8 million samples obtained from a global E3SMv2 simulation with the horizontal resolution of about one degree under cloud-free conditions. Incorporating domain spatial and temporal coordinates, as well as principle components extracted from training data, the dual-net surrogate model effectively captures the intricate representations of aerosol and the relationship with atmospheric state variables, achieving an R-squared score over $$95.7\%$$ for all the lognormal aerosol modes in the extrapolated regime. The validated model provides feature importance of input variables and their impact on the predictive capacity of the surrogate model in relation to the E3SM. The computational cost of online inference time deployed on CPUs and GPUs with lower precisions highlights ADON’s efficiency and potential in robust predictive modeling for large-scale Earth system computations.

Bai, Zhe↗
DOE OSTI2015

In-flight observations of low-mode ρR asymmetries in NIF implosions

Charged-particle spectroscopy is used to assess implosion symmetry in ignition-scale indirect-drive implosions for the first time. Surrogate D{sup 3}He gas-filled implosions at the National Ignition Facility produce energetic protons via D+{sup 3}He fusion that are used to measure the implosion areal density (ρR) at the shock-bang time. By using protons produced several hundred ps before the main compression bang, the implosion is diagnosed in-flight at a convergence ratio of 3–5 just prior to peak velocity. This isolates acceleration-phase asymmetry growth. For many surrogate implosions, proton spectrometers placed at the north pole and equator reveal significant asymmetries with amplitudes routinely ≳10%, which are interpreted as ℓ=2 Legendre modes. With significant expected growth by stagnation, it is likely that these asymmetries would degrade the final implosion performance. X-ray self-emission images at stagnation show asymmetries that are positively correlated with the observed in-flight asymmetries and comparable in magnitude, contradicting growth models; this suggests that the hot-spot shape does not reflect the stagnated shell shape or that significant residual kinetic energy exists at stagnation. More prolate implosions are observed when the laser drive is sustained (“no-coast”), implying a significant time-dependent asymmetry in peak drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

3D printer generated thorax phantom with mobile tumor for radiation dosimetry

This article describes the design, construction, and properties of an anthropomorphic thorax phantom with a moving surrogate tumor. This novel phantom permits detection of dose both inside and outside a moving tumor and within the substitute lung tissue material. A 3D printer generated the thorax shell composed of a chest wall, spinal column, and posterior regions of the phantom. Images of a computed tomography scan of the thorax from a patient with lung cancer provided the template for the 3D printing. The plastic phantom is segmented into two materials representing the muscle and bones, and its geometry closely matches a patient. A surrogate spherical plastic tumor controlled by a 3D linear stage simulates a lung tumor’s trajectory during normal breathing. Sawdust emulates the lung tissue in terms of average and distribution in Hounsfield numbers. The sawdust also provides a forgiving medium that permits tumor motion and sandwiching of radiochromic film inside the mobile surrogate plastic tumor for dosimetry. A custom cork casing shields the film and tumor and eliminates film bending during extended scans. The phantom, lung tissue surrogate, and radiochromic film are exposed to a seven field plan based on an ECLIPSE plan for 6 MV photons from a Trilogy machine delivering 230 cGy to the isocenter. The dose collected in a sagittal plane is compared to the calculated plan. Gamma analysis finds 8.8% and 5.5% gamma failure rates for measurements of large amplitude trajectory and static measurements relative to the large amplitude plan, respectively. These particular gamma analysis results were achieved using parameters of 3% dose and 3 mm, for regions receiving doses >150 cGy. The plan assumes a stationary detection grid unlike the moving radiochromic film and tissues. This difference was experimentally observed and motivated calculated dose distributions that incorporated the phase of the tumor periodic motion. These calculations modestly improve agreement between the measured and intended doses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Examining astrophysical gas cloud collapse using an optical depth-scaled, x-ray-irradiated, carbon-foam sphere

When stellar radiation interacts with a molecular cloud, the cloud's fate depends on the strength of the incident radiation and the radiation's mean-free-path within the cloud [F. Bertoldi, Astrophys. J. 346, 735–755 (1989)]. Under the right conditions, the radiation compresses the cloud and a star formation may occur. Where and when the stellar formation occurs in the cloud's collapse are open questions. Direct observation of the complete star–cloud lifecycle is nearly impossible due to the immense timescales and distances over which the interaction occurs. Laboratory astrophysics offers a way to investigate such a system by scaling the important astrophysical parameters to the laboratory. This work describes laboratory experiments to study the radiation-driven implosion of clouds, using x rays from a laser-irradiated, thin, gold foil as a surrogate star and a carbon-foam sphere as a surrogate cloud. An optically thick system, theoretically corresponding to a star-forming regime, was selected by choice of the foam density. Gold foil and sphere motions were imaged by x-ray radiography. Radiographic images show the formation of an interface between rarefied gold and carbon plasmas, a shock moving into the sphere, and a blunting of the initial sphere's shape. Measurements show that the shock moved linearly around 64 μm/ns into the sphere, and the gold–carbon interface formed by 2 ns at the sphere edge remained stationary. The deformation of the sphere was driven by the incident radiation and not by mechanical pressures applied by gold plasma. The blunting of the sphere was likely due to the geometric reduction of flux near the sphere's poles. Higher x-ray flux near the sphere's equator caused high compression and a faster shock, which flattened the sphere. We will discuss the results and implications of our observations.

VanDervort, R. W. [University of Michigan 1 , Ann ↗
DOE OSTI2022

Neural network acceleration of large-scale structure theory calculations

Here, we make use of neural networks to accelerate the calculation of power spectra required for the analysis of galaxy clustering and weak gravitational lensing data. For modern perturbation theory codes, evaluation time for a single cosmology and redshift can take on the order of two seconds. In combination with the comparable time required to compute linear predictions using a Boltzmann solver, these calculations are the bottleneck for many contemporary large-scale structure analyses. Here, we construct neural network-based surrogate models for Lagrangian perturbation theory (LPT) predictions of matter power spectra, real and redshift space galaxy power spectra, and galaxy-matter cross power spectra that attain ~ 0.1% (at one sigma) accuracy over a broad range of scales in a ωCDM parameter space. The neural network surrogates can be evaluated in approximately one millisecond, a factor of 1000 times faster than the full Boltzmann code and LPT computations. In a simulated full-shape redshift space galaxy power spectrum analysis, we demonstrate that the posteriors obtained using our surrogates are accurate compared to those obtained using the full LPT model. We make our surrogate models public at https://github.com/sfschen/EmulateLSS, so that others may take advantage of the speed gains they provide to enable rapid iteration on analysis settings, something that is essential in complex contemporary large-scale structure analyses.

79 ASTRONOMY AND ASTROPHYSICS↗