Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Experimental particle physics”

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 163 records · Page 9

Detailed Studies on Flame Extinction by Inert Particles in Normal- and Micro-gravity

The combustion of dusty flows has been studied to lesser extent than pure gas phase flows and sprays. Particles can have a strong effect by modifying the dynamic response and detailed structure of flames through the dynamic, thermal, and chemical couplings between the two phases. A rigorous understanding of the dynamics and structure of two-phase flows can be attained in stagnation flow configurations, which have been used by others to study spray combustion as well as reacting dusty flows. In earlier studies on reacting dusty flows, the thermal coupling between the two phases as well as the effect of gravity on the flame response were not considered. However, in Ref. 6, the thermal coupling between chemically inert particles and the gas was addressed in premixed flames. The effects of gravity was also studied showing that it can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature. The results showed a strong dynamic and thermal dependence of reacting dusty flows to particle number density. However, the work was only numerical and limited to twin-flames, stagnation, premixed flames. In Ref. 7 the effects of chemically inert particle clouds on the extinction of strained premixed and non-premixed flames were studied both experimentally and numerically at 1-g. It was shown and explained that large particles can cause more effective flame cooling compared to smaller particles. The effects of flame configuration and particle injection orientation were also addressed. The complexity of the coupling between the various parameters in such flows was demonstrated and it was shown that it was impossible to obtain a simple and still meaningful scaling that captured all the pertinent physics.

Andac, M. G.↗

Penetration experiments in aluminum 1100 targets using soda-lime glass projectiles

The cratering and penetration behavior of annealed aluminum 1100 targets, with thickness varied from several centimeters to ultra-thin foils less than 1 micrometer thick, were experimentally investigated using 3.2 mm diameter spherical soda-lime glass projectiles at velocities from 1 to 7 km/s. The objective was to establish quantitative, dimensional relationships between initial impact conditions (impact velocity, projectile diameter, and target thickness) and the diameter of the resulting crater or penetration hole. Such dimensional relationships and calibration experiments are needed to extract the diameters and fluxes of hypervelocity particles from space-exposed surfaces and to predict the performance of certain collisional shields. The cratering behavior of aluminum 1100 is fairly well predicted. However, crater depth is modestly deeper for our silicate impactors than the canonical value based on aluminum projectiles and aluminum 6061-T6 targets. The ballistic-limit thickness was also different. These differences attest to the great sensitivity of detailed crater geometry and penetration behavior on the physical properties of both the target and impactor. Each penetration experiment was equipped with a witness plate to monitor the nature of the debris plume emanating from the rear of the target. This plume consists of both projectile fragments and target debris. Both penetration hole and witness-plate spray patterns systematically evolve in response to projectile diameter/target thickness. The relative dimensions of the projectile and target totally dominate the experimental products documented in this report; impact velocity is an important contributor as well to the evolution of penetration holes, but is of subordinate significance for the witness-plate spray patterns.

Horz, Friedrich↗

Monodisperse Single-Material Granular Tribocharging Modeling and Experimental Validation

Charge transfer between insulating grains has been a topic of interest for many years as this phenomenon is extremely important to many areas in industry. Of particular interest to NASA is the behavior of electrically insulative dust grains such as those found on the lunar surface. Whether poured from a scoop during sample collection activities, agitated inside a drum via mining robots, or fluidized by a gas plume, particle-particle interactions between similar granular materials will be widespread on the Moon. Designers hope to have a predictive model for how dust grains are charged, transported, and deposited on flight systems to better understand issues electrostatics may cause in future missions. While difference in work function is often seen as a driver for charge transfer between two materials, this is an incomplete representation for insulative particles. Other asymmetries in the system can also lead to charge being transferred: one such being the difference between static vs dynamic particles. Experiments performed under vacuum in the Electrostatics and Surface Physics Laboratory (ESPL) clearly show that particles having more contacts with other particles tend to charge positively and particles having fewer contacts tend to charge negatively, despite being made from the same material. During these experiments, monodisperse (population of single diameter) spherical particles were held in a reservoir which was then elevated at one end to slope into a Faraday cup used to measure bulk charge of the grains. Care was taken to limit the number of spheres bouncing out of the cup, to isolate the beads from the container itself so that only grain-to-grain interactions occurred, and to assure neutrality of the beads before the experiment began. A discrete element method (DEM) modeling package already incorporating many granular mechanics interactions was augmented to include the electrostatic behaviors seen in experiments. The model can reproduce the experimental results above by tracking the number of contacts between pairs of individual grains and applying a charge transfer condition related to this contact parameter. The experimental results as well as their incorporation into the improved modeling suite will be discussed. Future work includes adding additional asymmetries into the experiments/model such as polydisperse populations of spherical particles and irregularly shaped grains.

granular↗

Monodisperse Single-Material Granular Tribocharging Modeling and Experimental Validation

Charge transfer between insulating grains has been a topic of interest for many years as this phenomenon is extremely important to many areas in industry. Of particular interest to NASA is the behavior of electrically insulative dust grains such as those found on the lunar surface. Whether poured from a scoop during sample collection activities, agitated inside a drum via mining robots, or fluidized by a gas plume, particle-particle interactions between similar granular materials will be widespread on the Moon. Designers hope to have a predictive model for how dust grains are charged, transported, and deposited on flight systems to better understand issues electrostatics may cause in future missions. While difference in work function is often seen as a driver for charge transfer between two materials, this is an incomplete representation for insulative particles. Other asymmetries in the system can also lead to charge being transferred: one such being the difference between static vs dynamic particles. Experiments performed under vacuum in the Electrostatics and Surface Physics Laboratory (ESPL) clearly show that particles having more contacts with other particles tend to charge positively and particles having fewer contacts tend to charge negatively, despite being made from the same material. During these experiments, monodisperse (population of single diameter) spherical particles were held in a reservoir which was then elevated at one end to slope into a Faraday cup used to measure bulk charge of the grains. Care was taken to limit the number of spheres bouncing out of the cup, to isolate the beads from the container itself so that only grain-to-grain interactions occurred, and to assure neutrality of the beads before the experiment began. A discrete element method (DEM) modeling package already incorporating many granular mechanics interactions was augmented to include the electrostatic behaviors seen in experiments. The model can reproduce the experimental results above by tracking the number of contacts between pairs of individual grains and applying a charge transfer condition related to this contact parameter. The experimental results as well as their incorporation into the improved modeling suite will be discussed. Future work includes adding additional asymmetries into the experiments/model such as polydisperse populations of spherical particles and irregularly shaped grains.

electrostatics↗

Traveling-Wave Tube Amplifier Model to Predict High-Order Modulation Intersymbol Interference

Demands for increased data rates in satellite communications necessitate higher order modulation schemes, larger system bandwidth, and minimum distortion of the modulated signal as it is passed through the traveling wave tube amplifier (TWTA). One type of distortion that the TWTA contributes to is intersymbol interference (ISI), and this becomes particularly disruptive with wide-band, complex modulation schemes. It is suspected that in addition to the dispersion of the TWT, frequency dependent reflections due to mismatches within the TWT are a significant contributor to ISI. To experimentally investigate the effect of these mismatches within the physical TWT on ISI would be prohibitively expensive, as it would require manufacturing numerous amplifiers in addition to the acquisition of the required digital hardware. In an attempt to develop a more accurate model to correlate IS1 with the TWTA and the operational signal, a fully three-dimensional (3D), time-dependent, TWT interaction model has been developed using the electromagnetic particle-in-cell (PIC) code MAFIA (solution of Maxwell's equations by the Finite-Integration-Algorithm). The model includes a user defined slow-wave circuit with a spatially tapered region of loss to implement a sever, and spatially varied geometry (such as helical pitch) to implement a phase velocity taper. The model also includes user defined input/output coupling and an electron beam contained by solenoidal, electrostatic, or periodic permanent magnet (PPM) focusing allowing standard or novel TWTs to be investigated. This model comprehensively takes into account the effects of frequency dependent nonlinear distortions (MAM and AMPM); gain ripple due to frequency dependent reflections at the input/output coupling, severs, and mismatches from dynamic pitch variations; drive induced oscillations; harmonic generation; intermodulation products; and backward waves.

Kory, Carol L.↗

Plans for Extreme Energy Cosmic Ray Observations from Space

Cosmic rays have been detected at energies beyond 10(exp 20) eV, where Universe is predicted to become opaque to protons. The acceleration of cosmic rays to such extreme energies in known astrophysical objects has also proven difficult to understand, leading to many suggestions that new physics may be required to explain their existence. This has prompted the construction of new experiments designed to detect cosmic rays with fluxes below 1 particle/km/century and follow their spectrum to even higher energies. To detect large numbers of these particles, the next generation of these experiments must be performed on space-based platforms that look on very large detection volumes in the Earth's atmosphere. The talk will review the experimental and theoretical investigations of extreme energy cosmic rays and discuss the present and planned experiments to extend measurements beyond 10(exp 21) eV.

Adams, James H., Jr.↗

Lunary Dust Grain Charging by Electron Impact: Complex Role of Secondary Electron Emissions in Space Environments

Dust grains in various astrophysical environments are generally charged electrostatically by photoelectric emissions with radiation from nearby sources, or by electron/ion collisions by sticking or secondary electron emissions (SEES). The high vacuum environment on the lunar surface leads to some unusual physical and dynamical phenomena involving dust grains with high adhesive characteristics, and levitation and transportation over long distances. Knowledge of the dust grain charges and equilibrium potentials is important for understanding a variety of physical and dynamical processes in the interstellar medium, and heliospheric, interplanetary/ planetary, and lunar environments. It has been well recognized that the charging properties of individual micron-/submicron-size dust grains are expected to be substantially different from the corresponding values for bulk materials. In this paper, we present experimental results on the charging of individual 0.2-13 m size dust grains selected from Apollo 11 and 17 dust samples, and spherical silica particles by exposing them to mono-energetic electron beams in the 10-200 eV energy range. The dust charging process by electron impact involving the SEES discussed is found to be a complex charging phenomenon with strong particle size dependence. The measurements indicate substantial differences between the polarity and magnitude of the dust charging rates of individual small-size dust grains, and the measurements and model properties of corresponding bulk materials. A more comprehensive plan of measurements of the charging properties of individual dust grains for developing a database for realistic models of dust charging in astrophysical and lunar environments is in progress.

Abbas, M. M.↗

Numerical and Analytical Model of an Electrodynamic Dust Shield for Solar Panels on Mars

Masuda and collaborators at the University of Tokyo developed a method to confine and transport particles called the electric curtain in which a series of parallel electrodes connected to an AC source generates a traveling wave that acts as a contactless conveyor. The curtain electrodes can be excited by a single-phase or a multi-phase AC voltage. A multi-phase curtain produces a non-uniform traveling wave that provides controlled transport of those particles [1-6]. Multi-phase electric curtains from two to six phases have been developed and studied by several research groups [7-9]. We have developed an Electrodynamic Dust Shield prototype using threephase AC voltage electrodes to remove dust from surfaces. The purpose of the modeling work presented here is to research and to better understand the physics governing the electrodynamic shield, as well as to advance and to support the experimental dust shield research.

Calle, C. I.↗

Overview of the Graphical User Interface for the GERMcode (GCR Event-Based Risk Model)

The descriptions of biophysical events from heavy ions are of interest in radiobiology, cancer therapy, and space exploration. The biophysical description of the passage of heavy ions in tissue and shielding materials is best described by a stochastic approach that includes both ion track structure and nuclear interactions. A new computer model called the GCR Event-based Risk Model (GERM) code was developed for the description of biophysical events from heavy ion beams at the NASA Space Radiation Laboratory (NSRL). The GERMcode calculates basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at NSRL for the purpose of simulating space radiobiological effects. For mono-energetic beams, the code evaluates the linear-energy transfer (LET), range (R), and absorption in tissue equivalent material for a given Charge (Z), Mass Number (A) and kinetic energy (E) of an ion. In addition, a set of biophysical properties are evaluated such as the Poisson distribution of ion or delta-ray hits for a specified cellular area, cell survival curves, and mutation and tumor probabilities. The GERMcode also calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle. The contributions from primary ion and nuclear secondaries are evaluated. The GERMcode accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections, and has been used by the GERMcode for application to thick target experiments. The GERMcode provides scientists participating in NSRL experiments with the data needed for the interpretation of their experiments, including the ability to model the beam line, the shielding of samples and sample holders, and the estimates of basic physical and biological outputs of the designed experiments. We present an overview of the GERMcode GUI, as well as providing training applications.

Kim, Myung-Hee Y.↗

Rate equation for desorbing particles

A rate equation describing the desorption of particles from substrates has been formulated using a simple classical model in terms of the particulate escape probabilities and escape velocities. Numerical results for desorption rates as a function of temperature using simple two-body potentials are calculated for various metallic systems. These results were obtained for the low coverage limits and are found to be in good agreement with experimental findings. The present formulation also provides an improved picture of the physical nature of the commonly discussed parameters, such as the frequency factor, desorption energy, and desorption temperature in relation to observed desorption phenomena.

Halicioglu, T.↗

Overview of the Graphical User Interface for the GERM Code (GCR Event-Based Risk Model

The descriptions of biophysical events from heavy ions are of interest in radiobiology, cancer therapy, and space exploration. The biophysical description of the passage of heavy ions in tissue and shielding materials is best described by a stochastic approach that includes both ion track structure and nuclear interactions. A new computer model called the GCR Event-based Risk Model (GERM) code was developed for the description of biophysical events from heavy ion beams at the NASA Space Radiation Laboratory (NSRL). The GERM code calculates basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at NSRL for the purpose of simulating space radiobiological effects. For mono-energetic beams, the code evaluates the linear-energy transfer (LET), range (R), and absorption in tissue equivalent material for a given Charge (Z), Mass Number (A) and kinetic energy (E) of an ion. In addition, a set of biophysical properties are evaluated such as the Poisson distribution of ion or delta-ray hits for a specified cellular area, cell survival curves, and mutation and tumor probabilities. The GERM code also calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle. The contributions from primary ion and nuclear secondaries are evaluated. The GERM code accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections, and has been used by the GERM code for application to thick target experiments. The GERM code provides scientists participating in NSRL experiments with the data needed for the interpretation of their experiments, including the ability to model the beam line, the shielding of samples and sample holders, and the estimates of basic physical and biological outputs of the designed experiments. We present an overview of the GERM code GUI, as well as providing training applications.

Kim, Myung-Hee↗

Measurement Capabilities of Planar Doppler Velocimetry in Large-Scale Wind Tunnels

Over the past few years, Planar Doppler Velocimetry (PDT) has been shown by several laboratories to offer an attractive means for measuring three-dimensional velocity vectors everywhere in a light sheet placed in a flow. Unlike some other optical means of measuring flow velocities, PDT is particularly attractive for use in large wind tunnels where distances to the sample region may be several meters, because it does not require the spatial resolution and tracking of individual scattering particles or the alignment of crossed beams at large distances. To date, demonstrations of PDT (also called Doppler Global Velocimetry by some authors) have been made either in low speed flows without quantitative comparison to other measurements, or in supersonic flows where the Doppler shift is large and its measurement is relatively insensitive to instrumental errors. Moreover, most reported applications have relied on the use of continuous-wave lasers, which limit the measurement to time-averaged velocity fields. This work summarizes the results of two previous studies of PDT in which the use of pulsed lasers to obtain instantaneous velocity vector fields is evaluated. The objective has been to quantitatively define and demonstrate PDT capabilities for applications in large-scale wind tunnels that are intended primarily for the testing of rotorcraft and subsonic aircraft at speeds typically less than 100 m/s. For such applications, the adequate resolution of low-speed flow fields requires accurate measurements of small Doppler shifts that are obtained at distances of several meters from the sample region and with a field of view that is sufficient to encompass the entire region of interest. The use of pulsed lasers provides the unique capability to obtain not only time-averaged fields, but also their statistical fluctuation amplitudes and the spatial excursions of unsteady flow regions such as wakes, separations, and rotor-tip vortices. To accomplish the objectives of these studies, the PDT measurement process was first modeled and its performance evaluated computationally. The noise sources considered included those related to the optical and electronic properties of Charge-Coupled Device (CCD) arrays and to speckle effects associated with the coherent illumination of aerosols from pulsed lasers. The signal noise estimates were incorporated into the PDT signal analysis process and combined with computed scattering signals using a Mae scattering theory for a distributed range of aerosol particle sizes. The results are used to define the necessary instrument configuration, to show that the expected signal levels from a practical PDT system are sufficiently large to allow its useful application in large facilities, and to show that the expected velocity measurement uncertainties are small compared to the mean velocities of interest for most subsonic, large-scale wind tunnel testing, Experimental studies using several bench-top setups are then described that validate the physics of the PDT model and demonstrate the ability to obtain accurate PDT measurements using procedures that are compatible with large-scale wind tunnel operations. The PDT measurement capabilities found in this study lead to the conclusion that PDT offers significant advantages compared to other means of measuring velocity fields in large-scale wind tunnels, including those operating at low speeds.

McKenzie, Robert L.↗

Experimenter's data package for the descending layers rocket

In response to a proposal from Science Applications International Corporation (SAIC), NASA Headquarters has approved a sounding rocket mission designed to study the physics of intermediate layers in the Earth's ionosphere at middle latitudes. The experiment will be carried out by a team of scientists and engineers from the NASA Wallops Flight Facility, SAIC, the NASA Goddard Space Flight Center, and the Millstone Hill radar observatory. The mission will involve the launch of an instrumented sounding rocket from the Wallops Island rocket range in the summer of 1994, with the objective of penetrating a descending ionized layer in the E-region between altitudes of 115 and 140 km. Instrumentation aboard the rocket will measure the ion and neutral composition of the layer, its plasma density, driving wind and electric field forces, the thermal ion distribution function, and electron temperature. Depending on payload weight constraints and subject to availability, a particle detector to measure energetic ion and/or electron fluxes near the layer may also be included. This document was prepared as a reference for the NASA payload development and experiment teams, for distribution at the Project Initiation Conference (PIC). The design specifications discussed herein are therefore of a preliminary nature; the intent is to promote open discussions between experimenters and NASA engineers that will lead to a final design capable of achieving the experiment objectives.

Earle, Greg↗

Search For Effects of an Electrostatic Potential on Clocks in the Frame of Reference of a Charged Particle

The Program Planning meeting was held on Jan 27-28 at the Washington University, Dept. of Physics. Experiment strategies and details were discussed and milestones mapped. Dr. Ringermacher, the program technical leader and Co-Principal Investigator, presented a Physics Department seminar on the theory and experiment in progress. Following the meeting Professor Conradi (CO-PI) and his students designed the NMR probe and electronics required for the experiment and submitted the drawings to the shop for fabrication. Dr. Cassenti, at UTRC, initiated analytic and computer analysis of the experimental conditions in order to anticipate possible classical effects. Dr. Ringermacher, at KRONOTRAN, re-evaluated the experiment interpretation from the theory and continued working to secure its underpinnings.

Ringermacher, Harry I.↗

Volcanic aerosols: Chemistry, evolution, and effects

Stratospheric aerosols have been the subject of scientific speculation since the 1880s, when the powerful eruption of Krakatoa attracted worldwide attention to the upper atmosphere through spectacular optical displays. The presence of a permanent tenuous dust layer in the lower stratosphere was postulated in the 1920s following studies of the twilight glow. Junge collected the first samples of these 'dust' particles and demonstrated that they were actually composed of sulfates, most likely concentrated sulfuric acid (Junge and Manson, 1961; Junge, 1963). Subsequent research has been spurred by the realization that stratospheric particles can influence the surface climate of earth through their effects on atmospheric radiation. Such aerosols can also influence, through chemical and physical effects, the trace composition of the atmosphere, ozone concentrations, and atmospheric electrical properties. The properties of stratospheric aerosols (both the background particles and those enhanced by volcanic eruptions) were measured in situ by balloon ascents and high altitude aircraft sorties. The aerosols were also observed remotely from the ground and from satellites using both active (lidar) and passive (solar occultation) techniques (remote sensing instruments were carried on aircraft and balloon platforms as well). In connection with the experimental work, models were developed to test theories of particle formation and evolution, to guide measurement strategies, to provide a means of connecting laboratory and field data, and to apply the knowledge gained to answer practical questions about global changes in climate, depletion of the ozone layer, and related environmental problems.

Turco, Richard↗

Positron Interactions with Atoms and Ions

Dirac, in 1928, combining the ideas of quantum mechanics and the ideas of relativity invented the well-known relativistic wave equation. In his formulation, he predicted an antiparticle of the electron of spin n-bar/2. He thought that this particle must be a proton. Dirac published his interpretation in a paper 'A theory of electrons and protons.' It was shown later by the mathematician Hermann Weyl that the Dirac theory was completely symmetric between negative and positive particles and the positive particle must have the same mass as that of the electron. In his J. Robert Oppenheimer Memorial Prize Acceptance Speech, Dirac notes that 'Blackett was really the first person to obtain hard evidence for the existence of a positron but he was afraid to publish it. He wanted confirmation, he was really over cautious.' Positron, produced by the collision of cosmic rays in a cloud chamber, was detected experimentally by Anderson in 1932. His paper was published in Physical Review in 1933. The concept of the positron and its detection were the important discoveries of the 20th century. I have tried to discuss various processes involving interactions of positrons with atoms and ions. This includes scattering, bound states and resonances. It has not been possible to include the enormous work which has been carried out during the last 40 or 50 years in theory and measurements.

Bhatia, Anand K.↗

3-D Flow Field Diagnostics and Validation Studies using Stereoscopic Tracking Velocimetry

The measurement of 3-D three-component velocity fields is of great importance in both ground and space experiments for understanding materials processing and fluid physics. Here, we present the investigation results of stereoscopic tracking velocimetry (STV) for measuring 3-D velocity fields. The effort includes diagnostic technology development, experimental velocity measurement, and comparison with analytical and numerical computation. The advantages of STV stems from the system simplicity for building compact hardware and in software efficiency for continual near-real-time process monitoring. It also has illumination flexibility for observing volumetric flow fields from arbitrary directions. STV is based on stereoscopic CCD observations of particles seeded in a flow. Neural networks are used for data analysis. The developed diagnostic tool is tested with a simple directional solidification apparatus using Succinonitrile. The 3-D velocity field in the liquid phase is measured and compared with results from detailed numerical computations. Our theoretical, numerical, and experimental effort has shown STV to be a viable candidate for reliably quantifying the 3-D flow field in materials processing and fluids experiments.

Cha, Soyoung Stephen↗

Predictive Modeling in Plasma Reactor and Process Design

Research continues toward the improvement and increased understanding of high-density plasma tools. Such reactor systems are lauded for their independent control of ion flux and energy enabling high etch rates with low ion damage and for their improved ion velocity anisotropy resulting from thin collisionless sheaths and low neutral pressures. Still, with the transition to 300 mm processing, achieving etch uniformity and high etch rates concurrently may be a formidable task for such large diameter wafers for which computational modeling can play an important role in successful reactor and process design. The inductively coupled plasma (ICP) reactor is the focus of the present investigation. The present work attempts to understand the fundamental physical phenomena of such systems through computational modeling. Simulations will be presented using both computational fluid dynamics (CFD) techniques and the direct simulation Monte Carlo (DSMC) method for argon and chlorine discharges. ICP reactors generally operate at pressures on the order of 1 to 10 mTorr. At such low pressures, rarefaction can be significant to the degree that the constitutive relations used in typical CFD techniques become invalid and a particle simulation must be employed. This work will assess the extent to which CFD can be applied and evaluate the degree to which accuracy is lost in prediction of the phenomenon of interest; i.e., etch rate. If the CFD approach is found reasonably accurate and bench-marked with DSMC and experimental results, it has the potential to serve as a design tool due to the rapid time relative to DSMC. The continuum CFD simulation solves the governing equations for plasma flow using a finite difference technique with an implicit Gauss-Seidel Line Relaxation method for time marching toward a converged solution. The equation set consists of mass conservation for each species, separate energy equations for the electrons and heavy species, and momentum equations for the gas. The sheath is modeled by imposing the Bohm velocity to the ions near the walls. The DSMC method simulates each constituent of the gas as a separate species which would be analogous in CFD to employing separate species mass, momentum, and energy equations. All particles including electrons are moved and allowed to collide with one another with the stipulation that the electrons remain tied to the ions consistent with the concept of ambipolar diffusion. The velocities of the electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. These benefits come at a price in terms of computational time and memory. The DSMC and CFD are made as consistent as possible by using similar chemistry and power deposition models. Although the comparison of CFD and DSMC is interesting, the main goal of this work is the increased understanding of high-density plasma flowfields that can then direct improvements in both techniques. This work is unique in the level of the physical models employed in both the DSMC and CFD for high-density plasma reactor applications. For example, the electrons are simulated in the present DSMC work which has not been done before for low temperature plasma processing problems. In the CFD approach, for the first time, the charged particle transport (discharge physics) has been self-consistently coupled to the gas flow and heat transfer.

Hash, D. B.↗