Search NASA⌕ Search

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.

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

At least 109 records · Page 6

Preliminary results of a Monte Carlo simulation of burst directions by BATSE

The effect of varying BATSE trigger parameters on the modeled BATSE detection rate is investigated numerically. Quantitative estimates are obtained of the systematic effects generated by the complex interplay of trigger times, variations in the signal-to-noise ratio for burst detection as well as the time scale over which backgrounds are determined, and the number of required detections above threshold. It is found, in particular that employing a third detection, at a reduced signal-to-noise ratio, to resolve a burst does not increase the detection rate. It is also shown that the BATSE detection rate is insensitive to variations in the time scale over which the background is determined.

Higdon, J. C.↗

Seasonal Variations in Mercury's Dayside Calcium Exosphere

The Mercury Atmospheric and Surface Composition Spectrometer on the MESSENGER spacecraft has observed calcium emission in Mercury's exosphere on a near-daily basis since March 2011. During MESSENGER's primary and first extended missions (March 2011 - March 2013) the dayside calcium exosphere was measured over eight Mercury years. We have simulated these data with a Monte Carlo model of exospheric source processes to show that (a) there is a persistent source of energetic calcium located in the dawn equatorial region, (b) there is a seasonal dependence in the calcium source rate, and (c) there are no obvious year-to-year variations in the near-surface dayside calcium exosphere. Although the precise mechanism responsible for ejecting the calcium has not yet been determined, the most likely process is the dissociation of Ca-bearing molecules produced in micrometeoroid impact plumes to form energetic, escaping calcium atoms.

Exosphere↗

L2-Charged Particle Environment (L2-CPE)Low Energy Radiation Fluence Model

The L2 Charged Particle Environment (L2-CPE) model provides estimates of number flux and fluence for the low energy electron, proton, and alpha particle populations in the near Earth solar wind and the Earth’s distant magnetosheath and magnetotail. The model is an engineering tool for specifying radiation environments over an energy range from a few eV to a few MeV of importance to surface dose and radiation damage to thin space exposed materials and is intended for use in space system design applications. Mission fluences are obtained by simulating a spacecraft flight trajectory through time-dependent bow shock and magnetopause boundaries with dimensions and orientations driven by solar wind parameters. Monte Carlo sampling of flux environments within individual plasma regimes and/or fluence accumulated along a fight trajectory through multiple regimes is used to determine statistical variations (means and extremes) of the differential number flux and fluence environments for each of the three charged particle species. Model output is differential (in energy) number fluence and flux for the three particle species in units of particles/cm2-keV and particles/cm2-sec-keV, respectively. Users can select the energy range of interest but the current version of the code (L2-CPE Version 1.4.2d) is limited to an energy range of 1 eV to 10 MeV. Flux is integrated over angle to give flux and fluence to surfaces in the ±XGSE, ±YGSE, and/or ±ZGSE directions. Figure 1 shows the opening screen from the L2-CPE graphical user interface (GUI) with options for flux, fluence, plotting model output. The flux scene generate is not currently implemented in the GUI.

Joseph I Minow↗

Dynamic performance of an aero-assist spacecraft - AFE

Dynamic performance of the Aero-assist Flight Experiment (AFE) spacecraft was investigated using a high-fidelity 6-DOF simulation model. Baseline guidance logic, control logic, and a strapdown navigation system to be used on the AFE spacecraft are also modeled in the 6-DOF simulation. During the AFE mission, uncertainties in the environment and the spacecraft are described by an error space which includes both correlated and uncorrelated error sources. The principal error sources modeled in this study include navigation errors, initial state vector errors, atmospheric variations, aerodynamic uncertainties, center-of-gravity off-sets, and weight uncertainties. The impact of the perturbations on the spacecraft performance is investigated using Monte Carlo repetitive statistical techniques. During the Solid Rocket Motor (SRM) deorbit phase, a target flight path angle of -4.76 deg at entry interface (EI) offers very high probability of avoiding SRM casing skip-out from the atmosphere. Generally speaking, the baseline designs of the guidance, navigation, and control systems satisfy most of the science and mission requirements.

Chang, Ho-Pen↗

Global-scale modes of surface temperature variability on interannual to century timescales

Using 100 years of global temperature anomaly data, we have performed a singluar value decomposition of temperature variations in narrow frequency bands to isolate coherent spatio-temporal modes of global climate variability. Statistical significance is determined from confidence limits obtained by Monte Carlo simulations. Secular variance is dominated by a globally coherent trend; with nearly all grid points warming in phase at varying amplitude. A smaller, but significant, share of the secular variance corresponds to a pattern dominated by warming and subsequent cooling in the high latitude North Atlantic with a roughly centennial timescale. Spatial patterns associated with significant peaks in variance within a broad period range from 2.8 to 5.7 years exhibit characteristic El Nino-Southern Oscillation (ENSO) patterns. A recent transition to a regime of higher ENSO frequency is suggested by our analysis. An interdecadal mode in the 15-to-18 years period and a mode centered at 7-to-8 years period both exhibit predominantly a North Atlantic Oscillation (NAO) temperature pattern. A potentially significant decadal mode centered on 11-to-12 years period also exhibits an NAO temperature pattern and may be modulated by the century-scale North Atlantic variability.

Mann, Michael E.↗

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↗

Measurement and Simulation of the Variation in Proton-Induced Energy Deposition in Large Silicon Diode Arrays

The proton induced charge deposition in a well characterized silicon P-i-N focal plane array is analyzed with Monte Carlo based simulations. These simulations include all physical processes, together with pile up, to accurately describe the experimental data. Simulation results reveal important high energy events not easily detected through experiment due to low statistics. The effects of each physical mechanism on the device response is shown for a single proton energy as well as a full proton space flux.

Howe, Christina L.↗

Simulations of sonic boom ray tube area fluctuations for propagation through atmospheric turbulence including caustics via a Monte Carlo method

A theory which gives statistical predictions for how often sonic booms propagating through the earth's turbulent boundary layer will encounter caustics, given the spectral properties of the atmospheric turbulence, is outlined. The theory is simple but approximately accounts for the variation of ray tube areas along ray paths. This theory predicts that the variation of ray tube areas is determined by the product of two similar area factors, psi (x) and phi (x), each satisfying a generic harmonic oscillator equation. If an area factor increases the peak acoustic pressure decreases, and if the factor decreases the peak acoustic pressure increases. Additionally, if an area factor decreases to zero and becomes negative, the ray has propagated through a caustic, which contributes a phase change of 90 degrees to the wave. Thus, it is clear that the number of times that a sonic boom wave passes through a caustic should be related to the distorted boom waveform received on the ground. Examples are given based on a characterization of atmospheric turbulence due to the structure function of Tatarski as modified by Crow.

Sparrow, Victor W.↗

A probabilistic model of a porous heat exchanger

This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

Agrawal, O. P.↗

A Sample/Jitter Monte Carlo Technique for Main Parachute Loads Predictions

Models for Orion parachute performance are based on reconstructions of the Capsule Parachute Assembly System (CPAS) drop test campaign and were documented in the CPAS “Model Memo.” Experience with similar Commercial Crew Program (CCP) parachute systems resulted in some updates to the Orion models in preparation for Artemis missions. The reefing cutter dispersion model for the drogues and mains had been overly-conservative by producing wide timing differences within clusters. A higher-fidelity timing model was generated by separating out in-lot variation and temperature effects. The main parachute inflation model had accounted for some correlations between parameters using complicated 2-D geometric bounding, but the results tended to exaggerate individual peak loads from fast (leading) inflations and under-emphasize actual lagging experience. Several flight tests were reconstructed again with an emphasis on matching peak load magnitudes using a search algorithm. A simpler method for generating inflation parameters uses the 3-D correlated reconstructed “samples” with some random “jitter” applied. Dispersed Monte Carlo inputs are then checked against flight test data to evaluate whether they represent reality.

parachutes↗

Independent pixel and Monte Carlo estimates of stratocumulus albedo

Monte Carlo radiative transfer methods are employed here to estimate the plane-parallel albedo bias for marine stratocumulus clouds. This is the bias in estimates of the mesoscale-average albedo, which arises from the assumption that cloud liquid water is uniformly distributed. The authors compare such estimates with those based on a more realistic distribution generated from a fractal model of marine stratocumulus clouds belonging to the class of 'bounded cascade' models. In this model the cloud top and base are fixed, so that all variations in cloud shape are ignored. The model generates random variations in liquid water along a single horizontal direction, forming fractal cloud streets while conserving the total liquid water in the cloud field. The model reproduces the mean, variance, and skewness of the vertically integrated cloud liquid water, as well as its observed wavenumber spectrum, which is approximately a power law. The Monte Carlo method keeps track of the three-dimensional paths solar photons take through the cloud field, using a vectorized implementation of a direct technique. The simplifications in the cloud field studied here allow the computations to be accelerated. The Monte Carlo results are compared to those of the independent pixel approximation, which neglects net horizontal photon transport. Differences between the Monte Carlo and independent pixel estimates of the mesoscale-average albedo are on the order of 1% for conservative scattering, while the plane-parallel bias itself is an order of magnitude larger. As cloud absorption increases, the independent pixel approximation agrees even more closely with the Monte Carlo estimates. This result holds for a wide range of sun angles and aspect ratios. Thus, horizontal photon transport can be safely neglected in estimates of the area-average flux for such cloud models. This result relies on the rapid falloff of the wavenumber spectrum of stratocumulus, which ensures that the smaller-scale variability, where the radiative transfer is more three-dimensional, contributes less to the plane-parallel albedo bias than the larger scales, which are more variable. The lack of significant three-dimensional effects also relies on the assumption of a relatively simple geometry. Even with these assumptions, the independent pixel approximation is accurate only for fluxes averaged over large horizontal areas, many photon mean free paths in diameter, and not for local radiance values, which depend strongly on the interaction between neighboring cloud elements.

Cahalan, Robert F.↗

Using probabilistic analysis to assess the reliability of predicted SRB aft-skirt stresses

Probabilistic failure analysis is a tool to predict the reliability of a part or system. Probabalistic techniques were used to predict critical stresses which occur in the solid rocket booster aft-skirt during main engine buildup, immediately prior to lift-off. More than any other hold down post (HDP) load component, the Z loads are sensitive to variations in strains and calibration constants. Also, predicted aft-skirt stresses are strongly affected by HDP load variations. Therefore, the instrumented HDP are not effective load transducers for Z loads, and, when used with aft skirt stress indicator equations, yield estimates with large uncertainty. Monte Carlo simulation proved to be a straight forward way of studying the overlapping effects of multiple parameters on predicted equipment performance. An advantage of probabilistic analysis is the degree of uncertainty of each parameter as stated explicitly by its probability distribution. It was noted, however, that the choice of parameter distribution had a large effect on the simulation results. Many times these distributions must be assumed. The engineer who is designing the part should be responsible for the choice of parameter distribution.

Richardson, James A.↗

Robust Stability Analysis of the Space Launch System Control Design: A Singular Value Approach

Classical stability analysis consists of breaking the feedback loops one at a time and determining separately how much gain or phase variations would destabilize the stable nominal feedback system. For typical launch vehicle control design, classical control techniques are generally employed. In addition to stability margins, frequency domain Monte Carlo methods are used to evaluate the robustness of the design. However, such techniques were developed for Single-Input-Single-Output (SISO) systems and do not take into consideration the off-diagonal terms in the transfer function matrix of Multi-Input-Multi-Output (MIMO) systems. Robust stability analysis techniques such as H(sub infinity) and mu are applicable to MIMO systems but have not been adopted as standard practices within the launch vehicle controls community. This paper took advantage of a simple singular-value-based MIMO stability margin evaluation method based on work done by Mukhopadhyay and Newsom and applied it to the SLS high-fidelity dynamics model. The method computes a simultaneous multi-loop gain and phase margin that could be related back to classical margins. The results presented in this paper suggest that for the SLS system, traditional SISO stability margins are similar to the MIMO margins. This additional level of verification provides confidence in the robustness of the control design.

Pei, Jing↗

Application of the Probabilistic Dynamic Synthesis Method to the Analysis of a Realistic Structure

The Probabilistic Dynamic Synthesis method is a new technique for obtaining the statistics of a desired response engineering quantity for a structure with non-deterministic parameters. The method uses measured data from modal testing of the structure as the input random variables, rather than more "primitive" quantities like geometry or material variation. This modal information is much more comprehensive and easily measured than the "primitive" information. The probabilistic analysis is carried out using either response surface reliability methods or Monte Carlo simulation. A previous work verified the feasibility of the PDS method on a simple seven degree-of-freedom spring-mass system. In this paper, extensive issues involved with applying the method to a realistic three-substructure system are examined, and free and forced response analyses are performed. The results from using the method are promising, especially when the lack of alternatives for obtaining quantitative output for probabilistic structures is considered.

Brown, Andrew M.↗

Application of the Probabilistic Dynamic Synthesis Method to Realistic Structures

The Probabilistic Dynamic Synthesis method is a technique for obtaining the statistics of a desired response engineering quantity for a structure with non-deterministic parameters. The method uses measured data from modal testing of the structure as the input random variables, rather than more "primitive" quantities like geometry or material variation. This modal information is much more comprehensive and easily measured than the "primitive" information. The probabilistic analysis is carried out using either response surface reliability methods or Monte Carlo simulation. In previous work, the feasibility of the PDS method applied to a simple seven degree-of-freedom spring-mass system was verified. In this paper, extensive issues involved with applying the method to a realistic three-substructure system are examined, and free and forced response analyses are performed. The results from using the method are promising, especially when the lack of alternatives for obtaining quantitative output for probabilistic structures is considered.

Brown, Andrew M.↗

A Monte Carlo analysis of the effects of instrumentation errors on aircraft parameter identification

An output error estimation algorithm was used to evaluate the effects of both static and dynamic instrumentation errors on the estimation of aircraft stability and control parameters. A Monte Carlo analysis, using simulated cruise flight data, was performed for a high performance military aircraft, a large commercial transport, and a small general-aviation aircraft. The effects of variations in the information content of the flight data, resulting from two different choices of control input maneuvers, were also determined. The results indicate that unmodeled instrumentation errors can cause inaccuracies in the estimated parameters which are comparable to their nominal values. Control input errors and angular accelerometer lags were found to be most significant of the instrumentation errors evaluated, and the perturbations they produce are much larger than those arising from the combined effects of static errors and white noise in the output response measurements.

Bryant, W. H.↗

Quantifying the Sensitivity of Condition Incidence Parameters in the Evidence Library

One approach to quantifying spaceflight risk at NASA makes use event driven probabilistic techniques. The Medical Extensible Dynamic Probabilistic Risk Assessment Tool (MEDPRAT) is such a tool that estimates medical risk metrics via simulation and enables optimization of medical resources subject to mission constraints [1]. Previous analyses have informed medical set composition, exercise countermeasures, and water intake, where each analysis quantifies the risk associated with proposed variations in system design. As future mission profiles extend beyond Low-Earth Orbit (LEO) and lengthen in duration, understanding these risks and contributing factors is critical. MEDPRAT employs Monte Carlo sampling techniques to simulate missions and track the occurrence of medical events. These events follow fault-tree-like progressions through levels of severity and mitigation via medical treatment to many possible outcomes and these are reported throughout the mission. Making this possible, are the medical databases that contain evidence gathered by the Human Research Program (HRP). Quantifying the impact of uncertainty or variability in the input data is an important step in evaluating the credibility of modeling and simulation results. In this work, we investigate the sensitivity of medical risk metrics with respect to the condition incidence parameters within the Evidence Library (EL) [2] as the medical database input for MEDPRAT. The medical conditions, contained in the EL, are equipped with incidence rates that describe the likelihood that the condition will occur. These incidence rates reflect historical spaceflight data or when appropriate, terrestrial data. In this presentation, we will explore how uncertainty in these rates propagate to the medical risk described by MEDPRAT. These results identify the conditions and parameters with the largest contribution to medical risks.

Ian Lim↗