Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical error”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

High-Resolution Regional Wave Hindcast for Hawaii

This report summarizes modeling efforts for hindcasting of wave climate within the exclusive economic zone (EEZ) around the state of Hawaii. Specifically, the report provides detailed description of model mesh development, sensitivity analysis of the configuration of model parameters, data used for model forcing and validation, model skill assessment, and model outputs for resource characterization. The unstructured, nested-grid modeling approach incorporates a regional high-resolution SWAN (Simulating WAves Nearshore) model within the global-regional WAVEWATCH III® (WW3) wave model. The models were configured to comply with the International Electrotechnical Commission (IEC) Technical Standards for wave energy resource characterization. The SWAN model domain covers the entire EEZ with a spatial resolution of approximately 300 m in the nearshore region. Two levels of WW3 nested grids, from 0.5 arc-degree at global to 7.5 arc-minute at regional scale, were used to provide boundary conditions that drive the high-resolution nearshore SWAN model. The WW3 and SWAN models were forced by 0.5 arc-degree Climate Forecast System Reanalysis wind and 5 km regional Weather Research and Forecasting wind, respectively. Wave resource and bulk wave parameters were simulated for a 32-year period from 1979 to 2010, and were subsequently validated using wave buoy data within the model domain. Model performance was evaluated using a set of error statistics. Overall, the model shows good skills in reproducing the wave climate in the Hawaii EEZ in both spatial and temporal domains. Comparisons of time series and scatter plots for the six IEC resource parameters and the associated error statistics tables, and horizontal two-dimensional plots of annually averaged and monthly averaged IEC parameters are provided.

16 TIDAL AND WAVE POWER↗
DOE OSTI2016

Controlling excited-state contamination in nucleon matrix elements

We present a detailed analysis of methods to reduce statistical errors and excited-state contamination in the calculation of matrix elements of quark bilinear operators in nucleon states. All the calculations were done on a 2+1-flavor ensemble with lattices of size 32 3 × 64 generated using the rational hybrid Monte Carlo algorithm at 𝑎 = 0.081 fm and with 𝑀 𝜋 = 312 MeV. The statistical precision of the data is improved using the all-mode-averaging method. We compare two methods for reducing excited-state contamination: a variational analysis and a 2-state fit to data at multiple values of the source-sink separation 𝑡 sep . We show that both methods can be tuned to significantly reduce excited-state contamination and discuss their relative advantages and cost effectiveness. Here a detailed analysis of the size of source smearing used in the calculation of quark propagators and the range of values of 𝑡 sep needed to demonstrate convergence of the isovector charges of the nucleon to the 𝑡 sep → ∞ estimates is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Adiabatic expansion cooling of antihydrogen

Magnetically trapped antihydrogen atoms can be cooled by expanding the volume of the trap in which they are confined. We report a proof-of-principle experiment in which antiatoms are deliberately released from expanded and static traps. Antiatoms escape at an average trap depth of 0.08 ± 0.01 K (statistical errors only) from the expanded trap while they escape at average depths of 0.22 ± 0.01 and 0.17 ± 0.01 K from two different static traps. (We employ temperature-equivalent energy units.) Detailed simulations qualitatively agree with the escape times measured in the experiment and show a decrease of 38 % (statistical error < 0.2 % ) in the mean energy of the population after the trap expansion without significantly increasing antiatom loss compared to typical static confinement protocols. This change is bracketed by the predictions of one-dimensional and three-dimensional semianalytic adiabatic expansion models. These experimental, simulational, and model results are consistent with obtaining an adiabatically cooled population of antihydrogen atoms that partially exchanged energy between axial and transverse degrees of freedom during the trap expansion. This result is important for future antihydrogen gravitational experiments which rely on adiabatic cooling, and it will enable antihydrogen cooling beyond the fundamental limits of laser cooling. Published by the American Physical Society 2024

Ahmadi, M.↗
DOE OSTI2015

Variance reduction for Fokker–Planck based particle Monte Carlo schemes

Recently, Fokker–Planck based particle Monte Carlo schemes have been proposed and evaluated for simulations of rarefied gas flows [1–3]. In this paper, the variance reduction for particle Monte Carlo simulations based on the Fokker–Planck model is considered. First, deviational based schemes were derived and reviewed, and it is shown that these deviational methods are not appropriate for practical Fokker–Planck based rarefied gas flow simulations. This is due to the fact that the deviational schemes considered in this study lead either to instabilities in the case of two-weight methods or to large statistical errors if the direct sampling method is applied. Motivated by this conclusion, we developed a novel scheme based on correlated stochastic processes. The main idea here is to synthesize an additional stochastic process with a known solution, which is simultaneously solved together with the main one. By correlating the two processes, the statistical errors can dramatically be reduced; especially for low Mach numbers. To assess the methods, homogeneous relaxation, planar Couette and lid-driven cavity flows were considered. For these test cases, it could be demonstrated that variance reduction based on parallel processes is very robust and effective.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Biases on cosmological parameter estimators from galaxy cluster number counts

Sunyaev-Zel'dovich (SZ) surveys are promising probes of cosmology — in particular for Dark Energy (DE) —, given their ability to find distant clusters and provide estimates for their mass. However, current SZ catalogs contain tens to hundreds of objects and maximum likelihood estimators may present biases for such sample sizes. In this work we study estimators from cluster abundance for some cosmological parameters, in particular the DE equation of state parameter w{sub 0}, the amplitude of density fluctuations σ{sub 8}, and the Dark Matter density parameter Ω{sub c}. We begin by deriving an unbinned likelihood for cluster number counts, showing that it is equivalent to the one commonly used in the literature. We use the Monte Carlo approach to determine the presence of bias using this likelihood and study its behavior with both the area and depth of the survey, and the number of cosmological parameters fitted. Our fiducial models are based on the South Pole Telescope (SPT) SZ survey. Assuming perfect knowledge of mass and redshift some estimators have non-negligible biases. For example, the bias of σ{sub 8} corresponds to about 40% of its statistical error bar when fitted together with Ω{sub c} and w{sub 0}. Including a SZ mass-observable relation decreases the relevance of the bias, for the typical sizes of current SZ surveys. Considering a joint likelihood for cluster abundance and the so-called ''distance priors'', we obtain that the biases are negligible compared to the statistical errors. However, we show that the biases from SZ estimators do not go away with increasing sample sizes and they may become the dominant source of error for an all sky survey at the SPT sensitivity. Finally, we compute the confidence regions for the cosmological parameters using Fisher matrix and profile likelihood approaches, showing that they are compatible with the Monte Carlo ones. The results of this work validate the use of the current maximum likelihood methods for present SZ surveys, but highlight the need for further studies for upcoming experiments. To perform the analyses of this work, we developed fast, accurate, and adaptable codes for cluster counts in the framework of the Numerical Cosmology Library.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Dark Energy Survey Year 1 Results: Multi-Probe Methodology and Simulated Likelihood Analyses

We present the methodology for and detail the implementation of the Dark Energy Survey (DES) 3x2pt DES Year 1 (Y1) analysis, which combines configuration-space two-point statistics from three different cosmological probes: cosmic shear, galaxy-galaxy lensing, and galaxy clustering, using data from the first year of DES observations. We have developed two independent modeling pipelines and describe the code validation process. We derive expressions for analytical real-space multi-probe covariances, and describe their validation with numerical simulations. We stress-test the inference pipelines in simulated likelihood analyses that vary 6-7 cosmology parameters plus 20 nuisance parameters and precisely resemble the analysis to be presented in the DES 3x2pt analysis paper, using a variety of simulated input data vectors with varying assumptions. We find that any disagreement between pipelines leads to changes in assigned likelihood $\Delta \chi^2 \le 0.045$ with respect to the statistical error of the DES Y1 data vector. We also find that angular binning and survey mask do not impact our analytic covariance at a significant level. We determine lower bounds on scales used for analysis of galaxy clustering (8 Mpc$~h^{-1}$) and galaxy-galaxy lensing (12 Mpc$~h^{-1}$) such that the impact of modeling uncertainties in the non-linear regime is well below statistical errors, and show that our analysis choices are robust against a variety of systematics. These tests demonstrate that we have a robust analysis pipeline that yields unbiased cosmological parameter inferences for the flagship 3x2pt DES Y1 analysis. We emphasize that the level of independent code development and subsequent code comparison as demonstrated in this paper is necessary to produce credible constraints from increasingly complex multi-probe analyses of current data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The need for accurate redshifts in supernova cosmology

Recent papers have shown that a small systematic redshift shift (Δ z ∼ 10{sup −5}) in measurements of type Ia supernovae can cause a significant bias (∼1%) in the recovery of cosmological parameters. Such a redshift shift could be caused, for example, by a gravitational redshift due to the density of our local environment. The sensitivity of supernova data to redshift shifts means supernovae make excellent probes of inhomogeneities. We therefore invert the analysis, and try to diagnose the nature of our local gravitational environment by fitting for Δ z as an extra free parameter alongside the usual cosmological parameters. Using the Joint Light-curve SN Ia dataset we find the best fit includes a systematic redshift shift of Δ z = (2.6{sup +2.7}{sub −2.8}) × 10{sup −4}. This is a larger shift than would be expected due to gravitational redshifts in a standard Λ-Cold Dark Matter universe (though still consistent with zero), and would correspond to a monopole Doppler shift of about 100 km s{sup −1} moving away from the Milky-Way. However, since most supernova measurements are made to a redshift precision of no better than 10{sup −3}, it is possible that a systematic error smaller than the statistical error remains in the data and is responsible for the shift; or that it is an insignificant statistical fluctuation. We find that when Δ z is included as a free parameter while fitting to the JLA SN Ia data, the constraints on the matter density shifts to Ω {sub m} = 0.313{sup +0.042}{sub −0.040}, bringing it into better agreement with the CMB cosmological parameter constraints from Planck. A positive Δ z ∼ 2.6×10{sup −4} would also cause us to overestimate the supernova measurement of Hubble's constant by Δ H {sub 0} ∼ 1 kms{sup −1}Mpc{sup −1}. However this overestimation should diminish as one increases the low-redshift cutoff, and this is not seen in the most recent data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Insights From Dayflow: A Historical Streamflow Reanalysis Dataset for the Conterminous United States

Abstract Reconstructed historical streamflow time series can supplement limited streamflow gauge observations. However, there are common challenges of typical modeling approaches: process‐based hydrologic models can be data/computation‐intensive, and statistics‐based models can be region/stream‐specific. Here we present a nationally scalable modeling framework integrating the simulated runoff from the Variable Infiltration Capacity (VIC) model with the Routing Application for Parallel computatIon of Discharge (RAPID) routing model leveraging high‐performance computing. We demonstrate an efficient method of assimilating streamflow at US Geological Survey (USGS) streamflow monitoring sites using a simple hierarchical approach in the VIC‐RAPID framework. The result is a reconstructed 36‐year (1980–2015) daily and monthly streamflow dataset (Dayflow) at ∼2.7 million NHDPlusV2 stream reaches in the conterminous US (CONUS). We perform a comprehensive evaluation at 7,526 USGS sites and characterize their error statistics. The results demonstrate that 49% of the USGS sites demonstrate Kling–Gupta Efficiency (KGE) > 0.5 and 58% of the sites show percentage bias within ±20% for the daily naturalized streamflow. Streamflow data assimilation across CONUS shows an overall improvement over naturalized streamflow, notably in the western semiarid‐to‐arid regions. Comparison to other national and global streamflow reanalysis datasets such as the National Water Model and Global Reach‐scale A priori Discharge Estimates for SWOT demonstrates improved KGE, reduced bias, and directions for Dayflow improvements. Investigations of error statistics with key hydrologic, hydroclimatic, and geomorphologic basin characteristics reveal region‐specific patterns which may help improve future framework applications. Overall, Dayflow may enable a better understanding of hydrologic conditions in a changing environment, especially in locations currently not represented by streamflow monitoring networks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Blinded challenge for precision cosmology with large-scale structure: Results from effective field theory for the redshift-space galaxy power spectrum

An accurate theoretical template for the galaxy power spectrum is key for the success of ongoing and future spectroscopic surveys. We examine to what extent the effective field theory (EFT) of large-scale structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a 100 times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge is to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and to pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within subpercent statistical errors corresponding to the total simulation volume.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Baryon acoustic oscillations in the Ly α forest of BOSS DR11 quasars

We report a detection of the baryon acousticoscillation (BAO) feature in the flux-correlation function of the Lyα forest of high-redshift quasars with a statistical significance of five standard deviations. The study uses 137,562 quasars in the redshift range 2.1 ≤ z ≤ 3.5 from the data release 11 (DR11) of the Baryon Oscillation Spectroscopic Survey (BOSS) of SDSS-III. This sample contains three times the number of quasars used in previous studies. The measured position of the BAO peak determines the angular distance, D A (z = 2.34) and expansion rate, H(z = 2.34), both on a scale set by the sound horizon at the drag epoch, r d . We find D A /r d = 11.28 ± 0.65(1σ)$+2.8\atop{-1.2}$(2σ) and D H /r d = 9.18 ± 0.28(1σ) ± 0.6(2σ) where D H = c/H. The optimal combination, ~D$0.7\atop{H}$ D$0.3\atop{A}/r d is determined with a precision of ~2%. For the value r d = 147.4 Mpc, consistent with the cosmic microwave background power spectrum measured by Planck, we find D A (z = 2.34) = 1662 ± 96(1σ) Mpc and H(z = 2.34) = 222 ± 7(1σ) km s -1 Mpc -1 . Tests with mock catalogs and variations of our analysis procedure have revealed no systematic uncertainties comparable to our statistical errors. Our results agree with the previously reported BAO measurement at the same redshift using the quasar-Lyα forest cross-correlation. The autocorrelation and cross-correlation approaches are complementary because of the quite different impact of redshift-space distortion on the two measurements. The combined constraints from the two correlation functions imply values of D A /r d that are 7% lower and 7% higher for D H /r d than the predictions of a flat ΛCDM cosmological model with the best-fit Planck parameters. With our estimated statistical errors, the significance of this discrepancy is ≈2.5σ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Clustering of quasars in SDSS-IV eBOSS: study of potential systematics and bias determination

Here, we study the first year of the eBOSS quasar sample in the redshift range 0.9 Q = 2.45 ± 0.05 at $\bar{z}$=1.55, together with its evolution with redshift. We use this information to determine the minimum mass of the halo hosting the quasars and the characteristic halo mass, which we find to be both independent of redshift within statistical error. Using a recently-measured quasar-luminosity-function we also determine the quasar duty cycle. Lastly, the size of this first year sample is insufficient to detect any luminosity dependence to quasar clustering and this issue should be further studied with the final ~500,000 eBOSS quasar sample.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

1–2.4 μm Near-IR Spectrum of the Giant Planet β Pictoris b Obtained with the Gemini Planet Imager

Using the Gemini Planet Imager located at Gemini South, we measured the near-infrared (1.0–2.4 μm) spectrum of the planetary companion to the nearby, young star β Pictoris. We compare the spectrum obtained with currently published model grids and with known substellar objects and present the best matching models as well as the best matching observed objects. Comparing the empirical measurement of the bolometric luminosity to evolutionary models, we find a mass of 12.9 ± 0.2 M{sub Jup}, an effective temperature of 1724 ± 15 K, a radius of 1.46 ± 0.01 R{sub Jup}, and a surface gravity of logg=4.18±0.01 [dex] (cgs). The stated uncertainties are statistical errors only, and do not incorporate any uncertainty on the evolutionary models. Using atmospheric models, we find an effective temperature of 1700–1800 K and a surface gravity of logg=3.5–4.0 [dex] depending upon the model. These values agree well with other publications and with “hot-start” predictions from planetary evolution models. Further, we find that the spectrum of β Pic b best matches a low surface gravity L2 ± 1 brown dwarf. Finally, comparing the spectrum to field brown dwarfs, we find the the spectrum best matches 2MASS J04062677–381210 and 2MASS J03552337+1133437.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Integrating a reservoir regulation scheme into a spatially distributed hydrological model

During the past several decades, numerous reservoirs have been built across the world for a variety of purposes such as flood control, irrigation, municipal/industrial water supplies, and hydropower generation. Consequently, natural streamflow timing and magnitude have been altered significantly by reservoir operations. In addition, the hydrological cycle can be modified by land use/land cover and climate changes. To understand the fine scale feedback between hydrological processes and water management decisions, a distributed hydrological model embedded with a reservoir component is of desire. In this study, a multi-purpose reservoir module with predefined complex operational rules was integrated into the Distributed Hydrology Soil Vegetation Model (DHSVM). Conditional operating rules, which are designed to reduce flood risk and enhance water supply reliability, were adopted in this module. The performance of the integrated model was tested over the upper Brazos River Basin in Texas, where two U.S. Army Corps of Engineers reservoirs, Lake Whitney and Aquilla Lake, are located. The integrated DHSVM model was calibrated and validated using observed reservoir inflow, outflow, and storage data. The error statistics were summarized for both reservoirs on a daily, weekly, and monthly basis. Using the weekly reservoir storage for Lake Whitney as an example, the coefficients of determination (R2) and the Nash-Sutcliff Efficiency (NSE) are 0.85 and 0.75, respectively. These results suggest that this reservoir module has promise for use in sub-monthly hydrological simulations. Enabled with the new reservoir component, the DHSVM model provides a platform to support adaptive water resources management under the impacts of evolving anthropogenic activities and substantial environmental changes.

Zhao, Gang↗
DOE OSTI2022

Angular systematics-free cosmological analysis of galaxy clustering in configuration space

ABSTRACT Galaxy redshift surveys are subject to incompleteness and inhomogeneous sampling due to the various constraints inherent to spectroscopic observations. This can introduce systematic errors on the summary statistics of interest, which need to be mitigated in cosmological analysis to achieve high accuracy. Standard practices involve applying weighting schemes based on completeness estimates across the survey footprint, possibly supplemented with additional weighting schemes accounting for density-dependent effects. In this work, we concentrate on pure angular systematics and describe an alternative approach consisting in analysing the galaxy two-point correlation function where angular modes are nulled. By construction, this procedure removes all possible known and unknown sources of angular observational systematics, but also part of the cosmological signal. We use a modified Landy–Szalay estimator for the two-point correlation function that relies on an additional random catalogue where angular positions are randomly drawn from the galaxy catalogue, and provide an analytical model to describe this modified statistic. We test the model by performing an analysis of the full anisotropic clustering in mock catalogues of luminous red and emission-line galaxies at 0.43 < z < 1.1. We find that the model fully accounts for the modified correlation function in redshift space, without introducing new nuisance parameters. The derived cosmological parameters from the analysis of baryon acoustic oscillations and redshift-space distortions display slightly larger statistical uncertainties, mostly for the growth rate of structure parameter fσ8 that exhibits a $50{{\ \rm per\ cent}}$ statistical error increase, but free from angular systematic error.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Precision cosmology with time delay lenses: High resolution imaging requirements

Lens time delays are a powerful probe of cosmology, provided that the gravitational potential of the main deflector can be modeled with sufficient precision. Recent work has shown that this can be achieved by detailed modeling of the host galaxies of lensed quasars, which appear as ``Einstein Rings'' in high resolution images. The distortion of these arcs and counter-arcs, as measured over a large number of pixels, provides tight constraints on the difference between the gravitational potential between the quasar image positions, and thus on cosmology in combination with the measured time delay. We carry out a systematic exploration of the high resolution imaging required to exploit the thousands of lensed quasars that will be discovered by current and upcoming surveys with the next decade. Specifically, we simulate realistic lens systems as imaged by the Hubble Space Telescope (HST), James Webb Space Telescope (JWST), and ground based adaptive optics images taken with Keck or the Thirty Meter Telescope (TMT). We compare the performance of these pointed observations with that of images taken by the Euclid (VIS), Wide-Field Infrared Survey Telescope (WFIRST) and Large Synoptic Survey Telescope (LSST) surveys. We use as our metric the precision with which the slope γ' of the total mass density profile ρ tot ∝ r–γ' for the main deflector can be measured. Ideally, we require that the statistical error on γ' be less than 0.02, such that it is subdominant to other sources of random and systematic uncertainties. We find that survey data will likely have sufficient depth and resolution to meet the target only for the brighter gravitational lens systems, comparable to those discovered by the SDSS survey. For fainter systems, that will be discovered by current and future surveys, targeted follow-up will be required. Furthermore, the exposure time required with upcoming facilitites such as JWST, the Keck Next Generation Adaptive Optics System, and TMT, will only be of order a few minutes per system, thus making the follow-up of hundreds of systems a practical and efficient cosmological probe.

79 ASTRONOMY AND ASTROPHYSICS↗