Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
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↗
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↗
DOE OSTI2012

MCNP6 Cosmic-Source Option

MCNP is a Monte Carlo radiation transport code that has been under development for over half a century. Over the last decade, the development team of a high-energy offshoot of MCNP, called MCNPX, has implemented several physics and algorithm improvements important for modeling galactic cosmic-ray (GCR) interactions with matter. In this presentation, we discuss the latest of these improvements, a new Cosmic-Source option, that has been implemented in MCNP6.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2012

Preliminary Simulations for Geometric Optimization of a High-Energy Delayed Gamma Spectrometer for Direct Assay of Pu in Spent Nuclear Fuel

High-energy, beta-delayed gamma-ray spectroscopy is under investigation as part of the Next Generation Safeguard Initiative effort to develop non-destructive assay instruments for plutonium mass quantification in spent nuclear fuel assemblies. Results obtained to date indicate that individual isotope-specific signatures contained in the delayed gamma-ray spectra can potentially be used to quantify the total fissile content and individual weight fractions of fissile and fertile nuclides present in spent fuel. Adequate assay precision for inventory analysis can be obtained using a neutron generator of sufficient strength and currently available detection technology. In an attempt to optimize the geometric configuration and material composition for a delayed gamma measurement on spent fuel, the current study applies MCNPX, a Monte Carlo radiation transport code, in order to obtain the best signal-to-noise ratio. Results are presented for optimizing the neutron spectrum tailoring material, geometries to maximize thermal or fast fissions from a given neutron source, and detector location to allow an acceptable delayed gamma-ray signal while achieving a reasonable detector lifetime while operating in a high-energy neutron field. This work is supported in part by the Next Generation Safeguards Initiative, Office of Nuclear Safeguards and Security, National Nuclear Security Administration.

Delayed Gamma Spectroscopy, active interrogation, ↗
DOE OSTI2013

OBJECT KINETIC MONTE CARLO SIMULATIONS OF MICROSTRUCTURE EVOLUTION

The objective is to report the development of the flexible object kinetic Monte Carlo (OKMC) simulation code KSOME (kinetic simulation of microstructure evolution) which can be used to simulate microstructure evolution of complex systems under irradiation. In this report we briefly describe the capabilities of KSOME and present preliminary results for short term annealing of single cascades in tungsten at various primary-knock-on atom (PKA) energies and temperatures.

Nandipati, Giridhar↗
DOE OSTI2017

Using Whisper-1.1 to Guide Improvements to Nuclear Data Evaluations

Along with the latest release of the general purpose Monte Carlo radiation transport code, MCNP6.2, came the first officially released version of Whisper, version 1.1, a statistical analysis code used for nuclear criticality safety (NCS) validation. The Whisper software package is primarily used to calculate the upper subcritical limit (USL) of a particular application case along with components of the margin of subcriticality (MOS). In the calculation of the MOS due to nuclear data uncertainties, Whisper uses the generalized linear-least-squares (GLLS) approach to obtain adjusted nuclear data uncertainties which are then propagated to estimate the nuclear data uncertainties in the k eff of the application. Due to the fact that Whisper contains the simulated and measured k eff values for over 1100 benchmarks, nuclear data sensitivity profiles for each benchmark along with the nuclear data covariance matrices, all of which are used in the GLLS method, it is possible to easily utilize these data available to help guide improvements to nuclear data evaluations. In the present work, the Whisper data and the GLLS method are used together to suggest an adjusted nuclear data cross section evaluation specifically for 233 U such that the agreement between the calculated and experimental k eff measurements of the relevant benchmarks are improved. With this adjusted nuclear data available, feedback to the theoretical and experimental nuclear data community is provided.

97 MATHEMATICS AND COMPUTING↗