Search NASA⌕ Search

SEARCH · Search NASA

Results for “['Monte”

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 127 records · Page 7

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

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

Simulations of particle acceleration in parallel shocks: Direct comparison between Monte Carlo and one-dimensional hybrid codes

We have made a direct comparison between two different computer simulations of a plane, parallel, collisionless shock including particle acceleration to energies typical of those of diffuse ions observed at the earth bow shock. Despite the fact that the one-dimensional hybrid and Monte Carlo techniques employ entirely different algorithms, they give surprisingly close agreement in the overall shapes of the complete distribution functions for protons as well as heavier ions. Both methods show that energetic ions emerge smoothly from the background thermal plasma with approximately the same relative injection rate and that the fraction of the incoming plasma's energy flux that is converted into downstream enthalpy flux of the accelerated population (i.e., the acceleration efficiency) is similar in the two cases. The fraction of the downstream proton distribution made up of superthermal particles is quite large, with at least 10% of the energy flux going into protons with energies above 10 keV. In addition, an upstream precursor, produced by backstreaming energetic particles, is present in both shocks, although the Monte Carlo precursor is considerably longer than that produced in the hybrid shock. These results offer convincing evidence that, at least in these ways, the two simulations are consistent in their description of parallel shock structure and particle acceleration, and they lay the groundwork for development of shock models employing a combination of both methods.

Ellison, Donald C.↗

A Monte Carlo simulation of the effect of ion self-collisions on the ion velocity distribution function in the high-latitude F-region

Non-Maxwellian ion velocity distribution functions have been theoretically predicted and confirmed by observations, to occur at high latitudes. These distributions deviate from Maxwellian due to the combined effect of the E x B drift and ion-neutral collisions. At high altitude and/or for solar maximum conditions, the ion-to-neutral density ratio increases and, hence, the role of ion self-collisions becomes appreciable. A Monte Carlo simulation was used to investigate the behavior of O(+) ions that are E x B-drifting through a background of neutral O, with the effect of O(+) (Coulomb) self-collisions included. Wide ranges of the ion-to-neutral density ratio n(sub i)/n(sub n) and the electrostatic field E were considered in order to investigate the change of ion behavior with solar cycle and with altitude. For low altitudes and/or solar minimum (n(sub i)/n(sub n) less than or equal to 10(exp -5)), the effect of self-collisions is negligible. For higher values of n(sub i)/n(sub n), the effect of self-collisions becomes significant and, hence, the non-Maxwellian features of the O(+) distribution are reduced. The Monte Carlo results were compared to those that used simplified collision models in order to assess their validity. In general, the simple collision models tend to be more accurate for low E and for high n(sub i)/n(sub n).

Barghouthi, I. A.↗

Parallel Monte Carlo Simulation for control system design

The research during the 1993/94 academic year addressed the design of parallel algorithms for stochastic robustness synthesis (SRS). SRS uses Monte Carlo simulation to compute probabilities of system instability and other design-metric violations. The probabilities form a cost function which is used by a genetic algorithm (GA). The GA searches for the stochastic optimal controller. The existing sequential algorithm was analyzed and modified to execute in a distributed environment. For this, parallel approaches to Monte Carlo simulation and genetic algorithms were investigated. Initial empirical results are available for the KSR1.

Schubert, Wolfgang M.↗

A Large-Particle Monte Carlo Code for Simulating Non-Linear High-Energy Processes Near Compact Objects

High-energy radiation processes in compact cosmic objects are often expected to have a strongly non-linear behavior. Such behavior is shown, for example, by electron-positron pair cascades and the time evolution of relativistic proton distributions in dense radiation fields. Three independent techniques have been developed to simulate these non-linear problems: the kinetic equation approach; the phase-space density (PSD) Monte Carlo method; and the large-particle (LP) Monte Carlo method. In this paper, we present the latest version of the LP method and compare it with the other methods. The efficiency of the method in treating geometrically complex problems is illustrated by showing results of simulations of 1D, 2D and 3D systems. The method is shown to be powerful enough to treat non-spherical geometries, including such effects as bulk motion of the background plasma, reflection of radiation from cold matter, and anisotropic distributions of radiating particles. It can therefore be applied to simulate high-energy processes in such astrophysical systems as accretion discs with coronae, relativistic jets, pulsar magnetospheres and gamma-ray bursts.

Stern, Boris E.↗

Development of a Space Radiation Monte Carlo Computer Simulation

The ultimate purpose of this effort is to undertake the development of a computer simulation of the radiation environment encountered in spacecraft which is based upon the Monte Carlo technique. The current plan is to adapt and modify a Monte Carlo calculation code known as FLUKA, which is presently used in high energy and heavy ion physics, to simulate the radiation environment present in spacecraft during missions. The initial effort would be directed towards modeling the MIR and Space Shuttle environments, but the long range goal is to develop a program for the accurate prediction of the radiation environment likely to be encountered on future planned endeavors such as the Space Station, a Lunar Return Mission, or a Mars Mission. The longer the mission, especially those which will not have the shielding protection of the earth's magnetic field, the more critical the radiation threat will be. The ultimate goal of this research is to produce a code that will be useful to mission planners and engineers who need to have detailed projections of radiation exposures at specified locations within the spacecraft and for either specific times during the mission or integrated over the entire mission. In concert with the development of the simulation, it is desired to integrate it with a state-of-the-art interactive 3-D graphics-capable analysis package known as ROOT, to allow easy investigation and visualization of the results. The efforts reported on here include the initial development of the program and the demonstration of the efficacy of the technique through a model simulation of the MIR environment. This information was used to write a proposal to obtain follow-on permanent funding for this project.

Pinsky, Lawrence S.↗

Study of the Transition Flow Regime using Monte Carlo Methods

This NASA Cooperative Agreement presents a study of the Transition Flow Regime Using Monte Carlo Methods. The topics included in this final report are: 1) New Direct Simulation Monte Carlo (DSMC) procedures; 2) The DS3W and DS2A Programs; 3) Papers presented; 4) Miscellaneous Applications and Program Modifications; 5) Solution of Transitional Wake Flows at Mach 10; and 6) Turbulence Modeling of Shock-Dominated Fows with a k-Enstrophy Formulation.

Hassan, H. A.↗

Monte Carlo Radiative Transfer Modeling of Lightning Observed in Galileo Images of Jupiter

We study lightning on Jupiter and the clouds illuminated by the lightning using images taken by the Galileo orbiter. The Galileo images have a resolution of 25 km/pixel and axe able to resolve the shape of the single lightning spots in the images, which have full widths at half the maximum intensity in the range of 90-160 km. We compare the measured lightning flash images with simulated images produced by our ED Monte Carlo light-scattering model. The model calculates Monte Carlo scattering of photons in a ED opacity distribution. During each scattering event, light is partially absorbed. The new direction of the photon after scattering is chosen according to a Henyey-Greenstein phase function. An image from each direction is produced by accumulating photons emerging from the cloud in a small range (bins) of emission angles. Lightning bolts are modeled either as points or vertical lines. Our results suggest that some of the observed scattering patterns axe produced in a 3-D cloud rather than in a plane-parallel cloud layer. Lightning is estimated to occur at least as deep as the bottom of the expected water cloud. For the six cases studied, we find that the clouds above the lightning are optically thick (tau > 5). Jovian flashes are more regular and circular than the largest terrestrial flashes observed from space. On Jupiter there is nothing equivalent to the 30-40-km horizontal flashes which axe seen on Earth.

Dyudine, U. A.↗

Oversight and Development of a Community Monte Carlo Radiative Transfer Model

Under this grant we have developed a Monte Carlo radiative transfer code that will act as the nucleus for the 13RC Community Monte Carlo Model. All code is written in ANSI-compliant Fortran-95. Many modules define public type and procedures to manipulate them, but do not allow access to the type's internal components. This allows each module to do its own exhaustive error checking up-front, then proceed in a streamlined way. Many modules can read and write the state of their objects to persistent files. The code has been tested on a Macintosh running OS 10.2.4 and the Absoft Fortran compiler, and on Sun UltraSparcs running Solaris 5.8 and Forte V8 compilers. The code exposed bugs in the Intel Fortran Compiler (ifc) on the 13RC Linux host, and we are waiting for a resolution of these bugs before finished the port. the code base is under CVS versions control. The code base consists of the core code (nine modules providing the infrastructure), example integrators, and a suite of utilities and examples.

Source record↗

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.↗

Monte Carlo Simulation of Sudden Death Bearing Testing

Monte Carlo simulations combined with sudden death testing were used to compare resultant bearing lives to the calculated hearing life and the cumulative test time and calendar time relative to sequential and censored sequential testing. A total of 30 960 virtual 50-mm bore deep-groove ball bearings were evaluated in 33 different sudden death test configurations comprising 36, 72, and 144 bearings each. Variations in both life and Weibull slope were a function of the number of bearings failed independent of the test method used and not the total number of bearings tested. Variation in L10 life as a function of number of bearings failed were similar to variations in lift obtained from sequentially failed real bearings and from Monte Carlo (virtual) testing of entire populations. Reductions up to 40 percent in bearing test time and calendar time can be achieved by testing to failure or the L(sub 50) life and terminating all testing when the last of the predetermined bearing failures has occurred. Sudden death testing is not a more efficient method to reduce bearing test time or calendar time when compared to censored sequential testing.

Vlcek, Brian L.↗

From Simple Plane-Parallel to Complex Monte Carlo Calculations of Solar Fluxes and Radiances for Cloudy Atmospheres

In my presentation, I will describe several approximation methods with different level of complexity; they will be gradually applied to simple examples of horizontally inhomogeneous clouds. Understanding of photon horizontal transport and radiative smoothing can help to improve accuracy of the methods The accuracy of the methods will be compared with the full Monte Carlo calculations. The specifics of Monte Carlo in cloudy atmospheres will be also discussed. A special emphasis will be put on the strong forward scattering peak in the phase functions.

Marshak, Alexander↗

Monte Carlo Simulation of THz Multipliers

Schottky Barrier diode frequency multipliers are critical components in submillimeter and Thz space based earth observation systems. As the operating frequency of these multipliers has increased, the agreement between design predictions and experimental results has become poorer. The multiplier design is usually based on a nonlinear model using a form of harmonic balance and a model for the Schottky barrier diode. Conventional voltage dependent lumped element models do a poor job of predicting THz frequency performance. This paper will describe a large signal Monte Carlo simulation of Schottky barrier multipliers. The simulation is a time dependent particle field Monte Carlo simulation with ohmic and Schottky barrier boundary conditions included that has been combined with a fixed point solution for the nonlinear circuit interaction. The results in the paper will point out some important time constants in varactor operation and will describe the effects of current saturation and nonlinear resistances on multiplier operation.

East, J.↗

Monte Carlo Technique Used to Model the Degradation of Internal Spacecraft Surfaces by Atomic Oxygen

Atomic oxygen is one of the predominant constituents of Earth's upper atmosphere. It is created by the photodissociation of molecular oxygen (O2) into single O atoms by ultraviolet radiation. It is chemically very reactive because a single O atom readily combines with another O atom or with other atoms or molecules that can form a stable oxide. The effects of atomic oxygen on the external surfaces of spacecraft in low Earth orbit can have dire consequences for spacecraft life, and this is a well-known and much studied problem. Much less information is known about the effects of atomic oxygen on the internal surfaces of spacecraft. This degradation can occur when openings in components of the spacecraft exterior exist that allow the entry of atomic oxygen into regions that may not have direct atomic oxygen attack but rather scattered attack. Openings can exist because of spacecraft venting, microwave cavities, and apertures for Earth viewing, Sun sensors, or star trackers. The effects of atomic oxygen erosion of polymers interior to an aperture on a spacecraft were simulated at the NASA Glenn Research Center by using Monte Carlo computational techniques. A two-dimensional model was used to provide quantitative indications of the attenuation of atomic oxygen flux as a function of the distance into a parallel-walled cavity. The model allows the atomic oxygen arrival direction, the Maxwell Boltzman temperature, and the ram energy to be varied along with the interaction parameters of the degree of recombination upon impact with polymer or nonreactive surfaces, the initial reaction probability, the reaction probability dependence upon energy and angle of attack, degree of specularity of scattering of reactive and nonreactive surfaces, and the degree of thermal accommodation upon impact with reactive and non-reactive surfaces to be varied to allow the model to produce atomic oxygen erosion geometries that replicate actual experimental results from space. The degree of erosion of various interior locations was compared with the erosion that would occur external to the spacecraft. Results of one cavity model indicate that, at depths into a two-dimensional cavity that are equal to 10 cavity widths, the erosion on the walls of the cavity is less than that on the top surface by over 2 orders of magnitude. Wall erosion near the surface of a cavity depends on which wall is receiving direct atomic oxygen attack. However, deep in the cavity little difference is present. Testing of various cavity models such as these gives spacecraft designers an indication of the level of threat to sensitive interior surfaces for different geometries. Even though the Monte Carlo model is two-dimensional, it can be used to provide qualitative information about spacecraft openings that are three-dimensional by offering reasonable insight as to the nature of the attenuation of damage that occurs within a spacecraft in low Earth orbit. As shown, there is more erosion on the side seeing direct atomic oxygen attack until a depth of approximately 5 times the width of the opening, where the erosion is the same on both sides.

Banks, Bruce A.↗

Monte Carlo Analysis of Pion Contribution to Absorbed Dose from Galactic Cosmic Rays

Accurate knowledge of the physics of interaction, particle production and transport is necessary to estimate the radiation damage to equipment used on spacecraft and the biological effects of space radiation. For long duration astronaut missions, both on the International Space Station and the planned manned missions to Moon and Mars, the shielding strategy must include a comprehensive knowledge of the secondary radiation environment. The distribution of absorbed dose and dose equivalent is a function of the type, energy and population of these secondary products. Galactic cosmic rays (GCR) comprised of protons and heavier nuclei have energies from a few MeV per nucleon to the ZeV region, with the spectra reaching flux maxima in the hundreds of MeV range. Therefore, the MeV - GeV region is most important for space radiation. Coincidentally, the pion production energy threshold is about 280 MeV. The question naturally arises as to how important these particles are with respect to space radiation problems. The space radiation transport code, HZETRN (High charge (Z) and Energy TRaNsport), currently used by NASA, performs neutron, proton and heavy ion transport explicitly, but it does not take into account the production and transport of mesons, photons and leptons. In this paper, we present results from the Monte Carlo code MCNPX (Monte Carlo N-Particle eXtended), showing the effect of leptons and mesons when they are produced and transported in a GCR environment.

Aghara, S.K.↗

Proton Upset Monte Carlo Simulation

The Proton Upset Monte Carlo Simulation (PROPSET) program calculates the frequency of on-orbit upsets in computer chips (for given orbits such as Low Earth Orbit, Lunar Orbit, and the like) from proton bombardment based on the results of heavy ion testing alone. The software simulates the bombardment of modern microelectronic components (computer chips) with high-energy (.200 MeV) protons. The nuclear interaction of the proton with the silicon of the chip is modeled and nuclear fragments from this interaction are tracked using Monte Carlo techniques to produce statistically accurate predictions.

O'Neill, Patrick M.↗

Geology of Libya Montes and the Interbasin Plains of Northern Tyrrhena Terra, Mars: First Year Results and Second Year Work Plan

The Libya Montes-Tyrrhena Terra highland-lowland transitional zone of Mars is a complex tectonic and erosional region that contains some of the oldest exposed materials on the Martian surface as well as aqueous mineral signatures that may be potential chemical artifacts of early highland formational processes. Our 1:1M scale mapping project includes the geologic materials and landforms contained within MTMs 00282, -05282, -10282, 00277, - 05277, and -10277, which cover the highland portion of the transitional zone. The map region extends from the Libya Montes southward into Tyrrhena Terra and to the northern rim of Hellas basin and includes volcanic rocks of Syrtis Major Planum and a broad lowlying plain (palus) that forms a topographic divide between Isidis and Hellas basins. The objective of this project is to describe the geologic history of regional massif and plains materials by combining geomorphological and compositional mapping observations. This abstract summarizes the technical approaches and interim scientific results of Year 1 efforts and the expected work plan for Year 2 efforts.

Skinner, J. A., Jr.↗