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 973 records · Page 54

Monte Carlo simulation of EAS generated by 10(14) - 10(16) eV protons

Detailed Monte Carlo simulations of extensive air showers to be detected by the Homestake Surface Underground Telescope and other similar detectors located at sea level and mountain altitudes have been performed for 10 to the 14th power to 10 to the 16th power eV primary energies. The results of these Monte Carlo calculations will provide an opportunity to compare the experimental data with different models for the composition and spectra of primaries and for the development of air showers. The results obtained for extensive air showers generated by 10 to the 14th power to 10 to the 16th power eV primary protons are reported.

Fenyves, E. J.↗

A potential energy scaling Monte Carlo simulation of thin film nucleation and growth

The initial growth of thin Ge fims on the (100) surface of an Fe substrate is investigated theoretically by means of Monte Carlo simulations based on a potential-energy-scaling technique. The substrate is modeled as a 20 x 20-square array with periodic boundary conditions, as described by Heinbockel et al. (1983), and the movement of surface atoms under the influence of the substrate interaction potential and the lateral interaction of neighboring atoms is explored via continuous updating (on the time scale of single events) of the potential energy at each site in the array. Results for the clustering of nine dispersed atoms over 1.0 s at 600 K and for deposition at 5 x 10 to the -14th/sq cm s over 2.0 s at 500 K are presented graphically.

Outlaw, R. A.↗

Probabilistic low cycle fatigue failure analysis with application to liquid propellant rocket engines

A probabilistic Low Cycle Fatigue (LCF) failure analysis of a candidate turbine disk for use in a turbopump of a rocket engine of the Space Shuttle Main Engine class is described. A state-of-the-art LCF failure prediction method was used in a Monte Carlo simulation to generate a distribution of failure lives. A stochastic Stress/Life (S/N) model was used for materials characterization. The LCF failure model expresses fatigue life as a function of stochastic parameters including environmental parameters, loads, material properties, structural parameters, and model specification errors. The rationale for the particular characterization of each stochastic input parameter is described. The results and interpretation of the failure analysis are given.

Newlin, L.↗

Ion and relativistic electron transport in solar flares

Recent models for the ion and relativistic electron transport in solar flare magnetic loops are reviewed, focusing on the ways in which the models handle loop structure, particle acceleration, and particle transport. Results are presented from Monte Carlo simulations of ion and relativistic ion transport. Consideration is given to results concerning the depth distributions of gamma-ray production, the attenuation of gamma-ray lines, the time dependences of the various emissions and the angular distribution of the bremsstrahlung.

Ramaty, R.↗

Hypersonic flow - Physical issues and computational simulation

Computational results are presented on hypersonic flows, with emphasis on thermal and chemical nonequilibrium effects and the role played by grids. The continuum equations and their boundary conditions are presented, and the direct simulation Monte Carlo method is described. The importance of using the correct computational and physical models in describing hypersonic nonequilibrium flows is demonstrated.

Hassan, H. A.↗

Scalar dissipation and mixing in turbulent reacting flows

The mixing term that appears in probability density function-based formulations of the statistical behavior of turbulent reacting flows was studied with the use of direct numerical simulations of an irreversible, second-order, isothermal chemical reaction between initially segregated reactants in decaying, homogeneous turbulence. Three-dimensional, time-dependent pseudospectral calculations with 64-cubed Fourier modes were made as well as were Monte Carlo simulations of the joint concentration pdf equations for coalescence-dispersion (C/D) models. It was found that the development of the scalar microscales depends largely on the initial concentration distributions (both isotropic and nonisotropic distributions were used) and is relatively independent of the reaction rate constant. Joint concentration probability density functions were compared to predictions of D/D models that make use of parameters evaluated from the nonreacting or inert system.

Leonard, Andy D.↗

A probabilistic fracture mechanics approach for structural reliability assessment of space flight systems

A probabilistic fracture mechanics approach for predicting the failure life distribution due to subcritical crack growth is presented. A state-of-the-art crack propagation method is used in a Monte Carlo simulation to generate a distribution of failure lives. The crack growth failure model expresses failure life as a function of stochastic parameters including environment, loads, material properties, geometry, and model specification errors. A stochastic crack growth rate model that considers the uncertainties due to scatter in the data and mode misspecification is proposed. The rationale for choosing a particular type of probability distribution for each stochastic input parameter and for specifying the distribution parameters is presented. The approach is demonstrated through a probabilistic crack growth failure analysis of a welded tube in the Space Shuttle Main Engine. A discussion of the results from this application of the methodology is given.

Sutharshana, S.↗

Oriented Scintillation Spectrometer Experiment observations of Co-57 in SN 1987A

The Oriented Scintillation Spectrometer Experiment (OSSE) on the Compton Gamma Ray Observatory has observed SN 1987A for two 2 week periods during the first 9 months of the mission. Evidence for gamma-ray line and continuum emission from Co-57 is observed with an intensity of about 10 exp -4 gamma/sq cm/s. This photon flux between 50 and 136 keV is demonstrated by Monte Carlo calculations to be independent of the radial distribution of Co-57 for models of low optical depth, viz., models having photoelectric absorption losses of 122 keV photons no greater than several percent. For such models the observed Co-57 flux indicates that the ratio Ni-57/Ni-56 produced in the explosion was about 1.5 times the solar system ratio of Fe-57/Fe-56. When compared with nearly contemporaneous bolometric estimates of the luminosity for SN 1987A, our observations imply that Co-57 radioactivity does not account for most of the current luminosity of the supernova remnant in low optical depth models. We suggest alternatives, including a large optical depth model that is able to provide the SN 1987A luminosity and is consistent with the OSSE flux. It requires a 57/56 production ratio about twice solar.

Kurfess, J. D.↗

Spacecraft contamination investigation by direct simulation Monte Carlo - Contamination on UARS/HALOE

A three-dimensional version of the direct simulation Monte Carlo method is adapted to study the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite, with the ultimate goal of estimating the cumulative level of molecular deposits on critical optical surfaces of the Halogen Occultation Experiment. Emphasis is placed on describing the development of preprocessing codes that simplify the incorporation of geometry and outgassing data into the main run code. A novel use of parallel processing software allows use of multiple processors to solve different regions of the overall flowfield. Initial results relating to contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the Halogen Occultation telescope's aperture plane are presented, along with data related to code performance.

Rault, Didier F. G.↗

On predicting contamination levels of HALOE optics aboard UARS using direct simulation Monte Carlo

A three-dimensional version of the direct simulation Monte Carlo method is adapted to assess the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite. Emphasis is placed on simulating a realistic, worst-case set of flowfield and surface conditions and geometric orientations in order to estimate an upper limit for the cumulative level of volatile organic molecular deposits at the aperture of the Halogen Occultation Experiment. Problems resolving species outgassing and vent flux rates that varied over many orders of magnitude were handled using species weighting factors. Results relating to contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the target surface are presented, along with data related to code performance. Using procedures developed in standard contamination analyses, the cumulative level of volatile organic deposits on HALOE's aperture over the instrument's 35-month nominal data collection period is estimated to be about 2700A.

Woronowicz, Michael S.↗

Design optimization and probabilistic analysis of a hydrodynamic journal bearing

A nonlinear constrained optimization of a hydrodynamic bearing was performed yielding three main variables: radial clearance, bearing length to diameter ratio, and lubricating oil viscosity. As an objective function a combined model of temperature rise and oil supply has been adopted. The optimized model of the bearing has been simulated for population of 1000 cases using Monte Carlo statistical method. It appeared that the so called 'optimal solution' generated more than 50 percent of failed bearings, because their minimum oil film thickness violated stipulated minimum constraint value. As a remedy change of oil viscosity is suggested after several sensitivities of variables have been investigated.

Liniecki, Alexander G.↗

Direct simulation Monte Carlo prediction of on-orbit contaminant deposit levels for HALOE

A three-dimensional version of the direct simulation Monte Carlo method is adapted to assess the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite. Emphasis is placed on simulating a realistic, worst-case set of flow field and surface conditions and geometric orientations for the satellite in order to estimate an upper limit for the cumulative level of volatile organic molecular deposits at the aperture of the Halogen Occultation Experiment. A detailed description of the adaptation of this solution method to the study of the satellite's environment is also presented. Results pertaining to the satellite's environment are presented regarding contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the target surface, along with data related to code performance. Using procedures developed in standard contamination analyses, along with many worst-case assumptions, the cumulative upper-limit level of volatile organic deposits on HALOE's aperture over the instrument's 35-month nominal data collection period is estimated at about 13,350 A.

Woronowicz, Michael S.↗

Sensor Based Engine Life Calculation: A Probabilistic Perspective

It is generally known that an engine component will accumulate damage (life usage) during its lifetime of use in a harsh operating environment. The commonly used cycle count for engine component usage monitoring has an inherent range of uncertainty which can be overly costly or potentially less safe from an operational standpoint. With the advance of computer technology, engine operation modeling, and the understanding of damage accumulation physics, it is possible (and desirable) to use the available sensor information to make a more accurate assessment of engine component usage. This paper describes a probabilistic approach to quantify the effects of engine operating parameter uncertainties on the thermomechanical fatigue (TMF) life of a selected engine part. A closed-loop engine simulation with a TMF life model is used to calculate the life consumption of different mission cycles. A Monte Carlo simulation approach is used to generate the statistical life usage profile for different operating assumptions. The probabilities of failure of different operating conditions are compared to illustrate the importance of the engine component life calculation using sensor information. The results of this study clearly show that a sensor-based life cycle calculation can greatly reduce the risk of component failure as well as extend on-wing component life by avoiding unnecessary maintenance actions.

Guo, Ten-Huei↗

A Hierarchical Approach to Fracture Mechanics

Recent research conducted under NASA LaRC's Creativity and Innovation Program has led to the development of an initial approach for a hierarchical fracture mechanics. This methodology unites failure mechanisms occurring at different length scales and provides a framework for a physics-based theory of fracture. At the nanoscale, parametric molecular dynamic simulations are used to compute the energy associated with atomic level failure mechanisms. This information is used in a mesoscale percolation model of defect coalescence to obtain statistics of fracture paths and energies through Monte Carlo simulations. The mathematical structure of predicted crack paths is described using concepts of fractal geometry. The non-integer fractal dimension relates geometric and energy measures between meso- and macroscales. For illustration, a fractal-based continuum strain energy release rate is derived for inter- and transgranular fracture in polycrystalline metals.

Saether, Erik↗

A Deterministic Electron, Photon, Proton and Heavy Ion Radiation Transport Suite for the Study of the Jovian System

A deterministic suite of radiation transport codes, developed at NASA Langley Research Center (LaRC), which describe the transport of electrons, photons, protons, and heavy ions in condensed media is used to simulate exposures from spectral distributions typical of electrons, protons and carbon-oxygen-sulfur (C-O-S) trapped heavy ions in the Jovian radiation environment. The particle transport suite consists of a coupled electron and photon deterministic transport algorithm (CEPTRN) and a coupled light particle and heavy ion deterministic transport algorithm (HZETRN). The primary purpose for the development of the transport suite is to provide a means for the spacecraft design community to rapidly perform numerous repetitive calculations essential for electron, proton and heavy ion radiation exposure assessments in complex space structures. In this paper, the radiation environment of the Galilean satellite Europa is used as a representative boundary condition to show the capabilities of the transport suite. While the transport suite can directly access the output electron spectra of the Jovian environment as generated by the Jet Propulsion Laboratory (JPL) Galileo Interim Radiation Electron (GIRE) model of 2003; for the sake of relevance to the upcoming Europa Jupiter System Mission (EJSM), the 105 days at Europa mission fluence energy spectra provided by JPL is used to produce the corresponding dose-depth curve in silicon behind an aluminum shield of 100 mils ( 0.7 g/sq cm). The transport suite can also accept ray-traced thickness files from a computer-aided design (CAD) package and calculate the total ionizing dose (TID) at a specific target point. In that regard, using a low-fidelity CAD model of the Galileo probe, the transport suite was verified by comparing with Monte Carlo (MC) simulations for orbits JOI--J35 of the Galileo extended mission (1996-2001). For the upcoming EJSM mission with a potential launch date of 2020, the transport suite is used to compute the traditional aluminum-silicon dose-depth calculation as a standard shield-target combination output, as well as the shielding response of high charge (Z) shields such as tantalum (Ta). Finally, a shield optimization algorithm is used to guide the instrument designer with the choice of graded-Z shield analysis.

Norman, Ryan B.↗

Use of Polarization Lidar for Investigation of Meteorological Formations

This paper presents the results of theoretical and experimental, investigations of depolarization characteristics of different meteorological formations. Experimental investigations are carried out with a monostatic lidar. The ruby laser radiation is polarized in a vertical plane. The radiation reflected is accepted by a lens system of 150 rom in diameter and a viewing angle of 12' and further it is divided by Wollaston prism into the components polarized orthogonally. In this case the principal plane of the prism is exposed parallel with the laser polarization plane. Investigations show the degree of radiation polarization, reflected from water clouds, to be changed within 1/0.7 (seldom up to 0.6) depending on their density. In most cases a signal reflected from the cloud leading edge is polarized completely. The time shift is observed between polarized and crosspolarized components of a Fignal, reflected from a cloud, depending on the density of a meteorological object. While penetrating into the cloud depth a degree of polarization decreases up to 0.8-0.7, and the character of this decrease is different for various types of clouds. For crystal clouds the shift between the components of the reflected signal is not observed and the magnitude of polarization degree amounts to 0.1/0.3 in comparison with water clouds. The polarization degree of radiation reflected by fog is not less than 0.6, and that in the rains of average intensity (about 5 mm/h) is always about 1. The authors have suggested an algorithm of numerical solution of nonstationary transfer equation in the vector form to forecast the influence of multiple scattering effects on polarization characteristics of the lidar light signal. The method of statistical simulation (Monte-Carlo technique) forms the basis of the algorithm. Numerical estimates obtained for a model of stratocumulus at lambda = 0.6943 microns under boundary conditions close to the conditions of natural experiment being discussed proved to be in a good agreement with the results of observation. Specifically, Fig. 1 shows the profiles of polarization (p) versus depth (L) of the following drop formations: fog (curve 1) with horizontal meteorological visibility of 4 km two stratocumulus at a height of 1100 m with the attenuation factors delta = 0.01m(exp -1) (curve 2) and delta = 0.05m(exp -1) (curve 3). Curve 3 shows the results of numerical estimates and the value of their statistical error.

Balin, Yu. S.↗

Landing Gear Noise Prediction and Analysis for Tube-and-Wing and Hybrid-Wing-Body Aircraft

Improvements and extensions to landing gear noise prediction methods are developed. New features include installation effects such as reflection from the aircraft, gear truck angle effect, local flow calculation at the landing gear locations, gear size effect, and directivity for various gear designs. These new features have not only significantly improved the accuracy and robustness of the prediction tools, but also have enabled applications to unconventional aircraft designs and installations. Systematic validations of the improved prediction capability are then presented, including parametric validations in functional trends as well as validations in absolute amplitudes, covering a wide variety of landing gear designs, sizes, and testing conditions. The new method is then applied to selected concept aircraft configurations in the portfolio of the NASA Environmentally Responsible Aviation Project envisioned for the timeframe of 2025. The landing gear noise levels are on the order of 2 to 4 dB higher than previously reported predictions due to increased fidelity in accounting for installation effects and gear design details. With the new method, it is now possible to reveal and assess the unique noise characteristics of landing gear systems for each type of aircraft. To address the inevitable uncertainties in predictions of landing gear noise models for future aircraft, an uncertainty analysis is given, using the method of Monte Carlo simulation. The standard deviation of the uncertainty in predicting the absolute level of landing gear noise is quantified and determined to be 1.4 EPNL dB.

Guo, Yueping↗