Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Monte Carlo Simulations of Crystal Defects in Open Ensembles

Zero- and two-dimensional crystal defects form in open statistical ensembles, such as the grand canonical, that are usually inaccessible with conventional simulation techniques. This longstanding challenge is overcome with a new Hamiltonian Monte Carlo method that samples energy-biased gradual transformations. In conclusion, the method enables free energy calculations for nonideal point defects and the direct prediction of finite-temperature interface structures.

Grain boundaries↗
NASA NTRS2010

Accelerated Monte Carlo Simulation for Safety Analysis of the Advanced Airspace Concept

Safe separation of aircraft is a primary objective of any air traffic control system. An accelerated Monte Carlo approach was developed to assess the level of safety provided by a proposed next-generation air traffic control system. It combines features of fault tree and standard Monte Carlo methods. It runs more than one order of magnitude faster than the standard Monte Carlo method while providing risk estimates that only differ by about 10%. It also preserves component-level model fidelity that is difficult to maintain using the standard fault tree method. This balance of speed and fidelity allows sensitivity analysis to be completed in days instead of weeks or months with the standard Monte Carlo method. Results indicate that risk estimates are sensitive to transponder, pilot visual avoidance, and conflict detection failure probabilities.

Thipphavong, David↗
NASA NTRS2011

Kinetic Monte Carlo Simulation of Oxygen and Cation Diffusion in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia (YSZ) is of interest to the aerospace community, notably for its application as a thermal barrier coating for turbine engine components. In such an application, diffusion of both oxygen ions and cations is of concern. Oxygen diffusion can lead to deterioration of a coated part, and often necessitates an environmental barrier coating. Cation diffusion in YSZ is much slower than oxygen diffusion. However, such diffusion is a mechanism by which creep takes place, potentially affecting the mechanical integrity and phase stability of the coating. In other applications, the high oxygen diffusivity of YSZ is useful, and makes the material of interest for use as a solid-state electrolyte in fuel cells. The kinetic Monte Carlo (kMC) method offers a number of advantages compared with the more widely known molecular dynamics simulation method. In particular, kMC is much more efficient for the study of processes, such as diffusion, that involve infrequent events. We describe the results of kinetic Monte Carlo computer simulations of oxygen and cation diffusion in YSZ. Using diffusive energy barriers from ab initio calculations and from the literature, we present results on the temperature dependence of oxygen and cation diffusivity, and on the dependence of the diffusivities on yttria concentration and oxygen sublattice vacancy concentration. We also present results of the effect on diffusivity of oxygen vacancies in the vicinity of the barrier cations that determine the oxygen diffusion energy barriers.

Good, Brian↗
NASA NTRS2003

Monte Carlo Simulations of the Response of the MARIE Instrument

The MARIE instrument aboard Mars Odyssey functions as a telescope for the detection of charged, energetic, nuclei. The directionality that leads to the telescope description is achieved by requiring coincident signals in two designated detectors in MARIE s silicon detector stack for the instrument to trigger. Because of this, MARIE is actually a bi directional telescope. Triggering particles can enter the detector stack by passing through the lightly shielded front of the instrument, but can also enter the back of the instrument by passing through the bulk of Odyssey. Because of this, understanding how to relate the signals recorded by MARIE to astrophysically important quantities such as particle fluxes or spectra exterior to the spacecraft clearly requires detailed modeling of the physical interactions that occur as the particles pass through the spacecraft and the instrument itself. In order to facilitate in the calibration of the MARIE data, we have begun a program to simulate the response of MARIE using the FLUKA [1] [2] Monte Carlo radiation transport code.

Andersen, V.↗
NASA NTRS1990

Hypersonic flow simulations using DSMC (direct simulation Monte Carlo)

A review of the DSMC method of Bird is presented. The DSMC method provides the capability of simulating real gas flows in the rarefied flow regime. Recent developments and applications of the method for hypersonic flows are reported for both ground-based tests and during entry. Results obtained using both axisymmetric and 3D codes are included.

Moss, James N.↗
NASA NTRS2018

High Fidelity Collision Probabilities Estimated Using Brute Force Monte Carlo Simulations

The NASA Conjunction Assessment Risk Analysis team has implemented new software to estimate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm employs a brute force Monte Carlo (BFMC) method which differs from most other methods because it uses orbital states and covariances propagated from their orbit determination epoch times using the full set of Astrodynamics Support Workstation higher order theory models, including the High Accuracy Satellite Drag Model. This paper de-scribes the BFMC algorithm, presents comparisons of BFMC P (sub c) estimates to those calculated using other methods, and discusses the implications for conjunction risk assessment.

Probabilities↗
NASA NTRS2002

Using Monte-Carlo Simulations to Study the Disk Structure in Cygnus X-1

As the first dynamically determined black hole X-ray binary system, Cygnus X-1 has been studied extensively. However, its broad-band spectra in hard state with BeppoSAX is still not well understood. Besides the soft excess described by the multi-color disk model (MCD), the power- law component and a broad excess feature above 10 keV (disk reflection component), there is also an additional soft component around 1 keV, whose origin is not known currently.We propose that the additional soft component is due to the thermal Comptonization process between the s oft disk photon and the warm plasma cloud just above the disk.i.e., a warm layer. We use Monte-Carlo technique t o simulate this Compton scattering process and build several table models based on our simulation results.

Yao, Y.↗
NASA NTRS2003

A Study of the Impact of Variations on Aerodynamic Flow in Gas Turbine Engines via Monte-Carlo Simulations

The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented with a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, under the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions that motivate this research are: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance. This paper presents an initial attempt to the flow behavior in terms (perturbation) using simulation outputs from a customer-deck model obtained from Pratt&Whitney, (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine. Flow behavior is analyzed for the high pressure compressor and combustion chamber employing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000-sample simulation. Furthermore, the characteristic waves and eave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures; which are clearly the case in the combustion chamber and high-pressure turbine. Finally a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Ngo, Khiem Viet↗
NASA NTRS2004

Impact of Variations on 1-D Flow in Gas Turbine Engines via Monte Carlo Simulations

The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented using a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, to determine the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions motivate this research: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance? This paper presents an initial attempt to model the flow behavior in terms of their eigenfrequencies subject to the assumption of the uncertainty or variation (perturbation). The flow behavior is explored using simulation outputs from a customer-deck model obtained from Pratt & Whitney. Variations of the main variables (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine are modeled. Flow behavior is analyzed for the high-pressure compressor and combustion chamber utilizing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high-pressure and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000 sample simulation. Furthermore, the characteristic waves and wave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures in the case of the combustion chamber and high-pressure turbine. Finally, a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Ngo, Khiem Viet↗
NASA NTRSDate not supplied

Applying Monte Carlo Simulation to Launch Vehicle Design and Requirements Verification

- This presentation applies statistics to launch vehicle design, but the methods may be used for other engineering applications - Vehicle Models: What do we know when? How to separate parameters for making Monte Carlo runs - Known when assembling the vehicle but not during design - Known prior to committing to flight - Unknown at lift off - How do we correctly model the vehicle during the design phases? - Number of Monte Carlo samples - Requirements success - Design parameter values

Monte Carlo↗
NASA NTRS1990

Electron heating in quasi-perpendicular shocks - A Monte Carlo simulation

To study the problem of electron heating in quasi-perpendicular shocks, under the combined effects of 'reversible' motion, in the shock electric potential and magnetic field, and wave-particle interactions a diffusion equation was derived, in the drift (adiabatic) approximation and it was solved by using a Monte Carlo method. The results show that most of the observations can be explained within this framework. The simulation has also definitively shown that the electron parallel temperature is determined by the dc electromagnetic field and not by any wave particle induced heating. Wave-particle interactions are effective in smoothing out the large gradients in phase space produced by the 'reversible' motion of the electrons, thus producing a 'cooling' of the electrons. Some constraints on the wave-particle interaction process may be obtained from a detailed comparison between the simulation and observations. In particular, it appears that the adiabatic approximation must be violated in order to explain the observed evolution of the perpendicular temperature.

Veltri, Pierluigi↗
NASA NTRS1998

Monte Carlo Simulations of Background Spectra in Integral Imager Detectors

Predictions of the expected gamma-ray backgrounds in the ISGRI (CdTe) and PiCsIT (Csl) detectors on INTEGRAL due to cosmic-ray interactions and the diffuse gamma-ray background have been made using a coupled set of Monte Carlo radiation transport codes (HETC, FLUKA, EGS4, and MORSE) and a detailed, 3-D mass model of the spacecraft and detector assemblies. The simulations include both the prompt background component from induced hadronic and electromagnetic cascades and the delayed component due to emissions from induced radioactivity. Background spectra have been obtained with and without the use of active (BGO) shielding and charged particle rejection to evaluate the effectiveness of anticoincidence counting on background rejection.

Armstrong, T. W.↗
NASA NTRS1992

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

Turbulent combustion can not be simulated adequately by conventional moment closure turbulent models. The probability density function (PDF) method offers an attractive alternative: in a PDF model, the chemical source terms are closed and do not require additional models. Because the number of computational operations grows only linearly in the Monte Carlo scheme, it is chosen over finite differencing schemes. A grid dependent Monte Carlo scheme following J.Y. Chen and W. Kollmann has been studied in the present work. It was found that in order to conserve the mass fractions absolutely, one needs to add further restrictions to the scheme, namely alpha(sub j) + gamma(sub j) = alpha(sub j - 1) + gamma(sub j + 1). A new algorithm was devised that satisfied this restriction in the case of pure diffusion or uniform flow problems. Using examples, it is shown that absolute conservation can be achieved. Although for non-uniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Hsu, Andrew T.↗
NASA NTRS2019

Assessing GEO and LEO Repeating Conjunctions Using High Fidelity Brute Force Monte Carlo Simulations

Probability of collision (P(sub c)) estimates for Earth-orbiting satellites typically assume a temporally-isolated conjunction event. However, under certain conditions two objects may experience multiple high-risk close approach events over the course of hours or days. In these repeating conjunction cases, the P(sub c) accumulates as each successive encounter occurs. The NASA Conjunction Assessment Risk Analysis team has updated its “brute force Monte Carlo” (BFMC) software to estimate such accumulating P(sub c) values for repeating conjunctions. This study describes the updated BFMC algorithm and discusses the implications for conjunction risk assessment.

Baars, Luis↗
NASA NTRS2019

Assessing Geo and Leo Repeating Conjunctions Using High Fidelity Brute Force Monte Carlo Simulations

Probability of collision (P(sub c)) estimates for Earth-orbiting satellites typically assume a temporally-isolated conjunction event. However, under certain conditions two objects may experience multiple high-risk close approach events over the course of hours or days. In these repeating conjunction cases, the P(sub c) accumulates as each successive encounter occurs. The NASA Conjunction Assessment Risk Analysis team has updated its “brute force Monte Carlo” (BFMC) software to estimate such accumulating P(sub c) values for repeating conjunctions. This study describes the updated BFMC algorithm and discusses the implications for conjunction risk assessment.

Baars, Luis↗