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 1,009 records · Page 56

Modeling human tracking error in several different anti-tank systems

An optimal control model for generating time histories of human tracking errors in antitank systems is outlined. Monte Carlo simulations of human operator responses for three Army antitank systems are compared. System/manipulator dependent data comparisons reflecting human operator limitations in perceiving displayed quantities and executing intended control motions are presented. Motor noise parameters are also discussed.

Kleinman, D. L.↗

A spatial model of wind shear and turbulence for flight simulation

A three dimensional model which combines measurements of wind shear in the real atmosphere with three dimensional Monte Carlo simulated turbulence was developed. The wind field over the body of an aircraft can be simulated and all aerodynamic loads and moments calculated.

Campbell, C. W.↗

Statistical mechanics of light elements at high pressure. VII - A perturbative free energy for arbitrary mixtures of H and He

A model free energy is presented which accurately represents results from 45 high-precision Monte Carlo calculations of the thermodynamics of hydrogen-helium mixtures at pressures of astrophysical and planetophysical interest. The free energy is calculated using free-electron perturbation theory (dielectric function theory), and is an extension of the expression given in an earlier paper in this series. However, it fits the Monte Carlo results more accurately, and is valid for the full range of compositions from pure hydrogen to pure helium. Using the new free energy, the phase diagram of mixtures of liquid metallic hydrogen and helium is calculated and compared with earlier results. Sample results for mixing volumes are also presented, and the new free energy expression is used to compute a theoretical Jovian adiabat and compare the adiabat with results from three-dimensional Thomas-Fermi-Dirac theory. The present theory gives slightly higher densities at pressures of about 10 megabars.

Hubbard, W. B.↗

Rule-based analysis of pilot decisions

The application of the rule identification technique to the analysis of human performance data is proposed. The relation between the language and identifiable consistencies is discussed. The advantages of production system models for the description of complex human behavior are studied. The use of a Monte Carlo significance testing procedure to assure the validity of the rule identification is examined. An example of the rule-based analysis of Palmer's (1983) data is presented.

Lewis, C. M.↗

Simulation and analyses of the aeroassist flight experiment attitude update method

A method which will be used to update the alignment of the Aeroassist Flight Experiment's Inertial Measuring Unit is simulated and analyzed. This method, the Star Line Maneuver, uses measurements from the Space Shuttle Orbiter star trackers along with an extended Kalman filter to estimate a correction to the attitude quaternion maintained by an Inertial Measuring Unit in the Orbiter's payload bay. This quaternion is corrupted by on-orbit bending of the Orbiter payload bay with respect to the Orbiter navigation base, which is incorporated into the payload quaternion when it is initialized via a direct transfer of the Orbiter attitude state. The method of updating this quaternion is examined through verification of baseline cases and Monte Carlo analysis using a simplified simulation, The simulation uses nominal state dynamics and measurement models from the Kalman filter as its real world models, and is programmed on Microvax minicomputer using Matlab, and interactive matrix analysis tool. Results are presented which confirm and augment previous performance studies, thereby enhancing confidence in the Star Line Maneuver design methodology.

Carpenter, J. R.↗

Shock interference prediction using direct simulation Monte Carlo

The shock interaction produced when an incident shock impinges on an inlet cowl lip is investigated for a cowl lip radius of 0.1 in. and flight conditions of Mach 15 and 35-km altitude. This problem is of interest in the design of hypersonic flight vehicles because the interference heating at the cowl lip surface from such a shock interaction is expected to be very high and the design of adequate thermal protection in this region will be challenging. Noncontinuum effects may be significant at this combination of altitude and small dimensions. If this is so, a computational solution based on the Navier-Stokes equations could significantly overpredict the surface heating on the cowl lip. Therefore, a kinetic theory solution technique, direct simulation Monte Carlo (DSMC), is used to predict the flowfield and surface heating. The details of modeling this problem and some preliminary results are presented. The same qualitative flow patterns are seen with DSMC and a continuum solution of the shock interference case, but the peak heating predicted with DSMC is somewhat lower.

Carlson, Ann B.↗

Aerodynamics of Shuttle Orbiter at high altitudes

The high-altitude/high-Knudsen number aerodynamics of the Shuttle Orbiter are computed from Low-Earth Orbit down to 100 km using three-dimensional direct simulation Monte Carlo and free molecule codes. Results are compared with Blanchard's latest Shuttle aerodynamic model, which is based on in-flight accelerometer measurements, and bridging formula models. Good comparison is observed, except for the normal force and pitching moment coefficients. The present results were obtained for a generic Shuttle geometry configuration corresponding to a zero deflection for all control surfaces.

Rault, Didier F. G.↗

Development of a Probabilistic Dynamic Synthesis Method for the Analysis of Nondeterministic Structures

Accounting for the statistical geometric and material variability of structures in analysis has been a topic of considerable research for the last 30 years. The determination of quantifiable measures of statistical probability of a desired response variable, such as natural frequency, maximum displacement, or stress, to replace experience-based "safety factors" has been a primary goal of these studies. There are, however, several problems associated with their satisfactory application to realistic structures, such as bladed disks in turbomachinery. These include the accurate definition of the input random variables (rv's), the large size of the finite element models frequently used to simulate these structures, which makes even a single deterministic analysis expensive, and accurate generation of the cumulative distribution function (CDF) necessary to obtain the probability of the desired response variables. The research presented here applies a methodology called probabilistic dynamic synthesis (PDS) to solve these problems. The PDS method uses dynamic characteristics of substructures measured from modal test as the input rv's, rather than "primitive" rv's such as material or geometric uncertainties. These dynamic characteristics, which are the free-free eigenvalues, eigenvectors, and residual flexibility (RF), are readily measured and for many substructures, a reasonable sample set of these measurements can be obtained. The statistics for these rv's accurately account for the entire random character of the substructure. Using the RF method of component mode synthesis, these dynamic characteristics are used to generate reduced-size sample models of the substructures, which are then coupled to form system models. These sample models are used to obtain the CDF of the response variable by either applying Monte Carlo simulation or by generating data points for use in the response surface reliability method, which can perform the probabilistic analysis with an order of magnitude less computational effort. Both free- and forced-response analyses have been performed, and the results indicate that, while there is considerable room for improvement, the method produces usable and more representative solutions for the design of realistic structures with a substantial savings in computer time.

Brown, A. M.↗

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗

Filtered Mass Density Function for Design Simulation of High Speed Airbreathing Propulsion Systems

The objective of this research is to develop and implement new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. We have just completed two (2) years of Phase I of this research. This annual report provides a brief and up-to-date summary of our activities during the period: September 1, 2000 through August 31, 2001. In the work within the past year, a methodology termed "velocity-scalar filtered density function" (VSFDF) is developed and implemented for large eddy simulation (LES) of turbulent flows. In this methodology the effects of the unresolved subgrid scales (SGS) are taken into account by considering the joint probability density function (PDF) of all of the components of the velocity and scalar vectors. An exact transport equation is derived for the VSFDF in which the effects of the unresolved SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source terms appear in closed form. The remaining unclosed terms in this equation are modeled. A system of stochastic differential equations (SDEs) which yields statistically equivalent results to the modeled VSFDF transport equation is constructed. These SDEs are solved numerically by a Lagrangian Monte Carlo procedure. The consistency of the proposed SDEs and the convergence of the Monte Carlo solution are assessed by comparison with results obtained by an Eulerian LES procedure in which the corresponding transport equations for the first two SGS moments are solved. The unclosed SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source in the Eulerian LES are replaced by corresponding terms from VSFDF equation. The consistency of the results is then analyzed for a case of two dimensional mixing layer.

Drozda, T. G.↗

Remote-Sensing Reflectance and Inherent Optical Properties for Optically Deep Waters: A Revisit

Remote-sensing reflectance (r(rs)) is defined as the ratio of upwelling radiance to downwelling irradiance. Relationships between remote-sensing reflectance and inherent optical properties serve as the basis for ocean-color modeling, as well as for spectral deduction of oceanic constituents through analytical/semi-analytical models of ocean color. A decade ago, a simple and concise formula based on Monte Carlo simulations was developed by relating rrs to a property u, the ratio of backscattering (b(b)) to the sum of absorption (a) and backscattering (u = b(b)/(a+b(b))). This relationship generally ignored the shape differences in phase functions between molecular scattering and particle scattering. In this study, the relationship is updated with separate parameters for molecular and particle scattering, based on the Radiative Transfer Equation through use of Hydrolight numerical solutions. The new approach fits r(rs) better than an earlier traditional formula, for both clear and turbid waters.

Lee, Zhong-Ping↗

Periodicity in extinction and the problem of catastrophism in the history of life

The hypothesis that extinction events have recurred periodically over the last quarter billion years is greatly strengthened by new data on the stratigraphic ranges of marine animal genera. In the interval from the Permian to Recent, these data encompass some 13,000 generic extinctions, providing a more sensitive indicator of species-level extinctions than previously used familial data. Extinction time series computed from the generic data display nine strong peaks that are nearly uniformly spaced at 26 Ma intervals over the last 270 Ma. Most of these peaks correspond to extinction events recognized in more detailed, if limited, biostratigraphic studies. These new data weaken or negate most arguments against periodicity, which have involved criticisms of the taxonomic data base, sampling intervals, chronometric time scales, and statistical methods used in previous analyses. The criticisms are reviewed in some detail and various new calculations and simulations, including one assessing the effects of paraphyletic taxa, are presented. Although the new data strengthen the case for periodicity, they offer little new insight into the deriving mechanism behind the pattern. However, they do suggest that many of the periodic events may not have been catastrophic, occurring instead over several stratigraphic stages or substages.

NASA Discipline Number 52-40↗

Explore Full Range of QSO/AGN Properties

The goal of the proposal is to perform ISO spectroscopic studies, including data analysis and modeling, of star formation regions using an ensemble of archival space-based data from the Infrared Space Observatory s Long Wavelength Spectrometer and Short Wavelength Spectrometer, but including as well some other spectroscopic databases. Four kinds of regions are considered in the studies: (1) disks around more evolved objects; (2) young, low or high mass pre-main sequence stars in star formation regions; (3) star formation in external, bright IR galaxies; and (4) the galactic center. One prime focus of the program is the OH lines in the far infrared. The program had the following goals: 1) Refine the data analysis of IS0 observations to obtain deeper and better SNR results on selected sources. The IS0 data itself underwent "pipeline 10" reductions in early 2001, and additional "hands-on data reduction packages" were supplied by the IS0 teams in 2001. The Fabry-Perot database is particularly sensitive to noise and slight calibration errors. 2) Model the atomic and molecular line shapes, in particular the OH lines, using revised monte- carlo techniques developed by the SWAS team at the Center for Astrophysics; 3) Attend scientific meetings and workshops; 4) Do E&PO activities related to infrared astrophysics and/or spectroscopy.

Oliversen, Ronald↗

Development of the Sasquatch Drop Test Footprint Tool

The Crew Exploration Vehicle Parachute Assembly System (CPAS) is the parachute system for NASA s Orion spacecraft. CPAS is currently in the design and testing phase of development. The test program consists of numerous drop tests, wherein a test article rigged with parachutes is extracted or released from an aircraft. During such tests, range safety is paramount, as is the recoverability of the parachutes and test article. It is crucial to establish an aircraft release point that will ensure that the article and all items released from it will land in safe locations. Early in the CPAS project, a legacy tool (previously used on the X-38 project) was used to determine a safe release point and to predict the landing locations (the footprint) of the payload and all released objects. Due to increasing test complexity and the need for a more flexible tool, a new footprint predictor tool, called Sasquatch, was created in MATLAB. This tool takes in a simulated trajectory for the test article, information about all released objects, and atmospheric wind data (simulated or actual) to calculate the trajectories of the released objects. Dispersions are applied to the landing locations of those objects, taking into account the variability of winds, aircraft release point, and object descent rate. A safe aircraft release point is determined based on the landing locations of the payload and released objects. The release point, landing locations, and dispersions are plotted on a simple map of the drop zone for easy reference. To date, Sasquatch has been used for thirteen drop tests. Comparing the predictions with actual test results has allowed for significant improvements in the tool s predictive capabilities, especially the incorporation of a well-correlated horizontal throw model. Intended future improvements to the tool include tighter dispersions on the landing locations, Monte Carlo capability, direct input from trajectory simulations, and a graphical user interface.

Bledsoe, Kristin J.↗

Ares I-X Separation and Reentry Trajectory Analyses

The Ares I-X Flight Test Vehicle was launched on October 28, 2009 and was the first and only test flight of NASA s two-stage Ares I launch vehicle design. The launch was successful and the flight test met all of its primary and secondary objectives. This paper discusses the stage separation and reentry trajectory analysis that was performed in support of the Ares I-X test flight. Pre-flight analyses were conducted to assess the risk of stage recontact during separation, to evaluate the first stage flight dynamics during reentry, and to define the range safety impact ellipses of both stages. The results of these pre-flight analyses were compared with available flight data. On-board video taken during flight showed that the flight test vehicle successfully separated without any recontact. Reconstructed trajectory data also showed that first stage flight dynamics were well characterized by pre-flight Monte Carlo results. In addition, comparisons with flight data indicated that the complex interference aerodynamic models employed in the reentry simulation were effective in capturing the flight dynamics during separation. Finally, the splash-down locations of both stages were well within predicted impact ellipses.

Tartabini, Paul V.↗

An Approach for the Assessment of System Upset Resilience

This report describes an approach for the assessment of upset resilience that is applicable to systems in general, including safety-critical, real-time systems. For this work, resilience is defined as the ability to preserve and restore service availability and integrity under stated conditions of configuration, functional inputs and environmental conditions. To enable a quantitative approach, we define novel system service degradation metrics and propose a new mathematical definition of resilience. These behavioral-level metrics are based on the fundamental service classification criteria of correctness, detectability, symmetry and persistence. This approach consists of a Monte-Carlo-based stimulus injection experiment, on a physical implementation or an error-propagation model of a system, to generate a system response set that can be characterized in terms of dimensional error metrics and integrated to form an overall measure of resilience. We expect this approach to be helpful in gaining insight into the error containment and repair capabilities of systems for a wide range of conditions.

Torres-Pomales, Wilfredo↗

Wear-Out Sensitivity Analysis Project Abstract

During the course of the Summer 2015 internship session, I worked in the Reliability and Maintainability group of the ISS Safety and Mission Assurance department. My project was a statistical analysis of how sensitive ORU's (Orbital Replacement Units) are to a reliability parameter called the wear-out characteristic. The intended goal of this was to determine a worst case scenario of how many spares would be needed if multiple systems started exhibiting wear-out characteristics simultaneously. The goal was also to determine which parts would be most likely to do so. In order to do this, my duties were to take historical data of operational times and failure times of these ORU's and use them to build predictive models of failure using probability distribution functions, mainly the Weibull distribution. Then, I ran Monte Carlo Simulations to see how an entire population of these components would perform. From here, my final duty was to vary the wear-out characteristic from the intrinsic value, to extremely high wear-out values and determine how much the probability of sufficiency of the population would shift. This was done for around 30 different ORU populations on board the ISS.

Harris, Adam↗

Comparison of the Integrated Medical Model Predictions to Real World ISS and STS Observations

The Human Research Program funded the development of the integrated medical model (IMM) to quantify the medical component of overall mission risk. The IMM uses Monte Carlo methodology to integrate space flight and ground medical data to assess the probability of mission medical outcomes and resource utilization. To determine the credibility of IMM output the IMM project team completed two validation studies that compare IMM output to observed medical events from a selection of Shuttle Transportation System (STS) and International Space Station (ISS) missions.

validation↗