Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Modelling”

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

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↗

Uncertainty quantification in elastic constants of SiC f /SiC m tubular composites using global sensitivity analysis

Silicon carbide fiber and silicon carbide matrix (SiC f /SiC m ) tubes produced through the chemical vapor infiltration process have become a candidate cladding material in nuclear applications. The performance of this composite is influenced by many variables such as braiding angle, porosity, material properties, etc., which vary over a range of values due to the inherent fluctuations in the manufacturing process. In this study, the variability in elastic constants of SiC f /SiC m composite has been quantified through multiscale finite element (FE) simulations, variable screening, and high-fidelity surrogate modeling. The key variables dominantly affecting the elastic constants of SiC f /SiC m tubes were identified using global sensitivity analysis. A surrogate to the high-fidelity FE-based model was used in Monte Carlo simulations to generate a hundred thousand samples from which the uncertainty in elastic constants was assessed. It turned out that the coefficient of variation was less than 10%.

Materials Science↗

Construction of Supersonic Flight Dynamics Test (SFDT-2) Vehicle Monte Carlo Splashdown Footprints for Use in Range Safety and Recovery Planning

Overview of Low Density Supersonic Decelerator (LDSD) SFDT-2 ( Supersonic Flight Dynamics Test-2) flight profile; Description of Monte Carlo flight dynamics simulations performed; Modeling and uncertainties; Results with statistical prediction of test vehicle splashdown location; Use of Monte Carlo results for range safety and recovery planning.

Strauss, Bill↗