SEARCH · Search NASA

Results for “modified Monte Carlo”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Carbon Back Sputter Modeling for Hall Thruster Testing

Lifetime requirements for electric propulsion devices, including Hall Effect thrusters, are continually increasing, driven in part by NASA's inclusion of this technology in it's exploration architecture. NASA will demonstrate high-power electric propulsion system on the Solar Electric Propulsion Technology Demonstration Mission (SEP TDM). The Asteroid Redirect Robotic mission is one candidate SEP TDM, which is projected to require tens of thousands of thruster life. As thruster life is increased, for example through the use of improved magnetic field designs, the relative influence of facility effects increases. One such effect is the sputtering and redeposition, or back sputter, of facility materials by the high energy thruster plumes. In support of wear testing for the Hall Effect Rocket with Magnetic Shielding (HERMeS) project, the back sputter from a Hall effect thruster plume has been modeled for the NASA Glenn Research Center's Vacuum Facility 5. The predicted wear at a near-worst case condition of 600 V, 12.5 kW was found to be on the order of 1 micron/kh in a fully carbon-lined chamber. A more detailed numerical Monte Carlo code was also modified to estimate back sputter for a detailed facility and pumping configuration. This code demonstrated similar back sputter rate distributions, but is not yet accurately modeling the magnitudes. The modeling has been benchmarked to recent HERMeS wear testing, using multiple microbalance measurements. These recent measurements have yielded values on the order of 1.5 - 2 micron/kh at 600 V and 12.5 kW.

Propulsion↗

NASA NTRS2016

Carbon Back Sputter Modeling for Hall Thruster Testing

In support of wear testing for the Hall Effect Rocket with Magnetic Shielding (HERMeS) program, the back sputter from a Hall effect thruster plume has been modeled for the NASA Glenn Research Centers Vacuum Facility 5. The predicted wear at a near-worst case condition of 600 V, 12.5 kW was found to be on the order of 3 4 mkhour in a fully carbon-lined chamber. A more detailed numerical monte carlo code was also modified to estimate back sputter for a detailed facility and pumping configuration. This code demonstrated similar back sputter rate distributions, but is not yet accurately modeling the magnitudes. The modeling has been benchmarked to recent HERMeS wear testing, using multiple microbalance measurements. These recent measurements have yielded values, on the order of 1.5- 2 microns/khour.

Facility↗

NASA NTRS2022

Coupling Carbon Oxidation and Surface Recession in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that are often made out of carbon-based reinforcement and a polymeric matrix. They endure high-temperature oxidation and surface recession when re-entering Earth’s atmosphere. Ablation is the result of many coupled and competing thermal, mechanical, and chemical phenomena, and it is difficult to isolate the role of each on the overall degradation of the TPS. Here we develop an ablation model for material recession coupled explicitly to finite rate carbon oxidation in complex microstructures. In this work, Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al, was adapted to tally surface reactions and other surface data on a per-grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also develop robust algorithms which handle the evolution of the flow cells and solid material regions, including split cells (flow cell divided in two by a solid surface). Finally, we demonstrate our implicit chemistry model for 2D and 3D geometries by producing reaction statistics and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗

NASA NTRSDate not supplied

A Parallelized Oxidation-Driven Surface Recession Framework in DSMC Code, SPARTA

Spacecrafts rely on ablative thermal protection systems (TPS) made of composites consisting of a carbon-based reinforcement and a polymeric matrix. These materials are designed to withstand high-temperature oxidation and surface recession during re-entry into the Earth's atmosphere. However, ablation occurs due to a complex interplay of thermal, mechanical, and chemical factors, making it challenging to determine the individual impact of each on the TPS's overall degradation. In this study, we have developed an ablation model that can leverage a finite rate carbon oxidation model to predict material recession and surface states more accurately. Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al can perform both gas-surface and pure-surface reactions and is now adapted to tally surface data on a per grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also briefly discuss some unique challenges associated with parallelizing this dynamic surface state and geometry. Finally, we analyze the performance of this parallelized implicit chemistry model with simple 2D and 3D benchmark cases by producing surface state statistics, area changes over time, and visualization across a range of surface temperatures and processors with and without load-balancing.

DSMC↗

NASA NTRS2019

Characterizing the System Impulse Response Function from Photon-Counting LiDAR Data

NASA's Multiple Altimeter Beam Experimental LiDAR (MABEL) is an aircraft-based photon-counting laser altimeter designed as a simulator to test measurement techniques and algorithms for Advanced Topographic Laser Altimeter System (ATLAS), the sole instrument on NASA's Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) mission. By measuring the time of flight, pointing angle, and absolute position for individual photons, ICESat-2 provides detailed elevation measurements of earth's surface. Calculating accurate and precise elevations requires an understanding of how photons interact with surfaces, and characterization of the photon distribution after returning from surfaces. Neither MABEL nor ATLAS records the transmitted laser pulse shape, relying instead on aggregating several pulses worth of photons, often using histograms, to characterize the pulse shape. In this paper, we assess the limitations of using histograms and propose a more robust method to describe MABEL's system impulse-response function using an exponentially modified Gaussian distribution. We also provide standard error estimates for the arithmetic mean and standard deviation calculations, and for exponentially modified Gaussian parameters using a Monte Carlo sensitivity analysis. We apply this method to photon returns from a sea ice lead and from a dry salt lake bed as case studies for estimating the standard error associated with sample size for the arithmetic mean and standard deviation, and for the exponentially modified Gaussian parameters. We use these standard errors to calculate the minimum number of photons required to find both Gaussian and exponentially modified Gaussian distribution parameters within 3 cm of their parent population values.

photoncounting↗

NASA NTRS1983

Computer simulation study of multiple germanium gamma-ray sensor arrays

Design studies of large high resolution gamma-ray spectrometer systems employing multiple sensors are greatly facilitated by the use of computer simulation techniques. A study has been made by interfacing a detector geometry code with the EGS-3 version of the Stanford Linear Accelerator Center Monte-Carlo Code which has been modified to run on our VAX computer. The detector responses have been simulated for input photons with energies between 0.1 and 10 MeV. Examples are given to illustrate the performance of complex sensor arrays with different configurations of the anticoincidence system. Of particular interest is a design composed of an array of 7 detector elements each consisting of a thin planar Ge (1 cm thick) detector mounted in front of and in tandem with a large coaxial Ge detector.

Nakano, G. H.↗

NASA NTRS1989

Computation of electron density in high altitude re-entry flows

The direct simulation Monte Carlo (DSMC) method has been modified to deal with the very small degree of ionization that occurs in the early stages of the 'communications blackout' phase of re-entry. The simplified model of the process assumes that the ionization is entirely a consequence of the reaction N + O yields NO(+) + e-. The measurements from the RAM C-II test vehicle at an altitude of 81 km have been chosen as the test case. It was found that the form of the electron distribution along the shock layer was in agreement with the measurements. That normal to the surface was also in agreement with the measurements if the surface is regarded as fully catalytic. However, quantitative agreement for the fully catalytic case requires a reaction rate about an order of magnitude greater than those generally associated with this reaction.

Bird, G. A.↗

NASA NTRS1991

Analysis of vibrational-translational energy transfer using the direct simulation Monte Carlo method

A new model is proposed for energy transfer between the vibrational and translational modes for use in the direct simulation Monte Carlo method (DSMC). The model modifies the Landau-Teller theory for a harmonic oscillator and the rate transition is related to an experimental correlation for the vibrational relaxation time. Assessment of the model is made with respect to three different computations: relaxation in a heat bath, a one-dimensional shock wave, and hypersonic flow over a two-dimensional wedge. These studies verify that the model achieves detailed balance, and excellent agreement with experimental data is obtained in the shock wave calculation. The wedge flow computation reveals that the usual phenomenological method for simulating vibrational nonequilibrium in the DSMC technique predicts much higher vibrational temperatures in the wake region.

Boyd, Iain D.↗

NASA NTRS1997

Entry Dispersion Analysis for the Stardust Comet Sample Return Capsule

Stardust will be the first mission to return samples from beyond the Earth-Moon system. The sample return capsule, which is passively controlled during the fastest Earth entry ever, will land by parachute in Utah. The present study analyzes the entry, descent, and landing of the returning sample capsule. The effects of two aerodynamic instabilities are revealed (one in the high altitude free molecular regime and the other in the transonic/subsonic flow regime). These instabilities could lead to unacceptably large excursions in the angle-of-attack near peak heating and main parachute deployment, respectively. To reduce the excursions resulting from the high altitude instability, the entry spin rate of the capsule is increased. To stabilize the excursions from the transonic/subsonic instability, a drogue chute with deployment triggered by an accelerometer and timer is added prior to main parachute deployment. A Monte Carlo dispersion analysis of the modified entry (from which the impact of off-nominal conditions during the entry is ascertained) shows that the capsule attitude excursions near peak heating and drogue chute deployment are within Stardust program limits. Additionally, the size of the resulting 3-sigma landing ellipse is 83.5 km in downrange by 29.2 km in crossrange, which is within the Utah Test and Training Range boundaries.

Desai, Prasun N.↗

NASA NTRS1999

Entry Trajectory Issues for the Stardust Sample Return Capsule

The Stardust mission was successfully launched on February 7, 1999. It will be the first mission to return samples from a comet. The sample return capsule, which is passively controlled during the fastest Earth entry ever, will land by parachute in Utah. The present study describes the analysis of the entry, descent, and landing of the returning sample capsule utilizing the final, launch configuration capsule mass properties. The effects of two aerodynamic instabilities are revealed (one in the high altitude free molecular regime and the other in the transonic/subsonic flow regime). These instabilities could lead to unacceptably large excursions in the angle-of-attack near peak heating and main parachute deployment, respectively. To reduce the excursions resulting from the high altitude instability, the entry spin rate of the capsule is increased. To stabilize the excursions from the transonic/subsonic instability, a drogue chute with deployment triggered by a gravity-switch and timer is added prior to main parachute deployment. A Monte Carlo dispersion analysis of the modified entry (from which the impact of off-nominal conditions during the entry is ascertained) predicts that the capsule attitude excursions near peak heating and drogue chute deployment are within Stardust mission limits. Additionally, the size of the resulting 3-sigma landing ellipse is 60.8 km in downrange by 19.9 km in crossrange, which is within the Utah Test and Training Range boundaries.

Desai, Prasun N.↗

NASA NTRS1999

Entry Trajectory Issues for the Stardust Sample Return Capsule

Desai, Prasun N.↗

NASA NTRS1992

Particle acceleration in modified oblique non-relativistic shocks

Previous applications of the Monte Carlo technique at the quasiparallel earth bow shock has motivated the extension of this technique to oblique shock geometries typical of those found in most astrophysical shock environments. In addition, such a generalization will permit the thorough examination of theoretical predictions of rapid acceleration times at quasi-perpendicular shocks. Therefore, we have embarked on the modification of our existing Monte Carlo code and in this paper outline the major technical aspects involved in developing a simulation of cosmic-ray acceleration at modified oblique nonrelativistic shocks.

Baring, Matthew G.↗

NASA NTRS1975

The structure and evolution of Jupiter - The fluid contraction stage

The complete evolution of a contracting star of Jovian mass consisting of a convective adiabatic homogeneous fluid is determined using stellar structure methods, improved model atmosphere calculations, and substantially improved thermodynamic properties for hydrogen and hydrogen-helium fluids. The model atmospheres are calculated in the form of time-averaged vertical temperature structures, including all relevant sources of opacity and a solar energy deposition component, and the thermodynamic properties are modified to obtain better agreement with Monte Carlo results for metallic fluids. The resultant gravitationally contracting evolutionary models are found to have two phases: an early stellar phase similar to a typical low-mass pre-main-sequence body and a later phase constituting an approach to a degenerate-dwarf cooling curve. The first phase is shown to have high luminosities and internal temperatures, while the second gives excellent agreement with the observed radius and luminosity of Jupiter. Analysis indicates that the equation of state and superadiabaticity have the strongest influence on evolution over planetary time scales.

Graboske, H. C., Jr.↗

NASA NTRS1993

Monte Carlo modeling of atomic oxygen attack of polymers with protective coatings on LDEF

Characterization of the behavior of atomic oxygen interaction with materials on the Long Duration Exposure Facility (LDEF) assists in understanding of the mechanisms involved. Thus the reliability of predicting in-space durability of materials based on ground laboratory testing should be improved. A computational model which simulates atomic oxygen interaction with protected polymers was developed using Monte Carlo techniques. Through the use of an assumed mechanistic behavior of atomic oxygen interaction based on in-space atomic oxygen erosion of unprotected polymers and ground laboratory atomic oxygen interaction with protected polymers, prediction of atomic oxygen interaction with protected polymers on LDEF was accomplished. However, the results of these predictions are not consistent with the observed LDEF results at defect sites in protected polymers. Improved agreement between observed LDEF results and predicted Monte Carlo modeling can be achieved by modifying of the atomic oxygen interactive assumptions used in the model. LDEF atomic oxygen undercutting results, modeling assumptions, and implications are presented.

Banks, Bruce A.↗

NASA NTRS2016

MAVEN Transition-To-Science Maneuvers and Comet Siding Spring

After the September 2014 arrival of the Mars Atmosphere Volatile Evolution mission (MAVEN) at Mars, the spacecraft was to perform five maneuvers to reduce the orbit period from 35 to 4.5 hours, and lower the periapsis altitude into the target in-situ science density corridor. While in route, however, the potential threat — and science opportunity — from Comet C/2013 Siding Spring crystallized, and the plan had to be modified to shield the spacecraft behind Mars and support new observations. This paper describes the plan modifications, Monte Carlo analyses necessary to ensure success, and results of executing the modified plan.

Young, Brian↗

NASA NTRS2004

Aerocapture Performance Analysis for a Neptune-Triton Exploration Mission

A systems analysis has been conducted for a Neptune-Triton Exploration Mission in which aerocapture is used to capture a spacecraft at Neptune. Aerocapture uses aerodynamic drag instead of propulsion to decelerate from the interplanetary approach trajectory to a captured orbit during a single pass through the atmosphere. After capture, propulsion is used to move the spacecraft from the initial captured orbit to the desired science orbit. A preliminary assessment identified that a spacecraft with a lift to drag ratio of 0.8 was required for aerocapture. Performance analyses of the 0.8 L/D vehicle were performed using a high fidelity flight simulation within a Monte Carlo executive to determine mission success statistics. The simulation was the Program to Optimize Simulated Trajectories (POST) modified to include Neptune specific atmospheric and planet models, spacecraft aerodynamic characteristics, and interplanetary trajectory models. To these were added autonomous guidance and pseudo flight controller models. The Monte Carlo analyses incorporated approach trajectory delivery errors, aerodynamic characteristics uncertainties, and atmospheric density variations. Monte Carlo analyses were performed for a reference set of uncertainties and sets of uncertainties modified to produce increased and reduced atmospheric variability. For the reference uncertainties, the 0.8 L/D flatbottom ellipsled vehicle achieves 100% successful capture and has a 99.87 probability of attaining the science orbit with a 360 m/s V budget for apoapsis and periapsis adjustment. Monte Carlo analyses were also performed for a guidance system that modulates both bank angle and angle of attack with the reference set of uncertainties. An alpha and bank modulation guidance system reduces the 99.87 percentile DELTA V 173 m/s (48%) to 187 m/s for the reference set of uncertainties.

Starr, Brett R.↗

NASA NTRS1994

A probabilistic fatigue analysis of multiple site damage

The variability in initial crack size and fatigue crack growth is incorporated in a probabilistic model that is used to predict the fatigue lives for unstiffened aluminum alloy panels containing multiple site damage (MSD). The uncertainty of the damage in the MSD panel is represented by a distribution of fatigue crack lengths that are analytically derived from equivalent initial flaw sizes. The variability in fatigue crack growth rate is characterized by stochastic descriptions of crack growth parameters for a modified Paris crack growth law. A Monte-Carlo simulation explicitly describes the MSD panel by randomly selecting values from the stochastic variables and then grows the MSD cracks with a deterministic fatigue model until the panel fails. Different simulations investigate the influences of the fatigue variability on the distributions of remaining fatigue lives. Six cases that consider fixed and variable conditions of initial crack size and fatigue crack growth rate are examined. The crack size distribution exhibited a dominant effect on the remaining fatigue life distribution, and the variable crack growth rate exhibited a lesser effect on the distribution. In addition, the probabilistic model predicted that only a small percentage of the life remains after a lead crack develops in the MSD panel.

Rohrbaugh, S. M.↗

NASA NTRS2002

Evaluation of the Effect of Water Vapor On the Performance of NASA's NMRO Catalysts for Carbon Monoxide Oxidation

The Noble Metal Reducible Oxide (NMRO) catalysts for the low temperature oxidation of carbon monoxide were developed by NASA for the reoxidation of carbon monoxide which forms by the dissociation of carbon dioxide during the operation of sealed carbon dioxide lasers. The NMRO catalyst, which consists of a noble metal in conjunction with a reducible metal oxide, was evaluated under conditions that will be encountered in a CO2 laser operation, namely temperatures in the range 298 to 373 K and no significant reaction gas components other than CO, CO2 and O2. The NMRO catalysts may have significant potential for spin-off applications such as the prevention of carbon monoxide build-up in closed spaces as in space vehicle cabins or submarines, and the elimination of the cold start-up problem of automobile exhaust catalysts. The most significant difference in the conditions of these possible future applications is the high moisture content of the gases to be processed. Lack of understanding of the effects of water vapor and high temperature on catalyst activity and operation for extended periods are currently the main stumbling blocks for the transfer of this NASA technology to be used for commercial purposes. In the original proposal the following objectives were stated: To obtain experimental data on the adsorption, desorption and reaction characteristics of CO and CO2 the catalysts under high moisture conditions; to collect evidence on the presence of carbonate and hydroxyl surface species and their involvement in the CO oxidation mechanism; and to model the reaction system using a Monte-Carlo simulation to gain insight on the various steps involved. After the work has commenced the NASA technical monitor Mr. David Scheyer informed us that there was increased interest in the possible use of the NMRO catalysts as automobile exhaust catalysts and therefore NASA wanted to know whether the catalysts can operate at high temperatures as well as with high moisture gases. At that meeting it was decided that investigation of the high temperature performance of the NMRO catalysts should be given priority and replace the Monte-Carlo simulation objective. As a result, the modified objectives of the investigation were taken as the investigation of the high-temperature activity of the NMRO catalysts, and the effect of water vapor on the performance of these catalysts.

Akyurtlu, Ates↗