Search NASA⌕ Search

SEARCH · Search NASA

Results for “reaction kinetics and dynamics”

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 91 records · Page 5

Optimization and analysis of large chemical kinetic mechanisms using the solution mapping method - Combustion of methane

A method of systematic optimization, solution mapping, as applied to a large-scale dynamic model is presented. The basis of the technique is parameterization of model responses in terms of model parameters by simple algebraic expressions. These expressions are obtained by computer experiments arranged in a factorial design. The developed parameterized responses are then used in a joint multiparameter multidata-set optimization. A brief review of the mathematical background of the technique is given. The concept of active parameters is discussed. The technique is applied to determine an optimum set of parameters for a methane combustion mechanism. Five independent responses - comprising ignition delay times, pre-ignition methyl radical concentration profiles, and laminar premixed flame velocities - were optimized with respect to thirteen reaction rate parameters. The numerical predictions of the optimized model are compared to those computed with several recent literature mechanisms. The utility of the solution mapping technique in situations where the optimum is not unique is also demonstrated.

Frenklach, Michael↗

Kinetics of DSB rejoining and formation of simple chromosome exchange aberrations

PURPOSE: To investigate the role of kinetics in the processing of DNA double strand breaks (DSB), and the formation of simple chromosome exchange aberrations following X-ray exposures to mammalian cells based on an enzymatic approach. METHODS: Using computer simulations based on a biochemical approach, rate-equations that describe the processing of DSB through the formation of a DNA-enzyme complex were formulated. A second model that allows for competition between two processing pathways was also formulated. The formation of simple exchange aberrations was modelled as misrepair during the recombination of single DSB with undamaged DNA. Non-linear coupled differential equations corresponding to biochemical pathways were solved numerically by fitting to experimental data. RESULTS: When mediated by a DSB repair enzyme complex, the processing of single DSB showed a complex behaviour that gives the appearance of fast and slow components of rejoining. This is due to the time-delay caused by the action time of enzymes in biomolecular reactions. It is shown that the kinetic- and dose-responses of simple chromosome exchange aberrations are well described by a recombination model of DSB interacting with undamaged DNA when aberration formation increases with linear dose-dependence. Competition between two or more recombination processes is shown to lead to the formation of simple exchange aberrations with a dose-dependence similar to that of a linear quadratic model. CONCLUSIONS: Using a minimal number of assumptions, the kinetics and dose response observed experimentally for DSB rejoining and the formation of simple chromosome exchange aberrations are shown to be consistent with kinetic models based on enzymatic reaction approaches. A non-linear dose response for simple exchange aberrations is possible in a model of recombination of DNA containing a DSB with undamaged DNA when two or more pathways compete for DSB repair.

DNA Damage↗

Rarefied gas dynamics; Proceedings of the Eighth International Symposium, Stanford University, Stanford, Calif., July 10-14, 1972

Recent studies in the theory of rarefied gases (kinetic theory, shock structure, molecular forces and interactions, evaporation and condensation, expansions and jets, flow past bodies, gas-solid interactions) are presented along with descriptions of numerous experiments in these areas. Some of the topics covered include: singular limits of solutions of Boltzmann's equations, the Boltzmann and Klimontovich formalisms with reference to turbulent resolution, higher-order kinetic theory for the Lorentz gas, intermolecular potentials from shock structure experiments, cooling of particles in aerosol beams, hypersonic-subsonic transition in spherically symmetric expanding flows, absolute velocity measurements in a rarefied gas flow by an ion time-of-flight technique, the product distribution of gas-solid reactions from absolute rate theory, and relaxation of a gas with internal degrees of freedom. Individual items are announced in this issue.

Karamcheti, K.↗

A nonequilibrium model for a moderate pressure hydrogen microwave discharge plasma

This document describes a simple nonequilibrium energy exchange and chemical reaction model to be used in a computational fluid dynamics calculation for a hydrogen plasma excited by microwaves. The model takes into account the exchange between the electrons and excited states of molecular and atomic hydrogen. Specifically, electron-translation, electron-vibration, translation-vibration, ionization, and dissociation are included. The model assumes three temperatures, translational/rotational, vibrational, and electron, each describing a Boltzmann distribution for its respective energy mode. The energy from the microwave source is coupled to the energy equation via a source term that depends on an effective electric field which must be calculated outside the present model. This electric field must be found by coupling the results of the fluid dynamics and kinetics solution with a solution to Maxwell's equations that includes the effects of the plasma permittivity. The solution to Maxwell's equations is not within the scope of this present paper.

Scott, Carl D.↗

Reactive Infiltration of Silicon Melt Through Microporous Amorphous Carbon Preforms

The kinetics of unidirectional capillary infiltration of silicon melt into microporous carbon preforms have been investigated as a function of the pore morphology and melt temperature. The infiltrated specimens showed alternating bands of dark and bright regions, which corresponded to the unreacted free carbon and free silicon regions, respectively. The decrease in the infiltration front velocity for increasing infiltration distances, is in qualitative agreement with the closed-form solution of capillarity driven fluid flow through constant cross section cylindrical pores. However, drastic changes in the thermal response and infiltration front morphologies were observed for minute differences in the preforms microstructure. This suggests the need for a dynamic percolation model that would account for the exothermic nature of the silicon-carbon chemical reaction and the associated pore closing phenomenon.

Sangsuwan, P.↗

The Supercritical Pile GRB Model: The Prompt to Afterglow Evolution

The 'Supercritical Pile' is a very economical gamma ray burst (GRB) model that provides for the efficient conversion of the energy stored in the protons of a Relativistic Blast Wave (RBW) into radiation and at the same time produces - in the prompt GRB phase, even in the absence of any particle acceleration - a spectral peak at an energy sim 1 MeV. We extend this model to include also the evolution of the RBW Lorentz factor Gamma and thus follow the spectral and temporal features of this model into the early GRB afterglow stage. One of the novel features of the present treatment is the inclusion of the feedback of the GRB produced radiation on the evolution of Gamma with radius. This feedback and the presence of kinematic and dynamic thresholds in the model can be the sources of rich time evolution which we have begun to explore. In particular, one can this way obtain afterglow light curves with steep decays followed by the more conventional flatter afterglow slopes, while at the same time preserving the desirable features of the model, i.e. the well defined relativistic electron source and radiative processes that produce the proper peak in the nu F spectra. Furthermore, the existence of a kinematic threshold in this model provides for a operational distinction of the prompt and afterglow GRB stages; in fact, the afterglow stage sets in when the RBW Lorentz factor cannot anymore fulfill the kinematic condition for pair formation in the photon - proton pair production reactions that constitute the fundamental process for the dissipation of the blast wave kinetic energy. We present the results of a specific set of parameters of this model with emphasis on the multiwavelength prompt emission and transition to the early afterglow.

Kazanas, Demosthenes↗

A Model of Reduced Kinetics for Alkane Oxidation Using Constituents and Species for N-Heptane

The reduction of elementary or skeletal oxidation kinetics to a subgroup of tractable reactions for inclusion in turbulent combustion codes has been the subject of numerous studies. The skeletal mechanism is obtained from the elementary mechanism by removing from it reactions that are considered negligible for the intent of the specific study considered. As of now, there are many chemical reduction methodologies. A methodology for deriving a reduced kinetic mechanism for alkane oxidation is described and applied to n-heptane. The model is based on partitioning the species of the skeletal kinetic mechanism into lights, defined as those having a carbon number smaller than 3, and heavies, which are the complement of the species ensemble. For modeling purposes, the heavy species are mathematically decomposed into constituents, which are similar but not identical to groups in the group additivity theory. From analysis of the LLNL (Lawrence Livermore National Laboratory) skeletal mechanism in conjunction with CHEMKIN II, it is shown that a similarity variable can be formed such that the appropriately non-dimensionalized global constituent molar density exhibits a self-similar behavior over a very wide range of equivalence ratios, initial pressures and initial temperatures that is of interest for predicting n-heptane oxidation. Furthermore, the oxygen and water molar densities are shown to display a quasi-linear behavior with respect to the similarity variable. The light species ensemble is partitioned into quasi-steady and unsteady species. The reduced model is based on concepts consistent with those of Large Eddy Simulation (LES) in which functional forms are used to replace the small scales eliminated through filtering of the governing equations; in LES, these small scales are unimportant as far as the overwhelming part of dynamic energy is concerned. Here, the scales thought unimportant for recovering the thermodynamic energy are removed. The concept is tested by using tabular information from the LLNL skeletal mechanism in conjunction with CHEMKIN II utilized as surrogate ideal functions replacing the necessary functional forms. The test reveals that the similarity concept is indeed justified and that the combustion temperature is well predicted, but that the ignition time is over-predicted, a fact traced to neglecting a detailed description of the processes leading to the heavies chemical decomposition. To palliate this deficiency, functional modeling is incorporated into this conceptual reduction in addition to the modeling the evolution of the global constituent molar density, the enthalpy evolution of the heavies, the contribution to the reaction rate of the unsteady lights from other light species and from the heavies, the molar density evolution of oxygen and water, and the mole fractions of the quasisteady light species. The model is compact in that there are only nine species-related progress variables. Results are presented showing the performance of the model for predicting the temperature and species evolution. The model reproduces the ignition time over a wide range of equivalence ratios, initial pressure, and initial temperature.

Harstad, Kenneth G.↗

Combustion of Metals in Reduced-Gravity and Extraterrestrial Environment

As a result of the ongoing exploration of Mars and the several unmanned and manned missions planned for the future, increased attention has been given to the use of the natural resources of the planet for rocket propellant production and energy generation. Since the atmosphere of Mars consists of approximately 95% carbon dioxide (CO2), this gas is the resource of choice to be employed for these purposes. Unfortunately, CO2 is also a final product in most combustion reactions, requiring further processing to extract useful reactants such as carbon monoxide (CO), oxygen (O2), and hydrocarbons. An exception is the use Of CO2 as an oxidizer reacting directly with metal fuel. Since many metals burn vigorously with CO2, these may be used as an energy source and as propellants for an ascent/descent vehicle in sample-collection missions on Mars. In response to NASA's Human Exploration and Development of Space (HEDS) Enterprise to search for appropriate in-situ resource utilization techniques, this investigation will study the burning characteristics of promising metal/CO2 combinations. The use of reduced gravity is essential to eliminate the intrusive buoyant flows that plague the high-temperature metal reactions, to remove the destructive effect of gravity on the shape of molten metal samples, and to study the influence of radiative heat transfer from solid oxides undisturbed by natural convection. In studies with large metal specimens, the burning process is invariably influenced by strong convective currents that accelerate the reaction and shorten the burning times. Although these currents are nearly absent from small burning particles, the high emissivity of the flames, rapid reaction, small length scales, and intermittent explosions make the gathering of any useful information on burning rates and flame structure very difficult. This investigation has the ultimate goal of providing a careful probing of flame structure and dynamics by taking advantage of large, free-floating spherical metal samples and their corresponding long burning times available in reduced gravity. The first set of experiments has been conducted with magnesium (Mg) samples burning in the low-gravity environment generated by an aircraft flying parabolic trajectories. Owing to its high adiabatic flame temperature, oxidizer/fuel ratio, and heat per unit mass of fuel, as well as its low toxicity and low ignition temperature, Mg has been identified as a promising metal fuel with CO2 as oxidizer. The experimental effort is complemented by the development of a numerical model combining gas-phase chemical kinetics and transport mechanisms.

Abbud-Madrid, A.↗

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗

Mathematical description of complex chemical kinetics and application to CFD modeling codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗

Mathematical Description of Complex Chemical Kinetics and Application to CFD Modeling Codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗

Model of ASTM Flammability Test in Microgravity: Iron Rods

There is extensive qualitative results from burning metallic materials in a NASA/ASTM flammability test system in normal gravity. However, this data was shown to be inconclusive for applications involving oxygen-enriched atmospheres under microgravity conditions by conducting tests using the 2.2-second Lewis Research Center (LeRC) Drop Tower. Data from neither type of test has been reduced to fundamental kinetic and dynamic systems parameters. This paper reports the initial model analysis for burning iron rods under microgravity conditions using data obtained at the LERC tower and modeling the burning system after ignition. Under the conditions of the test the burning mass regresses up the rod to be detached upon deceleration at the end of the drop. The model describes the burning system as a semi-batch, well-mixed reactor with product accumulation only. This model is consistent with the 2.0-second duration of the test. Transient temperature and pressure measurements are made on the chamber volume. The rod solid-liquid interface melting rate is obtained from film records. The model consists of a set of 17 non-linear, first-order differential equations which are solved using MATLAB. This analysis confirms that a first-order rate, in oxygen concentration, is consistent for the iron-oxygen kinetic reaction. An apparent activation energy of 246.8 kJ/mol is consistent for this model.

Steinberg, Theodore A↗

Propellant Chemistry for CFD Applications

Current concepts for reusable launch vehicle design have created renewed interest in the use of RP-1 fuels for high pressure and tri-propellant propulsion systems. Such designs require the use of an analytical technology that accurately accounts for the effects of real fluid properties, combustion of large hydrocarbon fuel modules, and the possibility of soot formation. These effects are inadequately treated in current computational fluid dynamic (CFD) codes used for propulsion system analyses. The objective of this investigation is to provide an accurate analytical description of hydrocarbon combustion thermodynamics and kinetics that is sufficiently computationally efficient to be a practical design tool when used with CFD codes such as the FDNS code. A rigorous description of real fluid properties for RP-1 and its combustion products will be derived from the literature and from experiments conducted in this investigation. Upon the establishment of such a description, the fluid description will be simplified by using the minimum of empiricism necessary to maintain accurate combustion analyses and including such empirical models into an appropriate CFD code. An additional benefit of this approach is that the real fluid properties analysis simplifies the introduction of the effects of droplet sprays into the combustion model. Typical species compositions of RP-1 have been identified, surrogate fuels have been established for analyses, and combustion and sooting reaction kinetics models have been developed. Methods for predicting the necessary real fluid properties have been developed and essential experiments have been designed. Verification studies are in progress, and preliminary results from these studies will be presented. The approach has been determined to be feasible, and upon its completion the required methodology for accurate performance and heat transfer CFD analyses for high pressure, tri-propellant propulsion systems will be available.

Farmer, R. C.↗

Arcjet thruster research and technology. Phase 1: Executive summary

The principal objective of this two phase program is to conduct the development research required to make the low power arcjet a flight ready technology. Many important results were obtained during Phase 1 to move closer to this objective. Fundamental analyses were performed of the arcjet nozzle, the gas kinetic reaction effects, the thermal environment, and the arc stabilizing vortex. These aided the conceptual understanding of the arcjet and guided design work. A hydrazine (N2H4) arcjet was designed that combined a flight qualified catalyst bed with a modular arcjet. Extensive testing was performed which demonstrated the feasibility of using this propellant in an arcjet for the first time. Startup techniques were developed, stability maintained, material compatibility tests conducted, and performance mapping tests performed. Specific impulse values from 400 to 730 seconds were produced with a non-optimized design. These levels are higher than were originally thought possible and proved that extremely high enthalpy values can be obtained with constricted arc technology. Erosion rate data are promising for lifetime extensions to meet flight application requirements. Power control unit (PCU) development was started with the design and fabrication of a laboratory high switching frequency supply. Valuable data were obtained on PCU operation and on the interaction with the dynamic arc. Phase 2 efforts presently underway are resolving key issues for multi-hundred hour lifetimes, are continuing to investigate arcjet/PCU interactions, and will demonstrate duty cycle N2H4 arcjet/PCU operation in a simulated flight mode for lifetimes consistent with initial applications.

Source record↗

Status of Middle Atmosphere-Climate Models: Results SPARC-GRIPS

The middle atmosphere is an important component of the climate system, primarily because of the radiative forcing of ozone. Middle atmospheric ozone can change, over long times, because of changes in the abundance of anthropogenic pollutants which catalytically destroy it, and because of the temperature sensitivity of kinetic reaction rates. There is thus a complex interaction between ozone, involving chemical and climatic mechanisms. One question of interest is how ozone will change over the next decades , as the "greenhouse-gas cooling" of the middle atmosphere increases but the concentrations of chlorine species decreases (because of policy changes). concerns the climate biases in current middle atmosphere-climate models, especially their ability to simulate the correct seasonal cycle at high latitudes, and the existence of temperature biases in the global mean. A major obstacle when addressing this question This paper will present a summary of recent results from the "GCM-Reality Intercomparison Project for SPARC" (GRIPS) initiative. A set of middle atmosphere-climate models has been compared, identifying common biases. Mechanisms for these biases are being studied in some detail, including off-line assessments of the radiation transfer codes and coordinated studies of the impacts of gravity wave drag due to sub-grid-scale processes. ensemble of models will be presented, along with numerical experiments undertaken with one or more models, designed to investigate the mechanisms at work in the atmosphere. The discussion will focus on dynamical and radiative mechanisms in the current climate, but implications for coupled ozone chemistry and the future climate will be assessed.

Pawson, Steven↗

Update on the Status of the FLUKA Monte Carlo Transport Code

The FLUKA Monte Carlo transport code is a well-known simulation tool in High Energy Physics. FLUKA is a dynamic tool in the sense that it is being continually updated and improved by the authors. Here we review the progresses achieved in the last year on the physics models. From the point of view of hadronic physics, most of the effort is still in the field of nucleus--nucleus interactions. The currently available version of FLUKA already includes the internal capability to simulate inelastic nuclear interactions beginning with lab kinetic energies of 100 MeV/A up the the highest accessible energies by means of the DPMJET-II.5 event generator to handle the interactions for greater than 5 GeV/A and rQMD for energies below that. The new developments concern, at high energy, the embedding of the DPMJET-III generator, which represent a major change with respect to the DPMJET-II structure. This will also allow to achieve a better consistency between the nucleus-nucleus section with the original FLUKA model for hadron-nucleus collisions. Work is also in progress to implement a third event generator model based on the Master Boltzmann Equation approach, in order to extend the energy capability from 100 MeV/A down to the threshold for these reactions. In addition to these extended physics capabilities, structural changes to the programs input and scoring capabilities are continually being upgraded. In particular we want to mention the upgrades in the geometry packages, now capable of reaching higher levels of abstraction. Work is also proceeding to provide direct import into ROOT of the FLUKA output files for analysis and to deploy a user-friendly GUI input interface.

Pinsky, L.↗

Unsteady numerical simulations of the stability and dynamics of flames

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years (from Feb. 1992) with emphasis on the work performed since the last microgravity combustion workshop. The primary objective of our research is to develop an understanding of the differences in the structure, stability, dynamics and extinction of flames in earth gravity and in microgravity environments. Numerical simulations, in which the various physical and chemical processes can be independently controlled, can significantly advance our understanding of these differences. Therefore, our approach is to use detailed time-dependent, multi-dimensional, multispecies numerical models to perform carefully designed computational experiments. The basic issues we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein. Some of the basic issues we have addressed recently are (1) the relative importance of wall losses and gravity on the extinguishment of downward-propagating flames; (2) the role of hydrodynamic instabilities in the formation of cellular flames; (3) effects of gravity on burner-stabilized flames, and (4) effects of radiative losses and chemical-kinetics on flames near flammability limits. We have also expanded our efforts to include hydrocarbon flames in addition to hydrogen flames and to perform simulations in support of other on-going efforts in the microgravity combustion sciences program. Modeling hydrocarbon flames typically involves a larger number of species and a much larger number of reactions when compared to hydrogen. In addition, more complex radiation models may also be needed. In order to efficiently compute such complex flames recent developments in parallel computing have been utilized to develop a state-of-the-art parallel flame code. This is discussed below in some detail after a brief discussion of the numerical models.

Kailasanath, K.↗

Chemical Reaction and Flow Modeling in Fullerene and Nanotube Production

The development of processes to produce fullerenes and carbon nanotubes has largely been empirical. Fullerenes were first discovered in the soot produced by laser ablation of graphite [1]and then in the soot of electric arc evaporated carbon. Techniques and conditions for producing larger and larger quantities of fullerenes depended mainly on trial and error empirical variations of these processes, with attempts to scale them up by using larger electrodes and targets and higher power. Various concepts of how fullerenes and carbon nanotubes were formed were put forth, but very little was done based on chemical kinetics of the reactions. This was mainly due to the complex mixture of species and complex nature of conditions in the reactors. Temperatures in the reactors varied from several thousand degrees Kelvin down to near room temperature. There are hundreds of species possible, ranging from atomic carbon to large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to complexes of metals and carbon. Most of the chemical kinetics of the reactions and the thermodynamic properties of clusters and complexes have only been approximated. In addition, flow conditions in the reactors are transient or unsteady, and three dimensional, with steep spatial gradients of temperature and species concentrations. All these factors make computational simulations of reactors very complex and challenging. This article addresses the development of the chemical reaction involved in fullerene production and extends this to production of carbon nanotubes by the laser ablation/oven process and by the electric arc evaporation process. In addition, the high-pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. The first one addresses the thermochemical aspects of modeling; and considers the development of chemical rate equations, estimates of reaction rates, and thermodynamic properties where they are available. The second part addresses modeling of the arc process for fullerene and carbon nanotube production using O-D, 1-D and 2-D fluid flow models. The third part addresses simulations of the pulsed laser ablation process using time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser ablation process. The fourth part addresses steady state modeling in O-D and 2-D of the HiPco process. In each of the simulations, there is a variety of simplifications that are made that enable one to concentrate on one aspect or another of the process. There are simplifications that can be made to the chemical reaction models , e.g. reduction in number of species by lumping some of them together in a representative species. Other simulations are carried out by eliminating the chemistry altogether in order to concentrate on the fluid dynamics. When solving problems with a large number of species in more than one spatial dimension, it is almost imperative that the problem be decoupled by solving for the fluid dynamics to find the fluid motion and temperature history of "particles" of fluid moving through a reactor. Then one can solve the chemical rate equations with complex chemistry following the temperature and pressure history. One difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to take dilution and mixing into account. This changes the ratio of carbon species to background gas. Commercially available codes may have no provision for including dilution as part of the input. One must the write special solvers for including dilution in decoupled problems. The article addresses both ful1erene production and single-walled carbon nanotube (SWNT) production. There are at least two schemes or concepts of SWNT growth. This article will only address growth in the gas phase by carbon and catalyst cluster growth and SW T formation by the addition of carbon. There are other models that conceive of SWNT growth as a phase separation process from clusters me up carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We will not deal with that concept in this article. Further research is needed to determine the rates at which these composite clusters form, evaporate, and segregate.

Scott, Carl D.↗