Engineering PapersSearch

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 55 records · Page 3

OGLE-2013-BLG-1761Lb: A Massive Planet around an MK Dwarf

We report the discovery and the analysis of the planetary microlensing event, OGLE-2013-BLG-1761. There are some degenerate solutions in this event because the planetary anomaly is only sparsely sampled. However, the detailed light curve analysis ruled out all stellar binary models and shows the lens to be a planetary system. There is the so-called close wide degeneracy in the solutions with the planet host mass ratio of q approx.(7.0+/-2.0) x 10(exp -3) and q approx.(8.1+/-2.6) x 10(exp -3) with the projected separation in Einstein radius units of s = 0.95 (close) and s = 1.18(wide), respectively. The microlens parallax effect is not detected, but the finite source effect is detected. Our Bayesian analysis indicates that the lens system is located -D(sub L) = 6.9(+ 1.0 -1.2)kpc away from us and the host star is an M/K dwarf with amass of M(sub L) = 0.33(+ 0.32- 1.9)Stellar Mass orbited by a super-Jupiter mass planet with a mass of m(sub p) = 2.7(+ 2.5 - 1.5) M(sub Jup) at the projected separation of a(sub l) = 1.8(+ 0.5 -0.5)au. The preference of the large lens distance in the Bayesian analysis is due to the relatively large observed source star radius. The distance and other physical parameters may be constrained by the future high-resolution imaging by large ground telescopes or HST. If the estimated lens distance is correct, then this planet provides another sample for testing the claimed deficit of planets in the Galactic bulge.

Hirao, Y.

Variability Timescale and Spectral Index of Sgr A* in the Near Infrared: Approximate Bayesian Computation Analysis of the Variability of the Closest Supermassive Black Hole

Sagittarius A* (Sgr A* ) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability data set of Sgr A* to date: eight 24 hr epochs of continuous monitoring of Sgr A* at 4.5 μm with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 μm data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 μm data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current fast Fourier transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is b 243 57 82 t = - + minutes (90% credible interval). The PSD has no detectable features on timescales down to 8.5 minutes (95% credible level), which is the ISCO orbital frequency for a dimensionless spin parameter a = 0.92. One light curve measured simultaneously at 2.12 and 4.5 μm during a low flux-density phase gave a spectral index αs = 1.6 ± 0.1 (F µ n s n -a ). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined data sets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a nonthermal electron population from below 20 GHz through the NIR.

accretion

Bayesian Statistical Analysis for Mass Spectrometric Data Processing

• Thermal ionization mass spectrometry (TIMS) is a “gold standard’ technique for actinide isotope amount ratio and assay measurements. • Ubiquitously used in: • Nuclear Nonproliferation • Nuclear Safeguards • Nuclear Forensics • Basic Science • Savannah River National Laboratory (SRNL) installed a new Thermo Scientific TRITON Plus TIMS in early 2024.

McLarty, Ellis C. [Savannah River National Laborat

Bayesian Statistical Analysis for Mass Spectrometric Data Processing

Thermal ionization mass spectrometry (TIMS) is a “gold standard’ technique for actinide isotope amount ratio and assay measurements. • Ubiquitously used in: • Nuclear Nonproliferation • Nuclear Safeguards • Nuclear Forensics • Basic Science • Savannah River National Laboratory (SRNL) installed a new Thermo Scientific TRITON Plus TIMS in early 2024.

McLarty, Ellis C. [Savannah River National Laborat

A Survey of Eight Hot Jupiters in Secondary Eclipse Using WIRCam at CFHT

We present near-infrared high-precision photometry for eight transiting hot Jupiters observed during their predicted secondary eclipses. Our observations were carried out using the staring mode of the WIRCam instrument on the Canada-France-Hawaii Telescope (CFHT).We present the observing strategies and data reduction methods which delivered time series photometry with statistical photometric precision as low as 0.11 percent. We performed a Bayesian analysis to model the eclipse parameters and systematics simultaneously. The measured planet-to-star flux ratios allowed us to constrain the thermal emission from the day side of these hot Jupiters, as we derived the planet brightness temperatures. Our results combined with previously observed eclipses reveal an excess in the brightness temperatures relative to the blackbody prediction for the equilibrium temperatures of the planets for a wide range of heat redistribution factors. We find a trend that this excess appears to be larger for planets with lower equilibrium temperatures. This may imply some additional sources of radiation, such as reflected light from the host star and/or thermal emission from residual internal heat from the formation of the planet.

Bayesian Analysis

A Bayesian modification to the Jelinski-Moranda software reliability growth model

The Jelinski-Moranda (JM) model for software reliability was examined. It is suggested that a major reason for the poor results given by this model is the poor performance of the maximum likelihood method (ML) of parameter estimation. A reparameterization and Bayesian analysis, involving a slight modelling change, are proposed. It is shown that this new Bayesian-Jelinski-Moranda model (BJM) is mathematically quite tractable, and several metrics of interest to practitioners are obtained. The BJM and JM models are compared by using several sets of real software failure data collected and in all cases the BJM model gives superior reliability predictions. A change in the assumption which underlay both models to present the debugging process more accurately is discussed.

Littlewood, B.

New approaches to Bayesian uncertainty quantification for Nuclear Science (Final Technical Report)

Inverse problems play a central role in experimentation and theory/data comparisons for many areas of modern Nuclear Physics (NP) and High-Energy Physics (HEP). Bayes’s Theorem is a powerful tool for solving Inverse Problems, providing conceptually transparent and unbiased constraints on theoretical parameters and their uncertainties (“Bayesian Inference”) and enabling the quantification of agreement or tension between models and data. However, analyses based on Bayesian Inference are often challenging for NP and HEP applications, either because of the large number of parameters in the problem, the high computational cost, or both. We propose a multi-institutional collaboration to develop and deploy novel Bayesian analysis tools that advance the scientific scope of a broad range of current and future NP experiments. This project brings together NP domain scientists working on several high-profile NP projects for which new, high-performance Bayesian Uncertainty Quantification (“Bayesian UQ”) methods are essential to carry out the science, and data scientists who are developing state-of-the-art methods applicable to these problems. The NP projects in this proposal comprise measurements of the mass and fundamental nature of the neutrino; study of the Quark-Gluon Plasma that filled the early universe; and mapping of natural and anthropogenic radiation environments. While these NP projects have very different scientific goals, with datasets and analysis approaches that differ significantly, they share common requirements for improving computationally intensive Bayesian analyses using advanced Machine Learning algorithms and will benefit strongly from a coherent effort to develop general solutions. This proposal brings together these projects and forefront ML-based data science algorithms to develop such general solutions. The methods developed in this project will also be more widely applicable, thereby advancing science in the larger Nuclear Physics portfolio.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Bayesian Linear Regression for Hugoniot Data

This repository provides the code and datasets used in the paper Bayesian Analysis of Linear Shock Compression Data. This paper analyzes publicly available shock compression datasets on copper, argon, and nickel from Marsh (1980) using Bayesian linear regression, and compares the results with those obtained using bootstrapping methods. References: - Marsh, S. P. (1980). LASL shock Hugoniot data (Vol. 5). Univ of California Press.

Bernstein, JasonA [Lawrence Livermore National Lab

The Challenges of Credible Thermal Protection System Reliability Quantification

The paper discusses several of the challenges associated with developing a credible reliability estimate for a human-rated crew capsule thermal protection system. The process of developing such a credible estimate is subject to the quantification, modeling and propagation of numerous uncertainties within a probabilistic analysis. The development of specific investment recommendations, to improve the reliability prediction, among various potential testing and programmatic options is then accomplished through Bayesian analysis.

Green, Lawrence L.

A Second Look at 12 Candidate Dual AGNs Using BAYMAX

We present an analysis of 12 optically selected dual active galactic nucleus (AGN) candidates at z < 0.34. Each candidate was originally identified via double-peaked [O III] λ5007 emission lines and received follow-up Chandra and Hubble Space Telescope (HST) observations. Because the X-ray data are low-count (<100 counts) with small separations (<1''), a robust analysis is necessary for classifying each source. Pairing long-slit [O III] observations with existing Chandra observations, we re-analyze the X-ray observations with Bayesian AnalYsis of AGNs in X-rays to determine whether the X-ray emission from each system is more likely a single or dual point source. We find that 4 of the 12 sources are likely dual X-ray point-source systems. We examine each point source's spectra via a Monte Carlo method that probabilistically identifies the likely origin of each photon. When doing so, we find that (i) the secondary X-ray point sources in two of the systems have L(X)< 10(exp 40) erg/s, such that we cannot rule out a non-AGN origin, (ii) one source has a secondary with L(X)> 10(exp 40) erg/s but a spectrum that is too soft to definitively preclude being X-ray emitting diffuse gas that was photoionized by the primary AGN, and (iii) one system (SDSS J1126+2944) is a dual AGN. Additionally, using complementary HST observations, we analyze a subsample of systems that are visually identified as merging. Our results suggest that dual AGNs may preferentially reside in mergers with small separations, consistent with both simulations and observations.

Adi Foord

An automated land-use mapping comparison of the Bayesian maximum likelihood and linear discriminant analysis algorithms

The Bayesian maximum likelihood parametric classifier has been tested against the data-based formulation designated 'linear discrimination analysis', using the 'GLIKE' decision and "CLASSIFY' classification algorithms in the Landsat Mapping System. Identical supervised training sets, USGS land use/land cover classes, and various combinations of Landsat image and ancilliary geodata variables, were used to compare the algorithms' thematic mapping accuracy on a single-date summer subscene, with a cellularized USGS land use map of the same time frame furnishing the ground truth reference. CLASSIFY, which accepts a priori class probabilities, is found to be more accurate than GLIKE, which assumes equal class occurrences, for all three mapping variable sets and both levels of detail. These results may be generalized to direct accuracy, time, cost, and flexibility advantages of linear discriminant analysis over Bayesian methods.

Tom, C. H.

Visible core spectroscopy at Wendelstein 7-X

This paper presents an overview of recent hardware extensions and data analysis developments to the Wendelstein 7-X visible core spectroscopy systems. These include upgrades to prepare the in-vessel components for long-pulse operation, nine additional spectrometers, a new line of sight array for passive spectroscopy, and a coherence imaging charge exchange spectroscopy diagnostic. Progress in data analysis includes ion temperatures and densities from multiple impurity species, a statistical comparison with x-ray crystal spectrometer measurements, neutral density measurements from thermal passive Balmer-alpha emission, and a Bayesian analysis of active hydrogen emission, which is able to infer electron density and main ion temperature profiles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND

16 O ⁡(𝑒,𝑒′⁢𝛼)⁢ 12 C measurements and the 12 C ⁡(𝛼,𝛾)⁢ 16 O astrophysical reaction rate

The 12 C ⁡(𝛼,𝛾)⁢ 16 O reaction, an important component of stellar helium burning, has a key role in nuclear astrophysics. It has significant impact on the evolution and final state of heavy to low mass stars, maximum mass of stellar formed black holes and also shapes the elemental abundances resulting from nucleosynthesis in such stars. Providing a reliable estimate for the energy dependence of this reaction at stellar helium burning temperatures has been a longstanding and important goal. In this work, we study the role of potential new 𝐸⁢1 and 𝐸⁢2 measurements of the 16 O ⁡(𝑒,𝑒′⁢𝛼)⁢ 12 C reaction in reducing the overall uncertainty in the astrophysical 𝑆 factor for the 12 C ⁡(𝛼,𝛾)⁢ 16 O 𝐸⁢1 and 𝐸⁢2 ground state capture extrapolated to a stellar energy of 300 keV. A multilevel 𝑅-matrix analysis is used to make extrapolations of the 𝑆 𝐸⁢1 (300 keV) and 𝑆 𝐸⁢2 (300 keV) factors for the 12 C ⁡(𝛼,𝛾)⁢ 16 O reaction from existing ground state capture data. Bayesian analysis is used to quantify the uncertainties in the extrapolations for both the existing data alone and also when possible new experimental data are included. In particular, we consider a new experiment that would make use of a high-intensity low-energy electron beam that impinges on a windowless oxygen gas target as a means to determine the total 𝐸⁢1 and 𝐸⁢2 ground state cross sections for this reaction. Here we find that the new data could significantly reduce the 𝑆⁡(300) uncertainties. Splitting the new data into high and low energy regions shows that both low and high energy data are effective in reducing the uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

A Preliminary Detection of Arcminute Scale Cosmic Microwave Background Anisotropy with the BIMA Array

We have used the Berkeley-Illinois-Maryland-Association (BIMA) array outfitted with sensitive cm-wave receivers to expand our search for minute scale anisotropy of the Cosmic Microwave Background (CMB). The interferometer was placed in a compact configuration to obtain high brightness sensitivity on arcminute scales over its 6.6' FWHM field of view. The sensitivity of this experiment to flat band power peaks at a multipole of 1 = 5530 which corresponds to an angular scale of -2'. We present the analysis of a total of 470 hours of on-source integration time on eleven independent fields which were selected based on their low IR contrast and lack of bright radio sources. Applying a Bayesian analysis to the visibility data, we find CMB anisotropy flat band power Q_flat = 6.1(+2.8/-4.8) microKelvin at 68% confidence. The confidence of a nonzero signal is 76% and we find an upper limit of Q_flat < 12.4 microKelvin at 95% confidence. We have supplemented our BIMA observations with concurrent observations at 4.8 GHz with the VLA to search for and remove point sources. We find the point sources make an insignificant contribution to the observed anisotropy.

Dawson, K. S.

Bayesian estimation - A sensitivity analysis

The robustness of the assigned prior distribution in a Bayesian estimation problem is examined. A Bayesian analysis for a stochastic intensity parameter of a Poisson distribution is summarized in which the natural conjugate is assigned as the prior distribution of the random parameter. The sensitivity analysis is carried out by assuming the existence of a true prior which is different in form from that of the assigned prior distribution. By using mean-squared error as a measure of performance, the ensuing Bayes decision function is compared to the corresponding minimum variance unbiased estimator. Results indicate that the Bayes estimator is largely robust to deviations from the assigned prior and remains squared-error superior to the MVU type within a broad region.

Canavos, G. C.