Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Method”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

SHDOM

SHDOM is a general purpose, publicly available, three-dimensional atmospheric radiative transfer model. SHDOM is an explicit method, which means it solves for the whole radiation field, as distinct from Monte Carlo methods which solve for particular radiative outputs. SHDOM is particularly well suited for remote sensing applications, where it can compute outgoing radiances at many angles from a cloud field at virtually no extra cost. SHDOM is not appropriate for calculating domain average quantities for which Monte Carlo methods excel. The I3RC intercomparison offers an opportunity to explore the pros and cons of SHDOM and Monte Carlo models on some real world inhomogeneous cloud fields. Specifically, we wish to determine the computer resources required to achieve a particular accuracy for a certain number of outputs using SHDOM and Monte Carlo models. This will help guide modelers on the appropriate choice of SHDOM or Monte Carlo for their applications. To emphasize the importance of this accuracy versus CPU time tradeoff, we are submitting two SHDOM entries (low and high resolution) in the I3RC.

Evans, K. Franklin↗
NASA NTRS1994

Radiative interactions in multi-dimensional chemically reacting flows using Monte Carlo simulations

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The amount and transfer of the emitted radiative energy in a finite volume element within a medium are considered in an exact manner. The spectral correlation between transmittances of two different segments of the same path in a medium makes the statistical relationship different from the conventional relationship, which only provides the non-correlated results for nongray methods is discussed. Validation of the Monte Carlo formulations is conducted by comparing results of this method of other solutions. In order to further establish the validity of the MCM, a relatively simple problem of radiative interactions in laminar parallel plate flows is considered. One-dimensional correlated Monte Carlo formulations are applied to investigate radiative heat transfer. The nongray Monte Carlo solutions are also obtained for the same problem and they also essentially match the available analytical solutions. the exact correlated and non-correlated Monte Carlo formulations are very complicated for multi-dimensional systems. However, by introducing the assumption of an infinitesimal volume element, the approximate correlated and non-correlated formulations are obtained which are much simpler than the exact formulations. Consideration of different problems and comparison of different solutions reveal that the approximate and exact correlated solutions agree very well, and so do the approximate and exact non-correlated solutions. However, the two non-correlated solutions have no physical meaning because they significantly differ from the correlated solutions. An accurate prediction of radiative heat transfer in any nongray and multi-dimensional system is possible by using the approximate correlated formulations. Radiative interactions are investigated in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The governing equations are based on the fully elliptic Navier-Stokes equations. Chemical reaction mechanisms were described by a finite rate chemistry model. The correlated Monte Carlo method developed earlier was employed to simulate multi-dimensional radiative heat transfer. Results obtained demonstrate that radiative effects on the flowfield are minimal but radiative effects on the wall heat transfer are significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and nozzle size on the radiative and conductive wall fluxes.

Liu, Jiwen↗
NASA NTRS1981

Direct simulation Monte Carlo technique for modeling of the environment in the vicinity of the Space Shuttle Orbiter

A program for analyzing the flowfield parameters in the neighborhood of the Space Shuttle Orbiter has been developed. The program uses the direct simulation Monte Carlo method, which is a completely probabilistic Monte Carlo technique capable of analyzing 3-dimensional steady or unsteady flow with prescribed internal and external boundary conditions. The freestream flux densities incident on the external flowfield boundaries are calculated from the drifting Maxwellian gas properties of the freestream. The flux entering the flowfield from the internal boundary is calculated from the outgassing flux density distribution over the Shuttle external surface and the prescribed discrete source fluxes. This technique produces a numerical flowfield solution which is the probabilistic equivalent of a complete solution of the time-dependent, 3-dimensional Boltzmann equation. Flowfield results are presented for the following configurations: (1) Shuttle angles of attack of 0 and 90 deg, (2) freestream density 10 to the 9th - 10 to the 11th per cu cm, (3) Shuttle outgassing flux density, and (4) operation of the aft downfiring vernier Reaction Control System (RCS) engine. Results are presented for column density distribution of outgassed and engine species and for the flux density of outgassed and engine species incident on the Shuttle bay and the windshield.

Melfi, L. T., Jr.↗
NASA NTRS1986

Composite load spectra for select space propulsion structural components

A multiyear program is performed with the objective to develop generic load models with multiple levels of progressive sophistication to simulate the composite (combined) load spectra that are induced in space propulsion system components, representative of Space Shuttle Main Engines (SSME), such as transfer ducts, turbine blades, and liquid oxygen (LOX) posts. Progress of the first year's effort includes completion of a sufficient portion of each task -- probabilistic models, code development, validation, and an initial operational code. This code has from its inception an expert system philosophy that could be added to throughout the program and in the future. The initial operational code is only applicable to turbine blade type loadings. The probabilistic model included in the operational code has fitting routines for loads that utilize a modified Discrete Probabilistic Distribution termed RASCAL, a barrier crossing method and a Monte Carlo method. An initial load model was developed by Battelle that is currently used for the slowly varying duty cycle type loading. The intent is to use the model and related codes essentially in the current form for all loads that are based on measured or calculated data that have followed a slowly varying profile.

Newell, J. F.↗
NASA NTRS1997

Efficient algorithms for robust feature matching

One of the basic building blocks in any point-based registration scheme involves matching feature points that are extracted from the sensed image to their counterparts in the reference image. This leads to the fundamental problem of point matching: given two sets of points, find the affine transformation that transforms one point set so that its distance from the other point set is minimized. Because of measurement errors and the presence of outlying data points, it is important that the distance measure between two point sets be robust to these effects. We measure distances using the generalized Hausdorff distance. Point matching can be a computationally intensive task, and there have been a number of algorithms and approaches proposed for solving this problem both theoretical and applied. We present two approaches to the point matching problem, in an attempt to reduce the computational complexity of the problem, while still providing guarantees on the quality of the final match. Our first method is an approximation algorithm, which is loosely based on a branch-and-bound approach due to Huttenlocher and Rucklidge. We show that by varying the approximation error bounds, it is possible to achieve a tradeoff between the quality of the match and the running time of the algorithm. Our second method involves a Monte Carlo method for accelerating the search process used in the first algorithm. With high probability this method succeeds in finding an approximately optimal match. We establish the efficiency of our approaches empirically.

Mount, David M.↗
NASA NTRS2001

Transport Phenomena of Off-Axis Sputtering Deposition

Various high quality epitaxial films, especially oxides, have been synthesized using off-axis sputtering deposition. In this presentation, we report the experiment results of ZnO films grown by the off-axis sputtering deposition. Films were synthesized in temperatures ranged from room temperature to 600 C, and pressures from 5 mTorr to 150 mTorr. Film growth rate was measured by surface profilometer, ellipsometer, and wavelength dispersive spectrometry. Due to the collisions between the sputtered species and the residue gases, the kinetic energy of species was reduced and the transport of depositing species changed from a ballistic movement for low pressure to a diffuse drift for high pressure in which the transport species were almost thermalized. The measurements show an increase of growth rates along the gravity vector when the Knodson (Knudsen??) number of transport species is less than 0.05, which suggests that gravity affected the transport characterization in off-axis sputtering deposition. Because the product of pressure (p) and travel distance (d) of sputtered species, p exceeds several mTorr-cm during film deposition, the classical simulations for sputtering process in high vacuum system may not be applied. Based on these experimental measurements, a transport process of the off-axis sputtering deposition is proposed. Several methods including the Monte Carlo method and gravity-driven flow dynamics simulation will be discussed.

Zhu, S.↗
NASA NTRS2021

Measurement Uncertainty Analysis of 6x9 Icing Research Tunnel (IRT)

This paper presents the results of the measurement uncertainty analysis that was performed on the aero-thermal characteristics of the 6- by 9-Foot Icing Research Tunnel (IRT) at the NASA Glenn Research Center. The engineering judgments and statistical methods used to determine the elemental uncertainties are described. Instrumentation uncertainty was quanti ed using MANTUS (Measurement Analysis Tool for Uncertainty in Systems) and instrument speci cation sheets. The Monte Carlo method was used to propagate systematic components of uncertainty in order to quantify the uncertainty of the Variables of Interest (VOI). A detailed description of the Monte Carlo method and the MANTUS tool can be found in the Uncertainty Analysis of the NASA Glenn 8x6 Supersonic Wind Tunnel report [1]. Detailed uncertainty results for test section airspeed and temperature as well as other variables of interest are described within this document.

6x9 Icing Research Tunnel↗
NASA NTRS1994

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗
NASA NTRS1973

LP-search and its use in analysis of the accuracy of control systems with acoustical models

The LP-search is proposed as an analog of the Monte Carlo method for finding values in nonlinear statistical systems. It is concluded that: To attain the required accuracy in solution to the problem of control for a statistical system in the LP-search, a considerably smaller number of tests is required than in the Monte Carlo method. The LP-search allows the possibility of multiple repetitions of tests under identical conditions and observability of the output variables of the system.

Sergeyev, V. I.↗
NASA NTRS1994

A Monte Carlo study of photochemical model uncertainties affecting the calculation of subsonic perturbations in the upper troposphere

Monte Carlo methods are frequently applied to the evaluation of uncertainties in models with multiple inputs that themselves have associated imprecisions. In the case of photochemical models used to evaluate changes in O3 or OH, inputs analyzed include rate coefficients measured in the laboratory, and chemical and physical constituents measured in the atmosphere. The Monte Carlo method was used with the 1-dimensional GSFC tropospheric photochemical model to examine uncertainty propagation to calculation of NO(x) and the major oxidants (O3, OH) in the upper troposphere. In all cases chemical kinetics inputs are varied and the NO(x) perturbation of a subsonic fleet is simulated. A series of model runs is used to explore sensitivities of model-computed parameters to other parameters: heterogeneous processes, the uncertainty in upper tropospheric H2O vapor measurements, aircraft emissions at different latitudes.

Thompson, Anne M.↗
NASA NTRS1978

Three-dimensional radiative exchange factors for nongray, nondiffuse surfaces

A computer program using a variation of the Monte Carlo method was developed to predict radiative exchange factors in three-dimensional configurations with curved surfaces. The model allows for arbitrary emission, absorption, and reflection characteristics and for radiation exchange between surfaces as well as external irradiation through openings. It was shown that exchange factors can be determined efficiently by using the normal Monte Carlo method for closed configurations and/or high surface emissivities, while energy partitioning is vastly superior in the case of open configurations and/or low surface emissivities. Comparison with some experiments performed on the Space Shuttle heat rejection system demonstrates the model's accuracy as well as its superiority over a simplified diffuse-surface analysis, which requires similar amounts of computer time.

Modest, M. F.↗
NASA NTRS1992

Solution of the Burnett equations for hypersonic flows near the continuum limit

The INCA code, a three-dimensional Navier-Stokes code for analysis of hypersonic flowfields, was modified to analyze the lower reaches of the continuum transition regime, where the Navier-Stokes equations become inaccurate and Monte Carlo methods become too computationally expensive. The two-dimensional Burnett equations and the three-dimensional rotational energy transport equation were added to the code and one- and two-dimensional calculations were performed. For the structure of normal shock waves, the Burnett equations give consistently better results than Navier-Stokes equations and compare reasonably well with Monte Carlo methods. For two-dimensional flow of Nitrogen past a circular cylinder the Burnett equations predict the total drag reasonably well. Care must be taken, however, not to exceed the range of validity of the Burnett equations.

Imlay, Scott T.↗
NASA NTRS1982

Transport of energetic electrons in a fully ionized hydrogen plasma

A Monte Carlo method for calculating energetic electron transport in a plasma is presented. The energy loss and angular deflection due to Coulomb collisions as a function of travel distance and energy are derived for energetic electrons in a fully ionized plasma with a uniform magnetic field. Formulas which include the effect of nonuniform B fields on the angular deflection are derived. The Monte Carlo method is applied to the thick-target beam model in which the energetic electrons are injected vertically downward.

Bai, T.↗
NASA NTRS1986

Monte Carlo study of a model of diffusion-controlled reactions

The Monte Carlo method is used to perform averages over sink configurations in the present study of diffusion-controlled reactions occurring between solute particles and immobile spherical sinks. In order to determine the average steady state solute concentration profile in a locally perturbed solution for sink volume fractions phi of less than 0.3, the diffusion equation in the monopolar plus dipolar approximation of diffusive couplings between the sinks is solved numerically. The Monte Carlo method is shown to be the most accurate and efficient in the phi = 0.001-0.1 region, where a system of only 25 sinks suffices and the monopolar approximation alone is sufficiently accurate.

Beenakker, C. W. J.↗
NASA NTRS1972

Monte Carlo investigation of transient acoustic fields in partially or completely bounded medium

A simple repetitive calculation was used to investigate what happens to the field in terms of the signal paths of disturbances originating from the energy source. The computation allowed the field to be reconstructed as a function of space and time on a statistical basis. The suggested Monte Carlo method is in response to the need for a numerical method to supplement analytical methods of solution which are only valid when the boundaries have simple shapes, rather than for a medium that is bounded. For the analysis, a suitable model was created from which was developed an algorithm for the estimation of acoustic pressure variations in the region under investigation. The validity of the technique was demonstrated by analysis of simple physical models with the aid of a digital computer. The Monte Carlo method is applicable to a medium which is homogeneous and is enclosed by either rectangular or curved boundaries.

Thanedar, B. D.↗