Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bayesian 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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Systematic KMTNet Planetary Anomaly Search. I. OGLE-2019-BLG-1053Lb, a BuriedTerrestrial Planet

In order to exhume the buried signatures of “missing planetary caustics” in Korea Microlensing Telescope Network (KMTNet) data, we conducted a systematic anomaly search of the residuals from point-source point-lens fits, based on a modified version of the KMTNet Event Finder algorithm. This search revealed the lowest-mass-ratio planetary caustic to date in the microlensing event OGLE-2019-BLG-1053, for which the planetary signal had not been noticed before. The planetary system has a planet–host mass ratio ofq= (1.25±0.13) × 10−5. A Bayesian analysis yielded estimates of the mass of the host star, Mhost =-0.61+0.29 -0.24 Mo, the mass of its planet, Mplanet =-2.48 +1.19 -0.98 Mo, the projected planet – host separation, a^= 3.4 +0.5/-0.5 au, and the lens distance, DL =-6.8 +0.6 -0.90kpc.The discovery of this very-low-mass-ratio planet illustrates the utility of our method and opens a new window for a large and homogeneous sample to study the microlensing planet–host mass ratio function down to q∼ 10−5.

Exoplanet detection methods↗
NASA NTRS2021

KMT-2021-BLG-0322: Severe Degeneracy Between Triple-lens and Higher-order Binary-lens Interpretations

We investigate the microlensing event KMT-2021-BLG-0322, for which the light curve exhibits three distinctive sets of caustic-crossing features. It is found that the overall features of the light curve are approximately described by a binary-lens (2L1S) model, but the model leaves substantial residuals. We test various interpretations with the aim of explaining the residuals. Methods. We find that the residuals can be explained either by considering a nonrectilinear lens-source motion caused by the microlens-parallax and lens-orbital effects or by adding a low-mass companion to the binary lens (3L1S model). The degeneracy between the higher-order 2L1S model and the 3L1S model is very severe, making it difficult to single out a correct solution based on the photometric data. This degeneracy was known before for two previous events (MACHO-97-BLG-41 and OGLE-2013-BLG-0723),which led to the false detections of planets in binary systems, and thus the identification of the degeneracy for KMT-2021-BLG-0322illustrates that the degeneracy can be not only common but also very severe, emphasizing the need to check both interpretations of deviations from 2L1S models. Results. From the Bayesian analysis conducted with the measured lensing observables of the event timescale, angular Einstein radius, and microlens parallax, it was estimated that the binary lens components have masses(M1,M2)=(0.62+0.25−0.26M,0.07+0.03−0.03M), for both2L1S and 3L1S solutions, and the mass of the tertiary lens component according to the 3L1S solution isM3=6.40+2.64−2.78MJ.

Gravitational lensing: micro↗
DOE OSTI2012

Global Fits of the Minimal Universal Extra Dimensions Scenario

In theories with Universal Extra-Dimensions (UED), the {gamma}{sub 1} particle, first excited state of the hypercharge gauge boson, provides an excellent Dark Matter (DM) candidate. Here we use a modified version of the SuperBayeS code to perform a Bayesian analysis of the minimal UED scenario, in order to assess its detectability at accelerators and with DM experiments. We derive in particular the most probable range of mass and scattering cross sections off nucleons, keeping into account cosmological and electroweak precision constraints. The consequences for the detectability of the {gamma}{sub 1} with direct and indirect experiments are dramatic. The spin-independent cross section probability distribution peaks at {approx} 10{sup -11} pb, i.e. below the sensitivity of ton-scale experiments. The spin-dependent cross-section drives the predicted neutrino flux from the center of the Sun below the reach of present and upcoming experiments. The only strategy that remains open appears to be direct detection with ton-scale experiments sensitive to spin-dependent cross-sections. On the other hand, the LHC with 1 fb{sup -1} of data should be able to probe the current best-fit UED parameters.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

A Bayesian Meta-Analysis of the Effect of Alcohol Use on HCV-Treatment Outcomes with a Comparison of Resampling Methods to Assess Uncertainty in Parameter Estimates.

There is mounting evidence that alcohol use is significantly linked to lower HCV treatment response rates in interferon-based therapies, though some of the evidence is conflicting. Furthermore, although health care providers recommend reducing or abstaining from alcohol use prior to treatment, many patients do not succeed in doing so. The goal of this meta-analysis was to systematically review and summarize the Englishlanguage literature up through January 30, 2015 regarding the relationship between alcohol use and HCV treatment outcomes, among patients who were not required to abstain from alcohol use in order to receive treatment. Seven pertinent articles studying 1,751 HCV-infected patients were identified. Log-ORs of HCV treatment response for heavy alcohol use and light alcohol use were calculated and compared. We employed a hierarchical Bayesian meta-analytic model to accommodate the small sample size. The summary estimate for the log-OR of HCV treatment response was -0.775 with a 95% credible interval of (-1.397, -0.236). The results of the Bayesian meta-analysis are slightly more conservative compared to those obtained from a boot-strapped, random effects model. We found evidence of heterogeneity (Q = 14.489, p = 0.025), accounting for 60.28% of the variation among log-ORs. Meta-regression to capture the sources of this heterogeneity did not identify any of the covariates investigated as significant. This meta-analysis confirms that heavy alcohol use is associated with decreased HCV treatment response compared to lighter levels of alcohol use. Further research is required to characterize the mechanism by which alcohol use affects HCV treatment response.

alcohol↗
DOE OSTI2011

EVALUATION OF SPRING OPERATED RELIEF VALVE MAINTENANCE INTERVALS AND EXTENSION OF MAINTENANCE TIMES USING A WEIBULL ANALYSIS WITH MODIFIED BAYESIAN UPDATING

The Savannah River Site (SRS) spring operated pressure relief valve (SORV) maintenance intervals were evaluated using an approach provided by the American Petroleum Institute (API RP 581) for risk-based inspection technology (RBI). In addition, the impact of extending the inspection schedule was evaluated using Monte Carlo Simulation (MCS). The API RP 581 approach is characterized as a Weibull analysis with modified Bayesian updating provided by SRS SORV proof testing experience. Initial Weibull parameter estimates were updated as per SRS's historical proof test records contained in the Center for Chemical Process Safety (CCPS) Process Equipment Reliability Database (PERD). The API RP 581 methodology was used to estimate the SORV's probability of failing on demand (PFD), and the annual expected risk. The API RP 581 methodology indicates that the current SRS maintenance plan is conservative. Cost savings may be attained in certain mild service applications that present low PFD and overall risk. Current practices are reviewed and recommendations are made for extending inspection intervals. The paper gives an illustration of the inspection costs versus the associated risks by using API RP 581 Risk Based Inspection (RBI) Technology. A cost effective maintenance frequency balancing both financial risk and inspection cost is demonstrated.

02 PETROLEUM↗
DOE OSTI2014

Developing a Comprehensive Risk Assessment Framework for Geological Storage CO 2

The operational risks for CCS projects include: risks of capturing, compressing, transporting and injecting CO₂; risks of well blowouts; risk that CO 2 will leak into shallow aquifers and contaminate potable water; and risk that sequestered CO 2 will leak into the atmosphere. This report examines these risks by using information on the risks associated with analogue activities such as CO 2 based enhanced oil recovery (CO 2 -EOR), natural gas storage and acid gas disposal. We have developed a new analysis of pipeline risk based on Bayesian statistical analysis. Bayesian theory probabilities may describe states of partial knowledge, even perhaps those related to non-repeatable events. The Bayesian approach enables both utilizing existing data and at the same time having the capability to adsorb new information thus to lower uncertainty in our understanding of complex systems. Incident rates for both natural gas and CO 2 pipelines have been widely used in papers and reports on risk of CO 2 pipelines as proxies for the individual risk created by such pipelines. Published risk studies of CO 2 pipelines suggest that the individual risk associated with CO2 pipelines is between 10 -3 and 10 -4 , which reflects risk levels approaching those of mountain climbing, which many would find unacceptably high. This report concludes, based on a careful analysis of natural gas pipeline failures, suggests that the individual risk of CO 2 pipelines is likely in the range of 10-6 to 10-7, a risk range considered in the acceptable to negligible range in most countries. If, as is commonly thought, pipelines represent the highest risk component of CCS outside of the capture plant, then this conclusion suggests that most (if not all) previous quantitative- risk assessments of components of CCS may be orders of magnitude to high. The potential lethality of unexpected CO 2 releases from pipelines or wells are arguably the highest risk aspects of CO 2 enhanced oil recovery (CO2-EOR), carbon capture, and storage (CCS). Assertions in the CCS literature, that CO 2 levels of 10% for ten minutes, or 20 to 30% for a few minutes are lethal to humans, are not supported by the available evidence. The results of published experiments with animals exposed to CO 2 , from mice to monkeys, at both normal and depleted oxygen levels, suggest that lethal levels of CO 2 toxicity are in the range 50 to 60%. These experiments demonstrate that CO 2 does not kill by asphyxia, but rather is toxic at high concentrations. It is concluded that quantitative risk assessments of CCS have overestimated the risk of fatalities by using values of lethality a factor two to six lower than the values estimated in this paper. In many dispersion models of CO 2 releases from pipelines, no fatalities would be predicted if appropriate levels of lethality for CO 2 had been used in the analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

On the Methodology to Calculate the Covariance of Estimated Resonance Parameters

Principles to determine resonance parameters and their covariance from experimental data are discussed. Different methods to propagate the covariance of experimental parameters are compared. A full Bayesian statistical analysis reveals that the level to which the initial uncertainty of the experimental parameters propagates, strongly depends on the experimental conditions. For high precision data the initial uncertainties of experimental parameters, like a normalization factor, has almost no impact on the covariance of the parameters in case of thick sample measurements and conventional uncertainty propagation or full Bayesian analysis. The covariances derived from a full Bayesian analysis and least-squares fit are derived under the condition that the model describing the experimental observables is perfect. When the quality of the model can not be verified a more conservative method based on a renormalization of the covariance matrix is recommended to propagate fully the uncertainty of experimental systematic effects. Finally, neutron resonance transmission analysis is proposed as an accurate method to validate evaluated data libraries in the resolved resonance region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Statistical Significance Testing for Mixed Priors: A Combined Bayesian and Frequentist Analysis

In many hypothesis testing applications, we have mixed priors, with well-motivated informative priors for some parameters but not for others. The Bayesian methodology uses the Bayes factor and is helpful for the informative priors, as it incorporates Occam’s razor via the multiplicity or trials factor in the look-elsewhere effect. However, if the prior is not known completely, the frequentist hypothesis test via the false-positive rate is a better approach, as it is less sensitive to the prior choice. We argue that when only partial prior information is available, it is best to combine the two methodologies by using the Bayes factor as a test statistic in the frequentist analysis. We show that the standard frequentist maximum likelihood-ratio test statistic corresponds to the Bayes factor with a non-informative Jeffrey’s prior. We also show that mixed priors increase the statistical power in frequentist analyses over the maximum likelihood test statistic. We develop an analytic formalism that does not require expensive simulations and generalize Wilks’ theorem beyond its usual regime of validity. In specific limits, the formalism reproduces existing expressions, such as the p-value of linear models and periodograms. We apply the formalism to an example of exoplanet transits, where multiplicity can be more than 10 7 . We show that our analytic expressions reproduce the $p$-values derived from numerical simulations. We offer an interpretation of our formalism based on the statistical mechanics. We introduce the counting of states in a continuous parameter space using the uncertainty volume as the quantum of the state. We show that both the $p$-value and Bayes factor can be expressed as an energy versus entropy competition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Improved estimates of the Milky Way’s stellar mass and star formation rate from hierarchical Bayesian meta-analysis

We present improved estimates of several global properties of the Milky Way, including its current star formation rate (SFR), the stellar mass contained in its disk and bulge+bar components, as well as its total stellar mass. We do so by combining previous measurements from the literature using a hierarchical Bayesian (HB) statistical method that allows us to account for the possibility that any value may be incorrect or have underestimated errors. We show that this method is robust to a wide variety of assumptions about the nature of problems in individual measurements or error estimates. Ultimately, our analysis yields an SFR for the Galaxy of M-dot {sub ⋆}=1.65±0.19 M{sub ⊙} yr{sup −1}, assuming a Kroupa initial mass function (IMF). By combining HB methods with Monte Carlo simulations that incorporate the latest estimates of the Galactocentric radius of the Sun, R{sub 0}, the exponential scale length of the disk, L{sub d}, and the local surface density of stellar mass, Σ{sub ⋆}(R{sub 0}), we show that the mass of the Galactic bulge+bar is M{sub ⋆}{sup B}=0.91±0.07×10{sup 10} M{sub ⊙}, the disk mass is M{sub ⋆}{sup D}=5.17±1.11×10{sup 10} M{sub ⊙}, and their combination yields a total stellar mass of M{sub ⋆}=6.08±1.14×10{sup 10} M{sub ⊙} (assuming a Kroupa IMF and an exponential disk profile). This analysis is based upon a new compilation of literature bulge mass estimates, normalized to common assumptions about the stellar IMF and Galactic disk properties, presented herein. We additionally find a bulge-to-total mass ratio for the Milky Way of B/T=0.150{sub −0.019}{sup +0.028} and a specific SFR of M-dot {sub ⋆}/M{sub ⋆}=2.71±0.59×10{sup −11} yr{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

High-Temperature Sapphire Pressure Sensors for Harsh Environments

The project aimed to develop and characterize an optical pressure sensor capable of withstanding high temperatures, up to 1000 °C. Survivability in harsh environments was enabled using sapphire for both the substrate and the optical fibers, therefore mitigating thermal expansion mismatch. Such material cannot be machined through conventional microfabrication techniques, and unique additive and subtractive methods had to be studied and implemented. Subtractive methods included laser-machining the sapphire substrate using a picosecond laser machining tool. At first, the parameters of importance were theoretically and experimentally evaluated such that the fluence, the pulse spacing, and the number of passes were found to strongly influence the laser ablation process. A simulation tool was designed to predict the final pattern when inputting a set of parameters and a profile. This model was complemented with laser ablation material physics models and Bayesian uncertainty analysis for one dimensional and three dimensional ablation processes. An additive process called thermocompression bonding was studied in order to seal the diaphragm substrate to the back cavity. At first, bonding strength using this technique was examined using the chevron-notch method that allow to determine the fracture toughness of this bond. This parameter is independent of the geometry used and therefore can be extended to other geometries used in actual sensor bonding. A theoretical and experimental study of internal damage caused by the laser micromachining of the sapphire material was carried out. The research identified enhanced strength and fracture toughness along laser ablated surfaces. This was determined using bend bar specimens and Vicker’s microindentations. Softening of the plasticity behavior (as indicated by nanoindentation experiments) was also investigated using finite deformation plasticity modeling via a contact mechanics finite element simulation. Good correlations between experiments and the model were achieved to determine material parameters that can be used in future MEMS based sensor design and modeling to minimize sensor failure loads. The enhanced strength and fracture toughness was determined to be due to compressive residual stress near laser ablated surfaces. This was confirmed through x-ray spectra taken along different crystal structure planes in combination and Bayesian uncertainty analysis. The Bayesian analysis inferred the entire residual strain tensor (six components) and the residual stress tensor. The results illustrated compressive in-plane stress was the most likely cause of enhanced strength and toughness after laser machining. Due to major equipment failures, a parallel route for fabricating the sensor was taken, using commercially available sapphire removal method (subtractive) and ceramic adhesive bonding (additive). The fully-packaged sensor utilized laser micromachining for the sapphire substrate fabrication and was comprised only of high-temperature-compatible materials. Calibration tests of the sensor in both room-temperature and high-temperature plane-wave tubes showed good linearity and consistent frequency response, in the tested frequency range. The sensor also showed good thermal stability in tests up to 300 °C, proving the fiber-optic lever and its packaging process to be a promising approach for high-temperature measurements. The calibration from 23 °C to 300 °C gave a raw sensitivity of 3.9 × 10 -7 V/Pa and a normalized sensitivity of 9.8 × 10 -8 Pa -1 which were consistent within the temperature range. The sensor had a resonance at 16.3 kHz and a flat band (±3 dB) up to 10.3 kHz. The noise floor of the sensor at 1 kHz was 316 nV/Hz 1/2 , resulting in a MDP of 0.8 Pa. Previously, diaphragm delamination due to CTE mismatch was found at 500 °C, causing packaging failure. Unfortunately high temperature validation in hot jet environment could not be made in time. In the future, it is suggested that the two sapphire substrates be bonded by thermal compression bonding with an intermediate platinum layer instead of using an alumina-based adhesive, and therefore achieve higher temperature capacity up to 1000 °C.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Goal-oriented a-posteriori estimation of model error as an aid to parameter estimation

In this work, a Bayesian model calibration framework is presented that utilizes goal-oriented a-posterior error estimates in quantities of interest (QoIs) for classes of high-fidelity models characterized by PDEs. It is shown that for a large class of computational models, it is possible to develop a computationally inexpensive procedure for calibrating parameters of high-fidelity models of physical events when the parameters of low-fidelity (surrogate) models are known with acceptable accuracy. The main ingredients in the proposed model calibration scheme are goal-oriented a-posteriori estimates of error in QoIs computed using a so-called lower fidelity model compared to those of an uncalibrated higher fidelity model. The estimates of error in QoIs are used to define likelihood functions in Bayesian inversion analysis. A standard Bayesian approach is employed to compute the posterior distribution of model parameters of high-fidelity models. As applications, parameters in a quasi-linear second-order elliptic boundary-value problem (BVP) are calibrated using a second-order linear elliptic BVP. In a second application, parameters of a tumor growth model involving nonlinear time-dependent PDEs are calibrated using a lower fidelity linear tumor growth model with known parameter values.

A-posterior estimates↗
DOE OSTI2012

A BAYESIAN MONTE CARLO ANALYSIS OF THE M-{sigma} RELATION

We present an analysis of selection biases in the M{sub bh}-{sigma} relation using Monte Carlo simulations including the sphere of influence resolution selection bias and a selection bias in the velocity dispersion distribution. We find that the sphere of influence selection bias has a significant effect on the measured slope of the M{sub bh}-{sigma} relation, modeled as {beta}{sub intrinsic} = -4.69 + 2.22{beta}{sub measured}, where the measured slope is shallower than the model slope in the parameter range of {beta} > 4, with larger corrections for steeper model slopes. Therefore, when the sphere of influence is used as a criterion to exclude unreliable measurements, it also introduces a selection bias that needs to be modeled to restore the intrinsic slope of the relation. We find that the selection effect due to the velocity dispersion distribution of the sample, which might not follow the overall distribution of the population, is not important for slopes of {beta} {approx} 4-6 of a logarithmically linear M{sub bh}-{sigma} relation, which could impact some studies that measure low (e.g., {beta} < 4) slopes. Combining the selection biases in velocity dispersions and the sphere of influence cut, we find that the uncertainty of the slope is larger than the value without modeling these effects and estimate an intrinsic slope of {beta} = 5.28{sup +0.84}{sub -0.55}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

UNSUPERVISED TRANSIENT LIGHT CURVE ANALYSIS VIA HIERARCHICAL BAYESIAN INFERENCE

Historically, light curve studies of supernovae (SNe) and other transient classes have focused on individual objects with copious and high signal-to-noise observations. In the nascent era of wide field transient searches, objects with detailed observations are decreasing as a fraction of the overall known SN population, and this strategy sacrifices the majority of the information contained in the data about the underlying population of transients. A population level modeling approach, simultaneously fitting all available observations of objects in a transient sub-class of interest, fully mines the data to infer the properties of the population and avoids certain systematic biases. We present a novel hierarchical Bayesian statistical model for population level modeling of transient light curves, and discuss its implementation using an efficient Hamiltonian Monte Carlo technique. As a test case, we apply this model to the Type IIP SN sample from the Pan-STARRS1 Medium Deep Survey, consisting of 18,837 photometric observations of 76 SNe, corresponding to a joint posterior distribution with 9176 parameters under our model. Our hierarchical model fits provide improved constraints on light curve parameters relevant to the physical properties of their progenitor stars relative to modeling individual light curves alone. Moreover, we directly evaluate the probability for occurrence rates of unseen light curve characteristics from the model hyperparameters, addressing observational biases in survey methodology. We view this modeling framework as an unsupervised machine learning technique with the ability to maximize scientific returns from data to be collected by future wide field transient searches like LSST.

79 ASTRONOMY AND ASTROPHYSICS↗