Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Disentangling the Black Hole Mass Spectrum with Photometric Microlensing Surveys

Abstract From the formation mechanisms of stars and compact objects to nuclear physics, modern astronomy frequently leverages surveys to understand populations of objects to answer fundamental questions. The population of dark and isolated compact objects in the Galaxy contains critical information related to many of these topics, but is only practically accessible via gravitational microlensing. However, photometric microlensing observables are degenerate for different types of lenses, and one can seldom classify an event as involving either a compact object or stellar lens on its own. To address this difficulty, we apply a Bayesian framework that treats lens type probabilistically and jointly with a lens population model. This method allows lens population characteristics to be inferred despite intrinsic uncertainty in the lens class of any single event. We investigate this method’s effectiveness on a simulated ground-based photometric survey in the context of characterizing a hypothetical population of primordial black holes (PBHs) with an average mass of 30 M ⊙ . On simulated data, our method outperforms current black hole (BH) lens identification pipelines and characterizes different subpopulations of lenses while jointly constraining the PBH contribution to dark matter to ≈25%. Key to robust inference, our method can marginalize over population model uncertainty. We find the lower mass cutoff for stellar origin BHs, a key observable in understanding the BH mass gap, particularly difficult to infer in our simulations. This work lays the foundation for cutting-edge PBH abundance constraints to be extracted from current photometric microlensing surveys.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Diverse properties of interstellar medium embedding gamma-ray bursts at the epoch of reionization

Analysis is performed on ultra-high-resolution large-scale cosmological radiation-hydrodynamic simulations to quantify, for the first time, the physical environment of long-duration gamma-ray bursts (GRBs) at the epoch of reionization. We find that, on parsec scales, 13% of GRBs remain in high-density (≥10{sup 4} cm{sup –3}) low-temperature star-forming regions, whereas 87% of GRBs occur in low-density (∼10{sup –2.5} cm{sup –3}) high-temperature regions heated by supernovae. More importantly, the spectral properties of GRB afterglows, such as the neutral hydrogen column density, total hydrogen column density, dust column density, gas temperature, and metallicity of intervening absorbers, vary strongly from sight line to sight line. Although our model explains extant limited observationally inferred values with respect to circumburst density, metallicity, column density, and dust properties, a substantially larger sample of high-z GRB afterglows would be required to facilitate a statistically solid test of the model. Our findings indicate that any attempt to infer the physical properties (such as metallicity) of the interstellar medium (ISM) of the host galaxy based on a very small number (usually one) of sight lines would be precarious. Utilizing high-z GRBs to probe the ISM and intergalactic medium should be undertaken properly, taking into consideration the physical diversities of the ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Estimates of Active Region Area Coverage through Simultaneous Measurements of the He i λλ 5876 and 10830 Lines

Simultaneous, high-quality measurements of the neutral helium triplet features at 5876 Å and 10830 Å in a sample of solar-type stars are presented. The observations were made with ESO telescopes at the La Silla Paranal Observatory under program ID 088.D-0028(A) and MPG Utility Run for Fiber Extended-range Optical Spectrograph 088.A-9029(A). The equivalent widths of these features combined with chromospheric models are utilized to infer the fractional area coverage, or filling factor, of magnetic regions outside of spots. We find that the majority of the sample is characterized by filling factors less than unity. However, discrepancies occur among the coolest K-type and the warmest and most rapidly rotating F-type dwarf stars. We discuss these apparently anomalous results and find that in the case of K-type stars, they are an artifact of the application of chromospheric models best suited to the Sun than to stars with significantly lower T {sub eff}. The case of the F-type rapid rotators can be explained by the measurement uncertainties of the equivalent widths, but they may also be due to a non-magnetic heating component in their atmospheres. With the exceptions noted above, preliminary results suggest that the average heating rates in the active regions are the same from one star to the other, differing in the spatially integrated, observed level of activity due to the area coverage. Hence, differences in activity in this sample are mainly due to the filling factor of active regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Kangaroo’s First Hop: The Early Fast Cooling Phase of EP250108a/SN 2025kg

Fast X-ray transients are a rare and poorly understood population of events. Previously difficult to detect in real time, the launch of the Einstein Probe with its Wide-field X-ray Telescope has led to a rapid expansionof the sample and allowed the exploration of these transients across the electromagnetic spectrum. EP250108a is a recently detected example linked to an optical counterpart, SN 2025kg, or “the kangaroo.” Together with a companion Letter we present our observing campaign and analysis of this event. In this letter, we focus on the early evolution of the optical counterpart over the first 6 days, including our measurement of the redshift of z = 0.17641. We compare to other supernovae and fast transients showing similar features, finding significant similarities with SN 2006aj and SN 2020bvc, and show that the source is well modelled by a rapidly expanding cooling blackbody. We show the observed X-ray and radio properties are consistent with a collapsar-powered jet that is low energy (≲10 51 erg) and/or fails to break out of the dense material surrounding it. While we examine the possibility that the optical emission emerges from the shock produced as the supernova ejecta expand into a dense shell of circumstellar material, due to our X-ray and radio inferences, we favour a model where it arises from a shocked cocoon resulting from a trapped jet. This makes SN 2025 one of the few examples of this currently observationally rare event.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1978

Stratospheric NO and HNO3 observations in the Northern Hemisphere for three seasons

NO, HNO3, and O3 levels and air temperature were measured as a function of latitude in the 18 to 21 km region of the stratosphere, and the sum of odd nitrogen, equal to NO + NO2 + HNO3, was calculated and compared with model predictions. (NO2 values were inferred from photochemical equilibrium characteristics.) The data show that NO measurements generally exhibit good agreement with model predictions for low and midlatitudes but poor agreement at high latitudes. The experimental sum of odd nitrogen mixing ratios and model predictions agree within a factor of 2-1/2 or better at both 20 and 40 deg N, and show excellent agreement for latitudinal dependence.

Loewenstein, M.↗
DOE OSTI2014

Likelihood-informed dimension reduction for nonlinear inverse problems

The intrinsic dimensionality of an inverse problem is affected by prior information, the accuracy and number of observations, and the smoothing properties of the forward operator. From a Bayesian perspective, changes from the prior to the posterior may, in many problems, be confined to a relatively low-dimensional subspace of the parameter space. We present a dimension reduction approach that defines and identifies such a subspace, called the 'likelihood-informed subspace' (LIS), by characterizing the relative influences of the prior and the likelihood over the support of the posterior distribution. This identification enables new and more efficient computational methods for Bayesian inference with nonlinear forward models and Gaussian priors. In particular, we approximate the posterior distribution as the product of a lower-dimensional posterior defined on the LIS and the prior distribution marginalized onto the complementary subspace. Markov chain Monte Carlo sampling can then proceed in lower dimensions, with significant gains in computational efficiency. We also introduce a Rao–Blackwellization strategy that de-randomizes Monte Carlo estimates of posterior expectations for additional variance reduction. We demonstrate the efficiency of our methods using two numerical examples: inference of permeability in a groundwater system governed by an elliptic PDE, and an atmospheric remote sensing problem based on Global Ozone Monitoring System (GOMOS) observations.

42 ENGINEERING↗
DOE OSTI2026

First Estimation of Model Parameters for Neutrino-Induced Nucleon Knockout Using Simulation-Based Inference

To enable an accurate determination of oscillation parameters, accelerator-based neutrino experiments require detailed simulations of nuclear interaction physics in the GeV regime. While substantial effort from both theory and experiment is currently being invested to improve the fidelity of these simulations, their present deficiencies typically oblige experimental collaborations to resort to empirical tuning of simulation model parameters. As the precision requirements of the field continue to become more stringent, machine learning techniques may provide a powerful means of handling corresponding growth in the complexity of future neutrino interaction model tuning exercises. To study the suitability of simulation-based inference (SBI) for this physics application, in this paper we revisit a tuned configuration of the GENIE neutrino event generator that was originally developed by the MicroBooNE collaboration. Despite closely reproducing the adopted values of four physics parameters when confronted with the tuned cross-section predictions as input, we find that our trained SBI algorithm prefers modestly different values (within MicroBooNE's assigned uncertainties) and achieves slightly better goodness-of-fit when inference is run on the experimental data set originally used by MicroBooNE. We also find that our trained algorithm can create a fair approximation of an alternative neutrino scattering simulation, NuWro, that shares only a subset of its physics model parameters with GENIE.

Tame-Narvaez, Karla [Fermilab] (ORCID:000000022249↗
DOE OSTI2014

WE-A-17A-11: Implanted Brachytherapy Seed Movement Due to Transrectal Ultrasound Probe-Induced Prostate Deformation

Purpose: To characterize the movement of implanted brachytherapy seeds due to transrectal ultrasound probe-induced prostate deformation and to estimate the effects on prostate dosimetry. Methods: Implanted probe-in and probe-removed seed distributions were reconstructed for 10 patients using C-arm fluoroscopy imaging. The prostate was delineated on ultrasound and registered to the fluoroscopy seeds using a visible subset of seeds and residual needle tracks. A linear tensor and shearing model correlated the seed movement with position. The seed movement model was used to infer the underlying prostate deformation and to simulate the prostate contour without probe compression. Changes in prostate and surrogate urethra dosimetry were calculated. Results: Seed movement patterns reflecting elastic decompression, lateral shearing, and rectal bending were observed. Elastic decompression was characterized by anterior-posterior expansion and superior-inferior and lateral contractions. For lateral shearing, anterior movement up to 6 mm was observed for extraprostatic seeds in the lateral peripheral region. The average intra-prostatic seed movement was 1.3 mm, and the residual after linear modeling was 0.6 mm. Prostate D90 increased by 4 Gy on average (8 Gy max) and was correlated with elastic decompression. For selected patients, lateral shearing resulted in differential change in D90 of 7 Gy between anterior and posterior quadrants, and increase in whole prostate D90 of 4 Gy. Urethra D10 increased by 4 Gy. Conclusion: Seed movement upon probe removal was characterized. The proposed model captured the linear correlation between seed movement and position. Whole prostate dose coverage increased slightly, due to the small but systematic seed movement associated with elastic decompression. Lateral shearing movement increased dose coverage in the anterior-lateral region, at the expense of the posterior-lateral region. The effect on whole prostate D90 was smaller due to the subset of peripheral seeds involved, but lateral shearing movement can have greater consequences for local dose coverage.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1983

A proposed origin for palimpsests and anomalous pit craters on Ganymede and Callisto

The hypothesis that palimpsests and anomalous pit craters are essentially pristine crater forms derived from high-velocity impacts and/or impacts into an ice crust with preimpact temperatures near melting is explored. The observational data are briefly reviewed, and an impact model is proposed for the direct formation of a palimpsest from an impact when the modification flow which produces the final crater is dominated by 'wet' fluid flow, as opposed to the 'dry' granular flow which produces normal craters. Conditions of 'wet' modification occur when the volume of impact melt remaining in the transient crater attains a volume comparable to the transient crater. The normal crater-palimpsest transition is found to occur for sufficiently large impacts or sufficiently fast impactors. The range of crater diameters and morphological characteristics inferred from the impact model is consistent with the observed characteristics of palimpsests and anomalous pit craters.

Croft, S. K.↗
NASA NTRS1979

Comparison of 2.8- and 21-cm microwave radiometer observations over soils with emission model calculations

An airborne experiment was conducted under NASA auspices to test the feasibility of detecting soil moisture by microwave remote sensing techniques over agricultural fields near Phoenix, Arizona at midday of April 5, 1974 and at dawn of the following day. Extensive ground data were obtained from 96 bare, sixteen hectare fields. Observations made using a scanning (2.8 cm) and a nonscanning (21 cm) radiometer were compared with the predictions of a radiative transfer emission model. It is shown that (1) the emitted intensity at both wavelengths correlates best with the near surface moisture, (2) surface roughness is found to more strongly affect the degree of polarization than the emitted intensity, (3) the slope of the intensity-moisture curves decreases in going from day to dawn, and (4) increased near surface moisture at dawn is characterized by increased polarization of emissions. The results of the experiment indicate that microwave techniques can be used to observe the history of the near surface moisture. The subsurface history must be inferred from soil physics models which use microwave results as boundary conditions.

Burke, W. J.↗
DOE OSTI2013

INTERACTING COSMIC RAYS WITH MOLECULAR CLOUDS: A BREMSSTRAHLUNG ORIGIN OF DIFFUSE HIGH-ENERGY EMISSION FROM THE INNER 2 Degree-Sign Multiplication-Sign 1 Degree-Sign OF THE GALACTIC CENTER

The high-energy activity in the inner few degrees of the Galactic center is traced by diffuse radio, X-ray, and {gamma}-ray emission. The physical relationship between different components of diffuse gas emitting at multiple wavelengths is a focus of this work. We first present radio continuum observations using the Green Bank Telescope and model the nonthermal spectrum in terms of a broken power-law distribution of {approx}GeV electrons emitting synchrotron radiation. We show that the emission detected by Fermi is primarily due to nonthermal bremsstrahlung produced by the population of synchrotron emitting electrons in the GeV energy range interacting with neutral gas. The extrapolation of the electron population measured from radio data to low and high energies can also explain the origin of Fe I 6.4 keV line and diffuse TeV emission, as observed with Suzaku, XMM-Newton, Chandra, and the H.E.S.S. observatories. The inferred physical quantities from modeling multiwavelength emission in the context of bremsstrahlung emission from the inner {approx}300 Multiplication-Sign 120 pc of the Galactic center are constrained to have the cosmic-ray ionization rate {approx}1-10 Multiplication-Sign 10{sup -15} s{sup -1}, molecular gas heating rate elevating the gas temperature to 75-200 K, fractional ionization of molecular gas 10{sup -6}-10{sup -5}, large-scale magnetic field 10-20 {mu}G, the density of diffuse and dense molecular gas {approx}100 and {approx}10{sup 3} cm{sup -3} over 300 pc and 50 pc path lengths, and the variability of Fe I K{alpha} 6.4 keV line emission on yearly timescales. Important implications of our study are that GeV electrons emitting in radio can explain the GeV {gamma}-rays detected by Fermi and that the cosmic-ray irradiation model, like the model of the X-ray irradiation triggered by past activity of Sgr A*, can also explain the origin of the variable 6.4 keV emission from Galactic center molecular clouds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Monte Carlo Control Loops for Cosmic Shear Cosmology with DES Year 1

Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. The Monte Carlo Control Loops (MCCL) is a forward modelling method designed to tackle this problem. It relies on the Ultra Fast Image Generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. We present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_{\ell}$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, Point Spread Function (PSF), an Approximate Bayesian Computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and the fast estimation of the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a non-tomographic setting, we derive a set of $C_\ell$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $\Omega_m$, $\sigma_8$, and intrinsic alignment amplitude $A_{\rm{IA}}$, defined as $S_8D_{\rm{IA}} = \sigma_8(\Omega_m/0.3)^{0.5}D_{\rm{IA}}$, where $D_{\rm{IA}}=1-0.11(A_{\rm{IA}}-1)$. We find $S_8D_{\rm{IA}}=0.895^{+0.054}_{-0.039}$, where systematics are at the level of roughly 60\% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Monte Carlo control loops for cosmic shear cosmology with DES Year 1 data

Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. Monte Carlo control loops (MCCL) is a forward modeling method designed to tackle this problem. It relies on the ultra fast image generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. Here, we present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_ℓ$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, point spread function (PSF), an approximate Bayesian computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and a fast method to estimate the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a nontomographic setting, we derive a set of $C_ℓ$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $Ω_m$, $σ_8$, and intrinsic alignment amplitude $A_{IA}$, defined as $S_8D_{IA}=σ_8(Ω_m/0.3)^{0.5}D_{IA}$, where $D_{IA}=1-0.11(A_{IA}-1)$. We find $S_8D_{IA}=0.8954_{-0.039}^{+0.054}$, where systematics are at the level of roughly 60% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Analysis of tipping-curve measurements performed at the DSS-13 beam-waveguide antenna at 32.0 and 8.45 GigaHertz

This article reports on the analysis of the Ka-band Antenna Performance Experiment tipping-curve data acquired at the DSS-13 research and development beam-waveguide (BWG) antenna. By measuring the operating system temperatures as the antenna is moved form zenith to low-elevation angles and fitting a model to the data, one can obtain information on how well the overall temperature model behaves at zenith and approximate the contribution due to the atmosphere. The atmospheric contribution estimated from the data can be expressed in the form of (1) atmospheric noise temperatures that can provide weather statistic information and be compared against those estimated from other methods and (2) the atmospheric loss factor used to refer efficiency measurements to zero atmosphere. This article reports on an analysis performed on a set of 68 8.4-GHz and 67 32-GHz tipping-curve data sets acquired between December 1993 and May 1995 and compares the results with those inferred from a surface model using input meteorological data and from water vapor radiometer (WVR) data. The general results are that, for a selected subset of tip curves, (1) the BWG tipping-curve atmospheric temperatures are in good agreement with those determined from WVR data (the average difference is 0.06 +/- 0.64 K at 32 GHz) and (2) the surface model average values are biased 3.6 K below those of the BWG and WVR at 32 GHz.

Morabito, D. D.↗
DOE OSTI2023

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Laser pulse-length dependent ablation and shock generation in silicon at 5 × 10 14 W/ cm 2 intensities

The effect of laser pulse duration on energy coupling into a planar silicon target is investigated in experiments at the OMEGA-EP facility by varying the laser pulse length τ —spanning 3 orders of magnitude from 100 ps to 10 ns—while maintaining a constant peak laser intensity, I 0 = 5 × 10 14 W/ cm 2 . In theoretical models, the ablation pressure primarily scales for a given material with laser intensity and wavelength, which are all fixed variables here, allowing us to explore the specific role of laser pulse duration. Two-dimensional radiation-hydrodynamics simulations benchmarked with optical probing of the expanding plasma show that the pulse duration is critical for the ablation pressure to reach a steady state. Moreover, the pulse duration impacts shock decay and multiple wave effects, which strongly dictate the evolving shock profile that propagates within the laser-shocked target as ultimately measured by rear-surface diagnostics. The shock velocities inferred from the theoretical model, after considering shock decay, impedance matching, and shock Hugoniot, are found to be in good agreement with velocimetry measurements. However, discrepancies are observed with simulations for the shorter (0.1 ns) and longer (10 ns) pulse durations, which are respectively attributed to unaccounted contributions of kinetic absorption mechanisms and instabilities in simulations. Published by the American Physical Society 2024

Bailly-Grandvaux, M. (ORCID:0000000175294013)↗
DOE OSTI2018

At the Nexus of History, Ecology, and Hydrobiogeochemistry: Improved Predictions across Scales through Integration

To improve predictions of ecosystem function in future environments, we need to integrate the ecological and environmental histories experienced by microbial communities with hydrobiogeochemistry across scales. A key issue is whether we can derive generalizable scaling relationships that describe this multiscale integration. There is a strong foundation for addressing these challenges. We have the ability to infer ecological history with null models and reveal impacts of environmental history through laboratory and field experimentation. Recent developments also provide opportunities to inform ecosystem models with targeted omics data. A major next step is coupling knowledge derived from such studies with multiscale modeling frameworks that are predictive under non-steady-state conditions. This is particularly true for systems spanning dynamic interfaces, which are often hot spots of hydrobiogeochemical function. We can advance predictive capabilities through a holistic perspective focused on the nexus of history, ecology, and hydrobiogeochemistry.

Stegen, James C. (ORCID:0000000191357424)↗