Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct simulation 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.

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

At least 379 records · Page 21

NTR plume modeling

Viewgraphs on nuclear thermal propulsion are presented. Topics covered include computational fluid dynamics (CFD) for plume analysis; molecular fluid dynamics; molecular CFD characteristics; direct-simulation Monte-Carlo (DSMC) method; integration of DSMC and Navier-Stokes computations; and density profiles.

Byers, D.↗

Reliability assessment of thrust chamber cooling concepts using probabilistic analysis techniques

The reliability of OFHC (Oxygen Free High Conductivity) copper and NARloy-Z thrust chambers is assessed by applying probabilistic structural analysis techniques to incorporate design parameter variability and uncertainty. Thrust chambers specifically evaluated are the cylindrical test fixtures employed in a plug-nozzle configuration at the NASA Lewis Research Center. Direct sampling Monte Carlo simulations based on a simplified life prediction methodology established probability densities of firing cycles to structural failure. Simulated cyclic lives demonstrated modest agreement to experiment. Similarly, regions of high structural failure probability were determined using a limit state approach employing calculated cumulative distribution functions for effective stress response and an assumed material strength distribution. A probability of failure of 0.012 was calculated at the center of the coolant channel hot-gas-side wall for an OFHC milled channel. Structural response was found to be sensitive to the uncertainties in the thrust chamber thermal environment and the material's thermal expansion coefficient.

Rapp, Douglas C.↗

High-performance computing for flight vehicles; Proceedings of the Symposium, Washington, Dec. 7-9, 1992

The present conference discusses high-performance computing systems for flight vehicles, large-scale simulations on high-performance flight computers and software, multidisciplinary and design/optimization applications of computers, computational electromagnetics and acoustics, the simulation of aircraft powerplant turbomachinery and reacting flows, and flow calculations on parallel machines. Also discussed are direct flow simulation Monte Carlo methods, structural mechanics sensitivity and fracture calculations on parallel machines, grid-generation and advanced algorithms for CFD, advanced solid-mechanics and structures applications, and advancements in flow visualization technology and neural networks.

Noor, Ahmed K.↗

Reliability assessment of thrust chamber cooling concepts using probabilistic analysis techniques

The reliability of OFHC (Oxygen Free High Conductivity) copper and NARloy-Z thrust chambers is assessed by applying probabilistic structural analysis techniques to incorporate design parameter variability and uncertainty. Thrust chambers specifically evaluated are the cylindrical test fixtures employed in a plug-nozzle configuration at the NASA Lewis Research Center. Direct sampling Monte Carlo simulations based on a simplified life prediction methodology established probability densities of firing cycles to structural failure. Simulated cyclic lives demonstrated modest agreement to experiment. Similarly, regions of high structural failure probability were determined using a limit state approach employing calculated cumulative distribution functions for effective stress response and an assumed material strength distribution. A probability of failure of 0.012 was calculated at the center of the coolant channel hot-gas-side wall for an OFHC milled channel. Structural response was found to be sensitive to the uncertainties in the thrust chamber thermal environment and the material's thermal expansion coefficient.

Rapp, Douglas C.↗

Simulation of Low-density Nozzle Plumes in Non-zero Ambient Pressures

The direct simulation Monte-Carlo (DSMC) method was applied to the analysis of low-density nitrogen plumes exhausting from a small converging-diverging nozzle into finite ambient pressures. Two cases were considered that simulated actual test conditions in a vacuum facility. The numerical simulations readily captured the complicated flow structure of the overexpanded plumes adjusting to the finite ambient pressures, including Mach disks and barrel shaped shocks. The numerical simulations compared well to experimental data of Rothe.

Chung, Chan-Hong↗

Meteors: A Delivery Mechanism of Organic Matter to The Early Earth

All potential exogenous pre-biotic matter arrived to Earth by ways of our atmosphere, where much material was ablated during a luminous phase called 1. meteors" in rarefied flows of high (up to 270) Mach number. The recent Leonid showers offered a first glimpse into the elusive physical conditions of the ablation process and atmospheric chemistry associated with high-speed meteors. Molecular emissions were detected that trace a meteor's brilliant light to a 4,300 K warm wake rather than to the meteor's head. A new theoretical approach using the direct simulation by Monte Carlo technique identified the source-region and demonstrated that the ablation process is critical in the heating of the meteor's wake. In the head of the meteor, organic carbon appears to survive flash heating and rapid cooling. The temperatures in the wake of the meteor are just right for dissociation of CO and the formation of more complex organic compounds. The resulting materials could account for the bulk of pre-biotic organic carbon on the early Earth at the time of the origin of life.

Jenniskens, Peter↗

Testing of Error-Correcting Sparse Permutation Channel Codes

A computer program performs Monte Carlo direct numerical simulations for testing sparse permutation channel codes, which offer strong error-correction capabilities at high code rates and are considered especially suitable for storage of digital data in holographic and volume memories. A word in a code of this type is characterized by, among other things, a sparseness parameter (M) and a fixed number (K) of 1 or "on" bits in a channel block length of N.

Shcheglov, Kirill, V.↗

Monte Carlo Simulation to Estimate Likelihood of Direct Lightning Strikes

A software tool has been designed to quantify the lightning exposure at launch sites of the stack at the pads under different configurations. In order to predict lightning strikes to generic structures, this model uses leaders whose origins (in the x-y plane) are obtained from a 2D random, normal distribution.

Mata, Carlos↗

Applications of Low Density Flow Techniques and Catalytic Recombination at the Johnson Space Center

The talk presents a brief background on defInitions of catalysis and effects associated with chemically nonequilibrium and low-density flows of aerospace interest. Applications of catalytic recombination on surfaces in dissociated flow are given, including aero heating on reentry spacecraft thermal protection surfaces and reflection of plume flow on pressure distributions associated with the space station. Examples include aero heating predictions for the X-38 test vehicle, the inlet of a proposed gas-sampling probe used in high enthalpy test facilities, and a parabolic body at angle of attack. The effect of accommodation coefficients on thruster induced pressure distributions is also included. Examples of tools used include simple aero heating formulas based on boundary layer solutions, an engineering approximation that uses axisymmetric viscous shock layer flow to simulate full three dimensional flow, full computational fluid dynamics, and direct simulation Monte-Carlo calculations. Methods of determining catalytic recombination rates in arc jet flow are discus ed. An area of catalysis not fully understood is the formation of single-wall carbon nanotubes (SWNT) with gas phase or nano-size metal particles. The Johnson Space Center is making SWNTs using both a laser ablation technique and an electric arc vaporization technique.

Scott, Carl D.↗

Towards Plume Impingement Modeling in Space Environments

After 30 years human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to act as a home base for Lunar exploration. A large space station must be assembled while in orbit, using a “piecemeal” approach. In this context, it means that the different modules will arrive at different times and attach to what is already in service. The ISS provides an excellent example of this approach, and the proposed Lunar Gateway will undergo a similar assembly process. This assembly is achieved via “docking” maneuvers between modules, which are made possible by sequential firings of the onboard reaction control system (RCS) thrusters. They work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking approach. The issue is that the gas from these thrusters' forms flow structures described as “plumes” and can impinge onto the outer surfaces of the space station, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion. All mechanisms that can damage the space station and must be avoided. Both permanent and visiting modules will have these RCS thruster exhaust impingement problems. Accurately and efficiently modeling these plume is an involved multi-physics calculation but also an important tool when designing the control algorithms of approaching modules. This poster presents progress towards this simulation on two fronts. First is the verification of OpenFOAM for rarefied plume impingement calculations by direct comparison to published DAC cases. Second is the estimation of plume impingement strikes over the time scale of an entire docking event. This is done by using a simple plume source flow model and a prescribed motion visualizer—coded in Python. Together these tools help push NASA's capabilities for simulating these plume impingement effects.

Rarefied Flows↗

The Venus Atmospheric Sample Return (VATMOS-SR) Mission Concept

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small space- craft mission concept that would return several gas samples from the upper atmosphere of Venus (< 110 km) to Earth for scientific analysis. This could be the first sample return mission for an extraterrestrial atmosphere and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and ac- quire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen (e.g. trajectory shown in Fig. 1).

direct simulation Monte Carlo↗

The Venus Atmospheric Sample Return (VATMOS-SR) Mission Concept

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small space- craft mission concept that would return several gas samples from the upper atmosphere of Venus (< 110 km) to Earth for scientific analysis. This could be the first sample return mission for an extraterrestrial atmosphere and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and ac- quire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen.

direct simulation Monte Carlo↗

Realistic Covariance Generation for the GPM Spacecraft

We present several different methods for generating realistic predictive covariance, including Monte-Carlo simulations and more direct linear methods which require the addition of process noise. The Monte-Carlo simulation starts with an epoch uncertainty sample basis and propagates each trial to a time of interest in the future. The variance-covariance of the state elements as well as other higher order sample statistics can be readily computed from the propagated sample. While this method preserves the nonlinear effects on the propagated uncertainty, it is computationally intensive as a statistically significant sample size must be considered in the propagation process. Moreover, if the epoch covariance is optimistic, this can result in an underestimation of the prediction error. Another method is to propagate a state sensitivity matrix simultaneously with the satellite state, which allows the epoch covariances to be propagated forward in a linear fashion. This method does not preserve the non-linearity of the satellite state uncertainty but is much less computationally intensive. The propagated state covariance is the scaled to represent the realistic level of GPM state uncertainty via a "e-tuning process." The tuning process generates an inflation factor based on the observed error statistics of the predictive satellite trajectories when compared to the definitive ones. Difference tuning strategies are considered and compared via Goodness-of-Fit method testing for the Gaussian properties of the scaled covariance.

prediction error↗

An Application of the Direct Coulomb Electron Pair Production Process to the Energy Measurement of the "VH-Group" in the "Knee" Region of the "All-Particle" Energy Spectrum

The "all-particle" cosmic ray energy spectrum appears to be exhibiting a significant change in the spectral index just above approximately 3000 TeV. This could indicate (1) a change in the propagation of the cosmic rays in the galactic medium, and/or (2) the upper limit of the supernova shock wave acceleration mechanism, and/or (3) a new source of high-energy cosmic rays. Air shower and JACEE data indicate the spectral change is associated with a composition change to a heavier element mixture whereas DICE does not indicate this. A detector concept will be presented that utilizes the energy dependence of the production of direct Coulomb electron-positron pairs by energetic heavy ions. Monte Carlo simulations of a direct electron pair detector consisting of Pb target foils interleaved with planes of 1-mm square scintillating optical fibers will be discussed. The goal is to design a large area, non-saturating instrument to measure the energy spectrum of the individual cosmic ray elements in the "VH-group" for energies greater than 10 TeV/nucleon.

Derrickson, J. H.↗

Simulations of particle acceleration in parallel shocks: Direct comparison between Monte Carlo and one-dimensional hybrid codes

We have made a direct comparison between two different computer simulations of a plane, parallel, collisionless shock including particle acceleration to energies typical of those of diffuse ions observed at the earth bow shock. Despite the fact that the one-dimensional hybrid and Monte Carlo techniques employ entirely different algorithms, they give surprisingly close agreement in the overall shapes of the complete distribution functions for protons as well as heavier ions. Both methods show that energetic ions emerge smoothly from the background thermal plasma with approximately the same relative injection rate and that the fraction of the incoming plasma's energy flux that is converted into downstream enthalpy flux of the accelerated population (i.e., the acceleration efficiency) is similar in the two cases. The fraction of the downstream proton distribution made up of superthermal particles is quite large, with at least 10% of the energy flux going into protons with energies above 10 keV. In addition, an upstream precursor, produced by backstreaming energetic particles, is present in both shocks, although the Monte Carlo precursor is considerably longer than that produced in the hybrid shock. These results offer convincing evidence that, at least in these ways, the two simulations are consistent in their description of parallel shock structure and particle acceleration, and they lay the groundwork for development of shock models employing a combination of both methods.

Ellison, Donald C.↗

Probalistic Finite Elements (PFEM) structural dynamics and fracture mechanics

The purpose of this work is to develop computationally efficient methodologies for assessing the effects of randomness in loads, material properties, and other aspects of a problem by a finite element analysis. The resulting group of methods is called probabilistic finite elements (PFEM). The overall objective of this work is to develop methodologies whereby the lifetime of a component can be predicted, accounting for the variability in the material and geometry of the component, the loads, and other aspects of the environment; and the range of response expected in a particular scenario can be presented to the analyst in addition to the response itself. Emphasis has been placed on methods which are not statistical in character; that is, they do not involve Monte Carlo simulations. The reason for this choice of direction is that Monte Carlo simulations of complex nonlinear response require a tremendous amount of computation. The focus of efforts so far has been on nonlinear structural dynamics. However, in the continuation of this project, emphasis will be shifted to probabilistic fracture mechanics so that the effect of randomness in crack geometry and material properties can be studied interactively with the effect of random load and environment.

Liu, Wing-Kam↗

Estimation variance bounds of importance sampling simulations in digital communication systems

In practical applications of importance sampling (IS) simulation, two basic problems are encountered, that of determining the estimation variance and that of evaluating the proper IS parameters needed in the simulations. The authors derive new upper and lower bounds on the estimation variance which are applicable to IS techniques. The upper bound is simple to evaluate and may be minimized by the proper selection of the IS parameter. Thus, lower and upper bounds on the improvement ratio of various IS techniques relative to the direct Monte Carlo simulation are also available. These bounds are shown to be useful and computationally simple to obtain. Based on the proposed technique, one can readily find practical suboptimum IS parameters. Numerical results indicate that these bounding techniques are useful for IS simulations of linear and nonlinear communication systems with intersymbol interference in which bit error rate and IS estimation variances cannot be obtained readily using prior techniques.

Lu, D.↗