Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo code”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

X-Ray Astronomy

Dr. S. N. Zhang has lead a seven member group (Dr. Yuxin Feng, Mr. XuejunSun, Mr. Yongzhong Chen, Mr. Jun Lin, Mr. Yangsen Yao, and Ms. Xiaoling Zhang). This group has carried out the following activities: continued data analysis from space astrophysical missions CGRO, RXTE, ASCA and Chandra. Significant scientific results have been produced as results of their work. They discovered the three-layered accretion disk structure around black holes in X-ray binaries; their paper on this discovery is to appear in the prestigious Science magazine. They have also developed a new method for energy spectral analysis of black hole X-ray binaries; four papers on this topics were presented at the most recent Atlanta AAS meeting. They have also carried Monte-Carlo simulations of X-ray detectors, in support to the hardware development efforts at Marshall Space Flight Center (MSFC). These computation-intensive simulations have been carried out entirely on the computers at UAH. They have also carried out extensive simulations for astrophysical applications, taking advantage of the Monte-Carlo simulation codes developed previously at MSFC and further improved at UAH for detector simulations. One refereed paper and one contribution to conference proceedings have been resulted from this effort.

Wu, S. T.↗
NASA NTRS2001

Orbits of Saturn's Inner Moons and Other Observations Connected within the 1995-1996 Saturnian Ring Plane Crossing

We reduced and analyzed radio images of Saturn taken in 1990-1995 at 0.35, 2.0, 3.6 and 6.1 cm wavelength. The latitudinal brightness distribution varies substantially over time. Detailed modeling of these data shows that NH4SH condenses at the thermochemical equilibrium temperature of 235.5 K only at Northern mid-latitudes. Over most of the planet, condensation does not occur until the gas cools to 190 +/- 5 K. Supersaturation may also cause the dark equatorial region seen in 1995 at 6.1 cm. Observations of the rings show that the West (dusk) ansa is brighter than the East (dawn) ansa by factors up to 2. The magnitude of the asymmetry increases with increasing wavelength and with decreasing distance to the planet. This East-West asymmetry may be due to multiple scattering in gravitational wakes. We improved the Monte Carlo scattering code used to model the microwave emission from Saturn's rings at multiple wavelengths. The basic code outlined in Dunn, Molnar, and Fix describes the rings as a uniform, infinite layer of ice spheres with a finite vertical extent. The particle sizes follow a power law distribution and scatter light according to a linear combination of Mie and isotropic phase functions. These parameters may be adjusted to reflect the physical properties in the scattering properties of the ring particles. This is used to construct a model which can be directly compared to Very Large Array data. The code has been extended by including anisotropic particle distributions, primarily in the form of the wake structures that are likely to exist in the enhancements. The vertical extent of the wakes equals that of the ring layer. Between the wakes are regions of diminished density so that the optical depth integrated over many wakes matches that of the Voyager occultation experiments. This allows us to model the brightness temperature as a function of azimuth and look for asymmetric scattering behavior in the rings. By changing the width and relative density of the wakes, we can alter the nature of the gross scattering properties of the rings to match what is seen in our data.

Lissauer, Jack J.↗
NASA NTRS2001

Chemical and Flowfield Modeling for Enhanced Analysis of Contamination Experiments

This paper describes the application of a new Direct Simulation Monte Carlo (DSMC) code, the Molecular Beam Simulator (MBS), which is designed to analyze laboratory scale molecular beam-surface (and crossed-beam) experiments. The MBS is primarily intended to model experiments associated with spacecraft contamination effects, but it can also be used to simulate a variety of surface chemistry and reactive flow measurements. The MBS code is fully three-dimensional, includes a wide-range of chemical processes, and can model one or multiple pulsed (non-steady) sources. As an example application of the MBS code, a fast, pulsed, oxygen atom-surface experiment which examines the chemistry behind erosion of graphite by oxygen atoms is analyzed. Unsteady DSMC simulations show that experimental observations of excited molecular states after the pulse has hit the surface are consistent with two distinct chemical mechanisms: a direct one where the excited molecules are formed on the surface, and a two-step mechanism where ground state molecules formed on the surface are collisionally excited after they leave the surface by trailing oxygen atoms in the pulse. Further DSMC calculations suggest experiments which can distinguish between these mechanisms.

Braunstein, Matthew↗
NASA NTRS2003

Monte Carlo Methodology Serves Up a Software Success

Widely used for the modeling of gas flows through the computation of the motion and collisions of representative molecules, the Direct Simulation Monte Carlo method has become the gold standard for producing research and engineering predictions in the field of rarefied gas dynamics. Direct Simulation Monte Carlo was first introduced in the early 1960s by Dr. Graeme Bird, a professor at the University of Sydney, Australia. It has since proved to be a valuable tool to the aerospace and defense industries in providing design and operational support data, as well as flight data analysis. In 2002, NASA brought to the forefront a software product that maintains the same basic physics formulation of Dr. Bird's method, but provides effective modeling of complex, three-dimensional, real vehicle simulations and parallel processing capabilities to handle additional computational requirements, especially in areas where computational fluid dynamics (CFD) is not applicable. NASA's Direct Simulation Monte Carlo Analysis Code (DAC) software package is now considered the Agency s premier high-fidelity simulation tool for predicting vehicle aerodynamics and aerothermodynamic environments in rarified, or low-density, gas flows.

Source record↗
NASA NTRS2003

Modeling the Extended Neutral Atmosphere and Plasma Environment near Saturn

In the three years of this study we have published five papers in refereed journals. We have first examined satellite sources and their contribution to the observed neutral cloud. Based on the total calculated satellite sources and the spatial distribution of sputtered H20, we concluded that they cannot produce observed OH cloud. The E-ring contribution has been also studied in detail.In order to produce observed OH cloud we suggested that the E-ring might be the dominant source in inner Saturnian magnetosphere. We proposed a possible resupply mechanism which is needed to keep both E-ring and OH cloud in the present state: collisions between E-ring grains and remains of a disrupted satellite near Enceladus Lagrangian points. In this scenario a large amount of material, ranging from a few micrometers to hundred of meters, which is optically invisible at present, is likely to exist there. The fourth paper compares the magnetosheaths of the outer planets. A surprising result is that the hot proton component comprises about 40% of the total density, much larger than predicted by shock theory. Gas dynamic models of the boundaries show that the magnetospheres of Jupiter and Saturn are flattened at the poles. The last paper was published in GRL and is the first based of the model of neutrals developed as a main goal of this grant and which is now operational.This Monte Carlo collision code self- consistently determines the neutral distribution. from the rings and satellites until they are lost by ionization, by collisions with rings, moons, or Saturn, or by escape from Saturn. Our model is unique in that it includes the effects of plasma chemistry and both plasma-neutral and neutral neutral collisions to determine the dynamical evolution of the water group neutrals in Saturn's magnetosphere. The lifetimes of the neutrals against loss to photoionization, charge exchange, electron dissociation and electron impact dissociation are based on the model given by (Richardson et al. 1998) and vary with position in the magnetosphere. The dominant neutral dissociation channels H20->OH+H, H20-> O+H2, and OH->O+H2 are considered.

Richardson, John D.↗
NASA NTRS2004

The Role of Beam Geometry in Population Statistics and Pulse Profiles of Radio and Gamma-ray Pulsars

We present results of a pulsar population synthesis study that incorporates a number of recent developments and some significant improvements over our previous study. We have included the results of the Parkes multi-beam pulsar survey in our select group of nine radio surveys, doubling our sample of radio pulsars. More realistic geometries for the radio and gamma-ray beams are included in our Monte Carlo computer code that simulates the characteristics of the Galactic population of radio and gamma-ray pulsars. We adopted with some modifications the radio beam geometry of Arzoumanian, Chernoff & Cordes (2002). For the gamma-ray beam, we have assumed the slot gap geometry described in the work of Muslimov & Harding (2003). To account for the shape of the distribution of radio pulsars in the P(dot) - P diagram, we continue to find that decay of the magnetic field on a timescale of 2.8 Myr is needed. With all nine surveys, our model predicts that EGRET should have seen 7 radio-quiet (below the sensitivity of these radio surveys) and 19 radio-loud gamma-ray pulsars. AGILE (nominal sensitivity map) is expected to detect 13 radio-quiet and 37 radio-loud gamma-ray pulsars, while GLAST, with greater sensitivity is expected to detect 276 radio-quiet and 344 radio-loud gamma-ray pulsars. When the Parkes multi-beam pulsar survey is excluded, the ratio of radio-loud to radio-quiet gamma-ray pulsars decreases, especially for GLAST. The decrease for EGRET is 45%, implying that some fraction of EGRET unidentified sources are radio-loud gamma-ray pulsars. In the radio geometry adopted, short period pulsars are core dominated. Unlike the EGRET gamma-ray pulsars, our model predicts that when two gamma-ray peaks appear in the pulse profile, a dominant radio core peak appears in between the gamma-ray peaks. Our findings suggest that further improvements are required in describing both the radio and gamma-ray geometries.

Gonthier, Peter L.↗
NASA NTRS2004

Radiolytic Model for Chemical Composition of Europa's Atmosphere and Surface

The overall objective of the present effort is to produce models for major and selected minor components of Europa s neutral atmosphere in 1-D versus altitude and in 2-D versus altitude and longitude or latitude. A 3-D model versus all three coordinates (alt, long, lat) will be studied but development on this is at present limited by computing facilities available to the investigation team. In this first year we have focused on 1-D modeling with Co-I Valery Shematovich s Direct Simulation Monte Carlo (DSMC) code for water group species (H2O, O2, O, OH) and on 2-D with Co-I Mau Wong's version of a similar code for O2, O, CO, CO2, and Na. Surface source rates of H2O and O2 from sputtering and radiolysis are used in the 1-D model, while observations for CO2 at the Europa surface and Na detected in a neutral cloud ejected from Europa are used, along with the O2 sputtering rate, to constrain source rates in the 2-D version. With these separate approaches we are investigating a range of processes important to eventual implementation of a comprehensive 3-D atmospheric model which could be used to understand present observations and develop science requirements for future observations, e.g. from Earth and in Europa orbit. Within the second year we expect to merge the full water group calculations into the 2-D version of the DSMC code which can then be extended to 3-D, pending availability of computing resources. Another important goal in the second year would be the inclusion of sulk and its more volatile oxides (SO, SO2).

Cooper, John F.↗
NASA NTRS2004

DSMC Simulations of Hypersonic Flows With Shock Interactions and Validation With Experiments

The capabilities of a relatively new direct simulation Monte Carlo (DSMC) code are examined for the problem of hypersonic laminar shock/shock and shock/boundary layer interactions, where boundary layer separation is an important feature of the flow. Flow about two model configurations is considered, where both configurations (a biconic and a hollow cylinder-flare) have recent published experimental measurements. The computations are made by using the DS2V code of Bird, a general two-dimensional/axisymmetric time accurate code that incorporates many of the advances in DSMC over the past decade. The current focus is on flows produced in ground-based facilities at Mach 12 and 16 test conditions with nitrogen as the test gas and the test models at zero incidence. Results presented highlight the sensitivity of the calculations to grid resolution, sensitivity to physical modeling parameters, and comparison with experimental measurements. Information is provided concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Moss, James N.↗
NASA NTRS2004

DSMC Simulations of Hypersonic Flows With Shock Interactions and Validation With Experiments

The capabilities of a relatively new direct simulation Monte Carlo (DSMC) code are examined for the problem of hypersonic laminar shock/shock and shock/boundary layer interactions, where boundary layer separation is an important feature of the flow. Flow about two model configurations is considered, where both configurations (a biconic and a hollow cylinder-flare) have recent published experimental measurements. The computations are made by using the DS2V code of Bird, a general two-dimensional/axisymmetric time accurate code that incorporates many of the advances in DSMC over the past decade. The current focus is on flows produced in ground-based facilities at Mach 12 and 16 test conditions with nitrogen as the test gas and the test models at zero incidence. Results presented highlight the sensitivity of the calculations to grid resolutions, sensitivity to physical modeling parameters, and comparison with experimental measurements. Information is provided concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Moss, James N.↗
NASA NTRS2005

The Application of FLUKA to Dosimetry and Radiation Therapy

Monte Carlo transport codes like FLUKA are useful for many purposes, and one of those is the simulation of the effects of radiation traversing the human body. In particular, radiation has been used in cancer therapy for a long time, and recently this has been extended to include heavy ion particle beams. The advent of this particular type of therapy has led to the need for increased capabilities in the transport codes used to simulate the detailed nature of the treatment doses to the Y ~ ~ O U S tissues that are encountered. This capability is also of interest to NASA because of the nature of the radiation environment in space.[l] While in space, the crew members bodies are continually being traversed by virtually all forms of radiation. In assessing the risk that this exposure causes, heavy ions are of primary importance. These arise both from the primary external space radiation itself, as well as fragments that result from interactions during the traversal of that radiation through any intervening material including intervening body tissue itself. Thus the capability to characterize the details of the radiation field accurately within a human body subjected to such external 'beams" is of critical importance.

Wilson, Thomas L.↗
NASA NTRS2006

Blunt Body Aerodynamics for Hypersonic Low Density Flows

Numerical simulations are performed for the Apollo capsule from the hypersonic rarefied to the continuum regimes. The focus is on flow conditions similar to those experienced by the Apollo 6 Command Module during the high altitude portion of its reentry. The present focus is to highlight some of the current activities that serve as a precursor for computational tool assessments that will be used to support the development of aerodynamic data bases for future capsule flight environments, particularly those for the Crew Exploration Vehicle (CEV). Results for aerodynamic forces and moments are presented that demonstrate their sensitivity to rarefaction; that is, free molecular to continuum conditions. Also, aerodynamic data are presented that shows their sensitivity to a range of reentry velocities, encompassing conditions that include reentry from low Earth orbit, lunar return, and Mars return velocities (7.7 to 15 km/s). The rarefied results obtained with direct simulation Monte Carlo (DSMC) codes are anchored in the continuum regime with data from Navier-Stokes simulations.

Moss, James N.↗
NASA NTRS1996

3-D Simulations of NSTAR Ion Thruster Plasma Interactions

Described is a Particle-in-Cell with Monte Carlo Collision code developed to perform detailed three-dimensional ion thruster simulations. To capture the full kinetic behavior of ion thruster plumes, both the electrons and ions are treated as test particles. Simulation results are given of the NSTAR ion thruster under ground test and in space conditions. Numerical results are compared.

NSTAR Ion Thruster Plasma Ionization Particles↗
NASA NTRS2006

Orion Aerodynamics for Hypersonic Free Molecular to Continuum Conditions

Numerical simulations are performed for the Orion Crew Module, previously known as the Crew Exploration Vehicle (CEV) Command Module, to characterize its aerodynamics during the high altitude portion of its reentry into the Earth's atmosphere, that is, from free molecular to continuum hypersonic conditions. The focus is on flow conditions similar to those that the Orion Crew Module would experience during a return from the International Space Station. The bulk of the calculations are performed with two direct simulation Monte Carlo (DSMC) codes, and these data are anchored with results from both free molecular and Navier-Stokes calculations. Results for aerodynamic forces and moments are presented that demonstrate their sensitivity to rarefaction, that is, for free molecular to continuum conditions (Knudsen numbers of 111 to 0.0003). Also included are aerodynamic data as a function of angle of attack for different levels of rarefaction and results that demonstrate the aerodynamic sensitivity of the Orion CM to a range of reentry velocities (7.6 to 15 km/s).

Moss, James N.↗
NASA NTRS2009

Simulation of TGF-Beta Activation by Low-Dose HZE Radiation in a Cell Culture

High charge (Z) and energy (E) (HZE) nuclei comprised in the galactic cosmic rays are main contributors to space radiation risk. They induce many lesions in living matter such as non-specific oxidative damage and the double-strand breaks (DSBs), which are considered key precursors of early and late effects of radiation. There is increasing evidence that cells respond collectively rather than individually to radiation, suggesting the importance of cell signaling1. The transforming growth factor (TGF ) is a signaling peptide that is expressed in nearly all cell type and regulates a large array of cellular processes2. TGF have been shown to mediate cellular response to DNA damage3 and to induce apoptosis in non-irradiated cells cocultured with irradiated cells4. TFG molecules are secreted by cells in an inactive complex known as the latency-associated peptide (LAP). TGF is released from the LAP by a conformational change triggered by proteases, thrombospondin-1, integrins, acidic conditions and .OH radical5. TGF then binds to cells receptors and activates a cascade of events mediated by Smad proteins6, which might interfere with the repair of DNA. Meanwhile, increasingly sophisticated Brownian Dynamics (BD) algorithms have appeared recently in the literature7 and can be applied to study the interaction of molecules with receptors. These BD computer models have contributed to the elucidation of signal transduction, ligand accumulation and autocrine loops in the epidermal growth factor (EGF) and its receptor (EFGR) system8. To investigate the possible roles of TGF in an irradiated cell culture, our Monte-Carlo simulation codes of the radiation track structure9 will be used to calculate the activation of TFG triggered by .OH produced by low doses of HZE ions. The TGF molecules will then be followed by a BD algorithm in a medium representative of a cell culture to estimate the number of activated receptors.

Plante, Ianik↗
NASA NTRS2010

Monte-Carlo Simulation of Radiation Track Structure and Calculation of Dose Deposition in Nanovolumes

INTRODUCTION: The radiation track structure is of crucial importance to understand radiation damage to molecules and subsequent biological effects. Of a particular importance in radiobiology is the induction of double-strand breaks (DSBs) by ionizing radiation, which are caused by clusters of lesions in DNA, and oxidative damage to cellular constituents leading to aberrant signaling cascades. DSB can be visualized within cell nuclei with gamma-H2AX experiments. MATERIAL AND METHODS: In DSB induction models, the DSB probability is usually calculated by the local dose obtained from a radial dose profile of HZE tracks. In this work, the local dose imparted by HZE ions is calculated directly from the 3D Monte-Carlo simulation code RITRACKS. A cubic volume of 5 micron edge (Figure 1) is irradiated by a (Fe26+)-56 ion of 1 GeV/amu (LET approx.150 keV/micron) and by a fluence of 450 H+ ions, 300 MeV/amu (LET approx. 0.3 keV/micron). In both cases, the dose deposited in the volume is approx.1 Gy. The dose is then calculated into each 3D pixels (voxels) of 20 nm edge and visualized in 3D. RESULTS AND DISCUSSION: The dose is deposited uniformly in the volume by the H+ ions. The voxels which receive a high dose (orange) corresponds to electron track ends. The dose is deposited differently by the 56Fe26+ ion. Very high dose (red) is deposited in voxels with direct ion traversal. Voxels with electron track ends (orange) are also found distributed around the path of the track. In both cases, the appearance of the dose distribution looks very similar to DSBs seen in gammaH2AX experiments, particularly when the visualization threshold is applied. CONCLUSION: The refinement of the dose calculation to the nanometer scale has revealed important differences in the energy deposition between high- and low-LET ions. Voxels of very high dose are only found in the path of high-LET ions. Interestingly, experiments have shown that DSB induced by high-LET radiation are more difficult to repair. Therefore, this new approach may be useful to understand the nature of DSB and oxidative damage induced by ionizing radiation.

Plante, I.↗
NASA NTRS2011

New Modeling Approaches to Investigate Cell Signaling in Radiation Response

Ionizing radiation damages individual cells and tissues leading to harmful biological effects. Among many radiation-induced lesions, DNA double-strand breaks (DSB) are considered the key precursors of most early and late effects [1] leading to direct mutation or aberrant signal transduction processes. In response to damage, a flow of information is communicated to cells not directly hit by the radiation through signal transduction pathways [2]. Non-targeted effects (NTE), which includes bystander effects and genomic instability in the progeny of irradiated cells and tissues, may be particularly important for space radiation risk assessment [1], because astronauts are exposed to a low fluence of heavy ions and only a small fraction of cells are traversed by an ion. NTE may also have important consequences clinical radiotherapy [3]. In the recent years, new simulation tools and modeling approaches have become available to study the tissue response to radiation. The simulation of signal transduction pathways require many elements such as detailed track structure calculations, a tissue or cell culture model, knowledge of biochemical pathways and Brownian Dynamics (BD) propagators of the signaling molecules in their micro-environment. Recently, the Monte-Carlo simulation code of radiation track structure RITRACKS was used for micro and nano-dosimetry calculations [4]. RITRACKS will be used to calculate the fraction of cells traversed by an ion and delta-rays and the energy deposited in cells in a tissue model. RITRACKS also simulates the formation of chemical species by the radiolysis of water [5], notably the .OH radical. This molecule is implicated in DNA damage and in the activation of the transforming growth factor beta (TGF), a signaling molecule involved in NTE. BD algorithms for a particle near a membrane comprising receptors were also developed and will be used to simulate trajectories of signaling molecules in the micro-environment and characterize autocrine and paracrine cell communication and signal transduction.

Plante, Ianik↗
NASA NTRS2011

Ares I-X Ascent Base Environments

Plume induced base heating environments were measured during the flight of the NASA Constellation Ares I-X developmental launch vehicle, successfully flown on October 28, 2009. The Ares IX first stage is a four segment Space Shuttle derived booster with base consisting of a flared aft skirt, deceleration and tumble motors, and a thermal curtain surrounding the first stage 7.2 area ratio nozzle. Developmental Flight Instrumentation (DFI) consisted of radiometers, calorimeters, pressure transducers and gas temperature probes installed on the aft skirt and nozzle to measure the base environments. In addition, thermocouples were also installed between the layers of the flexible thermal curtain to provide insight into the curtain response to the base environments and to assist in understanding curtain failure during reentry. Plume radiation environment predictions were generated by the Reverse Monte Carlo (RMC) code and the convective base heating predictions utilized heritage MSFC empirical methods. These predictions were compared to the DFI data and results from the flight videography. Radiation predictions agreed with the flight measured data early in flight but gauge failures prevented high altitude comparisons. The convective environment comparisons demonstrated the need to improve the prediction methodology; particularly for low altitude, local plume recirculation. The convective comparisons showed relatively good agreement at altitudes greater than 50,000 feet.

Mobley, B. L.↗
NASA NTRS2012

Calculation of the Frequency Distribution of the Energy Deposition in DNA Volumes by Heavy Ions

Radiation quality effects are largely determined by energy deposition in small volumes of characteristic sizes less than 10 nm representative of short-segments of DNA, the DNA nucleosome, or molecules initiating oxidative stress in the nucleus, mitochondria, or extra-cellular matrix. On this scale, qualitatively distinct types of molecular damage are possible for high linear energy transfer (LET) radiation such as heavy ions compared to low LET radiation. Unique types of DNA lesions or oxidative damages are the likely outcome of the energy deposition. The frequency distribution for energy imparted to 1-20 nm targets per unit dose or particle fluence is a useful descriptor and can be evaluated as a function of impact parameter from an ions track. In this work, the simulation of 1-Gy irradiation of a cubic volume of 5 micron by: 1) 450 (1)H(+) ions, 300 MeV; 2) 10 (12)C(6+) ions, 290 MeV/amu and 3) (56)Fe(26+) ions, 1000 MeV/amu was done with the Monte-Carlo simulation code RITRACKS. Cylindrical targets are generated in the irradiated volume, with random orientation. The frequency distribution curves of the energy deposited in the targets is obtained. For small targets (i.e. <25 nm size), the probability of an ion to hit a target is very small; therefore a large number of tracks and targets as well as a large number of histories are necessary to obtain statistically significant results. This simulation is very time-consuming and is difficult to perform by using the original version of RITRACKS. Consequently, the code RITRACKS was adapted to use multiple CPU on a workstation or on a computer cluster. To validate the simulation results, similar calculations were performed using targets with fixed position and orientation, for which experimental data are available [5]. Since the probability of single- and double-strand breaks in DNA as function of energy deposited is well know, the results that were obtained can be used to estimate the yield of DSB, and can be extended to include other targeted or non-target effects.

Plante, Ianik↗