Results for “Monte Carlo Modelling” — Page 64

NASA NTRS2011

Calculating Launch Vehicle Flight Performance Reserve

This paper addresses different methods for determining the amount of extra propellant (flight performance reserve or FPR) that is necessary to reach orbit with a high probability of success. One approach involves assuming that the various influential parameters are independent and that the result behaves as a Gaussian. Alternatively, probabilistic models may be used to determine the vehicle and environmental models that will be available (estimated) for a launch day go/no go decision. High-fidelity closed-loop Monte Carlo simulation determines the amount of propellant used with each random combination of parameters that are still unknown at the time of launch. Using the results of the Monte Carlo simulation, several methods were used to calculate the FPR. The final chosen solution involves determining distributions for the pertinent outputs and running a separate Monte Carlo simulation to obtain a best estimate of the required FPR. This result differs from the result obtained using the other methods sufficiently that the higher fidelity is warranted.

Hanson, John M.↗

NASA NTRS2012

Pulsar Pair Cascades in Magnetic Fields with Offset Polar Caps

Neutron star magnetic fields may have polar caps (PC) that are offset from the dipole axis, through field-line sweepback near the light cylinder or non-symmetric currents within the star. The effects of such offsets on electron-positron pair cascades are investigated, using simple models of dipole magnetic fields with small distortions that shift the PCs by different amounts or directions. Using a Monte Carlo pair cascade simulation, we explore the changes in the pair spectrum, multiplicity and energy flux across the PC, as well as the trends in pair flux and pair energy flux with spin-down luminosity, L(sub sd). We also give an estimate of the distribution of heating flux from returning positrons on the PC for different offsets. We find that even modest offsets can produce significant increases in pair multiplicity, especially for pulsars that are near or beyond the pair death lines for centered PCs, primarily because of higher accelerating fields. Pair spectra cover several decades in energy, with the spectral range of millisecond pulsars (MSPs) two orders of magnitude higher than for normal pulsars, and PC offsets allow significant extension of all spectra to lower pair energies. We find that the total PC pair luminosity L(sub pair) is proportional to L(sub sd), with L(sub pair) approximates 10(exp -3) L(sub sd) for normal pulsars and L(sub pair) approximates 10(exp -2) L(sub sd) for MSPs. Remarkably, the total PC heating luminosity for even large offsets increases by less than a factor of two, even though the PC area increases by much larger factors, because most of the heating occurs near the magnetic axis.

Harding, Alice K.↗

NASA NTRS2017

A Space-Based Observational Strategy for Characterizing the First Stars and Galaxies Using the Redshifted 21cm Global Spectrum

The redshifted 21 cm monopole is expected to be a powerful probe of the epoch of the first stars and galaxies(10 less than z less than 35). The global 21 cm signal is sensitive to the thermal and ionization state of hydrogen gas and thusprovides a tracer of sources of energetic photonsprimarily hot stars and accreting black holeswhich ionize andheat the high redshift intergalactic medium (IGM). This paper presents a strategy for observations of the globalspectrum with a realizable instrument placed in a low-altitude lunar orbit, performing night-time 40120 MHzspectral observations, while on the farside to avoid terrestrial radio frequency interference, ionospheric corruption,and solar radio emissions. The frequency structure, uniformity over large scales, and unpolarized state of theredshifted 21 cm spectrum are distinct from the spectrally featureless, spatially varying, and polarized emissionfrom the bright foregrounds. This allows a clean separation between the primordial signal and foregrounds. Forsignal extraction, we model the foreground, instrument, and 21 cm spectrum with eigenmodes calculated viaSingular Value Decomposition analyses. Using a Markov Chain Monte Carlo algorithm to explore the parameterspace defined by the coefficients associated with these modes, we illustrate how the spectrum can be measured andhow astrophysical parameters (e.g., IGM properties, first star characteristics) can be constrained in the presence offoregrounds using the Dark Ages Radio Explorer (DARE).

Burns, Jack O.↗

NASA NTRS2018

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan‐1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017, https://doi.org/10.1126/sciadv.1701471). The model reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U(sub 0) = 0.5 eV, U(sub W) = 0.078–0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

Tucker, O. J.↗

NASA NTRS2018

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan-1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017). Themodel reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U(o) = 0.5eV, U(w) = 0.078 – 0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

O. J. Tucker↗

NASA NTRS2018

Solar Wind Implantation Into the Lunar Regolith: Monte Carlo Simulations of H Retention in a Surface With Defects and the H2 Exosphere

The solar wind implants protons into the top 20–30 nm of lunar regolith grains, and the implanted hydrogen will diffuse out of the regolith but also interact with oxygen in the regolith oxides. We apply a statistical approach to estimate the diffusion of hydrogen in the regolith hindered by forming temporary bonds with regolith oxygen atoms. A Monte Carlo simulation was used to track the temporal evolution of bound OH surface content and the H2 exosphere. The model results are consistent with the interpretation of the Chandrayaan‐1 M3 observations of infrared absorption spectra by surface hydroxyls as discussed in Li and Milliken (2017, https://doi.org/10.1126/sciadv.1701471). The model reproduced the latitudinal concentration of OH by using a Gaussian energy distribution of f(U0 = 0.5 eV, UW = 0.078–0.1 eV) to characterize the activation energy barrier to the diffusion of hydrogen in space weathered regolith. In addition, the model results of the exospheric content of H2 are consistent with observations by the Lyman Alpha Mapping Project on the Lunar Reconnaissance Orbiter. Therefore, we provide support for hydroxyl formation by chemically trapped solar wind protons.

O J Tucker↗

NASA NTRS2022

Numerical phase space optics methods and applications to the analysis of fiber coupling efficiency in atmospheric turbulence

A fundamental requirement of free-space optical communication is the ability to efficiently couple atmospherically distorted light from a telescope to a detector. A numerical method is presented for modeling fiber-based receiver performance in atmospheric conditions based on phase space optics which does not rely on Monte Carlo methods. This method is employed to analyze the waveguide insertion loss and optimal coupling geometry in atmospheric conditions for step-index fibers, graded-index fibers, and photonic lanterns with and without tilt compensation and central obscurations in the telescope.

free-space optical communication↗

NASA NTRSDate not supplied

Numerical phase space optics methods and applications to the analysis of fiber coupling efficiency in atmospheric turbulence

A fundamental requirement of free-space optical communication is the ability to efficiently couple atmospherically distorted light from a telescope to a detector. A numerical method is presented for modeling fiber-based receiver performance in atmospheric conditions based on phase space optics which does not rely on Monte Carlo methods. This method is employed to analyze the waveguide insertion loss and optimal coupling geometry in atmospheric conditions for step-index fibers, graded-index fibers, and photonic lanterns with and without tilt compensation and central obscurations in the telescope.

free-space optical communication↗

NASA NTRSDate not supplied

Simulation of etch pit formation through active sites in carbon fiber micro-structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation as a result of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gasses such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

NASA NTRSDate not supplied

Simulation of Etch Pit Formation in DSMC Through Active Sites in Carbon Fiber Micro-Structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation because of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gases such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various detailed surface reaction mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

NASA NTRSDate not supplied

DSMC Simulation of Etch Pit Formation Through Active Sites in Carbon Fiber Micro-Structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation because of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gases such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various detailed surface reaction mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

NASA NTRS1975

UV stars and the interstellar medium - A statistical time-dependent model

A statistical time-dependent model for the effects on the interstellar gas of the ionizing radiation from UV stars is presented. The radiation from the stars is assumed to be emitted as monochromatic bursts at random times and locations. The thermal and ionization history of the gas is followed in a Monte Carlo simulation. Probability densities as well as mean values for the ionization fraction, temperature, and other observable quantities are derived. The fluctuations in observed quantities due to the stochastic nature of the model are estimated. Models of the interstellar medium perpendicular to the galactic plane are developed and compared with observations. The action of the UV stars can explain the bulk properties of the interstellar gas toward the galactic poles and is consistent with the Copernicus satellite observations toward lambda Scorpii.

Lyon, J.↗

NASA NTRS2008

Radiation Transport of Heliospheric Lyman-alpha from Combined Cassini and Voyager Data Sets

Heliospheric neutral hydrogen scatters solar Lyman-alpha radiation from the Sun with '27-day' intensity modulations observed near Earth due to the Sun's rotation combined with Earth's orbital motion. These modulations are increasingly damped in amplitude at larger distances from the Sun due to multiple scattering in the heliosphere, providing a diagnostic of the interplanetary neutral hydrogen density independent of instrument calibration. This paper presents Cassini data from 2003-2004 obtained downwind near Saturn at approximately 10 AU that at times show undamped '27-day' waves in good agreement with the single-scattering models of Pryor et al., 1992. Simultaneous Voyager 1 data from 2003- 2004 obtained upwind at a distance of 88.8-92.6 AU from the Sun show waves damped by a factor of -0.21. The observed degree of damping is interpreted in terms of Monte Carlo multiple-scattering calculations (e.g., Keller et al., 1981) applied to two heliospheric hydrogen two-shock density distributions (discussed in Gangopadhyay et al., 2006) calculated in the frame of the Baranov-Malama model of the solar wind interaction with the two-component (neutral hydrogen and plasma) interstellar wind (Baranov and Malama 1993, Izmodenov et al., 2001, Baranov and Izmodenov, 2006). We conclude that multiple scattering is definitely occurring in the outer heliosphere. Both models compare favorably to the data, using heliospheric neutral H densities at the termination shock of 0.085 cm(exp -3) and 0.095 cm(exp -3). This work generally agrees with earlier discussions of Voyager data in Quemerais et al., 1996 showing the importance of multiple scattering but is based on Voyager data obtained at larger distances from the Sun (with larger damping) simultaneously with Cassini data obtained closer to the Sun.

Pryor, W.↗

NASA NTRSDate not supplied

Modeling Space-Radiation Induced Chromosome Aberrations

Exposure to ionizing radiation is identified as one of the main hazards of space flight. Galactic cosmic rays (GCRs) are composed of energetic protons (87%), helium nuclei (12%) and high-charge and energy (HZE) ions. At the cellular scale, HZE ions have a pattern of energy deposition that depends on the ion linear energy transfer (LET), with high LET ions favoring the formation of complex DNA damages that induce chromosome aberrations (CAs). Understanding how physical (e.g., ion type and energy) and biological (e.g., cell type) parameters impact CA formation is required to assess space radiation risks such as carcinogenesis. The Monte Carlo tool RITCARD/RITRACKS [1,2] developed by NASA allows simulation of the transport of HZE ions and subsequent energy deposition at the cellular scale and models the creation of DNA double strand breaks (DSBs), DNA repair and formation of CAs. We used this radiation transport code to further understand how radiation quality and cell size affects CA formation.

Floriane Poignant↗

NASA NTRS2018

High Fidelity Collision Probabilities Estimated Using Brute Force Monte Carlo Simulations

The NASA Conjunction Assessment Risk Analysis team has implemented new software to estimate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm employs a brute force Monte Carlo (BFMC) method which differs from most other methods because it uses orbital states and covariances propagated from their orbit determination epoch times using the full set of Astrodynamics Support Workstation higher order theory models, including the High Accuracy Satellite Drag Model. This paper de-scribes the BFMC algorithm, presents comparisons of BFMC P (sub c) estimates to those calculated using other methods, and discusses the implications for conjunction risk assessment.

Probabilities↗

NASA NTRS1996

Clusters of DNA induced by ionizing radiation: formation of short DNA fragments. I. Theoretical modeling

We have developed a general theoretical model for the interaction of ionizing radiation with chromatin. Chromatin is modeled as a 30-nm-diameter solenoidal fiber comprised of 20 turns of nucleosomes, 6 nucleosomes per turn. Charged-particle tracks are modeled by partitioning the energy deposition between primary track core, resulting from glancing collisions with 100 eV or less per event, and delta rays due to knock-on collisions involving energy transfers >100 eV. A Monte Carlo simulation incorporates damages due to the following molecular mechanisms: (1) ionization of water molecules leading to the formation of OH, H, eaq, etc.; (2) OH attack on sugar molecules leading to strand breaks: (3) OH attack on bases; (4) direct ionization of the sugar molecules leading to strand breaks; (5) direct ionization of the bases. Our calculations predict significant clustering of damage both locally, over regions up to 40 bp and over regions extending to several kilobase pairs. A characteristic feature of the regional damage predicted by our model is the production of short fragments of DNA associated with multiple nearby strand breaks. The shapes of the spectra of DNA fragment lengths depend on the symmetries or approximate symmetries of the chromatin structure. Such fragments have subsequently been detected experimentally and are reported in an accompanying paper (B. Rydberg, Radiat, Res. 145, 200-209, 1996) after exposure to both high- and low-LET radiation. The overall measured yields agree well quantitatively with the theoretical predictions. Our theoretical results predict the existence of a strong peak at about 85 bp, which represents the revolution period about the nucleosome. Other peaks at multiples of about 1,000 bp correspond to the periodicity of the particular solenoid model of chromatin used in these calculations. Theoretical results in combination with experimental data on fragmentation spectra may help determine the consensus or average structure of the chromatin fibers in mammalian DNA.

Non-NASA Center↗

NASA NTRS1973

Statistical Evaluation of Time Series Analysis Techniques

The performance of a modified version of NASA's multivariate spectrum analysis program is discussed. A multiple regression model was used to make the revisions. Performance improvements were documented and compared to the standard fast Fourier transform by Monte Carlo techniques.

Benignus, V. A.↗

NASA NTRS1985

Composition of primary cosmic rays near the bend from a study of hadrons in air showers at sea level

Data on hadrons in air showers arriving at sea level were studied to find sensitivity to primary cosmic ray composition. The rate of showers which satisfy minimum shower density and hadron energy requirements as well as the rate of showers containing hadrons delayed with respect to the electron shower front are compared to Monte Carlo simulations. The data on the rate of total triggers and delayed hadrons are compared to predicted rates for two models of primary composition. The data are consistent with models which require an increasing heavy nuclei fraction near 10 to the 15th power eV. The spectra which are consistent with the observed rate are also compared to the observed shower size spectrum at sea level and mountain level.

Mincer, A. I.↗