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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Monte Carlo solution of Boltzmann equation for a simple model of highly nonequilibrium diatomic gases: Translational rotational energy relaxation

The semiclassical transition probability was incorporated in the simulation for energy exchange between rotational and translational energy. The results provide details on the fundamental mechanisms of gas kinetics where analytical methods were impractical. The validity of the local Maxwellian assumption and relaxation time, rotational-translational energy transition, and a velocity analysis of the inelastic collision were discussed in detail.

Yoshikawa, K. K.↗
NASA NTRS1990

Sampling errors for satellite-derived tropical rainfall - Monte Carlo study using a space-time stochastic model

Estimates of monthly average rainfall based on satellite observations from a low earth orbit will differ from the true monthly average because the satellite observes a given area only intermittently. This sampling error inherent in satellite monitoring of rainfall would occur even if the satellite instruments could measure rainfall perfectly. The size of this error is estimated for a satellite system being studied at NASA, the Tropical Rainfall Measuring Mission (TRMM). First, the statistical description of rainfall on scales from 1 to 1000 km is examined in detail, based on rainfall data from the Global Atmospheric Research Project Atlantic Tropical Experiment (GATE). A TRMM-like satellite is flown over a two-dimensional time-evolving simulation of rainfall using a stochastic model with statistics tuned to agree with GATE statistics. The distribution of sampling errors found from many months of simulated observations is found to be nearly normal, even though the distribution of area-averaged rainfall is far from normal. For a range of orbits likely to be employed in TRMM, sampling error is found to be less than 10 percent of the mean for rainfall averaged over a 500 x 500 sq km area.

Bell, Thomas L.↗
NASA NTRS1976

Secondary cosmic-ray electrons and positrons from 1 to 100 GeV in the upper atmosphere and interstellar space, and interpretation of a recent positron flux measurement

Secondary electron and positron fluxes generated in interstellar space and in the atmosphere from the decays of pions and kaons in inelastic nuclear interactions are calculated by Monte Carlo techniques for lepton energies in the range from 1 to 100 GeV and an assumed thickness of 10 g/sq cm or less for the interstellar or atmospheric material. A simple and accurate analytical model which summarizes the Monte Carlo results and identifies the essential parameters involved is developed and used to interpret a previous positron measurement. It is found that the thickness of interstellar and source material is about 4.3 g/sq cm for cosmic-ray positrons with energies exceeding 4 GeV, a result that is difficult to reconcile with recently proposed two-containment-volume propagation models which predict a thickness of 1.8 g/sq cm for the same energies on the basis of the energy dependence of the measured (Li+Be+B)/(C+O) ratio. It is shown that single-containment-volume (galactic) models invoking an energy-dependent leakage lifetime are compatible with the positron data, but lack a mechanism to explain the energy dependence.

Orth, C. D.↗
NASA NTRS2014

Object-Oriented/Data-Oriented Design of a Direct Simulation Monte Carlo Algorithm

Over the past decade, there has been much progress towards improved phenomenological modeling and algorithmic updates for the direct simulation Monte Carlo (DSMC) method, which provides a probabilistic physical simulation of gas Rows. These improvements have largely been based on the work of the originator of the DSMC method, Graeme Bird. Of primary importance are improved chemistry, internal energy, and physics modeling and a reduction in time to solution. These allow for an expanded range of possible solutions In altitude and velocity space. NASA's current production code, the DSMC Analysis Code (DAC), is well-established and based on Bird's 1994 algorithms written in Fortran 77 and has proven difficult to upgrade. A new DSMC code is being developed in the C++ programming language using object-oriented and data-oriented design paradigms to facilitate the inclusion of the recent improvements and future development activities. The development efforts on the new code, the Multiphysics Algorithm with Particles (MAP), are described, and performance comparisons are made with DAC.

Liechty, Derek S.↗
NASA NTRS1998

Influence of Subpixel Scale Cloud Top Structure on Reflectances from Overcast Stratiform Cloud Layers

Recent observational studies have shown that satellite retrievals of cloud optical depth based on plane-parallel model theory suffer from systematic biases that depend on viewing geometry, even when observations are restricted to overcast marine stratus layers, arguably the closest to plane parallel in nature. At moderate to low sun elevations, the plane-parallel model significantly overestimates the reflectance dependence on view angle in the forward-scattering direction but shows a similar dependence in the backscattering direction. Theoretical simulations are performed that show that the likely cause for this discrepancy is because the plane-parallel model assumption does not account for subpixel, scale variations in cloud-top height (i.e., "cloud bumps"). Monte Carlo simulation, comparing ID model radiances to radiances from overcast cloud field with 1) cloud-top height variation, but constant cloud volume extinction; 2) flat tops but horizontal variations in cloud volume extinction; and 3) variations in both cloud top height and cloud extinction are performed over a approximately equal to 4 km x 4 km domain (roughly the size of an individual GAC AVHRR pixel). The comparisons show that when cloud-top height variations are included, departures from 1D theory are remarkably similar (qualitatively) to those obtained observationally. In contrast, when clouds are assumed flat and only cloud extinction is variable, reflectance differences are much smaller and do not show any view-angle dependence. When both cloud-top height and cloud extinction variations are included, however, large increases in cloud extinction variability can enhance reflectance difference. The reason 3D-1D reflectance differences are more sensitive to cloud-top height variations in the forward-scattering direction (at moderate to low, sun elevations) is because photons leaving the cloud field in that direction experience fewer scattering events (low-order scattering) and are restricted to the topmost portions of the cloud. While reflectance deviations from 1D theory are much larger for bumpy clouds than for flat clouds with variable cloud extinction, differences in cloud albedo are comparable for these two cases.

Loeb, N. G.↗
NASA NTRS1993

Gravitational microlensing of gamma-ray bursts

A Monte Carlo code is developed to calculate gravitational microlensing in three dimensions when the lensing optical depth is low or moderate (not greater than 0.25). The code calculates positions of microimages and time delays between the microimages. The majority of lensed gamma-ray bursts should show a simple double-burst structure, as predicted by a single point mass lens model. A small fraction should show complicated multiple events due to the collective effects of several point masses (black holes). Cosmological models with a significant fraction of mass density in massive compact objects can be tested by searching for microlensing events in the current BATSE data. Our catalog generated by 10,000 Monte Carlo models is accessible through the computer network. The catalog can be used to take realistic selection effects into account.

Mao, Shude↗
NASA NTRS2017

Development of a Coherent Bistatic Vegetation Model for Signal of Opportunity Applications at VHF UHF-Bands

A coherent bistatic vegetation scattering model, based on a Monte Carlo simulation, is being developed to simulate polarimetric bi-static reflectometry at VHF/UHF-bands (240-270 MHz). The model is aimed to assess the value of geostationary satellite signals of opportunity to enable estimation of the Earth's biomass and root-zone soil moisture. An expression for bistatic scattering from a vegetation canopy is derived for the practical case of a ground-based/low altitude platforms with passive receivers overlooking vegetation. Using analytical wave theory in conjunction with distorted Born approximation (DBA), the transmit and receive antennas effects (i.e., polarization, orientation, height, etc.) are explicitly accounted for. Both the coherency nature of the model (joint phase and amplitude information) and the explicit account of system parameters (antenna, altitude, polarization, etc) enable one to perform various beamforming techniques to evaluate realistic deployment configurations. In this paper, several test scenarios will be presented and the results will be evaluated for feasibility for future biomass and root-zone soil moisture application using geostationary communication satellite signals of opportunity at low frequencies.

UFH↗
NASA NTRS1993

Numerical experiments with model monophyletic and paraphyletic taxa

The problem of how accurately paraphyletic taxa versus monophyletic (i.e., holophyletic) groups (clades) capture underlying species patterns of diversity and extinction is explored with Monte Carlo simulations. Phylogenies are modeled as stochastic trees. Paraphyletic taxa are defined in an arbitrary manner by randomly choosing progenitors and clustering all descendants not belonging to other taxa. These taxa are then examined to determine which are clades, and the remaining paraphyletic groups are dissected to discover monophyletic subgroups. Comparisons of diversity patterns and extinction rates between modeled taxa and lineages indicate that paraphyletic groups can adequately capture lineage information under a variety of conditions of diversification and mass extinction. This suggests that these groups constitute more than mere "taxonomic noise" in this context. But, strictly monophyletic groups perform somewhat better, especially with regard to mass extinctions. However, when low levels of paleontologic sampling are simulated, the veracity of clades deteriorates, especially with respect to diversity, and modeled paraphyletic taxa often capture more information about underlying lineages. Thus, for studies of diversity and taxic evolution in the fossil record, traditional paleontologic genera and families need not be rejected in favor of cladistically-defined taxa.

Non-NASA Center↗
NASA NTRS1976

Calculations of polarization and radiance in the atmosphere

Two different methods for the solution of radiative transfer problems, matrix operator and Monte Carlo, are discussed. As an example of the use of the method, results are given for the radiance and polarization of the radiation scattered from haze layers as well as from models of the real atmosphere. The variation of the radiance, polarization, and ellipticity with the aerosol amount in the real atmosphere model is presented. The Monte Carlo method is applied to the problem of calculation of the radiance and polarization of the photons in the atmosphere when there is an ocean as the lower surface. It is found that the turbidity of the ocean can be deduced from the upwelling radiance measured at some height in the atmosphere. The downwelling radiation just beneath the ocean surface is elliptically polarized at those angles where it is derived from the total internal reflection of the upwelling radiation at the ocean surface.

Plass, G. N.↗
NASA NTRSDate not supplied

Accelerating SNP Development with Radiation Source Approximations for Far-Field Environment Modeling

The successful development and deployment of space nuclear power and propulsion technology requires a high-fidelity, efficient means of modeling ex-core radiation fields. Monte Carlo particle transport codes provide enable the highest level of fidelity in radiation field modeling but are inefficient for modeling ex-core radiation without modification. This paper details the development of an approach to ex-core radiation field modeling that maintains the fidelity of a Monte Carlo transport approach without sacrificing computational efficiency. This approach is based on recent improvements to an existing Monte Carlo transport acceleration technique known as surface source banking. As the transport of many particles is required to ensure adequate uncertainty in far field transport, analytical reconstruction techniques are employed to represent surface source banks as sampleable series of distributions in particle phase from which an arbitrarily large number of particles can be generated and simulated. All steps of this approach are incorporated into an AMA model to test both the underlying mathematics of the analytical source reconstruction process, as well as the relevance of the technique as a whole towards SNP applications. These tests demonstrated that this approach is capable of accelerating far-field radiation modeling in SNP-relevant scenarios by factors of at least nine over purely eigenvalue-based scenarios without significantly sacrificing simulation accuracy.

Surface Source Banking↗
NASA NTRS2005

A Stochastic Cratering Model for Asteroid Surfaces

The observed cratering records on asteroid surfaces (four so far: Gaspra, Ida, Mathilde, and Eros [1-4]) provide us with important clues to their past bombardment histories. Previous efforts toward interpreting these records have led to two basic modeling styles for reproducing the statistics of the observed crater populations. The first, and most direct, method is to use Monte Carlo techniques [5] to stochastically populate a matrix-model test surface with craters as a function of time [6,7]. The second method is to use a more general, parameterized approach to duplicate the statistics of the observed crater population [8,9]. In both methods, several factors must be included beyond the simple superposing of circular features: (1) crater erosion by subsequent impacts, (2) infilling of craters by impact ejecta, and (3) crater degradation and era- sure due to the seismic effects of subsequent impacts. Here we present an updated Monte Carlo (stochastic) modeling approach, designed specifically with small- to medium-sized asteroids in mind.

Richardson, J. E.↗
NASA NTRS1998

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗
NASA NTRS2010

Post2 End-to-End Descent and Landing Simulation for ALHAT Design Analysis Cycle 2

The ALHAT project is an agency-level program involving NASA centers, academia, and industry, with a primary goal to develop a safe, autonomous, precision-landing system for robotic and crew-piloted lunar and planetary descent vehicles. POST2 is used as the 6DOF descent and landing trajectory simulation for determining integrated system performance of ALHAT landing-system models and lunar environment models. This paper presents updates in the development of the ALHAT POST2 simulation, as well as preliminary system performance analysis for ALDAC-2 used for the testing and assessment of ALHAT system models. The ALDAC-2 POST2 Monte Carlo simulation results have been generated and focus on HRN model performance with the fully integrated system, as well performance improvements of AGNC and TSAR model since the previous design analysis cycle

Davis, Jody L.↗
NASA NTRS1981

Monte Carlo simulation of charged particles upstream of the earth's bow shock

A simple model of the earth's bow shock is presented and compared with ISEE-1 observations of diffuse ions upstream of the shock. Ion behavior is modeled in a test particle simulation using Monte Carlo techniques for a plane parallel shock at which cold solar wind ions are heated sufficiently to allow diffusion back upstream. Fitting of the model to the observations indicates that the diffusion coefficient must increase faster than linearly with velocity, and provides an upper limit to the acceleration efficiency defined as the proportion of incoming solar wind ions accelerated and observed upstream. Results suggest that astrophysical shocks are capable of accelerating ions directly from the thermal plasma without the need for a distinct preenergizing injection stage.

Ellison, D. C.↗
NASA NTRS1994

Radiative interactions in chemically reacting compressible nozzle flows using Monte Carlo simulations

The two-dimensional spatially elliptic Navier-Stokes equations have been used to investigate the radiative interactions in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The radiative heat transfer term in the energy equation is simulated using the Monte Carlo method (MCM). The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The spectral correlation has been considered in the Monte Carlo formulations. Results obtained demonstrate that the effect of radiation on the flow field is minimal but its effect on the wall heat transfer is significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and the nozzle size on the radiative and conductive wall fluxes.

Liu, J.↗
NASA NTRS2004

The FLUKA code for space applications: recent developments

The FLUKA Monte Carlo transport code is widely used for fundamental research, radioprotection and dosimetry, hybrid nuclear energy system and cosmic ray calculations. The validity of its physical models has been benchmarked against a variety of experimental data over a wide range of energies, ranging from accelerator data to cosmic ray showers in the earth atmosphere. The code is presently undergoing several developments in order to better fit the needs of space applications. The generation of particle spectra according to up-to-date cosmic ray data as well as the effect of the solar and geomagnetic modulation have been implemented and already successfully applied to a variety of problems. The implementation of suitable models for heavy ion nuclear interactions has reached an operational stage. At medium/high energy FLUKA is using the DPMJET model. The major task of incorporating heavy ion interactions from a few GeV/n down to the threshold for inelastic collisions is also progressing and promising results have been obtained using a modified version of the RQMD-2.4 code. This interim solution is now fully operational, while waiting for the development of new models based on the FLUKA hadron-nucleus interaction code, a newly developed QMD code, and the implementation of the Boltzmann master equation theory for low energy ion interactions. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

Non-NASA Center↗
NASA NTRS1972

Monte Carlo inversion of seismic data

The analytic solution to the linear inverse problem provides estimates of the uncertainty of the solution in terms of standard deviations of corrections to a particular solution, resolution of parameter adjustments, and information distribution among the observations. It is shown that Monte Carlo inversion, when properly executed, can provide all the same kinds of information for nonlinear problems. Proper execution requires a relatively uniform sampling of all possible models. The expense of performing Monte Carlo inversion generally requires strategies to improve the probability of finding passing models. Such strategies can lead to a very strong bias in the distribution of models examined unless great care is taken in their application.

Wiggins, R. A.↗
NASA NTRS2008

Stochastic Analysis of Orbital Lifetimes of Spacecraft

A document discusses (1) a Monte-Carlo-based methodology for probabilistic prediction and analysis of orbital lifetimes of spacecraft and (2) Orbital Lifetime Monte Carlo (OLMC)--a Fortran computer program, consisting of a previously developed long-term orbit-propagator integrated with a Monte Carlo engine. OLMC enables modeling of variances of key physical parameters that affect orbital lifetimes through the use of probability distributions. These parameters include altitude, speed, and flight-path angle at insertion into orbit; solar flux; and launch delays. The products of OLMC are predicted lifetimes (durations above specified minimum altitudes) for the number of user-specified cases. Histograms generated from such predictions can be used to determine the probabilities that spacecraft will satisfy lifetime requirements. The document discusses uncertainties that affect modeling of orbital lifetimes. Issues of repeatability, smoothness of distributions, and code run time are considered for the purpose of establishing values of code-specific parameters and number of Monte Carlo runs. Results from test cases are interpreted as demonstrating that solar-flux predictions are primary sources of variations in predicted lifetimes. Therefore, it is concluded, multiple sets of predictions should be utilized to fully characterize the lifetime range of a spacecraft.

Sasamoto, Washito↗