SEARCH · Engineering Papers

Results for “Markov model”

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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Probabilistic Prognosis of Non-Planar Fatigue Crack Growth

Quantifying the uncertainty in model parameters for the purpose of damage prognosis can be accomplished utilizing Bayesian inference and damage diagnosis data from sources such as non-destructive evaluation or structural health monitoring. The number of samples required to solve the Bayesian inverse problem through common sampling techniques (e.g., Markov chain Monte Carlo) renders high-fidelity finite element-based damage growth models unusable due to prohibitive computation times. However, these types of models are often the only option when attempting to model complex damage growth in real-world structures. Here, a recently developed high-fidelity crack growth model is used which, when compared to finite element-based modeling, has demonstrated reductions in computation times of three orders of magnitude through the use of surrogate models and machine learning. The model is flexible in that only the expensive computation of the crack driving forces is replaced by the surrogate models, leaving the remaining parameters accessible for uncertainty quantification. A probabilistic prognosis framework incorporating this model is developed and demonstrated for non-planar crack growth in a modified, edge-notched, aluminum tensile specimen. Predictions of remaining useful life are made over time for five updates of the damage diagnosis data, and prognostic metrics are utilized to evaluate the performance of the prognostic framework. Challenges specific to the probabilistic prognosis of non-planar fatigue crack growth are highlighted and discussed in the context of the experimental results.

Leser, Patrick E.↗

NASA NTRS1991

Frequency weighted system identification and linear quadratic controller design

Application of filters for frequency weighting of Markov parameters (pulse response functions) is described in relation to system/observer identification. The time domain identification approach recovers a model which has a pulse response weighted according to frequency. The identified model is composed of the original system and filters. The augmented system is in a form which can be used directly for frequency weighted linear quadratic controller design. Data from either single or multiple experiments can be used to recover the Markov parameters. Measured acceleration signals from a truss structure are used for system identification and the model obtained is used for frequency weighted controller design. The procedure makes the identification and controler design complementary problems.

Horta, Lucas G.↗

NASA NTRS1980

Spectral characteristics of convolutionally coded digital signals

The power spectral density of the output symbol sequence of a convolutional encoder is computed for two different input symbol stream source models, namely, an NRZ signaling format and a first-order Markov source. The effects of alternate symbol inversion on this spectrum are also considered. The mathematical results are illustrated with many examples corresponding to optimal performance codes. It is demonstrated that only for the case of a purely random input source (e.g., NRZ data with equiprobable symbols) and a particular class of codes is the output spectrum identical to the input spectrum except for a frequency scaling (expansion) by the reciprocal of the code rate.

Divsalar, D.↗

NASA NTRS1980

An application of extremal field guidance in statistical low thrust guidance analysis

The subject of this paper is the application of extremal field guidance to the covariance analysis of low thrust retargeting guidance. The general procedure for low thrust navigation analysis is briefly reviewed. The basic equations of trajectory optimization and extremal field guidance are presented with emphasis on the low thrust problem. The adopted formulation includes thruster on-off times as control variables and accounts for the effects of model parameter variations. A method of solution of the guidance equations is outlined; and the application of the solution in a covariance analysis, including a procedure for the treatment of model parameter variations as either biases or first order Gauss-Markov stochastic processes, is discussed.

Jacobson, R. A.↗

NASA NTRS1980

The program of the Techniques Development Laboratory in aviation weather forecasting

Some operational and experimental products developed for aviation weather forecasting are briefly described. Experimental products include surface dew point, obstructions to vision, boundary layer model, computer-worded terminal forecasts, terminal alerting procedure, generalized equivalent Markov, and radar forecasts (0-2 hours).

Klein, W. H.↗

NASA NTRS1993

The effect of dissipative inhomogeneous medium on the statistics of the wave intensity

One of the main theoretical points in the theory of wave propagation in random medium is the derivation of closed form equations to describe the statistics of the propagating waves. In particular, in one dimensional problems, the closed form representation of the multiple scattering effects is important since it contributes in understanding such problems like wave localization, backscattering enhancement, and intensity fluctuations. In this the propagation of plane waves in a layer of one-dimensional dissipative random medium is considered. The medium is modeled by a complex permittivity whose real part is a constant representing the absorption. The one dimensional problem is mathematically equivalent to the analysis of a transmission line with randomly perturbed distributed parameters and a single mode lossy waveguide and the results can be used to study the propagation of radio waves through atmosphere and the remote sensing of geophysical media. It is assumed the scattering medium consists of an ensemble of one-dimensional point scatterers randomly positioned in a layer of thickness L with diffuse boundaries. A Poisson impulse process with density lambda is used to model the position of scatterers in the medium. By employing the Markov properties of this process an exact closed form equation of Kolmogorov-Feller type was obtained for the probability density of the reflection coefficient. This equation was solved by combining two limiting cases: (1) when the density of scatterers is small; and (2) when the medium is weakly dissipative. A two variable perturbation method for small lambda was used to obtain solutions valid for thick layers. These solutions are then asymptotically evaluated for small dissipation. To show the effect of dissipation, the mean and fluctuations of the reflected power are obtained. The results were compared with a lossy homogeneous medium and with a lossless inhomogeneous medium and the regions where the effect of absorption is not essential were discussed.

Saatchi, Sasan S.↗

NASA NTRS2011

Filtering Methods for Error Reduction in Spacecraft Attitude Estimation Using Quaternion Star Trackers

Precision attitude determination for recent and planned space missions typically includes quaternion star trackers (ST) and a three-axis inertial reference unit (IRU). Sensor selection is based on estimates of knowledge accuracy attainable from a Kalman filter (KF), which provides the optimal solution for the case of linear dynamics with measurement and process errors characterized by random Gaussian noise with white spectrum. Non-Gaussian systematic errors in quaternion STs are often quite large and have an unpredictable time-varying nature, particularly when used in non-inertial pointing applications. Two filtering methods are proposed to reduce the attitude estimation error resulting from ST systematic errors, 1) extended Kalman filter (EKF) augmented with Markov states, 2) Unscented Kalman filter (UKF) with a periodic measurement model. Realistic assessments of the attitude estimation performance gains are demonstrated with both simulation and flight telemetry data from the Lunar Reconnaissance Orbiter.

Calhoun, Philip C.↗

NASA NTRS2022

Alma Images the Eccentric HD 53143 Debris Disk

We present ALMA 1.3 mm observations of the HD 53143 debris disk—the first infrared or millimeter image produced of this ∼1 Gyr old solar analog. Previous HST STIS coronagraphic imaging did not detect flux along the minor axis of the disk, which could suggest a face-on geometry with two clumps of dust. These ALMA observations reveal a disk with a strikingly different structure. In order to fit models to the millimeter visibilities and constrain the uncertainties on the disk parameters, we adopt a Markov Chain Monte Carlo approach. This is the most eccentric debris disk observed to date with a forced eccentricity of 0.21 ± 0.02, nearly twice that of the Fomalhaut debris disk, and also displays an apocenter glow. Although this eccentric model fits the outer debris disk well, significant interior residuals remain, which may suggest a possible edge-on inner disk, which remains unresolved in these observations. Combined with the observed structure difference between HST and ALMA, these results suggest a potential previous scattering event or dynamical instability in this system. We also note that the stellar flux changes considerably over the course of our observations, suggesting flaring at millimeter wavelengths. Using simultaneous TESS observations, we determine the stellar rotation period to be 9.6 ± 0.1 days.

Radio astronomy↗

NASA NTRS1988

Reliability Model Generator for fault-tolerant systems

This paper discusses an analysis tool, the Reliability Model Generator, that reasons from structural and functional system design specifications to generate a reliability model for the system under investigation. The resultant model defines a system state space sufficient to characterize the effects of single and multiple component failures. This model may then be examined by the analyst or used as input to an existing reliability evaluation tool, the Semi-Markov Unreliability Range Evaluator (SURE), to compute numeric bounds for system reliability (i.e., safety, mission success, availability, etc.). In defining the input to the Reliability Model Generator, a separation of the component function from structural specification is proposed to allow easy modification for analysis of alternative architectures. A hierarchical system description paradigm promotes multiple abstractions, thereby enabling analysis at all phases of the design process. The work described in this paper has been supported under NASA contract NAS1-10899, Integrated Airframe/Propulsion Control System Architecture (IAPSA II). IAPSA II promotes a system engineering methodology and supporting analytical tools to evaluate flight control architecture configuration alternatives during early phases of the design cycle when the cost and schedule impact of design revisions is minimal. Emphasis is placed on traceability to ensure that the requirements drive the resulting design.

Catherine M McCann↗

NASA NTRS1990

Average reflected power from a one-dimensional slab of discrete scatterers

Reflection from a one-dimensional random medium of discrete scatterers is considered. The discrete scattering medium is modeled by a Poisson impulse process with concentration lambda. By employing the Markov property of the Poisson impulse process, an exact functional integro-differential equation of the Kolmogorov-Feller type is found for the average reflected power. Approximate solutions to this equation are obtained by regular perturbation and two variable expansion techniques in the limit of small lambda. The regular perturbation results is valid for small slab thicknesses, while the two-variable result is uniformly valid for any thickness. The two-variable result shows that as the slab size becomes infinite all of the incident power is reflected on the average.

Saatchi, Sasan S.↗

NASA NTRS2019

Atmospheric Correction of Satellite Ocean-Color Imagery During the PACE Era

The Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission will carry into space the Ocean Color Instrument (OCI), a spectrometer measuring at 5 nm spectral resolution in the ultraviolet (UV) to near infrared (NIR) with additional spectral bands in the shortwave infrared (SWIR), and two multi-angle polarimeters that will overlap the OCI spectral range and spatial coverage, i. e., the Spectrometer for Planetary Exploration (SPEXone) and the Hyper-Angular Rainbow Polarimeter (HARP2). These instruments, especially when used in synergy, have great potential for improving estimates of water reflectance in the post Earth Observing System (EOS) era. Extending the top-of-atmosphere (TOA) observations to the UV, where aerosol absorption is effective, adding spectral bands in the SWIR, where even the most turbid waters are black and sensitivity to the aerosol coarse mode is higher than at shorter wavelengths, and measuring in the oxygen A-band to estimate aerosol altitude will enable greater accuracy in atmospheric correction for ocean color science. The multi-angular and polarized measurements, sensitive to aerosol properties (e.g., size distribution, index of refraction), can further help to identify or constrain the aerosol model, or to retrieve directly water reflectance. Algorithms that exploit the new capabilities are presented, and their ability to improve accuracy is discussed. They embrace a modern, adapted heritage two-step algorithm and alternative schemes (deterministic, statistical) that aim at inverting the TOA signal in a single step. These schemes, by the nature of their construction, their robustness, their generalization properties, and their ability to associate uncertainties, are expected to become the new standard in the future. A strategy for atmospheric correction is presented that ensures continuity and consistency with past and present ocean-color missions while enabling full exploitation of the new dimensions and possibilities. Despite the major improvements anticipated with the PACE instruments, gaps/issues remain to be filled/tackled. They include dealing properly with whitecaps, taking into account Earth-curvature effects, correcting for adjacency effects, accounting for the coupling between scattering and absorption, modeling accurately water reflectance, and acquiring a sufficiently representative dataset of water reflectance in the UV to SWIR. Dedicated efforts, experimental and theoretical, are in order to gather the necessary information and rectify inadequacies. Ideas and solutions are put forward to address the unresolved issues. Thanks to its design and characteristics, the PACE mission will mark the beginning of a new era of unprecedented accuracy in ocean-color radiometry from space.

ocean color↗

NASA NTRS1993

Q-Markov covariance equivalent realization and its application to flexible structure identification

The Q-Markov covariance equivalent realization algorithm is applied to NASA's Minimast structure to identify state-space models for the purpose of designing closed-loop controllers, and laboratory test results are shown for the identification and for the closed-loop performance. The paper also presents for the first time a deterministic formulation of covariance parameters from a pulse response, a stochastic formulation of Markov parameters from a white-noise response, and a simple Q-Markov covariance equivalent realization formulation. This requires only pulse laboratory tests or only white-noise laboratory tests to generate Q-Markov covariance equivalent realizations.

Liu, Ketao↗

NASA NTRS1971

Bayesian recursive image estimation.

A procedure for recursively estimating images that are characterized statistically by the mean and correlation functions associated with the random process representing the brightness level is proposed for the case where the images are corrupted by additive noise. First, a dynamic model is developed with a response characteristic which matches that of the scanner output (the input of the estimator is the output of a horizontal line scanner) in a statistical sense. Such models have the form of an ordinary differential or difference equation with white noise input. An insignificant approximation is introduced by using a constant-coefficient model. The appropriate model is a vector valued difference equation with the solution representing a vector Markov process. The next step is to obtain the minimum mean square estimate of the image by using a Kalman filter. Since the image estimation is an interpolation problem, two successive runs over the observation are performed in opposite directions and the resultant estimates are averaged. Examples are included for illustration.

Nahi, N. E.↗

NASA NTRS1995

Interdecadal Variations in the Alaska Gyre

Climatic dynamic topography variations in the Alaska gyre during the period 1968-1990 are described with an objective analysis of more than 12000 STD and XBT stations, and COADS wind stress data Interannual the dynamic height and SST variations were correlated and were consistent with recently described large-scale climatic shifts in the North Pacific. The gyre was centered more to the east, circulation appeared stronger, and SST was lower during the early to mid-1970s than during the 1980s. The Aleutian low (NP and PNA indices) intensified during the interim, but the response did not appear as a gyre spinup. Instead, the associated wind stress anomalies forced a slowly varying dynamic height anomaly across the eastern and northern part of the gyre through Ekman convergence, which had the effect of displacing the gyre's low somewhat to the WSW in the 1980s. The wind curl spectrum was white, and the slow oceanic response was modeled as stochastic-forced climate variability with a simple first-order Markov autoregression process. Forcing was assumed to be Ekman pumping of the pycnocline, and the damping coefficient was estimated from the data to be approx. 1 yr. A hindcast with observed winds gave estimated dynamic height patterns similar to those observed, with a canonical correlation of 0.79 at 99% confidence. This response was weak in the western half of the gyre, where slow baroclinic variability may have been influenced by long Rossby wave propagation. A simple autoregression simulation using artificial white noise forcing shows the evolution of decadal variations similar in nature to those observed. This result, along with the low frequency correlation between dynamic height and SST, suggests that the upper-ocean climatic variability in this region is primarily wind forced.

Lagerloef, Gary S.E.↗

NASA NTRS2011

Unsupervised Detection of Planetary Craters by a Marked Point Process

With the launch of several planetary missions in the last decade, a large amount of planetary images is being acquired. Preferably, automatic and robust processing techniques need to be used for data analysis because of the huge amount of the acquired data. Here, the aim is to achieve a robust and general methodology for crater detection. A novel technique based on a marked point process is proposed. First, the contours in the image are extracted. The object boundaries are modeled as a configuration of an unknown number of random ellipses, i.e., the contour image is considered as a realization of a marked point process. Then, an energy function is defined, containing both an a priori energy and a likelihood term. The global minimum of this function is estimated by using reversible jump Monte-Carlo Markov chain dynamics and a simulated annealing scheme. The main idea behind marked point processes is to model objects within a stochastic framework: Marked point processes represent a very promising current approach in the stochastic image modeling and provide a powerful and methodologically rigorous framework to efficiently map and detect objects and structures in an image with an excellent robustness to noise. The proposed method for crater detection has several feasible applications. One such application area is image registration by matching the extracted features.

Troglio, G.↗

NASA NTRS1988

The inclusion of semi-Markov reconfiguration transitions into the computer-aided Markov evaluator (CAME) program

The modifications to the rule-based CAME program which allow it to more accurately model the fault-handling processes of fault-tolerant systems are described. This new capability is added to the CAME program by modeling the fault-handling processes of fault-tolerant systems. The integrated airframe/propulsion control system architecture (IAPSA II) reference configuration currently under development is detailed.

Rosch, Gene↗

NASA NTRS1991

Comparison of several system identification methods for flexible structures

In the last few years various methods of identifying structural dynamics models from modal testing data have appeared. This paper presents a comparison of four of these algorithms: the Eigensystem Realization Algorithm (ERA), the modified version ERA/DC where DC indicates that it makes use of data correlations, the Q-Markov Cover algorithm, and an algorithm due to Moonen, DeMoor, Vandenberghe and Vandewalle. The comparison is made using a five mode computer model of the 20 meter Mini-Mast truss structure at NASA Langley Research Center, and various noise levels are superimposed to produce simulated data. The results show that for the example considered ERA/DC generally gives the best results; that ERA/DC is always at least as good as ERA which is shown to be a special case of ERA/DC; that Q-Markov requires the use of significantly more data than ERA/DC to produce comparable results; and that in some situations Q-Markov cannot produce comparable results.

Lew, J.-S.↗

NASA NTRS2016

Real-time Identification and Control of Satellite Signal Impairments Solution and Application of the Stratonovich Equation Part 1. Theoretical Development

As satellite communications systems become both more complex and reliant with respect to their operating environment, it has become imperative to be able to identify, during real-time operation, the onset of one or more impairments to the quality of overall communications system integrity. One of the most important aspects to monitor of a satellite link operating within the Earth's atmosphere is the signal fading due to the occurrence of rain and/or phase scintillations. This, of course, must be done in the presence of the associated measurement uncertainty or potentially faulty measurement equipment such as in the Advanced Communication Technology Satellite (ACTS) experiment. In the present work, an approach originally suggested in 1991, and apparently still considered iconoclastic, will be significantly developed and applied to the satellite communications link on which the deleterious composite signal fade is the result of one or many component fade mechanisms. Through the measurement (with the attendant uncertainty or 'error' in the measurement) of such a composite fading satellite signal, it is desired to extract the level of each of the individual fading mechanisms so they can be appropriately mitigated before they impact the overall performance of the communications network. Rather than employing simple-minded deterministic filtering to the real-time fading, the present approach is built around all the models and/or descriptions used to describe the individual fade components, including their dynamic evolution. The latter is usually given by a first-order Langevin equation. This circumstance allows the description of the associated temporal transition probability densities of each of the component processes. By using this description, along with the real-time measurements of the composite fade (along with the measurement errors), one can obtain statistical estimates of the levels of each of the component fading mechanisms as well as their predicted values into the future. This is all accomplished by the use of the well-known Stratonovich integro-differential equation that results from the model of the measured signal fade that is also tailored to adaptively adjust the values of the parameters used in the statistical models of the individual fade mechanisms. Three examples of increasing complexity are addressed and solved for the iterative determination of fade component levels from the measured composite signal fade in the presence of measurement error and, in the last case, with uncertainty in the model parameters.

Countermeasure↗