Search NASA⌕ Search

SEARCH · Search NASA

Results for “Profile likelihood”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 55 records · Page 3

Spacecraft Charging and the Microwave Anisotropy Probe Spacecraft

The Microwave Anisotropy Probe (MAP), a MIDEX mission built in partnership between Princeton University and the NASA Goddard Space Flight Center (GSFC), will study the cosmic microwave background. It will be inserted into a highly elliptical earth orbit for several weeks and then use a lunar gravity assist to orbit around the second Lagrangian point (L2), 1.5 million kilometers, anti-sunward from the earth. The charging environment for the phasing loops and at L2 was evaluated. There is a limited set of data for L2; the GEOTAIL spacecraft measured relatively low spacecraft potentials (approx. 50 V maximum) near L2. The main area of concern for charging on the MAP spacecraft is the well-established threat posed by the "geosynchronous region" between 6-10 Re. The launch in the autumn of 2000 will coincide with the falling of the solar maximum, a period when the likelihood of a substorm is higher than usual. The likelihood of a substorm at that time has been roughly estimated to be on the order of 20% for a typical MAP mission profile. Because of the possibility of spacecraft charging, a requirement for conductive spacecraft surfaces was established early in the program. Subsequent NASCAP/GEO analyses for the MAP spacecraft demonstrated that a significant portion of the sunlit surface (solar cell cover glass and sunshade) could have nonconductive surfaces without significantly raising differential charging. The need for conductive materials on surfaces continually in eclipse has also been reinforced by NASCAP analyses.

Timothy, VanSant J.↗

Transition in a disturbed environment

The title of this presentation is the title of our research grant. While transition study is the objective of the work, the results to date are principally on the properties of turbulent boundary layers at low Reynolds numbers are discussed. Testing was done in a small return wind tunnel. Mean boundary layer development is given. The skin friction behavior of the turbulent points are considered. No standard laminar flow was observed. Furthermore, the turbulent mean flow data seem reasonable for the elevated disturbance levels of the tests in the sense that there is no discernible wake component to any of the profiles and that the variation of skin friction with R sub theta is consistent with zero wake strength. The no-grid data are in all likelihood transitional.

Reshotko, E.↗

Detecting ocean surface wind using the TRMM precipitation radar

Although designed specifically for the measurement of precipitation in the atmosphere, the Precipitation Radar (PR) onboard the Tropical Rainfall Measuring Mission (TRMM) satellite alsomeasures the normalized radar cross-section at the Earth's surface. As such, this instrument providesan interesting opportunity to explore an alternative radar configuration in the satellite remote sensingnf ocean winds. In particular, it can measure the strength of the ocean winds in relatively higherspatial resolution as compared to the conventional scatterometer systems presently in use. The addition of wind sensing capability to precipitation radar also complements its rain profiling capabilitynicely and allows coincident wind and rain measurements. In this paper, we discuss the new ocean wind algorithm developed specifically for TRMM PR using maximum likelihood estimation. A uniquefeature of this wind algorithm is its capability to incorporate measurement sensitivity and noiseinformation consistently in both the along-track and cross-track directions. This PR wind algorithm hasbeen tested using data from the TRMM Microwave Imager (TMI) and NASA QuikSCAT scatterometer.Excellent agreement on the retrieved wind strengths is achieved among three sensors.

precipitation radar TRMM GPM CAMEX-4 ocean surface↗

Fan Beam Emission Tomography for Estimating Scalar Properties in Laminar Flames

A new method of estimating temperatures and gas species concentrations (CO2 and H2O) in a laminar flame is reported. The path-integrated, spectral radiation intensities emitted from a laminar flame at multiple wavelengths and view angles are calculated using a narrow band radiation model. Synthetic data, in the form of radial profiles of temperature and gas concentrations, are used in these calculations. The calculations mimic measurements that would theoretically be obtained using a mid-infrared spectrometer with a scanner. The path integrated spectral radiation intensities are deconvoluted using a maximum likelihood estimation method in conjunction with an iterative scheme. The deconvolution algorithm accounts for the self-absorption of radiation by the intervening gases, and provides the local temperature and gas species concentrations. The deconvoluted temperatures and gas concentrations are compared with the synthetic data used for calculating the spectral radiation intensities. The deconvoluted temperatures and gas species concentrations are within 0.5 % of the synthetic data. The deconvolution algorithm is expected to provide combustion researchers with an easy method of obtaining the radial profiles of major gas species concentrations and temperatures in laminar flames non-intrusively using a mid-infrared spectrometer with a scanner.

Lim, Jongmook↗

Uncertainty quantification of graph convolution neural network models of evolving processes

The application of neural network models to scientific machine learning tasks has proliferated in recent years. In particular, neural networks have proved to be adept at modeling processes with spatial–temporal complexity. Nevertheless, these highly parameterized models have garnered skepticism in their ability to produce outputs with quantified error bounds over the regimes of interest. Hence there is a need to find uncertainty quantification methods that are suitable for neural networks. In this work we present comparisons of the parametric uncertainty quantification of neural networks modeling complex spatial–temporal processes with Hamiltonian Monte Carlo and Stein variational gradient descent and its projected variant. Specifically we apply these methods to graph convolutional neural network models of evolving systems modeled with recurrent neural network and neural ordinary differential equations architectures. We show that Stein variational inference is a viable alternative to Monte Carlo methods with some clear advantages for complex neural network models. For our exemplars, Stein variational interference gave similar pushed forward uncertainty profiles through time compared to Hamiltonian Monte Carlo, albeit with generally more generous variance. As a result, projected Stein variational gradient descent also produced similar uncertainty profiles to the non-projected counterpart, but large reductions in the active weight space were confounded by the stability of the neural network predictions and the convoluted likelihood landscape.

36 MATERIALS SCIENCE↗

HDSense: An efficient method for ranking observable sensitivity

Identifying which observables most effectively constrain model parameters can be computationally prohibitive when considering full likelihoods of many correlated observables. This is especially important for, e.g., hadronization models, where high precision is required to interpret the results of collider experiments. We introduce the High-Dimensional Sensitivity (HDSense) score, a computationally efficient metric for ranking observable sets using only one-dimensional histograms. Derived by profiling over unknown correlations in the Fisher information framework, the score balances total information content against redundancy between observables. We apply HDSense to rank a set observables in terms of their constraining power with respect to five parameters of the Lund string model of hadronization implemented in Pythia using simulated leptonic collider events at the $Z$ pole. Validation against machine-learning--based full-likelihood approximations demonstrates that HDSense successfully identifies near-optimal observable subsets. The framework naturally handles data from multiple experiments with different acceptances and incorporates detector effects. While demonstrated on hadronization models, the methodology applies broadly to generic parameter estimation problems where correlations are unknown or difficult to model.

Assi, Benoît [Cincinnati U.] (ORCID:00000003092433↗

Diurnal Cycles of Intense Thunderstorms As Seen By Satellites

Spaceborne sensors cannot directly measure aspects of severe weather occurring at ground-level. Tools from Geostationary Earth Orbit (GEO) and Low Earth Orbit (LEO) can be related to storm intensity and likelihood of severe weather. Visible and Infrared (GEO and LEO) measure cloud-top properties and horizontal structure. Some properties have empirical correlations with severe weather occurrence. Precipitation Radar (LEO) measures vertical profile of reflectivity related to particle size and concentration, but attenuates (2.2-cm and shorter wavelengths) well above the surface in intense thunderstorms. Lightning mappers (GEO and LEO) measure lightning flash rates and flash properties, which relate to updraft properties including size and vigor. SAR and optical sensors can detect damage paths on the surface after a storm has passed. (Depends on the vegetation type and maturity.)

Daniel J Cecil↗

The Atacama Cosmology Telescope (ACT): Beam Profiles and First SZ Cluster Maps

The Atacama Cosmology Telescope (ACT) is currently observing the cosmic microwave background with arcminute resolution at 148 GHz, 218 GHz, and 277 GHz, In this paper, we present ACT's first results. Data have been analyzed using a maximum-likelihood map-making method which uses B-splines to model and remove the atmospheric signal. It has been used to make high-precision beam maps from which we determine the experiment's window functions, This beam information directly impacts all subsequent analyses of the data. We also used the method to map a sample of galaxy clusters via the Sunyaev-Ze1'dovich (SZ) effect, and show five clusters previously detected with X-ray or SZ observations, We provide integrated Compton-y measurements for each cluster. Of particular interest is our detection of the z = 0.44 component of A3128 and our current non-detection of the low-redshift part, providing strong evidence that the further cluster is more massive as suggested by X-ray measurements. This is a compelling example of the redshift-independent mass selection of the SZ effect.

Hincks, A. D.↗

Modeling of solar oscillation power spectra

To produce accurate estimates of the line-profile parameters of a model used to represent the spectral features in a solar oscillation power spectrum, it is necessary to (1) select the appropriate probability density function when deriving the maximum-likelihood function to be employed for the parameter estimation and (2) allow for the redistribution of spectral power caused by gaps in the data string. This paper describes a maximum-likelihood method for estimating the model parameters (based on the observed power spectrum statistics) that accounts for redistribution of spectral power caused by gaps in the data string, by convolving the model with the power spectrum of the observed window function. The accuracy and reliability of the method were tested using both artificial and authentic solar oscillation power spectrum data. A comparison of this method with various least-squares techniques is also presented.

Anderson, Edwin R.↗

Profiles of Gamma-Ray Bursts and Their Component Pulses

One physically informative regularity of their otherwise heterogeneous ensemble, is that many Gamma-Ray Bursts consist of well defined pulses. To objectively quantify the temporal structure of BATSE bursts, we have developed an automatic modeling procedure that separates overlapping pulses and determines the energy-dependence of the pulse-shape parameters. No binning of photon arrival times is needed, so when applied to time-tagged events (TTE) the procedure captures variability information down to the shortest time scales present in the raw data. Maximizing the Bayesian likelihood function Pr(data/model) yields estimates of the model parameters, including the number of pulses present, and allows intercomparison of models of different forms. As with any nonlinear optimization, good initial guesses are crucial to avoid convergence to undesirable local minima. We find excellent initial pulse decompositions by wavelet-denoising a cumulative distribution of the raw photon arrival data; differentiation then gives a time profile mostly free of the systematic effects of degraded resolution (as in ordinary Fourier smoothing) and binning. We present statistical information on pulse rise-time, decay-time, peakedness, and amplitudes, plus their energy dependences - both within a single burst and for a large ensemble of bursts.

Scargle, Jeff D.↗

Modeling personnel turnover in the parametric organization

A model is developed for simulating the dynamics of a newly formed organization, credible during all phases of organizational development. The model development process is broken down into the activities of determining the tasks required for parametric cost analysis (PCA), determining the skills required for each PCA task, determining the skills available in the applicant marketplace, determining the structure of the model, implementing the model, and testing it. The model, parameterized by the likelihood of job function transition, has demonstrated by the capability to represent the transition of personnel across functional boundaries within a parametric organization using a linear dynamical system, and the ability to predict required staffing profiles to meet functional needs at the desired time. The model can be extended by revisions of the state and transition structure to provide refinements in functional definition for the parametric and extended organization.

Dean, Edwin B.↗

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Aerosol Observability and Predictability: From Research to Operations for Chemical Weather Forecasting. Lagrangian Displacement Ensembles for Aerosol Data Assimilation

A challenge common to many constituent data assimilation applications is the fact that one observes a much smaller fraction of the phase space that one wishes to estimate. For example, remotely sensed estimates of the column average concentrations are available, while one is faced with the problem of estimating 3D concentrations for initializing a prognostic model. This problem is exacerbated in the case of aerosols because the observable Aerosol Optical Depth (AOD) is not only a column integrated quantity, but it also sums over a large number of species (dust, sea-salt, carbonaceous and sulfate aerosols. An aerosol transport model when driven by high-resolution, state-of-the-art analysis of meteorological fields and realistic emissions can produce skillful forecasts even when no aerosol data is assimilated. The main task of aerosol data assimilation is to address the bias arising from inaccurate emissions, and Lagrangian misplacement of plumes induced by errors in the driving meteorological fields. As long as one decouples the meteorological and aerosol assimilation as we do here, the classic baroclinic growth of error is no longer the main order of business. We will describe an aerosol data assimilation scheme in which the analysis update step is conducted in observation space, using an adaptive maximum-likelihood scheme for estimating background errors in AOD space. This scheme includes e explicit sequential bias estimation as in Dee and da Silva. Unlikely existing aerosol data assimilation schemes we do not obtain analysis increments of the 3D concentrations by scaling the background profiles. Instead we explore the Lagrangian characteristics of the problem for generating local displacement ensembles. These high-resolution state-dependent ensembles are then used to parameterize the background errors and generate 3D aerosol increments. The algorithm has computational complexity running at a resolution of 1/4 degree, globally. We will present the result of assimilating AOD retrievals from MODIS (on both Aqua and TERRA satellites) from AERONET for validation. The impact on the GEOS-5 Aerosol Forecasting will be fully documented.

da Silva, Arlindo↗

Communications Blackout Predictions for Atmospheric Entry of Mars Science Laboratory

The Mars Science Laboratory (MSL) is expected to be a long-range, long-duration science laboratory rover on the Martian surface. MSL will provide a significant milestone that paves the way for future landed missions to Mars. NASA is studying options to launch MSL as early as 2009. There are three elements to the spacecraft; carrier (cruise stage), entry vehicle, and rover. The rover will have a UHF proximity link as the primary path for EDL communications and may have an X-band direct-to-Earth link as a back-up. Given the importance of collecting critical event telemetry data during atmospheric entry, it is important to understand the ability of a signal link to be maintained, especially during the period near peak convective heating. The received telemetry during entry (or played back later) will allow for the performance of the Entry-Descent-Landing technologies to be assessed. These technologies include guided entry for precision landing, a new sky-crane landing system and powered descent. MSL will undergo an entry profile that may result in a potential communications blackout caused by ionized particles for short periods near peak heating. The vehicle will use UHF and possibly X-band during the entry phase. The purpose of this rep0rt is to quantify or bound the likelihood of any such blackout at UHF frequencies (401 MHz) and X-band frequencies (8.4 GHz). Two entry trajectory scenarios were evaluated: a stressful entry trajectory to quantify an upper-bound for any possible blackout period, and a nominal trajectory to quantify likelihood of blackout for such cases.

communications↗

Communications Blackout Predictions for Atmospheric Entry of Mars Science Laboratory

The Mars Science Laboratory (MSL) is expected to be a long-range, long-duration science laboratory rover on the Martian surface. MSL will provide a significant milestone that paves the way for future landed missions to Mars. NASA is studying options to launch MSL as early as 2009. MSL will be the first mission to demonstrate the new technology of 'smart landers', which include precision landing and hazard avoidance in order to -land at scientifically interesting sites that would otherwise be unreachable. There are three elements to the spacecraft; carrier (cruise stage), entry vehicle, and rover. The rover will have an X-band direct-to-Earth (DTE) link as well as a UHF proximity link. There is also a possibility of an X-band proximity link. Given the importance of collecting critical event telemetry data during atmospheric entry, it is important to understand the ability of a signal link to be maintained, especially during the period near peak convective heating. The received telemetry during entry (or played back later) will allow for the performance of the Entry-Descent-Landing technologies to be assessed. These technologies include guided entry for precision landing, hazard avoidance, a new sky-crane landing system and powered descent. MSL will undergo an entry profile that may result in a potential communications blackout caused by ionized plasma for short periods near peak heating. The vehicle will use UHF and possibly X-band during the entry phase. The purpose of this report is to quantify or bound the likelihood of any such blackout at UHF frequencies (401 MHz) and X-band frequencies (8.4 GHz). Two entry trajectory scenarios were evaluated: a stressful entry trajectory to quantify an upper-bound for any possible blackout period, and a nominal likely trajectory to quantify likelihood of blackout for such cases.

UHF↗

An Ensemble Recentering Kalman Filter with an Application to Argo Temperature Data Assimilation into the NASA GEOS-5 Coupled Model

A two-step ensemble recentering Kalman filter (ERKF) analysis scheme is introduced. The algorithm consists of a recentering step followed by an ensemble Kalman filter (EnKF) analysis step. The recentering step is formulated such as to adjust the prior distribution of an ensemble of model states so that the deviations of individual samples from the sample mean are unchanged but the original sample mean is shifted to the prior position of the most likely particle, where the likelihood of each particle is measured in terms of closeness to a chosen subset of the observations. The computational cost of the ERKF is essentially the same as that of a same size EnKF. The ERKF is applied to the assimilation of Argo temperature profiles into the OGCM component of an ensemble of NASA GEOS-5 coupled models. Unassimilated Argo salt data are used for validation. A surprisingly small number (16) of model trajectories is sufficient to significantly improve model estimates of salinity over estimates from an ensemble run without assimilation. The two-step algorithm also performs better than the EnKF although its performance is degraded in poorly observed regions.

Particle Filter↗

Introduction to the IMPACT Probabilistic Risk and Tradespace Analysis Tool for Medical System Design

Background: Probabilistic risk analysis (PRA) is a method for estimating risk in complex engineered systems that, at a basic level, focuses on what can go wrong and the likelihood and consequences of those occurrences. NASA has used PRA as an integral component of medical system risk estimation and design for spaceflight. IMPACT (Informing Mission Planning via Analysis of Complex Tradespaces) is a novel tool to meet these goals for exploration missions. Overview: IMPACT performs hundreds of thousands of Monte Carlo simulations of missions to build aggregate pictures of medical risk. These simulations are based on 120 possible medical conditions (the IMPACT Condition List) selected in a consensus-based process because they are of highest likelihood and/or consequence for exploration spaceflight. The conditions are then tied to clinical capabilities which can be used for management (e.g., inserting an IV) and then to over 600 specific resources needed to deliver a capability (e.g., an angiocath or an ultrasound). Different mission profiles can be simulated with user-specified inputs such as mission duration, destination, number of crew and pre-existing medical conditions, and EVA frequency. While the IMPACT evidence base is designed for exploration environments, these user inputs allow the tool to be used across a broad range of missions. IMPACT’s primary outcome metrics include loss of crew life (LOCL, a measure of in-flight mortality due to medical conditions), need for evacuation (RTDC, return to definitive care), and crew disability (TTL, task time lost based on how medical conditions impact the ability to perform over 1000 specific exploration mission crew tasks). In addition to modeling medical risk, IMPACT also accepts user-specified constraints, such as limitations of mass or volume, and will output a recommended clinical capability set and specific medical resources that meet the mission constraints. Discussion: This abstract will provide an introduction to IMPACT and describe the nature of the underlying medical evidence. It will also detail potential use cases for how this tool can be utilized by NASA or commercial spaceflight providers.

Ben Easter↗

How Well Does the Latest Anthropomorphic Test Device Mimic Human Impact Responses?

One of the goals of the NASA Occupant Protection Group is to understand the human tolerance to dynamic loading. This knowledge has to come through indirect approaches such as existing human response databases, anthropometric test devices (ATD), animal testing, post‐mortem human subjects, and models. This study investigated the biofidelity of the National Highway Traffic Safety Administration's ATD named the THOR (test device for human occupant restraint). If THOR responds comparably to humans, then it could potentially be used as a human surrogate to help validate space vehicle requirements for occupant protection. The THOR responses to frontal and spinal impacts (ranging from 8 to 12 G with rise times of 40, 70, and 100 ms) were measured and compared to human volunteer responses (95 trials in frontal and 58 in spinal) previously collected by the U. S. Air Force on the same horizontal impact accelerator. The impact acceleration profiles tested are within the expected range of multi‐purpose crew vehicle (MPCV) landing dynamics. A correlation score was calculated for each THOR to human comparison using CORA (CORrelation and Analysis) software. A two‐parameter beta distribution model fit was obtained for each dependent variable using maximum likelihood estimation. For frontal impacts, the THOR head x‐acceleration peak response correlated with the human response at 8 and 10‐G 100 ms but not 10‐G 70 ms. The phase lagged the human response. Head z‐acceleration was not correlated. Chest x‐acceleration was in phase, had a higher peak response, and was well correlated with lighter subjects (Cora = 0.8 for 46 kg vs. Cora = 0.4 for 126 kg). Head x‐displacement had a leading phase. Several subjects responded with the same peak displacement but the mean of the group was lower. The shoulder x‐displacement was in phase but had higher peaks than the human response. For spinal impacts, the THOR head x‐acceleration was not well correlated. Head and chest z‐acceleration was in phase but had a higher peak response. Chest z‐acceleration was highly correlated with heavier subjects at lower G pulses (Cora = 0.86 for 125 kg at 8 G). The human response was variable in shoulder z‐displacement but the THOR was in phase and was comparable to the mean peak response. Head xand z‐displacement was in phase but had higher peaks. Seat pan forces were well correlated, were in phase, but had a larger peak response than most subjects. The THOR does not respond to frontal and spinal impacts exactly the same way that a human does. Some responses are well matched and others are not. Understanding the strengths and weaknesses of this ATD is an important first step in determining its usefulness in occupant protection at NASA

Newby, N.↗