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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Extended atmospheres of comets and outer planet-satellite systems

For the hydrogen coma of comet P/Halley, both a Lyman-alpha image and extensive Lyman-alpha scan data obtained by the Pioneer Venus Orbiter Ultraviolet spectrometer as well as H-alpha ground-based spectral observations obtained by the University of Wisconsin Space Physics Group were successfully interpreted and analyzed with our Monte Carlo particle trajectory model. The excellent fit of the model and the Halley data and the water production rate determined near perihelion (9 Feb. 1986) from 13 Dec. 1985 to 13 Jan. 1986 and from 1 Feb. to 7 Mar. 1986 are discussed. Studies for the circumplanetary distribution of atomic hydrogen in the Saturn and Neptune systems were undertaken for escape of H atoms from Titan and Triton, respectively. The discovery of a new mechanism which can dramatically change the normal cylindrically symmetric distribution of hydrogen about the planet is discussed. The implications for the Titan-Saturn and Triton-Neptune are summarized.

Smyth, William H.↗
NASA NTRS1993

On the evaluation of air mass factors for atmospheric near-ultraviolet and visible absorption spectroscopy

The interpretation of UV-visible twilight absorption measurements of atmospheric chemical constituents is dependent on how well the optical path, or air mass factor, of light collected by the spectrometer is understood. A simple single scattering model and a Monte Carlo radiative transfer scheme have been developed to study the effects of multiple scattering, aerosol scattering, surface albedo and refraction on air mass factors for scattered light observations. At fairly short visible wavelengths (less than about 450 nm), stratospheric air mass factors are found to be relatively insensitive to multiple scattering, surface albedo and refraction, as well as aerosol scattering by background aerosols. Longer wavelengths display greater sensitivity to refraction and aerosol scattering. Tropospheric air mass factors are found to be highly dependent on aerosol scattering, surface albedo and, at long visible wavelengths (about 650 nm), refraction. Absorption measurements of NO2 and O4 are shown to support these conclusions.

Perliski, Lori M.↗
NASA NTRS1993

Analysis of hydrogen H-alpha observations of the coma of Comet P/Halley

The Monte Carlo Particle Trajectory Model of Combi and Smyth (1988) is used here to analyze observations of the H-alpha coma of Comet Halley. The solar excitation mechanism for the H-alpha emissions line is described. The H2O production rates derived for the H-alpha brightness measurements are shown to be very consistent with the H2O production rates determined from other Comet Halley observations of the H, O, and OH comae. Revised H2O production rates determined from 6300 A brightness measurements are presented.

Smyth, William H.↗
NASA NTRS1993

PDF approach for turbulent scalar field: Some recent developments

The probability density function (PDF) method has been proven a very useful approach in turbulence research. It has been particularly effective in simulating turbulent reacting flows and in studying some detailed statistical properties generated by a turbulent field There are, however, some important questions that have yet to be answered in PDF studies. Our efforts in the past year have been focused on two areas. First, a simple mixing model suitable for Monte Carlo simulations has been developed based on the mapping closure. Secondly, the mechanism of turbulent transport has been analyzed in order to understand the recently observed abnormal PDF's of turbulent temperature fields generated by linear heat sources.

Gao, Feng↗
NASA NTRS1995

Observations and analysis of O(1D) and NH2 line profiles for the coma of comet P/Halley

A set of high-resolution Fabry-Perot measurements of the coma of comet P/Halley was acquired in the (O I) 6300 A and NH2 6298.62 A emission lines. These high-resolution measurements provide the first optical observations capable of studying directly the photochemical kinetics and dynamic outflow of the coma. The observations were analyzed by a Monte Carlo Particle Trajectory Model. The agreement of the model and observed line profiles was excellent and verified the underlying dynamics, exothermic photodissociative chemistry, and collisional thermalization in the coma. The somewhat wider intrinsic line profile width for the O(1D) emission in 1986 January compared to 1986 May, is, for example, produced by the larger outflow speeds and gas temperatures nearer perihelion in January. The January O(1D) profile, which is wider than the January NH2 profile, is indicative of the photochemical kinetics in the dissociation of the parent molecules H2O and OH in the coma. The absolute calibration of the observations in 1986 January allowed the production rates for H2O and the NH2-parent molecules to be determined. The average daily water production rates derived from the O(1D) emission data for January 16 and 17 are presented. These very large water production rates are consistent with the extrapolated (and 7.6 day time variable) water production rates determined from the analysis of lower spectral resolution observations for O(1D) and H-alpha emissions that covered the time period up to January 13. The large production rates on January 16 and 17 establish that the maximum water production rate for comet Halley accurred pre-perihelion in January. Implications drawn from comparison with 18 cm radio emission data in January suggest that the peak water production rate was even larger. The average production rate for NH3 determined from the NH2 emission data for January 17 was (1.48 +/- 0.10) x 10(exp 28) molecules/s, yielding an NH3/H2O production rate ratio of 0.55%.

Smyth, William H.↗
NASA NTRS1995

Analysis of hydrogen Lyman-alpha observations of the coma of Comet P/Halley near the perihelion

The pioneer Venus Orbiter Ultraviolet Spectrometer measurements of the Lyman-alpha intensity of atomic hydrogen excited by solar resonance scattering in the coma of Comet P/Halley acquired from December 28, 1985, to January 6, 1986, and from January 31, 1986, to March 6, 1986, are simulated with the Monte Carlo Particle Trajectory Model corrected for optical depth effects. Spatially detailed comparisons between data and model show excellent agreement and are used to infer that the highest cometary activity may not be at perihelion, but about 2 1/2 weeks before. An improved set of H2O production rates is presented for the period of time that the spectrometer was observing and is found to be consistent with the rates from other types of measurements. The apparent discrepancy between Stewart (1987) in early March and International Ultraviolet Explorer OH derived rates is resolved. The problem with the conversion of 18-cm OH radio brightness to H2O production rates is also discussed.

Smyth, William H.↗
NASA NTRS1998

Microwave Brightness Temperatures of Tilted Convective Systems

Aircraft and ground-based radar data from the Tropical Ocean and Global Atmosphere Coupled-Ocean Atmosphere Response Experiment (TOGA COARE) show that convective systems are not always vertical. Instead, many are tilted from vertical. Satellite passive microwave radiometers observe the atmosphere at a viewing angle. For example, the Special Sensor Microwave/Imager (SSM/I) on Defense Meteorological Satellite Program (DMSP) satellites and the Tropical Rainfall Measurement Mission (TRMM) Microwave Imager (TMI) on the TRMM satellite have an incident angle of about 50deg. Thus, the brightness temperature measured from one direction of tilt may be different than that viewed from the opposite direction due to the different optical depth. This paper presents the investigation of passive microwave brightness temperatures of tilted convective systems. To account for the effect of tilt, a 3-D backward Monte Carlo radiative transfer model has been applied to a simple tilted cloud model and a dynamically evolving cloud model to derive the brightness temperature. The radiative transfer results indicate that brightness temperature varies when the viewing angle changes because of the different optical depth. The tilt increases the displacements between high 19 GHz brightness temperature (Tb(sub 19)) due to liquid emission from lower level of cloud and the low 85 GHz brightness temperature (Tb(sub 85)) due to ice scattering from upper level of cloud. As the resolution degrades, the difference of brightness temperature due to the change of viewing angle decreases dramatically. The dislocation between Tb(sub 19) and Tb(sub 85), however, remains prominent.

Hong, Ye↗
NASA NTRS1998

Assay of the Martian Regolith with Neutrons

Different aspects of assaying Martian regolith using neutrons have been investigated. The epithermal portion of moderated neutrons spectra is dramatically effected by the presence of hydrogen (usually in the form of water). A simple analytic formula has been derived to describe the amplitude of this portion of the neutron spectrum as a function of water concentration. Several demonstration experiments have been performed and modeled with a Monte Carlo code. Results of these experiments generally agreed with the calculations to within 20%. In addition to He-3 detectors, lithium-glass scintillators and U-238 fission ion chambers were investigated to determine their applicability to space experiments.

Drake, Darrell M.↗
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.↗