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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Monte Carlo simulation of nonequilibrium shock fronts

A direct simulation of a shock-tube experiment carried by AVCO (pressure of 1 Torr and velocity of 6.4 km/sec), repeated by Sharma and Gillespie (1990), and used by Park (1988) to develop his two-temperature model, is presented. Results show that the electronic ground state of N2 is not in rotational nor vibrational equilibrium with that of N2(+). Moreover, a two-temperature model is inadequate to describe nonequilibrium flows behind shocks. The role of impurities is examined. It is shown that the effects of a small fraction of H2O are insignificant. Good agreement with the measurements of Sharma and Gillespie is indicated.

Olynick, David P.↗
NASA NTRS1991

Estimation variance bounds of importance sampling simulations in digital communication systems

In practical applications of importance sampling (IS) simulation, two basic problems are encountered, that of determining the estimation variance and that of evaluating the proper IS parameters needed in the simulations. The authors derive new upper and lower bounds on the estimation variance which are applicable to IS techniques. The upper bound is simple to evaluate and may be minimized by the proper selection of the IS parameter. Thus, lower and upper bounds on the improvement ratio of various IS techniques relative to the direct Monte Carlo simulation are also available. These bounds are shown to be useful and computationally simple to obtain. Based on the proposed technique, one can readily find practical suboptimum IS parameters. Numerical results indicate that these bounding techniques are useful for IS simulations of linear and nonlinear communication systems with intersymbol interference in which bit error rate and IS estimation variances cannot be obtained readily using prior techniques.

Lu, D.↗
NASA NTRS1988

Precision in the perception of direction of a moving pattern

The precision of the model of pattern motion analysis put forth by Adelson and Movshon (1982) who proposed that humans determine the direction of a moving plaid (the sum of two sinusoidal gratings of different orientations) in two steps is qualitatively examined. The volocities of the grating components are first estimated, then combined using the intersection of constraints to determine the velocity of the plaid as a whole. Under the additional assumption that the noise sources for the component velocities are independent, an approximate expression can be derived for the precision in plaid direction as a function of the precision in the speed and direction of the components. Monte Carlo simulations verify that the expression is valid to within 5 percent over the natural range of the parameters. The expression is then used to predict human performance based on available estimates of human precision in the judgment of single component speed. Human performance is predicted to deteriorate by a factor of 3 as half the angle between the wavefronts (theta) decreases from 60 to 30 deg, but actual performance does not. The mean direction discrimination for three human observers was 4.3 plus or minus 0.9 deg (SD) for theta = 60 deg and 5.9 plus or minus 1.2 for theta = 30 deg. This discrepancy can be resolved in two ways. If the noises in the internal representations of the component speeds are smaller than the available estimates or if these noises are not independent, then the psychophysical results are consistent with the Adelson-Movshon hypothesis.

Stone, Leland S.↗
NASA NTRS2010

Optical See-Through Head Mounted Display Direct Linear Transformation Calibration Robustness in the Presence of User Alignment Noise

Augmented Reality (AR) is a technique by which computer generated signals synthesize impressions that are made to coexist with the surrounding real world as perceived by the user. Human smell, taste, touch and hearing can all be augmented, but most commonly AR refers to the human vision being overlaid with information otherwise not readily available to the user. A correct calibration is important on an application level, ensuring that e.g. data labels are presented at correct locations, but also on a system level to enable display techniques such as stereoscopy to function properly [SOURCE]. Thus, vital to AR, calibration methodology is an important research area. While great achievements already have been made, there are some properties in current calibration methods for augmenting vision which do not translate from its traditional use in automated cameras calibration to its use with a human operator. This paper uses a Monte Carlo simulation of a standard direct linear transformation camera calibration to investigate how user introduced head orientation noise affects the parameter estimation during a calibration procedure of an optical see-through head mounted display.

Axholt, Magnus↗
NASA NTRS1994

Implications of the BATSE data for a helocentric origin of gamma-ray bursts

Several possible models have been suggested to explain the observed distribution of gamma-ray bursts: heliocentric distributions such as the Oort cloud, large galactic halos, and cosmological models. We report here on an investigation into the implications of the Burst and Transient Source Experiment (BATSE) gamma-ray burst distribution (Meegan et al. 1992a) data on the possible helocentric origin of gamma-ray bursts. We find no statistically significant anisotropy in the angular distribution of the bursts in a Sun-referenced coordinate system; there is no dipole moment in the direction of the Sun, and no quardrupole moment associated with the ecliptic plane. We have employed direct analytic calculations and Monte Carlo simulations of sources in the Oort cloud to constrain possible helicentric burst distributions. These can produce distributions consistent with the observed angular isotropy, the meal value of V/V(sub max), and the observed C/C(sub min) distribution of BATSE, and provide limits to burst energy of a few times approximately 10(exp 27) ergs. However, the agreement of the heliocentric C/C(sub min) distributions with the BATSE data is attributable to the relatively limited sampling of strong, nearby bursts. These bursts are known from observation to be homogeneously distributed, yet the density of sources in the Oort cloud is not constant in this region. Integral number-intensity distributions from the Oort cloud for larger numbers of bursts cannot reproduce the known homogeneity of the strong bursts without modification to the computed cometary number density and are therefore unlikely explanations of the gamma-ray burst distribution.

Horack, J. M.↗
NASA NTRS1978

Information content of sky intensity and polarization measurements at right angles to the solar direction

The paper presents the results of a Monte Carlo simulation study of the brightness and polarization at right angles to the solar direction both for ground-based observations (looking up) and for satellite-based systems (looking down). Calculations have been made for a solar zenith angle whose cosine was 0.6 and wavelengths ranging from 3500 A to 9500 A. A sensitivity of signatures to total aerosol loading, aerosol particle size distribution and refractive index, and the surface reflectance albedo has been demonstrated. For Lambertian-type surface reflection the albedo effects enter solely through the intensity sensitivity, and very high correlations have been found between the polarization term signatures for the ground-based and satellite-based systems. Potential applications of these results for local albedo predictions and satellite imaging systems recalibrations are discussed.

Holland, A. C.↗
NASA NTRS2011

Development of the Sasquatch Drop Test Footprint Tool

The Crew Exploration Vehicle Parachute Assembly System (CPAS) is the parachute system for NASA s Orion spacecraft. CPAS is currently in the design and testing phase of development. The test program consists of numerous drop tests, wherein a test article rigged with parachutes is extracted or released from an aircraft. During such tests, range safety is paramount, as is the recoverability of the parachutes and test article. It is crucial to establish an aircraft release point that will ensure that the article and all items released from it will land in safe locations. Early in the CPAS project, a legacy tool (previously used on the X-38 project) was used to determine a safe release point and to predict the landing locations (the footprint) of the payload and all released objects. Due to increasing test complexity and the need for a more flexible tool, a new footprint predictor tool, called Sasquatch, was created in MATLAB. This tool takes in a simulated trajectory for the test article, information about all released objects, and atmospheric wind data (simulated or actual) to calculate the trajectories of the released objects. Dispersions are applied to the landing locations of those objects, taking into account the variability of winds, aircraft release point, and object descent rate. A safe aircraft release point is determined based on the landing locations of the payload and released objects. The release point, landing locations, and dispersions are plotted on a simple map of the drop zone for easy reference. To date, Sasquatch has been used for thirteen drop tests. Comparing the predictions with actual test results has allowed for significant improvements in the tool s predictive capabilities, especially the incorporation of a well-correlated horizontal throw model. Intended future improvements to the tool include tighter dispersions on the landing locations, Monte Carlo capability, direct input from trajectory simulations, and a graphical user interface.

Bledsoe, Kristin J.↗
NASA NTRS1993

Transient motion of a confined rarefied gas due to wall heating or cooling

The direct Monte Carlo method is used to simulate the transient motion of a confined rarefied gas as a container wall is rapidly heated or cooled. Results show that, for the rapid wall temperature changes, the gas response is better measured by an acoustic, rather than conduction, timescale, with steady state conditions being achieved after about ten wave crossing times.

Wadsworth, Dean C.↗
NASA NTRS2010

Optical Imaging and Radiometric Modeling and Simulation

OPTOOL software is a general-purpose optical systems analysis tool that was developed to offer a solution to problems associated with computational programs written for the James Webb Space Telescope optical system. It integrates existing routines into coherent processes, and provides a structure with reusable capabilities that allow additional processes to be quickly developed and integrated. It has an extensive graphical user interface, which makes the tool more intuitive and friendly. OPTOOL is implemented using MATLAB with a Fourier optics-based approach for point spread function (PSF) calculations. It features parametric and Monte Carlo simulation capabilities, and uses a direct integration calculation to permit high spatial sampling of the PSF. Exit pupil optical path difference (OPD) maps can be generated using combinations of Zernike polynomials or shaped power spectral densities. The graphical user interface allows rapid creation of arbitrary pupil geometries, and entry of all other modeling parameters to support basic imaging and radiometric analyses. OPTOOL provides the capability to generate wavefront-error (WFE) maps for arbitrary grid sizes. These maps are 2D arrays containing digital sampled versions of functions ranging from Zernike polynomials to combination of sinusoidal wave functions in 2D, to functions generated from a spatial frequency power spectral distribution (PSD). It also can generate optical transfer functions (OTFs), which are incorporated into the PSF calculation. The user can specify radiometrics for the target and sky background, and key performance parameters for the instrument s focal plane array (FPA). This radiometric and detector model setup is fairly extensive, and includes parameters such as zodiacal background, thermal emission noise, read noise, and dark current. The setup also includes target spectral energy distribution as a function of wavelength for polychromatic sources, detector pixel size, and the FPA s charge diffusion modulation transfer function (MTF).

Ha, Kong Q.↗
NASA NTRS2012

Pulsar Pair Cascades in Magnetic Fields with Offset Polar Caps

Neutron star magnetic fields may have polar caps (PC) that are offset from the dipole axis, through field-line sweepback near the light cylinder or non-symmetric currents within the star. The effects of such offsets on electron-positron pair cascades are investigated, using simple models of dipole magnetic fields with small distortions that shift the PCs by different amounts or directions. Using a Monte Carlo pair cascade simulation, we explore the changes in the pair spectrum, multiplicity and energy flux across the PC, as well as the trends in pair flux and pair energy flux with spin-down luminosity, L(sub sd). We also give an estimate of the distribution of heating flux from returning positrons on the PC for different offsets. We find that even modest offsets can produce significant increases in pair multiplicity, especially for pulsars that are near or beyond the pair death lines for centered PCs, primarily because of higher accelerating fields. Pair spectra cover several decades in energy, with the spectral range of millisecond pulsars (MSPs) two orders of magnitude higher than for normal pulsars, and PC offsets allow significant extension of all spectra to lower pair energies. We find that the total PC pair luminosity L(sub pair) is proportional to L(sub sd), with L(sub pair) approximates 10(exp -3) L(sub sd) for normal pulsars and L(sub pair) approximates 10(exp -2) L(sub sd) for MSPs. Remarkably, the total PC heating luminosity for even large offsets increases by less than a factor of two, even though the PC area increases by much larger factors, because most of the heating occurs near the magnetic axis.

Harding, Alice K.↗
NASA NTRS2016

Progress of Aircraft System Noise Assessment with Uncertainty Quantification for the Environmentally Responsible Aviation Project

Aircraft system noise predictions have been performed for NASA modeled hybrid wing body aircraft advanced concepts with 2025 entry-into-service technology assumptions. The system noise predictions developed over a period from 2009 to 2016 as a result of improved modeling of the aircraft concepts, design changes, technology development, flight path modeling, and the use of extensive integrated system level experimental data. In addition, the system noise prediction models and process have been improved in many ways. An additional process is developed here for quantifying the uncertainty with a 95% confidence level. This uncertainty applies only to the aircraft system noise prediction process. For three points in time during this period, the vehicle designs, technologies, and noise prediction process are documented. For each of the three predictions, and with the information available at each of those points in time, the uncertainty is quantified using the direct Monte Carlo method with 10,000 simulations. For the prediction of cumulative noise of an advanced aircraft at the conceptual level of design, the total uncertainty band has been reduced from 12.2 to 9.6 EPNL dB. A value of 3.6 EPNL dB is proposed as the lower limit of uncertainty possible for the cumulative system noise prediction of an advanced aircraft concept.

Thomas, Russell H.↗
NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Physics-based modeling of hypersonic flows is predicated on the availability of chemical reaction rate coefficients and cross sections for the collisional processes. This approach has been built around the use of quantum mechanical calculations to describe the interaction between the colliding particles. In this approach a potential energy surface (PES) is computed by solving the electronic Schrödinger equation and collision cross sections are determined for that PES using classical, semiclassical or quantum mechanical scattering methods. The rate coefficients are computed by integrating the thermally weighted cross sections. State-to-state rate coefficients are determined by only integrating over a thermal distribution of collisional energies. Finally, thermal rate coefficients are determined by summation of the state-to-state rate coefficients for reactions of molecules in all relevant ro-vibrational energy levels. If the flow is in thermal non-equilibrium, the translational, vibrational and rotational energy modes can be represented in different ways: three unique temperatures can be used to describe the distributions, the populations of individual ro-vibrational energy levels can be determined by solving the Master Equation, or through the use of direct simulation in particle-based Monte Carlo sampling. The PES-to-rate coefficient approach had been proposed and attempted in the early days of digital computing, but it is only in the last 15 years that computer hardware and software have been up to the task of calculating accurate interatomic and intermolecular potentials.

Jaffe, Richard L.↗
NASA NTRS1978

Collisional evolution of mass-distribution spectrum of planetesimals

The early growth of planetesimals by mutual direct collisions is numerically simulated with a Monte Carlo technique to show how planetesimals with specific mass and velocity distributions evolve into a full-size planet. Four types of collisions are taken into account: rebound, erosion, catastrophic break-up, and coagulation. It is shown that evolution of mass-distribution spectrum is expressed by an inverse power relation. It is also suggested that whether or not a planetesimal can survive catastrophic collision is primarily dependent on mean relative velocity and mechanical properties of planetesimals. It is necessary for the early growth of rocky (basaltic) materials that mean relative velocity be much smaller than 0.1 km/s. Otherwise it is necessary to introduce something like nucleating agents (such as iron bodies which have plastic properties at temperatures higher than 200 K) for the formation of terrestrial planets.

Matsui, T.↗
NASA NTRS1980

Solar wind sputtering effects in the atmospheres of Mars and Venus

It is found through an investigation, combining Monte Carlo simulations and analytical techniques, of the direct collisional interaction of an energetic particle flux with the neutral components of a planetary atmosphere, that solar wind sputtering could provide an important exospheric mass sink on both Mars and Venus under appropriate conditions. The computed rates of helium loss in the Venusian atmosphere and of carbon, nitrogen and oxygen in the Martian atmosphere imply that sputtering would have a significant effect on the noble gas budget of Venus and dominate the chemical and photochemical loss processes of Mars. Because of diffusive separation of lighter elements and isotopes, and because the gravitational binding energy is proportional to the mass, the erosion process preferentially removes the lighter components of the atmosphere. Solar wind sputtering could therefore compete with other erosion mechanisms in generating substantial fractionation effects.

Watson, C. C.↗
NASA NTRS2009

KSC Tech Transfer News, Volume 2, No. 2

This issue contains articles about: (1) the Innovative Partnerships Program (IPP) and the manager of the program, Alexis Hongamen, (2) New Technology Report (NTR) on a Monte Carlo Simulation to Estimate the Likelihood of Direct Lightning Strikes, (3) Kennedy Space Center's Applied Physics Lab, (4) a virtual ruler that is used for many applications, (5) a portable device that finds low-level leaks, (6) a sun-shield, that supports in-space cryogenic propellant storage, (7) lunar dust modeling software, (8) space based monitoring of radiation damage to DNA, (9) the use of light-emitting diode (LED) arrays vegetable production system, (10) Dust Tolerant Intelligent Electrical Connection Systems, (11) Ice Detection Camera System Upgrade, (12) Repair Techniques for Composite Structures, (13) Cryogenic Orbital Testbed, and (14) copyright protection.

Makufka, David↗
NASA NTRS1989

Monte Carlo simulation of particle acceleration at astrophysical shocks

A Monte Carlo code was developed for the simulation of particle acceleration at astrophysical shocks. The code is implemented in Turbo Pascal on a PC. It is modularized and structured in such a way that modification and maintenance are relatively painless. Monte Carlo simulations of particle acceleration at shocks follow the trajectories of individual particles as they scatter repeatedly across the shock front, gaining energy with each crossing. The particles are assumed to scatter from magnetohydrodynamic (MHD) turbulence on both sides of the shock. A scattering law is used which is related to the assumed form of the turbulence, and the particle and shock parameters. High energy cosmic ray spectra derived from Monte Carlo simulations have observed power law behavior just as the spectra derived from analytic calculations based on a diffusion equation. This high energy behavior is not sensitive to the scattering law used. In contrast with Monte Carlo calculations diffusive calculations rely on the initial injection of supra-thermal particles into the shock environment. Monte Carlo simulations are the only known way to describe the extraction of particles directly from the thermal pool. This was the triumph of the Monte Carlo approach. The question of acceleration efficiency is an important one in the shock acceleration game. The efficiency of shock waves efficient to account for the observed flux of high energy galactic cosmic rays was examined. The efficiency of the acceleration process depends on the thermal particle pick-up and hence the low energy scattering in detail. One of the goals is the self-consistent derivation of the accelerated particle spectra and the MHD turbulence spectra. Presumably the upstream turbulence, which scatters the particles so they can be accelerated, is excited by the streaming accelerated particles and the needed downstream turbulence is convected from the upstream region. The present code is to be modified to include a better description of particle scattering (pitch-angle instead of hard-sphere) and as iterative procedure for treating the self-excitation of the MHD turbulence.

Campbell, Roy K.↗
NASA NTRS1991

Probabilistic simulation of the human factor in structural reliability

Many structural failures have occasionally been attributed to human factors in engineering design, analyses maintenance, and fabrication processes. Every facet of the engineering process is heavily governed by human factors and the degree of uncertainty associated with them. Factors such as societal, physical, professional, psychological, and many others introduce uncertainties that significantly influence the reliability of human performance. Quantifying human factors and associated uncertainties in structural reliability require: (1) identification of the fundamental factors that influence human performance, and (2) models to describe the interaction of these factors. An approach is being developed to quantify the uncertainties associated with the human performance. This approach consists of a multi factor model in conjunction with direct Monte-Carlo simulation.

Shah, Ashwin R.↗
NASA NTRS2000

Model for radial dependence of frequency distributions for energy imparted in nanometer volumes from HZE particles

This paper develops a deterministic model of frequency distributions for energy imparted (total energy deposition) in small volumes similar to DNA molecules from high-energy ions of interest for space radiation protection and cancer therapy. Frequency distributions for energy imparted are useful for considering radiation quality and for modeling biological damage produced by ionizing radiation. For high-energy ions, secondary electron (delta-ray) tracks originating from a primary ion track make dominant contributions to energy deposition events in small volumes. Our method uses the distribution of electrons produced about an ion's path and incorporates results from Monte Carlo simulation of electron tracks to predict frequency distributions for ions, including their dependence on radial distance. The contribution from primary ion events is treated using an impact parameter formalism of spatially restricted linear energy transfer (LET) and energy-transfer straggling. We validate our model by comparing it directly to results from Monte Carlo simulations for proton and alpha-particle tracks. We show for the first time frequency distributions of energy imparted in DNA structures by several high-energy ions such as cosmic-ray iron ions. Our comparison with results from Monte Carlo simulations at low energies indicates the accuracy of the method.

Non-NASA Center↗