Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Techniques”

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 253 records · Page 14

On Photospheric Fluorescence and the Nature of the 17.62 Angstrom Feature in Solar X-ray Spectra

The identification of the emission line feature at 17.62 Angstroms in solar x-ray spectra is re-examined. Using a Monte Carlo technique, we compute a realistic theoretical upper limit to the observed Fe L-alpha photospheric fluorescent line strength caused by irradiation from an overlying corona. These calculations demonstrate that the photospheric Fe L-alpha characteristic line is much too weak to account for the observed 17.62 Angstrom line flux. Instead, we identify this line with the configuration interaction 2s2p3p2P-2s2p6 2S transition in Fe XVIII seen in Electron Beam Ion Trap spectra and predicted in earlier theoretical work on the Fe XVIII x-ray spectrum.

Drake, Jeremy J.↗

Response of silicon-Based Linear Energy Transfer Spectrometers

Silicon-based linear energy transfer (LET) telescope,(e. g., DOSTEL and RRMD) have recently been flown in space. LET spectra measured using tissue equivalent proportional counters show differences that need to be fully understood. A Monte Carlo technique based on: 1. radiation transport cluster intra-cascade model. 2. Landau-Vavilov distribution, 3. telescope geometry and detector coincidence & discriminator settings, 4. spacecraft shielding geometry, and 5. the external free space radiation environment, including recent albedo measurements, was developed.

Aman, A.↗

The Use of the Satellite Breakup Risk Assessment Model (SBRAM) to Characterize Collision Risk to Manned Spacecraft

NASA uses environment models such as ORDEM96 to characterize the long-term orbital debris collision hazard for spacecraft in LEO. Occasionally, however, there are breakups of satellites or rocket bodies that create enhanced collision hazard for a period of time. This enhanced collision hazard can pose increased risks to space operations - especially those involving manned missions where the tolerance for risk is very low. NASA has developed SBRAM to simulate the enhanced debris environment in the days and weeks that follow such a breakup. This simulation provides the kind of risk probabilities that can be used by mission planners to consider if changes are warranted for the mission. Announcements of breakups come to NASA from US Space Command as soon as they are identified. The pre-breakup orbit and time of breakup are used to determine the initial conditions of the explosion. SBRAM uses the latest explosion models developed at NASA to simulate a debris cloud for the breakup. The model uses a Monte Carlo technique to create a random debris cloud from the probability distributions in the breakup model. Each piece of debris randomly created in the cloud is propagated in a deterministic manner to include the effects of drag and other orbital perturbations. The detailed geometry of each simulated close approach to the target spacecraft is noted and logged and the collision probability is computed using an estimated probability density in down-range and cross-range positions of both the target spacecraft and debris object. The collision probability is computed from the overlap of these probability densities for each close-approach geometry and summed over all computed conjunctions. Cloud propagation runs over the desired time interval are then repeated until the scale of the collision risk can be estimated to a desired precision. This paper presents an overview of the SBRAM model and a number of examples, both real and hypothetical, to demonstrate its use. In addition, a number of different examples are shown how the data can be used by decision makers on issues such as spacecraft orientation and timing of EVAs.

Matney, Mark J.↗

Leonid Shower Probe of Aerothermochemistry in Meteoric Plasmas and Implication for the Origin of Life

The rarefied and high Mach number (up to 270) of the flow field of a typical meteoroid as it enters the Earth's atmosphere implies conditions of ablation and atmospheric chemistry that have proven to be as difficult to grasp as the proverbial shooting star. An airborne campaign was organized to study these processes during an intense Leonid shower. A probe of molecular band emission now demonstrates that the flash of light from a common meteor originates in the wake of the object rather than in the meteor head. A new theoretical approach using the direct simulation Monte Carlo technique demonstrates that the ablation process is critical in heating the air in that wake. Air molecules impinge on a dense cloud of ablated material in front of the meteoroid head into an extended wake that has the observed excitation temperatures. These processes determine what extraterrestrial materials may have been delivered to Earth at the time of the origin of life.

Jenniskens, Peter S. I.↗

Meteors: A Delivery Mechanism of Organic Matter to The Early Earth

All potential exogenous pre-biotic matter arrived to Earth by ways of our atmosphere, where much material was ablated during a luminous phase called 1. meteors" in rarefied flows of high (up to 270) Mach number. The recent Leonid showers offered a first glimpse into the elusive physical conditions of the ablation process and atmospheric chemistry associated with high-speed meteors. Molecular emissions were detected that trace a meteor's brilliant light to a 4,300 K warm wake rather than to the meteor's head. A new theoretical approach using the direct simulation by Monte Carlo technique identified the source-region and demonstrated that the ablation process is critical in the heating of the meteor's wake. In the head of the meteor, organic carbon appears to survive flash heating and rapid cooling. The temperatures in the wake of the meteor are just right for dissociation of CO and the formation of more complex organic compounds. The resulting materials could account for the bulk of pre-biotic organic carbon on the early Earth at the time of the origin of life.

Jenniskens, Peter↗

Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm

One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous mixing. A similar difficulty may be encountered in any method relying on global modifications of phase space.

Darve, Eric↗

Investigation of Nitride Morphology After Self-Aligned Contact Etch

Self-Aligned Contact (SAC) etch has emerged as a key enabling technology for the fabrication of very large-scale memory devices. However, this is also a very challenging technology to implement from an etch viewpoint. The issues that arise range from poor oxide etch selectivity to nitride to problems with post etch nitride surface morphology. Unfortunately, the mechanisms that drive nitride loss and surface behavior remain poorly understood. Using a simple langmuir site balance model, SAC nitride etch simulations have been performed and compared to actual etched results. This approach permits the study of various etch mechanisms that may play a role in determining nitride loss and surface morphology. Particle trajectories and fluxes are computed using Monte-Carlo techniques and initial data obtained from double Langmuir probe measurements. Etched surface advancement is implemented using a shock tracking algorithm. Sticking coefficients and etch yields are adjusted to obtain the best agreement between actual etched results and simulated profiles.

Hwang, Helen H.↗

ISO Key Project: Exploring The Full Range of Quasar/AGN Properties

While most of the work on this program has been completed, as previously reported, the portion of the program dealing with the sub topic of ISO LWS data analysis and reduction for the LWS Extragalactic Science Team and its leader, Dr. Howard Smith, is still active. This program in fact continues to generate results, and newly available computer modeling has extended the value of the datasets, As a result the team has requested and been granted an obtained a no-cost extension to this program, through December 31, 2003. The essence of the proposal is to perform ISO spectroscopic studies, including data analysis and modeling, of star formation regions using an ensemble of archival space-based data from the Infrared Space Observatory's Long Wavelength Spectrometer and Short Wavelength Spectrometer, but including as well some other spectroscopic data bases. Four kinds of regions are considered in the studies: (1) disks around more evolved objects; (2) young, low or high mass pre-main sequence stars in star formation regions; (3) star formation in external, bright IR galaxies; and (4) the galactic center. One prime focus of the program is the OH lines in the far infrared. The program has the following goals: (1) refine the data analysis of ISO observations, to obtain deeper and better SNR results on selected sources. The ISO data itself underwent "pipeline 10" reductions in early 2001, and additional "hands-on data reduction packages" were supplied by the ISO teams in 2001. The Fabry-Perot database in particularly sensitive to noise can slight calibration errors. (2) model the atomic and molecular line shapes, in particular the OH lines, using revised Monte-Carlo techniques developed by the SWAS team at the Center for Astrophysics; (3) attend scientific meetings and workshops; (4) do E&PO activities related to infrared astrophysics and/or spectroscopy.

Wilkes, Belinda↗

Robotic Mars Sample Return: Risk Assessment and Analysis Report

A comparison of the risk associated with two alternative scenarios for a robotic Mars sample return mission was conducted. Two alternative mission scenarios were identified, the Jet Propulsion Lab (JPL) reference Mission and a mission proposed by Johnson Space Center (JSC). The JPL mission was characterized by two landers and an orbiter, and a Mars orbit rendezvous to retrieve the samples. The JSC mission (Direct/SEP) involves a solar electric propulsion (SEP) return to earth followed by a rendezvous with the space shuttle in earth orbit. A qualitative risk assessment to identify and characterize the risks, and a risk analysis to quantify the risks were conducted on these missions. Technical descriptions of the competing scenarios were developed in conjunction with NASA engineers and the sequence of events for each candidate mission was developed. Risk distributions associated with individual and combinations of events were consolidated using event tree analysis in conjunction with Monte Carlo techniques to develop probabilities of mission success for each of the various alternatives. The results were the probability of success of various end states for each candidate scenario. These end states ranged from complete success through various levels of partial success to complete failure. Overall probability of success for the Direct/SEP mission was determined to be 66% for the return of at least one sample and 58% for the JPL mission for the return of at least one sample cache. Values were also determined for intermediate events and end states as well as for the probability of violation of planetary protection. Overall mission planetary protection event probabilities of occurrence were determined to be 0.002% and 1.3% for the Direct/SEP and JPL Reference missions respectively.

Lalk, Thomas R.↗

ISO Key Project: Exploring the Full Range of Quasar/Agn Properties

While most of the work on this program has been completed, as previously reported, the portion of the program dealing with the subtopic of ISO LWS data analysis and reduction for the LWS Extragalactic Science Team and its leader, Dr. Howard Smith, is still active. This program in fact continues to generate results, and newly available computer modeling has extended the value of the datasets. As a result the team requests a one-year no-cost extension to this program, through 31 December 2004. The essence of the proposal is to perform ISO spectroscopic studies, including data analysis and modeling, of star-formation regions using an ensemble of archival space-based data from the Infrared Space Observatory's Long Wavelength Spectrometer and Short Wavelength Spectrometer, but including as well some other spectroscopic databases. Four kinds of regions are considered in the studies: (1) disks around more evolved objects; (2) young, low or high mass pre-main sequence stars in star-formation regions; (3) star formation in external, bright IR galaxies; and (4) the galactic center. One prime focus of the program is the OH lines in the far infrared. The program has the following goals: 1) Refine the data analysis of ISO observations to obtain deeper and better SNR results on selected sources. The ISO data itself underwent 'pipeline 10' reductions in early 2001, and additional 'hands-on data reduction packages' were supplied by the ISO teams in 2001. The Fabry-Perot database is particularly sensitive to noise and slight calibration errors; 2) Model the atomic and molecular line shapes, in particular the OH lines, using revised Monte-Carlo techniques developed by the SWAS team at the Center for Astrophysics; 3) Attend scientific meetings and workshops; 4) Perform E&PO activities related to infrared astrophysics and/or spectroscopy.

Wilkes, Belinda↗

Multivariate Error Covariance Estimates by Monte-Carlo Simulation for Assimilation Studies in the Pacific Ocean

One of the most difficult aspects of ocean state estimation is the prescription of the model forecast error covariances. The paucity of ocean observations limits our ability to estimate the covariance structures from model-observation differences. In most practical applications, simple covariances are usually prescribed. Rarely are cross-covariances between different model variables used. Here a comparison is made between a univariate Optimal Interpolation (UOI) scheme and a multivariate OI algorithm (MvOI) in the assimilation of ocean temperature. In the UOI case only temperature is updated using a Gaussian covariance function and in the MvOI salinity, zonal and meridional velocities as well as temperature, are updated using an empirically estimated multivariate covariance matrix. Earlier studies have shown that a univariate OI has a detrimental effect on the salinity and velocity fields of the model. Apparently, in a sequential framework it is important to analyze temperature and salinity together. For the MvOI an estimation of the model error statistics is made by Monte-Carlo techniques from an ensemble of model integrations. An important advantage of using an ensemble of ocean states is that it provides a natural way to estimate cross-covariances between the fields of different physical variables constituting the model state vector, at the same time incorporating the model's dynamical and thermodynamical constraints as well as the effects of physical boundaries. Only temperature observations from the Tropical Atmosphere-Ocean array have been assimilated in this study. In order to investigate the efficacy of the multivariate scheme two data assimilation experiments are validated with a large independent set of recently published subsurface observations of salinity, zonal velocity and temperature. For reference, a third control run with no data assimilation is used to check how the data assimilation affects systematic model errors. While the performance of the UOI and MvOI is similar with respect to the temperature field, the salinity and velocity fields are greatly improved when multivariate correction is used, as evident from the analyses of the rms differences of these fields and independent observations. The MvOI assimilation is found to improve upon the control run in generating the water masses with properties close to the observed, while the UOI failed to maintain the temperature and salinity structure.

Borovikov, Anna↗

Filtering in Hybrid Dynamic Bayesian Networks

We implement a 2-time slice dynamic Bayesian network (2T-DBN) framework and make a 1-D state estimation simulation, an extension of the experiment in (v.d. Merwe et al., 2000) and compare different filtering techniques. Furthermore, we demonstrate experimentally that inference in a complex hybrid DBN is possible by simulating fault detection in a watertank system, an extension of the experiment in (Koller & Lerner, 2000) using a hybrid 2T-DBN. In both experiments, we perform approximate inference using standard filtering techniques, Monte Carlo methods and combinations of these. In the watertank simulation, we also demonstrate the use of 'non-strict' Rao-Blackwellisation. We show that the unscented Kalman filter (UKF) and UKF in a particle filtering framework outperform the generic particle filter, the extended Kalman filter (EKF) and EKF in a particle filtering framework with respect to accuracy in terms of estimation RMSE and sensitivity with respect to choice of network structure. Especially we demonstrate the superiority of UKF in a PF framework when our beliefs of how data was generated are wrong. Furthermore, we investigate the influence of data noise in the watertank simulation using UKF and PFUKD and show that the algorithms are more sensitive to changes in the measurement noise level that the process noise level. Theory and implementation is based on (v.d. Merwe et al., 2000).

Andersen, Morten Nonboe↗

Studying the Warm Layer and the Hardening Factor in Cygnus X-1

As the first dynamically determined black hole X-ray binary system, Cygnus X-1 has been studied extensively. However, its broadband spectrum observed with BeppoSax is still not well understood. Besides the soft excess described by the multi-color disk model (MCD), the power-law hard component and a broad excess feature above 10 keV (a disk reflection component), there is also an additional soft component around 1 keV, whose origin is not known currently. Here we propose that the additional soft component is due to the thermal Comptonization between the soft disk photons and a warm plasma cloud just above the disk, i.e., a warm layer. We use the Monte-Carlo technique to simulate this Compton scattering process and build a table model based on our simulation results. With this table model, we study the disk structure and estimate the hardening factor to the MCD component in Cygnus X-1.

Yao, Yangsen↗

Using Monte-Carlo Simulations to Study the Disk Structure in Cygnus X-1

As the first dynamically determined black hole X-ray binary system, Cygnus X-1 has been studied extensively. However, its broad-band spectra in hard state with BeppoSAX is still not well understood. Besides the soft excess described by the multi-color disk model (MCD), the power- law component and a broad excess feature above 10 keV (disk reflection component), there is also an additional soft component around 1 keV, whose origin is not known currently.We propose that the additional soft component is due to the thermal Comptonization process between the s oft disk photon and the warm plasma cloud just above the disk.i.e., a warm layer. We use Monte-Carlo technique t o simulate this Compton scattering process and build several table models based on our simulation results.

Yao, Y.↗

Cosmological Constraints from Sunyaev-Zeldovich Effect and X-ray data for 37 Galaxy Clusters

We present a sample of 37 clusters of galaxies (0.14<=z<=0.89) for which we have Sunyaev-Zeldovich Effect data and Chandra X-ray imaging and spectroscopy. We use a Markov chain Monte Carlo technique to determine the spatial properties and surface brightness of the X-ray emitting gas, the X-ray temperature and abundance, and the SZE decrement for each cluster. From these observables, we infer each cluster's total mass, gas mass fraction, and angular diameter distance. The sample properties are used to constrain the Hubble parameter and Omega Matter, and scaling relations between integrated SZE flux and cluster mass and temperature are constructed. We also investigate the cosmological constraints that could be obtained if the X-ray/SZE sample were extended to a redshift of z approximately 2.

Bonamente, Massimiliano↗

Explore Full Range of QSO/AGN Properties

The goal of the proposal is to perform ISO spectroscopic studies, including data analysis and modeling, of star formation regions using an ensemble of archival space-based data from the Infrared Space Observatory s Long Wavelength Spectrometer and Short Wavelength Spectrometer, but including as well some other spectroscopic databases. Four kinds of regions are considered in the studies: (1) disks around more evolved objects; (2) young, low or high mass pre-main sequence stars in star formation regions; (3) star formation in external, bright IR galaxies; and (4) the galactic center. One prime focus of the program is the OH lines in the far infrared. The program had the following goals: 1) Refine the data analysis of IS0 observations to obtain deeper and better SNR results on selected sources. The IS0 data itself underwent "pipeline 10" reductions in early 2001, and additional "hands-on data reduction packages" were supplied by the IS0 teams in 2001. The Fabry-Perot database is particularly sensitive to noise and slight calibration errors. 2) Model the atomic and molecular line shapes, in particular the OH lines, using revised monte- carlo techniques developed by the SWAS team at the Center for Astrophysics; 3) Attend scientific meetings and workshops; 4) Do E&PO activities related to infrared astrophysics and/or spectroscopy.

Oliversen, Ronald↗

Managing Lunar and Mars Mission Radiation Risks: Cancer Risks, Uncertainties, and Shielding Effectiveness - Part 1

This document addresses calculations of probability distribution functions (PDFs) representing uncertainties in projecting fatal cancer risk from galactic cosmic rays (GCR) and solar particle events (SPEs). PDFs are used to test the effectiveness of potential radiation shielding approaches. Monte-Carlo techniques are used to propagate uncertainties in risk coefficients determined from epidemiology data, dose and dose-rate reduction factors, quality factors, and physics models of radiation environments. Competing mortality risks and functional correlations in radiation quality factor uncertainties are treated in the calculations. The cancer risk uncertainty is about four-fold for lunar and Mars mission risk projections. For short-stay lunar missins (<180 d), SPEs present the most significant risk, but one effectively mitigated by shielding. For long-duration (>180 d) lunar or Mars missions, GCR risks may exceed radiation risk limits. While shielding materials are marginally effective in reducing GCR cancer risks because of the penetrating nature of GCR and secondary radiation produced in tissue by relativisitc particles, polyethylene or carbon composite shielding cannot be shown to significantly reduce risk compared to aluminum shielding. Therefore, improving our knowledge of space radiobiology to narrow uncertainties that lead to wide PDFs is the best approach to ensure radiation protection goals are met for space exploration.

Cucinotta, Francis A.↗

Determination of the Cosmic Distance Scale from Sunyaev-Zel'dovich Effect and Chandra X-ray Measurements of High Redshift Galaxy Clusters

We determine the distance to 38 clusters of galaxies in the redshift range 0.14 less than or equal to z less than or equal to 0.89 using X-ray data from Chandra and Sunyaev-Zeldovich Effect data from the Owens Valley Radio Observatory and the Berkeley-Illinois-Maryland Association interferometric arrays. The cluster plasma and dark matter distributions are analyzed using a hydrostatic equilibrium model that accounts for radial variations in density, temperature and abundance, and, the statistical and systematic errors of this method are quantified. The analysis is performed via a Markov chain Monte Carlo technique that provides simultaneous estimation of all model parameters. W

Bonamente, Massimiliano↗