Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bias correction”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Aridity-Dependent Land Surface Skin Temperature Biases in CMIP5/6

Land surface skin temperature, a critical indicator of climate change, connects the water and energy cycles between the land and the atmosphere. Here, in this study, we evaluate the simulations of land surface skin temperature from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and CMIP6 models with satellite-based datasets and reanalysis. We find systematic cold skin temperature biases over arid regions in CMIP5/CMIP6 simulations. Over arid and semi-arid regions, latent heat biases drive skin temperature biases by evaporative cooling. Over humid regions, surface downward shortwave and albedo biases are relatively more critical. Spatial patterns of biases remain similar in the latest CMIP6 simulations, suggesting systematic biases in land-atmosphere interactions. These biases need to be corrected or considered while using models for future projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Liquid to crystal Si growth simulation using machine learning force field

Machine learning force field (ML-FF) has emerged as a potential promising approach to simulate various material phenomena for large systems with ab initio accuracy. However, most ML-FFs have been used to study the phenomena relatively close to the equilibrium ground states. In this work, we have studied a far from equilibrium system of liquid to crystal Si growth using ML-FF. Here, we found that our ML-FF based on ab initio decomposed atomic energy can reproduce all the aspects of ab initio simulated growth, from local energy fluctuations to transition temperatures, to diffusion constant, and growth rates. We have also compared the growth simulation with the Stillinger-Weber classical force field and found significant differences. A procedure is also provided to correct a systematic fitting bias in the ML-FF training process, which exists in all training models, otherwise critical results like transition temperature will be wrong.

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

Augmenting a Simulation Campaign for Hybrid Computer Model and Field Data Experiments

The Kennedy and O’Hagan (KOH) calibration framework uses coupled Gaussian processes (GPs) to meta-model an expensive simulator (first GP), tune its “knobs” (calibration inputs) to best match observations from a real physical/field experiment and correct for any modeling bias (second GP) when predicting under new field conditions (design inputs). There are well-established methods for placement of design inputs for data-efficient planning of a simulation campaign in isolation, that is, without field data: space-filling, or via criterion like minimum integrated mean-squared prediction error (IMSPE). Analogues within the coupled GP KOH framework are mostly absent from the literature. Here, in this study, we derive a closed form IMSPE criterion for sequentially acquiring new simulator data for KOH. We illustrate how acquisitions space-fill in design space, but concentrate in calibration space. Closed form IMSPE precipitates a closed-form gradient for efficient numerical optimization. We demonstrate that our KOH-IMSPE strategy leads to a more efficient simulation campaign on benchmark problems, and conclude with a showcase on an application to equilibrium concentrations of rare earth elements for a liquid–liquid extraction reaction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Private Tabular Survey Data Products through Synthetic Microdata Generation

We propose two synthetic microdata approaches to generate private tabular survey data products for public release. We adapt a pseudo posterior mechanism that downweights by-record likelihood contributions with weights ∈[0,1] based on their identification disclosure risks to producing tabular products for survey data. Our method applied to an observed survey database achieves an asymptotic global probabilistic differential privacy guarantee. Our two approaches synthesize the observed sample distribution of the outcome and survey weights, jointly, such that both quantities together possess a privacy guarantee. The privacy-protected outcome and survey weights are used to construct tabular cell estimates (where the cell inclusion indicators are treated as known and public) and associated standard errors to correct for survey sampling bias. Through a real data application to the Survey of Doctorate Recipients public use file and simulation studies motivated by the application, we demonstrate that our two microdata synthesis approaches to construct tabular products provide superior utility preservation as compared to the additive noise approach of the Laplace Mechanism. Moreover, our approaches allow the release of microdata to the public, enabling additional analyses at no extra privacy cost.

Mathematical Methods In Social Sciences↗
DOE OSTI2022

Optical and X-ray GRB Fundamental Planes as cosmological distance indicators

ABSTRACT Gamma-ray bursts (GRBs), can be employed as standardized candles, extending the distance ladder beyond Type Ia supernovae (SNe Ia, z = 2.26). We standardize GRBs using the three-dimensional (3D) Fundamental Plane relation (the Dainotti relation) among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity. Combining SNe Ia and GRBs, we constrain ΩM = 0.299 ± 0.009 assuming a flat Λ cold dark matter (ΛCDM) cosmology with and without correcting GRBs for selection biases and redshift evolution. Using a 3D optical Dainotti correlation, we find this sample is as efficacious in the determination of ΩM as the X-ray sample. We trimmed our GRB samples to achieve tighter planes to simulate additional GRBs. We determined how many GRBs are needed as stand-alone probes to achieve a comparable precision on ΩM to the one obtained by SNe Ia only. We reach the same error measurements derived using SNe Ia in 2011 and 2014 with 142 and 284 simulated optical GRBs, respectively, considering the error bars on the variables halved. These error limits will be reached in 2038 and in 2047, respectively. Using a doubled sample (obtained by future machine learning approaches allowing a light-curve reconstruction and the estimates of GRB redshifts when z is unknown) compared to the current sample, with error bars halved we will reach the same precision as SNe Ia in 2011 and 2014, now and in 2026, respectively. If we consider the current SNe precision, this will be reached with 390 optical GRBs by 2054.

Dainotti, M. G.↗
DOE OSTI2025

Accurate shear estimation with fourth-order moments

ABSTRACT As imaging surveys progress in exploring the large-scale structure of the Universe through the use of weak gravitational lensing, achieving sub-per cent accuracy in estimating shape distortions caused by lensing, or shear, is imperative for precision cosmology. In this paper, we extend the Fourier power function shapelets (FPFS) shear estimator using fourth-order shapelet moments and combine it with the original second-order shear estimator to reduce galaxy shape noise. We calibrate this novel shear estimator analytically to a sub-per cent level-accuracy using the AnaCal framework. This higher order shear estimator is tested with realistic image simulations, and after analytical correction for the detection/selection bias and noise bias, the multiplicative shear bias $|m|$ is below $3\times 10^{-3}$ (99.7 per cent confidence interval) for both isolated and blended galaxies. Once combined with the second-order FPFS shear estimator, the shape noise is reduced by $\sim 35~{{\ \rm per\ cent}}$ for isolated galaxies in simulations with Hyper Suprime-Cam and Vera C. Rubin Observatory Legacy Survey of Space and Time observational conditions. However, for blended galaxies, the effective number density does not significantly improve with the combination of the two estimators. Based on these results, we recommend exploration of how this framework can further reduce the systematic uncertainties in shear due to point spread function leakage and modelling error, and potentially provide improved precision in shear inference in high-resolution space-based images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Many-body electronic structure of LaScO 3 by real-space quantum Monte Carlo

We present real space quantum Monte Carlo (QMC) calculations of the scandate LaSc0 3 that proved to be challenging for traditional electronic structure approaches due to strong correlation effects resulting in inaccurate band gaps from DFT and GW methods when compared with existing experimental data. Besides calculating an accurate QMC band gap corrected for supercell size biases and in agreement with numerous experiments, we also predict the cohesive energy of the crystal using the standard fixed-node QMC without any empirical or non-variational parameters. We show that promotion (optical) gap and fundamental gap agree with each other illustrating a clear absence of significant excitonic effects in the ideal crystal. We obtained these results in perfect consistency in two independent tracks that employ different basis sets (plane wave vs. localized gaussians), different codes for generating orbitals (QUANTUM ESPRESSO vs. CRYSTAL), different QMC codes (QMCPACK vs. QWALK) and different high-accuracy pseudopotentials (ccECPs vs. Troullier-Martins) presenting the maturity and consistency of QMC methodology and tools for studies of strongly correlated problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Mixture density network estimation of continuous variable maximum likelihood using discrete training samples

Abstract Mixture density networks (MDNs) can be used to generate posterior density functions of model parameters $$\varvec{\theta }$$ θ given a set of observables $${\mathbf {x}}$$ x . In some applications, training data are available only for discrete values of a continuous parameter $$\varvec{\theta }$$ θ . In such situations, a number of performance-limiting issues arise which can result in biased estimates. We demonstrate the usage of MDNs for parameter estimation, discuss the origins of the biases, and propose a corrective method for each issue.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Measurements of higher-order cumulants of multiplicity and net-electric charge distributions in inelastic proton-proton interactions by NA61/SHINE

This paper presents the energy dependence of multiplicity and net-electric charge fluctuations in p+p interactions at beam momenta 20, 31, 40, 80, and 158 \(\text{ Ge }\hspace{-1.00006pt}\text{ V }\!/\!c\). Results are corrected for the experimental biases and quantified with the use of cumulants and factorial cumulants. Cumulant ratios are an essential tool in the search for the critical point of strongly interacting matter in heavy ion collisions. Measurements performed in p+p interactions provide a vital baseline estimation in these studies. The measured signals are compared with the string hadronic models Epos1.99 and FTFP-BERT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multiplicity and net-electric charge fluctuations in central Ar+Sc interactions at 13 A , 19 A , 30 A , 40 A , 75 A , and 150 A GeV/c beam momenta measured by NA61/SHINE at the CERN SPS

This paper presents results on multiplicity fluctuations of positively and negatively charged hadrons as well as net-electric charge fluctuations measured in central Ar+Sc interactions at beam momenta 13 A , 19 A , 30 A , 40 A , 75 A , and 150 A GeV/c. The fluctuation analysis is one of the tools to search for the predicted critical point of strongly interacting matter. Results are corrected for the experimental biases and quantified using cumulant ratios. In most instances, multiplicity and net-charge distributions appear narrower than the corresponding Poisson or Skellam distributions. Cumulant ratios are compared with the EPOS1.99 model predictions, which provide a qualitative description that aligns with observations for positively and negatively charged particles. The obtained results are also compared to earlier NA61/SHINE results from inelastic p+p interactions in the same analysis acceptance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude ( M w ) Estimation in Iraq’s Tectonically Diverse Zones

Robust estimation of moment magnitude (M w ) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of M w for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Therefore, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

58 GEOSCIENCES↗
DOE OSTI2020

AWCC Simulations

Neutron well counters provide a means to measure fissile mass by detecting coincident neutrons, which is a unique signature to fission. For coincidence counting, mass is determined using a calibration curve established with items that have the same material characteristics. Multiplicity analysis solves the point model equations for mass using known detector parameters. In both cases, mass measurements may exhibit a bias due to item characteristics and should be corrected for. Experiments can be done to estimate bias and material effects, but they require having both the well counter and the material well defined. In pursuit of estimating measurement bias and uncertainty in a High Efficiency Neutron Counter (HENC), an MCNP simulation was matched to measurements performed with an Active Well Coincidence Counter (AWCC) to de ne a neutron source term. The neutron source term will be used to predict measurement performance of similar material types in a HENC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

A Statistical Standard Siren Measurement of the Hubble Constant from the LIGO/Virgo Gravitational Wave Compact Object Merger GW190814 and Dark Energy Survey Galaxies

We present a measurement of the Hubble constant H 0 using the gravitational wave (GW) event GW190814, which resulted from the coalescence of a 23 M ⊙ black hole with a 2.6 M ⊙ compact object, as a standard siren. No compelling electromagnetic counterpart has been identified for this event; thus our analysis accounts for thousands of potential host galaxies within a statistical framework. The redshift information is obtained from the photometric redshift (photo- z ) catalog from the Dark Energy Survey. The luminosity distance is provided by the LIGO/Virgo gravitational wave sky map. Since this GW event has the second-smallest localization volume after GW170817, GW190814 is likely to provide the best constraint on cosmology from a single standard siren without identifying an electromagnetic counterpart. Furthermore, our analysis uses photo- z probability distribution functions and corrects for photo- z biases. We also reanalyze the binary black hole GW170814 within this updated framework. We explore how our findings impact the H 0 constraints from GW170817, the only GW merger associated with a unique host galaxy. From a combination of GW190814, GW170814, and GW170817, our analysis yields H 0 =72.0$_{+12}^{–8.2}$ km s -1 Mpc -1 (68% highest-density interval, HDI) for a prior in H 0 uniform between [20 and 140] km s -1 Mpc -1 . The addition of GW190814 and GW170814 to GW170817 improves the 68% HDI from GW170817 alone by ~18%, showing how well-localized mergers without counterparts can provide a significant contribution to standard siren measurements, provided that a complete galaxy catalog is available at the location of the event.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Improved laser ranging accuracy through dynamic threshold compensation

During a polar motion experiment (in which two laser ranging systems on a north-south line simultaneously tracked a satellite to determine motion of the pole), techniques for improving both the range noise and range bias were demonstrated. The corrections to the measured ranges (based on received signal levels relative to a fixed correction that would have been applied if pulse height information was not available) are given. It was shown that even existing long pulse systems (20 to 40 nsec duration) can be improved to the 20 to 30 cm level of precision without major or expensive changes to the system.

Johnson, T. S.↗
NASA NTRS1981

Equivalent source modeling of the main field using Magsat data

Magsat dipole solution models were obtained with both 32 deg and 21 deg resolution based on a data set extending over 4 months. Time dependence was modeled using first time derivatives for the dipole magnetization vector components. This doubles the n total number of parameters in the solution. The solutions displayed a very slow convergence in time derivatives although at each iteration the conversion of the dipole parameters to spherical harmonic coefficients showed close agreement with the MGST (12/80) spherical harmonic model to degree 13 in the constant terms and 8 in the secular variation terms. The program error in the option to simultaneously estimate observatory anomaly biases was found and corrected. A 32 deg resolution dipole model was generated using the selected magnetic observatory data from 1960-1977 used in the GSFC (9/80) spherical harmonic model.

Source record↗
NASA NTRS1983

The potential of SAR directional spectra in operational wave forecasting

Digitally processed Seasat SAR ocean wave imagery, when digitally transformed, smoothed, and corrected for known instrument biases, can accurately track multiple (at least three) wave systems in the open ocean, across major currents, and into shallow water. Furthermore, the relative backscatterer of the digitally processed SAR image, when sampled within a single pass and along a constant range, responds to the local wind through a simple power law relationship. The results indicate that SAR spectra, if obtained from a low altitude satellite, could be an important supplement to global winds and non-directional wave heights obtained through other methods.

Beal, R. C.↗
NASA NTRS1987

Relativistic beaming and the nuclei of double-lobed quasars

If the simple relativistic beaming model for extragalactic radio sources is at all correct, Doppler boosting may bias us toward choosing those objects pointing nearly at us in flux-density limited samples. In order to minimize any such effects in statistical tests of physical theories of these objects, an attempt is made to define a complete sample of sources with random orientations. For a sample of 26 double-lobed quasars, the dependence of a number of observable features on the angle to the line of sight and on the superluminal velocity is examined, and the distribution of superluminal velocities for the sources is measured. Two new superluminals are reported, and their apparent velocities and nuclear strengths are discussed.

Hough, David H.↗
NASA NTRS1989

A gyroscope calibration analysis for the Gamma Ray Observatory (GRO)

Current documentation and software do not adequately address the calculation and use of the optimal weight matrices involved in calibrating inertial reference units (IRU). Several facets of the GRO IRU calibration as it relates to the bias and misalignment weighting matrices are investigated. The physical meaning and use of the bias and misalignment weight matrices in IRU calibration are examined. The relation of the weighting and the final biases, misalignments, and their corrections are pursued. Methods for determining reliable, realistic weighting matrices to be used in the GRO IRU calibration (IRUCAL) utility are determined. Possible correlations among observation uncertainties are also explored. For the undetermined case where the maneuvers are insufficient to identify all calibration parameters, the weighting matrices allow as much information as possible to be extracted from the measurements. Finally, applicable simulated flight data are used, incorporating the appropriate calibration maneuvers, to test the weighting matrices in the IRUCAL utility, and examine correlation effects.

Ketchum, Eleanor A.↗