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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

An Analysis of Water Line Profiles in Star Formation Regions Observed by SWAS

We present spectral line profiles for the 557 GHz 1(sub 1,0) yields 1(sub 0,1) ground-state rotational transition of ortho-H2(16)O for 18 galactic star formation regions observed by SWAS. 2 Water is unambiguously detected in every source. The line profiles exhibit a wide variety of shapes, including single-peaked spectra and self-reversed profiles. We interpret these profiles using a Monte Carlo code to model the radiative transport. The observed variations in the line profiles can be explained by variations in the relative strengths of the bulk flow and small-scale turbulent motions within the clouds. Bulk flow (infall, outflow) must be present in some cloud cores, and in certain cases this bulk flow dominates the turbulent motions.

Ashby, Matthew L. N.↗
NASA NTRS2001

Transport Phenomena of Off-Axis Sputtering Deposition

Various high quality epitaxial films, especially oxides, have been synthesized using off-axis sputtering deposition. In this presentation, we report the experiment results of ZnO films grown by the off-axis sputtering deposition. Films were synthesized in temperatures ranged from room temperature to 600 C, and pressures from 5 mTorr to 150 mTorr. Film growth rate was measured by surface profilometer, ellipsometer, and wavelength dispersive spectrometry. Due to the collisions between the sputtered species and the residue gases, the kinetic energy of species was reduced and the transport of depositing species changed from a ballistic movement for low pressure to a diffusive drift for high pressure in which the transport species were almost thermalized. The measurements show an increase of growth rates along the gravity vector when the Knudsen number of transport species is less than 0.05, which suggests that gravity affected the transport characterization in off-axis sputtering deposition. Because the product of pressure (p) and travel distance (d) of sputtered species, pd, exceeds several mTorr-cm during film deposition, the classical simulations for sputtering process in high vacuum system may not be applied. Based on these experimental measurements, a transport process of the off-axis sputtering deposition is proposed. Several models including the Monte Carlo method and gravity-driven flow dynamics simulation will be discussed.

Zhu, S.↗
NASA NTRS2002

The Evolution of Disks and Winds in Dwarf Nova Outbursts - FUSE

This project was a project to study the FUV spectra of two proto-typical dwarf novae, U Gem and SS Cygni, through an outburst cycle. The luminosity of the boundary layer in the two systems, as evidenced by earlier EUVE observations, is different in the two systems. Our intensive study of the two systems was intended to (1) probe the ionization and kinematic structure of the wind as a function of system brightness, (2) isolate the contributions of the disk to the FUV spectra, and (3) examine physical conditions and abundances of material just being accreted onto the disk from the secondary. The U Gem and SS Cyg observations took place in March and October 2000, respectively. The data obtained with FUSE was of excellent quality. Analysis of the both observations is now essentially complete, although some modeling of the SS Cyg spectra is ongoing, as we complete an ApJ manuscript on this object. Our main results for U Gem are as follows: The plateau spectra have continuum shapes and fluxes that are approximated by steady state accretion disk model spectra with an accretion rate 7x10(exp 9) Msolar/yr. The spectra also show numerous absorption lines of H I, He II, and 2-5 times ionized transitions of C, N, O, P, S, and Si. There are no emission features in the spectra, with the possible exception of a weak feature on the red wing of the 0 VI doublet. The absorption lines are narrow (FWHM approx. 50 km/s), too narrow to arise from the disk photosphere, and at low velocities (less than or equal to 700 km/s). The S VI and O VI doublets are optically thick. The absorption lines in the plateau spectra show orbital variability: in spectra obtained at orbital phases between 0.53 and 0.79, low-ionization absorption lines appear and the central depths of the preexisting lines increase. The increase in line absorption occurs at the same orbital phases as previously observed EUV and X-ray light-curve dips. If the absorbing material is in (near-) Keplerian rotation around the disk, it must be located at large disk radii. The final observation occurred when U Gem was about 2 mag from optical quiescence. The spectra are dominated by emission from an approx. 43,000 K, metal-enriched white dwarf (WD). The inferred radius of the WD is 4.95x10(exp 8) cm, close to that observed in quiescence. Allowing for a hot heated region on the surface of the WD improves the fit to the spectrum at short (less than 960 A) wavelengths. Our main results for SS Cyg are as follows: The first two of four observations of SS Cyg show disk dominated spectra with accretion rates of order 10(exp -8) Msolar/yr. Except for narrow interstellar features (atomic and molecular H), the lines are all broad consistent with a disk or wind origin. The O VI line in the spectra is mostly of wind origin as detailed modeling with our Monte Carlo code (developed in part using funds from this project) show. The continua from spectra in observations 3 and 4, observed during the decline phase, are not well fit with steady-state disks, and show considerable resemblance to quiescent spectra obtained with HUT. The most probable interpretation for the emission features seen in the spectrum in the last two observations is that they arise from a photo-illuminated choronosphere above the disk, rather than a wind.

Long, Knox↗
NASA NTRS2002

Software for Simulation of Hyperspectral Images

A package of software generates simulated hyperspectral images 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 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 and 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.

Richtsmeier, Steven C.↗
NASA NTRS2004

Designing to Mitigate Food Growing Failures in Space

Future space life support systems may use crop plants to grow most of the crew s food. A harvest failure can reduce the food available for future consumption. If the previously stored food is insufficient to reach the next harvest, the crew may go hungry. This paper considers how the overall food supply system should be modified to cope with food production failures. The food supply concept for a mission will use grown food, or stored food, cIr both. The optimum food supply mix depends on the costs and failure probabilities of stored and grown food. A simple food system model assumes that either we obtain the nominal harvest or a failure occurs and no food is harvested. Given the probability that any particular harvest fails, it is easy to compute the expected number of failures and the total food shortfall over a mission. If some food is grown and the probability of harvest failure is high, a non-redundant system has an unacceptable likelihood that the crew will have no food for a full harvest period. Food supply reliability must be increased either by supplying more food initially or by increasing food production capacity. We can obtain a very reliable food supply even when the harvest failure rate is high. If the cost of growing food is much less than the cost of providing stored food, it is better to provide redundant food growing capacity than to increase initial storage. A more realistic biomass production failure model allows the harvest amount or time to vary around the nominal values, using stochastic modeling with repeated Monte Carlo simulation, but such failures have minor impact compared to a complete harvest failure.

Jones, Harry↗
NASA NTRS1998

A Novel Solution-Technique Applied to a Novel WAAS Architecture

The Federal Aviation Administration has embarked on an historic task of modernizing and significantly improving the national air transportation system. One system that uses the Global Positioning System (GPS) to determine aircraft navigational information is called the Wide Area Augmentation System (WAAS). This paper describes a reliability assessment of one candidate system architecture for the WAAS. A unique aspect of this study regards the modeling and solution of a candidate system that allows a novel cold sparing scheme. The cold spare is a WAAS communications satellite that is fabricated and launched after a predetermined number of orbiting satellite failures have occurred and after some stochastic fabrication time transpires. Because these satellites are complex systems with redundant components, they exhibit an increasing failure rate with a Weibull time to failure distribution. Moreover, the cold spare satellite build-time is Weibull and upon launch is considered to be a good-as-new system with an increasing failure rate and a Weibull time to failure distribution as well. The reliability model for this system is non-Markovian because three distinct system clocks are required: the time to failure of the orbiting satellites, the build time for the cold spare, and the time to failure for the launched spare satellite. A powerful dynamic fault tree modeling notation and Monte Carlo simulation technique with importance sampling are shown to arrive at a reliability prediction for a 10 year mission.

Bavuso, J.↗
NASA NTRS2005

Instrument Drift Uncertainties and the Long-Term TOMS/SBUV Total Ozone Record

Long-term climate records from satellites are often constructed from the measurements of a sequence of instruments launched at different times. Each of these instruments is calibrated prior to launch. After launch they are subjected to potential offsets and slow drifts in calibration. We illustrate these issues in the construction of a merged total ozone record from two TOMS and three SBUV instruments. This record extends from late 1978 through the present. The question is "How good are these records?". We have examined the uncertainty in determining the relative calibration of two instruments during an overlap period in their measurements. When comparing a TOMS instrument, such as that on Nimbus 7, with an SBUV instrument, also on Nimbus 7, we find systematic differences and random differences. We have combined these findings with estimates of individual instrument drift into a monte- carlo uncertainty propagation model. We estimate an instrument drift uncertainty of a little larger than 1 percent per decade over the 25-year history of the TOMS/SBUV measurements. We make an independent estimate of the drift uncertainty in the ground-based network of total ozone measurements and find it to be of similar, but slightly smaller magnitude. The implications of these uncertainties for trend and recovery determination will be discussed.

Solarski, Richard S.↗
NASA NTRS2006

The Effect of Cumulus Cloud Field Anisotropy on Domain-Averaged Solar Fluxes and Atmospheric Heating Rates

Cumulus clouds can become tilted or elongated in the presence of wind shear. Nevertheless, most studies of the interaction of cumulus clouds and radiation have assumed these clouds to be isotropic. This paper describes an investigation of the effect of fair-weather cumulus cloud field anisotropy on domain-averaged solar fluxes and atmospheric heating rate profiles. A stochastic field generation algorithm was used to produce twenty three-dimensional liquid water content fields based on the statistical properties of cloud scenes from a large eddy simulation. Progressively greater degrees of x-z plane tilting and horizontal stretching were imposed on each of these scenes, so that an ensemble of scenes was produced for each level of distortion. The resulting scenes were used as input to a three-dimensional Monte Carlo radiative transfer model. Domain-average transmission, reflection, and absorption of broadband solar radiation were computed for each scene along with the average heating rate profile. Both tilt and horizontal stretching were found to significantly affect calculated fluxes, with the amount and sign of flux differences depending strongly on sun position relative to cloud distortion geometry. The mechanisms by which anisotropy interacts with solar fluxes were investigated by comparisons to independent pixel approximation and tilted independent pixel approximation computations for the same scenes. Cumulus anisotropy was found to most strongly impact solar radiative transfer by changing the effective cloud fraction, i.e., the cloud fraction when the field is projected on a surface perpendicular to the direction of the incident solar beam.

Hinkelman, Laura M.↗
NASA NTRS2007

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.↗
NASA NTRS2010

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.↗
NASA NTRS2004

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.↗
NASA NTRS2012

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↗
NASA NTRS2012

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.↗
NASA NTRS2012

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.↗
NASA NTRS2013

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↗
NASA NTRS2012

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.↗
NASA NTRS2012

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.↗