Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed inference”

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,369 records · Page 76

Comparison of Precision of Biomass Estimates in Regional Field Sample Surveys and Airborne LiDAR-Assisted Surveys in Hedmark County, Norway

Airborne scanning LiDAR (Light Detection and Ranging) has emerged as a promising tool to provide auxiliary data for sample surveys aiming at estimation of above-ground tree biomass (AGB), with potential applications in REDD forest monitoring. For larger geographical regions such as counties, states or nations, it is not feasible to collect airborne LiDAR data continuously ("wall-to-wall") over the entire area of interest. Two-stage cluster survey designs have therefore been demonstrated by which LiDAR data are collected along selected individual flight-lines treated as clusters and with ground plots sampled along these LiDAR swaths. Recently, analytical AGB estimators and associated variance estimators that quantify the sampling variability have been proposed. Empirical studies employing these estimators have shown a seemingly equal or even larger uncertainty of the AGB estimates obtained with extensive use of LiDAR data to support the estimation as compared to pure field-based estimates employing estimators appropriate under simple random sampling (SRS). However, comparison of uncertainty estimates under SRS and sophisticated two-stage designs is complicated by large differences in the designs and assumptions. In this study, probability-based principles to estimation and inference were followed. We assumed designs of a field sample and a LiDAR-assisted survey of Hedmark County (HC) (27,390 km2), Norway, considered to be more comparable than those assumed in previous studies. The field sample consisted of 659 systematically distributed National Forest Inventory (NFI) plots and the airborne scanning LiDAR data were collected along 53 parallel flight-lines flown over the NFI plots. We compared AGB estimates based on the field survey only assuming SRS against corresponding estimates assuming two-phase (double) sampling with LiDAR and employing model-assisted estimators. We also compared AGB estimates based on the field survey only assuming two-stage sampling (the NFI plots being grouped in clusters) against corresponding estimates assuming two-stage sampling with the LiDAR and employing model-assisted estimators. For each of the two comparisons, the standard errors of the AGB estimates were consistently lower for the LiDAR-assisted designs. The overall reduction of the standard errors in the LiDAR-assisted estimation was around 40-60% compared to the pure field survey. We conclude that the previously proposed two-stage model-assisted estimators are inappropriate for surveys with unequal lengths of the LiDAR flight-lines and new estimators are needed. Some options for design of LiDAR-assisted sample surveys under REDD are also discussed, which capitalize on the flexibility offered when the field survey is designed as an integrated part of the overall survey design as opposed to previous LiDAR-assisted sample surveys in the boreal and temperate zones which have been restricted by the current design of an existing NFI.

Comparison↗

MODIS Aerosol Observations used to Constrain Dust Distributions and Lifecycle in the NASA GEOS-5 Model

Approximately 240 Tg of mineral dust aerosol are transported annually from Saharan Africa to the Atlantic Ocean. Dust affects the Earth radiation budget, and plays direct (through scattering and absorption of radiation) and indirect (through modification of cloud properties and environment) roles in climate. Deposition of dust to the surface provides an important nutrient source to terrestrial and oceanic ecosystems. Dust is additionally a contributor to adverse air quality. Among the tools toward understanding the lifecycle and impacts of mineral dust aerosols are numerical models. Important constraints on these models come from quantitative satellite observations, like those from the space-based Moderate Resolution Imaging Spectroradiometer (MODIS). In particular, Kauhan et al. [2005] used MODIS aerosol observations to infer transport and deposition fluxes of Saharan dust over the Atlantic, Caribbean, and Amazonian basins. Those observations are used here to constrain the transport of dust and its interannual variability simulated in the NASA GEOS-5 general circulation model and data assimilation system. Significant uncertainty exists in the MODIS-derived fluxes, however, due to uncertainty in the wind fields provided by meteorological analyses in this region. That same uncertainty in the wind fields is manifest in our GEOS-5 simulations of dust distributions. Here we use MODIS observations to investigate the seasonality and location of the Saharan dust plume and explore through sensitivity analysis of our model the meteorological controls on the dust distribution, including dust direct radiative effects and sub-gridscale source and sink processes.

Colarco, P.↗

C 13 ( n , 2 n γ ) C 12 γ -ray production in the 14–16 MeV incident neutron energy range

We report γ-ray emission from 12 C and 13 C samples irradiated with deuterium-tritium fusion neutrons was experimentally measured at the Omega Laser Facility and at the Ohio University Edwards Accelerator Laboratory. The intent of these measurements was to determine the feasibility of using 13 C-based plastic ablators with embedded 12 C layers for “dark mix” diagnosis of inertial confinement fusion implosions. Spectrally resolved measurements at Ohio University identified significant 4.44-MeV γ-ray emission from the 13 C(n, 2nγ) 12 C-L1 reaction channel. The recorded 4.44-MeV 13 C signal was compared against emission from an identically irradiated 12 C target with known 12 C(n, n'γ) 12 C-L1 cross section, which resulted in an average 13 C(n, 2nγ) 12 C-L1 cross section of 117 ± 17 mb over the incident neutron energy distribution range from 14.4 to 15.8 MeV. Integrated 13C γ -ray signals above 2.9 MeV recorded with the Gas Cherenkov Detector at Omega exceeded MCNP6.1 predictions by a factor of 3. The additional signal was attributed to 4.44-MeV γs resulting in an inferred 13 C(n, 2nγ) 12 C-L1 cross section of 95 ± 11 mb at 14.1-MeV average incident neutron energy. As a result, the 13 C-based “dark mix” diagnostic concept was deemed infeasible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Disentangling Beam Losses in The Fermilab Main Injector Enclosure Using Real-Time Edge AI

The Fermilab Main Injector enclosure houses two accelerators, the Main Injector and Recycler Ring. During normal operation, high intensity proton beams exist simultaneously in both. The two accelerators share the same beam loss monitors (BLM) and monitoring system. Deciphering the origin of any of the 260 BLM readings is often difficult. The (Accelerator) Real-time Edge AI for Distributed Systems project, or READS, has developed an AI/ML model, and implemented it on fast FPGA hardware, that disentangles mixed beam losses and attributes probabilities to each BLM as to which machine(s) the loss originated from in real-time. The model inferences are then streamed to the Fermilab accelerator controls network (ACNET) where they are available for operators and experts alike to aid in tuning the machines.

43 PARTICLE ACCELERATORS↗

Comparing Automated Posterior Estimation Techniques for Modeling Strong Lenses In Ground-based Survey Data

Current and future ground-based cosmological surveys, such as the Dark Energy Survey (DES), and the Vera Rubin Observatory Legacy Survey of Space and Time (LSST), are predicted to discover thousands to tens of thousands of strong gravitational lenses. The large number of strong lenses discoverable in future surveys will make strong lensing a highly competitive and complementary cosmic probe. However, conventional lens modeling techniques are unable to scale up to the sheer number of lenses that will be discovered through upcoming surveys. Therefore, the use of automated lens analysis techniques is necessary. We demonstrate that machine learning methods can be used to automate the inference of informative model posteriors of strong lensing systems in ground-based surveys with credible uncertainty estimation. We present two Simulation-Based Inference (SBI) approaches for lens parameter estimation of galaxy-galaxy lenses. We demonstrate applications of Neural Posteriors Estima tors (NPEs) and Bayesian Neural Network (BNNs) to automate the inference of a 12-parameter lensing system for DES-like ground-based imaging data. We apply a suite of diagnostics (e.g., posterior coverage and SBC) to validate the performance of our methods. We find that NPEs outperform the BNN, producing posterior distributions that are for the most part both more accurate and more precise; in particular, several source-light model parameters are systematically biased in the BNN implementation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hydrostatic equilibrium of gas, extent of cosmic ray confinement, and radio emission in the Galaxy

A detailed study of the hydrostatic equilibrium of the gas-field system perpendicular to the galactic plane in the neighborhood of the sun, assuming that the magnetic field is parallel to the plane, has been carried out. It is found that the distribution of pressure due to cosmic rays and the magnetic field is different from that due to the gas. The observed radio spectrum in the direction of the galactic pole can be understood if one includes in the hydrostatic equilibrium the distribution of gas in the halo. The results of this analysis show that in the galactic plane, the cosmic-ray density is approximately 2 to 3 eV/cu cm with a magnetic-field strength of about 5 to 6 microgauss. It is inferred that cosmic rays are confined to a more extended volume in the Galaxy (with an equivalent thickness of at least 2 kpc in the solar neighborhood) than the gas disk, and the magnetic-field strength slowly decreases to a value of 1 microgauss only beyond about 5 kpc from the plane. Further, there is no evidence for a quasi-spherical radio halo, although the radio emission extends beyond the gas disk with scale height of at least 1 kpc. A lower limit of about 6 million years for the confinement of cosmic rays in the Galaxy has been set.

Badhwar, G. D.↗

The Swift Burst Alert Telescope Detected Seyfert 1 Galaxies: X-Ray Broadband Properties and Warm Absorbers

We present results from an analysis of the broadband, 0.3-195 keV, X-ray spectra of 48 Seyfert 1-1.5 sources detected in the very hard X-rays with the Swift Burst Alert Telescope (BAT). This sample is selected in an all-sky survey conducted in the 14-195 keV band. Therefore, our sources are largely unbiased toward both obscuration and host galaxy properties. Our detailed and uniform model fits to Suzaku/BAT and XMM-Newton/BAT spectra include the neutral absorption, direct power-law, reflected emission, soft excess, warm absorption, and narrow Fe I K[alpha] emission properties for the entire sample. We significantly detect O VII and O VIII edges in 52% of our sample. The strength of these detections is strongly correlated with the neutral column density measured in the spectrum. Among the strongest detections, X-ray grating and UV observations, where available, indicate outflowing material. The ionized column densities of sources with O VII and O VIII detections are clustered in a narrow range with Nwarm [approx] 1021 cm-2, while sources without strong detections have column densities of ionized gas an order of magnitude lower. Therefore, we note that sources without strong detections likely have warm ionized outflows present but at low column densities that are not easily probed with current X-ray observations. Sources with strong complex absorption have a strong soft excess, which may or may not be due to difficulties in modeling the complex spectra of these sources. Still, the detection of a flat [Gamma] [approx] 1 and a strong soft excess may allow us to infer the presence of strong absorption in low signal-to-noise active galactic nucleus spectra. Additionally, we include a useful correction from the Swift BAT luminosity to bolometric luminosity, based on a comparison of our spectral fitting results with published spectral energy distribution fits from 33 of our sources.

Winter, Lisa M.↗

Alleviating prior dependencies for DESI DR1 clustering fits through reparameterization

Bayesian analyses of the full-shape clustering of Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) exhibit prior-volume projection effects, whereby weakly constrained nuisance parameters of the Effective Field Theory of Large Scale Structure (EFTofLSS) shift marginalized cosmological posteriors away from the posterior maximum. We reanalyze DESI DR1 power spectrum multipoles using two complementary mitigation strategies: (i) nonlinear orthogonalization to decorrelate nuisance and cosmological parameter priors, and (ii) a fully reparameterization-invariant Jeffreys prior over all EFTofLSS coefficients, evaluated on-the-fly via closed-form Jacobians. Including data from DESI, Big-Bang Nuclesynthesis and a constraint on $n_{\mathrm{s}}$, baseline priors lead to multi-$σ$ projection in the Hubble parameter $H_{0}$ and dark energy equation of state parameters $w_{0}$ and $w_{a}$; the Jeffreys prior successfully recenters these posteriors to enclose the maximum a posteriori estimate within the 68% credible regions, demonstrating clear mitigation of projection effects for these late-time expansion parameters. A hybrid Jeffreys+baseline-Gaussian configuration controls residual over-broad tails in the physical cold dark matter density $ω_{\mathrm{c}}$ while preserving the volume correction, and is our favoured approach. We compare the credible intervals derived using our methodology to those obtained using Halo Occupation Distribution (HOD)-informed priors and to confidence intervals derived using frequentist profile likelihood analyses, finding agreement in both central values and degeneracy directions in the $w_{0}$--$w_{a}$ plane. This demonstrates that, once projection effects are properly controlled, we can make robust inferences about the late-time cosmological expansion independent of the statistical framework adopted.

Bonici, M. [Waterloo U.; Perimeter Inst. Theor. Ph↗

Coordinated neutron and X-ray computed tomography of meteorites: Detection and distribution of hydrogen-bearing materials

The presence and distribution of hydrogen-bearing materials in meteorites are important constraints on processes in the early solar system, and the delivery of volatile constituents to growing planets. Here, in this paper, we show that coordinated neutron and X-ray computed tomography, NXCT, can reveal the presence and distributions of hydrogen-bearing materials in meteorites, and thus help constrain the presence and actions of water in the early solar system. NXCT is nearly nondestructive of meteorite samples. Neutron fluence in NXCT is approximately seven orders of magnitude less than in typical instrumental neutron activation analysis, and so produces little residual radioactivity and currently undetectable changes in isotope ratios. Heating during NXCT is minimal, but NXCT will overprint the record of cosmic ray exposure held in natural thermoluminescence. Two meteorites were examined. EET 87503 is a howardite, a regolith breccia inferred to be from the asteroid 4 Vesta, and contains fragments of eucrite basalt, diogenite pyroxenite, and H-rich carbonaceous chondrites. With NXCT, the chondrite fragments within the meteorite piece can be clearly located and characterized, in preparation for possible extraction and detailed analyses. Graves Nunataks (GRA) 06100 is a CR2 chondrite meteorite that contains abundant iron metal and H-bearing silicates from aqueous alteration. In NXCT, H-bearing altered material is clearly distinguished from metal, and its distribution in three dimensions is revealed as a constraint on the processes of alteration.

79 ASTRONOMY AND ASTROPHYSICS↗

Joint estimation of biogeochemical model parameters from multiple experiments: A bayesian approach applied to mercury methylation

Here, to characterize complex biogeochemical systems, results from multiple experiments, where each targets a specific subprocess, are commonly combined. The resulting datasets are interpreted through the calibration of biogeochemical models for process inference and predictions. Commonly used calibration approaches of fitting datasets from individual experiments to subprocess models one at a time is prone to missing information shared between datasets and incomplete uncertainty propagation. We propose a Bayesian joint-fitting scheme addressing the above-mentioned concerns by jointly fitting all the available datasets, thus calibrating the entire biogeochemical model in one go using Markov Chain Monte Carlo (MCMC). The identification of null spaces in the parameter distributions from MCMC guided the simplification of certain subprocess models. For example, fast kinetic sorption was replaced by equilibrium sorption, and Monod demethylation was replaced by first-order demethylation. Joint fitting of datasets resulted in complete uncertainty propagation with parameter estimates informed by all available data.

54 ENVIRONMENTAL SCIENCES↗

Analytic solutions for Asay foil trajectories with implications for ejecta source models and mass measurements

We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An investigation of several NACA 1-series inlets at Mach numbers from 0.4 to 1.29 for mass flow ratios near 1.0

An investigation to determine the performance of eight NACA 1-series inlets at massflow ratios near 1.0 was conducted in the Langley 16-foot transonic tunnel. The inlet diameter ratios (ratio of inlet diameter to maximum diameter) were 0.85 and 0.89 for an inlet length ratio (ratio of inlet length to maximum diameter) of 1.0. Inlet lip radius varied from 0.061 cm to 0.251 cm, and internal contraction area ratio (ratio of inlet area to throat area) varied from 1.006 to 1.201. Reynolds number based on model maximum diameter ranged from 3,600,000 at a Mach number of 400,000 to 5,900,000 at a Mach number of 1.29. The results indicate that nearly uniform pressure distributions on a given inlet were obtained over a limited range of mass-flow ratios and Mach numbers. When inlet lip thickness was increased by means of lip radius or contraction ratio, the inlet critical Mach number decreased. Drag-divergence Mach number inferred from forebody pressure integrations was above 0.94 for most of the inlets tested.

Re, R. J.↗

Patterns of seismic activity preceding large earthquakes

A mechanical model of seismic faults is employed to investigate the seismic activities that occur prior to major events. The block-and-spring model dynamically generates a statistical distribution of smaller slipping events that precede large events, and the results satisfy the Gutenberg-Richter law. The scaling behavior during a loading cycle suggests small but systematic variations in space and time with maximum activity acceleration near the future epicenter. Activity patterns inferred from data on seismicity in California demonstrate a regional aspect; increased activity in certain areas are found to precede major earthquake events. One example is given regarding the Loma Prieta earthquake of 1989 which is located near a fault section associated with increased activity levels.

Shaw, Bruce E.↗

Connection between ambient density fluctuations and clumpy Langmuir waves in type III radio sources

A recent stochastic-growth theory of clumpy Langmuir waves in type III sources is shown to imply that the clumps will have the same size distribution as the ambient low-frequency density fluctuations in the solar wind. Spectral analysis of Langmuir-wave time series from the ISEE 3 plasma wave instrument confirms this prediction to within the uncertainties in the spectra. The smallest Langmuir clump size is inferred to be in the range 0.4-30 km in general, and 2-30 km for beam-resonant waves, and it is concluded that the diffusion of waves in the source is anomalous.

Robinson, P. A.↗

AVIRIS calibration and application in coastal oceanic environments - Tracers of soluble and particulate constituents of the Tampa Bay coastal plume

AVIRIS is a testbed for future spacecraft sensors (such as HIRIS and MODIS) planned for the Earth Observing System. Model-derived absorption coefficients at 415 nm, a(415), and back-scattering coefficients at 671 nm, b sub b (671) for Tampa Bay waters were used to create images from AVIRIS data of the dissolved component of a(415) due to gelbstoff, a sub g (415), and salinity. Images of a sub g (415), salinity, and b sub b (671) were used to depict the distribution of dissolved and particulate constituents, respectively, for Tampa Bay plume during late, ebb-tidal conditions. Salinity covaried with a sub g (415), which provided a means of mapping salinity from the a sub g (415) imagery. The concentration of suspended particles, as inferred from b sub g (671), was extremely variable in the shallow regions where waves and currents interacted. Pollutants covarying with fresh water or suspended sediments can be mapped from a sub g (415) and b sub b (671) images, respectively.

Carder, K. L.↗

Tidal Disruptions of Main Sequence Stars - I. Observable Quantities and their Dependence on Stellar and Black Hole Mass

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses(MM0.1510)and six black hole masses( MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me;5×106 this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.

Taeho Ryu↗

Intensive disc-reverberation mapping of Fairall 9: first year of Swift and LCO monitoring

We present results of time-series analysis of the first year of the Fairall 9 intensive disc-reverberation campaign. We used Swift and the Las Cumbres Observatory global telescope network to continuously monitor Fairall 9 from X-rays to near-infrared at a daily to subdaily cadence. The cross-correlation function between bands provides evidence for a lag spectrum consistent with the τ ∝ λ^(4/3) scaling expected for an optically thick, geometrically thin blackbody accretion disc. Decomposing the flux into constant and variable components, the variable component’s spectral energy distribution is slightly steeper than the standard accretion disc prediction. We find evidence at the Balmer edge in both the lag and flux spectra for an additional bound-free continuum contribution that may arise from reprocessing in the broad-line region. The inferred driving light curve suggests two distinct components, a rapidly variable (<4 d) component arising from X-ray reprocessing, and a more slowly varying (>100 d) component with an opposite lag to the reverberation signal.

J V Hernández Santisteban↗

Parameter inference with non-linear galaxy clustering: accounting for theoretical uncertainties

ABSTRACT We implement euclidemulator (version 1), an emulator for the non-linear correction of the matter power spectrum, into the Markov chain Monte Carlo forecasting code montepython. We compare the performance of halofit, hmcode, and euclidemulator1, both at the level of power spectrum prediction and at the level of posterior probability distributions of the cosmological parameters, for different cosmological models and different galaxy power spectrum wavenumber cut-offs. We confirm that the choice of the power spectrum predictor has a non-negligible effect on the computed sensitivities when doing cosmological parameter forecasting, even for a conservative wavenumber cut-off of $0.2\, h\, {\rm Mpc}^{-1}$. We find that euclidemulator1 is on average up to 17 per cent more sensitive to the cosmological parameters than the other two codes, with the most significant improvements being for the Hubble parameter of up to 42 per cent and the equation of state of dark energy of up to 26 per cent, depending on the case. In addition, we point out that the choice of the power spectrum predictor contributes to the risk of computing a significantly biased mean cosmology when doing parameter estimations. For the four tested scenarios we find biases, averaged over the cosmological parameters, of between 0.5σ and 2σ (from below 1σ up to 6σ for individual parameters). This paper provides a proof of concept that this risk can be mitigated by taking a well-tailored theoretical uncertainty into account as this allows to reduce the bias by a factor of 2 to 5, depending on the case under consideration, while keeping posterior credibility contours small: the standard deviations are amplified by a factor of ≤1.4 in all cases.

Astronomy & Astrophysics↗