Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stochastic approximation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Effect of Head-Down Bed Rest and Artificial Gravity Countermeasure on Cardiac Autonomic and Advanced Electrocardiographic Function

To study the effects of 21 days of head-down bed rest (HDBR), with versus without an artificial gravity (AG) countermeasure, on cardiac autonomic and advanced electrocardiographic function. Fourteen healthy men participated in the study: seven experienced 21 days of HDBR alone ("HDBR controls") and seven the same degree and duration of HDBR but with approximately 1hr daily short-arm centrifugation as an AG countermeasure ("AG-treated"). Five minute supine high-fidelity 12-lead ECGs were obtained in all subjects: 1) 4 days before HDBR; 2) on the last day of HDBR; and 3) 7 days after HDBR. Besides conventional 12-lead ECG intervals and voltages, all of the following advanced ECG parameters were studied: 1) both stochastic (time and frequency domain) and deterministic heart rate variability (HRV); 2) beat-to-beat QT interval variability (QTV); 3) T-wave morphology, including signal-averaged T-wave residua (TWR) and principal component analysis ratios; 4) other SAECG-related parameters including high frequency QRS ECG and late potentials; and 5) several advanced ECG estimates of left ventricular (LV) mass. The most important results by repeated measures ANOVA were that: 1) Heart rates, Bazett-corrected QTc intervals, TWR, LF/HF power and the alpha 1 of HRV were significantly increased in both groups (i.e., by HDBR), but with no relevant HDBR*group differences; 2) All purely "vagally-mediated" parameters of HRV (e.g., RMSSD, HF power, Poincare SD1, etc.), PR intervals, and also several parameters of LV mass (Cornell and Sokolow-Lyon voltages, spatial ventricular activation times, ventricular gradients) were all significantly decreased in both groups (i.e., by HDBR), but again with no relevant HDBR*group differences); 3) All "generalized" or "vagal plus sympathetic" parameters of stochastic HRV (i.e., SDNN, total power, LF power) were significantly more decreased in the AG-treated group than in the HDBR-only group (i.e., here there was a relevant HDBR*group difference); and 4) QTV index was also significantly more changed (increased) in the AG-treated group than in the HDBR-only group, although this was clearly due to a greater decrease in generalized HRV and not to a greater increase in QTV proper because there was no relevant HDBR*group effect for either the SDNN or the RMSSD of QTV. Brief daily AG treatment by short-arm centrifuge during each of 21 days of HDBR does not appear to protect against HDBR-related losses of cardiac autonomic function or of LV mass as estimated by ECG.

Schlegel, T. T.↗
NASA NTRS1996

Simulation of Rotor Blade Element Turbulence

A turbulence model has been developed for blade-element helicopter simulation. This model, called Simulation of Rotor Blade Element Turbulence (SORBET), uses an innovative temporal and geometrical distribution algorithm that preserves the statistical characteristics of the turbulence spectra over the rotor disc, while providing velocity components in real time to each of five blade-element stations along each of four blades. An initial investigation of SORBET has been performed using a piloted, motion-based simulation of the Sikorsky UH60A Black Hawk. Although only the vertical component of stochastic turbulence was used in this investigation, vertical turbulence components induce vehicle responses in all translational and rotational degrees of freedom of the helicopter. The single-degree-of-freedom configuration of SORBET was compared to a conventional full 6-degrees-of-freedom baseline configuration, where translational velocity inputs are superimposed at the vehicle center of gravity, and rotational velocity inputs are created from filters that approximate the immersion rate into the turbulent field. For high-speed flight the vehicle responses were satisfactory for both models. Test pilots could not distinguish differences between the baseline configuration and SORBET. In low-speed flight the baseline configuration received criticism for its high frequency content, whereas the SORBET model elicited favorable pilot opinion. For this helicopter, which has fully articulated blades, results from SORBET show that vehicle responses to turbulent blade-station disturbances are severely attenuated. This is corroborated by in-flight observation of the rotor tip path plane as compared to vehicle responses.

McFarland, R. E.↗
NASA NTRS2004

Employing Sensitivity Derivatives for Robust Optimization under Uncertainty in CFD

A robust optimization is demonstrated on a two-dimensional inviscid airfoil problem in subsonic flow. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables), an approximate first-order statistical moment method is employed to represent the Computational Fluid Dynamics (CFD) code outputs as expected values with variances. These output quantities are used to form the objective function and constraints. The constraints are cast in probabilistic terms; that is, the probability that a constraint is satisfied is greater than or equal to some desired target probability. Gradient-based robust optimization of this stochastic problem is accomplished through use of both first and second-order sensitivity derivatives. For each robust optimization, the effect of increasing both input standard deviations and target probability of constraint satisfaction are demonstrated. This method provides a means for incorporating uncertainty when considering small deviations from input mean values.

Newman, Perry A.↗
NASA NTRS2021

Stochastic Guidance of Buoyancy Controlled Vehicles under Ice Shelves using Ocean Currents

We propose a novel technique for guidance ofbuoyancy-controlled vehicles in uncertain under-ice ocean flows.In-situ melt rate measurements collected at the grounding zoneof Antarctic ice shelves, where the ice shelf meets the underlyingbedrock, are essential to constrain models of future sea levelrise. Buoyancy-controlled vehicles, which control their verticalposition in the water column but have no means of horizontalpropulsion, offer an affordable and reliable platform for suchin-situ data collection. However, reaching the grounding zonerequires vehicles to traverse tens of kilometers under the iceshelf, with approximate position knowledge and no meansof communication, in highly variable and uncertain oceancurrents. To address this challenge, we propose a partiallyobservable MDP approach that exploits model-based knowledgeof the under-ice currents and, critically, of their uncertainty,to synthesize effective guidance policies. The approach usesapproximate dynamic programming to model uncertainty inthe currents, and QMDP to address localization uncertainty.Numerical experiments show that the policy can deliver upto 88.8% of underwater vehicles to the grounding zone – a33% improvement compared to state-of-the-art guidance techniques, and a 262% improvement over uncontrolled drifters.Collectively, these results show that model-based under-iceguidance is a highly promising technique for exploration ofunder-ice cavities, and has the potential to enable cost-effectiveand scalable access to these challenging and rarely observed environments.

Clark, Evan B.↗
NASA NTRS2011

Solar Cycle Fine Structure and Surface Rotation from Ca II K-Line Time Series Data

Analysis of three and a half decades of data from the NSO/AFRL/Sac Peak K-line monitoring program yields evidence for four components to the variation: (a) the solar cycle, with considerable fine structure and a quasi-periodicity of 122.4 days; (b) a stochastic process, faster than (a) and largely independent of it, (c) a quasi-periodic signal due to rotational modulation, and of course (d) observational errors (shown to be quite small). Correlation and power spectrum analyses elucidate periodic and aperiodic variation of these chromospheric parameters. Time-frequency analysis is especially useful for extracting information about differential rotation, and in particular elucidates the connection between its behavior and fine structure of the solar cycle on approximately one-year time scales. These results further suggest that similar analyses will be useful at detecting and characterizing differential rotation in stars from stellar light-curves such as those being produced at NASA's Kepler observatory.

Scargle, Jeff↗
NASA NTRS1990

Stochastic Fermi acceleration in solar flares

Proton spectra, valid from non- to ultra-relativistic energies, resulting from stochastic Fermi acceleration in solar flares are calculated. These spectra were obtained by numerically solving the Fokker-Planck equation, in which the escape of the particles from the acceleration region is characterized by an energy-independent escape time. In addition to equilibrium spectra, time-dependent energy spectra showing the approach to equilibrium are also presented. These numerical equilibrium spectra are compared with previous results which were obtained either by Monte Carlo simulations or approximate analytical treatments. There are no analytic solutions valid in the transrelativistic regime, which is very important for the production of pions and neutrons in solar flares. The acceleration efficiency is related to physical parameters, in particular the energy density in either magnetosonic or Alfven waves, and a lower limit is placed on either of these energy densities from acceleration times implied by gamma-ray observations. Also discussed is the physical interpretation of the escape time.

Miller, James A.↗
NASA NTRS2015

Physical Models of Layered Polar Firn Brightness Temperatures from 0.5 to 2 GHz

We investigate physical effects influencing 0.5-2 GHz brightness temperatures of layered polar firn to support the Ultra Wide Band Software Defined Radiometer (UWBRAD) experiment to be conducted in Greenland and in Antarctica. We find that because ice particle grain sizes are very small compared to the 0.5-2 GHz wavelengths, volume scattering effects are small. Variations in firn density over cm- to m-length scales, however, cause significant effects. Both incoherent and coherent models are used to examine these effects. Incoherent models include a 'cloud model' that neglects any reflections internal to the ice sheet, and the DMRT-ML and MEMLS radiative transfer codes that are publicly available. The coherent model is based on the layered medium implementation of the fluctuation dissipation theorem for thermal microwave radiation from a medium having a nonuniform temperature. Density profiles are modeled using a stochastic approach, and model predictions are averaged over a large number of realizations to take into account an averaging over the radiometer footprint. Density profiles are described by combining a smooth average density profile with a spatially correlated random process to model density fluctuations. It is shown that coherent model results after ensemble averaging depend on the correlation lengths of the vertical density fluctuations. If the correlation length is moderate or long compared with the wavelength (approximately 0.6x longer or greater for Gaussian correlation function without regard for layer thinning due to compaction), coherent and incoherent model results are similar (within approximately 1 K). However, when the correlation length is short compared to the wavelength, coherent model results are significantly different from the incoherent model by several tens of kelvins. For a 10-cm correlation length, the differences are significant between 0.5 and 1.1 GHz, and less for 1.1-2 GHz. Model results are shown to be able to match the v-pol SMOS data closely and predict the h-pol data for small observation angles.

layered polar firn emission↗
NASA NTRS1996

Stochastic Filtering of Rain Profiles Using Radar, Surface-Referenced Radar, or Combined Radar-Radiometer Measurements

This paper describes a computationally efficient nearly optimal Bayesian algorithm to estimate rain (and drop size distribution) profiles, given a radar reflectivity profile at a single attenuating wavelength. In addition to estimating the averages of all the mutually ambiguous combinations of rain parameters that can produce the data observed, the approach also calculates the n-ns uncertainty in its estimates (this uncertainty thus quantifies "the amount of ambiguity" in the "solution"). The paper also describes a more general approach that can make estimates based on a radar reflectivity profile together with an approximate measurement of the path-integrated attenuation, or a radar reflectivity profile and a set of passive microwave brightness temperatures. This more general "combined" algorithm is currently being adapted for the Tropical Rainfall Measuring Mission.

Haddad, Ziad S.↗
NASA NTRS2007

POD Model Reconstruction for Gray-Box Fault Detection

Proper orthogonal decomposition (POD) is the mathematical basis of a method of constructing low-order mathematical models for the "gray-box" fault-detection algorithm that is a component of a diagnostic system known as beacon-based exception analysis for multi-missions (BEAM). POD has been successfully applied in reducing computational complexity by generating simple models that can be used for control and simulation for complex systems such as fluid flows. In the present application to BEAM, POD brings the same benefits to automated diagnosis. BEAM is a method of real-time or offline, automated diagnosis of a complex dynamic system.The gray-box approach makes it possible to utilize incomplete or approximate knowledge of the dynamics of the system that one seeks to diagnose. In the gray-box approach, a deterministic model of the system is used to filter a time series of system sensor data to remove the deterministic components of the time series from further examination. What is left after the filtering operation is a time series of residual quantities that represent the unknown (or at least unmodeled) aspects of the behavior of the system. Stochastic modeling techniques are then applied to the residual time series. The procedure for detecting abnormal behavior of the system then becomes one of looking for statistical differences between the residual time series and the predictions of the stochastic model.

Park, Han↗
NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

residual stress↗
NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

processing↗
NASA NTRS1994

Electromagnetic and corpuscular emission from the solar flare of 1991 June 15: Continuous acceleration of relativistic particles

Data on X-, gamma-ray, optical and radio emission from the 1991 June 15 solar flare are considered. We have calculated the spectrum of protons that produces gamma-rays during the gradual phase of the flare. The primary proton spectrum can be described as a Bessel-function-type up to 0.8 GeV and a power law with the spectral index approximately equal to 3 from 0.8 up to 10 GeV or above. We have also analyzed data on energetic particles near the Earth. Their spectrum differed from that of primary protons producing gamma-ray line emission. In the gradual phase of the flare additional pulses of energy release occurred and the time profiles of cm-radio emission and gamma-rays in the 0.8-10 MeV energy band and above 50 MeV coincided. A continuous and simultaneous stochastic acceleration of the protons and relativistic electrons at the gradual phase of the flare is considered as a natural explanation of the data.

Kocharov, L. G.↗
NASA NTRS1994

Inter-Agency Consultative Group for Space Science (IACG): Handbook of Missions and Payloads

The ACE spacecraft design is based on the Charge Composition Explorer (CCE) built by Johns Hopkins University (JHU) and the Applied Physics Lab (APL) for the AMPTE program. ACE is designed as a spinning spacecraft with its spin axis aligned to the Earth-Sun axis. The ACE launch weight will be approx. 633 kg, including 105 kg of scientific instruments and 184 kg of propellant. Using a Delta-class expendable launch vehicle, ACE will be launched into an L1 libration point (240 R(sub e)) orbit. Telemetry will be 6.7 kbps average, using tape recorder storage with daily readout to DSN. The experiment power requirement is approximately 76 W nominal and 96 W peak. The prime objective of the ACE mission is: (1) to determine accurate elemental and isotropic abundances including solar matter, local interstellar matter and local galactic matter; (2) to study the origin of elements and evolutionary processing in galactic nucleosynthesis, galactic evolution, origin and evolution of the solar system; (3) to study coronal formation and solar-wind acceleration processes; and (4) to study particle acceleration and transport, including coronal shock acceleration, stochastic flare acceleration, interplanetary shock acceleration, and interstellar acceleration and propagation. To accomplish this objective, ACE will perform comprehensive and coordinated determinations of the elemental and isotopic composition of energetic nuclei accelerated on the Sun, in interplanetary space, and from galactic sources. These observations will span five decades in energy, from solar wind to galactic cosmic ray energies, and will cover the element range from H-1 to Zr-40. Comparison of these samples of matter will be used to study the origin and subsequent evolution of both solar system and galactic material by isolating the effects of fundamental processes that include nucleosynthesis, charged and neutral particle separation, bulk plasma acceleration, and the acceleration of suprathermal and high-energy particles.

Source record↗
NASA NTRS2009

The Principle of Energetic Consistency

A basic result in estimation theory is that the minimum variance estimate of the dynamical state, given the observations, is the conditional mean estimate. This result holds independently of the specifics of any dynamical or observation nonlinearity or stochasticity, requiring only that the probability density function of the state, conditioned on the observations, has two moments. For nonlinear dynamics that conserve a total energy, this general result implies the principle of energetic consistency: if the dynamical variables are taken to be the natural energy variables, then the sum of the total energy of the conditional mean and the trace of the conditional covariance matrix (the total variance) is constant between observations. Ensemble Kalman filtering methods are designed to approximate the evolution of the conditional mean and covariance matrix. For them the principle of energetic consistency holds independently of ensemble size, even with covariance localization. However, full Kalman filter experiments with advection dynamics have shown that a small amount of numerical dissipation can cause a large, state-dependent loss of total variance, to the detriment of filter performance. The principle of energetic consistency offers a simple way to test whether this spurious loss of variance limits ensemble filter performance in full-blown applications. The classical second-moment closure (third-moment discard) equations also satisfy the principle of energetic consistency, independently of the rank of the conditional covariance matrix. Low-rank approximation of these equations offers an energetically consistent, computationally viable alternative to ensemble filtering. Current formulations of long-window, weak-constraint, four-dimensional variational methods are designed to approximate the conditional mode rather than the conditional mean. Thus they neglect the nonlinear bias term in the second-moment closure equation for the conditional mean. The principle of energetic consistency implies that, to precisely the extent that growing modes are important in data assimilation, this term is also important.

Cohn, Stephen E.↗
NASA NTRS1981

Blue compact dwarf galaxies. I - Neutral hydrogen observations of 115 galaxies

HI observations are reported for a sample of 115 blue compact dwarf (M sub B greater than approximately -18) galaxies or 'extragalactic H II regions' chosen mostly from the objective prism surveys of Markarian (1967-1974) and Haro (1956), with a few objects from Zwicky (1971) and other investigators. Ninety-three galaxies are detected. H I profiles, neutral hydrogen masses, total masses, and all available optical data are given for the 115 galaxies in a consistent and homogeneous system and in a useful format for statistical studies. The data are used in a companion paper to study the stochastic mode of star formation in galaxies.

Thuan, T. X.↗
NASA NTRS1995

Filtering Image Records Using Wavelets and the Zakai Equation

Consider the problem of detecting and localizing a faint object moving In an "essentially stationary" background, using a sequence of two-dimensional low-SNR images of the scene. A natural approach consists of "digitizing" each snapshot into a discrete set of observations, sufficiently (perhaps not exactly) matched to the object In question, then tracking the object using an appropriate stochastic filter. The tracking would be expected to make up for the low signal-to-noise ratio, this allowing one to "coherently" process successive images in order to beat down the noise and localize the object. Thus, "tracking" here does not refer to the usual notion of detecting then tracking: rather, we track in order to detect The problem then becomes one of choosing the appropriate image representation as well as the optimal (and necessarily nonlinear filter. We propose exact and approximate solutions using wavelets and the Zakai equation. The smoothness of the wavelets used is required in the derivation of the evolution equation for the conditional density giving the filter, and their orthogonality makes it possible to carry out actual computations of the Ito- and change-of-gauge-terms in the algorithm effectively.

Haddad, Ziad S.↗
NASA NTRS2010

Navier-Stokes Dynamics by a Discrete Boltzmann Model

This work investigates the possibility of particle-based algorithms for the Navier-Stokes equations and higher order continuum approximations of the Boltzmann equation; such algorithms would generalize the well-known Pullin scheme for the Euler equations. One such method is proposed in the context of a discrete velocity model of the Boltzmann equation. Preliminary results on shock structure are consistent with the expectation that the shock should be much broader than the near discontinuity predicted by the Pullin scheme, yet narrower than the prediction of the Boltzmann equation. We discuss the extension of this essentially deterministic method to a stochastic particle method that, like DSMC, samples the distribution function rather than resolving it completely.

Rubinstein, Robet↗
NASA NTRS1978

Minimum energy control of a class of electrically driven vehicles

Theoretical and experimental results are obtained in the investigation of a minimum energy control problem associated with a class of electrically driven vehicles. Analytical results are obtained by making several justifiable approximations in the dynamical equations of the plant, the performance index of which is related to the minimization of system energy consumption for any required control action. The control problem of interest is simplified and solved by using Bellman's invariant imbedding technique. The optimality of the resulting control law is compared with those of the bang-bang control law, which is the exact solution for the case of speed-setting control action, and the classical control laws under identical conditions. A stochastic optimization problem is then formulated and solved by using Bellman's dynamic programming technique. The resulting control law is mechanized for a particular vehicle configuration which is simulated in the laboratory.

Sahinkaya, Y. E.↗