SEARCH · Search NASA

Results for “variances”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Constructing auxiliary dynamics for nonequilibrium stationary states by variance minimization

We present a strategy to construct guiding distribution functions (GDFs) based on variance minimization. Auxiliary dynamics via GDFs mitigates the exponential growth of variance as a function of bias in Monte Carlo estimators of large deviation functions. The variance minimization technique exploits the exact properties of eigenstates of the tilted operator that defines the biased dynamics in the nonequilibrium system. We demonstrate our techniques in two classes of problems. In the continuum, we show that GDFs can be optimized to study the interacting driven diffusive systems where the efficiency is systematically improved by incorporating higher correlations into the GDF. On the lattice, we use a correlator product state ansatz to study the 1D weakly asymmetric simple exclusion process. We show that with modest resources, we can capture the features of the susceptibility in large systems that mark the phase transition from uniform transport to a traveling wave state. Our work extends the repertoire of tools available to study nonequilibrium properties in realistic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

A case study: Comparing testers for means and variances in a balanced incomplete block design

In this article, we present a gauge R & R study that is different from the standard setup in several ways. A balanced incomplete block design (BIBD) was used to collect the data because the measurement is destructive. Besides addressing reproducibility, differences in repeatability (i.e., tester variances) were also of interest. Here we use Bayesian methods to analyze the study data to address whether the three tester means and variances are similar. We also consider what happens if the BIBD structure of the data and possibly different tester variances are ignored in analyzing the study data.

42 ENGINEERING↗

DOE OSTI2025

Suppressing the sample variance of DESI-like galaxy clustering with fast simulations

Ongoing and upcoming galaxy redshift surveys, such as the Dark Energy Spectroscopic Instrument (DESI) survey, will observe vast regions of sky and a wide range of redshifts. In order to model the observations and address various systematic uncertainties, N-body simulations are routinely adopted, however, the number of large simulations with sufficiently high mass resolution is usually limited by available computing time. Therefore, achieving a simulation volume with the effective statistical errors significantly smaller than those of the observations becomes prohibitively expensive. In this study, we apply the Convergence Acceleration by Regression and Pooling (CARPool) method to mitigate the sample variance of the DESI-like galaxy clustering in the AbacusSummit simulations, with the assistance of the quasi-N-body simulations FastPM. Based on the halo occupation distribution (HOD) models, we construct different FastPM galaxy catalogs, including the luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars, with their number densities and two-point clustering statistics well matched to those of AbacusSummit. We also employ the same initial conditions between AbacusSummit and FastPM to achieve high cross-correlation, as it is useful in effectively suppressing the variance. Our method of reducing noise in clustering is equivalent to performing a simulation with volume larger by a factor of 5 and 4 for LRGs and ELGs, respectively. We also mitigate the standard deviation of the LRG bispectrum with the triangular configurations k 2 = 2k 1 = 0.2 h Mpc -1 by a factor of 1.6. With smaller sample variance on galaxy clustering, we are able to constrain the baryon acoustic oscillations (BAO) scale parameters to higher precision. The CARPool method will be beneficial to better constrain the theoretical systematics of BAO, redshift space distortions (RSD) and primordial non-Gaussianity (NG).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Propagating sample variance uncertainties in redshift calibration: simulations, theory, and application to the COSMOS2015 data

ABSTRACT Cosmological analyses of galaxy surveys rely on knowledge of the redshift distribution of their galaxy sample. This is usually derived from a spectroscopic and/or many-band photometric calibrator survey of a small patch of sky. The uncertainties in the redshift distribution of the calibrator sample include a contribution from shot noise, or Poisson sampling errors, but, given the small volume they probe, they are dominated by sample variance introduced by large-scale structures. Redshift uncertainties have been shown to constitute one of the leading contributions to systematic uncertainties in cosmological inferences from weak lensing and galaxy clustering, and hence they must be propagated through the analyses. In this work, we study the effects of sample variance on small-area redshift surveys, from theory to simulations to the COSMOS2015 data set. We present a three-step Dirichlet method of resampling a given survey-based redshift calibration distribution to enable the propagation of both shot noise and sample variance uncertainties. The method can accommodate different levels of prior confidence on different redshift sources. This method can be applied to any calibration sample with known redshifts and phenotypes (i.e. cells in a self-organizing map, or some other way of discretizing photometric space), and provides a simple way of propagating prior redshift uncertainties into cosmological analyses. As a worked example, we apply the full scheme to the COSMOS2015 data set, for which we also present a new, principled SOM algorithm designed to handle noisy photometric data. We make available a catalogue of the resulting resamplings of the COSMOS2015 galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

A multishock model for the density variance of anisotropic, highly magnetized, supersonic turbulence

ABSTRACT Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $\sigma ^2_{\rho /\rho _0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parametrizes the star formation (SF) rate of a cloud. However, models for $\sigma ^2_{\rho /\rho _0}$ all share a common feature – the variance is assumed to be isotropic. This assumption does not hold when a trans-/sub-Alfvénic mean magnetic field, ${B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $\sigma _{\rho /\rho _0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $\sigma _{\rho /\rho _0}^2$ in the high Mach number $(\mathcal {M})$ limit as small-scale density fluctuations become suppressed by the strong ${B}_0$. The model reduces to the isotropic $\sigma _{\rho /\rho _0}^2\!-\!\mathcal {M}$ relation in the hydrodynamical limit. To validate our model, we calculate $\sigma _{\rho /\rho _0}^2$ from 12 high-resolution, three-dimensional, supersonic, sub-Alfvénic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented overdensities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvénic ${B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.

Beattie, James R.↗

DOE OSTI2023

Cosmic variance and the inhomogeneous UV luminosity function of galaxies during reionization

When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from z ≃ 6–15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same z. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by z = 6 (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, ‘rolling-over’ at M UV ≳ -17. CoDa II’s UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at z = 6–8, including the faint end, except for the faintest data point at z = 6, based on one galaxy at M UV = -12.5. According to CoDa II, the probability of observing the latter is $\sim 5~{{\ \rm per\ cent}}$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. Here, we find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2–4, primarily by underestimating the variance in halo abundance when the sample volume is small.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Variance reduction in lattice QCD observables via normalizing flows

Normalizing flows can be used to construct unbiased, reduced-variance estimators for lattice field theory observables that are defined by a derivative with respect to action parameters. This work implements the approach for observables involving gluonic operator insertions in the SU(3) Yang-Mills theory and two-flavor QCD in four space-time dimensions. Variance reduction by factors of 10–60 is achieved in glueball correlation functions and in gluonic matrix elements related to hadron structure, with demonstrated computational advantages. The observed variance reduction is found to be approximately independent of the lattice volume, so volume transfer can be utilized to minimize training costs.

Abbott, Ryan [Columbia U.; MIT, Cambridge, CTP; IA↗

DOE OSTI2025

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Two-Stage Estimation and Variance Modeling for Latency-Constrained Variational Quantum Algorithms

The quantum approximate optimization algorithm (QAOA) has enjoyed increasing attention in noisy, intermediate-scale quantum computing with its application to combinatorial optimization problems. QAOA has the potential to demonstrate a quantum advantage for NP-hard combinatorial optimization problems. As a hybrid quantum-classical algorithm, the classical component of QAOA resembles a simulation optimization problem in which the simulation outcomes are attainable only through a quantum computer. The simulation that derives from QAOA exhibits two unique features that can have a substantial impact on the optimization process: (i) the variance of the stochastic objective values typically decreases in proportion to the optimality gap, and (ii) querying samples from a quantum computer introduces an additional latency overhead. In this paper, we introduce a novel stochastic trust-region method derived from a derivative-free, adaptive sampling trust-region optimization method intended to efficiently solve the classical optimization problem in QAOA by explicitly taking into account the two mentioned characteristics. The key idea behind the proposed algorithm involves constructing two separate local models in each iteration: a model of the objective function and a model of the variance of the objective function. Exploiting the variance model allows us to restrict the number of communications with the quantum computer and also helps navigate the nonconvex objective landscapes typical in QAOA optimization problems. In conclusion, we numerically demonstrate the superiority of our proposed algorithm using the SimOpt library and Qiskit when we consider a metric of computational burden that explicitly accounts for communication costs.

Derivative-free Optimization↗

DOE OSTI2023

Cost-optimized Automated Variance Reduction for Highly Angle-dependent Radiation Transport Analyses [Dissertation]

Monte Carlo variance-reduction techniques that directly bias particle direction have historically received limited attention despite being useful for many radiation transport applications such as those in which radiation travels through large regions with a low probability of colliding. One such technique is known as DXTRAN (short for deterministic transport) in the MCNP Monte Carlo radiation transport code. Until now, effectively applying DXTRAN in calculations is largely based on empirical observations of computational performance. Optimal DXTRAN parameters are identified through manual iteration. This work develops new mathematical descriptions of the DXTRAN variance-reduction process and demonstrates a new automated variance-reduction method that applies these mathematical formulations to determine the optimal application of DXTRAN in a given problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

LLNS Request for Variance from 10 CFR 851 for Construction of LLNL Digital Infrastructure Capabilities Expansion Project

Lawrence Livermore National Security, LLC (LLNS) requests variance to 10 CFR 851 (the Rule) for the construction of the Digital Infrastructure Capabilities Expansion (DICE) project. The request isin continued support of the National Nuclear Security Administration (NNSA) “pilot” to streamline the delivery of commercial-like line-item construction projects under $50 million. This pilot is aligned with the approved NA-50 safety approach. Experience gained in the construction of the new Lawrence Livermore National Laboratory (LLNL) Emergency Operations Facility (EOC) under the variance granted in June 2020 is very encouraging. While detailed analysis will be forthcoming after that project is completed, there are strong indicators that a benefit has been realized by the subcontractors executing the work. The variance has permitted local/regional construction businesses contracted to Management & Operation contractor operated facilities and lower tier subcontractors to bid, operate, execute, and comply with safety and health standards they utilize in their business outside of the Department of Energy (DOE) Complex. While compliance to standards that are as safe and healthful as those in 10 CFR 851: confusion, additional cost, and inefficiencies associated with the burden of compliance to an unfamiliar regulatory spectrum has been avoided. Safety and health performance on the EOC construction project including the number of inspection and surveillance issues is comparable to other subcontractors and projects at LLNL executed under 10 CFR 851.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2023

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗

DOE OSTI2026

Allan Variance is Bootstrap Aggregation for Spectral Estimation

Characterization of clocks and inertial sensors, such as accelerometers and gyroscopes, typically includes Allan variance analysis. Allan variance is ubiquitous in timing and navigation communities which may appear niche compared with generalized spectral analysis. This note provides some motivation for Allan Variance for audiences more familiar with spectral analysis.

Walker, Michael Ray [Sandia National Laboratories ↗

DOE OSTI2020

Online Control of Process Variance Using Feedback

For many stochastic control problem formulations, design objectives typically include constraining the covariance of the outputs of the system. Recent engineering challenges have demonstrated an interest in controlling system output covariances to a setpoint. However, open-loop control of system output covariance can be exceptionally challenging when the system has significant plant uncertainties. This paper explores closed-loop control of the system’s output covariance using a windowed variance function. The stochastic nature of the windowed variance feedback means that root locus and other classical control design techniques cannot be immediately used. Expected value analysis is used to generate transfer functions for the average response of such systems, and the results are used to design controllers. Stochastic simulations are presented to verify the analytical results and make other observations concerning closed-loop variance control.

Bieniek, Mitchell↗

DOE OSTI2022

A method for transporting cloud-resolving model variance in a multiscale modeling framework

Abstract. An unphysical checkerboard pattern has recently been identified in the multiscale modeling framework configuration of the Energy Exascale Earth System Model (E3SM-MMF) that is hypothesized to be associated with the inability of large-scale dynamics to transport fluctuations within the embedded cloud-resolving model (CRM) on the global grid. To address this issue, a method is presented to facilitate the large-scale transport of CRM variance in E3SM-MMF. Simulation results show that the method is effective at reducing the occurrence of unphysical checkerboard patterns on a range of timescales from days to years. This result is confirmed both subjectively through visual inspection and quantitatively with a previously developed pattern categorization technique. The CRM variance transport does not significantly alter the model climate, although it does tend to reduce temporal variance on fields associated with convection on the global grid.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Global variance decomposition of downscaled and bias-corrected CMIP6 climate projections

This dataset provides the results of a global variance decomposition of downscaled and bias-corrected CMIP6 climate projections. The total projection variance for a set of climate metrics is partitioned into contributions from: scenario uncertainty, model/GCM uncertainty, downscaling and bias-correction uncertainty, and interannual variability. The contribution from each source is expressed as a percentage of the total variance. Seven climate metrics are analyzed: Annual average temperature (avg_tas.nc) Annual total precipitation (tot_pr.nc) Annual maximum of daily maximum temperature (max_tasmax.nc) Annual maximum 1-day precipitation (max_pr.nc) Annual number of extremely hot days (hot_days.nc) Annual number of extremely wet days (wet_days.nc) Annual number of dry days (dry_days.nc) Extremely hot/wet days are defined to occur when temperature/precipitation exceeds the local 99th percentile defined over 1980-2014. Dry days are defined to occur when daily precipitation is less than 1mm. all_metrics_timesliced.nc gives the results for all metrics averaged over three 20-year periods: 2020-2039, 2050-2069, 2080-2099. For more details on the methods, see: Lafferty & Sriver, Downscaling and bias-correction contribute considerable uncertainty to local climate projections in CMIP6, npj Climate & Atmospheric Science (2023) An interactive visualization of this data can be found at: https://lafferty-sriver-2023-downscaling-uncertainty.msdlive.org

Lafferty, David↗

DOE OSTI2020

The Effect Of Ionization Variance On Nuclear-recoil Dark Matter Searches

Direct detection dark matter experiments are increasingly interested in the low-mass dark matter parameter space, but zero-background low-mass searches require event separation between the electron and nuclear recoil bands, which requires a proper understanding of detector energy reconstruction. Previous simulations have shown that we do not entirely understand the ionization efficiency (yield) for electron and nuclear recoils, as the assumption that the distribution for the yield is normally distributed for a true recoil energy is violated. Since the yield distribution may directly affect dark matter low-mass limits, it is crucial we understand how the yield is distributed. A component to understanding the yield distribution is the variance in the number of electron-hole pairs produced or ionization variance. This effect has been studied relatively infrequently as experiments have been interested in large energy deposits (10 - 100 keV) and could accurately separate electron and nuclear recoil events. For electron recoils, the ionization variance is described by a “Fano” factor. For nuclear recoils the effect can be parameterized by an “effective” Fano factor, which has similar definition but a different physical origin. The nuclear recoil “effective” Fano factor is shown to be much larger than the electron-recoil version above around 10 keV deposited energy.

Matheny, Mitchell Douglas↗

DOE OSTI2024

Sequential Selection for Minimizing the Variance with Application to Crystallization Experiments

For many crystal-based products (e.g., pharmaceuticals, energy storage), the size uniformity is not only a key quality attribute, but sometimes also an indicator of other attributes such as solid purity. This article proposes a sequential selection approach to find a proper experimental setting that leads to high uniformity, or equivalently, small variance for crystal sizes, from the advanced slug flow reaction crystallization process of a model crystal, called manganese oxalate hydrate. The proposed sequential selection approach contains a Bayesian adaptive method to incorporate new uniformity measurements in each step and two design acquisition functions to improve the selection of the most promising experimental setting in terms of minimizing the variance. We study the performance of the proposed approach through multiple synthetic numerical studies, as well as a case study based on data from slug flow crystallization experiments. Throughout these studies, the proposed approach shows competitive performance in identifying the best experimental setting.

Expected improvement↗