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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Modelling of SAR polarisation phase difference from trees

The data for polarization phase difference Delta Phi between the HH- and VV-polarized backscattered waves from tree-covered fields were obtained with an airborne synthetic aperture radar at 1.225 GHz. The mean values over tree-covered fields were derived from the images of the phase difference and were examined as a function of incident beam angle from 15 to 55 deg. A theoretical model for simulating these data, based on the electromagnetic wave scatterings from the tree trunk and its branches, both of which are assumed as very long dielectric cylinders was developed. The radius and direction of a tree branch are taken as random variables and are chosen by a Monte Carlo method to encounter the incident waves in producing the scattering events. The Monte Carlo simulated results are in good agreement with the observations within experimental uncertainty.

Mo, Tsan↗
NASA NTRS1994

Microwave radiative transfer through horizontally inhomogeneous precipitating clouds

Recent advances in cloud microphysical models have led to realistic three-dimensional distributions of cloud constituents. Radiative transfer schemes can make use of this detailed knowledge in order to study the effects of horizontal as well as vertical inhomogeneities within clouds. This study looks specifically at the differences between three-dimensional radiative transfer results and those obtained by plane parallel, independent pixel approximations in the microwave spectrum. A three-dimensional discrete ordinates method as well as a backward Monte Carlo method are used to calculate realistic radiances emerging from the cloud. Analyses between these models and independent pixel approximations reveal that plane parallel approximations introduce two distinct types of errors. The first error is physical in nature and is related to the fact that plane parallel approximations do not allow energy to leak out of dense areas into surrouding areas. In general, it was found that these errors are quite small for emission-dominated frequencies (37 GHz and lower) and that physical errors are highly pronounced only at scattering frequencies (85 GHz) where large deviations and biases up to 8 K averaged over the entire cloud were found. The second error is more geometric in nature and is related to the fact that plane parallel approximations cannot accommodate physical boundaries in the horizontal dimension for off-nadir viewing angles. The geometric errors were comparable in magnitude for all frequencies. Their magnitude, however, depends on a number of factors including the scheme used to deal with the edge, the nature of the surface, and the viewing angle.

Roberti, Laura↗
NASA NTRS1995

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗
NASA NTRSDate not supplied

Improving Multi-Model Trajectory Simulation Estimators using Model Selection and Tuning

Multi-model Monte Carlo methods have been demonstrated to be an efficient and accurate alternative to standard Monte Carlo (MC) in the model-based propagation of uncertainty in entry, descent, and landing (EDL) applications. These multi-model MC methods fuse predictions from low-fidelity models with the high-fidelity EDL model of interest to produce unbiased statistics with a fraction of the computational cost. The accuracy and efficiency of the multi-model MC methods are dependent upon the magnitude of correlations of the low-fidelity models with the high-fidelity model, but also upon the correlation amongst the low-fidelity models, and their relative computational cost. Because of this layer of complexity, the question of how to optimally select the set of low-fidelity models has remained open. In this work, methods for optimal model construction and tuning are investigated as a means to increase the speed and precision of trajectory simulation for EDL. Specifically, the focus is on the inclusion of low-fidelity model tuning within the sample allocation optimization that accompanies multi-model MC methods. Preliminary results indicate that low-fidelity model tuning can significantly improve efficiency and precision of trajectory simulations and provide an increased edge to multi-model MC methods when compared to standard MC. The challenges and potential benefits to exploring a fully iterative and comprehensive optimization strategy in future work are highlighted.

uncertainty quantification↗
NASA NTRS1982

Backscattered UV radiation - Effects of multiple scattering and the lower boundary of the atmosphere

A method is proposed for the calculation of a multiple-scattering correction to the single-scattering calculation of the radiance of the terrestrial atmosphere resulting from backscattered ultraviolet solar radiation in the spectral region used in the ozone profile inversion. This method uses jointly the usual analytical and Monte Carlo methods. Effects of the lower boundary of the atmosphere, cloud tops, and ground surface are investigated both qualitatively and quantitatively. The ratio of multiple to single scattering is determined, and its importance in ozone profile inversion of backscattered UV solar radiation from the terrestrial atmosphere is evaluated. The polarization of the atmospheric radiance is treated briefly.

Aruga, T.↗
NASA NTRS2003

Towards a Probabilistic Preliminary Design Criterion for Buckling Critical Composite Shells

A probability-based analysis method for predicting buckling loads of compression-loaded laminated-composite shells is presented, and its potential as a basis for a new shell-stability design criterion is demonstrated and discussed. In particular, a database containing information about specimen geometry, material properties, and measured initial geometric imperfections for a selected group of laminated-composite cylindrical shells is used to calculate new buckling-load "knockdown factors". These knockdown factors are shown to be substantially improved, and hence much less conservative than the corresponding deterministic knockdown factors that are presently used by industry. The probability integral associated with the analysis is evaluated by using two methods; that is, by using the exact Monte Carlo method and by using an approximate First-Order Second- Moment method. A comparison of the results from these two methods indicates that the First-Order Second-Moment method yields results that are conservative for the shells considered. Furthermore, the results show that the improved, reliability-based knockdown factor presented always yields a safe estimate of the buckling load for the shells examined.

Arbocz, Johann↗
NASA NTRS2012

Computing Radiative Transfer in a 3D Medium

A package of software computes the time-dependent propagation of a narrow laser beam in an arbitrary three- dimensional (3D) medium with absorption and scattering, using the transient-discrete-ordinates method and a direct integration method. Unlike prior software that utilizes a Monte Carlo method, this software enables simulation at very small signal-to-noise ratios. The ability to simulate propagation of a narrow laser beam in a 3D medium is an improvement over other discrete-ordinate software. Unlike other direct-integration software, this software is not limited to simulation of propagation of thermal radiation with broad angular spread in three dimensions or of a laser pulse with narrow angular spread in two dimensions. Uses for this software include (1) computing scattering of a pulsed laser beam on a material having given elastic scattering and absorption profiles, and (2) evaluating concepts for laser-based instruments for sensing oceanic turbulence and related measurements of oceanic mixed-layer depths. With suitable augmentation, this software could be used to compute radiative transfer in ultrasound imaging in biological tissues, radiative transfer in the upper Earth crust for oil exploration, and propagation of laser pulses in telecommunication applications.

Von Allmen, Paul↗
NASA NTRS1972

The response of a 300 micron silicon detector to monoenergetic neutrons determined by the use of the Monte Carlo technique

The response of a 300 micron thick silicon detector to an incident monoenergetic neutron beam is evaluated by the Monte Carlo method for the cases of both a shielded and a bare detector. The result of Monte Carlo calculation, using elastic, inelastic, and absorption reactions indicates that the response of the silicon detector to neutrons is basically due to the elastic scattering. In addition, the gamma rays generated in the shield of the detector will result in a response which is 3 or 4 orders of magnitude smaller than response to incident photons. The response of a bare silicon detector is calculated for neutron energies up to 6 MeV and bias energies from 50 to 250 KeV. It is found that the maximum response for a 300 micron thick silicon detector is less than .004 c/n within this selected neutron and bias energy range. When the pulse height defect is introduced in the calculation the results at low energy neutrons were reduced.

Tahezadeh, M.↗
NASA NTRS1991

Monte Carlo study of vibrational relaxation processes

A new model is proposed for the computation of vibrational nonequilibrium in the direct simulation Monte Carlo method (DSMC). This model permits level to level vibrational transitions for the first time in a Monte Carlo flowfield simulation. The model follows the Landau-Teller theory for a harmonic oscillator in which the rates of transition are related to an experimental correlation for the vibrational relaxation time. The usual method for simulating such processes in the DSMC technique applies a constant exchange probability to each collision and the vibrational energy is treated as a continuum. A comparison of these two methods is made for the flow of nitrogen over a wedge. Significant differences exist for the vibrational temperatures computed. These arise as a consequence of the incorrect application of a constant exchange probability in the old method. It is found that the numerical performances of the two vibrational relaxation models are equal.

Boyd, Iain D.↗
NASA NTRS2016

A DSMC Surface Chemistry Model for Carbon-Based Ablators

A detailed molecular surface chemistry model for the DSMC (Direct Simulation Monte Carlo) method is proposed and implemented into the SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) DSMC solver. Molchanova et al. constructed a molecular model for surface recombination in DSMC that includes different surface processes (adsorption, desoprtion, Eley-Rideal and Langmuir-Hinshelwood). All surface processes can be divided into two groups: surface mechanisms, which involve only the particle adsorbed by the surface (desorption and Langmuir-Hinshelwood), and impact mechanisms, which also involve gas-phase particles (adsorption, Eley-Rideal). Using a similar approach, the 14-reaction kinetic model of oxygen-carbon interaction suggested by Zhlukhtov and Abe, as well as more recent models by Alba et al., Poovathinghal et al., and a new model developed in the scope of this work, are implemented in SPARTA. The computational results for the different oxidation models are compared with experimental results from Murray et al. (oxidation of a vitreous carbon surface due to a hyperthermal beam of O and O2), with a particular focus on fluxes, angular and Time-Of-Flight distributions of scattered particles.

oxidation↗
NASA NTRS1980

Assessment of risk due to the use of carbon fiber composites in commercial and general aviation

The development of a national risk profile for the total annual aircraft losses due to carbon fiber composite (CFC) usage through 1993 is discussed. The profile was developed using separate simulation methods for commercial and general aviation aircraft. A Monte Carlo method which was used to assess the risk in commercial aircraft is described. The method projects the potential usage of CFC through 1993, investigates the incidence of commercial aircraft fires, models the potential release and dispersion of carbon fibers from a fire, and estimates potential economic losses due to CFC damaging electronic equipment. The simulation model for the general aviation aircraft is described. The model emphasizes variations in facility locations and release conditions, estimates distribution of CFC released in general aviation aircraft accidents, and tabulates the failure probabilities and aggregate economic losses in the accidents.

Fiksel, J.↗
NASA NTRS1982

Effects of translational nonequilibrium on vacuum plume expansions

The plume flowfield of a helium vent nozzle exhausting into a vacuum is defined by two techniques: the method of characteristics and the direct-simulation Monte Carlo method. The method of characteristics is shown to severely underpredict gas static temperatures due to the assumption of translational equilibrium inherent in the continuum formulation. Mass flow per unit solid angle is shown to be accurately predicted for this flowfield by the method of characteristics until the flow angle approaches within 10 deg of the maximum Prandtl-Meyer expansion angle. Improved treatment of the noncontinuum flow within the nozzle (near the lip) is postulated to have minor effects on the accuracy of this continuum method in comparison with translational nonequilibrium in the external flowfield. Possible treatment of translational nonequilibrium by the method of characteristics is discussed.

Guernsey, C. S.↗
NASA NTRS2000

Comparison of Fatigue Life Estimation Using Equivalent Linearization and Time Domain Simulation Methods

The Monte Carlo simulation method in conjunction with the finite element large deflection modal formulation are used to estimate fatigue life of aircraft panels subjected to stationary Gaussian band-limited white-noise excitations. Ten loading cases varying from 106 dB to 160 dB OASPL with bandwidth 1024 Hz are considered. For each load case, response statistics are obtained from an ensemble of 10 response time histories. The finite element nonlinear modal procedure yields time histories, probability density functions (PDF), power spectral densities and higher statistical moments of the maximum deflection and stress/strain. The method of moments of PSD with Dirlik's approach is employed to estimate the panel fatigue life.

Mei, Chuh↗
NASA NTRS1991

Application of the DSMC method to hypersonic flow about a delta wing

The direct simulation Monte Carlo method was applied to a hypersonic rarefied flow about a delta wing at angle of incidence. Three-dimensional calculations were performed using a single-species gas model that accounts for energy exchange with the internal modes. Computed results were compared with the available experimental data obtained in a nitrogen wind-tunnel experiment. Discrepancies were shown to exist between the experimental and calculated results. These differences may be caused by the inadequately defined freestream conditions.

Celenligil, M. C.↗
NASA NTRS1999

Vertical Photon Transport in Cloud Remote Sensing Problems

Photon transport in plane-parallel, vertically inhomogeneous clouds is investigated and applied to cloud remote sensing techniques that use solar reflectance or transmittance measurements for retrieving droplet effective radius. Transport is couched in terms of weighting functions which approximate the relative contribution of individual layers to the overall retrieval. Two vertical weightings are investigated, including one based on the average number of scatterings encountered by reflected and transmitted photons in any given layer. A simpler vertical weighting based on the maximum penetration of reflected photons proves useful for solar reflectance measurements. These weighting functions are highly dependent on droplet absorption and solar/viewing geometry. A superposition technique, using adding/doubling radiative transfer procedures, is derived to accurately determine both weightings, avoiding time consuming Monte Carlo methods. Superposition calculations are made for a variety of geometries and cloud models, and selected results are compared with Monte Carlo calculations. Effective radius retrievals from modeled vertically inhomogeneous liquid water clouds are then made using the standard near-infrared bands, and compared with size estimates based on the proposed weighting functions. Agreement between the two methods is generally within several tenths of a micrometer, much better than expected retrieval accuracy. Though the emphasis is on photon transport in clouds, the derived weightings can be applied to any multiple scattering plane-parallel radiative transfer problem, including arbitrary combinations of cloud, aerosol, and gas layers.

Platnick, S.↗
NASA NTRS2004

Radiative Transfer Modeling of a Large Pool Fire by Discrete Ordinates, Discrete Transfer, Ray Tracing, Monte Carlo and Moment Methods

Five computational methods for solution of the radiative transfer equation in an absorbing-emitting and non-scattering gray medium were compared on a 2 m JP-8 pool fire. The temperature and absorption coefficient fields were taken from a synthetic fire due to the lack of a complete set of experimental data for fires of this size. These quantities were generated by a code that has been shown to agree well with the limited quantity of relevant data in the literature. Reference solutions to the governing equation were determined using the Monte Carlo method and a ray tracing scheme with high angular resolution. Solutions using the discrete transfer method, the discrete ordinate method (DOM) with both S(sub 4) and LC(sub 11) quadratures, and moment model using the M(sub 1) closure were compared to the reference solutions in both isotropic and anisotropic regions of the computational domain. DOM LC(sub 11) is shown to be the more accurate than the commonly used S(sub 4) quadrature technique, especially in anisotropic regions of the fire domain. This represents the first study where the M(sub 1) method was applied to a combustion problem occurring in a complex three-dimensional geometry. The M(sub 1) results agree well with other solution techniques, which is encouraging for future applications to similar problems since it is computationally the least expensive solution technique. Moreover, M(sub 1) results are comparable to DOM S(sub 4).

Jensen, K. A.↗