Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic modelling”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Probabilistic Structural Analysis Methods (PSAM) for Select Space Propulsion System Components

Probabilistic Structural Analysis Methods (PSAM) are described for the probabilistic structural analysis of engine components for current and future space propulsion systems. Components for these systems are subjected to stochastic thermomechanical launch loads. Uncertainties or randomness also occurs in material properties, structural geometry, and boundary conditions. Material property stochasticity, such as in modulus of elasticity or yield strength, exists in every structure and is a consequence of variations in material composition and manufacturing processes. Procedures are outlined for computing the probabilistic structural response or reliability of the structural components. The response variables include static or dynamic deflections, strains, and stresses at one or several locations, natural frequencies, fatigue or creep life, etc. Sample cases illustrates how the PSAM methods and codes simulate input uncertainties and compute probabilistic response or reliability using a finite element model with probabilistic methods.

Source record↗
NASA NTRS1987

Kinetic aspects of tail dynamics - Theory and simulation

Kinetic theories relevant to the geomagnetic tail are reviewed. The topics discussed include kinetic instabilities, simulations, and current-sheet particle acceleration. Tearing mode and reconnection theories are emphasized. Kinetic treatment is appropriate for these topics since the tail plasma is collisionless. Fluid calculations are appropriate when stochastic processes dominate and for studies where long wavelengths are important. However, fluid treatments of tearing modes and reconnection require a finite resistivity in the diffusion region. Thus, although 'anomalous resistivity' can be guessed or in some cases calculated, ideally the kinetic treatment is often to be preferred. Particle motion and acceleration in the current sheet can give rise to beam-like distributions in the plasma-sheet boundary layer. Studies of current-sheet particle motion have also been used as the basis for 'kinetic' tail equilibrium models. Furthermore, quite recently current-sheet particle motion is used directly in Coroniti's explosive tail reconnection model. The 'inertial conductivity' from the equilibrium models provides the 'dissipation' necessary for reconnection.

Speiser, T. W.↗
NASA NTRS2020

Calibration Probe Uncertainty and Validation for the Hypersonic Material Environmental Test System

This paper presents an uncertainty analysis of the stagnation-point calibration probe surface predictions for conditions that span the performance envelope of the Hypersonic Materials Environmental Test System facility located at NASA Langley Research Center. A second-order stochastic expansion was constructed over 47 uncertain parameters to evaluate the sensitivities, identify the most significant uncertain variables, and quantify the uncertainty in the stagnation-point heat flux and pressure predictions of the calibration probe for a low- and high-enthalpy test condition. A sensitivity analysis showed that measurement bias uncertainty is the most significant contributor to the stagnation-point pressure and heat flux variance for the low-enthalpy condition. For the high-enthalpy condition, a paradigm shift in sensitivities revealed the computational fluid dynamics model input uncertainty as the main contributor. A comparison between the prediction and measurement of the stagnation-point conditions under uncertainty showed that there was evidence of statistical disagreement. A validation metric was proposed and applied to the prediction uncertainty to account for the statistical disagreement when compared to the possible stagnation-point heat flux and pressure measurements.

Brune, Andrew J.↗
NASA NTRS1984

Reliability and maintainability assessment factors for reliable fault-tolerant systems

A long term goal of the NASA Langley Research Center is the development of a reliability assessment methodology of sufficient power to enable the credible comparison of the stochastic attributes of one ultrareliable system design against others. This methodology, developed over a 10 year period, is a combined analytic and simulative technique. An analytic component is the Computer Aided Reliability Estimation capability, third generation, or simply CARE III. A simulative component is the Gate Logic Software Simulator capability, or GLOSS. The numerous factors that potentially have a degrading effect on system reliability and the ways in which these factors that are peculiar to highly reliable fault tolerant systems are accounted for in credible reliability assessments. Also presented are the modeling difficulties that result from their inclusion and the ways in which CARE III and GLOSS mitigate the intractability of the heretofore unworkable mathematics.

Bavuso, S. J.↗
NASA NTRS1985

Numerical simulation of a turbulent flame stabilized behind a rearward-facing step

Flow of combustible mixtures in a plane channel past a smooth contraction followed by an abrupt expansion, in a typical dump combustor configuration, is modeled by a two-dimensional numerical technique based on the random vortex method. Both the inert and the reacting case are considered. In the latter, the flame is treated as an interface, self-advancing at a prescribed normal burning speed, while the dynamic effects of expansion due to the exothermicity of combustion are expressed by volumetric source lines delineated by its front. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and the reattachment length. The stochastic turbulent velocity components manifest interesting differences, especially near the walls where three-dimensional effects of turbulence are expected to be of importance.

Hsiao, C. C.↗
NASA NTRS2020

The NANOGrav 11 yr data set: Evolution of Gravitational-wave Background Statistics

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a gravitational-wave strain spectral index of −2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new “time slicing” analysis of the 11 yr data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential “false-positive” signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and we show how independently modeling noise across 30 Fourier components, spanning NANOGrav’s frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally, we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red-noise process with low spectral index noise, and we use a tailored noise model to remove these effects.

J S Hazboun↗
NASA NTRS2016

Characterizing X-Ray and Radio Emission in the Black Hole X-Ray Binary V404 Cygni During Quiescence

We present results from multi-wavelength simultaneous X-ray and radio observations of the black hole X-ray binary V404 Cyg in quiescence. Our coverage with NuSTAR provides the very first opportunity to study the X-ray spectrum of V404 Cyg at energies above 10 keV. The unabsorbed broadband (0.3-30 keV) quiescent luminosity of the source is 8.9 x 10(exp 32) erg per sec for a distance of 2.4 kpc. The source shows clear variability on short timescales (an hour to a couple of hours) in the radio, soft X-ray, and hard X-ray bands in the form of multiple flares. The broadband X-ray spectra obtained from XMM-Newton and NuSTAR can be characterized with a power-law model having a photon index of gamma = 2.12 +/- 0.07 (90% confidence errors); however, residuals at high energies indicate spectral curvature significant at a 3 sigma confidence level with the e-folding energy of the cutoff as 20(sub -7)(sup +20) keV. Such curvature can be explained using synchrotron emission from the base of a jet outflow. Radio observations using the VLA reveal that the spectral index evolves on very fast timescales (as short as 10 minutes), switching between optically thick and thin synchrotron emission, possibly due to instabilities in the compact jet or stochastic instabilities in the accretion rate. We explore different scenarios to explain this very fast variability.

accretion disks↗
NASA NTRS1991

Supercomputer optimizations for stochastic optimal control applications

Supercomputer optimizations for a computational method of solving stochastic, multibody, dynamic programming problems are presented. The computational method is valid for a general class of optimal control problems that are nonlinear, multibody dynamical systems, perturbed by general Markov noise in continuous time, i.e., nonsmooth Gaussian as well as jump Poisson random white noise. Optimization techniques for vector multiprocessors or vectorizing supercomputers include advanced data structures, loop restructuring, loop collapsing, blocking, and compiler directives. These advanced computing techniques and superconducting hardware help alleviate Bellman's curse of dimensionality in dynamic programming computations, by permitting the solution of large multibody problems. Possible applications include lumped flight dynamics models for uncertain environments, such as large scale and background random aerospace fluctuations.

Chung, Siu-Leung↗
NASA NTRS1993

Planet formation

Models of planetary formation are developed using the present single example of a planetary system, supplemented by limited astrophysical observations of star-forming regions and circumstellar disks. The solar nebula theory and the planetesimal hypothesis are discussed. The latter is found to provide a viable theory of the growth of the terrestrial planets, the cores of the giant planets, and the smaller bodies present in the solar system. The formation of solid bodies of planetary size should be a common event, at least around young stars which do not have binary companions orbiting at planetary distances. Stochastic impacts of large bodies provide sufficient angular momentum to produce the obliquities of the planets. The masses and bulk compositions of the planets can be understood in a gross sense as resulting from planetary growth within a disk whose temperature and surface density decreased with distance from the growing sun.

Lissauer, Jack J.↗
NASA NTRS1994

The CO/SiO radiative instability in cool star atmospheres revisited

We revisit the formation of radiative instabilities in cool star atmospheres and compare our results with those given by Muchmore, Nuth, & Stencel. We have considered the combined influence of CO and SiO molecules and have computed models for a grid of effective temperatures and geometrical dilution factors for the stellar radiation. Our results are based on the analysis of the energy balance of gas elements with prescribed thermodynamic properties. Our results show that radiative instabilities are most likely primarily caused by CO, whereas SiO is expected to play only a minor role, except when the CO density is reduced compared to Local Thermodynamic Equilibrium (LTE) values or the CO band can be assumed to be optically thick. The onset of radiative instabilities is expected to be strongly modified when dynamic phenomena such as stochastic shocks are present. Our results provide strong evidence that dust formation can most likely occur via a radiative instability alone. Therefore, we present a revised version of the Muchmore et al. dust formation paradigm, which also considers hydrodynamic cooling. The new paradigm is particularly relevant in cases where dust is formed relatively close to the stellar photosphere.

Cuntz, M.↗
NASA NTRS1995

Energization in regions of CIRs unconnected to shocks are probably not the result of cross-field transport

Corotating energetic ion populations (CEIPs) associated with the forward and reverse shocks of corotating interaction regions (CIRs) are observed in CIRs at places where models say are magnetically unconnected to either shock. Such disconnections between CEIPs and shocks are common and have been documented with data from Pioneers 10 and 11 and confirmed with data from Ulysses. They pose a problem for models that account for these CEIPs in terms of ion energization at the shocks followed by ion propagation along field lines. Two possible resolutions to this problem have been suggested: diffusion of the ions across field lines and extension of the ion energization process to regions beyond the shock waves. Here we quantitatively examine the first of these possibilities. We give the Green's function solution to the convection-diffusion equation applied to idealized CIR geometry, with a source at the reverse shock -- the main producer of CEIPs. Two kinds of diffusion are considered: resonant diffusion and stochastic field line diffusion. We find that for resonant diffusion the computed ratio is many orders of magnitudes below the observed ratio. For stochastic field line diffusion, the computed ratio approximately equals the observed ratio if a diffusion coefficient appropriate to the free solar wind is used. It is several orders of magnitude below the observed ratio, however, if a diffusion coefficient appropriate to CIRs is used. We conclude that cross-field diffusion probably does not account for the presence of energetic ions in regions of CIRs that are magnetically unconnected to its shock waves. We suggest that the alternative possibility -- the energetic ions in regions magnetically unconnected to shocks result from an acceleration process that is independent of shocks -- be pursued to the point where quantitative tests can be performed.

Intriligator, Devrie S.↗
NASA NTRS1979

VTOL controls for shipboard landing

The problem of landing a VTOL aircraft on a small ship in rough seas using an automatic controller is examined. The controller design uses the linear quadratic Gaussian results of modern control theory. Linear time invariant dynamic models are developed for the aircraft, ship, and wave motions. A hover controller commands the aircraft to track position and orientation of the ship deck using only low levels of control power. Commands for this task are generated by the solution of the steady state linear quadratic gaussian regulator problem. Analytical performance and control requirement tradeoffs are obtained. A landing controller commands the aircraft from stationary hover along a smooth, low control effort trajectory, to a touchdown on a predicted crest of ship motion. The design problem is formulated and solved as an approximate finite-time linear quadratic stochastic regulator. Performance and control results are found by Monte Carlo simulations.

Mcmuldroch, C. G.↗
NASA NTRS2000

A Probabilistic Framework for the Validation and Certification of Computer Simulations

The paper presents a methodology for quantifying, propagating, and managing the uncertainty in the data required to initialize computer simulations of complex phenomena. The purpose of the methodology is to permit the quantitative assessment of a certification level to be associated with the predictions from the simulations, as well as the design of a data acquisition strategy to achieve a target level of certification. The value of a methodology that can address the above issues is obvious, specially in light of the trend in the availability of computational resources, as well as the trend in sensor technology. These two trends make it possible to probe physical phenomena both with physical sensors, as well as with complex models, at previously inconceivable levels. With these new abilities arises the need to develop the knowledge to integrate the information from sensors and computer simulations. This is achieved in the present work by tracing both activities back to a level of abstraction that highlights their commonalities, thus allowing them to be manipulated in a mathematically consistent fashion. In particular, the mathematical theory underlying computer simulations has long been associated with partial differential equations and functional analysis concepts such as Hilbert spares and orthogonal projections. By relying on a probabilistic framework for the modeling of data, a Hilbert space framework emerges that permits the modeling of coefficients in the governing equations as random variables, or equivalently, as elements in a Hilbert space. This permits the development of an approximation theory for probabilistic problems that parallels that of deterministic approximation theory. According to this formalism, the solution of the problem is identified by its projection on a basis in the Hilbert space of random variables, as opposed to more traditional techniques where the solution is approximated by its first or second-order statistics. The present representation, in addition to capturing significantly more information than the traditional approach, facilitates the linkage between different interacting stochastic systems as is typically observed in real-life situations.

Ghanem, Roger↗
NASA NTRS2009

NASA Tech Briefs, October 2009

Topics covered include: Light-Driven Polymeric Bimorph Actuators; Guaranteeing Failsafe Operation of Extended-Scene Shack-Hartmann Wavefront Sensor Algorithm; Cloud Water Content Sensor for Sounding Balloons and Small UAVs; Pixelized Device Control Actuators for Large Adaptive Optics; T-Slide Linear Actuators; G4FET Implementations of Some Logic Circuits; Electrically Variable or Programmable Nonvolatile Capacitors; System for Automated Calibration of Vector Modulators; Complementary Paired G4FETs as Voltage-Controlled NDR Device; Three MMIC Amplifiers for the 120-to-200 GHz Frequency Band; Low-Noise MMIC Amplifiers for 120 to 180 GHz; Using Ozone To Clean and Passivate Oxygen-Handling Hardware; Metal Standards for Waveguide Characterization of Materials; Two-Piece Screens for Decontaminating Granular Material; Mercuric Iodide Anticoincidence Shield for Gamma-Ray Spectrometer; Improved Method of Design for Folding Inflatable Shells; Ultra-Large Solar Sail; Cooperative Three-Robot System for Traversing Steep Slopes; Assemblies of Conformal Tanks; Microfluidic Pumps Containing Teflon[Trademark] AF Diaphragms; Transparent Conveyor of Dielectric Liquids or Particles; Multi-Cone Model for Estimating GPS Ionospheric Delays; High-Sensitivity GaN Microchemical Sensors; On the Divergence of the Velocity Vector in Real-Gas Flow; Progress Toward a Compact, Highly Stable Ion Clock; Instruments for Imaging from Far to Near; Reflectors Made from Membranes Stretched Between Beams; Integrated Risk and Knowledge Management Program -- IRKM-P; LDPC Codes with Minimum Distance Proportional to Block Size; Constructing LDPC Codes from Loop-Free Encoding Modules; MMICs with Radial Probe Transitions to Waveguides; Tests of Low-Noise MMIC Amplifier Module at 290 to 340 GHz; and Extending Newtonian Dynamics to Include Stochastic Processes.

Source record↗
NASA NTRS1974

A relation between cosmic-ray fluctuations, gradient, and diffusion coefficient

The motion of charged particles in a stochastic magnetic field is considered via a generalized quasi-linear expansion of Liouville's equation. The result is an equation relating cosmic-ray scintillations to particle gradients and to magnetic-field fluctuations (or diffusion coefficient). The resulting theory may be regarded as an example of a fluctuation-dissipation phenomenon, in which the diffusion coefficient plays the role of the dissipative parameter. The resonant interaction between particles and the random interplanetary magnetic field is considered explicitly, and it is shown that observed scintillations of high-energy (about 1 GeV) cosmic rays may be reasonably explained by the model.

Owens, A. J.↗
NASA NTRS1992

The gravitational wave contribution to cosmic microwave background anisotropies and the amplitude of mass fluctuations from COBE results

A stochastic background of primordial gravitational waves may substantially contribute, via the Sachs-Wolfe effect, to the large-scale cosmic microwave background (CMB) anisotropies recently detected by COBE. This implies a bias in any resulting determination of the primordial amplitude of density fluctuations. We consider the constraints imposed on n is less than 1 ('tilted') power-law fluctuation spectra, taking into account the contribution from both scalar and tensor waves, as predicted by power-law inflation. The gravitational wave contribution to CMB anisotropies generally reduces the required rms level of mass fluctuation, thereby increasing the linear bias parameter, even in models where the spectral index is close to the Harrison-Zel'dovich value n = 1. This 'gravitational wave bias' helps to reconcile the predictions of CDM models with observations on pairwise galaxy velocity dispersion on small scales.

Lucchin, Francesco↗
NASA NTRS2021

Astrophysics Milestones for Pulsar Timing Array Gravitational-wave Detection

The NANOGrav Collaboration reported strong Bayesian evidence for a common-spectrum stochastic process in its12.5 yr pulsar timing array data set, with median characteristic strain amplitude at periods of a year of A(yr) = 1.92(+0.75,-0.55) x 10^(-15). However, evidence for the quadrupolar Hellings & Downs interpulsar correlations, which are characteristic of gravitational-wave signals, was not yet significant. We emulate and extend the NANOGrav data set, injecting a wide range of stochastic gravitational-wave background(GWB)signals that encompass a variety of amplitudes and spectral shapes, and quantify three key milestones.(I)Given the amplitude measured in the 12.5 yr analysis and assuming this signal is a GWB, we expect to accumulate robust evidence of an interpulsar-correlated GWB signal with 15–17 yr of data, i.e., an additional 2–5 yr from the 12.5 yr data set.(II)At the initial detection, we expect a fractional uncertainty of 40% on the power-law strain spectrum slope, which is sufficient to distinguish a GWB of supermassive black hole binary origin from some models predicting more exotic origins.(III)Similarly, the measured GWB amplitude will have an uncertainty of 44% upon initial detection, allowing us to arbitrate between some population models of supermassive black hole binaries. In addition, power-law models are distinguishable from those having low-frequency spectral turnovers once 20 yr of data are reached. Even though our study is based on the NANOGrav data, we also derive relations that allow for a generalization to other pulsar timing array data sets. Most notably, by combining the data of individual arrays into the International Pulsar Timing Array, all of these milestones can be reached significantly earlier.

Nihan S. Pol↗
NASA NTRS1996

Accurate Predictions of Mean Geomagnetic Dipole Excursion and Reversal Frequencies, Mean Paleomagnetic Field Intensity, and the Radius of Earth's Core Using McLeod's Rule

The geomagnetic spatial power spectrum R(sub n)(r) is the mean square magnetic induction represented by degree n spherical harmonic coefficients of the internal scalar potential averaged over the geocentric sphere of radius r. McLeod's Rule for the magnetic field generated by Earth's core geodynamo says that the expected core surface power spectrum (R(sub nc)(c)) is inversely proportional to (2n + 1) for 1 less than n less than or equal to N(sub E). McLeod's Rule is verified by locating Earth's core with main field models of Magsat data; the estimated core radius of 3485 kn is close to the seismologic value for c of 3480 km. McLeod's Rule and similar forms are then calibrated with the model values of R(sub n) for 3 less than or = n less than or = 12. Extrapolation to the degree 1 dipole predicts the expectation value of Earth's dipole moment to be about 5.89 x 10(exp 22) Am(exp 2)rms (74.5% of the 1980 value) and the expected geomagnetic intensity to be about 35.6 (mu)T rms at Earth's surface. Archeo- and paleomagnetic field intensity data show these and related predictions to be reasonably accurate. The probability distribution chi(exp 2) with 2n+1 degrees of freedom is assigned to (2n + 1)R(sub nc)/(R(sub nc). Extending this to the dipole implies that an exceptionally weak absolute dipole moment (less than or = 20% of the 1980 value) will exist during 2.5% of geologic time. The mean duration for such major geomagnetic dipole power excursions, one quarter of which feature durable axial dipole reversal, is estimated from the modern dipole power time-scale and the statistical model of excursions. The resulting mean excursion duration of 2767 years forces us to predict an average of 9.04 excursions per million years, 2.26 axial dipole reversals per million years, and a mean reversal duration of 5533 years. Paleomagnetic data show these predictions to be quite accurate. McLeod's Rule led to accurate predictions of Earth's core radius, mean paleomagnetic field intensity, and mean geomagnetic dipole power excursion and axial dipole reversal frequencies. We conclude that McLeod's Rule helps unify geo-paleomagnetism, correctly relates theoretically predictable statistical properties of the core geodynamo to magnetic observation, and provides a priori information required for stochastic inversion of paleo-, archeo-, and/or historical geomagnetic measurements.

Voorhies, Coerte V.↗