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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Monte Carlo Calculations of Polarized Microwave Radiation Emerging from Cloud Structures

The last decade has seen tremendous growth in cloud dynamical and microphysical models that are able to simulate storms and storm systems with very high spatial resolution, typically of the order of a few kilometers. The fairly realistic distributions of cloud and hydrometeor properties that these models generate has in turn led to a renewed interest in the three-dimensional microwave radiative transfer modeling needed to understand the effect of cloud and rainfall inhomogeneities upon microwave observations. Monte Carlo methods, and particularly backwards Monte Carlo methods have shown themselves to be very desirable due to the quick convergence of the solutions. Unfortunately, backwards Monte Carlo methods are not well suited to treat polarized radiation. This study reviews the existing Monte Carlo methods and presents a new polarized Monte Carlo radiative transfer code. The code is based on a forward scheme but uses aliasing techniques to keep the computational requirements equivalent to the backwards solution. Radiative transfer computations have been performed using a microphysical-dynamical cloud model and the results are presented together with the algorithm description.

Kummerow, Christian↗
NASA NTRS1998

Monte Carlo Simulations of Background Spectra in Integral Imager Detectors

Predictions of the expected gamma-ray backgrounds in the ISGRI (CdTe) and PiCsIT (Csl) detectors on INTEGRAL due to cosmic-ray interactions and the diffuse gamma-ray background have been made using a coupled set of Monte Carlo radiation transport codes (HETC, FLUKA, EGS4, and MORSE) and a detailed, 3-D mass model of the spacecraft and detector assemblies. The simulations include both the prompt background component from induced hadronic and electromagnetic cascades and the delayed component due to emissions from induced radioactivity. Background spectra have been obtained with and without the use of active (BGO) shielding and charged particle rejection to evaluate the effectiveness of anticoincidence counting on background rejection.

Armstrong, T. W.↗
NASA NTRS2000

Testing Models of Resonant Compton Scattering in X-Ray Pulsars

Over the performance period covered by the grant, the principal investigator modified a Monte Carlo Compton scattering code to model the propagation of x-rays through the magnetosphere of accreting neutron stars. These modifications were made to enable the author to compare the observations of x-ray pulsars to theoretical models of the system. The original code was designed to study relativistic plasmas with one of two geometries: a plane parallel plasma with a differential relativistic bulk velocity, and a static spherically symmetric plasma.- This code did not treat gravitational bending or bulk motion in the magnetosphere of a neutron star. Under the grant, the author incorporated code to trace light paths in a Schwarzschild metric. The code was modified to keep track of the photon polarization during propagati on. The investigator also modified the code so that bulk motion in an axisymmetric system is treated properly. An approximate treatment for resonant Compton scattering was added to the code. Finally, code was added that creates model observables that can be compared to observations, such as projected x-ray emission maps and energy-dependent light curves. Comparison to observations is now commencing.

Brainerd, Jerome J.↗
NASA NTRS2002

Bounding the Risk of Crew Loss Following Orbital Debris Penetration of the International Space Station at Assembly Stages 1J and 1E

Orbital debris impacts on the International Space Station occur frequently. To date, none of the impacting particles has been sufficiently large to penetrate manned pressurized volumes. We used the Manned Spacecraft Crew Survivability code to evaluate the risk to crew of penetrations of pressurized modules at two assembly stages: after Flight lJ, when the pressurized elements of Kibo, the Japanese Experiment Module, are present, and after Flight lE, when the European Columbus Module is present. Our code is a Monte Carlo simulation of impacts on the Station that considers several potential event types that could lead to crew loss. Among the statistics tabulated by the program is the probability of death of one or more crew members, expressed as the risk factor, R. This risk factor is dependent on details of crew operations during both ordinary circumstances and decompression emergencies, as well as on details of internal module configurations. We conducted trade studies considering these procedure and configuration details to determine the bounds on R at the 1J and 1E stages in the assembly sequence. Here we compare the R-factor bounds, and procedures and configurations that reduce R at these stages.

Evans, S.↗
NASA NTRS2003

Monte Carlo Simulations of the Response of the MARIE Instrument

The MARIE instrument aboard Mars Odyssey functions as a telescope for the detection of charged, energetic, nuclei. The directionality that leads to the telescope description is achieved by requiring coincident signals in two designated detectors in MARIE s silicon detector stack for the instrument to trigger. Because of this, MARIE is actually a bi directional telescope. Triggering particles can enter the detector stack by passing through the lightly shielded front of the instrument, but can also enter the back of the instrument by passing through the bulk of Odyssey. Because of this, understanding how to relate the signals recorded by MARIE to astrophysically important quantities such as particle fluxes or spectra exterior to the spacecraft clearly requires detailed modeling of the physical interactions that occur as the particles pass through the spacecraft and the instrument itself. In order to facilitate in the calibration of the MARIE data, we have begun a program to simulate the response of MARIE using the FLUKA [1] [2] Monte Carlo radiation transport code.

Andersen, V.↗
NASA NTRS2004

FUSE Observations of the Bright, Eclipsing Nova-like Cataclysmic Variable, UX UMa (FUSE 2000)

This was a project to study the disk and wind of the eclipsing nova-like variable UX UMa, in order to better define the wind geometry of the system, including the nature of the transition region between the disk photosphere and the supersonic wind. We proposed to use phase resolved spectroscopy of the system, taking advantage of the fact that UX UMa is an eclipsing system, to isolate different regions of the wind and to use a Monte Carlo radiative transfer code to simulate the spectra through the eclipse.

Long, Knox↗
NASA NTRS2000

A Radiation Shielding Code for Spacecraft and Its Validation

The HZETRN code, which uses a deterministic approach pioneered at NASA Langley Research Center, has been developed over the past decade to evaluate the local radiation fields within sensitive materials (electronic devices and human tissue) on spacecraft in the space environment. The code describes the interactions of shield materials with the incident galactic cosmic rays, trapped protons, or energetic protons from solar particle events in free space and low Earth orbit. The content of incident radiations is modified by atomic and nuclear reactions with the spacecraft and radiation shield materials. High-energy heavy ions are fragmented into less massive reaction products, and reaction products are produced by direct knockout of shield constituents or from de-excitation products. An overview of the computational procedures and database which describe these interactions is given. Validation of the code with recent Monte Carlo benchmarks, and laboratory and flight measurement is also included.

Shinn, J. L.↗
NASA NTRS1997

Modelling two dust features in comet Hale-Bopp (1995 O1)

A Monte Carlo computer simulation code for generating synthetic images of dust features in comets is applied to comet Hale-Bopp to model the diurnal evolution of a bright jet and the appearance of nearly concentric halos.

comet-Hale-Bopp diurnal evolution of a dust jet du↗
NASA NTRS2007

Quantifying Low Energy Proton Damage in Multijunction Solar Cells

An analysis of the effects of low energy proton irradiation on the electrical performance of triple junction (3J) InGaP2/GaAs/Ge solar cells is presented. The Monte Carlo ion transport code (SRIM) is used to simulate the damage profile induced in a 3J solar cell under the conditions of typical ground testing and that of the space environment. The results are used to present a quantitative analysis of the defect, and hence damage, distribution induced in the cell active region by the different radiation conditions. The modelling results show that, in the space environment, the solar cell will experience a uniform damage distribution through the active region of the cell. Through an application of the displacement damage dose analysis methodology, the implications of this result on mission performance predictions are investigated.

Messenger, Scott R.↗
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 NTRS2011

Probing the Accretion Geometry of Black Holes with X-Ray Polarization

In the coming years, new space missions will be able to measure X-ray polarization at levels of 1% or better in the approx.1-10 keV energy band. In particular, X-ray polarization is an ideal tool for determining the nature of black hole (BH) accretion disks surrounded by hot coronae. Using a Monte Carlo radiation transport code in full general relativity, we calculate the spectra and polarization features of these BH systems. At low energies, the signal is dominated by the thermal flux coming directly from the optically thick disk. At higher energies, the thermal seed photons have been inverse-Compton scattered by the corona, often reflecting back off the disk before reaching the observer, giving a distinctive polarization signature. By measuring the degree and angle of this X-ray polarization, we can infer the BH inclination, the emission geometry of the accretion flow, and also determine the spin of the black hole.

Schnitman, Jeremy D.↗
NASA NTRS2012

Performance Considerations for the SIMPL Single Photon, Polarimetric, Two-Color Laser Altimeter as Applied to Measurements of Forest Canopy Structure and Composition

The Slope Imaging Multi-polarization Photon-counting Lidar (SIMPL) is a multi-beam, micropulse airborne laser altimeter that acquires active and passive polarimetric optical remote sensing measurements at visible and near-infrared wavelengths. SIMPL was developed to demonstrate advanced measurement approaches of potential benefit for improved, more efficient spaceflight laser altimeter missions. SIMPL data have been acquired for wide diversity of forest types in the summers of 2010 and 2011 in order to assess the potential of its novel capabilities for characterization of vegetation structure and composition. On each of its four beams SIMPL provides highly-resolved measurements of forest canopy structure by detecting single-photons with 15 cm ranging precision using a narrow-beam system operating at a laser repetition rate of 11 kHz. Associated with that ranging data SIMPL provides eight amplitude parameters per beam unlike the single amplitude provided by typical laser altimeters. Those eight parameters are received energy that is parallel and perpendicular to that of the plane-polarized transmit pulse at 532 nm (green) and 1064 nm (near IR), for both the active laser backscatter retro-reflectance and the passive solar bi-directional reflectance. This poster presentation will cover the instrument architecture and highlight the performance of the SIMPL instrument with examples taken from measurements for several sites with distinct canopy structures and compositions. Specific performance areas such as probability of detection, after pulsing, and dead time, will be highlighted and addressed, along with examples of their impact on the measurements and how they limit the ability to accurately model and recover the canopy properties. To assess the sensitivity of SIMPL's measurements to canopy properties an instrument model has been implemented in the FLIGHT radiative transfer code, based on Monte Carlo simulation of photon transport. SIMPL data collected in 2010 over the Smithsonian Environmental Research Center, MD are currently being modelled and compared to other remote sensing and in situ data sets. Results on the adaptation of FLIGHT to model micropulse, single'photon ranging measurements are presented elsewhere at this conference. NASA's ICESat-2 spaceflight mission, scheduled for launch in 2016, will utilize a multi-beam, micropulse, single-photon ranging measurement approach (although non-polarimetric and only at 532 nm). Insights gained from the analysis and modelling of SIMPL data will help guide preparations for that mission, including development of calibration/validation plans and algorithms for the estimation of forest biophysical parameters.

Dabney, Philip W.↗
NASA NTRS2017

Estimating Consequences of MMOD Penetrations on ISS

The threat 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 briefly describes the overall 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, H.↗
NASA NTRS2016

Panchromatic Imaging of a Transitional Disk: The Disk of GM AUR in Optical and FUV Scattered Light

We have imaged GM Aurigae with the Hubble Space Telescope, detected its disk in scattered light at 1400 and 1650 A, and compared these with observations at 3300 A, 5550 A, 1.1 micron, and 1.6 micron. The scattered light increases at shorter wavelengths. The radial surface brightness profile at 3300 A shows no evidence of the 24 au radius cavity that has been previously observed in submillimeter observations. Comparison with dust grain opacity models indicates that the surface of the entire disk is populated with submicron grains. We have compiled a spectral energy distribution from 0.1 micron to 1 mm and used it to constrain a model of the star + disk system that includes the submillimeter cavity using the Monte Carlo radiative transfer code by Barbara Whitney. The best-fit model image indicates that the cavity should be detectable in the F330W bandpass if the cavity has been cleared of both large and small dust grains, but we do not detect it. The lack of an observed cavity can be explained by the presence of submicron grains interior to the submillimeter cavity wall. We suggest one explanation for this that could be due to a planet of mass less than 9M(sub J) interior to 24 au. A unique cylindrical structure is detected in the far- UV data from the Advanced Camera for Surveys/Solar Blind Channel. It is aligned along the system semiminor axis, but does not resemble an accretion-driven jet. The structure is limb brightened and extends 190 +/- 35 au above the disk midplane. The inner radius of the limb brightening is 40 +/- 10 au, just beyond the submillimeter cavity wall.

circumstellar matter aEuro" protoplanetary disks a↗
NASA NTRS2016

Constraining the Movement of the Spiral Features and the Locations of Planetary Bodies Within the AB Aur System

We present a new analysis of multi-epoch, H-band, scattered light images of the AB Aur system. We use a Monte Carlo radiative transfer code to simultaneously model the systems spectral energy distribution (SED) and H-band polarized intensity (PI) imagery. We find that a disk-dominated model, as opposed to one that is envelope dominated, can plausibly reproduce AB Aurs SED and near-IR imagery. This is consistent with previous modeling attempts presented in the literature and supports the idea that at least a subset of AB Aurs spirals originate within the disk. In light of this, we also analyzed the movement of spiral structures in multi-epoch H-band total light and PI imagery of the disk. We detect no significant rotation or change in spatial location of the spiral structures in these data, which span a 5.8-year baseline. If such structures are caused by disk planet interactions, the lack of observed rotation constrains the location of the orbit of planetary perturbers to be 47 au.

AB Aur system↗
NASA NTRS2018

Predicting the Consequences of MMOD Penetrations on the International Space Station

The threat 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 briefly describes the overall 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.

Hyde, James↗
NASA NTRS2017

The Shadow Knows: Using Shadows to Investigate the Structure of the Pretransitional Disk of HD 100453

We present Gemini Planet Imager polarized intensity imagery of HD 100453 in Y, J, and K1 bands that reveals an inner gap (9-18 au), an outer disk (18-39 au) with two prominent spiral arms, and two azimuthally localized dark features that are also present in Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) total intensity images. Spectral energy distribution fitting further suggests that the radial gap extends to 1 au. The narrow, wedge-like shape of the dark features appears similar to predictions of shadows cast by an inner disk that is misaligned with respect to the outer disk. Using the Monte Carlo radiative transfer code HOCHUNCK3D, we construct a model of the disk that allows us to determine its physical properties in more detail. From the angular separation of the features, we measure the difference in inclination between the disks (45deg) and their major axes, PA = 140deg east of north for the outer disk, and 100deg for the inner disk. We find an outer-disk inclination of 25deg +/- 10deg from face-on, in broad agreement with the Wagner et al. measurement of 34deg. SPHERE data in J and H bands indicate a reddish disk, which indicates that HD 100453 is evolving into a young debris disk.

Spectro-Polarimetric High-contrast Exoplanet REsea↗
NASA NTRS2018

Electromagnetic Chirps from Neutron Star-Black Hole Mergers

We calculate the electromagnetic signal of a gamma-ray flare coming from the surface of a neutron star shortly before merger with a black hole companion. Using a new version of the Monte Carlo radiation transport code Pandurata that incorporates dynamic spacetimes, we integrate photon geodesics from the neutron star surface until they reach a distant observer or are captured by the black hole. The gamma-ray light curve is modulated by a number of relativistic effects, including Doppler beaming and gravitational lensing. Because the photons originate from the inspiraling neutron star, the light curve closely resembles the corresponding gravitational waveform: a chirp signal characterized by a steadily increasing frequency and amplitude. We propose to search for these electromagnetic chirps using matched filtering algorithms similar to those used in LIGO data analysis.

gamma-ray burst↗