SEARCH · Search NASA

Results for “Variational Monte Carlo”

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

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Two-dimensional radiative transfer in cloudy atmospheres - The spherical harmonic spatial grid method

A new two-dimensional monochromatic method that computes the transfer of solar or thermal radiation through atmospheres with arbitrary optical properties is described. The model discretizes the radiative transfer equation by expanding the angular part of the radiance field in a spherical harmonic series and representing the spatial part with a discrete grid. The resulting sparse coupled system of equations is solved iteratively with the conjugate gradient method. A Monte Carlo model is used for extensive verification of outgoing flux and radiance values from both smooth and highly variable (multifractal) media. The spherical harmonic expansion naturally allows for different levels of approximation, but tests show that the 2D equivalent of the two-stream approximation is poor at approximating variations in the outgoing flux. The model developed here is shown to be highly efficient so that media with tens of thousands of grid points can be computed in minutes. The large improvement in efficiency will permit quick, accurate radiative transfer calculations of realistic cloud fields and improve our understanding of the effect of inhomogeneity on radiative transfer in cloudy atmospheres.

Evans, K. F.↗

NASA NTRS1990

Use of random Martian atmosphere to evaluate potential entry guidance schemes

A random Martian atmosphere was developed and was used with three guidance schemes to determine the effect of random density variations on the guidance. This random atmosphere was shown to be useful for testing the robustness of guidance schemes for vehicles encountering random disturbances during aerobraking for capture into planetary orbit. Levels of disturbance that could be tolerated and areas where performance could be improved were established. The need for Monte Carlo studies to define the excursion boundaries of capture orbit parameters was indicated.

Tomlinson, Barbara S.↗

NASA NTRS1985

Fine guidance sensor design optimization for Space Infrared Telescope Facility (SIRTF)

The present paper is concerned with the general issues of design optimization of a staring mosaic sensor for guide-star tracking, taking into account the choice of the proper interpolation algorithm. It was found that the final choice of algorithm depends on the severity of the positional accuracy requirement and the expected variation in the blur spot shape. The considered track algorithms estimate the position either of the peak or of the centroid of the image of a point source. Attention is given to aspects of problem definition, the total effect of error sources on algorithm accuracy, a parametric analysis of performance behavior, Monte Carlo results, and the design of the SIRTF Fine Guidance Sensor (FGS). The algorithm error, the signal-to-noise error, and the nonuniformity error are examined.

Grossman, S. B.↗

NASA NTRS1974

Monte Carlo analysis of the solar-wind modulation of galactic C-N-O at solar maximum

A previous (Moss and Giuli, 1971) Monte Carlo analysis of the solar-wind modulation of galactic cosmic-ray rigidity spectra is extended to include particles in the C-N-O group. A solar cavity 'radius' of 7 AU is found to be consistent with observations of the C-N-O component at solar maximum. The variation in modulation between solar maximum and solar minimum may be explained solely by variation of the solar cavity size, with no change in the scattering mean free path.

Moss, T. A.↗

NASA NTRS1977

Backscatter of hard X-rays in the solar atmosphere

The solar photosphere backscatters a substantial fraction of the hard X rays from solar flares incident upon it. This reflection was studied using a Monte Carlo simulation which takes into account Compton scattering and photo-electric absorption. Both isotropic and anisotropic X ray sources are considered. The bremsstrahlung from an anisotropic distribution of electrons are evaluated. By taking the reflection into account, the inconsistency is removed between recent observational data regarding the center-to-limb variation of solar X ray emission and the predictions of models in which accelerated electrons are moving down toward the photosphere.

Bai, T.↗

NASA NTRS1974

Model of an anisotropy solar proton event

The anisotropy of the particle distribution and its variation with time at 1 AU early in a solar cosmic-ray event can provide information on the pitch-angle scattering of the particles in the interplanetary medium. One proton event (7.6-55 MeV) is described in which the anisotropy remained large well into the decay of the omnidirectional intensity. A Monte Carlo technique which gives the pitch-angle distribution is employed to investigate two models put forward to explain the sustained anisotropy. It is shown that the observed event is consistent with a model in which the injection of particles at the sun reaches a peak 1.5 hours after flare maximum before decaying with an e-folding time of 7 hours. The uniform diffusion coefficient in this model corresponds to a classical mean free path of 1.6 AU.

Palmer, I. D.↗

NASA NTRS1987

Hydrogen Balmer alpha intensity distributions and line profiles from multiple scattering theory using realistic geocoronal models

The H Balmer alpha nightglow is investigated by using Monte Carlo models of asymmetric geocoronal atomic hydrogen distributions as input to a radiative transfer model of solar Lyman-beta radiation in the thermosphere and atmosphere. It is shown that it is essential to include multiple scattering of Lyman-beta radiation in the interpretation of Balmer alpha airglow data. Observations of diurnal variation in the Balmer alpha airglow showing slightly greater intensities in the morning relative to evening are consistent with theory. No evidence is found for anything other than a single sinusoidal diurnal variation of exobase density. Dramatic changes in effective temperature derived from the observed Balmer alpha line profiles are expected on the basis of changing illumination conditions in the thermosphere and exosphere as different regions of the sky are scanned.

Anderson, D. E., Jr.↗

NASA NTRSDate not supplied

Constraints on Triple Oxygen Isotope Kinetics

Triple oxygen isotope methods (e.g., using δ 18 O, and Δ′ 17 O) have been used to classify meteorites, to constrain atmospheric chemistry, and to estimate primary productivity, among other uses. Such utility relies on the mass dependency of most chemical reactions: δ 17 O is proportional to δ 18 O by roughly a factor of 0.5. However, in detail, mass-dependency factors vary. While constraints on mass dependency are well established for equilibrium isotope effects, similar constraints for kinetic isotope effects are lacking. Using an ab initio Monte Carlo approach (see envelopes in Figure 1), we find the range of kinetic isotope effects (i.e., resulting from unidirectional chemical reactions) contains and exceeds the range of accessible equilibrium isotope effects for the 16 O 17 O 18 O system (Hayles and Killingsworth, 2022). Accessible kinetic isotope effects can yield Δ′ 17 O variations as large 0.8 ‰ from a single kinetic fractionation step. Reactions that yield the largest deviations from equilibrium are those with the largest molecular mass difference between decomposition fragments, typically the two products such as a gas and solid generated from the thermal decomposition of an initial solid. Comparison of DFT transition state model results against previously published experiments on the thermal decomposition of calcite (CaCO 3 ) (Miller et al., 2002) and dehydroxylation of brucite (Mg(OH) 2 ) (Clayton and Mayeda, 2009) show they are matched with mass-dependent kinetic isotope effects. This approach can be applied to any system with more than two isotopes.

Oxygen isotopes↗

NASA NTRS1993

Monte Carlo analysis of the Titan III/Transfer Orbit Stage guidance system for the Mars Observer mission

An important part of space launch vehicle mission planning for a planetary mission is the integrated analysis of guidance and performance dispersions for both booster and upper stage vehicles. For the Mars Observer mission, an integrated trajectory analysis was used to maximize the scientific payload and to minimize injection errors by optimizing the energy management of both vehicles. This was accomplished by designing the Titan III booster vehicle to inject into a hyperbolic departure plane, and the Transfer Orbit Stage (TOS) to correct any booster dispersions. An integrated Monte Carlo analysis of the performance and guidance dispersions of both vehicles provided sensitivities, an evaluation of their guidance schemes and an injection error covariance matrix. The polynomial guidance schemes used for the Titan III variable flight azimuth computations and the TOS solid rocket motor ignition time and burn direction derivations accounted for a wide variation of launch times, performance dispersions, and target conditions. The Mars Observer spacecraft was launched on 25 September 1992 on the Titan III/TOS vehicle. The post flight analysis indicated that a near perfect park orbit injection was achieved, followed by a trans-Mars injection with less than 2sigma errors.

Bell, Stephen C.↗

NASA NTRS1997

Probabilistic Micromechanics and Macromechanics for Ceramic Matrix Composites

The properties of ceramic matrix composites (CMC's) are known to display a considerable amount of scatter due to variations in fiber/matrix properties, interphase properties, interphase bonding, amount of matrix voids, and many geometry- or fabrication-related parameters, such as ply thickness and ply orientation. This paper summarizes preliminary studies in which formal probabilistic descriptions of the material-behavior- and fabrication-related parameters were incorporated into micromechanics and macromechanics for CMC'S. In this process two existing methodologies, namely CMC micromechanics and macromechanics analysis and a fast probability integration (FPI) technique are synergistically coupled to obtain the probabilistic composite behavior or response. Preliminary results in the form of cumulative probability distributions and information on the probability sensitivities of the response to primitive variables for a unidirectional silicon carbide/reaction-bonded silicon nitride (SiC/RBSN) CMC are presented. The cumulative distribution functions are computed for composite moduli, thermal expansion coefficients, thermal conductivities, and longitudinal tensile strength at room temperature. The variations in the constituent properties that directly affect these composite properties are accounted for via assumed probabilistic distributions. Collectively, the results show that the present technique provides valuable information about the composite properties and sensitivity factors, which is useful to design or test engineers. Furthermore, the present methodology is computationally more efficient than a standard Monte-Carlo simulation technique; and the agreement between the two solutions is excellent, as shown via select examples.

Murthy, Pappu L. N.↗

NASA NTRS2005

Primal and Dual Integrated Force Methods Used for Stochastic Analysis

At the NASA Glenn Research Center, the primal and dual integrated force methods are being extended for the stochastic analysis of structures. The stochastic simulation can be used to quantify the consequence of scatter in stress and displacement response because of a specified variation in input parameters such as load (mechanical, thermal, and support settling loads), material properties (strength, modulus, density, etc.), and sizing design variables (depth, thickness, etc.). All the parameters are modeled as random variables with given probability distributions, means, and covariances. The stochastic response is formulated through a quadratic perturbation theory, and it is verified through a Monte Carlo simulation.

Patnaik, Surya N.↗

NASA NTRSDate not supplied

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

NASA NTRSDate not supplied

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

microstructure↗

NASA NTRS1997

Probabilistic Mesomechanical Fatigue Model

A probabilistic mesomechanical fatigue life model is proposed to link the microstructural material heterogeneities to the statistical scatter in the macrostructural response. The macrostructure is modeled as an ensemble of microelements. Cracks nucleation within the microelements and grow from the microelements to final fracture. Variations of the microelement properties are defined using statistical parameters. A micromechanical slip band decohesion model is used to determine the crack nucleation life and size. A crack tip opening displacement model is used to determine the small crack growth life and size. Paris law is used to determine the long crack growth life. The models are combined in a Monte Carlo simulation to determine the statistical distribution of total fatigue life for the macrostructure. The modeled response is compared to trends in experimental observations from the literature.

Tryon, Robert G.↗

NASA NTRS2003

The Development of Design Tools for Fault Tolerant Quantum Dot Cellular Automata Based Logic

We are developing software to explore the fault tolerance of quantum dot cellular automata gate architectures in the presence of manufacturing variations and device defects. The Topology Optimization Methodology using Applied Statistics (TOMAS) framework extends the capabilities of the A Quantum Interconnected Network Array Simulator (AQUINAS) by adding front-end and back-end software and creating an environment that integrates all of these components. The front-end tools establish all simulation parameters, configure the simulation system, automate the Monte Carlo generation of simulation files, and execute the simulation of these files. The back-end tools perform automated data parsing, statistical analysis and report generation.

Armstrong, Curtis D.↗

NASA NTRS2008

DSN Array Simulator

The DSN Array Simulator (wherein 'DSN' signifies NASA's Deep Space Network) is an updated version of software previously denoted the DSN Receive Array Technology Assessment Simulation. This software (see figure) is used for computational modeling of a proposed DSN facility comprising user-defined arrays of antennas and transmitting and receiving equipment for microwave communication with spacecraft on interplanetary missions. The simulation includes variations in spacecraft tracked and communication demand changes for up to several decades of future operation. Such modeling is performed to estimate facility performance, evaluate requirements that govern facility design, and evaluate proposed improvements in hardware and/or software. The updated version of this software affords enhanced capability for characterizing facility performance against user-defined mission sets. The software includes a Monte Carlo simulation component that enables rapid generation of key mission-set metrics (e.g., numbers of links, data rates, and date volumes), and statistical distributions thereof as functions of time. The updated version also offers expanded capability for mixed-asset network modeling--for example, for running scenarios that involve user-definable mixtures of antennas having different diameters (in contradistinction to a fixed number of antennas having the same fixed diameter). The improved version also affords greater simulation fidelity, sufficient for validation by comparison with actual DSN operations and analytically predictable performance metrics.

Tikidjian, Raffi↗

NASA NTRS1973

Reflection of a shock wave from a thermally accommodating wall - Molecular simulation.

Reflection of a plane shock wave from a wall has been simulated on a microscopic scale using a direct simulation Monte Carlo technique of the type developed by Bird. A monatomic gas model representing argon was used to describe the fluid medium and a simple one-parameter accommodation coefficient model was used to describe the gas-surface interaction. The influence of surface accommodation was studied parametrically by varying the accommodation coefficient from zero to one. Results are presented showing the temporal variations of flow field density, and mass, momentum, and energy fluxes to the wall during the shock wave reflection process. The energy flux was used to determine the wall temperature history. Comparisons with experiment are found to be satisfactory where data are available.

Deiwert, G. S.↗

NASA NTRS2007

The UDF05 Follow-up of the HUDF: I. The Faint-End Slope of the Lyman-Break Galaxy Population at zeta approx. 5

We present the UDF05 project, a HST Large Program of deep ACS (F606W, F775W, F850LP, and NICMOS (Fll0W, Fl60W) imaging of three fields, two of which coincide with the NICP1-4 NICMOS parallel observations of the Hubble Ultra Deep Field (HUDF). In this first paper we use the ACS data for the NICP12 field, as well as the original HUDF ACS data, to measure the UV Luminosity Function (LF) of z approximately 5 Lyman Break Galaxies (LBGs) down to very faint levels. Specifically, based on a V - i, i - z selection criterion, we identify a sample of 101 and 133 candidate z approximately 5 galaxies down to z(sub 850) = 28.5 and 29.25 magnitudes in the NICP12 field and in the HUDF, respectively. Using an extensive set of Monte Carlo simulations we derive corrections for observational biases and selection effects, and construct the rest-frame 1400 Angstroms LBG LF over the range M(sub 1400) = [-22.2, -17.1], i.e. down to approximately 0.04 L(sub *) at z = 5. We show that: (i) Different assumptions for the SED distribution of the LBG population, dust properties and intergalactic absorption result in a 25% variation in the number density of LBGs at z = 5 (ii) Under consistent assumptions for dust properties and intergalactic absorption, the HUDF is about 30% under-dense in z = 5 LBGs relative to the NICP12 field, a variation which is well explained by cosmic variance; (iii) The faint-end slope of the LF is independent of the specific assumptions for the input physical parameters, and has a value of alpha approximately -1.6, similar to the faint-end slope of the LF that has been measured for LBGs at z = 3 and z = 6. Our study therefore supports no variation in the faint-end of the LBG LF over the whole redshift range z = 3 to z = 6. The comparison with theoretical predictions suggests that (a,) the majority of the stars in the z = 5 LBG population are produced with a Top-Heavy IMF in merger-driven starbursts, and that (b) possibly, either the fraction of stellar mass produced in starburst, or the fraction of high mass stars in the bursts is increased towards the bright end of the LF.

Oesch, P. A.↗