Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Methods”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Direct simulation of typical AOTV entry flows

The direct simulation Monte Carlo method is applied to two flows that are relevant to the projected aeroassisted orbital transfer vehicle. The first is the axially symmetric flow about a 70 degree blunted cone shield with a cylindrical afterbody, and the second is the one-dimensional flow along the stagnation streamline of the first flow. The shock standoff distance for the stagnation streamline study is obtained from the full axially symmetric calculation. Advantage is then taken of the speed of the one-dimensional calculation to study the sensitivity of the stagnation point heat transfer to a number of physical variables. These include the chemical reaction rate coefficients, the surface catalycity, the electric field associated with ambipolar diffusion, and the nonequilibrium thermal radiation from the shock front.

Bird, G. A.↗
NASA NTRS1986

Bipropellant rocket exhaust plume analysis on the Galileo spacecraft

This paper describes efforts to quantify the contaminant flow field produced by 10 N thrust bipropellant rocket engines used on the Galileo spacecraft. The prediction of the composition of the rocket exhaust by conventional techniques is found to be inadequate to explain experimental observations of contaminant deposition on moderately cold (200 K) surfaces. It is hypothesized that low volatility contaminants are formed by chemical reactions which occur on the surfaces. The flow field calculations performed using the direct simulation Monte Carlo method give the expected result that the use of line-of-sight plume shields may have very little effect on the flux of vapor phase contaminant species to a surface, especially if the plume shields are located so close to the engine that the interaction of the plume with the shield is in the transition flow regime. It is shown that significant variations in the exhaust plume composition caused by nonequilibrium effects in the flow field lead to very low concentrations of species which have high molecular weights in the more rarefied regions of the flow field. Recommendations for the design of spacecraft plume shields and further work are made.

Guernsey, C. S.↗
NASA NTRS1986

Use of speckle for determining the response characteristics of Doppler imaging radars

An optical model is developed for imaging optical radars such as the SAR on Seasat and the Shuttle Imaging Radar (SIR-B) by analyzing the Doppler shift of individual speckles in the image. The signal received at the spacecraft is treated in terms of a Fresnel-Kirchhoff integration over all backscattered radiation within a Huygen aperture at the earth. Account is taken of the movement of the spacecraft along the orbital path between emission and reception. The individual points are described by integration of the point source amplitude with a Green's function scattering kernel. Doppler data at each point furnishes the coordinates for visual representations. A Rayleigh-Poisson model of the surface scattering characteristics is used with Monte Carlo methods to generate simulations of Doppler radar speckle that compare well with Seasat SAR data SIR-B data.

Tilley, D. G.↗
NASA NTRS1986

Nonlinear response - A time domain approach

The present paper reviews the basic concepts of nonlinear response of panels to surface flow and acoustic pressures, simulation of random processes, time domain solutions and the Monte Carlo Method. Applications of this procedure to the orbit-on-demand space vehicles, acoustic fatigue and composite materials are discussed. Numerical examples are included for a variety of nonlinear problems to illustrate the applicability of this method.

Vaicaitis, R.↗
NASA NTRS1986

Using a laser aureole to invert lidar return

An aureole generated by a laser beam was studied. The strength of the signal redirected towards a sensor high above the surface by a combination of one scattering event in the marine boundary layer (mbl) and one single reflection event from the ocean surface was estimated. A model of mbl aerosol size distributions was used to estimate Mie scattering for a wide range of meteorolocial conditions. The sea surface reflection was determined from a Gaussian model of the wave slopes. These laser aureoles which were estimated over the wide range of conditions and were normalized by the reflected laser light were found to be highly correlated with the optical depth of the boundary layer. By estimating optical depth from the aureole, the Bernoulli-Riccati inversion of lidar return could be constrained and the inversion accuracy improved. A Monte Carlo program was developed to study the laser aureole generated by up to 8 orders of reflection and scattering. The aureole was generated by a narrow, 10 nsec laser pulse at 1.06 microns and measured by a receiver 10 km above the ocean surface. The original theoretical computation compared well with the Monte Carlo method. When multiple scattereffects were included, the normalized aureole was still highly correlated with the mbl optical depth over the range of conditions.

Hooper, W. P.↗
NASA NTRS1987

Acoustic fatigue - A Monte Carlo approach

The present paper reviews the basic concepts of nonlinear structural response to surface flow and acoustic pressure inputs. A time domain solution and the Monte Carlo method are utilized to estimate the required statistics of the nonlinear response. Simulation procedures of multi-dimensional and multi-variate random processes are discussed. Application of these procedures to acoustic fatigue of flight structures are considered. Numerical examples are included to illustrate the applicability of this method for the solution of nonlinear problems.

Vaicaitis, R.↗
NASA NTRS1987

Nonequilibrium radiation during re-entry at 10 km/s

The direct simulation Monte Carlo method, including a real air model with thermal radiation, is applied to the flows associated with the two sets of measurements that are directly relevant to the projected aeroassisted orbital transfer vehicle. The first is a shock tube measurement of the radiation from a 10 km/s shock wave in air that was made at AVCO in 1962. The second is the flight data that was obtained from the Project Fire re-entry test vehicles in 1964. The calculations for both cases were made with a program that models the one-dimensional flow along a stagnation streamline. The shock standoff distance for the Fire vehicle was obtained from the theoretical studies that were associated with its launch. The simulation employed a partly phenomenological model for the nonequilibrium radiation. It was found that the results from the calculation were consistent with the measured radiation in each case, and also with the convective heat transfer data for the Fire vehicle. The uncertainties associated with the spectral absorptance and recombination probability at the surface appear to be as serious as those associated with the reaction rates.

Bird, G. A.↗
NASA NTRS1987

Estimating the arrival times of photon-limited laser pulses in the presence of shot and speckle noise

Maximum-likelihood (ML) receivers are frequently used to optimize the timing performance of laser-ranging and laser-altimetry systems in the presence of shot and speckle noise. Monte Carlo method was used to examine ML-receiver performance with return signals in the 10-5000-photoelectron (pe) range. The simulations were performed for shot noise only and for shot and speckle noise. The results agree with previous theory for signal strengths greater than about 100 pe's but show that the theory can significantly underestimate timing errors for weaker received signals. Sharp high-bandwidth features in the detected signals are shown to improve timing performance only if their signal levels are greater than 4-5 pe's.

Abshire, James B.↗
NASA NTRS1987

Effect of an electron scattering cloud on X-ray oscillations produced by beaming

The effect of a scattering cloud on the amplitude of oscillations produced by a rotating beam of X-rays is investigated using analytical and Monte Carlo methods. The scattering cloud was modeled as a uniform density sphere, and the source was represented as an anistropic distribution of radiation emerging from a point at the center of the scattering cloud. The intensity distribution produced by the source beam is examined as a function of optical depth. The relation between electron scattering optical depth and the forward-backward ratio is studied. It is observed that the scattering in a central corona of various optical depths reduces the amplitude of the oscillation. The data suggest that the quasi-periodic oscillations observed in the X-ray intensities of some luminous low-mass X-ray binaries are caused by oscillations in the luminosity of the X-ray star.

Brainerd, J.↗
NASA NTRS1988

Creep life prediction based on stochastic model of microstructurally short crack growth

A nondimensional model of microstructurally short crack growth in creep is developed based on a detailed observation of the creep fracture process of 304 stainless steel. In order to deal with the scatter of small crack growth rate data caused by microstructural inhomogeneity, a random variable technique is used in the model. A cumulative probability of the crack length at an arbitary time, G(bar a, bar t), and that of the time when a crack reaches an arbitary length, F(bar t, bar a), are obtained numerically by means of a Monte Carlo method. G(bar a, bar t), and F(bar t, bar a) are the probabilities for a single crack. However, multiple cracks generally initiate on the surface of a smooth specimen from the early stage of creep life to the final stage. TAking into account the multiple crack initiations, the actual crack length distribution observed on the surface of a specimen is predicted by the combination of probabilities for a single crack. The prediction shows a fairly good agreement with the experimental result for creep of 304 stainless steel at 923 K. The probability of creep life is obtained from an assumption that creep fracture takes place when the longest crack reaches a critical length. The observed and predicted scatter of the life is fairly small for the specimens tested.

Kitamura, Takayuki↗
NASA NTRS1986

The diffusion of individual molecules within a gas

The Direct Simulation Monte Carlo method is used to study the positional history of the individual molecules in a gas that is homogeneous at the macroscopic level and is in Maxwellian equilibrium at the microscopic level. The behavior at small times is characterized by 'persistence of velocity' effects, and a 'random walk' type of dispersal occurs over a longer timescale. It is shown that the rate of dispersal can be directly related to the self-diffusion coefficient. In addition, the diffusion coefficients are obtained directly from one-dimensional calculations, and the local Knudsen number at which the Chapman-Enskog theory breaks down is determined. Results are presented for both simple gases and gas mixtures.

Bird, G. A.↗
NASA NTRS1988

Nonequilibrium thermal radiation for an aeroassist flight experiment vehicle

The direct-simulation Monte Carlo method incorporating a dissociating and ionizing gas model for air with thermal radiation is used to characterize the hypersonic flow about an axisymmetric representation of an aeroassist flight experiment (AFE) vehicle, whose freestream conditions correspond to selected points along the entry, aerobraking, and exit phases of the trajectory. Calculations for two trajectory conditions indicate that the radiative heating of the AFE forebody is lower than the convective heating, but becomes significant as the maximum convective heating rate condition is approached.

Moss, James N.↗
NASA NTRS1988

Direct simulation of three-dimensional hypersonic flow about intersecting blunt wedges

The general three-dimensional direct simulation Monte Carlo method is used to study the hypersonic flow around two blunt wedges that intersect at a 90 deg angle. Results are obtained for the transitional flow regimes found at 85 and 100 km altitude with a reentry velocity of 7.5 km/s. The disturbance field in front of the double-wedge body is found to be larger than that produced by a single wedge. Surface pressures and flow densities are higher near the wedge intersection, whereas surface heating and shear streses are greater at locations removed from the corner. Results also show that three-dimensional flow structure occurs only near the wedge corner, and that the flow monotonically approaches the limiting two-dimensional wedge flow case in the spanwise direction.

Celenligil, M. Cevdet↗
NASA NTRS1987

Imaginary time path integral Monte Carlo route to rate coefficients for nonadiabatic barrier crossing

Nonadiabatic transitions are central to many areas of chemical and condensed matter physics, ranging from biological electron transfer to the optical properties of one-dimensional conductors. Here, a path integral Monte Carlo method is used to simulate such transitions, based on the observation that nonadiabatic rate coefficients are often dominated by saddle point trajectories that correspond to an imaginary time. Simple analytic theories can be used to continue these imaginary time correlation functions to determine rate coefficients. The advantages and drawbacks of this approach are discussed.

Wolynes, Peter G.↗
NASA NTRS1988

Interference effects on the hypersonic, rarefied flow about a flat plate

The Direct Simulation Monte Carlo method is used to study the hypersonic, rarified flow interference effects on a flat plate caused by nearby surfaces. Calculations focus on shock-boundary-layer and shock-lip interactions in hypersonic inlets. Results are presented for geometries consisting of a flat plate with different leading-edge shapes over a flat lower wall and a blunt-edge flat plate over a 5-degree wedge. The problems simulated correspond to a typical entry flight condition of 7.5 km/s at altitudes of 75 to 90 km. The results show increases in predicted local heating rates for shock-boundary-layer and shock-lip interactions that are quantitatively similar to those observed experimentally at much higher densities.

Wilmoth, Richard G.↗
NASA NTRS1988

Cyclotron line resonant transfer through neutron star atmospheres

Monte Carlo methods are used to study in detail the resonant radiative transfer of cyclotron line photons with recoil through a purely scattering neutron star atmosphere for both the polarized and unpolarized cases. For each case, the number of scatters, the path length traveled, the escape frequency shift, the escape direction cosine, the emergent frequency spectra, and the angular distribution of escaping photons are investigated. In the polarized case, transfer is calculated using both the cold plasma e- and o-modes and the magnetic vacuum perpendicular and parallel modes.

Wang, John C. L.↗
NASA NTRS1989

Influence of afterbodies on AOTV flows

The role of vibration-dissociation coupling in flows representative of AOTV's with afterbodies is discussed. A physical model designed to reflect the behavior of anharmonic diatomic molecules at high temperatures was developed and incorporated into the direct simulation Monte Carlo method of Bird (1976). The model leads to reduced dissociation and increased convective heating for flows with finite catalytic surfaces.

Olynick, David P.↗
NASA NTRS1989

Monte Carlo simulation of reentry plasmas

Attention is given to the treatment of ionization and plasma effects in the direct simulation Monte Carlo method. The requirements for accurate modeling of reentry plasmas are discussed along with the difficulties these requirements present. The current method for modeling such plasmas is reviewed and an alternative method is presented. Both methods are applied to the flow of a 10 km/s shock wave in air at 0.1 torr; a flowfield directly relevant to the projected aeroassisted orbital transfer vehicle. The results are compared and the differences between the methods are discussed.

Carlson, Ann B.↗