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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Using Black Hole Mergers to Explore Structure Formation

Observations of gravitational waves from massive black hole mergers will open a new window into the era of structure formation in the early universe. Past efforts have concentrated on calculating merger rates using different physical assumptions, resulting in merger rate estimates that span a wide range (0.1 - 10(exp 4) mergers/year). We develop a semi-analytical, phenomenological model of massive black hole mergers that includes plausible combinations of several physical parameters, which we then turn around to determine how well observations with the Laser Interferometer Space Antenna (LISA) will be able to enhance our understanding of the universe during the critical z approximately equal to 5-30 epoch. Our approach involves generating synthetic LISA observable data (total BH masses, BH mass ratios, redshifts, merger rates), which are then analyzed using a Markov Chain Monte Carlo method, thus finding constraints for the physical parameters of the mergers. We find that our method works well at estimating merger parameters and that the number of merger events is a key discriminant among models, therefore making our method robust against observational uncertainties. Our approach can also be extended to more physically-driven models and more general problems in cosmology. This work is supported in part by the Cooperative Education Program at NASA/GSFC.

Alicea-Munoz, E.↗
NASA NTRS2002

The response of tissue-equivalent proportional counters to heavy ions

The paper presents a theoretical model for the response of a tissue-equivalent proportional counter (TEPC) irradiated with charged particles. Heavy ions and iron ions in particular constitute a significant part of radiation in space. TEPCs are used for all space shuttle and International Space Station (ISS) missions to estimate the dose and radiation quality (in terms of lineal energy) inside spacecraft. The response of the tissue-equivalent proportional counters shows distortions at the wall/cavity interface. In this paper, we present microdosimetric investigation using Monte Carlo track structure calculations to simulate the response of a TEPC to charged particles of various LET (1 MeV protons, 2.4 MeV alpha particles, 46 MeV/nucleon 20Ne, 55 MeV/nucleon 20Ne, 45 MeV/nucleon 40Ar, and 1.05 GeV/nucleon 56Fe). Data are presented for energy lost and energy absorbed in the counter cavity and wall. The model calculations are in good agreement with the results of Rademacher et al. (Radiat. Res. 149, 387-389, 1998), including the study of the interface between the wall and the sensitive region of the counter. It is shown that the anomalous response observed at large event sizes in the experiment is due to an enhanced entry of secondary electrons from the wall into the gas cavity.

NASA Discipline Radiation Health↗
NASA NTRS1972

Degree and direction of polarization of multiple scattered light. I - Homogeneous cloud layers.

It is shown that the Monte Carlo method can provide useful information on the polarization and its direction for homogeneous layers corresponding to the haze C, nimbostratus, and ice crystals models. Results for various optical thicknesses and two solar zenith angles show the variation with these parameters. In most cases the direction of polarization with respect to the meridian plane (which contains the outgoing photon) does not change appreciably with the optical thickness of the scattering layer, so that the conclusion that the direction of polarization is determined mainly by the direction of the incident and scattered photons, obtained for single scattering, should apply approximately when multiple scattering is taken into account.

Kattawar, G. W.↗
NASA NTRS1985

Tolerance bounds for log gamma regression models

The present procedure for finding lower confidence bounds for the quantiles of Weibull populations, on the basis of the solution of a quadratic equation, is more accurate than current Monte Carlo tables and extends to any location-scale family. It is shown that this method is accurate for all members of the log gamma(K) family, where K = 1/2 to infinity, and works well for censored data, while also extending to regression data. An even more accurate procedure involving an approximation to the Lawless (1982) conditional procedure, with numerical integrations whose tables are independent of the data, is also presented. These methods are applied to the case of failure strengths of ceramic specimens from each of three billets of Si3N4, which have undergone flexural strength testing.

Jones, R. A.↗
NASA NTRS1986

Persistent time intervals between features in solar flare hard X-ray emission

Several solar hard X-ray events (greater than 100 keV) were observed simultaneously with identical instruments on the Venera 11, 12, 13, 14, and Prognoz spacecraft. High time resolution (= 2 ms) data were stored in memory when a trigger occurred. The observations of modulation are presented with a period of 1.6 s for the event on December 3, 1978. Evidence is also presented for fast time fluctuations from an event on November 6, 1979, observed from Venera 12 and another on September 6, 1981, observed from the Solar Maximum Mission. Power spectrum analysis, epoch folding, and Monte Carlo simulation were used to evaluate the statistical significance of persistent time delays between features. The results are discussed in light of the MHD model proposed by Zaitsev and Stepanov.

Desai, Upendra D.↗
NASA NTRS1988

Confusion-limited galaxy fields. I - Simulated optical and near-infrared images

Techniques for simulating images of galaxy fields are presented that extend to high redshifts and a surface density of galaxies high enough to produce overlapping images. The observed properties of galaxies and galaxy-ensembles in the 'local' universe are extrapolated to high redshifts using reasonable scenarios for the evolution of galaxies and their spatial distribution. This theoretical framework is then employed with Monte Carlo techniques to create fairly realistic two-dimensional distributions of galaxies plus optical and near-infrared sky images in a variety of model universes, using the appropriate density, luminosity, and angular size versus redshift relations.

Chokshi, Arati↗
NASA NTRS1991

New approach in direct-simulation of gas mixtures

Results are reported for an investigation of a new direct-simulation Monte Carlo method by which energy transfer and chemical reactions are calculated. The new method, which reduces to the variable cross-section hard sphere model as a special case, allows different viscosity-temperature exponents for each species in a gas mixture when combined with a modified Larsen-Borgnakke phenomenological model. This removes the most serious limitation of the usefulness of the model for engineering simulations. The necessary kinetic theory for the application of the new method to mixtures of monatomic or polyatomic gases is presented, including gas mixtures involving chemical reactions. Calculations are made for the relaxation of a diatomic gas mixture, a plane shock wave in a gas mixture, and a chemically reacting gas flow along the stagnation streamline in front of a hypersonic vehicle. Calculated results show that the introduction of different molecular interactions for each species in a gas mixture produces significant differences in comparison with a common molecular interaction for all species in the mixture. This effect should not be neglected for accurate DSMC simulations in an engineering context.

Chung, Chan-Hong↗
NASA NTRS1994

TDRS orbit determination by radio interferometry

In support of a NASA study on the application of radio interferometry to satellite orbit determination, MITRE developed a simulation tool for assessing interferometry tracking accuracy. The Orbit Determination Accuracy Estimator (ODAE) models the general batch maximum likelihood orbit determination algorithms of the Goddard Trajectory Determination System (GTDS) with the group and phase delay measurements from radio interferometry. ODAE models the statistical properties of tracking error sources, including inherent observable imprecision, atmospheric delays, clock offsets, station location uncertainty, and measurement biases, and through Monte Carlo simulation, ODAE calculates the statistical properties of errors in the predicted satellites state vector. This paper presents results from ODAE application to orbit determination of the Tracking and Data Relay Satellite (TDRS) by radio interferometry. Conclusions about optimal ground station locations for interferometric tracking of TDRS are presented, along with a discussion of operational advantages of radio interferometry.

Pavloff, Michael S.↗
NASA NTRS2004

3-D Radiative Transfer in Clouds

The development of the subject of 3D radiative transfer in clouds will be reviewed, starting from the 1970s when the first papers using simple cubic clouds were published. The two main lines of development, Monte Carlo and numerical-analytical, will be reviewed. The theme will be the extraction of simple results from all the complexity of 3D radiation modeling.

Wiscombe, Warren↗
NASA NTRS2004

Atomistic Modeling of Pd Site Preference in NiTi

An analysis of the site subsitution behavior of Pd in NiTi was performed using the BFS method for alloys. Through a combination of Monte Carlo simulations and detailed atom-by-atom energetic analyses of various computational cells, representing compositions of NiTi with up to 10 at% Pd, a detailed understanding of site occupancy of Pd in NiTi was revealed. Pd subsituted at the expense of Ni in a NiTi alloy will prefer the Ni-sites. Pd subsituted at the expense of Ti shows a very weak preference for Ti-sites that diminishes as the amount of Pd in the alloy increases and as the temperature increases.

Bozzolo, Guillermo↗
NASA NTRS2023

The Impact and Estimation of Uncertainty Correlation for Multi-Angle Polarimetric Remote Sensing of Aerosols and Ocean Color

Multi-angle polarimetric (MAP) measurements contain rich information for characterization of aerosol microphysical and optical properties that can be used to improve atmospheric correction in ocean color remote sensing. Advanced retrieval algorithms have been developed to obtain multiple geophysical parameters in the atmosphere-ocean system, although uncertainty correlation among measurements is generally ignored due to lack of knowledge on its strength and characterization. In this work, we provide a practical framework to evaluate the impact of the angular uncertainty correlation from retrieval results and a method to estimate correlation strength from retrieval fitting residuals. The Fast Multi-Angular Polarimetric Ocean coLor (FastMAPOL) retrieval algorithm, based on neural network forward models, is used to conduct the retrievals and uncertainty quantification. In addition, we also discuss a flexible approach to include a correlated uncertainty model in the retrieval algorithm. The impact of angular correlation on retrieval uncertainties is discussed based on synthetic AirHARP and HARP2 measurements using a Monte Carlo uncertainty estimation method. Correlation properties are estimated using auto-correlation functions based on the fitting residuals from both synthetic AirHARP and HARP2 data and real AirHARP measurement, with the resulting angular correlation parameters found to be larger than 0.9 and 0.8 for reflectance and DoLP, respectively, which correspond to correlation angles of 10° and 5°. Although this study focuses on angular correlation from HARP instruments, the methodology to study and quantify uncertainty correlation is also applicable to other instruments with angular, spectral, or spatial correlations, and can help inform laboratory calibration and characterization of the instrument uncertainty structure.

PACE↗
NASA NTRS1993

Parallelizing Timed Petri Net simulations

The possibility of using parallel processing to accelerate the simulation of Timed Petri Nets (TPN's) was studied. It was recognized that complex system development tools often transform system descriptions into TPN's or TPN-like models, which are then simulated to obtain information about system behavior. Viewed this way, it was important that the parallelization of TPN's be as automatic as possible, to admit the possibility of the parallelization being embedded in the system design tool. Later years of the grant were devoted to examining the problem of joint performance and reliability analysis, to explore whether both types of analysis could be accomplished within a single framework. In this final report, the results of our studies are summarized. We believe that the problem of parallelizing TPN's automatically for MIMD architectures has been almost completely solved for a large and important class of problems. Our initial investigations into joint performance/reliability analysis are two-fold; it was shown that Monte Carlo simulation, with importance sampling, offers promise of joint analysis in the context of a single tool, and methods for the parallel simulation of general Continuous Time Markov Chains, a model framework within which joint performance/reliability models can be cast, were developed. However, very much more work is needed to determine the scope and generality of these approaches. The results obtained in our two studies, future directions for this type of work, and a list of publications are included.

Nicol, David M.↗
NASA NTRS1993

Polarimetric thermal emission from rough surfaces

Recent theoretical works have suggested the potential of passive polarimetry in the remote sensing of geophysical media. It was shown that the third Stokes parameter U of the thermal emission may become larger for azimuthally asymmetric fields of observation. In order to investigate the potential applicability of passive polarimetry to the remote sensing of ocean surface, measurements of the polarimetric thermal emission from a sinusoidal water surface and a numerical study of the polarimetric thermal emission from randomly rough ocean surfaces were performed. Measurements of sinusoidal water surface thermal emission were performed using a sinusoidal water surface which was created by placing a thin sheet of fiberglass with a sinusoidal profile in two dimensions extended infinitely in the third dimension onto a water surface. The theory of thermal emission from a 'two-layer' periodic surface is derived and the exact solution is performed using both the extended boundary condition method (EBC) and the method of moments (MOM). The theoretical predictions are found to be in good agreement with the experimental results once the effects of the radiometer antenna pattern are included and the contribution of background noise to the measurements is modeled. The experimental results show that the U parameter indicates the direction of periodicity of the water surface and can approach values of up to 30 K for the surface observed. Next, a numerical study of polarimetric thermal emission from randomly rough surfaces was performed. A Monte Carlo technique utilizing an exact method for calculating thermal emission was chosen for the study to avoid any of the limitations of the commonly used approximate methods in rough surface scattering. In this Monte Carlo technique, a set of finite rough surface profiles in two dimensions with desired statistics was generated and extended periodically. The polarimetric thermal emission from each surface of the set was then calculated using both the EBC and the MOM and the results were averaged. The surface statistics chosen were intended to model a wind perturbed ocean surface in the X to Ku band microwave region. The results indicate that the U parameter is sensitive to the azimuthal angle between the surface periodicity and the looking angle and to the rms height of the surface, and that the U parameter is fairly insensitive to variations in polar angle, permittivity, surface power law spectrum, and surface spectrum high frequency cutoff. These properties give further strength to the idea of using the U parameter to detect wind direction over the ocean.

Johnson, J. T.↗
NASA NTRSDate not supplied

Modeling Radiation Sources for Experimental Validation

To build on previously completed experimental work, I have been working to implement a model of our lab’s Americium-Beryllium neutron source and Ludlum detection instruments computationally to provide secondary verification of our experimental results. This has been done using the radiation transport code MCNP6 (Monte Carlo N-Particle). I have developed files to calculate the mass ratio of various composite materials that we have exposed to the radiation source, as well as the components of the detection instruments and converted all of these into geometric coordinates that define the experimental setup. Identifying the proper energy distribution and detector efficiencies as well as the correct interpretation of the detector physics have proven to be the most difficult and time-consuming challenges of the work. It is currently too early for results based on the most updated and accurate models. With this experience, the next step is to develop another source model of a proton beam to simulate exposures on other composite targets. This will be developed based on existing experimental datasets, but not used to validate them. Each of these projects will further define the radiation shielding needs for humans in space radiation environments that astronauts will someday be exposed to beyond Low Earth Orbit.

Juliana Simon↗
NASA NTRS2014

James Webb Space Telescope (JWST) Stationkeeping Monte Carlo Simulations

The James Webb Space Telescope (JWST) will launch in 2018 into a Libration Point Orbit (LPO) around the Sun-EarthMoon (SEM) L2 point, with a planned mission lifetime of 11 years. This paper discusses our approach to Stationkeeping (SK) maneuver planning to determine an adequate SK delta-V budget. The SK maneuver planning for JWST is made challenging by two factors: JWST has a large Sunshield, and JWST will be repointed regularly producing significant changes in Solar Radiation Pressure (SRP). To accurately model SRP we employ the Solar Pressure and Drag (SPAD) tool, which uses ray tracing to accurately compute SRP force as a function of attitude. As an additional challenge, the future JWST observation schedule will not be known at the time of SK maneuver planning. Thus there will be significant variation in SRP between SK maneuvers, and the future variation in SRP is unknown. We have enhanced an earlier SK simulation to create a Monte Carlo simulation that incorporates random draws for uncertainties that affect the budget, including random draws of the observation schedule. Each SK maneuver is planned to optimize delta-V magnitude, subject to constraints on spacecraft pointing. We report the results of the Monte Carlo simulations and discuss possible improvements during flight operations to reduce the SK delta-V budget.

Monte Carlo↗
NASA NTRSDate not supplied

Aerial Vehicle Routing and Scheduling for UAS Traffic Management: A Hybrid Monte Carlo Tree Search Approach

We present the Multi-Route Weighted Package Delivery Problem (MRWPDP) and a scalable solution methodology as a major step towards enabling an airspace deconfliction service for drone delivery operations. The problem is motivated by Strategic deconfliction under the FAA’s “Unmanned Aircraft Systems Traffic Management” Concept of Operations. MRWPDP falls under a class of vehicle routing and scheduling problems, and as such is NP-Hard. In MRWPDP, a graph network is given which consists of depots, drop-off sites, and multiple routes connecting the two. In addition, routes are weighted by the associated ground risk and total travel distance for package delivery. The goal is to optimally schedule the departure time and assign routes to a known set of vehicles at the depot. We propose a heuristic solution to the problem by borrowing techniques from Mixed Integer Linear Programming (MILP), Constraint Programming, and Monte Carlo Tree Search (MCTS). The resulting hybrid framework is MCTS with Bound-and-Prune (BP) and rapid simulated updates (U), or MCTS-BP-U. This approach is able to quickly provide a feasible solution for MRWPDP, even for large problem instances up to 1000 vehicles. We provide a MILP formulation of MRWPDP and compare its performance against MCTS-BP-U in terms of solution quality. An agent-based model simulation is conducted as a final step to validate the efficacy of our approach.

air traffic scheduling↗
NASA NTRSDate not supplied

Aerial Vehicle Routing and Scheduling for UAS Traffic Management: A Hybrid Monte Carlo Tree Search Approach

We present the Multi-Route Weighted Package Delivery Problem (MRWPDP) and a scalable solution methodology as a major step towards enabling an airspace deconfliction service for drone delivery operations. The problem is motivated by Strategic deconfliction under the FAA’s “Unmanned Aircraft Systems Traffic Management” Concept of Operations. MRWPDP falls under a class of vehicle routing and scheduling problems, and as such is NP-Hard. In MRWPDP, a graph network is given which consists of depots, drop-off sites, and multiple routes connecting the two. In addition, routes are weighted by the associated ground risk and total travel distance for package delivery. The goal is to optimally schedule the departure time and assign routes to a known set of vehicles at the depot. We propose a heuristic solution to the problem by borrowing techniques from Mixed Integer Linear Programming (MILP), Constraint Programming, and Monte Carlo Tree Search (MCTS). The resulting hybrid framework is MCTS with Bound-and-Prune (BP) and rapid simulated updates (U), or MCTS-BP-U. This approach is able to quickly provide a feasible solution for MRWPDP, even for large problem instances up to 1000 vehicles. We provide a MILP formulation of MRWPDP and compare its performance against MCTS-BP-U in terms of solution quality. An agent-based model simulation is conducted as a final step to validate the efficacy of our approach.

air traffic scheduling↗
NASA NTRS1987

A dust scattered halo in starburst galaxy M82?

The source of the halo about M82 has been under discussion for several years. One explanation for it is the dust model of Solinger, Morrison and Markert in which they propose a diffuse cloud of dust through the M81 group, with M82 traveling through the group holding a denser cloud of dust around it. The feasibility of the dust theory is examined in the X-ray range, using the halo in the X-ray image of M82 taken by the Einstein Observatory. To this end the X-ray cross section for dust is presented, along with the single scattered image of an X-ray source surrounded by a dust cloud; multiply scattered images were simulated with a Monte Carlo program; profiles of the halo along the major and minor axes of M82 are presented. Also presented is an accounting for line spectrographs of M82 that show unusual splitting, using the dust model. The final model proposed for the X-ray image requires dust of radius 50 to 300 A, with density on the order of 10 to the -7th power cu cm, out to a distance of about 9 kpc for some regions.

Rohan, Michael↗