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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Consequences of MMOD Penetrations on the International Space Station

Risk from micrometeoroid and orbital debris (MMOD) impacts on space vehicles is often quantified in terms of the probability of no penetration (PNP). However, for large spacecraft, especially those with multiple compartments, a penetration may have a number of possible outcomes. The extent of the damage (diameter of hole, crack length or penetration depth), the location of the damage relative to critical equipment or crew, crew response, and even the time of day of the penetration are among the many factors that can affect the outcome. For the International Space Station (ISS), a Monte-Carlo style software code called Manned Spacecraft Crew Survivability (MSCSurv) is used to predict the probability of several outcomes of an MMOD penetration—broadly classified as loss of crew (LOC), crew evacuation (Evac), loss of escape vehicle (LEV), and nominal end of mission (NEOM). By generating large numbers of MMOD impacts (typically in the billions) and tracking the consequences, MSCSurv allows for the inclusion of a large number of parameters and models as well as enabling the consideration of uncertainties in the models and parameters. MSCSurv builds upon the results from NASA’s Bumper software (which provides the probability of penetration and critical input data to MSCSurv) to allow analysts to estimate the probability of LOC, Evac, LEV, and NEOM. This paper provides an overview of the methodology used by NASA to quantify LOC, Evac, LEV, and NEOM with particular emphasis on describing in broad terms how MSCSurv works and its capabilities and most significant models.

Evans, Hilary J.↗

NASA NTRS2019

Predicting the X-Ray Spectra of Stellar-mass Black Holes from Simulations

We describe results from a new technique for the prediction of complete, self-consistent X-ray spectra from three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations of black hole accretion flows. Density and cooling rate data from a harm3d GRMHD simulation are post-processed by an improved version of the Monte Carlo radiation transport code pandurata (in the corona) and the Feautrier solver PTRANSX (in the disk), with xstar subroutines. The codes are run in a sequential, iterative fashion to achieve globally energy-conserving and self-consistent radiation fields, temperature maps, and photoionization equilibria. The output is the X-ray spectrum as seen by a distant observer, including features, such as the Fe Kα emission line and corresponding K-edge absorption trough, due to disk reprocessing of coronal power. For the example cases we consider—a non-rotating 10 M(sub ⊙) black hole with solar abundances, accreting at 0.01, 0.03, 0.1, or 0.3 Eddington—we find spectra resembling actual observations of stellar-mass black holes in the soft or steep power-law state: broad thermal peaks (at 1–3 keV), steep power laws extending to high energy (Γ = 2.7–4.5), and prominent, asymmetric Fe Kα emission lines with equivalent widths in the range 40–400 eV (larger EW at lower accretion rates). By starting with simulation data, we obviate the need for parameterized descriptions of the accretion flow geometry—no a priori specification of the corona's shape or flux, or the disk temperature or density, etc., is needed. Instead, we apply the relevant physical principles to simulation output using appropriate numerical techniques; this procedure allows us to calculate inclination-dependent spectra after choosing only a small number of physically meaningful parameters: black hole mass and spin, accretion rate, and elemental abundances.

Kinch, Brooks E.↗

NASA NTRS2018

Variability of Disk Emission in Pre-Main-Sequence and Related Stars. IV. Investigating the Structural Changes in the Inner Disk Region of MWC 480

We present five epochs of near-IR observations of the protoplanetary disk around MWC 480 (HD 31648) obtained with the SpeX spectrograph on NASA’s Infrared Telescope Facility between 2007 and 2013, inclusive. Using the measured line fluxes in the Pa β and Br γ lines, we found the mass accretion rates to be (1.26–2.30) x 10(exp −7) Solar Mass/yr and (1.4–2.01) x 10(exp −7) Solar Mass/yr, respectively, but which varied by more than 50% from epoch to epoch. The spectral energy distribution reveals a variability of about 30% between 1.5 and 10 μm during this same period of time. We investigated the variability using of the continuum emission of the disk in using the Monte- Carlo Radiative Transfer Code HOCHUNK3D. We find that varying the height of the inner rim successfully produces a change in the NIR flux but lowers the far-IR emission to levels below all measured fluxes. Because the star exhibits bipolar flows, we utilized a structure that simulates an inner disk wind to model the variability in the near-IR, without producing flux levels in the far-IR that are inconsistent with existing data. For this object, variable near-IR emission due to such an outflow is more consistent with the data than changing the scale height of the inner rim of the disk.

Fernandes, Rachel B.↗

NASA NTRS2019

Estimating Drag and Heating Coefficients for Hollow Reentry Objects in Transitional Flow Using DSMC

In NASA’s Object Reentry Survival Analysis Tool (ORSAT), aerodynamic drag and aerothermal heating coefficients are computed for each of the free-molecular, continuum, and transitional flow regimes using analytical and semi-analytical methods. These methods are typically limited to convex, blunt objects (such as spheres) and are applied to other objects such as boxes and cylinders using multiplicative “shape factors” to account for the different behavior. Previous literature has analyzed the aerodynamic and aerothermodynamic properties of flow around sharp-edged objects like boxes and cylinders in transitional flow, though only those objects with solid external boundaries. However, many reentry objects we have encountered in real spacecraft have been hollow (i.e., with the potential to allow flow through them). We present here preliminary results from analyses performed using the NASA Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) on hollow cylinders and boxes (with varying wall thickness-diameter ratio).

Marichalar, Jeremiah J.↗

NASA NTRS1998

Dust Morphology of Comet Hale-Bopp (C/1995 01). II. Introduction of a Working Model

A Monte Carlo image simulation code for dust features in comets is applied to comet Hale-Bopp in order to model the object's persistent porcupine-like appearance on high-resolution images taken between May 11 and November 2, 1996. In this scenario, the dust-emission profile is dominated by several brief flare-ups, or puffs in the production of dust from one of the sources.

comet↗

NASA NTRS2021

Variability of Disk Emission in Pre-main Sequence and Related Stars. V. OccultationEvents from the Innermost Disk Region of the Herbig Ae Star HD 163296

HD 163296 is a Herbig Ae star that underwent a dramatic∼0.8 magnitude drop in brightness in the V photometric band in 2001 and a brightening in the near-IR in 2002. Because the star possesses Herbig–Haro objects traveling in out flowing bipolar jets, it was suggested that the drop in brightness was due to a clump of dust entrained in a disk wind, blocking the line of sight toward the star. In order to quantify this hypothesis, we investigated the brightness drop at visible wavelengths and the brightening at near-IR wavelengths of HD 163296 using the Monte Carlo Radiative Transfer Code,HOCHUNK3D. We created three models to understand the events. Model 1 describes the quiescent state of the system. Model 2 describes the change in structure that led to the drop in brightness in 2001.Model 3 describes the structure needed to produce the observed 2002 brightening of the near-IR wavelengths. Models 2 and 3 utilize a combination of a disk wind and central bipolar flow. By introducing a filled bipolar cavity in Models 2 and 3, we were able to successfully simulate a jet-like structure for the star with a disk wind and created the drop and subsequent increase in brightness of the system. On the other hand, when the bipolar cavity is not filled, Model 1 replicates the quiescent state of the system

Monika Pikhartova↗

NASA NTRS2015

Use of Existing CAD Models for Radiation Shielding Analysis

The utility of a radiation exposure analysis depends not only on the accuracy of the underlying particle transport code, but also on the accuracy of the geometric representations of both the vehicle used as radiation shielding mass and the phantom representation of the human form. The current NASA/Space Radiation Analysis Group (SRAG) process to determine crew radiation exposure in a vehicle design incorporates both output from an analytic High Z and Energy Particle Transport (HZETRN) code and the properties (i.e., material thicknesses) of a previously processed drawing. This geometry pre-process can be time-consuming, and the results are less accurate than those determined using a Monte Carlo-based particle transport code. The current work aims to improve this process. Although several Monte Carlo programs (FLUKA, Geant4) are readily available, most use an internal geometry engine. The lack of an interface with the standard CAD formats used by the vehicle designers limits the ability of the user to communicate complex geometries. Translation of native CAD drawings into a format readable by these transport programs is time consuming and prone to error. The Direct Accelerated Geometry -United (DAGU) project is intended to provide an interface between the native vehicle or phantom CAD geometry and multiple particle transport codes to minimize problem setup, computing time and analysis error.

Lee, K. T.↗

NASA NTRS2015

Space Radiation Transport Code Development: 3DHZETRN

The space radiation transport code, HZETRN, has been used extensively for research, vehicle design optimization, risk analysis, and related applications. One of the simplifying features of the HZETRN transport formalism is the straight-ahead approximation, wherein all particles are assumed to travel along a common axis. This reduces the governing equation to one spatial dimension allowing enormous simplification and highly efficient computational procedures to be implemented. Despite the physical simplifications, the HZETRN code is widely used for space applications and has been found to agree well with fully 3D Monte Carlo simulations in many circumstances. Recent work has focused on the development of 3D transport corrections for neutrons and light ions (Z < 2) for which the straight-ahead approximation is known to be less accurate. Within the development of 3D corrections, well-defined convergence criteria have been considered, allowing approximation errors at each stage in model development to be quantified. The present level of development assumes the neutron cross sections have an isotropic component treated within N explicit angular directions and a forward component represented by the straight-ahead approximation. The N = 1 solution refers to the straight-ahead treatment, while N = 2 represents the bi-directional model in current use for engineering design. The figure below shows neutrons, protons, and alphas for various values of N at locations in an aluminum sphere exposed to a solar particle event (SPE) spectrum. The neutron fluence converges quickly in simple geometry with N > 14 directions. The improved code, 3DHZETRN, transports neutrons, light ions, and heavy ions under space-like boundary conditions through general geometry while maintaining a high degree of computational efficiency. A brief overview of the 3D transport formalism for neutrons and light ions is given, and extensive benchmarking results with the Monte Carlo codes Geant4, FLUKA, and PHITS are provided for a variety of boundary conditions and geometries. Improvements provided by the 3D corrections are made clear in the comparisons. Developments needed to connect 3DHZETRN to vehicle design and optimization studies will be discussed. Future theoretical development will relax the forward plus isotropic interaction assumption to more general angular dependence.

Wilson, John W.↗

NASA NTRS1972

Radiation transport codes for potential applications related to radiobiology and radiotherapy using protons, neutrons, and negatively charged pions

Several Monte Carlo radiation transport computer codes are used to predict quantities of interest in the fields of radiotherapy and radiobiology. The calculational methods are described and comparisions of calculated and experimental results are presented for dose distributions produced by protons, neutrons, and negatively charged pions. Comparisons of calculated and experimental cell survival probabilities are also presented.

Armstrong, T. W.↗

NASA NTRS2022

Parallel Grand-Canonical Monte Carlo (ParaGrandMC) User’s Manual Version 2.0

This manual describes the commands and command line options for the Parallel Grand Canonical Monte Carlo version 2.0 (ParaGrandMC.2.0) simulation code. This is a highly scalable parallel FORTRAN 2003 code for simulating the thermodynamic evolution of materials at the atomic level, and predicting their thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is specifically designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces. In addition to performing Monte Carlo (MC) simulations, the code can also perform Molecular Dynamics (MD) and Langevin Dynamics (LD) simulations, which can be combined and interchanged with MC for faster and more efficient system evolution. A detailed description of the MC part of the code is provided in the NASA ParaGrandMC report: NASA/CR–2016-219202; http://www.sti.nasa.gov.

High performance computing↗

NASA NTRS1997

Earth and Planetary Science Letters

Systematic measurements of the concentrations of cosmogen Ca-41 (half-life = 1.04 x 10(exp 5) yr) in the Apollo 15 long core 15001-15006 were performed by accelerator mass spectroscopy. Earlier measurements of cosmogenic Be-10, C-14, Al-26, Cl-36, and Mn-53 in the same core have provided confirmation and improvement of theoretical models for predicting production profiles of nuclides by cosmic ray induced spallation in the Moon and large meteorites. Unlike these nuclides, Ca-40 in the lunar surface is produced mainly by thermal neutron capture reactions on Ca-40. The maximum production of Ca-41, about 1 dpm/g Ca, was observed at a depth in the Moon of about 150 g/sq cm. For depths below about 300 g/sq cm, Ca-41 production falls off exponentially with an e-folding length of 175 g/sq cm. Neutron production in the Moon was modeled with the Los Alamos High Energy Transport Code System, and yields of nuclei produced by low-energy thermal and epithermal neutrons were calculated with the Monte Carlo N-Particle code. The new theoretical calculations using these codes are in good agreement with our measured Ca-41 concentrations as well as with Co-60 and direct neutron fluence measurements in the Moon.

K. Nishiizumi↗

NASA NTRS2023

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante↗

NASA NTRS2014

A Monte Carlo Analysis of the Thrust Imbalance for the RSRMV Booster During Both the Ignition Transient and Steady State Operation

This paper presents the results of statistical analyses performed to predict the thrust imbalance between two solid rocket motor boosters to be used on the Space Launch System (SLS) vehicle. Two legacy internal ballistics codes developed for the Space Shuttle program were coupled with a Monte Carlo analysis code to determine a thrust imbalance envelope for the SLS vehicle based on the performance of 1000 motor pairs. Thirty three variables which could impact the performance of the motors during the ignition transient and thirty eight variables which could impact the performance of the motors during steady state operation of the motor were identified and treated as statistical variables for the analyses. The effects of motor to motor variation as well as variations between motors of a single pair were included in the analyses. The statistical variations of the variables were defined based on data provided by NASA's Marshall Space Flight Center for the upgraded five segment booster and from the Space Shuttle booster when appropriate. The results obtained for the statistical envelope are compared with the design specification thrust imbalance limits for the SLS launch vehicle

Foster, Winfred A., Jr.↗

NASA NTRS1993

Transport methods and interactions for space radiations

This report presents a brief history leading to the involvement of the Langley Research Center of the National Aeronautics and Space Administration (NASA) in space-radiation physics and protection. Indeed, a relatively complete summary of technical capability as of the summer of 1990 is given. The Boltzmann equations for coupled ionic and neutronic fields are presented and inversion techniques for the Boltzmann operator are discussed. Errors generated by the straight ahead approximation are derived and are shown to be negligible for most problems of space-radiation protection. A decoupling of projectile propagation from the target fields greatly simplifies the Boltzmann equations and allows an analytic solution of the target fragment transport. Analytic and numerical methods of solving the projectile transport equations are discussed. The nuclear physics underlying the coefficients in the Boltzmann equation is discussed. A coupled-channel optical model is found as a consequence of the loose binding of nuclear matter and closure of the nuclear states in high-energy reactions. Transport solutions with the developed data base are used with laboratory experiments to validate both the transport code and the data base. Numerical benchmarks and comparison with Monte Carlo calculations are also used for code validation.

Wilson, John W.↗

NASA NTRS2014

A Monte Carlo Analysis of the Thrust Imbalance for the Space Launch System Booster During Both the Ignition Transient and Steady State Operation

Foster, Winfred A., Jr.↗

NASA NTRS2002

Monte Carlo transport model comparison with 1A GeV accelerated iron experiment: heavy-ion shielding evaluation of NASA space flight-crew foodstuff

Deep-space manned flight as a reality depends on a viable solution to the radiation problem. Both acute and chronic radiation health threats are known to exist, with solar particle events as an example of the former and galactic cosmic rays (GCR) of the latter. In this experiment Iron ions of 1A GeV are used to simulate GCR and to determine the secondary radiation field created as the GCR-like particles interact with a thick target. A NASA prepared food pantry locker was subjected to the iron beam and the secondary fluence recorded. A modified version of the Monte Carlo heavy ion transport code developed by Zeitlin at LBNL is compared with experimental fluence. The foodstuff is modeled as mixed nuts as defined by the 71st edition of the Chemical Rubber Company (CRC) Handbook of Physics and Chemistry. The results indicate a good agreement between the experimental data and the model. The agreement between model and experiment is determined using a linear fit to ordered pairs of data. The intercept is forced to zero. The slope fit is 0.825 and the R2 value is 0.429 over the resolved fluence region. The removal of an outlier, Z=14, gives values of 0.888 and 0.705 for slope and R2 respectively. c2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

NASA Discipline Radiation Health↗

NASA NTRS1976

A Monte-Carlo simulation of galactic cosmic ray effects in the lunar regolith

A Monte-Carlo soil mixing computer code was used to evaluate isotopic effects induced by galactic cosmic rays in samples on the top surface of the lunar regolith and in lunar core tubes. For surface samples, good agreement was found between theoretical predictions of Ne-21 concentrations, neutron fluences, and Xe-131/Xe-126 ratios, and the corresponding observations.

Langevin, Y.↗

NASA NTRS1984

Investigation of the external flow analysis for density measurements at high altitude

The results of analysis performed on the external flow around the shuttle orbiter nose regions at the Shuttle Upper Atmosphere Mass Spectrometer (SUMS) inlet orifice are presented. The purpose of the analysis is to quantitatively characterize the flow conditions to facilitate SUMS flight data reduction and subsequent determination of orbiter aerodynamic force coefficients in the hypersonic rarefied flow regime. Experimental determination of aerodynamic force coefficients requires accurate simultaneous measurement of forces (or acceleration) and dynamic pressure along with independent knowledge of density and velocity. The SUMS provides independent measurement of dynamic pressure; however, it does so indirectly and requires knowledge of the relationship between measured orifice conditions and the dynamic pressure which can only be determined on the basis of molecule or theory for a winged configuration. Monte Carlo direct simulation computer codes were developed for both the flow field solution at the orifice and for the internal orifice flow. These codes were used to study issues associated with geometric modeling of the orbiter nose geometry and the modeling of intermolecular collisions including rotational energy exchange and a preliminary analysis of vibrational excitation and dissociation effects. Data obtained from preliminary simulation runs are presented.

Bienkowski, G. K.↗