Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo 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 433 records · Page 24

Analysis of a semiclassical model for rotational transition probabilities

A semiclassical model proposed by Pearson and Hansen (1974) for computing collision-induced transition probabilities in diatomic molecules is tested by the direct-simulation Monte Carlo method. Specifically, this model is described by point centers of repulsion for collision dynamics, and the resulting classical trajectories are used in conjunction with the Schroedinger equation for a rigid-rotator harmonic oscillator to compute the rotational energy transition probabilities necessary to evaluate the rotation-translation exchange phenomena. It is assumed that a single, average energy spacing exists between the initial state and possible final states for a given collision.

Deiwert, G. S.↗

A numerical model study of long-term planetary wave predictability

A numerical model was developed to study planetary wave behavior over periods up to 4 months under the influences of stochastic and diffusive forcing by synoptic scale motions and of sea surface temperature anomalies. The model is suitable for Monte Carlo experiments which provide ensemble mean and variance predictions and permit the response to heating anomalies to be partly separated from random variability. The model incorporates the linear-balanced dynamical equations for wind at two tropospheric levels and temperature at one level. The model atmosphere is forced by semiempirical January heating linearly dependent on local air temperature. Topography and the lower boundary layer are represented by vertical motion and drag functions, respectively, of the lower level wind. The linear phase speeds of Rossby waves are modified to reflect the effect of the stratosphere.

Wobus, R. L.↗

Analytic Guided-Search Model of Human Performance Accuracy in Target- Localization Search Tasks

Current models of human visual search have extended the traditional serial/parallel search dichotomy. Two successful models for predicting human visual search are the Guided Search model and the Signal Detection Theory model. Although these models are inherently different, it has been difficult to compare them because the Guided Search model is designed to predict response time, while Signal Detection Theory models are designed to predict performance accuracy. Moreover, current implementations of the Guided Search model require the use of Monte-Carlo simulations, a method that makes fitting the model's performance quantitatively to human data more computationally time consuming. We have extended the Guided Search model to predict human accuracy in target-localization search tasks. We have also developed analytic expressions that simplify simulation of the model to the evaluation of a small set of equations using only three free parameters. This new implementation and extension of the Guided Search model will enable direct quantitative comparisons with human performance in target-localization search experiments and with the predictions of Signal Detection Theory and other search accuracy models.

Eckstein, Miguel P.↗

Modeling and Control of Airport Queueing Dynamics under Severe Flow Restrictions

Based on field observations and interviews with controllers at BOS and EWR, we identify the closure of local departure fixes as the most severe class of airport departure restrictions. A set of simple queueing dynamics and traffic rules are developed to model departure traffic under such restrictions. The validity of the proposed model is tested via Monte Carlo simulation against 10 hours of actual operations data collected during a case-study at EWR on June 29,2000. In general, the model successfully reproduces the aggregate departure congestion. An analysis of the average error over 40 simulation runs indicates that flow-rate restrictions also significantly impact departure traffic; work is underway to capture these effects. Several applications and what-if scenarios are discussed for future evaluation using the calibrated model.

Carr, Francis↗

A simulation model for probabilistic analysis of Space Shuttle abort modes

A simulation model which was developed to provide a probabilistic analysis tool to study the various space transportation system abort mode situations is presented. The simulation model is based on Monte Carlo simulation of an event-tree diagram which accounts for events during the space transportation system's ascent and its abort modes. The simulation model considers just the propulsion elements of the shuttle system (i.e., external tank, main engines, and solid boosters). The model was developed to provide a better understanding of the probability of occurrence and successful completion of abort modes during the vehicle's ascent. The results of the simulation runs discussed are for demonstration purposes only, they are not official NASA probability estimates.

Hage, R. T.↗

A Monte Carlo approach to competition strategy

Variables taken into account in glider flight strategy decisions are modeled. These include height of clouds, distance between thermals, time of day, water ballast, present altitude, weather changes, lift organization, and distance to goal, as well as the strength of the next thermal. Results of the Monte Carlo atmospheric model are discussed.

Teter, M. P.↗

Advanced data assimilation in strongly nonlinear dynamical systems

Advanced data assimilation methods are applied to simple but highly nonlinear problems. The dynamical systems studied here are the stochastically forced double well and the Lorenz model. In both systems, linear approximation of the dynamics about the critical points near which regime transitions occur is not always sufficient to track their occurrence or nonoccurrence. Straightforward application of the extended Kalman filter yields mixed results. The ability of the extended Kalman filter to track transitions of the double-well system from one stable critical point to the other depends on the frequency and accuracy of the observations relative to the mean-square amplitude of the stochastic forcing. The ability of the filter to track the chaotic trajectories of the Lorenz model is limited to short times, as is the ability of strong-constraint variational methods. Examples are given to illustrate the difficulties involved, and qualitative explanations for these difficulties are provided. Three generalizations of the extended Kalman filter are described. The first is based on inspection of the innovation sequence, that is, the successive differences between observations and forecasts; it works very well for the double-well problem. The second, an extension to fourth-order moments, yields excellent results for the Lorenz model but will be unwieldy when applied to models with high-dimensional state spaces. A third, more practical method--based on an empirical statistical model derived from a Monte Carlo simulation--is formulated, and shown to work very well. Weak-constraint methods can be made to perform satisfactorily in the context of these simple models, but such methods do not seem to generalize easily to practical models of the atmosphere and ocean. In particular, it is shown that the equations derived in the weak variational formulation are difficult to solve conveniently for large systems.

Miller, Robert N.↗

Modeling and Simulation Techniques for the NASA SLS Service Module Panel Separation Event; from Loosely-Coupled Euler to Fully-Coupled 6-DOF, Time-Accurate, Navier-Stokes Methodologies

An aerodynamic database has been generated for use by the Orion Multi-Purpose Crew Vehicle (MPCV) Program to analyze Service Module (SM) panel jettison from the NASA SLS vehicle. The database is a combination of CFD data for the panel aerodynamic coefficients, and MATLAB code written to query the CFD data. The Cart3D inviscid CFD flow solver was used to generate the panel aerodynamic coefficients for static panel orientations and free stream conditions that can occur during the jettison event. The MATLAB code performs the multivariate interpolation to obtain aerodynamic coefficients. The MATLAB code uses input for SM panel parameters and returns the SM panel aerodynamic force and moment coefficients for use with a Six-Degree-of-Freedom (6-DOF) motion solver to model the jettison event. This paper examines the accuracy of the sequential-static database approach by modeling the panel jettison event with a fully-coupled, time-dependent, viscous, moving-body CFD simulation. The fully-coupled simulation is obtained using the Loci/Chem unstructured Navier-Stokes CFD solver. The results show that the fully-coupled approach agrees well with the loosely-coupled database/6-DOF approach, indicating that unsteady effects are minimal for the panel jettison event. These results suggest that the database/6-DOF approach is sufficient. In addition, this paper presents the development of an uncertainty model for use in Monte Carlo analysis of the panel jettison event. Here viscous CFD simulations are obtained with Loci/Chem and compared to the inviscid CFD forces and moments. An uncertainty model based on model-form error and numerical error is presented.

Hall, Leslie H.↗

Calculation of astrophysical spallation reactions using the RENO model

The RENO model for the Monte-Carlo treatment of astrophysical spallation reactions has been used to generate preliminary cross-sections for the purpose of illustrating the discrete-nucleon approach to spallation modeling and to exhibit differences between two versions of RENO. Comparisons with experimental, theoretical, and semiempirical data demonstrate the practicability of the discrete-nucleon approach.-

Ayres, C. L.↗

Highlights of Transient Plume Impingement Model Validation and Applications

This paper describes highlights of an ongoing validation effort conducted to assess the viability of applying a set of analytic point source transient free molecule equations to model behavior ranging from molecular effusion to rocket plumes. The validation effort includes encouraging comparisons to both steady and transient studies involving experimental data and direct simulation Monte Carlo results. Finally, this model is applied to describe features of two exotic transient scenarios involving NASA Goddard Space Flight Center satellite programs.

Woronowicz, Michael↗

Development of a Spectral Model Based on Charge Transport for the Swift/BAT 32K CdZnTe Detector Array

The properties of 32K CdZnTe (4 x 4 sq mm large, 2 mm thick) detectors have been studied in the pre-flight calibration of the Burst Alert Telescope (BAT) on-board the Swift Gamma-ray Burst Explorer (scheduled for launch in November 2004). In order to understand the energy response of the BAT CdZnTe array, we first quantify the mobility-lifetime (mu tau) products of carriers in individual CdZnTe detectors, which produce a position dependency in the charge induction efficiency and results in a low energy tail in the energy spectrum. Based on a new method utilizing (57)Co spectra obtained at different bias voltages, the mu tau for electrons ranges from 5.0 x 10(exp -4) to 1.0 x 10(exp -2) sq cm/V while the mu tau for holes ranges from 1.3 x 10(exp -5 to 1.8 x 10(exp -4) sq cm/V. We find that this wide distribution of mu tau products explains the large diversity in spectral shapes between CdZnTe detectors well. We also find that the variation of mu tau products can be attributed to the difference of crystal ingots or manufacturing harness. We utilize the 32K sets of extracted mu tau products to develop a spectral model of the detector. In combination with Monte Carlo simulations, we can construct a spectral model for any photon energy or any incident angle.

Sato, Goro↗

Monte Carlo approach in assessing damage in higher order structures of DNA

We have developed a computer monitor of nuclear DNA in the form of chromatin fibre. The fibres are modeled as a ideal solenoid consisting of twenty helical turns with six nucleosomes per turn. The chromatin model, in combination with are Monte Carlo theory of radiation damage induces by charged particles, based on general features of tack structure and stopping power theory, has been used to evaluate the influence of DNA structure on initial damage. An interesting has emerged from our calculations. Our calculated results predict the existence of strong spatial correlations in damage sites associated with the symmetries in the solenoidal model. We have calculated spectra of short fragments of double stranded DNA produced by multiple double strand breaks induced by both high and low LET radiation. The spectra exhibit peaks at multiples of approximately 85 base pairs (the nucleosome periodicity), and approximately 1000 base pairs (solenoid periodicity). Preliminary experiments to investigate the fragment distributions from irradiated DNA, made by B. Rydberg at Lawrence Berkeley Laboratory, confirm the existence of short DNA fragments and are in substantial agreement with the predictions of our theory.

Non-NASA Center↗

3-D Simulations of NSTAR Ion Thruster Plasma Environment

Full three-dimensional numerical models based on partical-in-cell with Monte Carlo collisions (PIC-MCC) simulations are developed to study ion thruster plasma interactions. 3-D simulation results for the NSTAR ion thruster plasma environment are presented.

NSTAR↗

Cometary atmospheres: Modeling the spatial distribution of observed neutral radicals

An algorithm for the random walk problem of multiple elastic collisions between newly formed non-thermal neutral cometary radicals and the outflowing cometary molecules was incorporated into the Monte Carlo particle-trajectory model. Preliminary model analysis has shown that the effects of collision on the observed spatial distribution of cometary radicals becomes important for the larger bright comets, especially at smaller values of the helicocentric distance. The model and early results are discussed herein.

Combi, M. R.↗

SUMS experiment flight results on STS-35

Calibrated pressure measurements for species with mass to charge ratios up to 50 amu/e(-) were obtained from the Shuttle Upper Atmosphere Mass Spectrometer (SUMS) experiment during reentry on the STS-35 mission. Data were collected from 180 km, when the signal rose above the background, to about 87 km, when the SUMS system automatically closed the gas inlet value. However, data above 115 km was contaminated from a source of gas emanating from pressure transducers connected in parallel to the mass spectrometer. At lower altitudes, the pressure transducer data is compared with the mass spectrometer total pressure with excellent agreement. The free-stream density in the rarefied flow flight regime is calculated using an orifice pressure coefficient model based upon direct simulation Monte Carlo results. This density, when compared with the 1976 U.S. standard atmosphere model, exhibits the wave-like nature seen on previous flights using accelerometry. In addition, selected spectra are presented at higher altitudes (320 km) showing the effects of the ingestion of gases from a forward fuselage fuel dump. An analysis of the spectra data from this event is presented to show that no significant permanent changes occurred which affected the data interpretation at lower altitudes. Further, the localized chemistry from the individual species during the onset of aerodynamic heating is examined to the extent possible for a closed source system, such as SUMS. Near the orifice entrance, a significant amount of CO2 was generated from chemical reactions with the carbon panels of the Orbiter and absorbed oxygen on the system tubing.

Blanchard, Robert C.↗

Shuttle Upper Atmosphere Mass Spectrometer Experimental Flight Results

Calibrated pressure measurements for species with mass-to-charge ratios up to 50 amu/e(-) were obtained trom the shuttle upper atmosphere mass spectrometer experiment during re-entry on the STS-35 mission. The principal experimental objective is to obtain measurements of freestream density in the hypersonic rarefied flow flight regime. Data were collected from 180 to about 87 km. However, data above 115 km were contaminated from a source of gas emanating from pressure transdueers connected in parallel to the mass spectrometer. At lower altitudes, the pressure transducer data are compared to the mass spectrometer total pressure with excellent agreement. Near the orifice entrance, a significant amount of CO2 was generated from chemical reactions. The freestream density in the rarefied flow flight regime is calculated using an orifice pressure coefficient model based upon direct simulation Monte Carlo results. This density, when compared with the 1976 U.S. Standard Atmosphere model, exhibits the wavelike nature seen on previous flights using accelerometry. Selected spectra are presented at higher altitudes (320 km) showing the effects of the ingestion of gases from a forward fuselage fuel dump.

Blanchard, R. C.↗

Evaluation of Fast-Time Wake Vortex Models using Wake Encounter Flight Test Data

This paper describes a methodology for the integration and evaluation of fast-time wake models with flight data. The National Aeronautics and Space Administration conducted detailed flight tests in 1995 and 1997 under the Aircraft Vortex Spacing System Program to characterize wake vortex decay and wake encounter dynamics. In this study, data collected during Flight 705 were used to evaluate NASA's fast-time wake transport and decay models. Deterministic and Monte-Carlo simulations were conducted to define wake hazard bounds behind the wake generator. The methodology described in this paper can be used for further validation of fast-time wake models using en-route flight data, and for determining wake turbulence constraints in the design of air traffic management concepts.

Ahmad, Nashat N.↗

A model of the 8-25 micron point source infrared sky

A model is developed that permits the representation of the Galactic disk and bulge, the molecular ring, extragalactic sky, and other Galactic structures in terms of the IR point-source sky. Full characterization is achieved for up to 87 types of Galactic sources in terms of scale heights and absolute magnitudes, and the model is guided by a parallel Monte Carlo simulation of the Galaxy at 12 microns. The model predicts differential and cumulative IR source counts for all bandpasses in the IRAS ranges and can therefore be employed to predict the character of the point-source sky related to IR space observations.

Wainscoat, Richard J.↗