Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution chemistry”

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 487 records · Page 27

Physical Principle for Generation of Randomness

A physical principle (more precisely, a principle that incorporates mathematical models used in physics) has been conceived as the basis of a method of generating randomness in Monte Carlo simulations. The principle eliminates the need for conventional random-number generators. The Monte Carlo simulation method is among the most powerful computational methods for solving high-dimensional problems in physics, chemistry, economics, and information processing. The Monte Carlo simulation method is especially effective for solving problems in which computational complexity increases exponentially with dimensionality. The main advantage of the Monte Carlo simulation method over other methods is that the demand on computational resources becomes independent of dimensionality. As augmented by the present principle, the Monte Carlo simulation method becomes an even more powerful computational method that is especially useful for solving problems associated with dynamics of fluids, planning, scheduling, and combinatorial optimization. The present principle is based on coupling of dynamical equations with the corresponding Liouville equation. The randomness is generated by non-Lipschitz instability of dynamics triggered and controlled by feedback from the Liouville equation. (In non-Lipschitz dynamics, the derivatives of solutions of the dynamical equations are not required to be bounded.)

Zak, Michail↗

Calibration of Axisymmetric and Quasi-1D Solvers for High Enthalpy Nozzles

The proposed paper will present a numerical investigation of the flow characteristics and boundary layer development in the nozzles of high enthalpy shock tunnel facilities used for hypersonic propulsion testing. The computed flow will be validated against existing experimental data. Pitot pressure data obtained at the entrance of the test cabin will be used to validate the numerical simulations. It is necessary to accurately model the facility nozzles in order to characterize the test article flow conditions. Initially the axisymmetric nozzle flow will be computed using a Navier Stokes solver for a range of reservoir conditions. The calculated solutions will be compared and calibrated against available experimental data from the DLR HEG piston-driven shock tunnel and the 16-inch shock tunnel at NASA Ames Research Center. The Reynolds number is assumed to be high enough at the throat that the boundary layer flow is assumed turbulent at this point downstream. The real gas affects will be examined. In high Mach number facilities the boundary layer is thick. Attempts will be made to correlate the boundary layer displacement thickness. The displacement thickness correlation will be used to calibrate the quasi-1D codes NENZF and LSENS in order to provide fast and efficient tools of characterizing the facility nozzles. The calibrated quasi-1D codes will be implemented to study the effects of chemistry and the flow condition variations at the test section due to small variations in the driver gas conditions.

Papadopoulos, P. E.↗

Monte Carlo transport model comparison with 1A GeV accelerated iron experiment: heavy-ion shielding evaluation of NASA space flight-crew foodstuff

Deep-space manned flight as a reality depends on a viable solution to the radiation problem. Both acute and chronic radiation health threats are known to exist, with solar particle events as an example of the former and galactic cosmic rays (GCR) of the latter. In this experiment Iron ions of 1A GeV are used to simulate GCR and to determine the secondary radiation field created as the GCR-like particles interact with a thick target. A NASA prepared food pantry locker was subjected to the iron beam and the secondary fluence recorded. A modified version of the Monte Carlo heavy ion transport code developed by Zeitlin at LBNL is compared with experimental fluence. The foodstuff is modeled as mixed nuts as defined by the 71st edition of the Chemical Rubber Company (CRC) Handbook of Physics and Chemistry. The results indicate a good agreement between the experimental data and the model. The agreement between model and experiment is determined using a linear fit to ordered pairs of data. The intercept is forced to zero. The slope fit is 0.825 and the R2 value is 0.429 over the resolved fluence region. The removal of an outlier, Z=14, gives values of 0.888 and 0.705 for slope and R2 respectively. c2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

NASA Discipline Radiation Health↗

Importance of formaldehyde in cloud chemistry

A physical-chemical model which is an extension of that of Hong and Carmichael (1983) is used to investigate the role of formaldehyde in cloud chemistry. This model takes into account the mass transfer of SO2, O3, NH3, HNO3, H2O2, CO2, HCl, HCHO, O2, OH and HO2 into cloud droplets and their subsequent chemical reactions. The model is used to assess the importance of S(IV)-HCHO adduct formation, the reduction of H2O2 by HCHO, HCHO-free radical interactions, and the formation of HCOOH in the presence of HCHO in cloud droplets. Illustrative calculations indicate that the presence of HCHO inhibits sulfate production rate in cloud droplets. The direct inhibition of sulfate production rate in cloudwater due to nucleophilic addition of HSO3(-) to HCHO(aq) to form hydroxymethanesulfonate is generally low for concentrations of HCHO typical of ambient air. However, inhibition of sulfate production due to formaldehyde-free radical interactions in solution can be important. These formaldehyde-free radical reactions can also generate appreciable quantities of formic acid.

Adewuyi, Y. G.↗

Biogeochemical Protocols and Diagnostics for the CMIP6 Ocean Model Intercomparison Project (OMIP)

The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948-2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF [subscript] 6) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and 14CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facilitate their implementation.

Carbon↗

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.

Sharma, Maitreyee↗

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.

Sharma, Maitreyee↗

Center for Isotope Cosmochemistry and Geochronology at NASA Johnson Space Center

Ten years in the making the Center for Isotope Cosmochemistry and Geochronology (CICG) at NASA Johnson Space Center is designed to be a cooperative analytical facility where novel techniques are developed and performed. Our team works collaboratively to operate this integrated lab suite comprising sample preparation, metal-free clean chemistry, and mass spectrometers to measure a wide variety of elements and isotopes in astromaterials. Contamination control is enacted throughout by a positive pressure cascade of HEPA filtered air, weekly surface cleaning and particle counting. Astromaterials processing is performed with application-specific mortar and pestles, heavy liquid separation, a Frantz magnetic separator, as well as hand-picking using optical microscopes in laminar flow boxes. Microsampling is also carried out using a New Wave MicroMill. Elements for isotopic analysis are isolated from acid-digested samples using ion-exchange chromatography. Cross-contamination is minimized with dedicated laminar flow boxes and labware for different purification protocols. Instrumentation includes a ThermoFisher Element XR inductively coupled plasma mass spectrometer (ICP-MS) for quantification of elemental concentrations in dissolved sample solutions. For in situ analysis, a Teledyne Photon Machines Analyte Excite+ laser ablation system is coupled to the Element XR. The ThermoFisher Triton thermal ionization mass spectrometer (TIMS) is used to measure the isotopic composition of elements such as CA, K, Cr, Rb, Sr, Sm, Nd at high precision in processed samples. For a more diverse range of isotopic analyses, CICG houses a Nu Sapphire 1700 multi-collector inductively coupled plasma mass spectrometer (MC-ICP-MS), which is a large geometry instrument with a multi-pole collision cell for interference removal. An Applied Spectra Instruments J200 femto-second laser is planned for install this Fall and will be coupled to the Nu Sapphire 1700. NASA’s Internal Scientist Funding Model (ISFM)supports this facility. Advanced analytical geochemistry techniques are applied to targeted astromaterials to address important problems in planetary sciences, with a longer-term goal of preparing for samples collected and returned by space missions. CICG is open to collaborations on projects that utilize the facility’s unique capabilities and are within its research scope. We are particularly eager to support and share expertise with the next generation of planetary scientists.

Justin I Simon↗

Advancements in Dual-Pump Broadband CARS for Supersonic Combustion Measurements

Space- and time-resolved measurements of temperature and species mole fractions of nitrogen, oxygen, and hydrogen were obtained with a dual-pump coherent anti-Stokes Raman spectroscopy (CARS) system in hydrogen-fueled supersonic combustion free jet flows. These measurements were taken to provide time-resolved fluid properties of turbulent supersonic combustion for use in the creation and verification of computational fluid dynamic (CFD) models. CFD models of turbulent supersonic combustion flow currently facilitate the design of air-breathing supersonic combustion ramjet (scramjet) engines. Measurements were made in supersonic axi-symmetric free jets of two scales. First, the measurement system was tested in a laboratory environment using a laboratory-scale burner (approx.10 mm at nozzle exit). The flow structures of the laboratory-burner were too small to be resolved with the CARS measurements volume, but the composition and temperature of the jet allowed the performance of the system to be evaluated. Subsequently, the system was tested in a burner that was approximately 6 times larger, whose length scales are better resolved by the CARS measurement volume. During both these measurements, weaknesses of the CARS system, such as sensitivity to vibrations and beam steering and inability to measure temperature or species concentrations in hydrogen fuel injection regions were indentified. Solutions were then implemented in improved CARS systems. One of these improved systems is a dual-pump broadband CARS technique called, Width Increased Dual-pump Enhanced CARS (WIDECARS). The two lowest rotational energy levels of hydrogen detectable by WIDECARS are H2 S(3) and H2 S(4). The detection of these lines gives the system the capability to measure temperature and species concentrations in regions of the flow containing pure hydrogen fuel at room temperature. WIDECARS is also designed for measurements of all the major species (except water) in supersonic combustion flows fueled with hydrogen and hydrogen/ethylene mixtures (N2, O2, H2, C2H4, CO, and CO2). This instrument can characterize supersonic combustion fueled with surrogate fuel mixtures of hydrogen and ethylene. This information can lead to a better understanding of the chemistry and performance of supersonic combustion fueled with cracked jet propulsion (JP)-type fuel.

Tedder, Sarah Augusta Umberger↗

Analysis of Fuel Vaporization, Fuel-Air Mixing, and Combustion in Integrated Mixer-Flame Holders

Requirements to limit pollutant emissions from the gas turbine engines for the future High-Speed Civil Transport (HSCT) have led to consideration of various low-emission combustor concepts. One such concept is the Integrated Mixer-Flame Holder (IMFH). This report describes a series of IMFH analyses performed with KIVA-II, a multi-dimensional CFD code for problems involving sprays, turbulence, and combustion. To meet the needs of this study, KIVA-II's boundary condition and chemistry treatments are modified. The study itself examines the relationships between fuel vaporization, fuel-air mixing, and combustion. Parameters being considered include: mixer tube diameter, mixer tube length, mixer tube geometry (converging-diverging versus straight walls), air inlet velocity, air inlet swirl angle, secondary air injection (dilution holes), fuel injection velocity, fuel injection angle, number of fuel injection ports, fuel spray cone angle, and fuel droplet size. Cases are run with and without combustion to examine the variations in fuel-air mixing and potential for flashback due to the above parameters. The degree of fuel-air mixing is judged by comparing average, minimum, and maximum fuel/air ratios at the exit of the mixer tube, while flame stability is monitored by following the location of the flame front as the solution progresses from ignition to steady state. Results indicate that fuel-air mixing can be enhanced by a variety of means, the best being a combination of air inlet swirl and a converging-diverging mixer tube geometry. With the IMFH configuration utilized in the present study, flashback becomes more common as the mixer tube diameter is increased and is instigated by disturbances associated with the dilution hole flow.

Deur, J. M.↗

Geochemical Properties of Rocks and Soils in Gusev Crater, Mars: APXS Results from Cumberland Ridge to Home Plate

The Mars Exploration Rover Spirit landed in Gusev crater on Jan. 4, 2004. Spirit has traversed the Gusev crater plains, ascended to the top of Husband Hill, and entered into the Inner Basin of the Columbia Hills. The Athena science payload onboard Spirit has recorded numerous measurements on the chemistry and mineralogy of materials encountered during nearly 2 Mars years of operation within the crater. Rocks and soils have been grouped into classes based upon their unique differences in mineralogy and chemistry [1-3]. Some of the most significant chemical discoveries include the composition of Adirondack class flood basalts [4-6]; high sulfur in Clovis and Peace Class rocks [7,2]; high P and Ti in Wishstone Class rocks [7,2]; composition of alkalic basalts [2,6]; very high S in Paso Robles class soils [7,2], and the possible occurrence of a smectite-like chemical composition in Independence class rocks [8]. Water has played a significant role in the alteration of rocks and soils in the Columbia Hills. The occurrence of goethite and ferric sulfate alone suggests that liquid water was involved in their formation [3]. The pervasively altered materials in Husband Hill outcrops and rocks may have formed by the aqueous alteration of basaltic rocks, volcaniclastic materials, and/or impact ejecta by solutions that were rich in acid-volatile elements [2]. The objective of this paper is to provide an update on the health of the Alpha Particle X-ray Spectrometer (APXS) and to expand the geochemical dataset from sol 470 to sol 1368. Specific objectives are to (1) update the rock and soil classifications, (2) characterize elemental relationships among the major rock and soil classes, and (3) evaluate the involvement of water in the formation or alteration of the materials in these classes.

Ming, D. W.↗

Effect of Surface Nonequilibrium Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using finite-rate gas/surface interaction model provides time-accurate solutions for multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal Response and Ablation Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas momentum conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of gas/surface interaction chemistry between air and carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was a Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kang↗

Non-Magnetic, Tough, Corrosion- and Wear-Resistant Knives From Bulk Metallic Glasses and Composites

Quality knives are typically fabricated from high-strength steel alloys. Depending on the application, there are different requirements for mechanical and physical properties that cause problems for steel alloys. For example, diver's knives are generally used in salt water, which causes rust in steel knives. Titanium diver's knives are a popular alternative due to their salt water corrosion resistance, but are too soft to maintain a sharp cutting edge. Steel knives are also magnetic, which is undesirable for military applications where the knives are used as a tactical tool for diffusing magnetic mines. Steel is also significantly denser than titanium (8 g/cu cm vs. 4.5 g/cu cm), which results in heavier knives for the same size. Steel is hard and wear-resistant, compared with titanium, and can keep a sharp edge during service. A major drawback of both steel and titanium knives is that they must be ground or machined into the final knife shape from a billet. Since most knives have a mirrored surface and a complex shape, manufacturing them is complex. It would be more desirable if the knife could be cast into a net or near-net shape in a single step. The solution to the deficiencies of titanium, steel, and ceramic knives is to fabricate them using bulk metallic glasses (or composites). These alloys can be cast into net or near-net shaped knives with a combination of properties that exceed both titanium and steel. A commercially viable BMG (bulk metallic glass) or composite knife is one that exhibits one or all of the following properties: It is based on titanium, has a self-sharpening edge, can retain an edge during service, is hard, is non-magnetic, is corrosion-resistant against a variety of corrosive environments, is tough (to allow for prying), can be cast into a net-shape with a mirror finish and a complex shape, has excellent wear resistance, and is low-density. These properties can be achieved in BMG and composites through alloy chemistry and processing. For each desired property for knife fabrication and performance, there is an alloy development strategy that optimizes behavior. Although BMG knives have been demonstrated as far back as 1995, they never found commercial success because they had to be ground (which presented problems because the alloys contained beryllium), they weren't low cost (because they weren't cast to a net-shape), they were brittle (because they were made with a low-quality commercial material), and they had extremely poor corrosion resistance (because corrosion was not well-understood in these materials). Ultimately, these shortcomings prevented the widespread commercialization. In the current work, the inventors have applied more than a decade of research on BMGs from Caltech and JPL to develop a better understanding of how to make BMG knives that exhibit an optimal combination of properties, processing and cost. Alloys have been developed based in titanium (and other metals), that exhibit high toughness, high hardness, excellent corrosion resistance, no ferromagnetism, edge-retaining selfsharpening, and the ability to be cast like a plastic using commercially available casting techniques (currently used by commercial companies such as Liquidmetal Technologies and Visser Precision Casting). The inventors argue that depending on the application (diving, military, tactical, utility, etc.) there is an optimal combination of design and alloy composition. Moreover, with new casting technologies not available at the inception of these materials, net-shaped knives can be cast into complex shapes that require no aftermarket forming, except for sharpening using water-cooled polishing wheel. These combinations of discoveries seek to make low-cost BMG knives commercially viable products that have no equal among metal or ceramic knives. Current work at JPL focuses on net-shape casting of these alloys and testing their mechanical properties versus commercially available knives to demonstrate their benefits.

Hoffman, Douglas C.↗

Renal Stone Risk During Space Flight

Space flight produces a number of metabolic and physiological changes in the crewmembers exposed to microgravity. Following launch, body fluid volumes, electrolyte levels, and bone and muscle undergo changes as the human body adapts to the weightless environment. Changes in the urinary chemical composition may lead to the potentially serious consequences of renal stone formation. Previous data collected immediately after space flight indicate changes in the urine chemistry favoring an increased risk of calcium oxalate and uric acid stone formation (n = 323). During short term Shuttle space flights, the changes observed include increased urinary calcium and decreased urine volume, pH and citrate resulting in a greater risk for calcium oxalate and brushite stone formation (n = 6). Results from long duration Shuttle/Mir missions (n = 9) followed a similar trend and demonstrated decreased fluid intake and urine volume and increased urinary calcium resulting in a urinary environment saturated with the calcium stone-forming salts. The increased risk occurs rapidly upon exposure to microgravity, continues throughout the space flight and following landing. Dietary factors, especially fluid intake, or pharmacologic intervention can significantly influence the urinary chemical composition. Increasing fluid intake to produce a daily urine output of 2 liters/day may allow the excess salts in the urine to remain in solution, crystals formation will not occur and a renal stone will not develop. Results from long duration crewmembers (n = 2) who had urine volumes greater than 2.5 L/day minimized their risk of renal stone formation. Also, comparisons of stone-forming risk in short duration crewmembers clearly identified greater risk in those who produced less than 2 liters of urine/day. However, hydration and increased urine output does not correct the underlying calcium excretion due to bone loss and only treats the symptoms and not the cause of the increased urinary salts. Dietary modification and promising pharmacologic treatments may also be used to reduce the potential risk for renal stone formation. Potassium citrate is being used clinically to increase the urinary inhibitor levels to minimize the development of crystals and the growth of renal stones. Bisphosphonates are a class of drugs recently shown to help in patients with osteoporosis by inhibiting the loss of bones in elderly patients. This drug could potentially prevent the bone loss observed in astronauts and thereby minimize the increase in urinary calcium and reduce the risk for renal stone development. Results of NASA's renal stone risk assessment program clearly indicate that exposure to microgravity changes the urinary chemical environment such that there is an increased risk for supersaturation of stone-forming salts, including calcium oxalaie and brushite. These studies have indicated specific avenues for development of countermeasures for the increased renal stone risk observed during and following space flight. Increased hydration and implementation of pharmacologic countermeasures should largely mitigate the in-flight risk of renal stones.

Whitson, Peggy A.↗

Chemical application of diffusion quantum Monte Carlo

The diffusion quantum Monte Carlo (QMC) method gives a stochastic solution to the Schroedinger equation. This approach is receiving increasing attention in chemical applications as a result of its high accuracy. However, reducing statistical uncertainty remains a priority because chemical effects are often obtained as small differences of large numbers. As an example, the single-triplet splitting of the energy of the methylene molecule CH sub 2 is given. The QMC algorithm was implemented on the CYBER 205, first as a direct transcription of the algorithm running on the VAX 11/780, and second by explicitly writing vector code for all loops longer than a crossover length C. The speed of the codes relative to one another as a function of C, and relative to the VAX, are discussed. The computational time dependence obtained versus the number of basis functions is discussed and this is compared with that obtained from traditional quantum chemistry codes and that obtained from traditional computer architectures.

Reynolds, P. J.↗

Effect of Non-Equilibrium Surface Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using non-equilibrium gas/surface interaction model provides time-accurate solutions for the multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and AblatioN Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas mass conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between the material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of nonequilibrium gas/surface interaction chemistry between air and the carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kanq↗

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.↗

Triaxial Burke-Schumann Flames with Applications to Flame Synthesis

The problem of a flame generated by three coaxial flows is solved by extending the Burke-Schumann methodology to include a third stream. The solution is particularly relevant to flame synthesis wherein multiple tubes are often employed either to introduce inert as a diffusion barrier or to introduce more than two reactants. The general problem is solved where the inner and outer tubes contain reactants and the middle tube contains either an inert or a third reactant. Relevant examples are considered and the results show that the triaxial Burke-Schumann flame can be substantially more complicated than the traditional Burke-Schumann flame. When the middle flow is inert the flame temperature is no longer constant but increases axially, reaching a maximum at the flame centerline. At the exit the flame does not sit on the tube exit but instead resides between the inner and outer tubes, resulting in an effective barrier for particle build-up on the burner rim. For the case of a third reactant in the middle flow, synthesis chemistry where the inner reaction is endothermic and the outer reaction is exothermic is considered. In addition to showing the flame temperature and flame shape, the results identify conditions wherein reaction is not possible due to insufficient heat transfer from the outer flame to support the inner flame reaction.

Chao, B. H.↗