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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Origin of Intense Electron Heating in Relativistic Blast Waves

The modeling of gamma-ray burst afterglow emission bears witness to strong electron heating in the precursor of Weibel-mediated, relativistic collisionless shock waves propagating in unmagnetized electron–ion plasmas. In this Letter, we propose a theoretical model, which describes electron heating via a Joule-like process caused by pitch-angle scattering in the decelerating, self-induced microturbulence and the coherent charge-separation field induced by the difference in inertia between electrons and ions. The emergence of this electric field across the precursor of electron–ion shocks is confirmed by large-scale particle-in-cell (PIC) simulations. Integrating the model using a Monte Carlo-Poisson method, we compare the main observables to the PIC simulations to conclude that the above mechanism can indeed account for the bulk of electron heating.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1974

Chain pooling to minimize prediction error in subset regression

Monte Carlo studies using population models intended to represent response surface applications are reported. Simulated experiments were generated by adding pseudo random normally distributed errors to population values to generate observations. Model equations were fitted to the observations and the decision procedure was used to delete terms. Comparison of values predicted by the reduced models with the true population values enabled the identification of deletion strategies that are approximately optimal for minimizing prediction errors.

Holms, A. G.↗
NASA NTRS1978

A general cratering-history model and its implications for the lunar highlands

Through analysis of a large number of Monte Carlo and Markov Chain simulations, a model for determining crater accumulation and crater obliteration histories has been derived. The model generally applies to populations of large craters. It predicts that the following relationships hold for subequilibrium-density crater populations: (1) the more negative the production function's exponent, alpha (N near D super alpha) the lower the crater density at which the population size-frequency distribution will significantly depart from its production function; (2) the more negative the production function's exponent, the less obliteration a crater population will sustain after a set number of impacts. Application of the model to the lunar highlands implies (1) the production function for the large craters is highly structured, resembling the observed size-frequency distribution and not the function N near D to the -2; (2) even the densely cratered highlands have not attained crater saturation or equilibrium. Direct simulations of the highlands' crater population supports the model's implications.

Woronow, A.↗
NASA NTRS1981

Verification of reflectance models in turbid waters

Inherent optical parameters of very turbid waters were used to evaluate existing water reflectance models. Measured upwelling radiance spectra and Monte Carlo simulations of the radiative transfer equations were compared with results from models based upon two flow, quasi-single scattering, augmented isotropic scattering, and power series approximation. Each model was evaluated for three separate components of upwelling radiance: (1) direct sunlight; (2) diffuse skylight; and (3) internally reflected light. Limitations of existing water reflectance models as applied to turbid waters and possible applications to the extraction of water constituent information are discussed.

Tanis, F. J.↗
NASA NTRS1983

Luminosity correlations in quasars

Simulations are conducted with and without flux thresholds in an investigation of quasar luminosity correlations by means of a Monte Carlo analysis, for various model distributions of quasars in X-rays and optical luminosity. For the case where the X-ray photons are primary, an anticorrelation between X-ray-to-optical luminosity ratio and optical luminosity arises as a natural consequence which resembles observations. The low optical luminosities of X-ray selected quasars can be understood as a consequence of the same effect, and similar conclusions may hold if the X-ray and optical luminosities are determined independently by a third parameter, although they do not hold if the optical photons are primary. The importance of such considerations is demonstrated through a reanalysis of the published X-ray-to-optical flux ratios for the 3CR sample.

Chanan, G. A.↗
NASA NTRS1984

Sensitivity of airborne fluorosensor measurements to linear vertical gradients in chlorophyll concentration

A semianalytic Monte Carlo radiative transfer simulation model for airborne laser fluorosensors has been extended to investigate the effects of inhomogeneities in the vertical distribution of phytoplankton concentrations in clear seawater. Simulation results for linearly varying step concentrations of chlorophyll are presented. The results indicate that statistically significant differences can be seen under certain conditions in the water Raman-normalized fluorescence signals between nonhomogeneous and homogeneous cases. A statistical test has been used to establish ranges of surface concentrations and/or verticl gradients in which calibration by surface samples would by inappropriate, and the results are discussed.

Venable, D. D.↗
NASA NTRS1986

Universal amplitude ratios and the interfacial tension near consolute points of binary liquid mixtures

The interfacial tension of the binary liquid mixtures triethylamine/heavy water, triethylamine/water, and cyclohexane/methanol near their consolute temperatures is determined experimentally by measuring the densities of the phases and the capillary lengths. The results are presented in tables and graphs and compared with published data obtained by other methods and with theoretical predictions. The findings are shown to be in good agreement with previous experiments, confirming the discrepancy between temperature-independent universal amplitude ratios calculated from such data and those predicted by renormalization-group calculations or Monte Carlo simulations of Ising models.

Moldover, M. R.↗
NASA NTRS1986

Neutral cometary atmospheres. V - C2 and CN in comets

The radial Haser scale lengths for the production and destruction of C2 and CN are calculated on the basis of brightness profiles for comets Bennett 1970 II and Kohoutek 1973 XII, applying the Monte Carlo particle-trajectory model of Combi and Delsemme (1980). Findings reported include C2/CN ratio independent of heliocentric distance, C2 ejection velocity from its parent molecule about 0.5 km/s, and 1-AU photochemical lifetime 31,000 s for the C2 parent molecule and 120,000 s for C2.

Combi, M. R.↗
NASA NTRS1988

Evaluation of concepts and methods for extraction of flow parameters from laser transit anemometry data

A Monte-Carlo type simulation system was developed to test concepts and associated algorithms proposed by Mayo and Smart (1984) for the extraction of two-dimensional flow parameters from laser transit anemometry data sets. The Monte-Carlo LTA simulation program models a flow field as well as the sample volume geometry of the LTA system, thereby permitting simulation of two-dimensional data acquisition in a three-dimensional flow field. The results of these tests on the proposed processing concepts suggest that the instrument is capable of measuring mean velocities to less than +0.15 percent error and flow angles to less then +0.3 percent for turbulence intensities of up to 15.0 percent.

Humphreys, W. M., Jr.↗
NASA NTRS1991

An Exploratory Exercise in Taguchi Analysis of Design Parameters: Application to a Shuttle-to-space Station Automated Approach Control System

The chief goals of the summer project have been twofold - first, for my host group and myself to learn as much of the working details of Taguchi analysis as possible in the time allotted, and, secondly, to apply the methodology to a design problem with the intention of establishing a preliminary set of near-optimal (in the sense of producing a desired response) design parameter values from among a large number of candidate factor combinations. The selected problem is concerned with determining design factor settings for an automated approach program which is to have the capability of guiding the Shuttle into the docking port of the Space Station under controlled conditions so as to meet and/or optimize certain target criteria. The candidate design parameters under study were glide path (i.e., approach) angle, path intercept and approach gains, and minimum impulse bit mode (a parameter which defines how Shuttle jets shall be fired). Several performance criteria were of concern: terminal relative velocity at the instant the two spacecraft are mated; docking offset; number of Shuttle jet firings in certain specified directions (of interest due to possible plume impingement on the Station's solar arrays), and total RCS (a measure of the energy expended in performing the approach/docking maneuver). In the material discussed here, we have focused on single performance criteria - total RCS. An analysis of the possibility of employing a multiobjective function composed of a weighted sum of the various individual criteria has been undertaken, but is, at this writing, incomplete. Results from the Taguchi statistical analysis indicate that only three of the original four posited factors are significant in affecting RCS response. A comparison of model simulation output (via Monte Carlo) with predictions based on estimated factor effects inferred through the Taguchi experiment array data suggested acceptable or close agreement between the two except at the predicted optimum point, where a difference outside a rule-of-thumb bound was observed. We have concluded that there is most likely an interaction effect not provided for in the original orthogonal array selected as the basis for our experimental design. However, we feel that the data indicates that this interaction is a mild one and that inclusion of its effect will not alter the location of the optimum.

Deal, Don E.↗
NASA NTRS2002

Dynamic Stability of Uncertain Laminated Beams Under Subtangential Loads

Because of the inherent complexity of fiber-reinforced laminated composites, it can be challenging to manufacture composite structures according to their exact design specifications, resulting in unwanted material and geometric uncertainties. In this research, we focus on the deterministic and probabilistic stability analysis of laminated structures subject to subtangential loading, a combination of conservative and nonconservative tangential loads, using the dynamic criterion. Thus a shear-deformable laminated beam element, including warping effects, is derived to study the deterministic and probabilistic response of laminated beams. This twenty-one degrees of freedom element can be used for solving both static and dynamic problems. In the first-order shear deformable model used here we have employed a more accurate method to obtain the transverse shear correction factor. The dynamic version of the principle of virtual work for laminated composites is expressed in its nondimensional form and the element tangent stiffness and mass matrices are obtained using analytical integration The stability is studied by giving the structure a small disturbance about an equilibrium configuration, and observing if the resulting response remains small. In order to study the dynamic behavior by including uncertainties into the problem, three models were developed: Exact Monte Carlo Simulation, Sensitivity Based Monte Carlo Simulation, and Probabilistic FEA. These methods were integrated into the developed finite element analysis. Also, perturbation and sensitivity analysis have been used to study nonconservative problems, as well as to study the stability analysis, using the dynamic criterion.

Goyal, Vijay K.↗
NASA NTRS2001

Predictions of secondary neutrons and their importance to radiation effects inside the International Space Station

As part of a study funded by NASA MSFC to assess thecontribution of secondary particles in producing radiation damage to optoelectronics devices located on the International Space Station (IS), Monte Carlo calculations have been made to predict secondary spectra vs. shielding inside ISS modules and in electronics boxes attached on the truss (Armstrong and Colborn, 1998). The calculations take into account secondary neutron, proton, and charged pion production from the ambient galactic cosmic-ray (GCR) proton, trapped proton, and neutron albedo environments. Comparisons of the predicted neutron spectra with measurments made on the Mir space station and other spacecraft have also been made (Armstrong and Colborn, 1998). In this paper, some initial results from folding the predicted neutron spectrum inside ISS modules from Armstrong and Colborn (1998) with several types of radiation effects response functions related to electronics damage and astronaut-dose are given. These results provide an estimate of the practical importance of neutrons compared to protons in assessing radiation effects for the ISS. Also, the important neutron energy ranges for producing these effects have been estimated, which provides guidance for onboard neutron measurement requirements.

Neutrons↗
NASA NTRS2005

Characterization of Space Shuttle Ascent Debris Aerodynamics Using CFD Methods

An automated Computational Fluid Dynamics process for determining the aerodynamic Characteristics of debris shedding from the Space Shuttle Launch Vehicle during ascent is presented. This process uses Cartesian fully-coupled, six-degree-of-freedom simulations of isolated debris pieces in a Monte Carlo fashion to produce models for the drag and crossrange behavior over a range of debris shapes and shedding scenarios. A validation of the Cartesian methods against ballistic range data for insulating foam debris shapes at flight conditions, as well as validation of the resulting models, are both contained. These models are integrated with the existing shuttle debris transport analysis software to provide an accurate and efficient engineering tool for analyzing debris sources and their potential for damage.

Murman, Scott M.↗
NASA NTRS2007

GCR-induced Photon Luminescence of the Moon: The Moon as a CR Detector

We report on the results of a preliminary study of the GCR-induced photon luminescence of the Moon using the Monte Carlo program FLUKA. The model of the lunar surface is taken to be the chemical composition of soils found at various landing sites during the Apollo and Luna programs, averaged over all such sites to define a generic regolith for the present analysis. This then becomes the target that is bombarded by Galactic Cosmic Rays (GCRs) in FLUKA to determine the photon fluence when there is no sunshine or Earthshine. From the photon fluence we derive the energy spectrum which can be utilized to design an orbiting optical instrument for measuring the GCR-induced luminescence. This is to be distinguished from the gamma-ray spectrum produced by the radioactive decay of its radiogenic constituents lying in the surface and interior. Also, we investigate transient optical flashes from high-energy CRs impacting the lunar surface (boulders and regolith). The goal is to determine to what extent the Moon could be used as a rudimentary CR detector. Meteor impacts on the Moon have been observed for centuries to generate such flashes, so why not CRs?

Wilson, Thomas L.↗
NASA NTRS2007

Particle Analysis Pitfalls

This viewgraph presentation reviews the use of particle analysis to assist in preparing for the 4th Hubble Space Telescope (HST) Servicing mission. During this mission the Space Telescope Imaging Spectrograph (STIS) will be repaired. The particle analysis consisted of Finite element mesh creation, Black-body viewfactors generated using I-DEAS TMG Thermal Analysis, Grey-body viewfactors calculated using Markov method, Particle distribution modeled using an iterative Monte Carlo process, (time-consuming); in house software called MASTRAM, Differential analysis performed in Excel, and Visualization provided by Tecplot and I-DEAS. Several tests were performed and are reviewed: Conformal Coat Particle Study, Card Extraction Study, Cover Fastener Removal Particle Generation Study, and E-Graf Vibration Particulate Study. The lessons learned during this analysis are also reviewed.

Hughes, David↗
NASA NTRS2008

Space-radiation-induced Photon Luminescence of the Moon

We report on the results of a study of the photon luminescence of the Moon induced by Galactic Cosmic Rays (GCRs) and space radiation from the Sun, using the Monte Carlo program FLUKA. The model of the lunar surface is taken to be the chemical composition of soils found at various landing sites during the Apollo and Luna programs, averaged over all such sites to define a generic regolith for the present analysis. This then becomes the target that is bombarded by Galactic Cosmic Rays (GCRs) and Solar Energetic Particles (SEPs) above 1 keV in FLUKA to determine the photon fluence albedo produced by the Moon's surface when there is no sunlight and Earthshine. This is to be distinguished from the gamma-ray spectrum produced by the radioactive decay of radiogenic constituents lying in the surface and interior of the Moon. From the photon fluence we derive the spectrum which can be utilized to examine existing lunar spectral data and to design orbiting instrumentation for measuring various components of the space-radiation-induced photon luminescence present on the Moon.

Wilson, Thomas↗