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 631 records · Page 35

Current Status on the use of Parallel Computing in Turbulent Reacting Flow Computations Involving Sprays, Monte Carlo PDF and Unstructured Grids

The state of the art in multidimensional combustor modeling as evidenced by the level of sophistication employed in terms of modeling and numerical accuracy considerations, is also dictated by the available computer memory and turnaround times afforded by present-day computers. With the aim of advancing the current multi-dimensional computational tools used in the design of advanced technology combustors, a solution procedure is developed that combines the novelty of the coupled CFD/spray/scalar Monte Carlo PDF (Probability Density Function) computations on unstructured grids with the ability to run on parallel architectures. In this approach, the mean gas-phase velocity and turbulence fields are determined from a standard turbulence model, the joint composition of species and enthalpy from the solution of a modeled PDF transport equation, and a Lagrangian-based dilute spray model is used for the liquid-phase representation. The gas-turbine combustor flows are often characterized by a complex interaction between various physical processes associated with the interaction between the liquid and gas phases, droplet vaporization, turbulent mixing, heat release associated with chemical kinetics, radiative heat transfer associated with highly absorbing and radiating species, among others. The rate controlling processes often interact with each other at various disparate time 1 and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and liquid phase evaporation in many practical combustion devices.

Raju, M. S.↗

Uncertainty Quantification Methodology for Sonic Boom Loudness of a Low-Boom Supersonic Concept

This paper presents a comprehensive analysis of uncertainty quantification for sonic boom loudness metrics associated with the X-59 aircraft during its acoustic validation phase. The study focuses on employing advanced methodologies to characterize uncertainties in key sonic boom parameters predicted from a database derived from a computational fluid dynamics (CFD) solver. These methodologies include creating a polynomial chaos expansion (PCE) surrogate model, which is then used to run a Monte Carlo analysis to obtain 95% uncertainty intervals for sonic boom loudness parameters across the full carpet. Atmospheric uncertainties are rigorously considered, leveraging weather models based on historical data near Edwards Air Force Base to simulate potential flight conditions. The findings include preliminary uncertainty results for the full carpet during the validation phase of the Quesst mission. Additionally, a sensitivity study reveals that the primary sources of uncertainty are humidity and cruise weight.

Laura M. White↗

Panchromatic Imaging of a Transitional Disk: The Disk of GM AUR in Optical and FUV Scattered Light

We have imaged GM Aurigae with the Hubble Space Telescope, detected its disk in scattered light at 1400 and 1650 A, and compared these with observations at 3300 A, 5550 A, 1.1 micron, and 1.6 micron. The scattered light increases at shorter wavelengths. The radial surface brightness profile at 3300 A shows no evidence of the 24 au radius cavity that has been previously observed in submillimeter observations. Comparison with dust grain opacity models indicates that the surface of the entire disk is populated with submicron grains. We have compiled a spectral energy distribution from 0.1 micron to 1 mm and used it to constrain a model of the star + disk system that includes the submillimeter cavity using the Monte Carlo radiative transfer code by Barbara Whitney. The best-fit model image indicates that the cavity should be detectable in the F330W bandpass if the cavity has been cleared of both large and small dust grains, but we do not detect it. The lack of an observed cavity can be explained by the presence of submicron grains interior to the submillimeter cavity wall. We suggest one explanation for this that could be due to a planet of mass less than 9M(sub J) interior to 24 au. A unique cylindrical structure is detected in the far- UV data from the Advanced Camera for Surveys/Solar Blind Channel. It is aligned along the system semiminor axis, but does not resemble an accretion-driven jet. The structure is limb brightened and extends 190 +/- 35 au above the disk midplane. The inner radius of the limb brightening is 40 +/- 10 au, just beyond the submillimeter cavity wall.

circumstellar matter aEuro" protoplanetary disks a↗

Proposal for a Sample Pilot Response Model

Regulations to establish operational and performance requirements for unmanned aircraft systems (UAS) are being developed by a consortium of government, industry and academic institutions. Those requirements will apply to the new detect and avoid (DAA) systems and other equipment necessary to integrate UAS with the National Airspace System (NAS) and are determined according to their contribution to the overall safety case for such an integration. In order to perform end-to-end verification and validation of DAA requirements, a committee accepted pilot response model is needed to emulate UAS pilots in closed-loop Monte Carlo fast-time simulations. This briefing focuses on background of the pilot model used in previous Airspace Concept Evaluation System studies, identification of strengths and weaknesses in that model, and proposals on model improvements. The briefing also looks to build consensus around a common architecture and assumptions for modeling the pilot for UAS DAA systems.

Santiago, Confesor↗

Monte Carlo Calculations of Polarized Microwave Radiation Emerging from Cloud Structures

The last decade has seen tremendous growth in cloud dynamical and microphysical models that are able to simulate storms and storm systems with very high spatial resolution, typically of the order of a few kilometers. The fairly realistic distributions of cloud and hydrometeor properties that these models generate has in turn led to a renewed interest in the three-dimensional microwave radiative transfer modeling needed to understand the effect of cloud and rainfall inhomogeneities upon microwave observations. Monte Carlo methods, and particularly backwards Monte Carlo methods have shown themselves to be very desirable due to the quick convergence of the solutions. Unfortunately, backwards Monte Carlo methods are not well suited to treat polarized radiation. This study reviews the existing Monte Carlo methods and presents a new polarized Monte Carlo radiative transfer code. The code is based on a forward scheme but uses aliasing techniques to keep the computational requirements equivalent to the backwards solution. Radiative transfer computations have been performed using a microphysical-dynamical cloud model and the results are presented together with the algorithm description.

Kummerow, Christian↗

Implications of a quadratic stream definition in radiative transfer theory.

An explicit definition of the radiation-stream concept is stated and applied to approximate the integro-differential equation of radiative transfer with a set of twelve coupled differential equations. Computational efficiency is enhanced by distributing the corresponding streams in three-dimensional space in a totally symmetric way. Polarization is then incorporated in this model. A computer program based on the model is briefly compared with a Monte Carlo program for simulation of horizon scans of the earth's atmosphere. It is found to be considerably faster.

Whitney, C.↗

Solar wind sputtering effects in the Martian atmosphere

A Monte Carlo simulation of the sputtering of the upper atmosphere of Mars by the solar wind was performed. The calculated sputtering yields imply loss rates (molecules/cm square - sec escaping the planet) for carbon dioxide, carbon, and oxygen of R(CO2) = 2.6 X 1000000/cm square - sec, R(C) = 6.6 X 1000000/cm square - sec, and R(O) = 7.7 X 1000000/cm - sec. The total mass loss by sputtering is only about 10% of that due to chemical and photo-chemical processes, but sputtering provides a major exospheric sink for carbon. The erosion process described here preferentially removes the lighter components of the atmosphere. Calculations based on a Monte Carlo simulation suggest that for a model atmosphere, 97% of the N2 and 33% of the CO2 originally present may have been sputtered away over 4.5 X 10 to the 9th power y. In the same length of time the (15)N/(14)N isotopic ratio for the bulk atmosphere would have increased by a factor 1.7.

Watson, C. C.↗

Statistical aspects of carbon fiber risk assessment modeling

The probabilistic and statistical aspects of the carbon fiber risk assessment modeling of fire accidents involving commercial aircraft are examined. Three major sources of uncertainty in the modeling effort are identified. These are: (1) imprecise knowledge in establishing the model; (2) parameter estimation; and (3)Monte Carlo sampling error. All three sources of uncertainty are treated and statistical procedures are utilized and/or developed to control them wherever possible.

Gross, D.↗

Solar reflection from interacting and shadowing cloud elements

An attempt is made to relate the reflection from a broken cloud field to that from a plane parallel cloud of the same optical properties by using Monte Carlo simulations for a very simple model of an array of regularly spaced, infinitely long bar clouds. Reflection from the array is represented by an effective cloud fraction, and it is shown that for many cases this may be obtained by treating the cloud elements to be noninteracting azimuthally symmetric reflectors. The effective cloud fraction is a function of solar zenith angle, the aspect ratio of the individual cloud elements, and the normal or zenith cloud fraction of the array. Model computations for reflection from an array of randomly distributed identical cylinders are also presented as an example of a possible application of this technique.

HARSHVARDHAN↗

High energy gamma-rays and hadrons at Mount Fuji

The energy spectra of high energy gamma-rays and hadrons were obtained by the emulsion chamber with 40 c.u. thickness at Mt. Fuji (3750 m). These results are compared with the Monte Carlo calculation based on the same model which is used in a family analysis. Our data are compatible with the model of heavy-enriched primary and scaling in the fragmentation region.

Amenomori, M.↗

Correlation of high energy muons with primary composition in extensive air shower

An experimental investigation of high energy muons above 200 GeV in extensive air showers has been made for studying high energy interaction and primary composition of cosmic rays of energies in the range 10 to the 14th power approx. 10 to the 15th power eV. The muon energies are estimated from the burst sizes initiated by the muons in the rock, which are measured by four layers of proportional counters, each of area 5 x 2.6 sq m, placed at 30 m.w.e. deep, Funasaka tunnel vertically below the air shower array. These results are compared with Monte Carlo simulations based on the scaling model and the fireball model for two primary compositions, all proton and mixed.

Chou, C.↗

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.↗

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.↗

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.↗

An improved correction algorithm for number density measurements made with the forward scattering spectrometer probe

A correction factor to the number density measured by the Forward Scattering Spectrometer Probe (FSSP) which compensates for dead time and coincidence errors was determined by calculating the probabilities of, and the average number of particles in, the six possible types of dead time and coincidence events. These probabilities and averages were calculated by means of a probabilistic model based on Poisson statistics. A Monte Carlo computer simulation of the FSSP operation was also carried out and the number density correction factor was compared with the Monte Carlo data. For an actual number density of 2000/cu cm, it was found that the measured number density was of the order of 300/cu cm.

Lock, James A.↗