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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Models for dynamic surface temperatures during rarefied aeropass maneuvers

New models for calculating dynamic surface temperatures in a direct simulation Monte Carlo particle method during atmospheric entry are described. The models take into account the radiative emissivity, thermal conductivity, and heat capacity of the material. They are applied to thin panels in a highly rarefied flows and yield the transient thermal response of individual surface facets. Simulation results for radiation panels with internal conductivity are verified by comparing them with free molecule theory.

Haas, Brian L.↗
NASA NTRS2003

Monte Carlo Methods in Materials Science Based on FLUKA and ROOT

A comprehensive understanding of mitigation measures for space radiation protection necessarily involves the relevant fields of nuclear physics and particle transport modeling. One method of modeling the interaction of radiation traversing matter is Monte Carlo analysis, a subject that has been evolving since the very advent of nuclear reactors and particle accelerators in experimental physics. Countermeasures for radiation protection from neutrons near nuclear reactors, for example, were an early application and Monte Carlo methods were quickly adapted to this general field of investigation. The project discussed here is concerned with taking the latest tools and technology in Monte Carlo analysis and adapting them to space applications such as radiation shielding design for spacecraft, as well as investigating how next-generation Monte Carlos can complement the existing analytical methods currently used by NASA. We have chosen to employ the Monte Carlo program known as FLUKA (A legacy acronym based on the German for FLUctuating KAscade) used to simulate all of the particle transport, and the CERN developed graphical-interface object-oriented analysis software called ROOT. One aspect of space radiation analysis for which the Monte Carlo s are particularly suited is the study of secondary radiation produced as albedoes in the vicinity of the structural geometry involved. This broad goal of simulating space radiation transport through the relevant materials employing the FLUKA code necessarily requires the addition of the capability to simulate all heavy-ion interactions from 10 MeV/A up to the highest conceivable energies. For all energies above 3 GeV/A the Dual Parton Model (DPM) is currently used, although the possible improvement of the DPMJET event generator for energies 3-30 GeV/A is being considered. One of the major tasks still facing us is the provision for heavy ion interactions below 3 GeV/A. The ROOT interface is being developed in conjunction with the CERN ALICE (A Large Ion Collisions Experiment) software team through an adaptation of their existing AliROOT (ALICE Using ROOT) architecture. In order to check our progress against actual data, we have chosen to simulate the ATIC14 (Advanced Thin Ionization Calorimeter) cosmic-ray astrophysics balloon payload as well as neutron fluences in the Mir spacecraft. This paper contains a summary of status of this project, and a roadmap to its successful completion.

Pinsky, Lawrence↗
NASA NTRS2016

Simulation of Satellite, Airborne and Terrestrial LiDAR with DART (I):Waveform Simulation with Quasi-Monte Carlo Ray Tracing

Light Detection And Ranging (LiDAR) provides unique data on the 3-D structure of atmosphere constituents and the Earth's surface. Simulating LiDAR returns for different laser technologies and Earth scenes is fundamental for evaluating and interpreting signal and noise in LiDAR data. Different types of models are capable of simulating LiDAR waveforms of Earth surfaces. Semi-empirical and geometric models can be imprecise because they rely on simplified simulations of Earth surfaces and light interaction mechanisms. On the other hand, Monte Carlo ray tracing (MCRT) models are potentially accurate but require long computational time. Here, we present a new LiDAR waveform simulation tool that is based on the introduction of a quasi-Monte Carlo ray tracing approach in the Discrete Anisotropic Radiative Transfer (DART) model. Two new approaches, the so-called "box method" and "Ray Carlo method", are implemented to provide robust and accurate simulations of LiDAR waveforms for any landscape, atmosphere and LiDAR sensor configuration (view direction, footprint size, pulse characteristics, etc.). The box method accelerates the selection of the scattering direction of a photon in the presence of scatterers with non-invertible phase function. The Ray Carlo method brings traditional ray-tracking into MCRT simulation, which makes computational time independent of LiDAR field of view (FOV) and reception solid angle. Both methods are fast enough for simulating multi-pulse acquisition. Sensitivity studies with various landscapes and atmosphere constituents are presented, and the simulated LiDAR signals compare favorably with their associated reflectance images and Laser Vegetation Imaging Sensor (LVIS) waveforms. The LiDAR module is fully integrated into DART, enabling more detailed simulations of LiDAR sensitivity to specific scene elements (e.g., atmospheric aerosols, leaf area, branches, or topography) and sensor configuration for airborne or satellite LiDAR sensors.

Gastellu-Etchegorry, Jean-Philippe↗
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 NTRS2010

Study of Plume Impingement Effects in the Lunar Lander Environment

Plume impingement effects from the descent and ascent engine firings of the Lunar Lander were analyzed in support of the Lunar Architecture Team under the Constellation Program. The descent stage analysis was performed to obtain shear and pressure forces on the lunar surface as well as velocity and density profiles in the flow field in an effort to understand lunar soil erosion and ejected soil impact damage which was analyzed as part of a separate study. A CFD/DSMC decoupled methodology was used with the Bird continuum breakdown parameter to distinguish the continuum flow from the rarefied flow. The ascent stage analysis was performed to ascertain the forces and moments acting on the Lunar Lander Ascent Module due to the firing of the main engine on take-off. The Reacting and Multiphase Program (RAMP) method of characteristics (MOC) code was used to model the continuum region of the nozzle plume, and the Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) was used to model the impingement results in the rarefied region. The ascent module (AM) was analyzed for various pitch and yaw rotations and for various heights in relation to the descent module (DM). For the ascent stage analysis, the plume inflow boundary was located near the nozzle exit plane in a region where the flow number density was large enough to make the DSMC solution computationally expensive. Therefore, a scaling coefficient was used to make the DSMC solution more computationally manageable. An analysis of the effectiveness of this scaling technique was performed by investigating various scaling parameters for a single height and rotation of the AM. Because the inflow boundary was near the nozzle exit plane, another analysis was performed investigating three different inflow contours to determine the effects of the flow expansion around the nozzle lip on the final plume impingement results.

Marichalar, Jeremiah↗
NASA NTRS2011

Calculation of Dose Deposition in 3D Voxels by Heavy Ions and Simulation of gamma-H2AX Experiments

The biological response to high-LET radiation is different from low-LET radiation due to several factors, notably difference in energy deposition and formation of radiolytic species. Of particular importance in radiobiology is the formation of double-strand breaks (DSB), which can be detected by -H2AX foci experiments. These experiments has revealed important differences in the spatial distribution of DSB induced by low- and high-LET radiations [1,2]. To simulate -H2AX experiments, models based on amorphous track with radial dose are often combined with random walk chromosome models [3,4]. In this work, a new approach using the Monte-Carlo track structure code RITRACKS [5] and chromosome models have been used to simulate DSB formation. At first, RITRACKS have been used to simulate the irradiation of a cubic volume of 5 m by 1) 450 1H+ ions of 300 MeV (LET 0.3 keV/ m) and 2) by 1 56Fe26+ ion of 1 GeV/amu (LET 150 keV/ m). All energy deposition events are recorded to calculate dose in voxels of 20 m. The dose voxels are distributed randomly and scattered uniformly within the volume irradiated by low-LET radiation. Many differences are found in the spatial distribution of dose voxels for the 56Fe26+ ion. The track structure can be distinguished, and voxels with very high dose are found in the region corresponding to the track "core". These high-dose voxels are not found in the low-LET irradiation simulation and indicate clustered energy deposition, which may be responsible for complex DSB. In the second step, assuming that DSB will be found only in voxels where energy is deposited by the radiation, the intersection points between voxels with dose > 0 and simulated chromosomes were obtained. The spatial distribution of the intersection points is similar to -H2AX foci experiments. These preliminary results suggest that combining stochastic track structure and chromosome models could be a good approach to understand radiation-induced DSB and chromosome aberrations.

Plante, I.↗
NASA NTRS1993

Thermal analysis of combinatorial solid geometry models using SINDA

Algorithms have been developed using Monte Carlo techniques to determine the thermal network parameters necessary to perform a finite difference analysis on Combinatorial Solid Geometry (CSG) models. Orbital and laser fluxes as well as internal heat generation are modeled to facilitate satellite modeling. The results of the thermal calculations are used to model the infrared (IR) images of targets and assess target vulnerability. Sample analyses and validation are presented which demonstrate code products.

Gerencser, Diane↗
NASA NTRS1991

ASRM radiation and flowfield prediction status

Existing and proposed methods for the prediction of plume radiation are discussed in terms of their application to the NASA Advanced Solid Rocket Motor (ASRM) and Space Shuttle Main Engine (SSME) projects. Extrapolations of the Solid Rocket Motor (SRM) are discussed with respect to preliminary predictions of the primary and secondary radiation environments. The methodology for radiation and initial plume property predictions are set forth, including a new code for scattering media and independent secondary source models based on flight data. The Monte Carlo code employs a reverse-evaluation approach which traces rays back to their point of absorption in the plume. The SRM sea-level plume model is modified to account for the increased radiation in the ASRM plume due to the ASRM's propellant chemistry. The ASRM cycle-1 environment predictions are shown to identify a potential reason for the shutdown spike identified with pre-SRM staging.

Reardon, J. E.↗
NASA NTRS2010

MIRO Computational Model

A computational model calculates the excitation of water rotational levels and emission-line spectra in a cometary coma with applications for the Micro-wave Instrument for Rosetta Orbiter (MIRO). MIRO is a millimeter-submillimeter spectrometer that will be used to study the nature of cometary nuclei, the physical processes of outgassing, and the formation of the head region of a comet (coma). The computational model is a means to interpret the data measured by MIRO. The model is based on the accelerated Monte Carlo method, which performs a random angular, spatial, and frequency sampling of the radiation field to calculate the local average intensity of the field. With the model, the water rotational level populations in the cometary coma and the line profiles for the emission from the water molecules as a function of cometary parameters (such as outgassing rate, gas temperature, and gas and electron density) and observation parameters (such as distance to the comet and beam width) are calculated.

Broderick, Daniel↗
NASA NTRS1988

Thermodynamics of dense molecular hydrogen-helium mixtures at high pressure

Effective intermolecular pair potentials derived from liquid hydrogen and helium shock wave experiments are presently used in Monte Carlo simulations of mixtures of hydrogen and helium, at densities of up to 1.2 g/cu cm. The model interaction Helmholtz free energy derived accurately reproduces both the Monte Carlo calculation results and the experimental data obtained for densities of up to about 0.6 g/cu cm. An equation of state is derived from the free energy expression that could be useful in interior models of such Jovian planets as Saturn.

Marley, Mark S.↗
NASA NTRS2013

Galactic Cosmic Ray Event-Based Risk Model (GERM) Code

This software describes the transport and energy deposition of the passage of galactic cosmic rays in astronaut tissues during space travel, or heavy ion beams in patients in cancer therapy. Space radiation risk is a probability distribution, and time-dependent biological events must be accounted for physical description of space radiation transport in tissues and cells. A stochastic model can calculate the probability density directly without unverified assumptions about shape of probability density function. The prior art of transport codes calculates the average flux and dose of particles behind spacecraft and tissue shielding. Because of the signaling times for activation and relaxation in the cell and tissue, transport code must describe temporal and microspatial density of functions to correlate DNA and oxidative damage with non-targeted effects of signals, bystander, etc. These are absolutely ignored or impossible in the prior art. The GERM code provides scientists data interpretation of experiments; modeling of beam line, shielding of target samples, and sample holders; and estimation of basic physical and biological outputs of their experiments. For mono-energetic ion beams, basic physical and biological properties are calculated for a selected ion type, such as kinetic energy, mass, charge number, absorbed dose, or fluence. Evaluated quantities are linear energy transfer (LET), range (R), absorption and fragmentation cross-sections, and the probability of nuclear interactions after 1 or 5 cm of water equivalent material. In addition, a set of biophysical properties is evaluated, such as the Poisson distribution for a specified cellular area, cell survival curves, and DNA damage yields per cell. Also, the GERM code calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle in a selected material. The GERM code makes the numerical estimates of basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the first option, properties of monoenergetic beams are treated. In the second option, the transport of beams in different materials is treated. Similar biophysical properties as in the first option are evaluated for the primary ion and its secondary particles. Additional properties related to the nuclear fragmentation of the beam are evaluated. The GERM code is a computationally efficient Monte-Carlo heavy-ion-beam model. It includes accurate models of LET, range, residual energy, and straggling, and the quantum multiple scattering fragmentation (QMSGRG) nuclear database.

Cucinotta, Francis A.↗
NASA NTRS2001

Response of Silicon-Based Linear Energy Transfer Spectrometers: Implication for Radiation Risk Assessment in Space Flights

There is considerable interest in developing silicon-based telescopes because of their compactness and low power requirements. Three such telescopes have been flown on board the Space Shuttle to measure the linear energy transfer spectra of trapped, galactic cosmic ray, and solar energetic particles. Dosimeters based on single silicon detectors have also been flown on the Mir orbital station. A comparison of the absorbed dose and radiation quality factors calculated from these telescopes with that estimated from measurements made with a tissue equivalent proportional counter show differences which need to be fully understood if these telescopes are to be used for astronaut radiation risk assessments. Instrument performance is complicated by a variety of factors. A Monte Carlo-based technique was developed to model the behavior of both single element detectors in a proton beam, and the performance of a two-element, wide-angle telescope, in the trapped belt proton field inside the Space Shuttle. The technique is based on: (1) radiation transport intranuclear-evaporation model that takes into account the charge and angular distribution of target fragments, (2) Landau-Vavilov distribution of energy deposition allowing for electron escape, (3) true detector geometry of the telescope, (4) coincidence and discriminator settings, (5) spacecraft shielding geometry, and (6) the external space radiation environment, including albedo protons. The value of such detailed modeling and its implications in astronaut risk assessment is addressed. c2001 Elsevier Science B.V. All rights reserved.

manned↗
NASA NTRS2014

Variations in Modeled Dengue Transmission over Puerto Rico Using a Climate Driven Dynamic Model

Dengue fever is a mosquito‐borne viral disease reemerging throughout much of the tropical Americas. Dengue virus transmission is explicitly influenced by climate and the environment through its primary vector, Aedes aegypti. Temperature regulates Ae. aegypti development, survival, and replication rates as well as the incubation period of the virus within the mosquito. Precipitation provides water for many of the preferred breeding habitats of the mosquito, including buckets, old tires, and other places water can collect. Because of variations in topography, ocean influences and atmospheric processes, temperature and rainfall patterns vary across Puerto Rico and so do dengue virus transmission rates. Using NASA's TRMM (Tropical Rainfall Measuring Mission) satellite for precipitation input, ground‐based observations for temperature input, and laboratory confirmed dengue cases reported by the Centers for Disease Control and Prevention for parameter calibration, we modeled dengue transmission at the county level across Puerto Rico from 2010‐2013 using a dynamic dengue transmission model that includes interacting vector ecology and epidemiological components. Employing a Monte Carlo approach, we performed ensembles of several thousands of model simulations for each county in order to resolve the model uncertainty arising from using different combinations of parameter values that are not well known. The top 1% of model simulations that best reproduced the reported dengue case data were then analyzed to determine the most important parameters for dengue virus transmission in each county, as well as the relative influence of climate variability on transmission. These results can be used by public health workers to implement dengue control methods that are targeted for specific locations and climate conditions.

Morin, Cory↗
NASA NTRS2013

Modeling Analysis for NASA GRC Vacuum Facility 5 Upgrade

A model of the VF5 test facility at NASA Glenn Research Center was developed using the direct simulation Monte Carlo Hypersonic Aerothermodynamics Particle (HAP) code. The model results were compared to several cold flow and thruster hot fire cases. The main uncertainty in the model is the determination of the effective sticking coefficient -- which sets the pumping effectiveness of the cryopanels and oil diffusion pumps including baffle transmission. An effective sticking coefficient of 0.25 was found to provide generally good agreement with the experimental chamber pressure data. The model, which assumes a cold diffuse inflow, also fared satisfactorily in predicting the pressure distribution during thruster operation. The model was used to assess other chamber configurations to improve the local effective pumping speed near the thruster. A new configuration of the existing cryopumps is found to show more than 2x improvement over the current baseline configuration.

cryopumping↗
NASA NTRS1992

Collisional simulations of particles in a narrow planetary ring

A model is used to investigate how particle-particle collisions might modify some consequences of satellite perturbations relevant to short-term features of narrow planetary rings. A Monte Carlo-like simulation ring-model particle velocity alteration allows the simulation of collisions while individual particles are tracked. Periodic features visible in Voyager images are reproduced in a 2D numerical model which incorporates the collisional simulation as well as the perturbations of adjacent satellites. Collisions are noted to wash out periodic features within one collisional relaxation time.

Kolvoord, Robert A.↗
NASA NTRS1979

Analysis of SRB reentry acoustic environments

Space shuttle solid rocket booster reentry aeroacoustic environments were estimated. Particular emphasis was given to the aft skirt/exit plane region for the Mach number regime 0.6 = or greater than M infinity = or less than 3.5. The analysis is based on the evaluation of wind tunnel model results in conjunction with Monte Carlo simulation of trajectory parameters. The experimental approach is described as well as the evaluation process utilized. Predicted environments are presented in terms of one-third octave band spectra representing space averaged values for critical regions on the solid rocket booster.

Coffin, T.↗
NASA NTRS1982

Lidar clear atmosphere multiple scattering dependence on receiver range

The degree of multiple scattering encountered by nadir-directed lidar in clear atmosphere conditions was assessed by the Monte Carlo method based on a model of the vertical distribution of aerosol scattering in the atmosphere. The lowest 3 km of the atmosphere were regarded as containing 90% of the aerosol optical thickness, with the unity normalized value for the forward peak of the aerosol phase function being 4.38 at a 694.3 nm laser wavelength. Results were obtained for two lidar receiver heights for three receiver field of view (FOV) angular halfwidths. An increase in receiver height was determined to cause a significant increase in the amount of angular scattering of the lidar signal. A factor of 20 change in receiver height produced an order of magnitude change in the single to multiple scattering ratio.

Spinhirne, J. D.↗