Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical 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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An Algorithm to Calibrate and Correct the Response to Unpolarized Radiation of the X-Ray Polarimeter Onboard IXPE

The Gas Pixel Detector (GPD) is an X-ray polarimeter to fly onboard IXPE and other missions. To correctly measure the source polarization, the response of IXPE’s GPDs to unpolarized radiation has to be calibrated and corrected. In this paper, we describe the way such response is measured with laboratory sources and the algorithm to apply such correction to the observations of celestial sources. The latter allows to correct the response to polarization of single photons, therefore allowing great flexibility in all the subsequent analysis. Our correction approach is tested against both monochromatic and nonmonochromatic laboratory sources and with simulations, finding that it correctly retrieves the polarization up to the statistical limits of the planned IXPE observations.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2022

Measuring the Speed of Sound Using Cumulants of Baryon Number

Here, we show that the values of the first three cumulants of the baryon number distribution can be used to calculate the isothermal speed of sound and its logarithmic derivative with respect to the baryon-number density. We discuss applications of this result to heavy-ion collision experiments and address possible challenges, including effects due to baryon-number conservation, differences between proton and baryon cumulants, and the influence of finite-number statistics on fluctuation observables in both ex periment and hadronic transport simulations. In particular, we investigate the relation between quantities calculated in infinite, continuous matter, and observables obtained in simulations using a finite number of particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

RTN-056: Study of the Photon Transfer Curve in the CCD detectors of the Vera C. Rubin Observatory

The RECA internship program provides Colombian students with an opportunity to enhance their research skills in Astronomy, Astrophysics, and Cosmology. During this three-month program, our main objective was to study the Photon Transfer Curves (PTC) of the Vera C. Rubin Observatory, specifically the gain, and to compare it with the gain obtained through pairs of flats. Overall, the study of PTCs is crucial in understanding the performance of detectors and instruments used in Astronomy. The Vera C. Rubin Observatory is an important facility that will enable researchers to carry out a wide range of studies in this field, making it essential to investigate its gain performance. We used run 13144 to construct the PTCs and 13186 to analyze the crosstalk. We employed the LSST Science Pipelines (also known as the DM stack), a software under development for this observatory, which performs all the necessary reductions for the construction of the PTCs. We also used simulations to replicate the observed effects. Initially, we found a 5% difference between the gain calculated by PTC and pairs of flats for a flow range between 5000 and 10000 ADU. Simulations showed that this difference was due to the handling of statistics and the assumption that the distribution following the Lupton equation is Gaussian. We found an error interval for this flow region based on the vendor, with (1.8 ± 0.7, 4.1 ± 0.9) % for E2V and (0.85 ± 0.7, 2.2 ± 0.9) % for ITL. From the PTC, we also obtained the average Full Well Capacity of LSSTCam as 130000 ± 10000$ electrons. We identified a list of segments where we found differences with the results obtained by SLAC National Acceleration Laboratory in PTC parameters, low saturation level, or other defects. We detected and corrected the effect of statistics in the gain calculation using pairs of flats and proposed a code change, which was implemented in the pipeline software. We do not recommend correcting for crosstalk as it does not significantly affect the parameters and does not change the shape of the PTC. However, the opposite is true for the nonlinearity correction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Methane Emissions Produced by Pneumatic Devices and Produced Water Tanks on Natural Gas Wellsites

In recent years, more scrutiny has been placed on the release of greenhouse gasses (GHGs) by the oil and gas industry in North America as climate concerns increase. In the coming years, stricter regulations proposed by the Environmental Protection Agency (EPA) concerning the release of GHGs will further direct attention to GHG emissions quantification, modeling, and overall reduction. Active natural gas wellsites employ production equipment that routinely emits methane during normal operation. Production equipment often includes, but is not limited to, gas production units (GPUs), pneumatic controllers (PCs) and actuators powered by produced natural gas, and produced water storage tanks. This work details direct methane emission measurements taken at active natural gas wellsites originating from intermittent pneumatic controllers and produced water storage tanks. Efforts to model methane flows originating from these sources are also detailed, and those results are compared to direct measurement. The direct methane emissions measurement campaign collected flow data from two active dry wellsites in the Marcellus shale region. These measurement campaigns captured, in total, several weeks of continuous data from these sites. Data included multiple pneumatic device vent flowrate channels, a single channel monitoring total emissions vented from the produced water storage tanks, and corresponding site weather data. Periodic composition measurements were taken from flows vented from produced water tanks. Computational modeling efforts attempted to mimic a given site’s methane emissions behavior given limited input parameters. Direct measurement data and simulation data yielded, through reduction and analysis, apparent methane emissions factors, statistics on pneumatic device actuation events, important relationships between GPU parameters and pneumatic device actuation events, produced water storage tank vent flows and flow compositions, and parameters most contributing to the magnitude and frequency of these fugitive methane emissions. Measurements, considering data from a 72hr period, produced PC emissions factors from 0.008 to 0.03 SCFH per device with an average device emitting 0.026 SCFH at the first site. The second site produced PC emissions factors spanning from 0.05 to 25.47 SCFH with an average device emitting 1.93 SCFH not including the 25.47 SCFH device, which was assumed to be malfunctioning. Two unique water tank emission measurements were made, 618 g CH4 per hour and 32.6 g CH4 per hour. On a water basis, the water tank emissions factors were 0.66 and 0.042 kg CH4 per bbl of water, respectively.

03 NATURAL GAS↗
DOE OSTI2021

Simulated Changes in Tropical Cyclone Size, Accumulated Cyclone Energy and Power Dissipation Index in a Warmer Climate

Detection, attribution and projection of changes in tropical cyclone intensity statistics are made difficult from the potentially decreasing overall storm frequency combined with increases in the peak winds of the most intense storms as the climate warms. Multi-decadal simulations of stabilized climate scenarios from a high-resolution tropical cyclone permitting atmospheric general circulation model are used to examine simulated global changes from warmer temperatures, if any, in estimates of tropical cyclone size, accumulated cyclonic energy and power dissipation index. Changes in these metrics are found to be complicated functions of storm categorization and global averages of them are unlikely to easily reveal the impact of climate change on future tropical cyclone intensity statistics.

Wehner, Michael (ORCID:0000000159910082)↗
DOE OSTI2022

Deep Learning for Subgrid‐Scale Turbulence Modeling in Large‐Eddy Simulations of the Convective Atmospheric Boundary Layer

Abstract In large‐eddy simulations, subgrid‐scale (SGS) processes are parameterized as a function of filtered grid‐scale variables. First‐order, algebraic SGS models are based on the eddy‐viscosity assumption, which does not always hold for turbulence. Here we apply supervised deep neural networks (DNNs) to learn SGS stresses from a set of neighboring coarse‐grained velocity from direct numerical simulations of the convective boundary layer at friction Reynolds numbers Re τ up to 1243 without invoking the eddy‐viscosity assumption. The DNN model was found to produce higher correlation between SGS stresses compared to the Smagorinsky model and the Smagorinsky‐Bardina mixed model in the surface and mixed layers and can be applied to different grid resolutions and various stability conditions ranging from near neutral to very unstable. The DNN model can capture key statistics of turbulence in a posteriori (online) tests when applied to large‐eddy simulations of the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Map-level baryonification: unified treatment of weak lensing two-point and higher-order statistics

Precision cosmology benefits from extracting maximal information from cosmic structures, motivating the use of higher-order statistics (HOS) at small spatial scales. However, predicting how baryonic processes modify matter statistics at these scales has been challenging. The baryonic correction model (BCM) addresses this by modifying dark-matter-only simulations to mimic baryonic effects, providing a flexible, simulation-based framework for predicting both two-point and HOS. We show that a 3-parameter version of the BCM can jointly fit weak lensing maps' two-point statistics, wavelet phase harmonics coefficients, scattering coefficients, and the third and fourth moments to within 2% accuracy across all scales ℓ < 2000 and tomographic bins for a DES-Y3-like redshift distribution ( z ≲ 2), using the FLAMINGO simulations. These results demonstrate the viability of BCM-assisted, simulation-based weak lensing inference of two-point and HOS, paving the way for robust cosmological constraints that fully exploit non-Gaussian information on small spatial scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Modeling regional precipitation over the Indus River basin of Pakistan using statistical downscaling

Complex processes govern spatiotemporal distribution of precipitation within the high-mountainous headwater regions (commonly known as the upper Indus basin (UIB)), of the Indus River basin of Pakistan. Reliable precipitation simulations particularly over the UIB present a major scientific challenge due to regional complexity and inadequate observational coverage. Here, we present a statistical downscaling approach to model observed precipitation of the entire Indus basin, with a focus on UIB within available data constraints. Taking advantage of recent high altitude (HA) observatories, we perform precipitation regionalization using K-means cluster analysis to demonstrate effectiveness of low-altitude stations to provide useful precipitation inferences over more uncertain and hydrologically important HA of the UIB. We further employ generalized linear models (GLM) with gamma and Tweedie distributions to identify major dynamic and thermodynamic drivers from a reanalysis dataset within a robust cross-validation framework that explain observed spatiotemporal precipitation patterns across the Indus basin. Final statistical models demonstrate higher predictability to resolve precipitation variability over wetter southern Himalayans and different lower Indus regions, by mainly using different dynamic predictors. The modeling framework also shows an adequate performance over more complex and uncertain trans-Himalayans and the northwestern regions of the UIB, particularly during the seasons dominated by the westerly circulations. However, the cryosphere-dominated trans-Himalayan regions, which largely govern the basin hydrology, require relatively complex models that contain dynamic and thermodynamic circulations. Furthermore, we also analyzed relevant atmospheric circulations during precipitation anomalies over the UIB, to evaluate physical consistency of the statistical models, as an additional measure of reliability. Overall, our results suggest that such circulation-based statistical downscaling has the potential to improve our understanding towards distinct features of the regional-scale precipitation across the upper and lower Indus basin. Additionally, such understanding should help to assess the response of this complex, data-scarce, and climate-sensitive river basin amid future climatic changes, to serve communal and scientific interests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Comparison of quantum advantage experiments using random circuit sampling

Random circuit sampling, the task of sampling bit strings from a random unitary operator, has been implemented to demonstrate quantum advantage on the Sycamore quantum processor with 53 qubits and on the Zuchongzhi quantum processor with 56 and 61 qubits. Recently, it was claimed that classical computers using tensor network simulation could catch on to current noisy quantum processors for random circuit sampling. While the linear cross-entropy benchmark fidelity was used to certify all these claims, it may not capture statistical properties of outputs in detail. Here, we compare the bit strings sampled from classical computers using tensor network simulation by Pan et al. [F. Pan, K. Chen, and P. Zhang, Phys. Rev. Lett. 129, 090502 (2022)] and by Kalachev et al. [G. Kalachev, P. Panteleev, P. Zhou, and M.-H. Yung, arXiv:2112.15083] with the bit strings from the Sycamore quantum processor. It is shown that all of Kalachev et al.'s samples passed the NIST random number tests. The heat maps of bit strings show that Pan et al.'s and Kalachev et al.'s samples are quite different from the Sycamore or Zuchongzhi samples. The analysis with the Marchenko-Pastur distribution and the Wasssertein distances demonstrates that Kalachev et al.'s samples are statistically closer to the Sycamore samples than Pan et al.'s while the three datasets have similar values for the linear cross-entropy fidelity. In conclusion, our finding implies that further study is needed to certify or beat the claims of quantum advantage using random circuit sampling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Large-Eddy Simulation of Flow Over Boeing Gaussian Bump Using Multiagent Reinforcement Learning Wall Model: Preprint

We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 x 10^6, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation.

boundary layer↗
DOE OSTI2023

zPerf: A Statistical Gray-Box Approach to Performance Modeling and Extrapolation for Scientific Lossy Compression

With the scaling up of simulation-based scientific discovery on high-performance computing systems, the disparity between compute and I/O has increased, forcing domain scientists to save only a small amount of simulation data to persistent storage. This can result in the loss of essential physics fields that are needed for data analysis. While error-bounded lossy compression has made tremendous progress in bridging the gap between compute and I/O, the lack of understanding of compression performance remains a key hurdle to its wide adoption. Here, in this work, we present zPerf, a statistical gray-box performance modeling approach for scientific lossy compression. Our contributions are threefold: 1) We develop zPerf to estimate the performance of lossy compression techniques, based on in-depth understanding and statistical modeling for data features and core compression metrics; 2) We demonstrate the in-detailed implementation of zPerf using two case studies, where we derive the performance modeling for SZ and ZFP, two leading lossy compressors; 3) We evaluate the effectiveness of zPerf on real-world datasets across various domains. Based on the evaluation, we demonstrate the efficacy of the zPerf performance model; 4) We further discuss three case studies where zPerf is applied to extrapolate the compression ratio of SZ and ZFP with alternative encoding schemes as well as ZFP with an alternative transform scheme. Through the case studies, we demonstrate the potential of zPerf for exploring the design space of lossy compression, which has hardly been studied in the literature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Stellar Activity F-statistic for Exoplanet Surveys (SAFE)

In the search for planets orbiting distant stars, the presence of stellar activity in the atmospheres of observed stars can obscure the radial velocity signal used to detect such planets. Furthermore, this stellar activity contamination is set by the star itself and cannot simply be avoided with better instrumentation. Various stellar activity indicators have been developed that may correlate with this contamination. We introduce a new stellar activity indicator called the Stellar Activity F-statistic for Exoplanet surveys (SAFE) that has higher statistical power (i.e., probability of detecting a true stellar activity signal) than many traditional stellar activity indicators in a simulation study of an active region on a Sun-like star with a moderate-to-high signal-to-noise ratio. Also through simulation, the SAFE is demonstrated to be associated with the projected area on the visible side of the star covered by active regions. We also demonstrate that the SAFE detects statistically significant stellar activity in most of the spectra for HD 22049, a star known to have high stellar variability. Additionally, the SAFE is calculated for recent observations of the three low-variability stars HD 34411, HD 10700, and HD 3651, the latter of which is known to have a planetary companion. As expected, the SAFE for these three only occasionally detects activity. Furthermore, initial exploration appears to indicate that the SAFE may be useful for disentangling stellar activity signals from planet-induced Doppler shifts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Turbulence-Resolving Simulations of Hurricane Laura (2020): Insights Into Extreme Winds and Eyewall Turbulence

Extreme weather events pose significant risks to coastal and offshore energy infrastructure. In this work, we investigate the structure of mean winds and turbulence near the surface ( z < 300 m) that occurred during Hurricane Laura in 2020 on the US Gulf Coast. To this end, we perform turbulence-resolving simulations (..delta..x = 33.33 m) that encompass the entire eyewall of Hurricane Laura by introducing a novel method that we will refer to as Moving-Mesoscale to Static-LES Integrated Coupling (MOSAIC). The simulation results are validated against surface, boundary-layer, and mid-tropospheric observations. Our analysis of the extreme wind conditions near the surface indicates that the mean wind and turbulence profiles vary greatly at and near the eyewall, resulting in extreme values of mean wind speed (U > 50 m * s-1) and turbulence kinetic energy (k ~ 50 m2 * s-2) at altitudes above 50 m . To investigate in detail the nature of turbulence within the eyewall of the storm, we conduct higher-resolution simulations (..delta..x = 11.11 m) of a portion of the hurricane. We provide a comprehensive statistical description of turbulence in the eyewall region, highlighting the need to employ ..delta..x ~ 10 m grid spacing to capture the integral length-scales near the surface, second- and higher-order central moments (i.e., skewness and kurtosis), and spectral coherence in the flow accurately. These numerical simulations provide the most detailed description of mean winds and turbulent conditions within the eyewall of a historical tropical cyclone to date, illustrating how high-resolution simulations can enhance limited turbulence measurements in hurricanes and giving insight into the extreme conditions threatening offshore and coastal infrastructure.

17 WIND ENERGY↗
DOE OSTI2022

Perturbation theory models for LSST-era galaxy clustering: Tests with subpercent mock catalog measurements in Fourier and configuration space

Here, we analyze the clustering of galaxies using the z = 1.006 snapshot of the CosmoDC2 simulation, a high-fidelity synthetic galaxy catalog designed to validate analysis methods for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We present sub-percent measurements of the galaxy auto-correlation and galaxy-dark matter cross correlations in Fourier space and configuration space for a magnitude-limited galaxy sample. At these levels of precision, the statistical errors of the measurement are comparable to the systematic effects present in the simulation and measurement procedure; nevertheless, using a hybrid-PT model, we are able to model non-linear galaxy bias with 0.5% precision up to scales of k max = 0.5 h/Mpc and r min = 4 Mpc/h. While the linear bias parameter is measured with 0.01% precision, other bias parameters are determined with considerably weaker constraints and sometimes bimodal posterior distributions. We compare our fiducial model with lower dimensional models where the higher-order bias parameters are fixed at their co-evolution values and find that leaving these parameters free provides significant improvements in our ability to model small scale information. We also compare bias parameters for galaxy samples defined using different magnitude bands and find agreement between samples containing equal numbers of galaxies. Finally, we compare bias parameters between Fourier space and configuration space and find moderate to significant tension between the two approaches. Although our model is often unable to fit the CosmoDC2 galaxy samples within the 0.1% precision of our measurements, our results suggest that the hybrid-PT model used in this analysis is capable of modeling non-linear galaxy bias within the percent level precision needed for upcoming galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Search for 2p2h Interactions in the NOνA Near Detector

The physics of 2p2h interactions and their contribution to the NO$\nu$A near detector data are not fully understood. This study attempts to shed some light on these interactions and the accuracy of the models used to simulate them through a search for a specific 2p2h interaction in the NO$\nu$A near detector. By performing an event selection algorithm based on particle identifier algorithms run over reconstructed data, a signal region is created to minimize the background while maximizing the number of 2p2h events where a muon neutrino interacts with a neutron and a proton coupled by a meson exchange current and produces two protons and one muon. In the signal region, separation is found between the signal events and the background in plots of the angles between the protons and the muon. Although a full statistical analysis is not completed in this study, comparing the angle plots for simulation and data shows that the model used to simulate the events reasonably approximates reality and that the near detector data likely includes signal events. Signal events are also identified in event displays, further indicating that there is some contribution of the signal to the overall NO$\nu$A near detector data.

Gable, Kyle↗
DOE OSTI2021

Multi-Scenario Extreme Weather Simulator

The Multi-scenario Extreme Weather Simulator (MEWS) takes existing historical weather data and inserts trends and extreme events such as heat waves, hurricanes, extreme precipitation, or winter storms. The code includes functionality for projecting increasing frequency and duration of such events over a large range of statistically sampled realizations of the future. These weather files can then be used in building energy simulations through the Energy Plus software created by DOE. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525 SAND2021-8658 O

Villa, Daniel↗