Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “combustion kinetics”

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 361 records · Page 20

International Symposium on Combustion, 14th, Pennsylvania State University, University Park, Pa., August 20-25, 1972, Proceedings.

Some elementary steps for the formation and destruction of nitric oxide in the H-N-O system, high-temperature oxidation of hydrogen by nitrous oxide in shock waves, and high-temperature fast-flow reactor studies of metal-atom oxidation kinetics are among the topics covered in papers concerned with elementary reactions in combustion. Other areas covered include oscillatory, turbulent, and supersonic combustion; flames in furnaces and combustors; and pollutant formation and destruction in flames. Individual items are announced in this issue.

Source record↗

Extinction of a Free Methanol Droplet in Microgravity

The combustion of isolated, unsupported methanol droplets under microgravity conditions has been studied in the NASA-Lewis 2.2-sec drop tower. The substitution of O2/He mixtures for O2/N2 mixtures results in a substantial increase in the oxygen index required for ignition, the rate of gasification, and the diameter at extinction. As a result, droplets with initial size larger than 1 mm can be produced, deployed, ignited, and burned to extinction, all under microgravity conditions in this facility. A comprehensive time dependent numerical model incorporating detailed multicomponent molecular transport mechanisms and complex elementary chemical kinetics was developed and applied to simulate the entire transient droplet combustion processes (ignition/quasi-steady burning/extinction). Model outputs show good agreement with the experimental measurements, without parameter adjustment. Analysis of modeling results shows that the thermal conductivity of the gas phase under combustion environments increased by a factor of two when the 50-percent O2/50 percent He mixture is used instead of air. The increased thermal conductivity accelerates the heat transfer from the flame front to the droplet surface, resulting in a factor of two increase in the burning rate.

Cho, S. Y.↗

The PDF method for turbulent combustion

Probability Density Function (PDF) methods provide a means of calculating the properties of turbulent reacting flows. They have been successfully applied to many turbulent flames, including some with finite rate kinetic effects. Here the methods are reviewed with an emphasis on computational issues and their application to turbulent combustion.

Pope, S. B.↗

Exploratory Investigation of Coal in Nonequilibrium Plasma

Coal is an abundant natural resource and there is motivation to find new uses for it that do not intrinsically involve combustion. One approach is to explore new ways of processing coal, and in this work, we focus on the transformation of coal in a nonequilibrium plasma generated from an equimolar mixture of nitrogen and hydrogen. The outcome of the nonequilibrium plasma reaction is fundamentally different than a thermal control reaction carried out using the same gas composition, pressure, and temperature range. The nonequilibrium plasma produces a gas mixture that is enriched in acetylene and its derivatives. Furthermore, when compared to the thermal control experiment, the solid char byproduct of the nonequilibrium plasma has a very reactive surface and is spontaneously combustible at ambient temperature. Experiments performed to characterize the reaction kinetics of coal in the plasma suggest that the mechanism proceeds through a sequential process by which the coal particle temperature rises to a point where devolatilization can occur, the devolatilization reaction happens, followed by parallel reactions of released organic vapors in the plasma phase and surface activation. In conclusion, the reaction rate appears to be limited by the time it takes for the coal particle temperature to rise, consistent with previous results reported for reactions of coal in thermal plasma.

Acetylene↗

Understanding how chemical structure affects ignition-delay-time $\phi$-sensitivity

$\phi$-sensitivity is the change in ignition delay time (IDT) with respect to the fuel-to-air equivalence ratio ($\phi$). High $\phi$-sensitivity is a desirable fuel property for applications in advanced compression ignition and multi-mode engine designs. Understanding how $\phi$-sensitivity depends on chemical structure is essential for selecting promising biofuels from the ever-growing list of proposed candidates. Here, we investigate the effect of chemical structure on $\phi$-sensitivity with experiment, simulation, and theory. Experimental Advanced Fuel Ignition Delay Analyzer (AFIDA) measurements for 2,4-dimethylpentane and diisopropyl ether provide evidence that branching and functional groups strongly impact $\phi$-sensitivity. Further insights into this dependence are obtained with 0-D kinetic simulations with existing mechanisms for n-pentane, diethyl ether, 3-pentanone, n-heptane, 2-methylhexane, 2,4-dimethylpentane, and 2,2,3-trimethylbutane. Quantum mechanical (QM) G4 calculations of low-temperature reactions help explain the observed experimental and simulation trends. Specifically, these QM calculations provide theoretical estimates of the ketohydroperoxide (KHP) dissociation rates, the HO 2 formation rates from peroxy radical (ROO), and the “cross-over” temperatures, i.e., the temperature at which ROO dissociation is favored compared to hydroperoxyl radical (QOOH) formation. Each of these reaction rates is compared to the n-alkane reference point to determine the impact of branching and different functional groups. Although kinetic mechanisms typically assume that KHP dissociation rates are invariant of chemical environment, our QM results suggest that this rate can span a range of roughly two orders of magnitude. We also discuss the importance of including the peroxy-hydroperoxy (OO-OOH) hydrogen transfer reaction for branched ethers. Finally, the insights gained assist in proposing a highly $\phi$-sensitive compound, namely, isopropyl propyl ether.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗

Pulse Detonation Engine for Advanced Oxy-Combustion of Coal-Based Fuel for Direct Power Extraction Applications

Pressure gain combustion and magnetohydrodynamic (MHD) systems have the potential to provide a step increase in the efficiency of combined-cycle power plants. Specific advantages include a net pressure increase to the system instead of a pressure drop, the high temperature of the detonation waves can increase the efficiency of power extraction compared to other processes, significant thermal energy can be released in a compact region, and the high velocities of the flow increase extraction of electrical power. In summary, a pressure gain combustor coupled with a MHD has the potential to be transformative. Despite the potential advantages, relatively little research has been conducted considering coupled pressure gain combustion systems with MHD systems. With this background and motivation, the overall goal of this effort was to advance the knowledge, technology, and computational tools associated with coupled detonation and MHD systems. A joint experimental and computational approach was used while seeking to accomplish the goals of this work. Specifically, two pulse detonation engines were developed and used for the experiments to produce detonations. Detonation speeds were measured for a variety of flow and fuel conditions (e.g., methane, propane, with coal particles). Preliminary electrical conductivity measurements were collected. An extensive amount of research was performed to identify sensitivities of detonation behaviors to the presence of combustion products. Computationally, a twofold approach was used in this work. First, a solver was developed for solving the governing equations for a reactive flow with coupled detonation and ionization chemistry. The solver was applied to study the impacts of seed material ionization on detonation. Second, a conservation element-solution element (CE-SE) based numerical solver for detonation studies with a reduced reaction mechanism for oxy-methane combustion was developed and verified on standard test cases. Key findings and contributions from this work are as follows. A system was developed for injecting powderized coal, or other seeding material, into a pulse-detonation system. The influence of a combustion product (i.e., CO 2 ) on detonation behavior was identified. Knowledge gained from this work is applicable to devices such as rotating detonation engines, where combustion products mix with fresh reactants. A system for measuring the electrical conductivity of the exhaust from a pulse-detonation engine was developed. The open-source solver, Clawpack, was extended to solve the reactive Euler equations for simulating detonations. A coupled combustion and ionization chemistry was developed in a single chemical kinetic model for methane oxidation. This model can be used to solve coupled MHD and detonation simulations. It was found that parasitic interactions from ionization chemistry with the magnetic field can reduce the detonation velocity by up to 8%, with a potential impact on power extraction of 15%. It is recommended that interactions between the detonation front and MHD field be considered. Higher gas temperatures and velocities were achieved owing to oxy-fuel detonations. Use of radical dissociation reactions in the reduced reaction mechanism, was found to be critical in predicting detonation temperature and velocity accurately.

01 COAL, LIGNITE, AND PEAT↗

Review of HxPyOz-Catalyzed H + OH Recombination in Scramjet Nozzle Expansions; and Possible Phosphoric Acid Enhancement of Scramjet Flameholding, from Extinction of H3PO4 + H2 - Air Counterflow Diffusion Flames

Recent detailed articles by Twarowski indicate that small quantities of phosphorus oxides and acids in the fuel-rich combustion products of H2 + phosphine (PH3) + air should significantly catalyze H, OH and O recombination kinetics during high-speed nozzle expansions -- to reform H2O, release heat, and approach equilibrium more rapidly and closely than uncatalyzed kinetics. This paper is an initial feasibility study to determine (a) if addition of phosphoric acid vapor (H3PO4) to a H2 fuel jet -- which is much safer than using PH3 -- will allow combustion in a high-speed scramjet engine test without adverse effects on localized flameholding, and (b) if phosphorus-containing exhaust emissions are environmentally acceptable. A well-characterized axisymmetric straight-tube opposed jet burner (OJB) tool is used to evaluate H3PO4 addition effects on the air velocity extinction limit (flame strength) of a H2 versus air counterflow diffusion flame. Addition of nitric oxide (NO), also believed to promote catalytic H-atom recombination, was evaluated for comparison. Two to five mass percent H3PO4 in the H2 jet increased flame strength 4.2%, whereas airside addition decreased it 1%. Adding 5% NO to the H2 caused a 2% decrease. Products of H-atom attack on H3PO4 produced an intense green chemiluminescence near the stagnation point. The resultant exothermic production of phosphorus oxides and acids, with accelerated H-atom recombination, released sufficient heat near the stagnation point to increase flame strength. In conclusion, the addition of H3PO4 vapor (or more reactive P sources) to hydrogen in scramjet engine tests may positively affect flameholding stability in the combustor and thrust production during supersonic expansion -- a possible dual benefit with system design / performance implications. Finally, a preliminary assessment of possible environmental effects indicates that scramjet exhaust emissions should consist of phosphoric acid aerosol, with gradual conversion to phosphate aerosol. This is compared to various natural abundances and sources.

Pellett, Gerald↗

Experimental and Modeling Studies of the Combustion Characteristics of Conventional and Alternative Jet Fuels. Final Report

The objectives of this project have been to develop a comprehensive set of fundamental data regarding the combustion behavior of jet fuels and appropriately associated model fuels. Based on the fundamental study results, an auxiliary objective was to identify differentiating characteristics of molecular fuel components that can be used to explain different fuel behavior and that may ultimately be used in the planning and design of optimal fuel-production processes. The fuels studied in this project were Fischer-Tropsch (F-T) fuels and biomass-derived jet fuels that meet certain specifications of currently used jet propulsion applications. Prior to this project, there were no systematic experimental flame data available for such fuels. One of the key goals has been to generate such data, and to use this data in developing and verifying effective kinetic models. The models have then been reduced through automated means to enable multidimensional simulation of the combustion characteristics of such fuels in real combustors. Such reliable kinetic models, validated against fundamental data derived from laminar flames using idealized flow models, are key to the development and design of optimal combustors and fuels. The models provide direct information about the relative contribution of different molecular constituents to the fuel performance and can be used to assess both combustion and emissions characteristics.

Meeks, Ellen↗

Photo-sensitized ignition of hydrogen/oxygen mixtures for hypersonic flight vehicles

An approach toward improving supersonic ignition and combustion efficiency was investigated under a feasibility program with NASA-Ames. The idea is to selectively irradiate targetted species in the combustion zone. The ensuing photodissociation reactions produce highly reactive radicals which modify the gas phase kinetics. The overall objective of this research was to demonstrate photosensitized ignition of H2/O2 mixtures, and to gain a preliminary understanding of the ignition mechanism, especially as it pertains to the concept of supersonic combustion. To this end, a series of experiments were performed with a model sensitizer, NO2. At subatmospheric pressures and moderate temperatures (above 100 F), mixtures of H2/O2/NO2 are ignitable by a pulsed eximer laser at 308 nm, 115 to 200 mJ. Ignitability was found to be a function of NO2 pressure, temperature, laser energy per pulse, and laser lens focal length (fluence at focal volume). Mixtures could not be ignited in the absence of the sensitizer. While increasing temperature improved sensitized ignition in the lower temperature range, a retarding effect was observed at the highest temperatures. To more thoroughly understand these thermal effects, the UV-Vis and IR spectra of NO2 and H2/O2/NO2 mixtures were studied over the temperature range 70 to 900 F. Initially, increasing temperature shifts the N2O4 in a reversible reaction to 2 NO2 equilibrium toward NO2, favoring photochemical ignition. However, at still higher temperature, NO2 converts to non-absorbing NO + 1/2 O2, thus preventing photosensitized ignition.

Lavid, Moshe↗

Parallel Simulation of Unsteady Turbulent Flames

Time-accurate simulation of turbulent flames in high Reynolds number flows is a challenging task since both fluid dynamics and combustion must be modeled accurately. To numerically simulate this phenomenon, very large computer resources (both time and memory) are required. Although current vector supercomputers are capable of providing adequate resources for simulations of this nature, the high cost and their limited availability, makes practical use of such machines less than satisfactory. At the same time, the explicit time integration algorithms used in unsteady flow simulations often possess a very high degree of parallelism, making them very amenable to efficient implementation on large-scale parallel computers. Under these circumstances, distributed memory parallel computers offer an excellent near-term solution for greatly increased computational speed and memory, at a cost that may render the unsteady simulations of the type discussed above more feasible and affordable.This paper discusses the study of unsteady turbulent flames using a simulation algorithm that is capable of retaining high parallel efficiency on distributed memory parallel architectures. Numerical studies are carried out using large-eddy simulation (LES). In LES, the scales larger than the grid are computed using a time- and space-accurate scheme, while the unresolved small scales are modeled using eddy viscosity based subgrid models. This is acceptable for the moment/energy closure since the small scales primarily provide a dissipative mechanism for the energy transferred from the large scales. However, for combustion to occur, the species must first undergo mixing at the small scales and then come into molecular contact. Therefore, global models cannot be used. Recently, a new model for turbulent combustion was developed, in which the combustion is modeled, within the subgrid (small-scales) using a methodology that simulates the mixing and the molecular transport and the chemical kinetics within each LES grid cell. Finite-rate kinetics can be included without any closure and this approach actually provides a means to predict the turbulent rates and the turbulent flame speed. The subgrid combustion model requires resolution of the local time scales associated with small-scale mixing, molecular diffusion and chemical kinetics and, therefore, within each grid cell, a significant amount of computations must be carried out before the large-scale (LES resolved) effects are incorporated. Therefore, this approach is uniquely suited for parallel processing and has been implemented on various systems such as: Intel Paragon, IBM SP-2, Cray T3D and SGI Power Challenge (PC) using the system independent Message Passing Interface (MPI) compiler. In this paper, timing data on these machines is reported along with some characteristic results.

Menon, Suresh↗

Autoignition study of iso -cetane/tetralin blends at low temperature

iso-Cetane and tetralin are the two fuel components commonly considered in literature diesel surrogates, and are hydrocarbon classes representative of iso-alkanes and naphthoaromatics, respectively. Since both surrogate components are involved in the key reactions controlling the autoignition reactivity of diesel surrogates, autoignition studies of iso-cetane, tetralin, and their blending behavior at low temperatures can be helpful in developing/refining chemical kinetic models of diesel surrogates, in order to better predict diesel ignition response under low temperature combustion conditions. In this investigation, autoignition experiments of iso-cetane, tetralin, and their binary mixtures at varying blending ratios have been conducted in a rapid compression machine (RCM) in the temperature/pressure range of 630–930 K and 10–20 bar. For the neat fuel components, the newly-acquired RCM data have been compared with the literature experimental results at some overlapping conditions. In addition, a chemical kinetic model describing the oxidation of iso-cetane and tetralin compiled/updated from the diesel surrogate model developed by Lawrence Livermore National Laboratory is used to simulate the experimental datasets. It is shown that this model exhibits good agreement with the experimental results of tetralin and also predicts the total ignition delay times of iso-cetane at lower temperatures reasonably well. However, the current model fails to fully capture the increase of total ignition delay time with increasing amount of tetralin in the binary blends observed in the present RCM experiments. Therefore, chemical kinetic analyses of the current model are conducted to help identify possible reasons leading to the discrepancies between simulations and experiments. Furthermore, the chemical kinetic interactions between iso-cetane and tetralin are explored by including additional cross-reactions among iso-cetane, tetralin, and their fuel fragments in the current model to assess their effects on model predictions. As a result, based on the chemical kinetic analysis results, future direction for model refinements is discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Fuel Properties and Kinetics

The Co-Optimization of Fuels and Engines (CoOptima) program is focused on a two-pronged approach to develop new, lower net-carbon fuels and advanced engine technologies in tandem. This effort brings nine national labs and several university and industry partners together to simultaneously identify new fuel candidates and to further engine technology development in the light-, medium-, and heavy-duty vehicle sectors with a focus on fuel economy and emission reduction. The fuel properties and kinetics project is aimed at achieving a better understanding of fuel property kinetics and how they impact engine performance such that more efficient and lower emitting engines can be made a reality. Several promising biofuel candidates were identified through tiered screening approaches and were evaluated utilizing unique bench-scale measurement techniques. Kinetic models were developed to describe autoignition and soot pre-cursor formation mechanisms to identify the most promising candidates that advanced various combustion strategies.

ADVANCED PROPULSION SYSTEMS↗

Tuning Process for the Modified Magnussen Combustion Model

In the application of CFD to turbulent reacting flows, one of the main limitations to predictive accuracy is the chemistry model. Using a full or skeletal kinetics model may provide good predictive ability, however, at considerable computational cost. Adding the ability to account for the interaction between turbulence and chemistry improves the overall fidelity of a simulation but adds to this cost. An alternative is the use of simple models, such as the Magnussen model, which has negligible computational overhead, but lacks general predictive ability except for cases that can be tuned to the flow being solved. In this paper, a technique will be described that allows the tuning of the Magnussen model for an arbitrary fuel and flow geometry without the need to have experimental data for a particular case. The tuning is based on comparing the results of the Magnussen model and full finite-rate chemistry when applied to perfectly and partially stirred reactor sim- ulations. In addition, a modification to the Magnussen model is proposed that allows the upper kinetic limit for the reaction rate to be set, giving better physical agreement with full kinetic mechanisms. In order to improve the agreement with flame temperatures, the thermal properties of the product species is adjusted to better match the mixture proper- ties of the full mechanism. The combustion model is then applied to the simulation of a representative scramjet flowpath, and the results compared to experimental data and other kinetic models. This procedure allows a simple reacting model to be used in a predictive manner, and affords significant savings in computational costs for CFD simulations.

Norris, A. T.↗

Application of mixing-controlled combustion models to gas turbine combustors

Gas emissions were studied from a staged Rich Burn/Quick-Quench Mix/Lean Burn combustor were studied under test conditions encountered in High Speed Research engines. The combustor was modeled at conditions corresponding to different engine power settings, and the effect of primary dilution airflow split on emissions, flow field, flame size and shape, and combustion intensity, as well as mixing, was investigated. A mathematical model was developed from a two-equation model of turbulence, a quasi-global kinetics mechanism for the oxidation of propane, and the Zeldovich mechanism for nitric oxide formation. A mixing-controlled combustion model was used to account for turbulent mixing effects on the chemical reaction rate. This model assumes that the chemical reaction rate is much faster than the turbulent mixing rate.

Nguyen, Hung Lee↗

A Theoretical Study of NH 2 Radical Reactions with Propane and Its Kinetic Implications in NH 3 -Propane Blends’ Oxidation

The reaction of NH 2 radicals with C 3 H 8 is crucial for understanding the combustion behavior of NH 3 /C 3 H 8 blends. In this study, we investigated the temperature dependence of the rate coefficients for the hydrogen abstraction reactions of C 3 H 8 by NH 2 radicals using high-level theoretical approaches. The potential energy surface was constructed at the CCSD(T)/cc-pV(T, Q)//M06-2X/aug-cc-pVTZ level of theory, and the rate coefficients were computed using conventional transition state theory, incorporating the corrections for quantum tunneling and hindered internal rotors (HIR). The computed rate coefficients showed a strong curvature in the Arrhenius behavior, capturing the experimental literature data well at low temperatures. However, at T > 1500 K, the theory severely overpredicted the experimental data. The available theoretical studies did not align with the experiment at high temperatures, and the possible reasons for this discrepancy are discussed. At 300 K, the reaction of NH 2 with C 3 H 8 predominantly occurs at the secondary C-H site, which accounts for approximately 95% of the total reaction flux. However, the hydrogen abstraction reaction at the primary C-H site becomes the dominant reaction above 1700 K. A composite kinetic model was built, which incorporated the computed rate coefficients for NH 2 + C 3 H 8 reactions. The importance of NH 2 + C 3 H 8 reactions in predicting the combustion behavior of NH 3 /C 3 H 8 blends was demonstrated by kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Opposed Jet Burner Extinction Limits: Simple Mixed Hydrocarbon Scramjet Fuels vs Air

Opposed Jet Burner tools have been used extensively by the authors to measure Flame Strength (FS) of laminar non-premixed H2 air and simple hydrocarbon (HC) air counterflow diffusion flames at 1-atm. FS represents a strain-induced extinction limit based on air jet velocity. This paper follows AIAA-2006-5223, and provides new HC air FSs for global testing of chemical kinetics, and for characterizing idealized flameholding potentials during early scramjet-like combustion. Previous FS data included six HCs, pure and N2-diluted; and three HC-diluted H2 fuels, where FS decayed very nonlinearly as HC was added to H2, due to H-atom scavenging. This study presents FSs on mixtures of (candidate surrogate) HCs, some with very high FS ethylene. Included are four binary gaseous systems at 300 K, and a hot ternary system at approx. 600 K. The binaries are methane + ethylene, ethane + ethylene, methane + ethane, and methane + propylene. The first three also form two ternary systems. The hot ternary includes both 10.8 and 21.3 mole % vaporized n-heptane and full ranges of methane + ethylene. Normalized FS data provide accurate means of (1) validating, globally, chemical kinetics for extinction of non-premixed flames, and (2) estimating (scaling by HC) the loss of incipient flameholding in scramjet combustors. The n-heptane is part of a proposed baseline simulant (10 mole % with 30% methane + 60% ethylene) that mimics the ignition of endothermically cracked JP-7 like kerosene fuel, as suggested by Colket and Spadaccini in 2001 in their shock tube Scramjet Fuels Autoignition Study. Presently, we use FS to gauge idealized flameholding, and define HC surrogates. First, FS was characterized for hot nheptane + methane + ethylene; then a hot 36 mole % methane + 64% ethylene surrogate was defined that mimics FS of the baseline simulant system. A similar hot ethane + ethylene surrogate can also be defined, but it has lower vapor pressure at 300 K, and thus exhibits reduced gaseous capacity. The new FS results refine our earlier idealized reactivity scale that shows wide ranging (50 x) diameter-normalized FSs for various HCs. These range from JP-10 and methane to H2 air, which produces an exceptionally strong flame that agrees within approx. 1% of recent 2-D numerically simulations. Finally, we continue advocating the FS approach as more direct and fundamental, for assessing idealized scramjet flameholding potentials, than measurements of unstrained laminar burning velocity or blowout in a Perfectly Stirred Reactor.

Pellett, Gerald L.↗

Liquid rocket combustor computer code development

The Advanced Rocket Injector/Combustor Code (ARICC) that has been developed to model the complete chemical/fluid/thermal processes occurring inside rocket combustion chambers are highlighted. The code, derived from the CONCHAS-SPRAY code originally developed at Los Alamos National Laboratory incorporates powerful features such as the ability to model complex injector combustion chamber geometries, Lagrangian tracking of droplets, full chemical equilibrium and kinetic reactions for multiple species, a fractional volume of fluid (VOF) description of liquid jet injection in addition to the gaseous phase fluid dynamics, and turbulent mass, energy, and momentum transport. Atomization and droplet dynamic models from earlier generation codes are transplated into the present code. Currently, ARICC is specialized for liquid oxygen/hydrogen propellants, although other fuel/oxidizer pairs can be easily substituted.

Liang, P. Y.↗