SEARCH · Engineering Papers

Results for “likelihood function”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Extraction of Lateral-Directional Stability and Control Derivatives for the Basic F-18 Aircraft at High Angles of Attack

The results of parameter identification to determine the lateral-directional stability and control derivatives of an F-18 research aircraft in its basic hardware and software configuration are presented. The derivatives are estimated from dynamic flight data using a specialized identification program developed at NASA Dryden Flight Research Center. The formulation uses the linearized aircraft equations of motions in their continuous/discrete form and a maximum likelihood estimator that accounts for both state and measurement noise. State noise is used to model the uncommanded forcing function caused by unsteady aerodynamics, such as separated and vortical flows, over the aircraft. The derivatives are plotted as functions of angle of attack between 3 deg and 47 deg and compared with wind-tunnel predictions. The quality of the derivative estimates obtained by parameter identification is somewhat degraded because the maneuvers were flown with the aircraft's control augmentation system engaged, which introduced relatively high correlations between the control variables and response variables as a result of control motions from the feedback control system.

Iliff, Kenneth W.↗

NASA NTRS2012

Achievable Precision for Optical Ranging Systems

Achievable RMS errors in estimating the phase, frequency, and intensity of a direct-detected intensity-modulated optical pulse train are presented. For each parameter, the Cramer-Rao-Bound (CRB) is derived and the performance of the Maximum Likelihood estimator is illustrated. Approximations to the CRBs are provided, enabling an intuitive understanding of estimator behavior as a function of the signaling parameters. The results are compared to achievable RMS errors in estimating the same parameters from a sinusoidal waveform in additive white Gaussian noise. This establishes a framework for a performance comparison of radio frequency (RF) and optical science. Comparisons are made using parameters for state-of-the-art deep-space RF and optical links. Degradations to the achievable errors due to clock phase noise and detector jitter are illustrated.

estimation theory↗

DOE OSTI2023

Fitting Matérn smoothness parameters using automatic differentiation

The Mat$\acute{e}$rn covariance function is ubiquitous in the application of Gaussian processes to spatial statistics and beyond. Perhaps the most important reason for this is that the smoothness parameter $\nu$ gives complete control over the mean-square differentiability of the process, which has significant implications for the behavior of estimated quantities such as interpolants and forecasts. Unfortunately, derivatives of the Mat$\acute{e}$rn covariance function with respect to $\nu$ require derivatives of the modified second-kind Bessel function $K$ $\nu$ with respect to $\nu$. While closed form expressions of these derivatives do exist, they are prohibitively difficult and expensive to compute. For this reason, many software packages require fixing $\nu$ as opposed to estimating it, and all existing software packages that attempt to offer the functionality of estimating $\nu$ use finite difference estimates for $\partial$ $\nu$ $K$ $\nu$ . In this work, we introduce a new implementation of $K$$\nu$ that has been designed to provide derivatives via automatic differentiation (AD), and whose resulting derivatives are significantly faster and more accurate than those computed using finite differences. Here, we provide comprehensive testing for both speed and accuracy and show that our AD solution can be used to build accurate Hessian matrices for second-order maximum likelihood estimation in settings where Hessians built with finite difference approximations completely fail.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1986

Process for detecting the likelihood of ice formation, ice warning system for carrying out the process, and utilization thereof

An ice warning system comprising a diaphragm set into vibration at one of its frequencies of resonance by a piezoelectric cell. The diaphragm is coupled thermally to a metallic plate alternately cooled and heated respectively below and above the ambient temperature by Peltier elements connected to a reversible DC current generator. A microprocessor measures any variation of the resonance frequency caused by a deposit of ice on the diaphragm during the cooling or heating periods and delivers an alarm signal if this variation of frequency reaches and/ or exceeds a predetermined value. The microprocessor controls the period, intensity and direction of the current delivered by the DC current generator as a function of the ambient temperature and the temperature of the diaphragm. The ice warning system is preferably utilized in the aeronautical field for indicating the likelihood of natural ice formation on the engines and the wings of aircraft.

Martin Lustenberger↗

NASA NTRS1963

STATISTICAL FILTERING OF SPACE NAVIGATION MEASUREMENTS

Statistical filter design theory indicates that if space navigation measurement uncertainties have gaussian distributions, optimum filter theory and maximum likelihood estimate are identical

SPACE NAVIGATION↗

DOE OSTI2026

Combined dark matter search towards dwarf spheroidal galaxies with Fermi -LAT, HAWC, H.E.S.S., MAGIC, and VERITAS

Dwarf spheroidal galaxies (dSphs) are excellent targets for indirect dark matter (DM) searches using gamma-ray telescopes because they are thought to have high DM content and a low astrophysical background. The sensitivity of these searches is improved by combining the observations of dSphs made by different gamma-ray telescopes. We present the results of a combined search by the most sensitive currently operating gamma-ray telescopes, namely: the satellite-borne Fermi -LAT telescope; the ground-based imaging atmospheric Cherenkov telescope arrays H.E.S.S., MAGIC, and VERITAS; and the HAWC water Cherenkov detector. Individual datasets were analyzed using a common statistical approach. Results were subsequently combined via a global joint likelihood analysis. We obtain constraints on the velocity-weighted cross section 〈σv〉 for DM self-annihilation as a function of the DM particle mass. This five-instrument combination allows the derivation of up to 2-3 times more constraining upper limits on 〈σv〉 than the individual results over a wide mass range spanning from 5 GeV to 100 TeV. Depending on the DM content modeling, the 95% confidence level observed limits reach 1.5×10 -24 cm 3 s -1 and 3.2×10 -25 cm 3 s -1 , respectively, in the τ + τ - annihilation channel for a DM mass of 2 TeV.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1977

Estimated probability density functions for the times between flashes in the storms of 12 September 1975, 26 August 1975, and 13 July 1976

A report is given to supplement the progress report of June 17, 1977. In that progress report gamma, lognormal, and Rayleigh probability density functions were fitted to the times between lightning flashes in the storms of 9/12/75, 8/26/75, and 7/13/76 by the maximum likelihood method. The goodness of fit is checked by the Kolmogoroff-Smirnoff test. Plots of the estimated densities along with normalized histograms are included to provide a visual check on the goodness of fit. The lognormal densities are the most peaked and have the highest tails. This results in the best fit to the normalized histogram in most cases. The Rayleigh densities have too broad and rounded peaks to give good fits. In addition, they have the lowest tails. The gamma densities fall inbetween and give the best fit in a few cases.

Tretter, S. A.↗

NASA NTRS1999

Failure Engineering Study and Accelerated Stress Test Results for the Mars Global Surveyor Spacecraft's Power Shunt Assemblies

Due to a post launch failure of a part a new plan for the Mars Global Surveyor was developed. This new plan involved the addition of many deep thermal cycles to the Power Shunt Assemblies (PSA's). This new plan exceeds the previous acceptance cold level, and fatigue life on packaging design. This presentation reviews the experiments that were used to test the capabilities of the PSA to function in the new situation. It also reviews the analyses preformed to verify the most likely failure mechanism, and the likelihood that these failures would impact the new mission requirements.

Gibbel, Mark↗

NASA NTRS2003

Software for Better Documentation of Other Software

The Literate Programming Extraction Engine is a Practical Extraction and Reporting Language- (PERL-)based computer program that facilitates and simplifies the implementation of a concept of self-documented literate programming in a fashion tailored to the typical needs of scientists. The advantage for the programmer is that documentation and source code are written side-by-side in the same file, reducing the likelihood that the documentation will be inconsistent with the code and improving the verification that the code performs its intended functions. The advantage for the user is the knowledge that the documentation matches the software because they come from the same file. This program unifies the documentation process for a variety of programming languages, including C, C++, and several versions of FORTRAN. This program can process the documentation in any markup language, and incorporates the LaTeX typesetting software. The program includes sample Makefile scripts for automating both the code-compilation (when appropriate) and documentation-generation processes into a single command-line statement. Also included are macro instructions for the Emacs display-editor software, making it easy for a programmer to toggle between editing in a code or a documentation mode.

Pinedo, John↗

NASA NTRS2022

Constraining Warm Dark Matter and Population III Stars with the Global 21 cm Signal

Upcoming ground- and space-based experiments may have sufficient accuracy to place significant constraints upon high-redshift star formation, Reionization, and dark matter (DM) using the global 21 cm signal of the intergalactic medium. In the early universe, when the relative abundance of low-mass DM halos was important, measuring the global signal would place constraints on the damping of structure formation caused by DM having a higher relic velocity (warm dark matter, or WDM) than in cold dark matter. Such damping, however, can be mimicked by altering the star formation efficiency (SFE) and can be difficult to detect because of the presence of Population III stars with unknown properties. We study these various cases and their degeneracies with the WDM mass parameter mX using a Fisher matrix analysis. We study the mX = 7 keV case and a star formation model that parameterizes the SFE as a strong function of halo mass and include several variations of this model along with three different input noise levels for the likelihood; we also use a minimum halo virial temperature for collapse near the molecular cooling threshold. We find that when the likelihood includes only Population II stars, mX is constrained to an uncertainty of ∼0.4 keV for all models and noise levels at the 68% confidence level. When the likelihood includes weak Population III stars, mX ∼ 0.3 keV, and if Population III star formation is relatively efficient, mX ∼ 0.1 keV uncertainty, with tight Population III star formation parameter constraints. Our results show that the global 21 cm signal is a promising test-bed for WDM models, even in the presence of strong degeneracies with astrophysical parameters.

Cosmology↗

DOE OSTI2025

Optimizers for stabilizing likelihood-free inference

A growing number of applications in particle physics and beyond use neural networks as unbinned likelihood ratio estimators applied to real or simulated data. Precision requirements on the inference tasks demand a high-level of stability from these networks, which are affected by the stochastic nature of training. We show how physics concepts can be used to stabilize network training through a physics-inspired optimizer. In particular, the energy conserving descent (ECD) optimization framework uses classical Hamiltonian dynamics on the space of network parameters to reduce the dependence on the initial conditions while also stabilizing the result near the minimum of the loss function. We develop a version of this optimizer known as , which has few free hyperparameters with limited ranges guided by physical reasoning. We apply to representative likelihood-ratio estimation tasks in particle physics and find on average that it out-performs the widely used Adam optimizer. We expect that ECD will be a useful tool for wide array of data-limited problems, where it is computationally expensive to exhaustively optimize hyperparameters and mitigate fluctuations with ensembling.

Monte Carlo methods↗

NASA NTRS1987

A new method of boundary parameter estimation for a two-dimensional diffusion system under noisy observations

The purpose of this paper is to establish a method for identifying unknown parameters involved in the boundary state of a class of diffusion systems under noisy observations. A mathematical model of the system dynamics is given by a two-dimensional diffusion equation. Noisy observations are made by sensors allocated on the system boundary. Starting with the mathematical model mentioned above, an online parameter estimation algorithm is proposed within the framework of the maximum likelihood estimation. Existence of the optimal solution and related necessary conditions are discussed. By solving a local variation of the cost functional with respect to the perturbation of parameters, the estimation mechanism is proposed in a form of recursive computations. Finally, the feasibility of the estimator proposed here is demonstrated through results of digital simulation experiments.

Sunahara, Y.↗

NASA NTRS1988

A method of boundary parameter estimation for a two-dimensional diffusion system under noisy observations

Sunahara, Y.↗

NASA NTRS1999

Flight-Determined, Subsonic, Lateral-Directional Stability and Control Derivatives of the Thrust-Vectoring F-18 High Angle of Attack Research Vehicle (HARV), and Comparisons to the Basic F-18 and Predicted Derivatives

The subsonic, lateral-directional, stability and control derivatives of the thrust-vectoring F-1 8 High Angle of Attack Research Vehicle (HARV) are extracted from flight data using a maximum likelihood parameter identification technique. State noise is accounted for in the identification formulation and is used to model the uncommanded forcing functions caused by unsteady aerodynamics. Preprogrammed maneuvers provided independent control surface inputs, eliminating problems of identifiability related to correlations between the aircraft controls and states. The HARV derivatives are plotted as functions of angles of attack between 10deg and 70deg and compared to flight estimates from the basic F-18 aircraft and to predictions from ground and wind tunnel tests. Unlike maneuvers of the basic F-18 aircraft, the HARV maneuvers were very precise and repeatable, resulting in tightly clustered estimates with small uncertainty levels. Significant differences were found between flight and prediction; however, some of these differences may be attributed to differences in the range of sideslip or input amplitude over which a given derivative was evaluated, and to differences between the HARV external configuration and that of the basic F-18 aircraft, upon which most of the prediction was based. Some HARV derivative fairings have been adjusted using basic F-18 derivatives (with low uncertainties) to help account for differences in variable ranges and the lack of HARV maneuvers at certain angles of attack.

Iliff, Kenneth W.↗

NASA NTRS2018

Random Variables with Moment-Matching Staircase Density Functions

This paper proposes a family of random variables for uncertainty modeling. The variables of interest have a bounded support set, and prescribed values for the first four moments. We present the feasibility conditions for the existence of any of such variables, and propose a class of variables that conforms to such constraints. This class is called staircase because the density of its members is a piecewise constant function. Convex optimization is used to calculate their distributions according to several optimality criteria, including maximal entropy and maximal log-likelihood. The flexibility and efficiency of staircases enable modeling phenomena having a possibly skewed and/or multimodal response at a low computational cost. Furthermore, we provide a means to account for the uncertainty in the distribution caused by estimating staircases from data. These ideas are illustrated by generating empirical staircase predictor models. We consider the case in which the predictor matches the sample moments exactly (a setting applicable to large datasets), as well as the case in which the predictor accounts for the sampling error in such moments (a setting applicable to sparse datasets). A predictor model for the dynamics of an aeroelastic airfoil subject to flutter instability is used as an example. The resulting predictor not only describes the system's response accurately, but also enables carrying out a risk analysis for safe flight.

Luis G. Crespo↗

DOE OSTI2024

Inoculated biocrust cover and functions diverged over a gradient of soil textures and water availability

Restoring biological crust (biocrust) in disturbed drylands is challenging due to the difficult environmental conditions, such as limited soil moisture, low soil nutrients, and extreme temperatures, that impede growth. Understanding how the key components of biocrust—mosses, lichens, and cyanobacteria—react to different environmental factors informs the optimal timing, locations, and species composition for biocrust reintroduction, thereby increasing the likelihood of establishment. Here, we inoculated soils with a diverse range of biocrust organisms, analogous to seeding an area with diverse vascular plant seeds, and varied environmental conditions to observe how these changes influenced the development and functions of reintroduced biocrust. We found that by manipulating soil texture and time spent wet, we can change the proportional cover of biocrust within a restoration‐like setting. Specifically, we found that 4 months after inoculation, finer textured soils that received more water become dominated by moss cover, while coarser textured soils with less water remained dominated by cyanobacteria cover, and the interactions between texture and time spent wet strongly influenced cover. We found biocrust morphological group cover had a small, but detectable, effect on ecosystem functions (soil stability and nitrogenase activity, a proxy for nitrogen fixation), but that environmental conditions had a stronger impact on the functions we measured. Manipulative experiments in controlled environments, like this one, can help elucidate the mechanisms underlying the establishment rate and patterns of biocrusts post‐inoculation, and inform implementation of inoculations in the field.

Environmental Sciences & Ecology↗

NASA NTRS1995

Gaussian statistics of the cosmic microwave background: Correlation of temperature extrema in the COBE DMR two-year sky maps

We use the two-point correlation function of the extrema points (peaks and valleys) in the Cosmic Background Explorer (COBE) Differential Microwave Radiometers (DMR) 2 year sky maps as a test for non-Gaussian temperature distribution in the cosmic microwave background anisotropy. A maximum-likelihood analysis compares the DMR data to n = 1 toy models whose random-phase spherical harmonic components a(sub lm) are drawn from either Gaussian, chi-square, or log-normal parent populations. The likelihood of the 53 GHz (A+B)/2 data is greatest for the exact Gaussian model. There is less than 10% chance that the non-Gaussian models tested describe the DMR data, limited primarily by type II errors in the statistical inference. The extrema correlation function is a stronger test for this class of non-Gaussian models than topological statistics such as the genus.

Kogut, A.↗

DOE OSTI2023

Measurement of the production of a $W$ boson in association with a charmed hadron in $pp$ collisions at $\sqrt{s}$ = $\mathrm{13 TeV}$ with the ATLAS detector

The production of a $W$ boson in association with a single charm quark is studied using 140 fb –1 of $\sqrt{s}$ = $13$ $\mathrm{TeV}$ proton-proton collision data collected with the ATLAS detector at the Large Hadron Collider. The charm quark is tagged by the presence of a charmed hadron reconstructed with a secondary-vertex fit. The $W$ boson is reconstructed from the decay to either an electron or a muon and the missing transverse momentum present in the event. The charmed mesons reconstructed are D + → K – π + π + and D *+ → D 0 π + → (K – π + ) π + and the charge conjugate decays in the fiducial regions where $p$ T ($e,μ$) > 30 GeV, |$η(e,μ)$| < 2.5, $p$ T (D (*) ) > 8 GeV, and |$η$(D (*) )| < 2.2. The integrated and normalized differential cross sections as a function of the pseudorapidity of the lepton from the $W$ boson decay, and of the transverse momentum of the charmed hadron, are extracted from the data using a profile likelihood fit. The measured total fiducial cross sections are ${σ}_{fid}^{OS – SS}$ (W – + D + ) = 50.2 ± 0.2 ${(stat)}_{–2.3}^{+2.4}$(syst) pb, ${σ}_{fid}^{OS – SS}$ (W + + D – ) = 48.5 ± 0.2 ${(stat)}_{–2.2}^{+2.3}$(syst) pb, ${σ}_{fid}^{OS – SS}$ (W – + D *+ ) = 51.1 ± 0.4 ${(stat)}_{–1.8}^{+1.9}$(syst) pb, ${σ}_{fid}^{OS – SS}$ (W + + D *– ) = 50.0 ± 0.4 ${(stat)}_{–1.8}^{+1.9}$(syst) pb. Results are compared with the predictions of next-to-leading-order quantum chromodynamics calculations performed using state-of-the-art parton distribution functions. Additionally, the ratio of charm to anticharm production cross sections is studied to probe the $s$-$\bar{s}$ quark asymmetry. The ratio is found to be ${R}_{c}^{±}$= 0.971 ± 0.006 (stat) ± 0.011 (syst). The ratio and cross-section measurements are consistent with the predictions obtained with parton distribution function sets that have a symmetric $s$-$\bar{s}$ sea, indicating that any $s$-$\bar{s}$ asymmetry in the Bjorken-x region relevant for this measurement is small.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗