Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variance analysis”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Landscape‐scale lidar analysis of aboveground biomass distribution in secondary Brazilian Atlantic Forest

Abstract Secondary forests account for more than half of tropical forests and represent a growing carbon sink, but rates of biomass accumulation vary by a factor of two or more even among plots in the same landscape. To better understand the drivers of this variability, we used airborne lidar to measure forest canopy height and estimate biomass over 4529 ha at Serra do Conduru Park in Southern Bahia, Brazil. We measured trees in 30 georeferenced field plots (0.25‐ha each) to estimate biomass using allometry. Then we estimated aboveground biomass density ( ABD ) across the lidar study area using a statistical model developed from our field plots. This model related the 95th percentile of the distribution of lidar return heights to ABD . We overlaid this map of ABD on a Landsat‐derived forest age map to determine rates of biomass accumulation. We found rapid initial biomass regeneration (~6 Mg/ha yr), which slowed as forests aged. We also observed high variability in both height and biomass across the landscape within forests of similar age. Nevertheless, a regression model that accounted for spatial autocorrelation and included forest age, slope, and distance to roads or open areas explained 62 and 77 percent of the landscape variation in ABD and canopy height, respectively. Thus, while there is high spatial heterogeneity in forest recovery, and the drivers of this heterogeneity warrant further investigation, we suggest that a relatively simple set of predictor variables is sufficient to explain the majority of variance in both height and ABD in this landscape.

Becknell, Justin M.↗
DOE OSTI2020

Experimental Study of the Effect of Start of Injection and Blend Ratio on Single Fuel Reformate RCCI

Abstract A new concept of single fuel reactivity-controlled compression ignition (RCCI) has been proposed through the catalytic partial oxidation (CPOX) reformation of diesel fuel. The reformed fuel mixture is then used as the low reactivity fuel and diesel itself is used as the high reactivity fuel. In this paper, two reformate mixtures from the reformation of diesel were selected for further analysis. Each reformate fuel mixture contained a significant fraction of inert gases (89% and 81%). The effects of the difference in the molar concentrations of the reformate mixtures were studied by experimenting with diesel as the direct injected fuel in RCCI over a varying start of injection timings and different blend ratios (i.e., the fraction of low and high reactivity fuels). The reformate mixture with the lower inert gas concentration had earlier combustion phasing and shorter combustion duration at any given diesel start of injection timing. The higher reactivity separation between reformate mixture and diesel, compared with gasoline and diesel, causes the combustion phasing of reformate-diesel RCCI to be more sensitive to the start of injection timing. The maximum combustion efficiency was found at a CA50 before top dead center (TDC), whereas the maximum thermal efficiency occurs at a CA50 after TDC. The range of energy-based blend ratios in which reformate-diesel RCCI is possible is between 25% and 45%, limited by ringing intensity (RI) at the low limit of blend ratios, and coefficient of variance (COV) of net indicated mean effective pressure (IMEPn) and combustion efficiency at the high limit. Intake boosting becomes necessary due to the oxygen deficiency caused by the low energy density of the reformate mixtures as it displaces intake air.

Engineering↗
DOE OSTI2026

Evaluating the limitations of Bayesian metabolic control analysis

Bayesian Metabolic Control Analysis (BMCA) is a promising framework for inferring metabolic control coefficients in data-limited scenarios, combining Bayesian inference with linear-logarithmic (lin-log) rate laws. These metabolic control coefficients quantify how changes in enzyme activities affect steady-state fluxes and metabolite concentrations across a metabolic network. However, its predictive accuracy and limitations remain underexplored. This study systematically evaluates BMCA’s ability to infer elasticity values, flux control coefficients (FCC), and concentration control coefficients (CCC) under varying data availability conditions using three synthetic metabolic network models. We demonstrate that BMCA predictions are highly dependent on the inclusion of flux and enzyme concentration data, with the omission of these datasets leading to severe inaccuracies. In our synthetic, enzyme-perturbation datasets, external metabolite concentrations had minimal impact and, in some cases, their exclusion improved predictions; when external-nutrient perturbations were introduced and those concentrations were observed, gains were at most modest. Additionally, we find that posterior estimation with both ADVI and HMC can underestimate large-magnitude elasticities in our synthetic settings, with ADVI showing somewhat higher variance under strong up-regulation; thus, recovering |elasticity| ≳ 1.5 remains challenging regardless of the inference engine. ADVI also fails to accurately infer allosteric interactions, even when regulatory effects are strong. While BMCA maintains reasonable accuracy in partially recovering the rankings of the highest FCC values, its estimates of absolute values remain constrained by prior assumptions and data limitations. Our findings reveal the BMCA algorithm’s strengths and weaknesses, providing guidance on its application in metabolic engineering, and highlighting the need for methodological refinements to enhance its predictive capabilities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Fisher for complements: extracting cosmology and neutrino mass from the counts-in-cells PDF

In this work, we comprehensively analyse the cosmology dependence of counts-in-cells statistics. Furthermore, we focus on the shape of the one-point probability distribution function (PDF) of the matter density field at mildly non-linear scales. Based on large-deviation statistics, we parametrize the cosmology dependence of the matter PDF in terms of the linear power spectrum, the growth factor, the spherical collapse dynamics, and the non-linear variance. We extend our formalism to include massive neutrinos, finding that the total matter PDF is highly sensitive to the total neutrino mass Mν and can disentangle it from the clustering amplitude σ 8 . Additionally, using more than a million PDFs extracted from the Quijote simulations, we determine the response of the matter PDF to changing parameters in the νΛCDM model and successfully cross-validate the theoretical model and the simulation measurements. We present the first νΛCDM Fisher forecast for the matter PDF at multiple scales and redshifts, and its combination with the matter power spectrum. We establish that the matter PDF and the matter power spectrum are highly complementary at mildly non-linear scales. The matter PDF is particularly powerful for constraining the matter density Ωm, clustering amplitude σ 8 and the total neutrino mass Mν. Adding the mildly non-linear matter PDF to the mildly non-linear matter power spectrum improves constraints on Ω m by a factor of 5 and σ 8 by a factor of 2 when considering the three lowest redshifts. In our joint analysis of the matter PDF and matter power spectrum at three redshifts, the total neutrino mass is constrained to better than 0.01 eV with a total volume of 6 (Gpc h -1 ) 3 . We discuss how density-split statistics can be used to translate those encouraging results for the matter PDF into realistic observables in galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

PRIMORDIAL GRAVITATIONAL WAVE DETECTABILITY WITH DEEP SMALL-SKY COSMIC MICROWAVE BACKGROUND EXPERIMENTS

We use the Bayesian estimation on direct T - Q - U cosmic microwave background (CMB) polarization maps to forecast errors on the tensor-to-scalar power ratio r, and hence on primordial gravitational waves, as a function of sky coverage f{sub sky}. This map-based likelihood filters the information in the pixel-pixel space into the optimal combinations needed for r detection for cut skies, providing enhanced information over a first-step linear separation into a combination of E, B, and mixed modes, and ignoring the latter. With current computational power and for typical resolutions appropriate for r detection, the large matrix inversions required are accurate and fast. Our simulations explore two classes of experiments, with differing bolometric detector numbers, sensitivities, and observational strategies. One is motivated by a long duration balloon experiment like Spider, with pixel noise {proportional_to}{radical}(f{sub sky}) for a specified observing period. This analysis also applies to ground-based array experiments. We find that, in the absence of systematic effects and foregrounds, an experiment with Spider-like noise concentrating on f{sub sky} {approx} 0.02-0.2 could place a 2{sigma}{sub r} Almost-Equal-To 0.014 boundary ({approx}95% confidence level), which rises to 0.02 with an l-dependent foreground residual left over from an assumed efficient component separation. We contrast this with a Planck-like fixed instrumental noise as f{sub sky} varies, which gives a Galaxy-masked (f{sub sky} = 0.75) 2{sigma}{sub r} Almost-Equal-To 0.015, rising to Almost-Equal-To 0.05 with the foreground residuals. Using as the figure of merit the (marginalized) one-dimensional Shannon entropy of r, taken relative to the first 2003 WMAP CMB-only constraint, gives -2.7 bits from the 2012 WMAP9+ACT+SPT+LSS data, and forecasts of -6 bits from Spider (+ Planck); this compares with up to -11 bits for CMBPol, COrE, and PIXIE post-Planck satellites and -13 bits for a perfectly noiseless cosmic variance limited experiment. We thus confirm the wisdom of the current strategy for r detection of deeply probed patches covering the f{sub sky} minimum-error trough with balloon and ground experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Long‐Term Observations of Turbulence Vertical Velocity Spectra in a Convective Mixed Layer: Dependence on Land‐Surface Forcing in the U.S. Southern Great Plains

Doppler lidar vertical velocity retrievals were analyzed for the scale and structure of mixed-layer turbulence over a 7-year period, on fair-weather warm season days (442 cases) in the U.S. Southern Great Plains. Based on the hypothesized influence of land-surface forcing and updraft size on convective boundary layer clouds, a spectral analysis was performed to quantify effects of surface forcing (surface buoyancy flux, friction velocity, and evaporative fraction) on turbulence scale. Significant (order-of-magnitude) variations in spectral density were found in the energy-production subrange and mesoscale regimes. Integral scale decreased with increasing buoyancy flux (and Monin-Obukhov stability parameter), while spectral density in the energy-production subrange increased, implying a transition to buoyancy-driven cellular structures with narrower updrafts. The influence of stability parameter was limited to the neutral to convective transition, and could not explain the wide variation in spectral density in the mesoscale regime. However, high friction velocity (u* > 0.5 m s –1 ) was associated with larger integral scale (updraft width), and a greater portion of spectral density in the mesoscale regime. Lidar profiles and radar reflectivity for individual cases revealed evidence of roll and wave-like structures contributing to mesoscale variability on high friction velocity days, suggesting shear instabilities on days with wetter land surface conditions and smaller buoyancy flux. The role of friction velocity emphasizes the multivariate nature of surface influences on turbulence, and implies that variance-based turbulence closures may not be adequate for capturing effects of surface forcing on updraft size.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

An a priori DNS study of the shadow-position mixing model

In this study, the modeling of mixing by molecular diffusion is a central aspect for transported probability density function (tPDF) methods. In this paper, the newly-proposed shadow position mixing model (SPMM) is examined, using a DNS database for a temporally evolving di-methyl ether slot jet flame. Two methods that invoke different levels of approximation are proposed to extract the shadow displacement (equivalent to shadow position) from the DNS database. An approach for a priori analysis of the mixing-model performance is developed. The shadow displacement is highly correlated with both mixture fraction and velocity, and the peak correlation coefficient of the shadow displacement and mixture fraction is higher than that of the shadow displacement and velocity. This suggests that the composition-space localness is reasonably well enforced by the model, with appropriate choices of model constants. The conditional diffusion of mixture fraction and major species from DNS and from SPMM are then compared, using mixing rates that are derived by matching the mixture fraction scalar dissipation rates. Good qualitative agreement is found, for the prediction of the locations of zero and maximum/minimum conditional diffusion locations for mixture fraction and individual species. Similar comparisons are performed for DNS and the IECM (interaction by exchange with the conditional mean) model. The agreement between SPMM and DNS is better than that between IECM and DNS, in terms of conditional diffusion iso-contour similarities and global normalized residual levels. It is found that a suitable value for the model constant c that controls the mixing frequency can be derived using the local normalized scalar variance, and that the model constant a controls the localness of the model. A higher-Reynolds-number test case is anticipated to be more appropriate to evaluate the mixing models, and stand-alone transported PDF simulations are required to more fully enforce localness and to assess model performance.

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

ON THE NON-KOLMOGOROV NATURE OF FLARE-PRODUCTIVE SOLAR ACTIVE REGIONS

A magnetic power spectral analysis is performed on 53 solar active regions, observed from 2011 August to 2012 July. Magnetic field data obtained from the Helioseismic and Magnetic Imager, inverted as Active Region Patches, are used to study the evolution of the magnetic power index as each region rotates across the solar disk. Active regions are classified based on the numbers and sizes of solar flares they produce in order to study the relationship between flare productivity and the magnetic power index. The choice of window size and inertial range plays a key role in determining the correct magnetic power index. The overall distribution of magnetic power indices has a range of 1.0–2.5. Flare-quiet regions peak at a value of 1.6. However, flare-productive regions peak at a value of 2.2. Overall, the histogram of the distribution of power indices of flare-productive active regions is well separated from flare-quiet active regions. Only 12% of flare-quiet regions exhibit an index greater than 2, whereas 90% of flare-productive regions exhibit an index greater than 2. Flare-quiet regions exhibit a high temporal variance (i.e., the index fluctuates between high and low values), whereas flare-productive regions maintain an index greater than 2 for several days. This shows the importance of including the temporal evolution of active regions in flare prediction studies, and highlights the potential of a 2–3 day prediction window for space weather applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

MOOSE ProbML: Parallelizable Probabilistic Machine Learning and Uncertainty Quantification Capabilities

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a widely used open- source finite element software for performing multiphysics multiscale simulations in a massively parallel fashion. Recently, the computational team at Idaho National Laboratory (INL) has implemented Probabilistic Machine Learning (ProbML) capabilities in MOOSE—in a parallelized fashion—and enable active learning with large-scale computational models for tasks such as surrogate model development, scale bridging, forward/inverse uncertainty quantification (UQ), Bayesian optimization, etc. This presentation summarizes these developments in MOOSE along with demonstrations on several real applications relevant to nuclear energy. At the fundamental level, samplers like Monte Carlo/Latin Hypercube, variance reduction, parallelized Markov Chain Monte Carlo (MCMC) support uncertainty propagation in both forward and inverse settings. These samplers can be integrated with the Gaussian processes (GP) suite in MOOSE, which offer several variants like scalar GPs, multi-output GPs, and deep GPs, to enable active learning. These GPs can be tuned using gradient-based optimization methods like Adam and its variants or gradient-free methods like the elliptical slice sampler (a variant of MCMC adept under Gaussian settings) for more complex covariance kernels or likelihoods whose gradient computations can be cumbersome. A variety of batch acquisition functions permit parallelized evaluation of the computational model and support different learning objectives with high efficiency like Bayesian inference, global surrogate development, optimization, etc. Furthermore, libtorch integration supports training, evaluation, and re-training of neural networks and other complex machine learning models in active learning settings. The impacts of these developments are shown on several real applications: (1) nuclear fuel inverse UQ and model inadequacy assessment using the Kennedy O’Hagan framework; (2) uncertainty aware surrogate modeling for additive manufacturing to predict field quantities; (3) nuclear reactor rare events analysis; and (4) complex fluid flow prediction using a global surrogate with quantified prediction uncertainty. Finally, the outlook of MOOSE ProbML is discussed for both outer-loop and inner-loop computations in the broad view to accelerate fuels and materials qualification, address gaps in knowledge and data, and assess new reactor/fuel systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Unexpected topology of the temperature fluctuations in the cosmic microwave background

We study the topology generated by the temperature fluctuations of the cosmic microwave background (CMB) radiation, as quantified by the number of components and holes, formally given by the Betti numbers, in the growing excursion sets. We compare CMB maps observed by the Planck satellite with a thousand simulated maps generated according to the $Λ$CDM paradigm with Gaussian distributed fluctuations. The comparison is multi-scale, being performed on a sequence of degraded maps with mean pixel separation ranging from 0.05 to 7.33°. The survey of the CMB over S 2 is incomplete due to obfuscation effects by bright point sources and other extended foreground objects like our own galaxy. To deal with such situations, where analysis in the presence of “masks” is of importance, we introduce the concept of relative homology. The parametric χ 2 -test shows differences between observations and simulations, yielding p -values at percent to less than permil levels roughly between 2 and 7°, with the difference in the number of components and holes peaking at more than 3 σ sporadically at these scales. The highest observed deviation between the observations and simulations for b 0 and b 1 is approximately between 3 σ and 4 σ at scales of 3–7°. There are reports of mildly unusual behaviour of the Euler characteristic at 3.66° in the literature, computed from independent measurements of the CMB temperature fluctuations by Planck ’s predecessor, the Wilkinson Microwave Anisotropy Probe (WMAP) satellite. The mildly anomalous behaviour of the Euler characteristic is phenomenologically related to the strongly anomalous behaviour of components and holes, or the zeroth and first Betti numbers, respectively. Further, since these topological descriptors show consistent anomalous behaviour over independent measurements of Planck and WMAP, instrumental and systematic errors may be an unlikely source. These are also the scales at which the observed maps exhibit low variance compared to the simulations, and approximately the range of scales at which the power spectrum exhibits a dip with respect to the theoretical model. Non-parametric tests show even stronger differences at almost all scales. Crucially, Gaussian simulations based on power-spectrum matching the characteristics of the observed dipped power spectrum are not able to resolve the anomaly. Understanding the origin of the anomalies in the CMB, whether cosmological in nature or arising due to late-time effects, is an extremely challenging task. Regardless, beyond the trivial possibility that this may still be a manifestation of an extreme Gaussian case, these observations, along with the super-horizon scales involved, may motivate the study of primordial non-Gaussianity. Alternative scenarios worth exploring may be models with non-trivial topology, including topological defect models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Megavoltage Computed Tomography Image Guidance With Helical Tomotherapy in Patients With Vertebral Tumors: Analysis of Factors Influencing Interobserver Variability

Purpose: To evaluate megavoltage computed tomography (MVCT)-based image guidance with helical tomotherapy in patients with vertebral tumors by analyzing factors influencing interobserver variability, considered as quality criterion of image guidance. Methods and Materials: Five radiation oncologists retrospectively registered 103 MVCTs in 10 patients to planning kilovoltage CTs by rigid transformations in 4 df. Interobserver variabilities were quantified using the standard deviations (SDs) of the distributions of the correction vector components about the observers' fraction mean. To assess intraobserver variabilities, registrations were repeated after {>=}4 weeks. Residual deviations after setup correction due to uncorrectable rotational errors and elastic deformations were determined at 3 craniocaudal target positions. To differentiate observer-related variations in minimizing these residual deviations across the 3-dimensional MVCT from image resolution effects, 2-dimensional registrations were performed in 30 single transverse and sagittal MVCT slices. Axial and longitudinal MVCT image resolutions were quantified. For comparison, image resolution of kilovoltage cone-beam CTs (CBCTs) and interobserver variability in registrations of 43 CBCTs were determined. Results: Axial MVCT image resolution is 3.9 lp/cm. Longitudinal MVCT resolution amounts to 6.3 mm, assessed as full-width at half-maximum of thin objects in MVCTs with finest pitch. Longitudinal CBCT resolution is better (full-width at half-maximum, 2.5 mm for CBCTs with 1-mm slices). In MVCT registrations, interobserver variability in the craniocaudal direction (SD 1.23 mm) is significantly larger than in the lateral and ventrodorsal directions (SD 0.84 and 0.91 mm, respectively) and significantly larger compared with CBCT alignments (SD 1.04 mm). Intraobserver variabilities are significantly smaller than corresponding interobserver variabilities (variance ratio [VR] 1.8-3.1). Compared with 3-dimensional registrations, 2-dimensional registrations have significantly smaller interobserver variability in the lateral and ventrodorsal directions (VR 3.8 and 2.8, respectively) but not in the craniocaudal direction (VR 0.75). Conclusion: Tomotherapy image guidance precision is affected by image resolution and residual deviations after setup correction. Eliminating the effect of residual deviations yields small interobserver variabilities with submillimeter precision in the axial plane. In contrast, interobserver variability in the craniocaudal direction is dominated by the poorer longitudinal MVCT image resolution. Residual deviations after image guidance exist and need to be considered when dose gradients ultimately achievable with image guided radiation therapy techniques are analyzed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Extreme Rainfall Variability in Australia: Patterns, Drivers, and Predictability*

Leading patterns of observed monthly extreme rainfall variability in Australia are examined using an empirical orthogonal teleconnection (EOT) method. Extreme rainfall variability is more closely related to mean rainfall variability during austral summer than in winter. The leading EOT patterns of extreme rainfall explain less variance in Australia-wide extreme rainfall than is the case for mean rainfall EOTs. The authors illustrate that, as with mean rainfall, the El Niño–Southern Oscillation (ENSO) has the strongest association with warm-season extreme rainfall variability, while in the cool season the primary drivers are atmospheric blocking and the subtropical ridge. The Indian Ocean dipole and southern annular mode also have significant relationships with patterns of variability during austral winter and spring. Leading patterns of summer extreme rainfall variability have predictability several months ahead from Pacific sea surface temperatures (SSTs) and as much as a year in advance from Indian Ocean SSTs. Predictability from the Pacific is greater for wetter-than-average summer months than for months that are drier than average, whereas for the Indian Ocean the relationship has greater linearity. Several cool-season EOTs are associated with midlatitude synoptic-scale patterns along the south and east coasts. These patterns have common atmospheric signatures denoting moist onshore flow and strong cyclonic anomalies often to the north of a blocking anticyclone. Tropical cyclone activity is observed to have significant relationships with some warm-season EOTs. This analysis shows that extreme rainfall variability in Australia can be related to remote drivers and local synoptic-scale patterns throughout the year.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

LEISURE: Load-Balanced Network-Wide Traffic Measurement and Monitor Placement

We report that network-wide traffic measurement is of interest to network operators to uncover global network behavior for the management tasks of traffic accounting, debugging or troubleshooting, security, and traffic engineering. Increasingly, sophisticated network measurement tasks such as anomaly detection and security forensic analysis are requiring in-depth fine-grained flow-level measurements. However, performing in-depth per-flow measurements (e.g., detailed payload analysis) is often an expensive process. Given the fast-changing Internet traffic landscape and large traffic volume, a single monitor is not capable of accomplishing the measurement tasks for all applications of interest due to its resource constraint. Moreover, uncovering global network behavior requires network-wide traffic measurements at multiple monitors across the network since traffic measured at any single monitor only provides a partial view and may not be sufficient or accurate. These factors call for coordinated measurements among multiple distributed monitors. In this paper, we present a centralized optimization framework, LEISURE (Load- EqualIzed meaSUREment), for load-balancing network measurement workloads across distributed monitors. Specifically, we consider various load-balancing problems under different objectives and study their extensions to support both fixed and flexible monitor deployment scenarios. We formulate the latter flexible monitor deployment case as an MILP (Mixed Integer Linear Programming) problem and propose several heuristic algorithms to approximate the optimal solution and reduce the computation complexity. We evaluate LEISURE via detailed simulations on Abilene and GEANT network traces to show that LEISURE can achieve much better load-balanced performance (e.g., 4.75 × smaller peak workload and 70 × smaller variance in workloads) across all coordinated monitors in comparison to a naive solution (uniform assignment) to accomplish network-wide traffic measurement tasks under the fixed monitor deployment scenario. Lastly, we also show that under the flexible monitor deployment setting, our heuristic solutions can achieve almost the same load-balancing performance as the optimal solution while reducing the computation times by a factor up to 22.5 × in Abilene and 800 × in GEANT.

42 ENGINEERING↗
DOE OSTI2012

The full squeezed CMB bispectrum from inflation

The small-scale CMB temperature we observe on the sky is modulated by perturbations that were super-horizon at recombination, giving differential focussing and lensing that generate a non-zero bispectrum even for single-field inflation where local physics is identical. Understanding this signal is important for primordial non-Gaussianity studies and also parameter constraints from the CMB lensing bispectrum signal. Because of cancellations individual effects can appear larger or smaller than they are in total, so a full analysis may be required to avoid biases. I relate angular scales on the sky to physical scales at recombination using the optical equations, and give full-sky results for the large-scale adiabatic temperature bispectrum from Ricci focussing (expansion of the ray bundle), Weyl lensing (convergence and shear), and temperature redshift modulations of small-scale power. The δN expansion of the beam is described by the constant temperature 3-curvature, and gives a nearly-observable version of the consistency relation prediction from single-field inflation. I give approximate arguments to quantify the likely importance of dynamical effects, and argue that they can be neglected for modulation scales l∼<100, which is sufficient for lensing studies and also allows robust tests of local primordial non-Gaussianity using only the large-scale modulation modes. For accurate numerical results early and late-time ISW effects must be accounted for, though I confirm that the late-time non-linear Rees-Sciama contribution is negligible compared to other more important complications. The total corresponds to f{sub NL} ∼ 7 for Planck-like temperature constraints and f{sub NL} ∼ 11 for cosmic-variance limited data to l{sub max} = 2000. Temperature lensing bispectrum estimates are affected at the 0.2σ level by Ricci focussing, and up to 0.5σ with polarization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Exploration of multifidelity UQ sampling strategies for computer network applications.

Network modeling is a powerful tool to enable rapid analysis of complex systems that can be challenging to study directly using physical testing. Here, two approaches are considered: emulation and simulation. The former runs real software on virtualized hardware, while the latter mimics the behavior of network components and their interactions in software. Although emulation provides an accurate representation of physical networks, this approach alone cannot guarantee the characterization of the system under realistic operative conditions. Operative conditions for physical networks are often characterized by intrinsic variability (payload size, packet latency, etc.) or a lack of precise knowledge regarding the network configuration (bandwidth, delays, etc.); therefore uncertainty quantification (UQ) strategies should be also employed. UQ strategies require multiple evaluations of the system with a number of evaluation instances that roughly increases with the problem dimensionality, i.e., the number of uncertain parameters. It follows that a typical UQ workflow for network modeling based on emulation can easily become unattainable due to its prohibitive computational cost. In this paper, a multifidelity sampling approach is discussed and applied to network modeling problems. The main idea is to optimally fuse information coming from simulations, which are a low-fidelity version of the emulation problem of interest, in order to decrease the estimator variance. By reducing the estimator variance in a sampling approach it is usually possible to obtain more reliable statistics and therefore a more reliable system characterization. Several network problems of increasing difficulty are presented. For each of them, the performance of the multifidelity estimator is compared with respect to the single fidelity counterpart, namely, Monte Carlo sampling. For all the test problems studied in this work, the multifidelity estimator demonstrated an increased efficiency with respect to MC.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Range cost-effectiveness of plug-in electric vehicle for heterogeneous consumers: An expanded total ownership cost approach

Plug-in electric vehicles (PEV) appears to have sales momentum in major personal vehicle markets but are still at the early market stage. Opportunities to accelerate PEV adoption can be discovered through comprehensive total cost of ownership (TCO) analysis. Understanding the cost-effective electric ranges of PEVs for consumers, manufacturers, and the society is critical for any discussion of PEV mass markets. This study expanded the traditional TCO approach by (1) fully considering heterogeneous consumer driving patterns, (2) quantifying the charging inconvenience and range anxiety cost of battery-electric vehicles (BEVs), and (3) monetizing both tangible and intangible PEV policies. Uncertainties were handled through Monte Carlo simulation. The results suggest that BEVs with an electric range of 250–350 km have the lowest TCO in cities with government-enacted purchase limitations, and internal combustion engine vehicles (ICEVs) have the lowest TCO in cities without purchase limitations, even when considering PEV subsidies. The lowest TCO for some consumer groups is obtained by BEVs with an electric range of 400–450 km, especially in northern China, where the weather is colder. The cost-effective all-electric range for BEVs in each city in 2025 will decrease due to improved battery performance in cold environments and an expanded charging infrastructure. Finally, based on TCO, plug-in hybrid electric vehicles (PHEVs) are currently more suitable for drivers with a high average daily mileage or a large mileage variance. However, by 2025, BEVs with a long driving range may become a more cost-effective choice for these drivers.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2014

Temporal and spatial turbulent spectra of MHD plasma and an observation of variance anisotropy

The nature of magnetohydrodynamic (MHD) turbulence is analyzed through both temporal and spatial magnetic fluctuation spectra. A magnetically turbulent plasma is produced in the MHD wind tunnel configuration of the Swarthmore Spheromak Experiment. The power of magnetic fluctuations is projected into directions perpendicular and parallel to a local mean field; the ratio of these quantities shows the presence of variance anisotropy which varies as a function of frequency. Comparisons among magnetic, velocity, and density spectra are also made, demonstrating that the energy of the turbulence observed is primarily seeded by magnetic fields created during plasma production. Direct spatial spectra are constructed using multi-channel diagnostics and are used to compare to frequency spectra converted to spatial scales using the Taylor hypothesis. Evidence for the observation of dissipation due to ion inertial length scale physics is also discussed, as well as the role laboratory experiments can play in understanding turbulence typically studied in space settings such as the solar wind. Finally, all turbulence results are shown to compare fairly well to a Hall-MHD simulation of the experiment.

79 ASTRONOMY AND ASTROPHYSICS↗