Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Modeling”

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 739 records · Page 41

Determination of Turboprop Reduction Gearbox System Fatigue Life and Reliability

Two computational models to determine the fatigue life and reliability of a commercial turboprop gearbox are compared with each other and with field data. These models are (1) Monte Carlo simulation of randomly selected lives of individual bearings and gears comprising the system and (2) two-parameter Weibull distribution function for bearings and gears comprising the system using strict-series system reliability to combine the calculated individual component lives in the gearbox. The Monte Carlo simulation included the virtual testing of 744,450 gearboxes. Two sets of field data were obtained from 64 gearboxes that were first-run to removal for cause, were refurbished and placed back in service, and then were second-run until removal for cause. A series of equations were empirically developed from the Monte Carlo simulation to determine the statistical variation in predicted life and Weibull slope as a function of the number of gearboxes failed. The resultant L(sub 10) life from the field data was 5,627 hr. From strict-series system reliability, the predicted L(sub 10) life was 774 hr. From the Monte Carlo simulation, the median value for the L(sub 10) gearbox lives equaled 757 hr. Half of the gearbox L(sub 10) lives will be less than this value and the other half more. The resultant L(sub 10) life of the second-run (refurbished) gearboxes was 1,334 hr. The apparent load-life exponent p for the roller bearings is 5.2. Were the bearing lives to be recalculated with a load-life exponent p equal to 5.2, the predicted L(sub 10) life of the gearbox would be equal to the actual life obtained in the field. The component failure distribution of the gearbox from the Monte Carlo simulation was nearly identical to that using the strict-series system reliability analysis, proving the compatibility of these methods.

Zaretsky, Erwin V.↗

Development of an Uncertainty Model for the National Transonic Facility

This paper introduces an uncertainty model being developed for the National Transonic Facility (NTF). The model uses a Monte Carlo technique to propagate standard uncertainties of measured values through the NTF data reduction equations to calculate the combined uncertainties of the key aerodynamic force and moment coefficients and freestream properties. The uncertainty propagation approach to assessing data variability is compared with ongoing data quality assessment activities at the NTF, notably check standard testing using statistical process control (SPC) techniques. It is shown that the two approaches are complementary and both are necessary tools for data quality assessment and improvement activities. The SPC approach is the final arbiter of variability in a facility. Its result encompasses variation due to people, processes, test equipment, and test article. The uncertainty propagation approach is limited mainly to the data reduction process. However, it is useful because it helps to assess the causes of variability seen in the data and consequently provides a basis for improvement. For example, it is shown that Mach number random uncertainty is dominated by static pressure variation over most of the dynamic pressure range tested. However, the random uncertainty in the drag coefficient is generally dominated by axial and normal force uncertainty with much less contribution from freestream conditions.

Walter, Joel A.↗

Simulation of Hyperspectral Images

A software package generates simulated hyperspectral imagery for use in validating algorithms that generate estimates of Earth-surface spectral reflectance from hyperspectral images acquired by airborne and spaceborne instruments. This software is based on a direct simulation Monte Carlo approach for modeling three-dimensional atmospheric radiative transport, as well as reflections from surfaces characterized by spatially inhomogeneous bidirectional reflectance distribution functions. In this approach, "ground truth" is accurately known through input specification of surface and atmospheric properties, and it is practical to consider wide variations of these properties. The software can treat both land and ocean surfaces, as well as the effects of finite clouds with surface shadowing. The spectral/spatial data cubes computed by use of this software can serve both as a substitute for, and a supplement to, field validation data.

Richsmeier, Steven C.↗

Radiative Transfer Through Clouds and Its Applications in Support of the Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) Mission

The Greenland and Antarctic ice sheets, which contain enough ice to raise sea level by about 7 and 60 m, respectively, are losing mass at an increasing rate. To acquire continuous information of the cryosphere, after the Ice, Cloud, and land Elevation Satellite (ICESat) (2003-2010), NASA is actively planning for the ICESat-2 mission. Both ICESat and ICESat-2 are space-borne lidar altimetry systems. The systems measure the time of flight of the arriving photons that are reflected by the surface to deduce the elevation of the underlying terrain. As one of NASA's top priority missions, ICESat-2 is scheduled to launch in 2016. One of the major science goals of ICESat-2 is to quantify the ice sheet mass balance to determine its contributions to the sea level change and its impacts on ocean circulation (Abdalati et al. 2010). Compared to ICESat, which operates at 40 Hz and records the reflected laser energy as a waveform, the significantly improved ICESat-2 lidar employs a 532 nm micro-pulse photon counting system that operates at a high frequency of 10kHz with single photon detectability (Yang et al. 2012). To achieve its science goals, ICESat-2 requires the ability of detecting the elevation change with an accuracy of 0.2 cm/year over the entire ice sheet. Since every photon emitted by the lidar system will travel through the atmosphere, clouds can certainly affect the flight time of the arriving photons. Forward scattering by cloud particles increases the photon path length, thus resulting in biases in ice sheet elevation measurements known as atmospheric path delay (Duta et al. 2001, Yang et al. 2010, 2011). To ensure the accuracy of ICESat-2 surface elevation measurements, it is critical to understand how clouds would affect the travel time of arriving photons. In this talk, we will first present a framework that simulates the behavior of a space-borne 532 mn micro-pulse photon counting lidar in cloudy and clear atmospheres. To investigate the process of laser propagation through clouds, a 3-D Monte Carlo radiative transfer model is used to simulate the photon path distribution and the Poisson distribution is adopted for the number of photon returns. Since the photon counting system only registers the time of the first arriving photon within the detector "dead time", the retrieved average surface elevation tends to bias towards higher values. This is known as the first photon bias. With the scenarios simulated here, the first photon bias for clear sky is about 6.5 cm. Clouds affect surface altimetry in two ways: (1) cloud attenuation lowers the average number of arriving photons and hence reduces the first photon bias; (2) cloud forward scattering increases the photon path length and makes the surface appear further away from the satellite. Compared to clear sky, the average surface elevation detected by the photon counting system for cloudy sky with optical depth 1.0 is 4.0 to 6.0 cm lower for the simulations conducted. The effect of surface roughness on the accuracy of elevation retrievals will also discussed.

Yang, Yuekui↗

Polarimetric Imaging of Large Cavity Structures in the Pre-transitional Protoplanetary Disk around PDS 70: Observations of the Disk

We present high resolution H-band polarized intensity (PI; FWHM = 0."1: 14 AU) and L'-band imaging data (FWHM = 0."11: 15 AU) of the circumstellar disk around the weak-lined T Tauri star PDS 70 in Centaurus at a radial distance of 28 AU (0."2) up to 210 AU (1."5). In both images, a giant inner gap is clearly resolvro for the first time, and the radius of the gap is approx 70 AU. Our data show that the geometric center of the disk shifts by approx 6 AU toward the minor axis. We confirm that the brown dwarf companion candidate to the north of PDS 70 is a background star based on its proper motion. As a result of SED fitting by Monte Carlo radiative transfer modeling, we infer the existence of an optically thick inner disk at a few AU. Combining our observations and modeling, we classify the disk of PDS 70 as a pre-transitional disk. Furthermore, based on the analysis of L'-band imaging data, we put an upper limit mass of companions at approx 30 to approx 50M(sub J) within the gap. Taking account of the presence of the large and sharp gap, we suggest that the gap could be formed by dynamical interactions of sub-stellar companions or multiple unseen giant planets in the gap.

Hashimoto, J.↗

Limits on Large Extra Dimensions Based on Observations of Neutron Stars with the Fermi-LAT

We present limits for the compactification scale in the theory of Large Extra Dimensions (LED) proposed by Arkani-Hamed, Dimopoulos, and Dvali. We use 11 months of data from the Fermi Large Area Telescope (Fermi-LAT) to set gamma ray flux limits for 6 gamma-ray faint neutron stars (NS). To set limits on LED we use the model of Hannestad and Raffelt (HR) that calculates the Kaluza-Klein (KK) graviton production in supernova cores and the large fraction subsequently gravitationally bound around the resulting NS. The predicted decay of the bound KK gravitons to should contribute to the flux from NSs. Considering 2 to 7 extra dimensions of the same size in the context of the HR model, we use Monte Carlo techniques to calculate the expected differential flux of gamma-rays arising from these KK gravitons, including the effects of the age of the NS, graviton orbit, and absorption of gamma-rays in the magnetosphere of the NS. We compare our Monte Carlo-based differential flux to the experimental differential flux using maximum likelihood techniques to obtain our limits on LED. Our limits are more restrictive than past EGRET-based optimistic limits that do not include these important corrections. Additionally, our limits are more stringent than LHC based limits for 3 or fewer LED, and comparable for 4 LED. We conclude that if the effective Planck scale is around a TeV, then for 2 or 3 LED the compactification topology must be more complicated than a torus.

Ferrara, E. C.↗

Characterization of Aerodynamic Interactions with the Mars Science Laboratory Reaction Control System Using Computation and Experiment

On August 5 , 2012, the Mars Science Laboratory (MSL) entry capsule successfully entered Mars' atmosphere and landed the Curiosity rover in Gale Crater. The capsule used a reaction control system (RCS) consisting of four pairs of hydrazine thrusters to fly a guided entry. The RCS provided bank control to fly along a flight path commanded by an onboard computer and also damped unwanted rates due to atmospheric disturbances and any dynamic instabilities of the capsule. A preliminary assessment of the MSL's flight data from entry showed that the capsule flew much as predicted. This paper will describe how the MSL aerodynamics team used engineering analyses, computational codes and wind tunnel testing in concert to develop the RCS system and certify it for flight. Over the course of MSL's development, the RCS configuration underwent a number of design iterations to accommodate mechanical constraints, aeroheating concerns and excessive aero/RCS interactions. A brief overview of the MSL RCS configuration design evolution is provided. Then, a brief description is presented of how the computational predictions of RCS jet interactions were validated. The primary work to certify that the RCS interactions were acceptable for flight was centered on validating computational predictions at hypersonic speeds. A comparison of computational fluid dynamics (CFD) predictions to wind tunnel force and moment data gathered in the NASA Langley 31-Inch Mach 10 Tunnel was the lynch pin to validating the CFD codes used to predict aero/RCS interactions. Using the CFD predictions and experimental data, an interaction model was developed for Monte Carlo analyses using 6-degree-of-freedom trajectory simulation. The interaction model used in the flight simulation is presented.

Schoenenberger, Mark↗

Polarimetric Imaging of Large Cavity Structures in the Pre-transitional Protoplanetary Disk Around PDS 70: Observations of the Disk

We present high resolution H-band polarized intensity (PI; FWHM = 0."1: 14 AU) and L'-band imaging data (FWHM = 0."11: 15 AU) of the circumstellar disk around the weak-lined T Tauri star PDS 70 in Centaurus at a radial distance of 28 AU (0."2) up to 210 AU (1."5). In both images, a giant inner gap is clearly resolved for the first time, and the radius of the gap is approx.70 AU. Our data show that the geometric center of the disk shifts by approx.6 AU toward the minor axis. We confirm that the brown dwarf companion candidate to the north of PDS 70 is a background star based on its proper motion. As a result of SED fitting by Monte Carlo radiative transfer modeling, we infer the existence of an optically thick inner disk at a few AU. Combining our observations and modeling, we classify the disk of PDS 70 as a pre-transitional disk. Furthermore, based on the analysis of L'-band imaging data, we put an upper limit mass of companions at approx.30 to approx.50MJ within the gap. Taking account of the presence of the large and sharp gap, we suggest that the gap could be formed by dynamical interactions of sub-stellar companions or multiple unseen giant planets in the gap. Subject headings: planetary systems - protoplanetary disks - stars: individual (PDS 70) - stars: pre-main sequence - polarization

Hashimoto, J.↗

Polarimetric Imaging of Large Cavity Structures in the Pre-transitional Protoplanetary Disk Around PDS 70: Observations of the Disk

We present high-resolution H-band polarized intensity (FWHM=0".1:14AU) and L'-band imaging data(FWHM= 0".11:15 AU) of the circumstellar disk around the weak-lined T Tauri star PDS 70 in Centaurus at a radial distance of 28 AU (0".2) up to 210 AU (1".5). In both images, a giant inner gap is clearly resolved for the first time, and the radius of the gap is approx.70 AU. Our data show that the geometric center of the disk shifts by approx.6 AU toward the minor axis. We confirm that the brown dwarf companion candidate to the north of PDS 70 is a background star based on its proper motion. As a result of spectral energy distribution fitting by Monte Carlo radiative transfer modeling, we infer the existence of an optically thick inner disk at a few AU. Combining our observations and modeling, we classify the disk of PDS 70 as a pre-transitional disk. Furthermore, based on the analysis of L'-band imaging data, we put an upper limit of approx.30 to approx.50 M(sub J) on the mass of companions within the gap. Taking into account the presence of the large and sharp gap, we suggest that the gap could be formed by dynamical interactions of sub-stellar companions or multiple unseen giant planets in the gap. Key words: planetary systems - polarization - protoplanetary disks - stars: individual (PDS 70) - stars: pre-main sequence.

Hashimoto, J.↗

Anisotropic Solar Wind Sputtering of the Lunar Surface Induced by Crustal Magnetic Anomalies

The lunar exosphere is generated by several processes each of which generates neutral distributions with different spatial and temporal variability. Solar wind sputtering of the lunar surface is a major process for many regolith-derived species and typically generates neutral distributions with a cosine dependence on solar zenith angle. Complicating this picture are remanent crustal magnetic anomalies on the lunar surface, which decelerate and partially reflect the solar wind before it strikes the surface. We use Kaguya maps of solar wind reflection efficiencies, Lunar Prospector maps of crustal field strengths, and published neutral sputtering yields to calculate anisotropic solar wind sputtering maps. We feed these maps to a Monte Carlo neutral exospheric model to explore three-dimensional exospheric anisotropies and find that significant anisotropies should be present in the neutral exosphere depending on selenographic location and solar wind conditions. Better understanding of solar wind/crustal anomaly interactions could potentially improve our results.

exosphere↗

Multi-Model Ensemble Wake Vortex Prediction

Several multi-model ensemble methods are investigated for predicting wake vortex transport and decay. This study is a joint effort between National Aeronautics and Space Administration and Deutsches Zentrum fuer Luft- und Raumfahrt to develop a multi-model ensemble capability using their wake models. An overview of different multi-model ensemble methods and their feasibility for wake applications is presented. The methods include Reliability Ensemble Averaging, Bayesian Model Averaging, and Monte Carlo Simulations. The methodologies are evaluated using data from wake vortex field experiments.

Koerner, Stephan↗

Multi-Model Ensemble Wake Vortex Prediction

Several multi-model ensemble methods are investigated for predicting wake vortex transport and decay. This study is a joint effort between National Aeronautics and Space Administration and Deutsches Zentrum fuer Luft- und Raumfahrt to develop a multi-model ensemble capability using their wake models. An overview of different multi-model ensemble methods and their feasibility for wake applications is presented. The methods include Reliability Ensemble Averaging, Bayesian Model Averaging, and Monte Carlo Simulations. The methodologies are evaluated using data from wake vortex field experiments.

Koerner, Stephan↗

Multimodel Ensemble Methods for Prediction of Wake-Vortex Transport and Decay Originating NASA

Several multimodel ensemble methods are selected and further developed to improve the deterministic and probabilistic prediction skills of individual wake-vortex transport and decay models. The different multimodel ensemble methods are introduced, and their suitability for wake applications is demonstrated. The selected methods include direct ensemble averaging, Bayesian model averaging, and Monte Carlo simulation. The different methodologies are evaluated employing data from wake-vortex field measurement campaigns conducted in the United States and Germany.

Korner, Stephan↗

Full Inversion of Solar Relativistic Electron Events Measured by the Helios Spacecraft

The Parker Solar Probe and the incoming Solar Orbiter mission will provide measurements of solar energetic particle (SEP) events at close heliocentric distances from the Sun. Up to present, the largest data set of SEP events in the inner heliosphere are the observations by the two Helios spacecraft. Aims.We re-visit a sample of 15 solar relativistic electron events measured by the Helios mission with the goal of better characterizing the injection histories of solar energetic particles and their interplanetary transport conditions at heliocentric distances <1AU. Methods. The measurements provided by the E6 instrument on board Helios provide us with the electron directional distributionsin eight di erent sectors that we use to infer the detailed evolution of the electron pitch-angle distributions. The results of a Monte Carlo interplanetary transport model, combined with a full inversion procedure, were used to fit the observed directional intensities in the 300–800 keV nominal energy channel. Unlike previous studies, we have considered both the energy and angular responses of the detector. This method allowed us to infer the electron release time profile at the source and determine the electron interplanetary transport conditions. Results. We discuss the duration of the release time profiles and the values of the radial mean free path, and compare them with thevalues reported previously in the literature using earlier approaches. Five of the events show short injection histories (<30 min) at the Sun and ten events show long-lasting (>30 min) injections. The values of mean free path range from 0.02AU to 0.27AU.Conclusions. The inferred injection histories match with the radio and soft X-ray emissions found in literature.We find no dependence of the radial mean free path on the radial distance. In addition, we find no apparent relation between the strength of interplanetary scattering and the size of the solar particle release.

Pacheco, D.↗

Computational Bayesian Methods Applied to Complex Problems in Bio and Astro Statistics

In this dissertation we apply computational Bayesian methods to three distinct problems. In the first chapter, we address the issue of unrealistic covariance matrices used to estimate collision probabilities. We model covariance matrices with a Bayesian Normal-Inverse-Wishart model, which we fit with Gibbs sampling. In the second chapter, we are interested in determining the sample sizes necessary to achieve a particular interval width and establish non-inferiority in the analysis of prevalences using two fallible tests. To this end, we use a third order asymptotic approximation. In the third chapter, we wish to synthesize evidence across multiple domains in measurements taken longitudinally across time, featuring a substantial amount of structurally missing data, and fit the model with Hamiltonian Monte Carlo in a simulation to analyze how estimates of a parameter of interest change across sample sizes.

Elrod, Chris↗

Optical Performance Prediction of the Thirty Meter Telescope after Initial Alignment Using Optical Modeling

We present an estimate of the optical performance of the Thirty Meter Telescope (TMT) after execution of the full telescope alignment plan. The TMT alignment is performed by the Global Metrology System (GMS) and the Alignment and Phasing System (APS). The GMS first measures the locations of the telescope optics and instruments as a function of elevation angle. These initial measurements will be used to adjust the optics positions and build initial elevation look-up tables. Then the telescope is aligned using starlight as the input for the APS at multiple elevation angles. APS measurements are used to refine the telescope alignment to build elevation and temperature dependent look-up tables. Due to the number of degrees of freedom in the telescope (over 10,000), the ability of the primary mirror to correct aberrations on other optics, the tight optical performance requirements and the multiple instrument locations, it is challenging to develop, test and validate these alignment procedures. In this paper, we consider several GMS and APS operational scenarios. We apply the alignment procedures to the model-generated TMT, which consists of various quasi-static errors such as polishing errors, passive supports errors, thermal and gravity deformations and installation position errors. Using an integrated optical model and Monte-Carlo framework, we evaluate the TMT’s aligned states using optical performance metrics at multiple instrument and field of view locations. The optical performance metrics include the Normalized Point Source Sensitivity (PSSN), RMS wavefront error before and after Adaptive Optics (AO) correction, pupil position change, and plate scale distortion.

Rogers, John↗

Development of the Ames Global Hyperspectral Synthetic Data Set: Surface Bidirectional Reflectance Distribution Function

This study introduces the Ames Global Hyperspectral Synthetic Data set (AGHSD), in particular the surface bidirectional reflectance distribution function (BRDF) product, to support the NASA Surface Biology and Geology (SBG) mission development. The data set is generated based on the corresponding multispectral BRDF products from NASA's MODIS satellite sensor. Based on theories of radiative transfer in vegetation canopies, we derive a simple but robust relationship that indicates that the hyperspectral surface BRDF can be accurately approximated as a weighted sum of the soil surface reflectance, the leaf single albedo, and the canopy scattering coefficient, where the weights or coefficients are spectrally invariant and thus readily estimated from the multispectral MODIS products. We validate the algorithm with simulations by a Monte Carlo Ray Tracing model and find the results highly consistent with the theoretic derivation. Using reflectance spectra of soil and vegetation derived from existing spectral libraries, we apply the algorithm to generate the AGHSD BRDF product at 1 km and 8-day resolutions for the year of 2019. The data set is biogeochemically and biogeophysically coherent and consistent, and serves the goal to support the SBG community in developing sciences and applications for the future global imaging spectroscopy mission.

Hyperspectral Remote Sensing↗