Monte Carlo simulation of an adaptive policy for multiple impulse correction of the trajectory of a spacecraft.
Monte Carlo simulation for adaptive policy for multiple impulse correction of trajectory of spacecraft
MIDCOURSE GUIDANCE↗
SEARCH · Search NASA
Results for “Monte Carlo Simulation”
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.
Applying Monte Carlo Simulation to Launch Vehicle Design and Requirements Verification
This paper is focused on applying Monte Carlo simulation to probabilistic launch vehicle design and requirements verification. The approaches developed in this paper can be applied to other complex design efforts as well. Typically the verification must show that requirement "x" is met for at least "y" % of cases, with, say, 10% consumer risk or 90% confidence. Two particular aspects of making these runs for requirements verification will be explored in this paper. First, there are several types of uncertainties that should be handled in different ways, depending on when they become known (or not). The paper describes how to handle different types of uncertainties and how to develop vehicle models that can be used to examine their characteristics. This includes items that are not known exactly during the design phase but that will be known for each assembled vehicle (can be used to determine the payload capability and overall behavior of that vehicle), other items that become known before or on flight day (can be used for flight day trajectory design and go/no go decision), and items that remain unknown on flight day. Second, this paper explains a method (order statistics) for determining whether certain probabilistic requirements are met or not and enables the user to determine how many Monte Carlo samples are required. Order statistics is not new, but may not be known in general to the GN&C community. The methods also apply to determining the design values of parameters of interest in driving the vehicle design. The paper briefly discusses when it is desirable to fit a distribution to the experimental Monte Carlo results rather than using order statistics.
Hanson, John M.↗
Radiation modeling with direct simulation Monte Carlo
Improvements in the modeling of radiation in low density shock waves with direct simulation Monte Carlo (DSMC) are the subject of this study. A new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. This scheme attempts to move the DSMC programs toward a more detailed modeling of the physics and more reliance on available experimental data. The new method is compared with the current modeling technique and both techniques are compared with available data. The differences in the results are evaluated. The test case is based on an AVCO-Everett shock tube experiment, a 10-km/s standing shock wave in air at 0.1 Torr. The new method agrees with the available data as well as the results from the earlier scheme and is more easily extrapolated to different flow conditions.
Carlson, Ann B.↗
Monte Carlo simulation of electron acceleration in modified relativistic shocks
We give a brief review of Monte Carlo simulations of nonlinear Fermi shock acceleration and then give new results on electron acceleration in SNRs and in relativistic parallel shocks. The acceleration of low energy electrons in shocks is poorly understood, but even when energetic electrons are considered, where electron and proton scattering should be qualitatively similar, dramatic differences result between electron and proton acceleration in relativistic shocks. If the shocked plasma is a mixture of electrons and protons, the electrons are accelerated much less efficiently than protons. We predict that only e(-)-e(+) pair dominated plasmas can produce significant radio emission in relativistic flows if the standard Fermi mechanism operates in parallel shocks.
Ellison, Donald C.↗
Radiation Modeling with Direct Simulation Monte Carlo
Improvements in the modeling of radiation in low density shock waves with direct simulation Monte Carlo (DSMC) are the subject of this study. A new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. This scheme attempts to move the DSMC programs toward a more detailed modeling of the physics and more reliance on available rate data. The new method is compared with the current modeling technique and both techniques are compared with available experimental data. The differences in the results are evaluated. The test case is based on experimental measurements from the AVCO-Everett Research Laboratory electric arc-driven shock tube of a normal shock wave in air at 10 km/s and .1 Torr. The new method agrees with the available data as well as the results from the earlier scheme and is more easily extrapolated to di erent ow conditions.
Carlson, Ann B.↗
Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys
In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.
Bentz, Daniel N.↗
A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data
The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.
DAC↗
Monte Carlo simulation of lunar megaregolith and implications
A realistic Monte Carlo model closely simulating the evolution of the lunar megaregolith over a large area of 67 million sq. km of the front surface of the moon is presented. Craters larger than 100 km in diameter observed over the entire surface of the moon and those less than 100 km lying in the referenced area are included in the simulation. A total of 21,664 craters are processed. The model predicts the average thickness of the megaregolith to be about 1.9-2.0 km. Curves for the variation of the regolith thickness across the simulated area are given and show that about 50% of the area is covered with regolith less than 1 km thick. The model produces crater structures similar to the ones observed in the lunar highlands, it partially supports the layering theory for crater structures that the variations in strength and density of target materials may be responsible for the observed differences in the morphologies of lunar craters, and rules out the possibility that all craters when formed are bowl-shaped with a fixed depth/diameter ratio characteristic of small craters.
Aggarwal, H. R.↗
Monte Carlo Simulation of a Knudsen Effusion Mass Spectrometer Sampling System
Knudsen flow is easily simulated with a Monte Carlo method. In this study we develop a Visual Basic for Excel (VBA) code to simulate the molecular beam from a vaporizing solid. The system at NASA Glenn uses the restricted collimation method of Chatillon and colleagues, which consists of two apertures between the effusion cell and the ionizer. The diameter of the first aperture is smaller than the diameter of the effusion cell orifice, so the ionizer effectively only sees inside the effusion cell. The code is able to calculate the transmission coefficient through the cell orifice, through the cell orifice and the first aperture, and through the cell orifice and first and second aperatures. Calculated transmission coefficients through the cell orifice are compared to tabulated values to validate the code. Then transmission coefficients are calculated through the cell orifice and both apertures to the ionizer. This allows the geometry (aperture spacing and diameters) of the sampling system the sampling system to be optimized. Calculated transmission factors are also compared to literature values calculated via an analytic method.
Radke, Michael J.↗
A Time-Slicing Method for Capturing Neutron Source Locations for Initializing Dynamic Monte Carlo Simulations
Explore the source record for details and available documents.
Monte Carlo transport, neutronics, numerical metho↗
Monte Carlo simulations of particle acceleration at oblique shocks
The Fermi shock acceleration mechanism may be responsible for the production of high-energy cosmic rays in a wide variety of environments. Modeling of this phenomenon has largely focused on plane-parallel shocks, and one of the most promising techniques for its study is the Monte Carlo simulation of particle transport in shocked fluid flows. One of the principal problems in shock acceleration theory is the mechanism and efficiency of injection of particles from the thermal gas into the accelerated population. The Monte Carlo technique is ideally suited to addressing the injection problem directly, and previous applications of it to the quasi-parallel Earth bow shock led to very successful modeling of proton and heavy ion spectra, as well as other observed quantities. Recently this technique has been extended to oblique shock geometries, in which the upstream magnetic field makes a significant angle Theta(sub B1) to the shock normal. Spectral resutls from test particle Monte Carlo simulations of cosmic-ray acceleration at oblique, nonrelativistic shocks are presented. The results show that low Mach number shocks have injection efficiencies that are relatively insensitive to (though not independent of) the shock obliquity, but that there is a dramatic drop in efficiency for shocks of Mach number 30 or more as the obliquity increases above 15 deg. Cosmic-ray distributions just upstream of the shock reveal prominent bumps at energies below the thermal peak; these disappear far upstream but might be observable features close to astrophysical shocks.
Baring, Matthew G.↗
Monte Carlo Simulations of 347H Stainless Steel Aging for the Synthetic Generation of Microstructures Under Creep Conditions
Here, a Monte Carlo simulation method capable of replicating the kinetics of M 23 C 6 precipitation in 347H stainless steels was developed for the purpose of producing synthetic microstructures that approximate its microstructural evolution under aging periods of up to 10,000 hours at temperatures between 600 °C and 750 °C. To accomplish this, experimental data from the literature was used to parameterize simulations and replicate the nucleation and growth kinetics of M 23 C 6 particles within 347H and similar austenitic stainless steel alloys. These simulations were found to have considerable fidelity to previous efforts to study the precipitation of M 23 C 6 in other 300 series stainless steel alloys. Synthetic 347H microstructures were then generated that accounted the effects of aging temperature, duration, dislocation density, and the presence of boron within the microstructure. These simulations predict several key trends, those being that (1) the size of M 23 C 6 precipitates decreased with aging temperature and (2) the growth rate of M 23 C 6 particles decreased with aging temperature. Further, while (3) the addition of dislocation density due to creep conditions resulted in increasing intragranular nucleation of M 23 C 6 precipitates with increasing dislocation density and (4) B additions within the microstructure led to modest increases in precipitate size above 700 °C, which indicates that more complex physics are necessary to account for the presence of B.
36 MATERIALS SCIENCE↗
Result of Monte-Carlo simulation of electron-photon cascades in lead and layers of lead-scintillator
Results of Monte-Carlo simulation of electromagnetic cascade development in lead and lead-scintillator sandwiches are analyzed. It is demonstrated that the structure function for core approximation is not applicable in the case in which the primary energy is higher than 100 GeV. The simulation data has shown that introducing an inhomogeneous chamber structure results in subsequent reduction of secondary particles.
Wasilewski, A.↗
Object-Oriented/Data-Oriented Design of a Direct Simulation Monte Carlo Algorithm
Over the past decade, there has been much progress towards improved phenomenological modeling and algorithmic updates for the direct simulation Monte Carlo (DSMC) method, which provides a probabilistic physical simulation of gas Rows. These improvements have largely been based on the work of the originator of the DSMC method, Graeme Bird. Of primary importance are improved chemistry, internal energy, and physics modeling and a reduction in time to solution. These allow for an expanded range of possible solutions In altitude and velocity space. NASA's current production code, the DSMC Analysis Code (DAC), is well-established and based on Bird's 1994 algorithms written in Fortran 77 and has proven difficult to upgrade. A new DSMC code is being developed in the C++ programming language using object-oriented and data-oriented design paradigms to facilitate the inclusion of the recent improvements and future development activities. The development efforts on the new code, the Multiphysics Algorithm with Particles (MAP), are described, and performance comparisons are made with DAC.
Liechty, Derek S.↗
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.↗
Direct simulation Monte Carlo with ionization and radiation
Improvements in the modeling of radiation in low density shock waves with Direct Simulation Monte Carlo (DSMC) are the subject of this study. The physical processes which determine the amount of radiation in a shock wave were investigated and the way in which they were modeled with DSMC was evaluated. Three physical processes were identified for which an improvement in the modeling technique could result in improved radiation predictions. New physical modeling schemes are introduced in this report for the three processes. First a method for determining the electric field and its effect on the flow is introduced. Second, a two step reaction process for electron impact ionization reactions is evaluated. Finally, a new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. Each new scheme attempts to move the DSMC method toward more accurate physics or more reliance on experimental data. The new schemes are all compared to the current modeling techniques and the differences in the results are evaluated. In all cases the results agree with the available data as well as, or better than the results from earlier schemes.
Carlson, Ann B.↗
Procedure for Adapting Direct Simulation Monte Carlo Meshes
A technique is presented for adapting computational meshes used in the G2 version of the direct simulation Monte Carlo method. The physical ideas underlying the technique are discussed, and adaptation formulas are developed for use on solutions generated from an initial mesh. The effect of statistical scatter on adaptation is addressed, and results demonstrate the ability of this technique to achieve more accurate results without increasing necessary computational resources.
Woronowicz, Michael S.↗
An Investigation of Parallel Programming Techniques Applied to Monte Carlo Simulations for Post-Flight Reconstruction of Spacecraft Trajectory
Parallelizing software to execute on multi-core central processing units (CPUs) and graphics processing units (GPUs) can be challenging. For some fields outside of Computer Science, this transition comes with new issues. For example, memory limitations can require modifications to code not initially developed to run on GPUs. This work applies the Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) directive-based parallelization strategies on a Monte Carlo simulation approach for trajectory reconstruction enabling it to run on multi-core CPUs and GPUs. Large matrix operations are the most common use of GPUs, which are not present in this algorithm; however, the natural parallelism of independent trajectories in Monte Carlo simulations is exploited. Benchmarking data are presented comparing execution times of the software for single-thread CPUs, multi-thread CPUs with OpenMP, and multi-thread GPUs using OpenACC. These data were collected using nodes with Intel® Xeon® E5-2670 (Sandy Bridge) CPUs enhanced with NVIDIA® Tesla® K40 GPUs on the Pleiades Supercomputer cluster at the National Aeronautics and Space Administration (NASA) Ames Research Center (ARC) and a local Intel® Xeon Phi™ node at NASA Langley Research Center (LaRC).
Williams, R. Anthony↗