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.

SourceResearch connectionFrom yearThrough year

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

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRS1980

Optimal regulation in systems with stochastic time sampling

An optimal control theory that accounts for stochastic variable time sampling in a distributed microprocessor based flight control system is presented. The theory is developed by using a linear process model for the airplane dynamics and the information distribution process is modeled as a variable time increment process where, at the time that information is supplied to the control effectors, the control effectors know the time of the next information update only in a stochastic sense. An optimal control problem is formulated and solved for the control law that minimizes the expected value of a quadratic cost function. The optimal cost obtained with a variable time increment Markov information update process where the control effectors know only the past information update intervals and the Markov transition mechanism is almost identical to that obtained with a known and uniform information update interval.

Montgomery, R. C.↗

NASA NTRS1991

Ion pickup, scattering, and stochastic acceleration in the cometary environment of P/Giacobini-Zinner

Observations and theory related to the scattering and acceleration of cometary pickup ions are reviewed with emphasis on Comet P/Giacobini-Zinner. A comparison of the regions upstream and downstream of the bow shock is made to assess the relative merits of each as a site for stochastic acceleration of ions above the pickup energy through interaction with low-frequency MHD waves. In the far upstream region the data are most consistent with a model where pickup ions generate a low level of MHD waves but remain relatively scatter-free. In the downstream region intense magnetic fluctuations gives rise to rapid isotropization of the ions and a second-order stochastic acceleration. The properties of the MHD power spectrum are related to the energetic ion spectrum in the framework of a leaky box model where the bulk of the acceleration occurs downstream of the shock throughout the cometosheath. Good agreement of the observations with theory is evident for both P/Giacobini-Zinner and P/Halley.

Barbosa, D. D.↗

NASA NTRS2014

Learning the Task Management Space of an Aircraft Approach Model

Validating models of airspace operations is a particular challenge. These models are often aimed at finding and exploring safety violations, and aim to be accurate representations of real-world behavior. However, the rules governing the behavior are quite complex: nonlinear physics, operational modes, human behavior, and stochastic environmental concerns all determine the responses of the system. In this paper, we present a study on aircraft runway approaches as modeled in Georgia Tech's Work Models that Compute (WMC) simulation. We use a new learner, Genetic-Active Learning for Search-Based Software Engineering (GALE) to discover the Pareto frontiers defined by cognitive structures. These cognitive structures organize the prioritization and assignment of tasks of each pilot during approaches. We discuss the benefits of our approach, and also discuss future work necessary to enable uncertainty quantification.

Validation↗

NASA NTRS1991

Hybrid state-space self-tuning control using dual-rate sampling

This paper presents a hybrid state-space self-tuning control scheme using dual-rate sampling for suboptimal digital adaptive control of linear time-invariant continuous-time multivariable stochastic systems with unknown parameters. An equivalent fast-rate discrete-time state-space innovation model (with estimated states) of the continuous-time system is constructed by using the estimated system parameters and Kalman gain. To utilize the existing optimal regional-pole assignment method developed in the continuous-time domain, the constructed fast-rate discrete-time model is converted into an equivalent continuous-time model for the development of a state-feedback optimal control law with pole placement in a specific region. The developed analog optimal control law is then converted into an equivalent pseudo-slow-rate digital control law via the proposed digital redesign technique, which can be realized via slow-rate digital electronics. The proposed method enables the development of a digitally implementable advanced control algorithm for digital adaptive control of continuous-time multivariable stochastic systems which may be unstable and/or have nonminimum phase.

Shieh, Leang S.↗

NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗

NASA NTRS1979

On reliable control system designs with and without feedback reconfigurations

This paper contains an overview of a theoretical framework for the design of reliable multivariable control systems, with special emphasis on actuator failures and necessary actuator redundancy levels. Using a linear model of the system, with Markovian failure probabilities and quadratic performance index, an optimal stochastic control problem is posed and solved. The solution requires the iteration of a set of highly coupled Riccati-like matrix difference equations; if these converge one has a reliable design; if they diverge, the design is unreliable, and the system design cannot be stabilized. In addition, it is shown that the existence of a stabilizing constant feedback gain and the reliability of its implementation is equivalent to the convergence properties of a set of coupled Riccati-like matrix difference equations. In summary, these results can be used for offline studies relating the open loop dynamics, required performance, actuator mean time to failure, and functional or identical actuator redundancy, with and without feedback gain reconfiguration strategies.

Birdwell, J. D.↗

NASA NTRS1985

The structure of dilute combusting sprays

An experimental and theoretical study of drop processes in a turbulent flame is described. The experiments involved a monodisperse (105 and 180 micro m initial diameter) stream of methanol drops injected at the base of a turbulent methane-fueled diffusion flame burning in still air. The following measurements were made: mean and fluctuating phase velocities, mean drop number flux, drop-size distributions and mean gas-phase temperatures. Measurements were compared with predictions of two separated flow models: (1) deterministic separated flow, where drop-turbulence interactions are ignored; and (2) stochastic separated flow, where drop-turbulence interactions are considered using random-walk computations. The stochastic separated flow analysis yielded best agreement with measurements, since it provides for turbulent dispersion of drops which was important for present test conditions (and probably for most combusting sprays as well). Distinguishing the presence or absence of envelope flames around the drops, however, was relatively unimportant for present test conditions, since the drops spent most of their lifetime in fuel-rich regions of the flow where this distinction is irrelevant.

Shuen, J. S.↗

NASA NTRS1986

Ion heating by waves with frequencies below the ion gyrofrequency

Ion heating by broadband and coherent waves with relatively small amplitude and frequencies below the ion gyrofrequency is examined. The theoretical results, verified by test particle calculations, show that such waves can heat ions to high energies, and in particular produce the ion conic distributions observed in the earth's magnetosphere. The model utilized differs from other approaches in that neither quasi-linear theory nor stochastic acceleration in a large amplitude wave are used to heat ions. This heating mechanism may be significant in other geophysical and astrophysical environments.

Temerin, M.↗

NASA NTRS1988

Particle-laden weakly swirling free jets: Measurements and predictions

A theoretical and experimental investigation of particle-laden, weakly swirling, turbulent free jets was conducted. Glass particles, having a Sauter mean diameter of 39 microns with a standard deviation of 15 microns, were used. A single loading ratio of 0.2 was used in the experiments. Measurements are reported for three swirl numbers, ranging from 0.0 to 0.3. The measurements included mean and fluctuating velocities of both phases, and particle mass flux distributions. Measurements were compared with predictions from three types of multiphase flow analysis: locally homogeneous flow (LHF); deterministic separated flow (DSF); and stochastic separated flow (SSF). For the particle-laden jets, the LHF and DSF models did not provide very satisfactory predictions. The LHF model generally overestimated the rate of decay of particle mean axial and angular velocities with streamwise distance, due to the neglect of particle inertia. The LHF model predictions of particle mass flux also showed poor agreement with measurements due to the assumption of no-slip between phases. The DSF model also performed quite poorly for predictions of particle mass flux, because turbulent dispersion of the particles was neglected. The SSF model, which accounts for both particle inertia and turbulent dispersion of the particles, yielded reasonably good predictions throughout the flow field for the particle-laden jets.

Bulzan, Daniel L.↗

NASA NTRS1988

Particle-laden weakly swirling free jets - Measurements and predictions

A theoretical and experimental investigation of particle-laden, weakly swirling, turbulent free jets was conducted. Glass particles, having a Sauter mean diameter of 39 microns with a standard deviation of 15 microns, were used. A single loading ratio of 0.2 was used in the experiments. Measurements are reported for three swirl numbers, ranging from 0.0 to 0.3. The measurements included mean and fluctuating velocities of both phases, and particle mass flux distributions. Measurements were compared with predictions from three types of multiphase flow analysis: locally homogeneous flow (LHF); deterministic separated flow (DSF); and stochastic separated flow (SSF). For the particle-laden jets, the LHF and DSF models did not provide very satisfactory predictions. The LHF model generally overestimated the rate of decay of particle mean axial and angular velocities with streamwise distance, due to the neglect of particle inertia. The LHF model predictions of particle mass flux also showed poor agreement with measurements due to the assumption of no-slip between phases. The DSF model also performed quite poorly for predictions of particle mass flux, because turbulent dispersion of the particles was neglected. The SSF model, which accounts for both particle inertia and turbulent disperison of the particles, yielded reasonably good predictions throughout the flow field for the particle-laden jets.

Bulzan, Daniel L.↗

NASA NTRS1989

Mechanical model of an earthquake fault

The dynamic behavior of a simple mechanical model of an earthquake fault is studied. This model, introduced originally by Burridge and Knopoff (1967), consists of an elastically coupled chain of masses in contact with a moving rough surface. The present version of the model retains the full Newtonian dynamics with inertial effects and contains no externally imposed stochasticity or spatial inhomogeneity. The only nonlinear feature is a velocity-weakening stick-slip friction force between the masses and the moving surface. This system is being driven persistently toward a slipping instability and therefore exhibits noisy sequences of earthquakelike events. These events are observed in numerical simulations, and many of their features can be predicted analytically.

Carlson, J. M.↗

NASA NTRS1990

An experimental and numerical study of particle-laden coaxial jet flows

A detailed experimental and numerical study of the developing region of coaxial jet flows with and without glass beads is performed. A two-component phase/Doppler interferometer is used to measure mean and fluctuating velocity components for each phase and particle-number density. The numerical calculation is based on a stochastic Lagrangian treatment for the particles and a recently proposed two-equation turbulence model for two-phase flows. Results show that the particle-number density profile becomes narrower than the corresponding profile for round jet flow and that the particles attain a uniform velocity across the jet radius. The particles attenuate the level of gas turbulence and increase their anisotropy level. The numerical calculations yield reasonable and encouraging agreement with the measurements.

Mostafa, A. A.↗

NASA NTRS1995

Simulation of rotor blade element turbulence

A piloted, motion-based simulation of Sikorsky's Black Hawk helicopter was used as a platform for the investigation of rotorcraft responses to vertical turbulence. By using an innovative temporal and geometrical distribution algorithm that preserved the statistical characteristics of the turbulence over the rotor disc, stochastic velocity components were applied at each of twenty blade-element stations. This model was implemented on NASA Ames' Vertical Motion Simulator (VMS), and ten test pilots were used to establish that the model created realistic cues. The objectives of this research included the establishment of a simulation-technology basis for future investigation into real-time turbulence modeling. This goal was achieved; our extensive additions to the rotor model added less than a 10 percent computational overhead. Using a VAX 9000 computer the entire simulation required a cycle time of less than 12 msec. Pilot opinion during this simulation was generally quite favorable. For low speed flight the consensus was that SORBET (acronym for title) was better than the conventional body-fixed model, which was used for comparison purposes, and was determined to be too violent (like a washboard). For high speed flight the pilots could not identify differences between these models. These opinions were something of a surprise because only the vertical turbulence component on the rotor system was implemented in SORBET. Because of the finite-element distribution of the inputs, induced outputs were observed in all translational and rotational axes. Extensive post-simulation spectral analyses of the SORBET model suggest that proper rotorcraft turbulence modeling requires that vertical atmospheric disturbances not be superimposed at the vehicle center of gravity but, rather, be input into the rotor system, where the rotor-to-body transfer function severely attenuates high frequency rotorcraft responses.

Mcfarland, R. E.↗

NASA NTRS2007

Stochastic Characterization of Flutter using Historical Wind Tunnel Data

Methods for predicting the onset of flutter during an experiment are traditionally applied treating the data as deterministic values. Uncertainty and variation in the data is often glossed over by using best-fit curves to represent the information. This paper applies stochastic treatments to wind tunnel data obtained for the Piezoelectric Aeroelastic Response Tailoring Investigation model. These methods include modal amplitude tracking, modal frequency tracking and several applications of the flutter margin method. The flutter margin method was developed by Zimmerman and Weissenburger, and extended by Poirel, Dunn and Porter to incorporate uncertainty. Much of the current work follows the future work recommendations of Poirel, Dunn and Porter.

Heeg, Jennifer↗

NASA NTRSDate not supplied

Simulation of etch pit formation through active sites in carbon fiber micro-structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation as a result of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gasses such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

NASA NTRSDate not supplied

Simulation of Etch Pit Formation in DSMC Through Active Sites in Carbon Fiber Micro-Structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation because of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gases such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various detailed surface reaction mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

NASA NTRSDate not supplied

DSMC Simulation of Etch Pit Formation Through Active Sites in Carbon Fiber Micro-Structures

DSMC↗

NASA NTRS2009

Spreadsheets in Team X: Preserving Order in an Inherently Chaotic Environment

JPL is NASA's prime center for deep space missions. In response to the need to reduce the cost and time to complete early concept studies and proposals JPL created the first concurrent engineering team in the aerospace industry: Team X. Started in 1995, Team X has carried out over 800 studies, dramatically reducing the time and cost involved, and has been the model for other concurrent engineering teams both within NASA and throughout the larger aerospace community. Since its inception, the software backbone of this highly successful design team - engaged in examining some of NASA's cutting edge concepts - has been the unassuming spreadsheet. Over the years the Team X spreadsheet-based tools have evolved from simple standalone engineering models into a networked spreadsheet intensive system with real time parameter updating. Recent new capabilities include stochastic cost estimation and a graphical drag and drop block diagram that automatically populates the related spreadsheet parameters of cost, mass and power. This paper describes how the spreadsheet functions within Team X: its history, architecture, current capabilities, enabling strengths and persistent weaknesses. In addition, the verification methods and institutional oversight that have evolved as the Team X products became increasingly critical to Laboratory success are also discussed.

concurrent engineering↗