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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

MCRadar: A Monte Carlo Solver for Cloud and Precipitation Radar

Multiple scattering produces anomalous echoes in observed radar profiles that cannot be explained by other phenomena. These effects are most obvious for spaceborne platforms, at shorter wavelengths, and in convection, as multiple scattering is governed by antenna beam width, optical depth, and albe do; however, multiple scattering has been observed in a range of precipitating conditions. To account for the effects of multiple scattering, Monte Carlo integration is employed in a flexible framework to enable arbitrary radar configurations with finite Gaussian beams for three-dimensional atmospheric scenarios. The three-dimensional nature of the model, coupled with the finite antenna response, also allows for consideration of nonuniform beam filling effects that often coincide with multiple scattering. Examples of airborne and spaceborne radars at various wavelengths are used to illustrate the effects of multiple scattering and nonuniform beam filling. The MCRadar code is currently available in the development version of ARTS.

Adams, Ian Stuart↗
NASA NTRS2019

The Effects of High-Temperature Fractional Crystallization on Calcium Isotopic Composition

We report Ca isotope fractionation among a co-genetic suite of samples from the Mauritanian Guelb el Azib ultramafic-mafic-anorthosite complex (GAC), which represents the fractional crystallization sequence of an igneous magma chamber. We investigate how the composition of an evolving crystallizing liquid and the resultant mineralogy of co-crystallizing assemblages affects inter-phase Ca isotope fractionation. Because the GAC is an Archean metamorphosed complex, we first investigate the resistance of the Ca isotope signature to secondary hydrothermal alteration and to metamorphism. While we cannot conclude here that Ca isotopes are undisturbed at the mineral scale, we show that they are preserved at the bulk rock scale. This is adequate for our study where several samples are almost monomineralic at the bulk scale level. The δ(exp 44/40)Ca (relative to NIST standard SRM 915a) of GAC layers range from −1.53 to 1.61‰, with the earliest-formed ultramafic cumulate layer being isotopically heaviest and the later, more felsic layers being isotopically lighter. For the first time, we show that plagioclase is much more enriched in light Ca isotopes compared to olivine and (Ca-) pyroxene. Monte Carlo simulations suggest that Ca isotopes are fractionated among co-existing silicate minerals during fractional crystallization, possibly in relation with the residual liquid composition. In qualitative agreement with computational models based on first principles lattice dynamics, we observe that Ca isotope fractionation is mineralogically controlled and importantly, that the degree of fractionation can vary according to the CaO composition of the residual liquid. While previous studies have aimed to understand partial melting as a source of Ca isotope fractionation, our results suggest that fractional crystallization is also a source of Ca isotope variability among co-genetic samples.

Maria C. Valdes↗
NASA NTRS2015

Space Radiation Transport Code Development: 3DHZETRN

The space radiation transport code, HZETRN, has been used extensively for research, vehicle design optimization, risk analysis, and related applications. One of the simplifying features of the HZETRN transport formalism is the straight-ahead approximation, wherein all particles are assumed to travel along a common axis. This reduces the governing equation to one spatial dimension allowing enormous simplification and highly efficient computational procedures to be implemented. Despite the physical simplifications, the HZETRN code is widely used for space applications and has been found to agree well with fully 3D Monte Carlo simulations in many circumstances. Recent work has focused on the development of 3D transport corrections for neutrons and light ions (Z < 2) for which the straight-ahead approximation is known to be less accurate. Within the development of 3D corrections, well-defined convergence criteria have been considered, allowing approximation errors at each stage in model development to be quantified. The present level of development assumes the neutron cross sections have an isotropic component treated within N explicit angular directions and a forward component represented by the straight-ahead approximation. The N = 1 solution refers to the straight-ahead treatment, while N = 2 represents the bi-directional model in current use for engineering design. The figure below shows neutrons, protons, and alphas for various values of N at locations in an aluminum sphere exposed to a solar particle event (SPE) spectrum. The neutron fluence converges quickly in simple geometry with N > 14 directions. The improved code, 3DHZETRN, transports neutrons, light ions, and heavy ions under space-like boundary conditions through general geometry while maintaining a high degree of computational efficiency. A brief overview of the 3D transport formalism for neutrons and light ions is given, and extensive benchmarking results with the Monte Carlo codes Geant4, FLUKA, and PHITS are provided for a variety of boundary conditions and geometries. Improvements provided by the 3D corrections are made clear in the comparisons. Developments needed to connect 3DHZETRN to vehicle design and optimization studies will be discussed. Future theoretical development will relax the forward plus isotropic interaction assumption to more general angular dependence.

Wilson, John W.↗
NASA NTRS1971

Nuclear electromagnetic cascades from nuclei with Z larger than or equal to 3

A Monte Carlo simulation method was developed for studying nuclear-electromagnetic cascades initiated by high energy nuclei with Z or = 3 incident on heavy absorbers. The calculations are based on a cascade model which was first adjusted until it agreed with measurements made with protons at an accelerator. Modifications of the model used for protons include the incorporation of the probabilities for fragmentation of heavy nuclei into lighter nuclei, alpha particles, and nucleons. Mean values and fluctuations of the equivalent numbers of particles in the cascades at various depths in an iron absorber are presented for protons, carbon, and iron nuclei over the 30 to 300 GeV/nucleon energy range.

Jones, W. V.↗
NASA NTRS1976

Remote sensing of chlorophyll concentration: State-of-the-art, 1975

Remote measurement of chlorophyll concentration of the world's oceans from satellite observations could potentially be extremely useful for assessments of productivity in large areas for which measurements by other means would be impractical. The basis of these measurements rests with the physics of the interaction of light with material dissolved and suspended in the water. It is theoretically possible to predict the nature of light upwelled from the ocean surface from a solution to the radiative transfer equation. Practically, however, this is difficult. Monte-Carlo methods presently are thought to be the most viable method to treat the general theoretical problem. With restrictive assumptions of the nature of scattering, it is possible to construct simpler models. Algorithms developed to relate chlorophyll concentration (or some other parameter, i.e., seechi depth) to the upwelled light spectrum are discussed.

Atwell, B. H.↗
NASA NTRS1979

Physical and dynamical evolution of long-period comets

The source of long-period comets was investigated using a Monte Carlo simulation of comet evolution under the influence of a combination of physical and dynamical processes. The perturbation of cometary orbits by major planets and by non-gravitational forces was modeled, as was physical loss of comets due to random disruption (splitting) and planetary collision; a model was also derived for loss of all volatiles. The importance of each of these processes was examined. The primary end states found for long-period comets were: ejection on hyperbolic orbit, 65.2%; random disruption, 27.6%; and formation of silicate crusts, 7.1%. The basic correctness of the Ort hypothesis was confirmed.

Weissman, P. R.↗
NASA NTRS1980

Assessment of risk due to the use of carbon fiber composites in commercial and general aviation

The development of a national risk profile for the total annual aircraft losses due to carbon fiber composite (CFC) usage through 1993 is discussed. The profile was developed using separate simulation methods for commercial and general aviation aircraft. A Monte Carlo method which was used to assess the risk in commercial aircraft is described. The method projects the potential usage of CFC through 1993, investigates the incidence of commercial aircraft fires, models the potential release and dispersion of carbon fibers from a fire, and estimates potential economic losses due to CFC damaging electronic equipment. The simulation model for the general aviation aircraft is described. The model emphasizes variations in facility locations and release conditions, estimates distribution of CFC released in general aviation aircraft accidents, and tabulates the failure probabilities and aggregate economic losses in the accidents.

Fiksel, J.↗
NASA NTRS1983

Highlights of the SEASAT-SASS program - A review

Some important concepts of the SEASAT-SASS program are described and some of the decisions made during the program as to methods for relating wind to backscatter are discussed. The radar scatterometer design is analyzed along with the model function, which is an empirical relationship between the backscatter value and the wind speed, wind direction, and incidence angle of the radar beam with the sea surface. The results of Monte Carlo studies of mesoscale turbulence and of studies of wind stress on the sea surface involving SASS are reviewed.

Pierson, W. J., Jr.↗
NASA NTRS1990

A numerical algorithm to propagate navigation error covariance matrices associated with generalized strapdown inertial measurement units

The design and operation of a Strapdown Navigation Analysis Program (SNAP) developed to perform covariance analysis on spacecraft inertial-measurement-unit (IMU) navigation errors are described and demonstrated. Consideration is given to the IMU modeling subroutine (with user-specified sensor characteristics), the data input procedures, state updates and the simulation of instrument failures, the determination of the nominal trajectory, the mapping-matrix and Monte Carlo covariance-matrix propagation methods, and aided-navigation simulation. Numerical results are presented in tables for sample applications involving (1) the Galileo/IUS spacecraft from its deployment from the Space Shuttle to a point 10 to the 8th ft from the center of the earth and (2) the TDRS-C/IUS spacecraft from Space Shuttle liftoff to a point about 2 h before IUS deployment. SNAP is shown to give reliable results for both cases, with good general agreement between the mapping-matrix and Monte Carlo predictions.

Weir, Kent A.↗
NASA NTRS2003

Monte Carlo Simulations and Generation of the SPI Response

In this paper we discuss the methods developed for the production of the INTEGRAL/SPI instrument response. The response files were produced using a suite of Monte Carlo simulation software developed at NASA/GSFC based on the GEANT-3 package available from CERN. The production of the INTEGRAL/SPI instrument response also required the development of a detailed computer mass model for SPI. We discuss our extensive investigations into methods to reduce both the computation time and storage requirements for the SPI response. We also discuss corrections to the simulated response based on our comparison of ground and inflight calibration data with MGEANT simulation.

Sturner, S. J.↗
NASA NTRS2003

Monte Carlo Simulations and Generation of the SPI Response

In this paper we discuss the methods developed for the production of the INTEGRAL/SPI instrument response. The response files were produced using a suite of Monte Carlo simulation software developed at NASA/GSFC based on the GEANT-3 package available from CERN. The production of the INTEGRAL/SPI instrument response also required the development of a detailed computer mass model for SPI. We discuss ow extensive investigations into methods to reduce both the computation time and storage requirements for the SPI response. We also discuss corrections to the simulated response based on our comparison of ground and infiight Calibration data with MGEANT simulations.

Sturner, S. J.↗
NASA NTRS1994

Transversely accelerated ions in the topside ionosphere

Data from the rocket campaigns Mechanism in the Auroral Region for Ion Energization (MARIE) and TOpside Probe of the Auroral Zone (TOPAZ) III, within regions of low-altitude transversely accelerated ions, are interpreted to explain the acceleration of the ions. Using the Monte Carlo kinetic technique to evaluate the ion heating produced by the simultaneously observed lower hybrid waves, we find that their observed electric field amplitudes are sufficient to explain the observed ion energies in the MARIE event. Much of the uncertainty in evaluating the efficiency of a plasma wave induced particle heating process which is dependent on a velocity resonance comes from the lack of information on the phase velocities of the waves. In the case of the MARIE observations, our modeling efforts show that features in the ion velocity distribution are consistent with the wave phase velocities inferred from interferometer measurements of wavelengths. The lower hybrid waves with which low-altitude transversely accelerated ions are associated are frequently observed to be concentrated in small-scale wave packets called 'spikelets'. We demonstrate through the scaling of the size of these wave packets that they are consistent with the theory of lower hybrid collapse. Using the Monte Carlo technique, we find that if the lower hybrid field energy is concentrated in these wave packets, it is still adequate to accelerate the ionospheric ions to the observed energies.

Retterer, John M.↗
NASA NTRS1998

LSPRAY: Lagrangian Spray Solver for Applications With Parallel Computing and Unstructured Gas-Phase Flow Solvers

Sprays occur in a wide variety of industrial and power applications and in the processing of materials. A liquid spray is a phase flow with a gas as the continuous phase and a liquid as the dispersed phase (in the form of droplets or ligaments). Interactions between the two phases, which are coupled through exchanges of mass, momentum, and energy, can occur in different ways at different times and locations involving various thermal, mass, and fluid dynamic factors. An understanding of the flow, combustion, and thermal properties of a rapidly vaporizing spray requires careful modeling of the rate-controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as other phenomena. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on sprays to unstructured grids and parallel computing. LSPRAY, which was developed by M.S. Raju of Nyma, Inc., is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo probability density function (PDF) solver. The LSPRAY solver accommodates the use of an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements in the gas-phase solvers. It is used specifically for fuel sprays within gas turbine combustors, but it has many other uses. The spray model used in LSPRAY provided favorable results when applied to stratified-charge rotary combustion (Wankel) engines and several other confined and unconfined spray flames. The source code will be available with the National Combustion Code (NCC) as a complete package.

Raju, Manthena S.↗
NASA NTRS1993

Comparison of remote sensing measurements with a two-scale polarimetric emission and scattering model for sea surfaces

Recently, it has been observed that the brightness temperatures of sea surfaces correlate with the azimuth angle of the ocean wind vector, indicating that ocean wind direction can not only be retrieved from the microwave backscatter, but also from the brightness temperature measurements. In this paper, comparison of the theoretical result calculated from a two-scale emission and scattering model with the Seasat and SSM/I model functions and aircraft measurements is presented and potential applications of microwave polarimetry are discussed. In our two-scale model, the modified reflectivities of large scale surfaces are calculated by extending the small perturbation method to the second order for small scale perturbation with anisotropic directional spectrum. It was found that the modified reflectivities derived from the second-order scattered field agree excellently well with the results obtained from a Monte Carlo simulation technique which numerically calculates the polarimetric reflectivities of one-dimensional random rough surfaces with a power-law spectrum. Without the second-order correction, the modified reflectivities of the rough surfaces are significantly over-estimated and sign errors are observed in the third and fourth Stokes parameters for thermal emissions. The surface spectrum parameters and two-scale cutoff are selected so that the calculated scattering and emission signatures agree with the reported model functions for Seasat and SSM/I. Subsequently, the polarimetric signatures of sea surfaces are illustrated to indicate the possibility of reducing the number of azimuthal looks required for spaceborne sensors in the remote sensing of ocean wind by using the polarimetric information. Furthermore, it is found that contrary to the dependence of backscattering coefficients on incidence angles, polarimetric brightness temperatures display a stronger wind direction dependence in the near nadir-looking direction than away from nadir. Finally, we discuss the potential of a single-look polarimetric radiometer for wind retrieval.

Yueh, S. H.↗
NASA NTRS2018

Aircraft System Noise Prediction Uncertainty Quantification for a Hybrid Wing Body Subsonic Transport Concept

Aircraft system level noise prediction for advanced, unconventional concepts has undergone significant improvement over the past two decades. The prediction modeling uncertainty must be quantified so that potential benefits of unconventional configurations, which are outside of the range of empirical models, can be reliably assessed. This paper builds on previous work in an effort to improve estimates of element prediction uncertainties where the prediction methodology has been improved, or new experimental validation data are available, to provide an estimate of the system level uncertainty in the prediction process. In general, the uncertainty of the prediction will be strongly dependent on the aircraft configuration as well as which technologies are integrated. While the quantitative uncertainty values contained here are specific to the hybrid wing body design presented, the underlying process is the same regardless of configuration. A refined process for determining the uncertainty for each element of the noise prediction is detailed in this paper. The system level uncertainty in the prediction of the aircraft noise is determined at the three certification points, using a Monte Carlo method. Comparisons with previous work show a reduction of 1 EPNdB in the 95%coverage interval of the cumulative noise level. The largest impediment for continued reduction in uncertainty for the hybrid wing body concept is the need for improved modeling and validation experiments for fan noise, propulsion airframe aeroacoustic effects, and the Krueger flap, which comprise the bulk of the uncertainty in the cumulative certification noise level.

June, Jason 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 NTRS1972

Monte Carlo calculations of high energy nucleon meson cascades and applications to galactic cosmic ray transport

Results obtained using a recently developed calculational method for determining the nucleon-meson cascade induced in thick materials by high-energy nucleons and charged pions are presented. The calculational method uses the intranuclear-cascade-evaporation model to treat nonelastic collisions by particles with energies approximately or smaller than GeV and an extrapolation model at higher energies. The following configurations are considered: (1) 19.2-GeV/c protons incident on iron; (2) 30.3-GeV/c protons incident on iron; (3) solar and galactic protons incident on the moon, and (4) galactic protons incident on tissue. For the first three configurations, experimental results are available and comparisons between the experimental and calculated results are given.

Armstrong, T. W.↗
NASA NTRS2008

Effect of the Inhomogeneity of Ice Crystals on Retrieving Ice Cloud Optical Thickness and Effective Particle Size

Spherical or spheroidal air bubbles are generally trapped in the formation of rapidly growing ice crystals. In this study the single-scattering properties of inhomogeneous ice crystals containing air bubbles are investigated. Specifically, a computational model based on an improved geometric-optics method (IGOM) has been developed to simulate the scattering of light by randomly oriented hexagonal ice crystals containing spherical or spheroidal air bubbles. A combination of the ray-tracing technique and the Monte Carlo method is used. The effect of the air bubbles within ice crystals is to smooth the phase functions, diminish the 22deg and 46deg halo peaks, and substantially reduce the backscatter relative to bubble-free particles. These features vary with the number, sizes, locations and shapes of the air bubbles within ice crystals. Moreover, the asymmetry factors of inhomogeneous ice crystals decrease as the volume of air bubbles increases. Cloud reflectance lookup tables were generated at wavelengths 0.65 m and 2.13 m with different air-bubble conditions to examine the impact of the bubbles on retrieving ice cloud optical thickness and effective particle size. The reflectances simulated for inhomogeneous ice crystals are slightly larger than those computed for homogenous ice crystals at a wavelength of 0.65 microns. Thus, the retrieved cloud optical thicknesses are reduced by employing inhomogeneous ice cloud models. At a wavelength of 2.13 microns, including air bubbles in ice cloud models may also increase the reflectance. This effect implies that the retrieved effective particle sizes for inhomogeneous ice crystals are larger than those retrieved for homogeneous ice crystals, particularly, in the case of large air bubbles.

Xie, Yu↗