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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Model scattering envelopes of young stellar objects. II - Infalling envelopes

We present scattered light images for models of young stellar objects surrounded by dusty envelopes. The envelopes are assumed to have finite angular momentum and are falling in steady flow onto a disk. The model envelopes include holes, such as might be created by energetic bipolar flows. We calculate images using the Monte Carlo method to follow the light scattered in the dusty envelope and circumstellar disk, assuming that the photons originate from the central source. Adopting typical interstellar medium dust opacities and expected mass infall rates for protostars of about 10 exp -6 solar mass/yr, we find that detectable amounts of optical radiation can escape from envelopes falling into a disk as small as about 10-100 AU, depending upon the viewing angle and the size of the bipolar flow cavity. We suggest that the extended optical and near-IR light observed around several young stars is scattered by dusty infalling envelopes rather than disks.

Whitney, Barbara A.↗

NASA NTRS1994

Computation of axisymmetric and ionized hypersonic flows using particle and continuum methods

Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo (DSMC) method, and the continuum approach solves the appropriate equations of fluid flow. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code that is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuum of a previous study that considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry model employed.

Boyd, Iain D.↗

NASA NTRS1996

Dust Around Herbig Ae Stars: Additional Constraints from their Photometric and Polarimetric Variability

For the Herbig Ae stars with Algol-like minima (UX Ori, WW Vul, etc), the effects of circumstellar dust include: excess infrared emission, anomalous ultraviolet extinction, the 'blueing' of the stars in minima accompanying by an increase of intrinsic polarization. Using a Monte-Carlo code for polarized radiation transfer we have simulated these effects and compared the results obtained for different models with the observational data available. We found that the photometric and polarimetric behavior of the stars provided essential additional constraints on the circumstellar dust models. The models with spheroidal shell geometry and compact (non-fluffy) dust grains do not appear to be able to explain all the data.

Krivova, N. A.↗

NASA NTRS2000

DSMC Simulations of Shock Interactions About Sharp Double Cones

This paper presents the results of a numerical study of shock interactions resulting from Mach 10 flow about sharp double cones. Computations are made by using the direct simulation Monte Carlo (DSMC) method of Bird. The sensitivity and characteristics of the interactions are examined by varying flow conditions, model size, and configuration. The range of conditions investigated includes those for which experiments have been or will be performed in the ONERA R5Ch low-density wind tunnel and the Calspan-University of Buffalo Research Center (CUBRC) Large Energy National Shock (LENS) tunnel.

Moss, James N.↗

NASA NTRS2002

Simulations of Ground and Space-Based Oxygen Atom Experiments

Fast, pulsed atomic oxygen sources are a key tool in ground-based investigations of spacecraft contamination and surface erosion effects. These technically challenging ground-based studies provide a before and after picture of materials under low-earth-orbit (LEO) conditions. It would be of great interest to track in real time the pulsed flux from the source to the surface sample target and beyond in order to characterize the population of atoms and molecules that actually impact the surface and those that make it downstream to any coincident detectors. We have performed simulations in order to provide such detailed descriptions of these ground-based measurements and to provide an assessment of their correspondence to the actual LEO environment. Where possible we also make comparisons to measured fluxes and erosion yields. To perform the calculations we use a detailed description of a measurement beam and surface geometry based on the W, pulsed apparatus at Montana State University. In this system, a short pulse (on the order of 10 microseconds) of an O/O2 beam impacts a flat sample about 40 cm downstream and slightly displaced &om the beam s central axis. Past this target, at the end of the beam axis is a quadrupole mass spectrometer that measures the relative in situ flux of 0102 to give an overall normalized erosion yield. In our simulations we use the Direct Simulation Monte Carlo (DSMC) method, and track individual atoms within the atomic oxygen pulse. DSMC techniques are typically used to model rarefied (few collision) gas-flows which occur at altitudes above approximately 110 kilometers. These techniques are well suited for the conditions here, and multi-collision effects that can only be treated by this or a similar technique are included. This simulation includes collisions with the surface and among gas atoms that have scattered from the surface. The simulation also includes descriptions of the velocity spread and spatial profiles of the O/O2 beam obtained from separate measurements. These computations use basic engineering models for the gas-gas and gas-surface scattering and focus on the influence of multi-collision effects. These simulations characterize many important quantities of interest including the actual flux of atoms that reach the surface, the energy distribution of this flux, as well as the direction of the velocity of the flux that strikes the surface. These quantities are important in characterizing the conditions which give rise to measured surface erosion. The calculations also yield time- snapshots of the pulse as it impacts and flows around the surface. These snapshots reveal the local environment of gas near the surface for the duration of the pulse. We are also able to compute the flux of molecules that travel downstream and reach the spectrometer, and we characterize their velocity distribution. The number of atoms that reach the spectrometer can in fact be influenced by the presence of the surface due to gas-gas collisions from atoms scattered h m the surface, and it will generally be less than that with the surface absent. This amounts to an overall normalization factor in computing erosion yields. We discuss these quantities and their relationship to the gas-surf$ce interaction parameters. We have also performed similar calculations corresponding to conditions (number densities, temperatures, and velocities) of low-earth orbit. The steady-state nature and lower overall flux of the actual space environment give rise to differences in the nature of the gas-impacts on the surface from those of the ground-based measurements using a pulsed source.

Minton, T. K.↗

NASA NTRS2003

"Hypothetical" Heavy Particles Dynamics in LES of Turbulent Dispersed Two-Phase Channel Flow

The extensive experimental study of dispersed two-phase turbulent flow in a vertical channel has been performed in Eaton's research group in the Mechanical Engineering Department at Stanford University. In Wang & Squires (1996), this study motivated the validation of LES approach with Lagrangian tracking of round particles governed by drag forces. While the computed velocity of the flow have been predicted relatively well, the computed particle velocity differed strongly from the measured one. Using Monte Carlo simulation of inter-particle collisions, the computation of Yamamoto et al. (2001) was specifically performed to model Eaton's experiment. The results of Yamamoto et al. (2001) improved the particle velocity distribution. At the same time, Vance & Squires (2002) mentioned that the stochastic simualtion of inter-particle collisions is too expensive, requiring significantly more CPU resources than one needs for the gas flow computation. Therefore, the need comes to account for the inter-particle collisions in a simpler and still effective way. To present such a model in the framework of LES/Lagrangian particle approach, and to compare the calculated results with Eaton's measurement and modeling of Yamamoto is the main objective of the present paper.

Gorokhovski, M.↗

NASA NTRS1991

A kill curve for Phanerozoic marine species

A kill curve for Phanerozoic species is developed from an analysis of the stratigraphic ranges of 17,621 genera, as compiled by Sepkoski. The kill curve shows that a typical species' risk of extinction varies greatly, with most time intervals being characterized by very low risk. The mean extinction rate of 0.25/m.y. is thus a mixture of long periods of negligible extinction and occasional pulses of much higher rate. Because the kill curve is merely a description of the fossil record, it does not speak directly to the causes of extinction. The kill curve may be useful, however, to li inverted question markmit choices of extinction mechanisms.

NASA Discipline Number 52-40↗

NASA NTRS2012

Two-Dimensional Distribution of Volatiles in the Lunar Regolith from Space Weathering Simulations

We present simulations of space weathering effects on ice deposits in regions of permanent shadow on the Moon. These Monte Carlo simulations follow the effects of space weathering processes on the distribution of the volatiles over time. The model output constrains the coherence of volatile deposits with depth, lateral separation, and time. The results suggest that ice sheets become broken and buried with time. As impacts begin to puncture an initially coherent surficial ice sheet, small areas with a deficit of ice compared to surrounding areas are formed first. As time progresses, holes become prevalent and the anomalous regions are local enhancements of ice concentration in a volume. The 3-D distribution is also heterogeneous because the ice is buried to varying depths in different locations. Analysis of the coherence of ice on 10 cm scales predicts that putative ice sheets in anomalous radar craters are < 100 Myr old. Surface frost becomes homogenized within 20 Myr. The simulations show the data from the LCROSS impact and surrounding region are consistent with the ice deposit in Cabeus being >1000 Myr old. For future in situ analysis of cold trap volatiles, a horizontal range of 10 m is sufficient to acquire surface-based measurements of heterogeneously distributed ice. These results also support previous analyses that Mercury's cold traps are young.

Hurley, Dana M.↗

NASA NTRS2013

The LF of TP-AGB stars in the LMC/SMC

We show that Monte Carlo simulations of the TP-AGB stellar population in the LMC and SMC galaxies using the CB. models produce LF and color distributions that are in closer agreement with observations than those obtained with the BC03 and CB07 models. This is a progress report of work that will be published elsewhere.

Bruzual, Gustavo↗

NASA NTRS2019

Benchmarking Simulations of the Compton Spectrometer and Imager with Calibrations

The Compton Spectrometer and Imager (COSI) is a balloon-borne 𝛾-ray (0.2-5 MeV) telescope designed to study astrophysical sources. COSI employs a compact Compton telescope design utilizing 12 high-purity germanium double-sided strip detectors and is inherently sensitive to polarization. In 2016, COSI was launched from Wanaka, New Zealand and completed a successful 46-day flight on NASA’s new Super Pressure Balloon. In order to perform imaging, spectral, and polarization analysis of the sources observed during the 2016 flight, we compute the detector response from well-benchmarked simulations. As required for accurate simulations of the instrument, we have built a comprehensive mass model of the instrument and developed a detailed detector effects engine which applies the intrinsic detector performance to Monte Carlo simulations. The simulated detector effects include energy, position, and timing resolution, thresholds, dead strips, charge sharing, charge loss, crosstalk, dead time, and detector trigger conditions. After including these effects, the simulations closely resemble the measurements, the standard analysis pipeline used for measurements can also be applied to the simulations, and the responses computed from the simulations are accurate. We have computed the systematic error that we must apply to measured fluxes at certain energies, which is 6.3% on a rage. Here we describe the detector effects engine and the benchmarking tests performed with calibrations.

Sleator, Clio C.↗

NASA NTRS2020

InSight Entry, Descent and Landing Pre-Flight Performance Predictions

On November 26, 2018, the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) lander successfully touched down on the surface of Mars. Over its seven-plus year development, NASA Langley Research Center’s (LaRC) Program to Optimize Simulated Trajectories II (POST2) was used to assess the mission’s Entry, Descent and Landing (EDL) vehicle system performance against related requirements across the full range of possible environmental and spacecraft conditions. Much of the simulation code was derived from the Phoenix mission, for which this vehicle is very similar. The InSight six degree-of-freedom simulation included models for Mars atmosphere, gravity and digital elevation maps of the landing location. Additionally, vehicle specific aerodynamic, parachute, engine, navigation sensor, flight software and landing radar models were also included. A set of dispersions for each model, as well as for additional simulation input parameters, were also included in order to provide a statistical, Monte Carlo prediction of the EDL system performance. An overview of the pre-flight performance assessments completed, including the various simulation campaigns used, will be provided. Ultimately, this work was critical in the assessment of readiness for InSight launch. A brief description of the use of this simulation in support of flight operations is also discussed.

Robert W Maddock↗

NASA NTRSDate not supplied

Comparison of Dynamic Stability Testing Techniques with Magnetic Suspension Wind Tunnel

Dynamic stability testing techniques currently utilized at NASA Langley Research Center (LaRC) are conducted in multiple facilities and consists of free flight, forced oscillation, and free-to-oscillate tests. The MIT/NASA/ODU Magnetic Suspension and Balance System (MSBS) has been recommissioned to explore its utility as an additional facility to expand the dynamic stability test capabilities currently available at NASA LaRC. Simulations were created to replicate each current test facility and method as closely as possible. Data collected from the simulated environments was corrupted with replicated noise sources of the different testing environments and then compared to real data collected during tests when such data was available. The corrupted data was then passed through data reduction and System Identification (SID) to estimate the accuracy of the results with the known aerodynamic model that was utilized within the simulation to generate the original data. Magnitudes of noise were varied utilizing Monte Carlo analysis to perform sensitivity analysis of each noise source on the extracted dynamic stability coefficients. Some preliminary results will be presented.

Otoniel Ramirez↗

NASA NTRSDate not supplied

Comparison of Dynamic Stability Testing Techniques with Magnetic Suspension Wind Tunnel

Otoniel Ramirez↗

NASA NTRS2003

Mesoscale Computational Investigation of Shocked Heterogeneous Materials with Application to Large Impact Craters

The propagation of shock waves through target materials is strongly influenced by the presence of small-scale structure, fractures, physical and chemical heterogeneities. Pre-existing fractures often create craters that appear square in outline (e.g. Meteor Crater). Reverberations behind the shock from the presence of physical heterogeneity have been proposed as a mechanism for transient weakening of target materials. Pre-existing fractures can also affect melt generation. In this study, we are attempting to bridge the gap in numerical modeling between the micro-scale and the continuum, the so-called meso-scale. To accomplish this, we are developing a methodology to be used in the shock physics hydrocode (CTH) using Monte-Carlo-type methods to investigate the shock properties of heterogeneous materials. By comparing the results of numerical experiments at the micro-scale with experimental results and by using statistical techniques to evaluate the performance of simple constitutive models, we hope to embed the effect of physical heterogeneity into the field variables (pressure, stress, density, velocity) allowing us to directly imprint the effects of micro-scale heterogeneity at the continuum level without incurring high computational cost.

Crawford, D. A.↗

NASA NTRS2014

Optimizing Medical Kits for Spaceflight

The Integrated Medical Model (IMM) is a probabilistic model that estimates medical event occurrences and mission outcomes for different mission profiles. IMM simulation outcomes describing the impact of medical events on the mission may be used to optimize the allocation of resources in medical kits. Efficient allocation of medical resources, subject to certain mass and volume constraints, is crucial to ensuring the best outcomes of in-flight medical events. We implement a new approach to this medical kit optimization problem. METHODS We frame medical kit optimization as a modified knapsack problem and implement an algorithm utilizing a dynamic programming technique. Using this algorithm, optimized medical kits were generated for 3 different mission scenarios with the goal of minimizing the probability of evacuation and maximizing the Crew Health Index (CHI) for each mission subject to mass and volume constraints. Simulation outcomes using these kits were also compared to outcomes using kits optimized..RESULTS The optimized medical kits generated by the algorithm described here resulted in predicted mission outcomes more closely approached the unlimited-resource scenario for Crew Health Index (CHI) than the implementation in under all optimization priorities. Furthermore, the approach described here improves upon in reducing evacuation when the optimization priority is minimizing the probability of evacuation. CONCLUSIONS This algorithm provides an efficient, effective means to objectively allocate medical resources for spaceflight missions using the Integrated Medical Model.

Optimization↗

NASA NTRS2015

Radiation Environment Model of Protons and Heavier Ions at Jupiter

We performed an in depth study of the methods used to review the geometric factors (GF) and sensitivity to charge particles of the Energetic Particle Detector instrument on board the Galileo Spacecraft. Monte Carlo simulations were performed to understand the interactions of electrons and ions (i. e., protons and alphas) with the sensitive regions of the instrument. The DC0 and B0 channels were studied with the intention of using them to update the jovian proton radiation model. The results proved that the B0 is a clean proton chanel without any concerns for contamination by heavier ions and electrons. In contrast, DC0 was found to be contaminated by electrons. Furthermore, we also found out that the B2 channel is a clean alpha particle channel (in other words, no contamination by electrons and/or protons).

Jupiter Radiation Environment↗

NASA NTRS2013

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species With Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for nearequilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion of the heating and is then compared to the total heating measured in flight.

Liechty, Derek S.↗

NASA NTRS2014

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure Induced by HZE Particles

The formation of double‐strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte‐Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high‐energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation‐induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex‐type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05‐0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta‐rays becomes more important.

Ponomarev, Artem↗