Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

At least 1,675 records · Page 93

Numerical Propulsion System Simulation (NPSS) 1999 Industry Review

The technologies necessary to enable detailed numerical simulations of complete propulsion systems are being developed at the NASA Glenn Research Center in cooperation with industry, academia, and other government agencies. Large scale, detailed simulations will be of great value to the nation because they eliminate some of the costly testing required to develop and certify advanced propulsion systems. In addition, time and cost savings will be achieved by enabling design details to be evaluated early in the development process before a commitment is made to a specific design. This concept is called the Numerical Propulsion System Simulation (NPSS). NPSS consists of three main elements: (1) engineering models that enable multidisciplinary analysis of large subsystems and systems at various levels of detail, (2) a simulation environment that maximizes designer productivity, and (3) a cost-effective, high-performance computing platform. A fundamental requirement of the concept is that the simulations must be capable of overnight execution on easily accessible computing platforms. This will greatly facilitate the use of large-scale simulations in a design environment. This paper describes the current status of the NPSS with specific emphasis on the progress made over the past year on air breathing propulsion applications. In addition, the paper contains a summary of the feedback received from industry partners in the development effort and the actions taken over the past year to respond to that feedback. The NPSS development was supported in FY99 by the High Performance Computing and Communications Program.

Lytle, John↗

Recommendations for Achieving Accurate Numerical Simulation of Tip Clearance Flows in Transonic Compressor Rotors

The tip clearance flows of transonic compressor rotors are important because they have a significant impact on rotor and stage performance. While numerical simulations of these flows are quite sophisticated. they are seldom verified through rigorous comparisons of numerical and measured data because these kinds of measurements are rare in the detail necessary to be useful in high-speed machines. In this paper we compare measured tip clearance flow details (e.g. trajectory and radial extent) with corresponding data obtained from a numerical simulation. Recommendations for achieving accurate numerical simulation of tip clearance flows are presented based on this comparison. Laser Doppler Velocimeter (LDV) measurements acquired in a transonic compressor rotor, NASA Rotor 35, are used. The tip clearance flow field of this transonic rotor was simulated using a Navier-Stokes turbomachinery solver that incorporates an advanced k-epsilon turbulence model derived for flows that are not in local equilibrium. Comparison between measured and simulated results indicates that simulation accuracy is primarily dependent upon the ability of the numerical code to resolve important details of a wall-bounded shear layer formed by the relative motion between the over-tip leakage flow and the shroud wall. A simple method is presented for determining the strength of this shear layer.

VanZante, Dale E.↗

CFD Simulation of Liquid Rocket Engine Injectors

Detailed design issues associated with liquid rocket engine injectors and combustion chamber operation require CFD methodology which simulates highly three-dimensional, turbulent, vaporizing, and combusting flows. The primary utility of such simulations involves predicting multi-dimensional effects caused by specific injector configurations. SECA, Inc. and Engineering Sciences, Inc. have been developing appropriate computational methodology for NASA/MSFC for the past decade. CFD tools and computers have improved dramatically during this time period; however, the physical submodels used in these analyses must still remain relatively simple in order to produce useful results. Simulations of clustered coaxial and impinger injector elements for hydrogen and hydrocarbon fuels, which account for real fluid properties, is the immediate goal of this research. The spray combustion codes are based on the FDNS CFD code' and are structured to represent homogeneous and heterogeneous spray combustion. The homogeneous spray model treats the flow as a continuum of multi-phase, multicomponent fluids which move without thermal or velocity lags between the phases. Two heterogeneous models were developed: (1) a volume-of-fluid (VOF) model which represents the liquid core of coaxial or impinger jets and their atomization and vaporization, and (2) a Blob model which represents the injected streams as a cloud of droplets the size of the injector orifice which subsequently exhibit particle interaction, vaporization, and combustion. All of these spray models are computationally intensive, but this is unavoidable to accurately account for the complex physics and combustion which is to be predicted, Work is currently in progress to parallelize these codes to improve their computational efficiency. These spray combustion codes were used to simulate the three test cases which are the subject of the 2nd International Workshop on-Rocket Combustion Modeling. Such test cases are considered by these investigators to be very valuable for code validation because combustion kinetics, turbulence models and atomization models based on low pressure experiments of hydrogen air combustion do not adequately verify analytical or CFD submodels which are necessary to simulate rocket engine combustion. We wish to emphasize that the simulations which we prepared for this meeting are meant to test the accuracy of the approximations used in our general purpose spray combustion models, rather than represent a definitive analysis of each of the experiments which were conducted. Our goal is to accurately predict local temperatures and mixture ratios in rocket engines; hence predicting individual experiments is used only for code validation. To replace the conventional JANNAF standard axisymmetric finite-rate (TDK) computer code 2 for performance prediction with CFD cases, such codes must posses two features. Firstly, they must be as easy to use and of comparable run times for conventional performance predictions. Secondly, they must provide more detailed predictions of the flowfields near the injector face. Specifically, they must accurately predict the convective mixing of injected liquid propellants in terms of the injector element configurations.

Farmer, Richard↗

Vertical Axis Rotational Motion Cues in Hovering Flight Simulation

A previous study that examined how yaw motion affected a pilot's ability to perform realistic hovering flight tasks indicated that any amount of pure yaw motion had little-to-no effect on pilot performance or opinion. In that experiment, pilots were located at the vehicle's center of rotation; thus lateral or longitudinal accelerations were absent. The purpose of the new study described here was to investigate further these unanticipated results for additional flight tasks, but with the introduction of linear accelerations associated with yaw rotations when the pilot is not at the center of rotation. The question of whether a yaw motion degree-of-freedom is necessary or not is important to government regulators who specify what simulator motions are necessary according to prescribed levels of simulator sophistication. Currently, specifies two levels of motion sophistication for flight simulators: full 6-degree-of-freedom and 3-degree-of-freedom. For the less sophisticated simulator, the assumed three degrees of freedom are pitch, roll, and heave. If other degrees of freedom are selected, which are different f rom these three, they must be qualified on a case-by-case basis. Picking the assumed three axes is reasonable and based upon experience, but little empirical data are available to support the selection of critical axes. Thus, the research described here is aimed at answering this question. The yaw and lateral degrees of freedom were selected to be examined first, and maneuvers were defined to uncouple these motions from changes in the gravity vector with respect to the pilot. This approach simplifies the problem to be examined. For this experiment, the NASA Ames Vertical Motion Simulator was used in a comprehensive investigation. The math model was an AH-64 Apache in hover, which was identified from flight test data and had previously been validated by several AH-64 pilots. The pilot's head was located 4.5 ft in front of the vehicle center of gravity, which is representative of the AH-64 pilot location. Six test pilots flew three tasks that were specifically designed to represent a broad class of situations in which both lateral and yaw motion cues may be useful. For the first task, the pilot controlled only the yaw axis and was required to rapidly acquire a North heading from 15 deg yaw offsets to either the East or West. This task allowed for full, or 1:1, motion to be used in all axes (yaw, lateral, and longitudinal). The second task was a 10 sec., 180 deg. pedal turn over a runway, but with the pilot only controlling the yaw degree-of-freedom. The position of the vehicle's center-of-mass remained fixed. This maneuver was taken from a current U.S. Army rotary wing design standard5 and is representative of a maneuver performed for acceptance of military helicopters; however, it does not allow for full 1:1 motion, since the simulator cab cannot rotate 180 deg. The third task required the pilot to perform a rapid 9 ft climb at a constant heading. This task was challenging, because rapid collective lever movement in the unaugmented AH64 results in a substantial yawing moment (due to engine torque) that must be countered by the pilot. This task also had full motion in all axes, but, in this case, the pilot had two axes to control simultaneously, rather than one as in the previous tasks. Four motion configurations were examined for each task: full motion (except for the 180 deg turn, for which the motion system was configured to provide as much motion as possible), full linear with no yaw motion, full yaw with no linear motion, and no motion. Each configuration was flown four times in a randomized test matrix, and the pilots were not informed of the configuration given. Vehicle state data were recorded for objective performance comparisons, and pilots provided subjective comments and ratings. As part of the pilots' evaluation, they were asked to rate the compensation required, the overall fidelity of the motion as compared to real flight, and whether motion was detected or not in each of the six degrees of freedom. In addition, the pilots provided a numerical level-of confidence rating, between 1 and 7, corresponding to how sure they were whether or not motion was present in each degree-of-freedom. The latter rating allow classical signal detection analysis to be performed.

Schroeder, Jeffrey A.↗

Textbook Multigrid Efficiency for Computational Fluid Dynamics Simulations

Considerable progress over the past thirty years has been made in the development of large-scale computational fluid dynamics (CFD) solvers for the Euler and Navier-Stokes equations. Computations are used routinely to design the cruise shapes of transport aircraft through complex-geometry simulations involving the solution of 25-100 million equations; in this arena the number of wind-tunnel tests for a new design has been substantially reduced. However, simulations of the entire flight envelope of the vehicle, including maximum lift, buffet onset, flutter, and control effectiveness have not been as successful in eliminating the reliance on wind-tunnel testing. These simulations involve unsteady flows with more separation and stronger shock waves than at cruise. The main reasons limiting further inroads of CFD into the design process are: (1) the reliability of turbulence models; and (2) the time and expense of the numerical simulation. Because of the prohibitive resolution requirements of direct simulations at high Reynolds numbers, transition and turbulence modeling is expected to remain an issue for the near term. The focus of this paper addresses the latter problem by attempting to attain optimal efficiencies in solving the governing equations. Typically current CFD codes based on the use of multigrid acceleration techniques and multistage Runge-Kutta time-stepping schemes are able to converge lift and drag values for cruise configurations within approximately 1000 residual evaluations. An optimally convergent method is defined as having textbook multigrid efficiency (TME), meaning the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in the discretized system of equations (residual equations). In this paper, a distributed relaxation approach to achieving TME for Reynolds-averaged Navier-Stokes (RNAS) equations are discussed along with the foundations that form the basis of this approach. Because the governing equations are a set of coupled nonlinear conservation equations with discontinuities (shocks, slip lines, etc.) and singularities (flow- or grid-induced), the difficulties are many. This paper summarizes recent progress towards the attainment of TME in basic CFD simulations.

Brandt, Achi↗

A Methodology for Evaluating the Fidelity of Ground-Based Flight Simulators

An analytical and experimental investigation was undertaken to model the manner in which pilots perceive and utilize visual, proprioceptive, and vestibular cues in a ground-based flight simulator. The study was part of a larger research effort which has the creation of a methodology for determining flight simulator fidelity requirements as its ultimate goal. The study utilized a closed-loop feedback structure of the pilot/simulator system which included the pilot, the cockpit inceptor, the dynamics of the simulated vehicle and the motion system. With the exception of time delays which accrued in visual scene production in the simulator, visual scene effects were not included in this study. The NASA Ames Vertical Motion Simulator was used in a simple, single-degree of freedom rotorcraft bob-up/down maneuver. Pilot/vehicle analysis and fuzzy-inference identification were employed to study the changes in fidelity which occurred as the characteristics of the motion system were varied over five configurations. The data from three of the five pilots that participated in the experimental study were analyzed in the fuzzy-inference identification. Results indicate that both the analytical pilot/vehicle analysis and the fuzzy-inference identification can be used to reflect changes in simulator fidelity for the task examined.

Zeyada, Y.↗

The Impact on Simulated Storm Structure and Intensity of Variations in the Lifted Condensation Level and the Level of Free Convection

The sensitivities of convective storm structure and intensity to changes in the altitudes of the prestorm environmental lifted condensation level and level of free convection axe studied using a full-physics three-dimensional cloud model. Matrices of simulations are conducted for a range of LCL=LFC altitudes, using a single moderately-sheared curved hodograph trace in conjunction with 1 convective available potential energy values of either 800 or 2000 J/kg, with the matrices consisting of all four combinations of two distinct choices of buoyancy and shear profile shape. For each value of CAPE, the LCL=LFC altitudes are also allowed to vary in a series of simulations based on the most highly compressed buoyancy and shear profiles for that CAPE, with the environmental buoyancy profile shape, subcloud equivalent potential temperature, subcloud lapse rates of temperature and moisture, and wind profile held fixed. For each CAPE, one final simulation is conducted using a near optimal LFC, but a lowered LCL, with a neutrally buoyant environmental thermal profile specified in between. Results show that, for the buoyancy-starved small-CAPE environments, the simulated storms are supercells and are generally largest and most intense when LCL=LFC altitudes lie in the approximate range 1.5-2.5 km above the surface. The simulations show similar trends for the shear-starved large-CAPE environments, except that conversion from supercell to multicell morphology frequently occurs when the LCL is high. For choices of LCL=LFC height within the optimal 1.5-2.5 km range, peak storm updraft overturning efficiency may approaches unity relative to parcel theory, while for lower LCL=LFC heights, overturning efficiency is reduced significantly. The enhancements of overturning efficiency and updraft diameter with increasing LFC height are shown to be the result of systematic increases in the mean equivalent potential temperature of the updraft at cloud base. For the shear-starved environments, the tendency for outflow dominance is eliminated, but a large overturning efficiency maintained, when a low LCL is used in conjunction with a high LFC. The result regarding outflow dominance at high LCL is consistent with expectations, but the beneficial effect of a high LFC on convective overturning efficiency has not previously been widely recognized. The simulation findings here also appear to be consistent with statistics from previous severe storm environment climatologies, but provide a new framework for interpreting those statistics.

McCaul, Eugene W., Jr.↗

Forced and Free Intra-Seasonal Variability Over the South Asian Monsoon Region Simulated by 10 AGCMs

This study examines intra-seasonal (20-70 day) variability in the South Asian monsoon region during 1997/98 in ensembles of 10 simulations with 10 different atmospheric general circulation models. The 10 ensemble members for each model are forced with the same observed weekly sea surface temperature (SST) but differ from each other in that they are started from different initial atmospheric conditions. The results show considerable differences between the models in the simulated 20-70 day variability, ranging from much weaker to much stronger than the observed. A key result is that the models do produce, to varying degrees, a response to the imposed weekly SST. The forced variability tends to be largest in the Indian and western Pacific Oceans where, for some models, it accounts for more than 1/4 of the 20-70 day intra-seasonal variability in the upper level velocity potential during these two years. A case study of a strong observed MJO (intraseasonal oscillation) event shows that the models produce an ensemble mean eastward propagating signal in the tropical precipitation field over the Indian Ocean and western Pacific, similar to that found in the observations. The associated forced 200 mb velocity potential anomalies are strongly phase locked with the precipitation anomalies, propagating slowly to the east (about 5 m/s) with a local zonal wave number two pattern that is generally consistent with the developing observed MJO. The simulated and observed events are, however, approximately in quadrature, with the simulated response 2 leading by 5-10 days. The phase lag occurs because, in the observations, the positive SST anomalies develop upstream of the main convective center in the subsidence region of the MJO, while in the simulations, the forced component is in phase with the SST. For all the models examined here, the intraseasonal variability is dominated by the free (intra-ensemble) component. The results of our case study show that the free variability has a predominately zonal wave number one pattern, and has propagation speeds (10 - 15 m/s) that are more typical of observed MJO behavior away from the convectively active regions. The free variability appears to be synchronized with the forced response, at least, during the strong event examined here. The results of this study support the idea that coupling with SSTs plays an important, though probably not dominant, role in the MJO. The magnitude of the atmospheric response to the SST appears to be in the range of 15% - 30% of the 20-70 day variability over much of the tropical eastern Indian and western Pacific Oceans. The results also highlight the need to use caution when interpreting atmospheric model simulations in which the prescribed SST resolve MJO time scales.

Wu, Man Li C.↗

A Unique Software System For Simulation-to-Flight Research

"Simulation-to-Flight" is a research development concept to reduce costs and increase testing efficiency of future major aeronautical research efforts at NASA. The simulation-to-flight concept is achieved by using common software and hardware, procedures, and processes for both piloted-simulation and flight testing. This concept was applied to the design and development of two full-size transport simulators, a research system installed on a NASA B-757 airplane, and two supporting laboratories. This paper describes the software system that supports the simulation-to-flight facilities. Examples of various simulation-to-flight experimental applications were also provided.

Chung, Victoria I.↗

The NASA Ames Hypersonic Combustor-Model Inlet CFD Simulations and Experimental Comparisons

Computations have been performed on a three-dimensional inlet associated with the NASA Ames combustor model for the hypersonic propulsion experiment in the 16-inch shock tunnel. The 3-dimensional inlet was designed to have the combustor inlet flow nearly two-dimensional and of sufficient mass flow necessary for combustion. The 16-inch shock tunnel experiment is a short duration test with test time of the order of milliseconds. The flow through the inlet is in chemical non-equilibrium. Two test entries have been completed and limited experimental results for the inlet region of the combustor-model are available. A number of CFD simulations, with various levels of simplifications such as 2-D simulations, 3-D simulations with and without chemical reactions, simulations with and without turbulent conditions, etc., have been performed. These simulations have helped determine the model inlet flow characteristics and the important factors that affect the combustor inlet flow and the sensitivity of the flow field to these simplifications. In the proposed paper, CFD modeling of the hypersonic inlet, results from the simulations and comparison with available experimental results will be presented.

Venkatapathy, E.↗

Simulation of Unsteady Combustion in a Ramjet Engine Using a Highly Parallel Computer

Combustion instability in ramjets is a complex phenomenon that involve nonlinear interaction between acoustic waves, vortex motion and unsteady heat release in the combustor. To numerically simulate this 3-D, transient phenomenon, enormous computer resources (time, memory and disk storage) are required. Although current generation vector supercomputers are capable of providing adequate resources for simulations of this nature, their high cost and limited availability, makes such machines less than satisfactory for routine use. The primary focus of this study is to assess the feasibility of using highly parallel computer systems as a cost-effective alternative for conducting such unsteady flow simulations. Towards this end, a large-eddy simulation model for combustion instability was implemented on the Intel iPSC/860 and a careful study was conducted to determine the benefits and the problems associated with the use of such machines for transient simulations. Details of this study along with the results obtained from the unsteady combustion simulations carried out on the iPSC/860 are discussed in this paper.

Menon, Suresh↗

A Comparison of Three-Dimensional Simulations of Traveling-Wave Tube Cold-Test Characteristics Using CST MICROWAVE STUDIO and MAFIA

Previously, it was shown that MAFIA (solutions of Maxwell's equations by the Finite Integration Algorithm), a three-dimensional simulation code, can be used to produce accurate cold-test characteristics including frequency-phase dispersion, interaction impedance, and attenuation for traveling-wave tube (TWT) slow-wave structures. In an effort to improve user-friendliness and simulation time, a model was developed to compute the cold-test parameters using the electromagnetic field simulation software package CST MICROWAVE STUDIO (MWS). Cold-test parameters were calculated for several slow-wave circuits including a ferruled coupled-cavity, a folded waveguide, and a novel finned-ladder circuit using both MWS and MAFIA. Comparisons indicate that MWS provides more accurate cold-test data with significantly reduced simulation times. Both MAFIA and MWS are based on the finite integration (FI) method; however, MWS has several advantages over MAFIA. First, it has a Windows based interface for PC operation, making it very user-friendly, whereas MAFIA is UNIX based. MWS uses a new Perfect Boundary Approximation (PBA), which increases the accuracy of the simulations by avoiding stair step approximations associated with MAFIA's representation of structures. Finally, MWS includes a Visual Basic for Applications (VBA) compatible macro language that enables the simulation process to be automated and allows for the optimization of user-defined goal functions, such as interaction impedance.

Chevalier, C. T.↗

A Distributed Simulation Facility to Support Human Factors Research in Advanced Air Transportation Technology

A distributed real-time simulation of the civil air traffic environment developed to support human factors research in advanced air transportation technology is presented. The distributed environment is based on a custom simulation architecture designed for simplicity and flexibility in human experiments. Standard Internet protocols are used to create the distributed environment, linking all advanced cockpit simulator, all Air Traffic Control simulator, and a pseudo-aircraft control and simulation management station. The pseudo-aircraft control station also functions as a scenario design tool for coordinating human factors experiments. This station incorporates a pseudo-pilot interface designed to reduce workload for human operators piloting multiple aircraft simultaneously in real time. The application of this distributed simulation facility to support a study of the effect of shared information (via air-ground datalink) on pilot/controller shared situation awareness and re-route negotiation is also presented.

Amonlirdviman, Keith↗

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. we consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise vorticity is generated and ejected due to the oscillatory motion of the lid. The present simulations complement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin (1994). The viscous effects are taken into account by modifying the strength Of the particles, whereas fast multipole schemes employing hundreds of thousands of particles allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies.

Mansour, N. N.↗

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. We consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise voracity is generated and ejected due to the oscillatory motion of the lid. The present simulations c Implement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University. Jacobson and Reynolds used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds. More recently, Lachowiez and Wlezien are investigating the flow generated by an electro-mechanically driven lid to be used for assertion control in aerodynamic applications. We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin. The viscous effects are taken into account by modifying the strength of the particles, whereas fast multipole schemes employing hundreds of thousands ol'particle's allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies. Estimates of the forces needed to drive the devices will also be presented.

Mansour, N. N.↗

A Study of the Electrostatic Interaction Between Insulators and Martian/Lunar Soil Simulants

Using our previous experience with the Mars Environmental Compatibility Assessment (MECA) electrometer, we have designed a new type of aerodynamic electrometer. The goal of the research was to measure the buildup of electrostatic surface charge on a stationary cylindrical insulator after windborne granular particles have collided with the insulator surface in a simulated dust storm. The experiments are performed inside a vacuum chamber. This allows the atmospheric composition and pressure to be controlled in order to simulate the atmospheric conditions near the equator on the Martian surface. An impeller fan was used to propel the dust particles at a cylindrically shaped insulator under low vacuum conditions. We tested the new electrometer in a 10 mbar CO2 atmosphere by exposing two types of cylindrical insulators, Teflon (1.9 cm diameter) and Fiberglass (2.5 cm diameter), to a variety of windborne granular particulate materials. The granular materials tested were JSC Mars-1 simulant, which is a mixture of coarse and fine (<5microns diameter) particle sizes, and some of the major mineral constituents of the Martian soil. The minerals included Ottawa sand (SiO2), iron oxide (Fe2O3), aluminum oxide (Al2O3) and magnesium oxide (MgO). We also constructed a MECA-like electrometer that contained an insulator capped planar electrode for measuring the amount of electrostatic charge produced by rubbing an insulator surface over Martian and lunar soil simulants. The results of this study indicate that it is possible to detect triboelectric charging of insulator surfaces by windborne Martian soil simulant, and by individual mineral constituents of the soil simulant. We have also found that Teflon and Fiberglass insulator surfaces respond in different ways by developing opposite polarity surface charge, which decays at different rates after the particle impacts cease.

Mantovani, James G.↗

High-Density Plasma Reactors: Simulations for Design

The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.

Hash, David B.↗

Measured and Simulated Dark J-V Characteristics of a-Si:H Single Junction p-i-n Solar Cells Irradiated with 40 keV Electrons

The effect of 40 keV electron irradiation on a-Si:H p-i-n single-junction solar cells was investigated using measured and simulated dark J-V characteristics. EPRI-AMPS and PC-1D simulators were explored for use in the studies. The EPRI-AMPS simulator was employed and simulator parameters selected to produce agreement with measured J-V characteristics. Three current mechanisms were evident in the measured dark J-V characteristics after electron irradiation, namely, injection, shunting and a term of the form CV(sup m). Using a single discrete defect state level at the center of the band gap, good agreement was achieved between measured and simulated J-V characteristics in the forward-bias voltage region where the dark current density was dominated by injection. The current mechanism of the form CV(sup m) was removed by annealing for two hours at 140 C. Subsequent irradiation restored the CV(sup m) current mechanism and it was removed by a second anneal. Some evidence of the CV(sup m) term is present in device simulations with a higher level of discrete density of states located at the center of the bandgap.

Lord, Kenneth↗