Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

PISCES two-detector covariance matrix fit for the NOvA Experiment

NOvA is a long-baseline neutrino oscillation experiment with two functionally identical detectors: a Near Detector (ND) at Fermilab, placed 1 km from the neutrino source, and a Far Detector (FD) located 810 km away from the ND in Minnesota. NOvA's primary physics goals are the precise measurements of neutrino oscillation parameters $\theta_{23}$ and $\Delta m^2_{32}$ , determine the neutrino mass ordering, and constrain the value of $\delta_{CP}$, via the study of muon neutrino to electron neutrino oscillation. In the standard NOvA three-flavor analysis, oscillation parameters are extracted using an extrapolation technique in which the ND data constrain the FD prediction through a ratio method. While this allows for systematic uncertainties sharing the same effects in both detectors to cancel, it remains an FD-only fit and does not fully leverage the constraining power of the high-statistics ND. This analysis proposes a simultaneous ND+FD fit using the PISCES method. PISCES (Parameter Inference with Systematic Covariance and Exact Statistics) is a framework designed to support complex configurations such as a joint ND+FD fit. This allows PISCES to take full advantage of the ND data to directly constrain systematic uncertainties across all samples. In PISCES, systematic uncertainties are encoded in a fractional covariance matrix, and statistical uncertainties are handled with a Poisson likelihood, making the approach well suited for low-statistics samples. For interpretability, we further use a Newton–Raphson + PCA method to recover per-systematic pulls from the covariance formulation. This poster presents the full PISCES joint ND+FD fit for the NOvA three-flavor analysis, describes its implementation and evaluates its performance through extensive robustness tests and fake data studies. It also provides a comparison between the PISCES joint ND+FD results and the standard NOvA extrapolation method.

Rajaoalisoa, Miriama [Cincinnati U.] (ORCID:000000↗
DOE OSTI2024

Toward Accelerated Nuclear-physics Parameter Estimation from Binary Neutron Star Mergers: Emulators for the Tolman–Oppenheimer–Volkoff Equations

Abstract Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doing so is the computational cost involved in solving the neutron-star structure equations, known as the Tolman–Oppenheimer–Volkoff (TOV) equations. In this paper, we explore approaches to emulating solutions for the TOV equations: multilayer perceptrons (MLPs), Gaussian processes, and a data-driven variant of the reduced basis method (RBM). We implement these emulators for three different parameterizations of the nuclear EOS, each with a different degree of complexity represented by the number of model parameters. We find that our MLP-based emulators are generally more accurate than the other two algorithms, whereas the RBM results in the largest speedup with respect to the full high-fidelity TOV solver. We employ these emulators for a simple parameter inference using a potentially loud BNS observation and show that the posteriors predicted by our emulators are in excellent agreement with those obtained from the full TOV solver.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

High-precision Galaxy Clustering Predictions from Small-volume Hydrodynamical Simulations via Control Variates

Abstract Cosmological simulations of galaxy formation are an invaluable tool for understanding galaxy formation and its impact on cosmological parameter inference from large-scale structures. However, their high computational cost is a significant obstacle for running simulations that probe cosmological volumes comparable to those analyzed by contemporary large-scale structure experiments. In this work, we explore the possibility of obtaining high-precision galaxy clustering predictions from small-volume hydrodynamical simulations such as MillenniumTNG and FLAMINGO via control variates. In this approach, the hydrodynamical full-physics simulation is paired with a matched low-resolution gravity-only simulation. By learning the galaxy–halo connection from the hydrodynamical simulation and applying it to the gravity-only counterpart, one obtains a galaxy population that closely mimics the one in the more expensive simulation. One can then construct an estimator of galaxy clustering that combines the clustering amplitudes in the small-volume hydrodynamical and gravity-only simulations with clustering amplitudes in a large-volume gravity-only simulation. Depending on the galaxy sample, clustering statistic, and scale, this galaxy clustering estimator can have an effective volume of up to around 100 times the volume of the original hydrodynamical simulation in the nonlinear regime. With this approach, we can construct galaxy clustering predictions from existing simulations that are precise enough for mock analyses of next-generation large-scale structure surveys such as the Dark Energy Spectroscopic Instrument and the Legacy Survey of Space and Time.

Doytcheva, Alexandra (ORCID:0009000111254888)↗
DOE OSTI2025

Can Neutron Star Tidal Effects Obscure Deviations from General Relativity?

Abstract One of the main goals of gravitational-wave astrophysics is to study gravity in the strong-field regime and constrain deviations from general relativity (GR). Any such deviation affects not only binary dynamics and gravitational-wave emission but also the structure and tidal properties of compact objects. In the case of neutron stars, masses, radii, and tidal deformabilities can all differ significantly between different theories of gravity. Currently, the measurement uncertainties in neutron star radii and tidal deformabilities are quite large. However, much less is known about how the large uncertainty in the nuclear equation of state (EOS) might affect tests of GR using binary neutron star mergers. Conversely, using the wrong theory of gravity might lead to incorrect constraints on the nuclear EOS. Here, we study this problem within scalar–tensor (ST) theory. We apply the recently derived ℓ = 2 tidal Love numbers in this theory to parameter estimation of GW170817. Correspondingly, we test if physics beyond GR could bias measurements of the nuclear EOS and neutron star radii. We find that parameter inference for both the GR and ST cases returns consistent component masses and tidal deformabilities. The radius and the EOS posteriors, however, differ between the two theories, but neither is excluded by current observational limits. This indicates that measurements of the nuclear EOS may be biased and that deviations from GR could go undetected when analyzing current binary neutron star mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Interstellar clouds containing optically thin H2

The theory of Black and Delgarno that the relative populations of the excited rotational levels of H2 can be understood in terms of cascading following absorption in the Lyman and Werner bands is employed to infer the gas densities and radiation fields within diffuse interstellar clouds containing H2 that is optically thin in those bands. The procedure is described for computing the populations of the different rotation levels, the relative distribution among the different rotation levels of newly formed H2 is determined on the basis of five simplified models, and the rate of H2 formation is estimated. The results are applied to delta Ori, two components of iota Ori, the second components of rho Leo and zeta Ori, tau Sco, gamma Vel, and zeta Pup. The inferred parameters are summarized for each cloud.

Jura, M.↗
NASA NTRS1983

Data processing techniques used with MST radars: A review

The data processing methods used in high power radar probing of the middle atmosphere are examined. The radar acts as a spatial filter on the small scale refractivity fluctuations in the medium. The characteristics of the received signals are related to the statistical properties of these fluctuations. A functional outline of the components of a radar system is given. Most computation intensive tasks are carried out by the processor. The processor computes a statistical function of the received signals, simultaneously for a large number of ranges. The slow fading of atmospheric signals is used to reduce the data input rate to the processor by coherent integration. The inherent range resolution of the radar experiments can be improved significant with the use of pseudonoise phase codes to modulate the transmitted pulses and a corresponding decoding operation on the received signals. Commutability of the decoding and coherent integration operations is used to obtain a significant reduction in computations. The limitations of the processors are outlined. At the next level of data reduction, the measured function is parameterized by a few spectral moments that can be related to physical processes in the medium. The problems encountered in estimating the spectral moments in the presence of strong ground clutter, external interference, and noise are discussed. The graphical and statistical analysis of the inferred parameters are outlined. The requirements for special purpose processors for MST radars are discussed.

Rastogi, P. K.↗
NASA NTRS1985

Discrete Random Media Techniques for Microwave Modeling of Vegetated Terrain

Microwave remote sensing of vegetated terrain has been studied. Vegetation is modeled so that backscattered radar signals can be used to infer parameters which characterize the vegetation and underlying ground. The vegetation is modeled by discrete lossy dielectric scatterers with prescribed characteristics. The goal of the modeling effort is to remotely sense vegetation type (classification), growth stage, and plant/ground moisture. This information can then be used as input into agricultural, forestry and global circulation models. The microwave frequency spectrum, particularly L and C bands, are especially appropriate for this purpose since the wavelength is comparable to plant leaf and stem size. The resulting resonant interaction leads to backscattered data highly depend on plant shape and orientation. In addition, the transparent nature of the atmosphere in this frequency regime allows for algorithm development which requires no atmospheric correction.

Lang, R. H.↗
NASA NTRS1986

Satellite 'wakes' and the orbit of the Encke Gap moonlet

A general model is developed for the quasi-periodic optical depth fluctuations observed in photopolarimeter and radio occultation profiles near the Encke Gap of Saturn's rings, which may be the gravitational 'wakes' of a moonlet orbiting within the same gap. The results obtained are used to estimate the moonlet's orbit and mass from the occultation data. In view of the present treatment's neglect of particle collisions and self-gravity, the longitude is consistent with other estimates based on a photopolarimeter scan exterior to the gap. The consistency of the orbit parameters inferred from photopolarimeter, radio, and wavy edge data strongly suggests that a single, dominant moonlet orbits within the gap.

Showalter, M. R.↗
NASA NTRS1993

BBXRT and Ginga observations of the Seyfert 1 galaxy Markarian 335

We present new X-ray observations of Mrk 335 in the 0.4-10 keV band using BBXRT and in the 2-20 keV band with Ginga. The BBXRT data show a break to a softer spectrum below about 1 keV. The Ginga data show evidence for variability in absorption, uncorrelated with flux changes of the continuum source. The inconsistency between the opacity observed at the iron K edge and at softer X-ray energies implies that the absorber has an ionization state about FeXV. We suggest that the nuclear emission can be explained by a self-consistent pair cascade model as developed by Zdziarski et al. (1990) and Zdziarski and Coppi (1991), where the compactness parameter inferred from the spectral fits to both data sets is in good agreement with the value calculated from the ratio of observed X-ray luminosity and variability time scale.

Turner, T. J.↗
NASA NTRS1994

Radiative acceleration in outflows from broad absorption line quasi-stellar objects. 2: Wind models

We investigate the dynamics of radiatively driven broad absorption-line (BAL) outflows in quasi-stellar objects (QSOs) by developing radial and time-independent numerical models. Two limits are explored. The first assumes that the absorbing matter is not forced to comove with the substrate, which provides pressure confinement. This assumption allows us to explore in detail a case in which the acceleration is entirely due to radiation pressure. Using the parameters inferred from observations, we find that under these conditions radiative acceleration (mainly due to resonance line scattering) can readily accelerate the flow to the observed velocities. An important feature of the noncoupled flow is that the line profiles tend to stay relatively flat throughout the velocity interval covered by the line. We discuss how relaxing the assumptions of radial symmetry and time independent may help to explain the structures observed in BALs. In the second class of models, the absorbing flow is assumed to be completely coupled to the substrate in which it is embedded. Aside from being more plausible physically, these models produce line profiles that trail off at higher velocities, a behavior observed in some BALs. We show that, even if the substrate is massless, we have to assume a starting radius very close to the inferred radius of the broad emission-line region (approximately 0.1 pc) in order to obtain a significant contribution from radiative acceleration, given a typical active galactic nucleus (AGN) spectrum. The reason is that the energy input needed to pressurize the substrate, allowing the flow to become supersonic and to retain a reasonable ionization equilibrium, at the same time contributes appreciably to the acceleration. A way to relax the small starting radius constraint is to use a softer ionizing spectrum.

Arav, Nahum↗
NASA NTRS2000

The nu Andromedae System: Models and Stability

Radial velocity observations of the F8 V star nu Andromedae taken at Lick and at Whipple Observatories have revealed evidence of three periodicities in the line-of-sight velocity of the star. These periodicities have been interpreted as evidence for at least three low-mass companions (LMCs) revolving around nu Andromedae. The mass and orbital parameters inferred for these companions raise questions about the dynamical stability of the system. We report here results from our independent analysis of the published radial velocity data, as well as new unpublished data taken at Lick Observatory. Our results confirm the finding of three periods in the data. Our best fits to the data, on the assumption that these periods arise from the gravitational perturbations of companions in Keplerian orbits, are also generally in agreement but with some differences from the earlier findings. We find that the available data do not constrain well the orbital eccentricity of the middle companion in a three-companion model of the data. We also find that in order for our best-fit model to the Lick data to be dynamically stable over the lifetime of the star (approximately 2 billion years), the system must have a mean inclination to the plane of the sky greater than 13 deg. The corresponding minimum inclination for the best fit to the Whipple data set is 19 deg. These values imply that the maximum mass for the outer companion can be no greater than about 20 Jupiter masses. Our analysis of the stability of the putative systems also places constraints on the relative inclinations of the orbital planes of the companions. We comment on global versus local (i.e., method of steepest descent) means of finding best-fit orbits from radial velocity data sets.

Stepinski, Tomasz F.↗
NASA NTRS1991

The Au Cathode in the System Li2CO3-CO2-CO at 800 to 900 C

Gold is one of several metals being evaluated at NASA Lewis Research Center as positive electrode catalysts for an alkali metal/molten alkali metal carbonate/carbon dioxide electrochemical cell. Such a cell is proposed for CO2-rich planetary atmospheres such as those of Mars and Venus. Its application could be as a primary power supply, as a secondary power supply recharged either 'chemically' by replenishment of the alkali metal or electrochemically from a central station power source, or as a converter of carbon dioxide to oxygen via a complete electrochemical cycle. For the work being reported, lithium was assumed to be the alkali metal of choice for the negative electrode of the cell, and therefore molten lithium carbonate was the electrolyte used in the Au electrode experiments. Cathodic linear sweep voltammetry (LSV) was the primary analytical technique for evaluating the performance of the Au cathode. interest comprised the cell temperature and the total pressure and composition of the reactant gas. In the absence of operational difficulties, the effect of bubbling the reactant gas through the melt was also determined. On the basis of the variation of electrode performance with changes in these parameters, inferences have been made concerning the electrochemical and chemical processes at and near the electrode. The results of post-test micrographic analyses of the Au cathode are also presented. An attempt is then made to project from the experimental results to some relevant conclusions pertaining to a gold cathode in a practical alkali metal - carbon dioxide cell.

Hagedorn, Norman H.↗
NASA NTRS2010

Combined Heat Transfer in High-Porosity High-Temperature Fibrous Insulations: Theory and Experimental Validation

Combined radiation and conduction heat transfer through various high-temperature, high-porosity, unbonded (loose) fibrous insulations was modeled based on first principles. The diffusion approximation was used for modeling the radiation component of heat transfer in the optically thick insulations. The relevant parameters needed for the heat transfer model were derived from experimental data. Semi-empirical formulations were used to model the solid conduction contribution of heat transfer in fibrous insulations with the relevant parameters inferred from thermal conductivity measurements at cryogenic temperatures in a vacuum. The specific extinction coefficient for radiation heat transfer was obtained from high-temperature steady-state thermal measurements with large temperature gradients maintained across the sample thickness in a vacuum. Standard gas conduction modeling was used in the heat transfer formulation. This heat transfer modeling methodology was applied to silica, two types of alumina, and a zirconia-based fibrous insulation, and to a variation of opacified fibrous insulation (OFI). OFI is a class of insulations manufactured by embedding efficient ceramic opacifiers in various unbonded fibrous insulations to significantly attenuate the radiation component of heat transfer. The heat transfer modeling methodology was validated by comparison with more rigorous analytical solutions and with standard thermal conductivity measurements. The validated heat transfer model is applicable to various densities of these high-porosity insulations as long as the fiber properties are the same (index of refraction, size distribution, orientation, and length). Furthermore, the heat transfer data for these insulations can be obtained at any static pressure in any working gas environment without the need to perform tests in various gases at various pressures.

Daryabeigi, Kamran↗
NASA NTRS2020

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↗
NASA NTRS2020

Tidal Disruptions of Main Sequence Stars - II. Simulation Methodology and Stellar Mass Dependence of the Character of Full Tidal Disruptions

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 (0.15 < M_/M<10)and six black hole masses( ́10<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/Me5×106this 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 of 0.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↗
NASA NTRS2023

CME Propagation Through the Heliosphere: Status and Future of Observations and Model Development

The ISWAT (International Space Weather Action Teams) heliosphere clusters H1 and H2 have a focus on interplanetary space and its characteristics, especially on the large-scale co-rotating and transient structures impacting Earth. Solar wind stream interaction regions, generated by the interaction between high-speed solar wind originating in large-scale open coronal magnetic fields and slower solar wind from closed magnetic fields, are regions of compressed plasma and magnetic field followed by high-speed streams that recur at the ~27 day solar rotation period. Short-term reconfigurations of the lower coronal magnetic field generate flare emissions and provide the energy to accelerate enormous amounts of magnetised plasma and particles in the form of coronal mass ejections into interplanetary space. The dynamic interplay between these phenomena changes the configuration of interplanetary space on various temporal and spatial scales which in turn influences the propagation of individual structures. While considerable efforts have been made to model the solar wind, we outline the limitations arising from the rather large uncertainties in parameters inferred from observations that make reliable predictions of the structures impacting Earth difficult. Moreover, the increased complexity of interplanetary space as solar activity rises in cycle 25 is likely to pose a challenge to these models. Combining observational and modeling expertise will extend our knowledge of the relationship between these different phenomena and the underlying physical processes, leading to improved models and scientific understanding and more-reliable space-weather forecasting. The current paper summarizes the efforts and progress achieved in recent years, identifies open questions, and gives an outlook for the next 5–10 years. It acts as basis for updating the existing COSPAR roadmap by Schrijver et al. (2015), as well as providing a useful and practical guide for peer-users and the next generation of space weather scientists.

Space weather↗
NASA NTRS1989

Inference of geophysical parameters from multifrequency polarimetric radar observations and model inversion

In summer 1988 an intensive field experiment was conducted in the vicinity of the Pisgah lava flow in the Mojave Desert. Physical properties such as microtopography, composition, soil moisture, and dielectric constant were measured at five sites representing surfaces with rms heights varying from less than one centimeter to tens of centimeters. In addition, polarimetric radar images at P-band, L-band, and C-band were acquired at three different incidence angles with the NASA/JPL airborne imaging radar polarimeter. Using trihedral corner reflectors deployed in the area prior to imaging, the L- and C-band images were calibrated to provide sigma0 values for each resolution element in the scene. The results of inferring geophysical parameters such as rms surface height and correlation length of the measured surfaces by fitting the observed signatures with those predicted by the small-perturbation model are presented: for smoother surfaces, the rms height values inferred are in good agreement with in situ measurements.

Van Zyl, J. J.↗
NASA NTRS1988

Inference of geologic surface parameters from polarimetric radar observations and model inversion

The results of inferring geologic parameters such as rms surface height, correlation length, and dielectric constant of rough surfaces by fitting observed polarization signatures with those predicted by the second order Rice model are discussed. The inferred results are compared to measured values of rms height and correlation length. The rms height values inferred are in good agreement with in situ measurements. The inferred correlation lengths generally do not agree with measured values. The results allow the separation of the effects of surface roughness and dielectric constant on the overall backscatter from rough surfaces.

Vanzyl, Jakob J.↗