Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo codes”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2012

Mapping the Milky Way Galaxy with LISA

Gravitational wave detectors in the mHz band (such as the Laser Interferometer Space Antenna, or LISA) will observe thousands of compact binaries in the galaxy which can be used to better understand the structure of the Milky Way. To test the effectiveness of LISA to measure the distribution of the galaxy, we simulated the Close White Dwarf Binary (CWDB) gravitational wave sky using different models for the Milky Way. To do so, we have developed a galaxy density distribution modeling code based on the Markov Chain Monte Carlo method. The code uses different distributions to construct realizations of the galaxy. We then use the Fisher Information Matrix to estimate the variance and covariance of the recovered parameters for each detected CWDB. This is the first step toward characterizing the capabilities of space-based gravitational wave detectors to constrain models for galactic structure, such as the size and orientation of the bar in the center of the Milky Way

McKinnon, Jose A.↗
NASA NTRS2014

Geant4 Predictions of Energy Spectra in Typical Space Radiation Environment

Accurate knowledge of energy spectra inside spacecraft is important for protecting astronauts as well as sensitive electronics from the harmful effects of space radiation. Such knowledge allows one to confidently map the radiation environment inside the vehicle. The purpose of this talk is to present preliminary calculations for energy spectra inside a spherical shell shielding and behind a slab in typical space radiation environment using the 3D Monte-Carlo transport code Geant4. We have simulated proton and iron isotropic sources and beams impinging on Aluminum and Gallium arsenide (GaAs) targets at energies of 0.2, 0.6, 1, and 10 GeV/u. If time permits, other radiation sources and beams (_, C, O) and targets (C, Si, Ge, water) will be presented. The results are compared to ground-based measurements where available.

Sabra, M. S.↗
NASA NTRS2015

Characterization of Vacuum Facility Background Gas Through Simulation and Considerations for Electric Propulsion Ground Testing

The background gas in a vacuum facility for electric propulsion ground testing is examined in detail through a series of cold flow simulations using a direct simulation Monte Carlo (DSMC) code. The focus here is on the background gas itself, its structure and characteristics, rather than assessing its interaction and impact on thruster operation. The background gas, which is often incorrectly characterized as uniform, is found to have a notable velocity within a test facility. The gas velocity has an impact on the proper measurement of pressure and the calculation of ingestion flux to a thruster. There are also considerations for best practices for tests that involve the introduction of supplemental gas flows to artificially increase the background pressure. All of these effects need to be accounted for to properly characterize the operation of electric propulsion thrusters across different ground test vacuum facilities.

Test Facilities↗
NASA NTRS2017

Scale Dependence of Cirrus Horizontal Heterogeneity Effects on TOA Measurements: MODIS Brightness Temperatures in the Thermal Infrared - Part I

This paper presents a study on the impact of cirrus cloud heterogeneities on MODIS simulated thermal infrared (TIR) brightness temperatures (BTs) at the top of the atmosphere (TOA) as a function of spatial resolution from 50 meters to 10 kilometers. A realistic 3-D (three-dimensional) cirrus field is generated by the 3DCLOUD model (average optical thickness of 1.4, cloudtop and base altitudes at 10 and 12 kilometers, respectively, consisting of aggregate column crystals of D (sub eff) equals 20 microns), and 3-D thermal infrared radiative transfer (RT) is simulated with the 3DMCPOL (3-D Monte Carlo Polarized) code. According to previous studies, differences between 3-D BT computed from a heterogenous pixel and 1-D (one-dimensional) RT computed from a homogeneous pixel are considered dependent at nadir on two effects: (i) the optical thickness horizontal heterogeneity leading to the plane-parallel homogeneous bias (PPHB); and the (ii) horizontal radiative transport (HRT) leading to the independent pixel approximation error (IPAE). A single but realistic cirrus case is simulated and, as expected, the PPHB mainly impacts the low-spatial resolution results (above approximately 250 meters), with averaged values of up to 5-7 K (thousand), while the IPAE mainly impacts the high-spatial resolution results (below approximately 250 meters) with average values of up to 1-2 K (thousand). A sensitivity study has been performed in order to extend these results to various cirrus optical thicknesses and heterogeneities by sampling the cirrus in several ranges of parameters. For four optical thickness classes and four optical heterogeneity classes, we have found that, for nadir observations, the spatial resolution at which the combination of PPHB and HRT effects is the smallest, falls between 100 and 250 meters. These spatial resolutions thus appear to be the best choice to retrieve cirrus optical properties with the smallest cloud heterogeneity-related total bias in the thermal infrared. For off-nadir observations, the average total effect is increased and the minimum is shifted to coarser spatial resolutions.

TOA↗
NASA NTRS2017

Investigation of the High-Energy Oxidation of FiberForm from DSMC Analysis of Molecular Beam Experiments

A collaborative effort between the University of Illinois at Urbana-Champaign (UIUC), NASA Ames Research Center (ARC) and Montana State University (MSU) succeeded at developing a new finite-rate carbon oxidation model from molecular beam scattering experiments on vitreous carbon (VC). We now aim to use the direct simulation Monte Carlo (DSMC) code SPARTA to apply the model to each fiber of the porous fibrous Thermal Protection Systems (TPS) material FiberForm (FF). The detailed micro-structure of FF was obtained from X-ray micro-tomography and then used in DSMC. Both experiments and simulations show that the CO/O products ratio increased at all temperatures from VC to FF. We postulate this is due to the larger number of collisions an O atom encounters inside the porous FF material compared to the flat surface of VC. For the simulations, we particularly focused on the lowest and highest temperatures studied experimentally, 1023 K and 1823 K, and found good agreement between the finite-rate DSMC simulations and experiments.

Oxidation↗
NASA NTRS2017

Development of a Detailed Surface Chemistry Framework in DSMC

Many of the current direct simulation Monte Carlo (DSMC) codes still employ only simple surface catalysis models. These include only basic mechanisms such as dissociation, recombination, and exchange reactions, without any provision for adsorption and finite rate kinetics. Incorporating finite rate chemistry at the surface is increasingly becoming a necessity for various applications such as high speed re-entry flows over thermal protection systems (TPS), micro-electro-mechanical systems (MEMS), surface catalysis, etc. In the recent years, relatively few works have examined finite-rate surface reaction modeling using the DSMC method.In this work, a generalized finite-rate surface chemistry framework incorporating a comprehensive list of reaction mechanisms is developed and implemented into the DSMC solver SPARTA. The various mechanisms include adsorption, desorption, Langmuir-Hinshelwood (LH), Eley-Rideal (ER), Collision Induced (CI), condensation, sublimation, etc. The approach is to stochastically model the various competing reactions occurring on a set of active sites. Both gas-surface (e.g., ER, CI) and pure-surface (e.g., LH, desorption) reaction mechanisms are incorporated. The reaction mechanisms could also be catalytic or surface altering based on the participation of the bulk-phase species (e.g., bulk carbon atoms). Marschall and MacLean developed a general formulation in which multiple phases and surface sites are used and we adopt a similar convention in the current work. Microscopic parameters of reaction probabilities (for gas-surface reactions) and frequencies (for pure-surface reactions) that are require for DSMC are computed from the surface properties and macroscopic parameters such as rate constants, sticking coefficients, etc. The energy and angular distributions of the products are decided based on the reaction type and input parameters. Thus, the user has the capability to model various surface reactions via user-specified reaction rate constants, surface properties and parameters.

Gas↗
NASA NTRS2018

Significance of DSMC Computed Aerothermal Environments in the Rarefied Regime for Atmospheric Entry Material Response

During Mars atmospheric entry, the Mars Science Laboratory (MSL) was protected by a 4.5 meters diameter ablative heatshield assembled in 113 tiles. The heatshield was made of NASA's flagship ablative material, the Phenolic Impregnated Carbon Ablator (PICA). Prior work compared the traditional one-dimensional and three-dimensional material response models at different locations in the heatshield. It was observed that the flow was basically one-dimensional in the nose and flank regions, but three-dimensional flow effects were observed in the outer flank. The objective of this work is to study the effects of the aerothermal environment on the material response. We extend prior work by computing aerothermal environments using the direct simulation Monte Carlo (DSMC) code SPARTA and the CFD code Data Parallel Line Relaxation (DPLR). SPARTA is used to compute environment in the rarefied regime prior to 48.4s of entry where the Knudsen number is such that the Navier-Stokes equations can be inaccurate. Similarly to previous work, the DPLR software is used to compute the hypersonic environment for laminar then turbulent boundary layer assumptions from 48.4 s up to 100 s after Entry Interface (EI) along the MSL 08-TPS-02/01a trajectory. We observe that extending the aerothermal environments to times prior to 48.4 s modifies the thermal response of the heat shield at the surface and in-depth; however the effects on the recession are minimal. Additionally, using the assumption of a turbulent boundary layer versus a laminar one leads to higher surface and in-depth temperatures, larger recession, and a displacement of the peak heating and peak recession location.

Aerothermal↗