Numerical Model Development for Supercritical-CO2 Oxy-Combustion
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
A simplified kinetic scheme for Jet-A, and methane fuels with water injection was developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) or even simple FORTRAN codes. The two time step method is either an initial time averaged value (step one) or an instantaneous value (step two). The switch is based on the water concentration in moles/cc of 1x10(exp -20). The results presented here results in a correlation that gives the chemical kinetic time as two separate functions. This two time step method is used as opposed to a one step time averaged method previously developed to determine the chemical kinetic time with increased accuracy. The first time averaged step is used at the initial times for smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used with higher water concentrations, gives the chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure and temperature (T4). The simple correlations would then be compared to the turbulent mixing times to determine the limiting rates of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide and NOx are obtained for Jet-A fuel and methane with and without water injection to water mass loadings of 2/1 water to fuel. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, and pressure.
A simplified kinetic scheme for Jet-A, and methane fuels with water injection was developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) or even simple FORTRAN codes that are being developed at Glenn. The two time step method is either an initial time averaged value (step one) or an instantaneous value (step two). The switch is based on the water concentration in moles/cc of 1x10(exp -20). The results presented here results in a correlation that gives the chemical kinetic time as two separate functions. This two step method is used as opposed to a one step time averaged method previously developed to determine the chemical kinetic time with increased accuracy. The first time averaged step is used at the initial times for smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used with higher water concentrations, gives the chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure and temperature (T4). The simple correlations would then be compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates were then used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide and NOx were obtained for Jet-A fuel and methane with and without water injection to water mass loadings of 2/1 water to fuel. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, and pressure.
This paper focuses on the development and validation of a combustion model for Computational Fluid Dynamics (CFD) modeling of Rotating Detonation Engines. A zero-dimensional Partially Stirred Reactor (PaSR) with a detailed chemical kinetic mechanism for hydrogen and air is used to model turbulent combustion. The model is computationally efficient and is based on the notion of partial mixing at the sub-grid level with turbulent exchange between mixed and unmixed regions. The ability of the PaSR model to accurately represent both detonative and deflagrative combustion is assessed by validating the results against experimental data. The effects of mesh resolution on the solution are also studied in order to determine if a mesh independent solution is obtainable with the Large Eddy Simulation (LES) approach to modeling turbulence. A comparison is made between the PaSR model and simply ignoring turbulence chemistry interactions which assumes that all species are perfectly mixed at the sub-grid level.
Dimethyl ether (DME) is a biofuel that has the potential to replace diesel in heavy-duty engines. A blend of propane (C3H8) and DME could reduce emissions when compared to diesel at heavy-duty compression engine-relevant conditions. Testing a potential new mixture in a compression engine requires a high-fidelity chemical kinetics model that accurately predicts autoignition delay times, which is especially necessary for a compression engine. In this work, autoignition data has been gathered at an equivalence ratio of 2.0 for pressures of 60 and 80 bar. DME and propane were combusted in synthetic air while excited hydroxyl (OH*) chemiluminescence was used to gather ignition delay times. Data was compared to recent chemical kinetic mechanisms for the two different pressures. Using the mechanism with the best fit of autoignition data, a sensitivity analysis was conducted to analyze the chemical kinetics of the DME/C3H8 combustion at elevated pressures.
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.
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.
A computational simulation of reacting 2-D and 3-D flowfields in a model inlet section of a Pre-Mixed, Shock-Induced Combustion (PM/SIC) engine concept was performed. LARCK, a multi-dimensional Navier-Stokes code with finite-rate kinetics chemistry developed at NASA LaRC by J.A. White, was adapted for this simulation. The flow conditions in the simulation match those envisioned for the PM/SIC engine experiments currently planned at LaRC. The reacting flowfields were Mach 6.3 freestream air and Mach 2 hydrogen at various pressure and temperature conditions injected through a slot injector at the base of the inlet section. In the PM/SIC engine, fuel is injected at the inlet section upstream of the combustor, and reaction is initiated by the shock wave at the inlet which increases the gas temperature and pressure beyond the kinetic limits for reaction. Many challenges exist prior to establishing shock-controlled combustion as a practical engine concept. These challenges include fuel injection schemes that can provide proper fuel-air mixing without creating large losses in the inlet section, and control of the combustion process so that early ignition or combustion propagation through the inlet boundary layer does not occur. For this project, a parametrics study was carried out to model the fuel injection of hydrogen at different flow conditions. It was found that, as the fuel temperature and pressure were increased, the potential for pre-ignition was high at a short distance downstream of the slot injector. The next stage of this work will investigate injection techniques for enhancing mixing of fuel and air in a manner that prevents or reduces the potential for premature ignition observed numerically.
Nitrogen-containing organic molecules play essential roles in nearly every aspect of chemistry. The importance of these species within gas-phase chemistry has increased in recent years, due to their emissions from wildfires, crude biofuels, ammonia combustion, and CO 2 capture facilities. However, there is a lack of detailed relevant mechanistic and kinetic studies of N-compounds, resulting in poor representation of these pathways within chemical models. In this report, we detail our exploration of the gas-phase oxidation mechanisms of pyrrole, imidazole, pyrrolidine, methylamine, and dimethylamine, done by employing multiplexed VUV photoionization mass spectrometry coupled with KinBot, a computational tool that automatically explores multi-well potential energy surfaces. Overall, we vastly expanded the known gas-phase oxidation mechanisms and chemical kinetics of these species. This work will improve the accuracy and completeness of atmospheric and combustion chemistry models, ultimately leading to better-informed decision-makers in the energy and environmental sectors.
At flight speeds, the residence time for atmospheric air ingested into a scramjet inlet and exiting from the engine nozzle is on the order of a millisecond. Therefore, fuel injected into the air must efficiently mix within tens of microseconds and react to release its energy in the combustor. The overall combustion process should be mixing controlled to provide a stable operating environment; in reality, however, combustion in the upstream portion of the combustor, particularly at higher Mach numbers, is kinetically controlled where ignition delay times are on the same order as the fluid scale. Both mixing and combustion time scales must be considered in a detailed study of mixing and reaction in a scramjet to understand the flow processes and to ultimately achieve a successful design. Although the geometric configuration of a scramjet is relatively simple compared to a turbomachinery design, the flow physics associated with the simultaneous injection of fuel from multiple injector configurations, and the mixing and combustion of that fuel downstream of the injectors is still quite complex. For this reason, many researchers have considered the more tractable problem of a spatially developing, primarily supersonic, chemically reacting mixing layer or jet that relaxes only the complexities introduced by engine geometry. All of the difficulties introduced by the fluid mechanics, combustion chemistry, and interactions between these phenomena can be retained in the reacting mixing layer, making it an ideal problem for the detailed study of supersonic reacting flow in a scramjet. With a good understanding of the physics of the scramjet internal flowfield, the designer can then return to the actual scramjet geometry with this knowledge and apply engineering design tools that more properly account for the complex physics. This approach will guide the discussion in the remainder of this section.
The fundamental physico-chemical mechanisms governing droplet vaporization and combustion are discussed. Specific topics include governing equations and simplifications, the classical d(2)-Law solution and its subsequent modification, finite-rate kinetics and the flame structure, droplet dynamics, near- and super-critical combustion, combustion of multicomponent fuel blends/emulsions/suspensions, and droplet interaction. Potential research topics are suggested.
Cool-flame-supported quasi-steady droplet combustion was first observed during droplet combustion experiments conducted onboard the International Space Station (ISS) in the Combustion Integrated Rack (CIR) with n-alkane fuels. These cool flames, not visible to the naked eye, result from the existence of the negative-temperature-coefficient (NTC) region during the combustion n-alkane fuels when the chemical-kinetic pathways experience a transition from the low-temperature chemistry to the normal, high-temperature-controlled chemistry. In this study, we present results for the cool-flame-supported droplet burning rates and extinction diameters for n-alkane fuels (n-heptane, n-octane, n-decane, and n-dodecane) in inert-diluent-substituted ambient environments. In these ISS experiments, some of the nitrogen in the ambient gas mixture is replaced by an inert gas, namely, carbon dioxide, helium, or xenon. Different diluents were observed to exert remarkably different influences on the cool-flame combustion of fuel droplets. The results show that the cool-flame extinction diameter is drastically increased with helium dilution compared to N2, CO2 or Xe, and Xe dilution leads to longer burning with a smaller extinction diameter. Unlike hot flames, the cool-flame burning rates are controlled by the cool-flame temperature established in the NTC region by the chemical kinetics, rather than by equilibrium thermodynamics, and the thermo-physical properties of the ambient gases, which control the heat feed back to the liquid droplet. The cool-flame extinction is triggered by the instability that occurs at the low-temperature end of the NTC region, below which cool diffusion flames become statically unstable. A quasi-steady theoretical model, presented earlier, which involves a six-step reduced-chemistry mechanism, valid around the lower end of the NTC region, is shown to explain the observed experimental trends qualitatively.
In order to develop an extensive envelope of test conditions for NASA's space-based Droplet Combustion Experiment (DCE) as well those droplet experiments which can be performed using a drop tower, the transient vaporization and combustion of methanol and n-heptane droplets were simulated using a recently developed fully time-dependent, spherically symmetric droplet combustion model. The transient vaporization of methanol and n-heptane was modeled to characterize the instantaneous gas phase composition surrounding the droplet prior to the introduction of an ignition source. The results for methanol/air showed that the entire gas phase surrounding a 2 mm methanol droplet deployed in zero-g .quickly falls outside the lean flammability limit. The gas phase surrounding an identically-sized n-heptane droplet, on the other hand, remains flammable. The combustion of methanol was then modeled considering a detailed gas phase chemical kinetic mechanism (168 steps, 26 species) and the effect of the dissolution of flame-generated water into the liquid droplet. These results were used to determine the critical ignition diameter required to achieve quasi-steady droplet combustion in a given oxidizing environment. For droplet diameters greater than the critical ignition diameter, the model predicted a finite diameter at which the flame would extinguish. These extinction diameters were found to vary significantly with initial droplet diameter. This phenomenon appears to be unique to the transient heat transfer, mass transfer and chemical kinetics of the system and thus has not been reported elsewhere to date. The extinction diameter was also shown to vary significantly with the liquid phase Lewis number since the amount of water present in the droplet at extinction is largely governed by the rate at which water is transported into the droplet via mass diffusion. Finally, the numerical results for n-heptane combustion were obtained using both 2 step and 96 step semi-emperical chemical kinetic mechanisms. Neither mechanism exhibited the variation of extinction diameter with initial diameter.
Work focussed on a large number of preliminary studies of supersonic combustion in a simple combustion duct - thrust nozzle combination, investigating effects of Mach number, equivalence ratio, combustor divergence, fuel injecting angle and other parameters with an influence on the combustion process. This phase lasted for some three or four years, during which strongest emphasis was placed on responding to the request for preliminary experimental information on high enthalpy effects, to support the technology maturation activities of the NASP program. As the need for preliminary data became less urgent, it was possible to conduct more systematic studies of high enthalpy combustion phenomena, and to initiate other projects aimed at improving the facilities and instrumentation used for studying scramjet phenomena at high enthalpies. The combustion studies were particularly directed towards hypersonic combustion, and to the effects of injecting fuel along the combustion chamber wall. A substantial effort was directed towards a study of the effect of scale on the supersonic combustion process. The influence of wave phenomena (both compression waves and expansion waves) on the realization of thrust from a supersonic combustion process was also investigated. The effect of chemical kinetics was looked into, particularly as it affected the composition of the test flow provided by a ground facility. The effect of injection of the fuel through wall orifices was compared with injection from a strut spanning the stream, and the effect of heating the fuel prior to injection was investigated. Studies of fuel-air mixing by shock impingement were also done, as well as mass spectrometer surveys of a combustion wake. The use of hypersonic nozzles with an expansion tube was investigated. A new method was developed for measuring the forces acting of a model in less than one millisecond. Also included in this report are listings of published journal papers and conference presentations.
Abstract Accurate, predictive reaction models are critical for the design and optimization of chemical looping combustion (CLC) reactors. The formulation and estimation of kinetic parameters for these reaction models using a first‐principles equation‐oriented (EO) approach is particularly beneficial as large amounts of experimental data spanning process‐relevant conditions can be used to estimate parameters in a computationally tractable way. This work demonstrates the application of a novel EO framework to develop reduction reaction kinetic models of an iron‐based CLC oxygen carrier (OC). An optimization problem is formulated to estimate kinetic parameters that provide the best fit to the experimental data. The model predicts the state of the OC with mean square error values of 2.5%–4.4% across the full range of validation data, including multiple reduction cycles.
Surrogate fuels that reproduce the characteristics of full-boiling range fuels are key tools to enable numerical simulations of fuel-related processes and ensure reproducibility of experiments by eliminating batch-to-batch variability. Within the PACE initiative, a surrogate fuel for regular-grade E10 (10%vol ethanol) gasoline representative of a U.S. market gasoline, termed PACE-20, was developed and adopted as baseline fuel for the consortium. Although extensive testing demonstrated that PACE-20 replicates the properties and combustion behavior of the full-boiling range gasoline, several concerns arose regarding the purity level required for the species that compose PACE-20. This is particularly important for cyclo-pentane, since commercial-grade cyclo-pentane typically shows 60%–85% purity. Here, in the present work, the effects of the purity level of cyclo-pentane on the properties and combustion characteristics of PACE-20 were studied. Chemical kinetic simulations were performed to predict the effects of cyclo-pentane impurities on the properties, octane rating, and autoignition reactivity under homogeneous charge compression-ignition conditions of PACE-20. From the numerical results, cyclo-pentane with 85% purity or higher is required to reasonably match both the research octane number and motor octane number of the target gasoline. Finally, homogeneous charge compression-ignition engine simulations show that impurities have only a modest effect on reactivity at naturally aspirated conditions, but cyclo-pentane purity is critical to properly replicate the pressure dependency of the reactivity.
Questions of the importance and feasibility of performing experiments on droplet burning at zero gravity in Spacelab were studied. Information on the physics and chemistry of droplet combustion, with attention directed specifically to the chemical kinetics, heat and mass transfer, and fluid mechanics of the phenomena involved, are presented. The work was divided into three phases, the justification, the feasibility, and the conceptual development of a preliminary design. Results from the experiments performed revealed a few new facts concerning droplet burning, notably burning rates in excess of theoretical prediction and a phenomenon of flash extinction, both likely traceable to accumulation of carbon produced by gas-phase pyrolysis in the fuel-rich zone enclosed by the reaction surface. These experiments also showed that they were primarily due to timing difficulties.
A scale-invariant model of statistical mechanics is applied to present invariant forms of mass, energy, linear, and angular momentum conservation equations in reactive fields. The resulting conservation equations at molecular-dynamic scale are solved by the method of large activation energy asymptotics to describe the hydro-thermo-diffusive structure of laminar premixed flames. The predicted temperature and velocity profiles are in agreement with the observations. Also, with realistic physico-chemical properties and chemical-kinetic parameters for a single-step overall combustion of stoichiometric methane-air premixed flame, the laminar flame propagation velocity of 42.1 cm/s is calculated in agreement with the experimental value.