Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct simulation Monte Carlo”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Mean sea surface height of the Antarctic Circumpolar Current from Geosat data: Method and application

The mean sea surface height across the Antarctic Circumpolar Current (ACC) has been reconstructed from height variability measured by the Geosat altimeter without assuming prior knowledge of the geoid. For this study, an automated technique has been developed to estimate mean sea surface height for each satellite ground track using a meandering Gaussian jet model, and errors have been estimated using Monte Carlo simulation. The results are objectively mapped to produce a picture of the mean Subantarctic and Polar Fronts, which together comprise the major components of the Antarctic Circumpolar Current. The meandering jet model explains between 40% and 70% of the height variance along the jet axes. The results show that the fronts are substantially steered by topography and that the jets have an average Gaussian width of about 44 km in the meridional direction and meander about 75 km to either side of their mean locations. The average height difference across the Subantarctic Front (SAF) is 0.7 m and across the Polar Front (PF) 0.6 m. The mean widths of the fronts are correlated with the size of the baroclinic Rossby radius.

Gille, Sarah T.↗
NASA NTRS2000

Analytic Guided-Search Model of Human Performance Accuracy in Target- Localization Search Tasks

Current models of human visual search have extended the traditional serial/parallel search dichotomy. Two successful models for predicting human visual search are the Guided Search model and the Signal Detection Theory model. Although these models are inherently different, it has been difficult to compare them because the Guided Search model is designed to predict response time, while Signal Detection Theory models are designed to predict performance accuracy. Moreover, current implementations of the Guided Search model require the use of Monte-Carlo simulations, a method that makes fitting the model's performance quantitatively to human data more computationally time consuming. We have extended the Guided Search model to predict human accuracy in target-localization search tasks. We have also developed analytic expressions that simplify simulation of the model to the evaluation of a small set of equations using only three free parameters. This new implementation and extension of the Guided Search model will enable direct quantitative comparisons with human performance in target-localization search experiments and with the predictions of Signal Detection Theory and other search accuracy models.

Eckstein, Miguel P.↗
NASA NTRS1973

Decision-directed detector for overlapping PCM/NRZ signals.

A decision-directed (DD) technique for the detection of overlapping PCM/NRZ signals in the presence of white Gaussian noise is investigated. The performance of the DD detector is represented by probability of error Pe versus input signal-to-noise ratio (SNR). To examine how much improvement in performance can be achieved with this technique, Pe's with and without DD feedback are evaluated in parallel. Further, analytical results are compared with those found by Monte Carlo simulations. The results are in good agreement.

Wang, C. D.↗
NASA NTRS1981

Monte Carlo simulation of charged particles upstream of the earth's bow shock

A simple model of the earth's bow shock is presented and compared with ISEE-1 observations of diffuse ions upstream of the shock. Ion behavior is modeled in a test particle simulation using Monte Carlo techniques for a plane parallel shock at which cold solar wind ions are heated sufficiently to allow diffusion back upstream. Fitting of the model to the observations indicates that the diffusion coefficient must increase faster than linearly with velocity, and provides an upper limit to the acceleration efficiency defined as the proportion of incoming solar wind ions accelerated and observed upstream. Results suggest that astrophysical shocks are capable of accelerating ions directly from the thermal plasma without the need for a distinct preenergizing injection stage.

Ellison, D. C.↗
NASA NTRS1982

Monte Carlo studies of ocean wind vector measurements by SCATT: Objective criteria and maximum likelihood estimates for removal of aliases, and effects of cell size on accuracy of vector winds

The scatterometer on the National Oceanic Satellite System (NOSS) is studied by means of Monte Carlo techniques so as to determine the effect of two additional antennas for alias (or ambiguity) removal by means of an objective criteria technique and a normalized maximum likelihood estimator. Cells nominally 10 km by 10 km, 10 km by 50 km, and 50 km by 50 km are simulated for winds of 4, 8, 12 and 24 m/s and incidence angles of 29, 39, 47, and 53.5 deg for 15 deg changes in direction. The normalized maximum likelihood estimate (MLE) is correct a large part of the time, but the objective criterion technique is recommended as a reserve, and more quickly computed, procedure. Both methods for alias removal depend on the differences in the present model function at upwind and downwind. For 10 km by 10 km cells, it is found that the MLE method introduces a correlation between wind speed errors and aspect angle (wind direction) errors that can be as high as 0.8 or 0.9 and that the wind direction errors are unacceptably large, compared to those obtained for the SASS for similar assumptions.

Pierson, W. J.↗
NASA NTRS2007

Extinction-to-Backscatter Ratios of Lofted Aerosol Layers Observed During the First Three Months of CALIPSO Measurements

Case studies from the first three months of the Cloud and Aerosol Lidar and Infrared Pathfinder Spaceborne Observations (CALIPSO) measurements of lofted aerosol layers are analyzed using transmittance [Young, 1995] and two-wavelength algorithms [Vaughan et al., 2004] to determine the aerosol extinction-to-backscatter ratios at 532 and 1064 nm. The transmittance method requires clear air below the layer so that the transmittance through the layer can be determined. Suitable scenes are selected from the browse images and clear air below features is identified by low 532 nm backscatter signal and confirmed by low depolarization and color ratios. The transmittance and two-wavelength techniques are applied to a number of lofted layers and the extinction-to-backscatter ratios are compared with values obtained from the CALIPSO aerosol models [Omar et al., 2004]. The results obtained from these studies are used to adjust the aerosol models and develop observations based extinction-to-backscatter ratio look-up tables and phase functions. Values obtained by these techniques are compared to Sa determinations using other independent methods with a goal of developing probability distribution functions of aerosol type-specific extinction to backscatter ratios. In particular, the results are compared to values determined directly by the High Spectral Resolution Lidar (HSRL) during the CALIPSO CloudSat Validation Experiments (CCVEX) and Sa determined by the application of the two-wavelength lidar Constrained Ratio Aerosol Model-fit (CRAM) retrieval approach [Cattrall et al., 2005; Reagan et al., 2004] to the HSRL data. The results are also compared to values derived using the empirical relationship between the multiple-scattering fraction and the linear depolarization ratio by using Monte Carlo simulations of water clouds [Hu et al., 2006].

Omar, Ali H.↗
NASA NTRS2010

Absorbed Dose and Dose Equivalent Calculations for Modeling Effective Dose

While in orbit, Astronauts are exposed to a much higher dose of ionizing radiation than when on the ground. It is important to model how shielding designs on spacecraft reduce radiation effective dose pre-flight, and determine whether or not a danger to humans is presented. However, in order to calculate effective dose, dose equivalent calculations are needed. Dose equivalent takes into account an absorbed dose of radiation and the biological effectiveness of ionizing radiation. This is important in preventing long-term, stochastic radiation effects in humans spending time in space. Monte carlo simulations run with the particle transport code FLUKA, give absorbed and equivalent dose data for relevant shielding. The shielding geometry used in the dose calculations is a layered slab design, consisting of aluminum, polyethylene, and water. Water is used to simulate the soft tissues that compose the human body. The results obtained will provide information on how the shielding performs with many thicknesses of each material in the slab. This allows them to be directly applicable to modern spacecraft shielding geometries.

Welton, Andrew↗
NASA NTRS1993

The cosmic Mach number - Direct comparisons of observations and models

A robust method is used to derive the cosmic Mach number, M, from existing peculiar velocity data sets. For the IR Tully-Fisher distances of spirals in the Local Supercluster of Aaronson et al. (1979, 1982), M of 1.03 and a characteristic distance from the Local Group of 1464 km/s are found. For the elliptical galaxy sample of Faber et al. (1987, 1988), M of 0.57, with a characteristic distance of 2572 km/s is found. These results are compared with Monte Carlo simulations of the observational realizations drawn from numerical simulations of the universe based on various scenarios, including the standard cold dark matter (CDM) scenario. The effect of velocity bias on the derived Mach number is found to be small. Only 5 percent of the CDM simulations have M as large as or larger than that derived for the Aaronson et al. sample.

Strauss, Michael A.↗
NASA NTRS1971

Study of automatic flare and decrab guidance and control system for the space shuttle

The automatic landing of an unpowered space shuttle vehicle was studied along with the design of the guidance and flight controls, selection of the ground reference component of the landing system (signals in space), and evaluation of the performance of the overall system. The study was carried out by means of an all-digital simulation which incorporated wind shears, gusts and error models for the hardware components. The performance of the landing system operating in winds and turbulence and with hardware errors was determined by Monte Carlo techniques. It was found that the interactions of turbulence and vehicle characteristics were the limiting factors on system performance. The touchdown errors due to the measurement errors of the ground reference system were less significant. The guidance designed was predictive in the longitudinal direction and fixed path in the lateral direction.

Cockayne, W.↗
NASA NTRS2014

Solar Wind Implantation into Lunar Regolith: Hydrogen Retention in a Surface with Defects

Solar wind protons are implanted directly into the top 100 nm of the lunar near-surface region, but can either quickly diffuse out of the surface or be retained, depending upon surface temperature and the activation energy, U, associated with the implantation site. In this work, we explore the distribution of activation energies upon implantation and the associated hydrogen-retention times; this for comparison with recent observation of OH on the lunar surface. We apply a Monte Carlo approach: for simulated solar wind protons at a given local time, we assume a distribution of U values with a central peak, U(sub c) and width, U(sub w), and derive the fraction retained for long periods in the near-surface. We find that surfaces characterized by a distribution with predominantly large values of U (greater than 1 eV) like that expected at defect sites will retain implanted H (to likely form OH). Surfaces with the distribution predominantly at small values of U (less than 0.2 eV) will quickly diffuse away implanted H. However, surfaces with a large portion of activation energies between 0.3 eV less than U less than 0.9 eV will tend to be H-retentive in cool conditions but transform into H-emissive surfaces when warmed (as when the surface rotates into local noon). These mid-range activation energies give rise to a diurnal effect with diffusive loss of H at noontime.

regolith↗
NASA NTRS1985

Shock acceleration of diffuse ions at the earth's bow shock Acceleration efficiency and A/Z enhancement

Observations of particle spectra, intensity, and enhancement of alpha particles over protons at diffuse ion events at the quasi-parallel earth bow shock are compared to a Monte Carlo simulation of diffusive shock acceleration. The simulation includes the back reaction of accelerated particles on the shock structure, particle escape at an upstream free escape boundary, and a low energy per nucleon threshold for thermal leakage of downstream, shock-heated particles into the upstream region. The simulation assumes that the same scattering operator that gives rise to shock acceleration can also describe a viscous shock governed by hydrodynamic turbulence. This implies that accelerated ions can be drawn directly from the thermal solar wind with no separate superthermal seed population. Good agreement between the simulation and observations made during nearly radial magnetic field configurations lends support to thermal leakage of downstream, shock-heated ions as the mode of injection for diffusion ion events.

Ellison, D. C.↗
NASA NTRS2011

Calcium in Mercury's Exosphere: Modeling MESSENGER Data

Mercury is surrounded by a surface-bounded exosphere comprised of atomic species including hydrogen, sodium, potassium, calcium, magnesium, and likely oxygen. Because it is collisionless. the exosphere's composition represents a balance of the active source and loss processes. The Mercury Atmospheric and Surface Composition Spectrometer (MASCS) on the MErcury Surface. Space ENvironment. GEochemistry. and Ranging (MESSENGER) spacecraft has made high spatial-resolution observations of sodium, calcium, and magnesium near Mercury's surface and in the extended, anti-sunward direction. The most striking feature of these data has been the substantial differences in the spatial distribution of each species, Our modeling demonstrates that these differences cannot be due to post-ejection dynamics such as differences in photo-ionization rate and radiation pressure. but instead point to differences in the source mechanisms and regions on the surface from which each is ejected. The observations of calcium have revealed a strong dawn/dusk asymmetry. with the abundance over the dawn hemisphere significantly greater than over the dusk. To understand this asymmetry, we use a Monte Carlo model of Mercury's exosphere that we developed to track the motions of exospheric neutrals under the influence of gravity and radiation pressure. Ca atoms can be ejected directly from the surface or produced in a molecular exosphere (e.g., one consisting of CaO). Particles are removed from the system if they stick to the surface or escape from the model region of interest (within 15 Mercury radii). Photoionization reduces the final weighting given to each particle when simulating the Ca radiance. Preliminary results suggest a high temperature ( I-2x 10(exp 4) K) source of atomic Ca concentrated over the dawn hemisphere. The high temperature is consistent with the dissociation of CaO in a near-surface exosphere with scale height <= 100 km, which imparts 2 eV to the freshly produced Ca atom. This source region and energy are consistent with data from the three MESSENGER flybys; whether this holds true for the data obtained in orbit is under investigation.

Burger, Matthew H.↗
NASA NTRS1998

Calculations of multiple scattering-induced errors in the GLAS mission

The Geoscience Lidar Altimeter System (GLAS) is a laser altimeter designed to measure temporal changes in the topography and mass balance of the Earth's ice sheets. GLAS is scheduled to be launched as part of the Earth Observing System (EOS) in 2001. Due to the large pulse width of GLAS and the slope (and roughness) of the ice sheets, the mean elevation of the laser's surface spot will be estimated from the centroid of the return pulse. While the accuracy of individual altimetry measurements is expected to be 15 cm, a cross-over technique that averages the elevation differences at selected points on the ice sheets will allow GLAS to detect mean ice elevation changes as small as 2 cm per year. One factor that may degrade the proposed accuracy of the altimeter measurements is the effect of multiple scattering by thin clouds and aerosols over the ice sheets. In the return pulse, some photons from the lidar are slightly deflected by cloud particles but still return to the instrument detector. Since these photons travel a longer path than photons that pass directly to and from the surface, the mean travel time of the pulse is lengthened and the centroid of the return pulse is shifted toward a later time . This study intends to determine the importance of this source of altimetry error in the GLAS measurements. Two methods will be used to estimate path delays in the GLAS lidar for a range of cloud conditions. An analytic expression of a doubly scattered return signal (Eloranta, 1972) will be compared to Monte Carlo simulations that estimate path delays by cloud, aerosol and molecular scattering. The results from both methods will be discussed in relationship to current knowledge of Arctic cloud climatology, and the impact of multiple scattering on GLAS altimeter measurements will be evaluated.

Duda, David P.↗
NASA NTRS2014

Path Toward a Unifid Geometry for Radiation Transport

The Direct Accelerated Geometry for Radiation Analysis and Design (DAGRAD) element of the RadWorks Project under Advanced Exploration Systems (AES) within the Space Technology Mission Directorate (STMD) of NASA will enable new designs and concepts of operation for radiation risk assessment, mitigation and protection. This element is designed to produce a solution that will allow NASA to calculate the transport of space radiation through complex computer-aided design (CAD) models using the state-of-the-art analytic and Monte Carlo radiation transport codes. Due to the inherent hazard of astronaut and spacecraft exposure to ionizing radiation in low-Earth orbit (LEO) or in deep space, risk analyses must be performed for all crew vehicles and habitats. Incorporating these analyses into the design process can minimize the mass needed solely for radiation protection. Transport of the radiation fields as they pass through shielding and body materials can be simulated using Monte Carlo techniques or described by the Boltzmann equation, which is obtained by balancing changes in particle fluxes as they traverse a small volume of material with the gains and losses caused by atomic and nuclear collisions. Deterministic codes that solve the Boltzmann transport equation, such as HZETRN [high charge and energy transport code developed by NASA Langley Research Center (LaRC)], are generally computationally faster than Monte Carlo codes such as FLUKA, GEANT4, MCNP(X) or PHITS; however, they are currently limited to transport in one dimension, which poorly represents the secondary light ion and neutron radiation fields. NASA currently uses HZETRN space radiation transport software, both because it is computationally efficient and because proven methods have been developed for using this software to analyze complex geometries. Although Monte Carlo codes describe the relevant physics in a fully three-dimensional manner, their computational costs have thus far prevented their widespread use for analysis of complex CAD models, leading to the creation and maintenance of toolkit-specific simplistic geometry models. The work presented here builds on the Direct Accelerated Geometry Monte Carlo (DAGMC) toolkit developed for use with the Monte Carlo N-Particle (MCNP) transport code. The workflow for achieving radiation transport on CAD models using MCNP and FLUKA has been demonstrated and the results of analyses on realistic spacecraft/habitats will be presented. Future work is planned that will further automate this process and enable the use of multiple radiation transport codes on identical geometry models imported from CAD. This effort will enhance the modeling tools used by NASA to accurately evaluate the astronaut space radiation risk and accurately determine the protection provided by as-designed exploration mission vehicles and habitats

Lee, Kerry↗
NASA NTRS1972

Martian doublet craters.

A large number of Mars craters are nearly tangential to other craters. They occur in clusters or as isolated crater doublets. Results of probability calculations and a Monte Carlo cratering simulation model show conclusively that many of the Mars craters could not have resulted from random single-body impact. The possibility that these craters are calderas is considered possible only if calderas on Mars form by mechanisms different from those on earth. However, clusters and doublets could be caused by meteoroid breakup resulting from stresses induced in the meteoroid by the gravitational field of Mars. It is concluded that, under certain conditions, doublets should be produced on Mars as a direct result of breakup of an impacting meteoroid. The impact process can yield nonrandom crater distributions that should be observed in different degrees of development on different planetary surfaces.

Oberbeck, V. R.↗
NASA NTRS2017

Kinetic Monte Carlo Simulations of Oxygen Diffusion in Environmental Barrier Coating Materials

Ceramic Matrix Composite (CMC) materials are of interest for use in next-generation turbine engine components, offering a number of significant advantages, including reduced weight and high operating temperatures. However, in the hot environment in which such components operate, the presence of water vapor can lead to corrosion and recession, limiting the useful life of the components. Such degradation can be reduced through the use of Environmental Barrier Coatings (EBCs) that limit the amount of oxygen and water vapor reaching the component. Candidate EBC materials include Yttrium and Ytterbium silicates. In this work we present results of kinetic Monte Carlo (kMC) simulations of oxygen diffusion, via the vacancy mechanism, in Yttrium and Ytterbium disilicates, along with a brief discussion of interstitial diffusion. An EBC system typically includes a bond coat located between the EBC and the component surface. Bond coat materials are generally chosen for properties other than low oxygen diffusivity, but low oxygen diffusivity is nevertheless a desirable characteristic, as the bond coat could provide some additional component protection, particularly in the case where cracks in the coating system provide a direct path from the environment to the bond coat interface. We have therefore performed similar kMC simulations of oxygen diffusion in this material.

Good, Brian S.↗
NASA NTRS1988

Power spectral density analysis of wind-shear turbulence for related flight simulations

Meteorological phenomena known as microbursts can produce abrupt changes in wind direction and/or speed over a very short distance in the atmosphere. These changes in flow characteristics have been labelled wind shear. Because of its adverse effects on aerodynamic lift, wind shear poses its most immediate threat to flight operations at low altitudes. The number of recent commercial aircraft accidents attributed to wind shear has necessitated a better understanding of how energy is transferred to an aircraft from wind-shear turbulence. Isotropic turbulence here serves as the basis of comparison for the anisotropic turbulence which exists in the low-altitude wind shear. The related question of how isotropic turbulence scales in a wind shear is addressed from the perspective of power spectral density (psd). The role of the psd in related Monte Carlo simulations is also considered.

Laituri, Tony R.↗
NASA NTRS1999

Diffuse Reflection of Laser Light From Clouds

Laser light reflected from an aqueous suspension of particles or "cloud" with known thickness and particle size distribution defines the "cloud radiative Green's function", G. G is sensitive to cloud thickness, allowing retrieval of that important quantity. We describe a laboratory simulation of G, useful in design of an offbeam Lidar instrument for remote sensing of cloud thickness. Clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The size distribution extends from 0.5 microns to 25 microns, roughly lognormal, similar to real clouds. Density of suspended spheres is adjusted so photon mean-free-path is about 10 cm, 1000 times smaller than in real clouds. The light source is a Nd:YAG laser at 530 nm. Detectors are flux and photon-counting PMTs, with a glass probe for precise positioning. A Labview 5 VI controls position and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider,and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns. Step size is selectable. Far from the beam, the rate of exponential increase in the beam direction scales as expected from diffusion theory, linearly with cloud thickness, and inversely as the square root of the reduced optical thickness, independent of particle size. Nearer the beam the signal increases faster than exponential and depends on particle size. Results verify 3D Monte Carlo simulations that demonstrate detectability of remotely sensed offbeam returns, without filters at night, with narrow bandpass filter in day.

Cahalan, R. F.↗