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 667 records · Page 37

Experimental and numerical investigation of ester droplet combustion: Application to butyl acetate

This paper presents an experimental and numerical study of the combustion of isolated n–butyl acetate droplets in the standard atmosphere. Numerical simulations are reported using a model that incorporates unsteady gas and liquid transport, variable properties, and radiation. Three skeletal mechanisms of n-butyl acetate, derived from a large detailed mechanism comprised of 819 species and 52,698 reactions, were used in the numerical simulations to evaluate the influence of the kinetic mechanism on burning. The reduced mechanisms comprised 212 species and 5413 reactions, 157 species and 3089 reactions, and 105 species and 1035 reactions. The numerical model did not include soot formation, though qualitatively mild sooting was noted only for droplets larger than 0.7 mm. The numerical predictions were in good agreement with experimental measurements of droplet and flame diameters. Flame extinction was numerically predicted which was attributed to a decrease of the characteristic diffusion time relative to the chemical time as droplet burned. Effects of initial droplet diameter on the evolution of maximum gas temperature (T max ) and peak mole fractions of CO 2 and CO are also examined numerically.

09 BIOMASS FUELS↗

Frhodo: A program for simulating chemical kinetic measurements and optimizing kinetic mechanisms

Chemical kinetic simulations are frequently used to extract kinetic and mechanistic information from experimental data by fitting simulations to the data. For simple systems of consecutive single channel reactions this is a relatively straightforward task. However, as the complexity of the system increases and particularly for systems involving competitive multi-channel reactions manually optimizing a mechanism to obtain the best fits to a range of experimental data quickly becomes a time-intensive and challenging process, especially when rate coefficients have pressure and temperature dependencies. There is considerable potential for exploiting automated optimization methods to rapidly screen reaction mechanisms and optimize the fit by adjusting rate coefficients within well-defined constraints. A new chemical kinetics simulation program, Frhodo, has been developed for this purpose. The program allows either manual simulation of individual experiments or fully automated optimization against a range of experimental data by adjusting user selected reactions within user defined constraints. Frhodo incorporates machine learning-based optimization routines that work by either minimizing standardized residuals or through Bayesian parameter estimation. Frhodo's optimization capabilities are demonstrated by reexamining two previously studied systems of relevance to combustion. Dissociation of diacetyl followed by recombination of methyl radicals is an example of a sequence of single channel, consecutive reactions. Pyrolysis of 2-methyl furan exemplifies multi-channel, unimolecular reactions that are often encountered and have temperature and pressure dependencies in the rate coefficients of the channels and branching ratios between them. For both systems the optimization routines resulted in solutions similar to the original studies. Finally, comments are made on the strengths and limitations of the approach.

0D simulation↗

Three Stage Cool Flame Droplet Burning Behavior of n-Alkane Droplets at Elevated Pressure Conditions: Hot, Warm and Cool Flame

Transient, isolated n-alkane droplet combustion is simulated at elevated pressure for helium-diluent substituted-air mixtures. We report the presence of unique quasi-steady, three-stage burning behavior of large sphero-symmetric n-alkane droplets at these elevated pressures and helium substituted ambient fractions. Upon initiation of reaction, hot-flame diffusive burning of large droplets is initiated that radiatively extinguishes to establish cool flame burning conditions in nitrogen/oxygen “air” at atmospheric and elevated pressures. However, at elevated pressure and moderate helium substitution for nitrogen ( X He > 20%), the initiated cool flame burning proceeds through two distinct, quasi-steady-state, cool flame burning conditions. The classical “Hot flame ” ( ∼1500 K) radiatively extinguishes into a “Warm flame ” burning mode at a moderate maximum reaction zone temperature ( ∼ 970 K), followed by a transition to a lower temperature ( ∼765 K), quasi-steady “Cool flame ” burning condition. The reaction zone (“flame”) temperatures are associated with distinctly different yields in intermediate reaction products within the reaction zones and surrounding near-field, and the flame-standoff ratios characterizing each burning mode progressively decrease. The presence of all three stages first appears with helium substitution near 20%, and the duration of each stage is observed to be strongly dependent on helium substitutions level between 20–60%. For helium substitution greater than 60%, the hot flame extinction is followed by only the lower temperature cool flame burning mode. In addition to the strong coupling between the diffusive loss of both energy and species and the slowly evolving degenerate branching in the low and negative temperature coefficient (NTC) kinetic regimes, the competition between the low-temperature chain branching and intermediate-temperature chain termination reactions control the “Warm” and “Cool ” flame quasi-steady conditions and transitioning dynamics.

Droplet Combustion↗

Investigation of Cycle-to-Cycle Variations in Internal Combustion Engine Using Proper Orthogonal Decomposition

The understanding, modeling and control of the cycle-to-cycle variation (CCV) in the modern internal combustion engine (ICE) is a key scientific challenge to achieve stable engine operation. High CCV in the engine combustion chamber may contribute to partial burn, misfire and knock, which adversely affects the engine performance and may potentially damage the engine. The objective of the current study is to leverage high-fidelity numerical simulations to improve the understanding of the causes of CCV. Using the massively parallel code, Nek5000, multi-cycle, wall-resolved large-eddy simulations (LES) were performed for the General Motors (GM), Transparent Combustion Chamber (TCC-III) optical engine under motored operating conditions. Further, the large-scale structures of the in-cylinder flow were investigated using a triple proper orthogonal decomposition (POD) technique to explore the characteristics of different parts of the flow and their contributions to CCV. The kinetic energy of the subset of flow structures were determined and correlated between the intake and compression strokes. The insights from the analysis of the large-scale flow structures will be used to assist the development of improved engine designs with reduced CCV and enhance the engine performance.

33 ADVANCED PROPULSION SYSTEMS↗

Characterization of a swept-strut hydrogen fuel-injector for scramjet applications

Results of an experimental investigation of a swept-strut hydrogen fuel-injector simulating the center strut of a three strut scramjet module at Mach 6 flight conditions are presented. Detailed wall pressure distributions from over 100 separate tests with overall fuel flow from 0.1 to 1.3 times stoichiometric and test gas stagnation temperature from 1100 to 2400 K were recorded. The distance for pressure rise from the point of injection was found to increase with increasing test stagnation temperature. This trend indicates that chemical kinetics in the immediate region of perpendicular injection are not likely to be the mechanism controlling the onset of pressure rise. A fluid dynamic mechanism is suggested involving separation of the boundary layer downstream of injection which is forced upstream from the trailing edge by pressure rise due to combustion occurring in the base region of the strut. The results obtained indicate that the swept-strut fuel-injector concept can be adapted to a wide range of flight conditions by varying the amount of perpendicular fuel injection.

Northam, G. B.↗

Combustion Characteristics of a Tier II Gasoline Certification Fuel and its Surrogate with Iso-butanol: experiments and detailed numerical modeling

This presentation will describe results of an experimental and numerical study of combustion of a Tier II gasoline certification fuel and its surrogate blended with iso-butanol as a potential biofuel additive. The experiments were carried out under conditions where gas transport originated entirely by fuel evaporation with external convection (both forced and natural) were minimized. The resulting near spherically symmetric droplet flames were amenable to detailed numerical modeling. The droplet initial diameters were in the narrow range of 0.55 mm to 0.62 mm while and burning is performed under standard atmosphere. The combustion characteristics such as evolution of droplet, flame and soot shell diameters for indolene/iso-butanol and surrogate/iso-butanol mixtures are compared. The presence of iso-butanol is experimentally found to reduce the formation of soot as the iso-butanol concentration was increased. Simulations are also reported using a detailed kinetic chemistry comprised of 324 species and 17608 reactions using the CRECK POLIMI mechanism. The detailed numerical model included unsteady gas and liquid transport, soot chemistry and radiative transport. Droplet burn rates were well predicted while flame diameters showed more variability depending on how the flame was defined in the numerical simulation (peak gas temperature and OH concentration). Simulations also show the influence of iso-butanol concentration on flame temperature, greenhouse gas emissions and particulate emissions.

Guo, Songtao↗

Editorial: Particulate Matter Emissions From Conventional and Reformulated Fuel Combustion: Advances in Experiments and Simulations

Particulate matter (PM or soot) emissions generated by combustion of conventional and reformulated fuels adversely impact the environment and human health (Bond et al., 2013; Landrigan et al., 2018; World Health Organization (WHO), 2016), generating public awareness and motivating efforts towards the mitigation of their harmful effects. Soot formation is one of the most complex phenomena in combustion, involving interactions between combustion chemistry, fluid mechanics, mass/heat transport, and particle dynamics, spanning different spatial and temporal scales. A fundamental understanding of soot formation process is thus necessary to achieve a strong reduction of PM emissions and design cleaner and more efficient combustion systems. This has motivated a longstanding and ongoing research activity aimed at improving our understanding of the physical and chemical processes involved in soot formation, well reviewed in these recent works (Wang and Chung, 2019; Michelsen et al., 2020; Martin et al., 2022). Despite the wide interest, the transition from gas-phase molecules to incipient soot particles is still elusive and the successive particle growth and oxidation processes are far from being fully understood, especially in conditions relevant to real-world applications. Recent advances in combustion PM emission diagnostic and computational capabilities helped in improving the predictability of fundamental chemical and aerosol models for practical applications, thus tackling some of the above-mentioned challenges. The aim of this research topic is to display the ongoing research efforts in addressing the existing gaps on particulate formation from various fuel sources and in conditions typical of practical combustion applications (e.g. flames, engines, pool fires), through both experimental and numerical approaches.

42 ENGINEERING↗

SHOCK TUBE IGNITION STUDIES OF RENEWABLE DIESEL FUELS FOR MEDIUM AND HEAVY-DUTY TRANSPORTATION

Currently extensive research on alternative fuels is being conducted due to their increasing demand to reduce greenhouse emissions. One renewable fuel studied in this work is dimethyl ether (DME) blended with propane(C3H8) as a potential mixture for heavy-duty engines used in semi-trucks. The blend has the potential to drastically reduce particulate and greenhouse gas emissions compared to a conventional diesel engine operating under similar conditions. To develop the use of mixture, one must conduct detailed conceptual and simulation studies before progressing to detail studies in CFD, engine modifications, and live testing. For simulations, accurate high-fidelity chemical kinetic models are necessary. However, the validity of the chemical kinetic mechanism for operating conditions of a heavy-duty mixing-controlled compression (MCCI) engine was widely unknown until recent work presented here and published. In this work, we studied the ignition of DME and propane blends in a shock tube under MCCI engine conditions. Ignition delay time (IDT) gathered behind the reflected shock for DME-propane mixtures for heavy-duty compression ignition (CI) engine parameters. Testing was conducted for undiluted varieties spanning from temperatures of 700 to 1100 K at pressures ranging from 55 to 84 bar for various blends (100% CH3OCH3, 100% C3H8, 60% CH3OCH3/ 40% C3H8) of DME and propane were combusted in synthetic air (21% O2/ 79% N2). Several experiments were conducted at higher pressures (90-120 bar) to improve the model performance and accuracy. The ignition delay times (IDTs) were compared to recent mechanisms, including Aramco3.0, NUIG, and Dames et al. A common trend among the mechanisms was overpredicted experimental IDTs. Further studies were conducted by a sensitivity analysis using the Dames et al. model, and critical reactions sensitive to IDTs of DME-propane mixture near 60 bar are outlined. Chemical analysis was conducted on the NTC region to explain chemical kinetics which is critical for developing MCCI heavy duty engines.

Mohammed, Zuhayr Pasha↗

The Effects of Buoyancy and Dilution on the Structure and Lift-off of Coflow Laminar Diffusion Flames

The ability to predict the coupled effects of complex transport phenomena with detailed chemical kinetics in diffusion flames is critical in the modeling of turbulent reacting flows and in understanding the processes by which soot formation and radiative transfer take place. In addition, an understanding of the factors that affect flame extinction in diffusion flames is critical in the suppression of fires and in improving engine efficiency. The goal of our characterizations of coflow laminar diffusion flames is to bring to microgravity the multidimensional diagnostic tools available in normal gravity, and in so doing provide a broader understanding of the successes and limitations of current combustion models. This will lead to a more detailed understanding of the interaction of convection, diffusion and chemistry in both buoyant and nonbuoyant environments. As a sensitive marker of changes in the flame shape, the number densities of excited-state CH (A(exp 2)delta, denoted CH*), and excited-state OH (A(exp 2)Sigma, denoted OH*) are measured in mu-g and normal gravity. Two-dimensional CH* and OH* number densities are deconvoluted from line-of-sight chemiluminescence measurements made on the NASA KC-135 reduced-gravity aircraft. Measured signal levels are calibrated, post-flight, with Rayleigh scattering. Although CH* and OH* kinetics are not well understood, the CH*, OH*, and ground-state CH distributions are spatially coincident in the flame anchoring region. Therefore, the ground-state CH distribution, which is easily computed, and the readily measured CH*/OH* distributions can be used to provide a consistent and convenient way of measuring lift-off height and flame shape in the diffusion flame under investigation. Given that the fuel composition affects flame chemistry and that buoyancy influences the velocity profile of the flow, we have the opportunity to computationally and experimentally study the roles of fluids and chemistry. In performing this microgravity study, improvements to the computational model have been made and new calculations performed for a range of gravity and flow conditions. Furthermore, modifications to the experimental approach were required as a consequence of the constraints imposed by existing microgravity facilities. Results from the computations and experiments are presented in the following sections.

Walsh, Kevin T.↗

Automated Reaction Kinetics of Gas-Phase Organic Species over Multiwell Potential Energy Surfaces

Automation of rate-coefficient calculations for gas-phase organic species became possible in recent years and has transformed how we explore these complicated systems computationally. Kinetics workflow tools bring rigor and speed and eliminate a large fraction of manual labor and related error sources. In this paper we give an overview of this quickly evolving field and illustrate, through five detailed examples, the capabilities of our own automated tool, KinBot. Here we bring examples from combustion and atmospheric chemistry of C-, H-, O-, and N-atom-containing species that are relevant to molecular weight growth and autoxidation processes. The examples shed light on the capabilities of automation and also highlight particular challenges associated with the various chemical systems that need to be addressed in future work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detonation structure in the presence of mixture stratification using reaction-resolved simulations

Many investigations of detonation-based combustors have identified reactant mixture inhomogeneity as having a leading-order impact on wave dynamics and combustion efficiency. To examine this phenomenon in a simplified context, an array of two- and three-dimensional channel detonation simulations are conducted in the present work. The reactant mixture consists of stratified fuel and air, wherein the randomly distributed equivalence ratio field features a characteristic stratification length scale. Detailed chemical kinetics are implemented in an adaptive mesh refinement solution framework where the region near the shock front is resolved with Ο (100) cells per representative ZND induction length. The results show that in comparison to baseline cases with uniform reactant mixtures, reactant stratification has a marked impact on the detonation structure. Increasing the stratification length scale increases the size and irregularity of the detonation cells, yielding larger variations in wave speed. Triple point collisions in fuel-rich regions lead to local wave speeds above the notional mean CJ speed, but wave passage through inert regions causes the local wave speed and strength to diminish. Further, conditional statistics show that increasing the stratification length scale increases the variance in pressure and temperature in the primary reaction zone, as well as the variance in heat release over a range of mixture conditions. In addition to the reactant mixture, the impact of the boundary condition behind the detonation is also investigated. The results show that an inflow boundary condition acts to over-drive the wave, leading to higher peak pressures, smaller detonation cells, and increased reactant consumption. On the other hand, cases with a wall behind the wave exhibit weaker waves with lower peak pressures and heat release rates, as well as greater variance in conditional quantities. Comparisons between complementary two- and three-dimensional simulations show reasonable qualitative agreement in wave structure, speed, and conditional statistics.

42 ENGINEERING↗

Particle Size Distribution and Its Impacts on Ash Deposition and Radiative Transfer during Oxy-Combustion of Rice Husk–Natural Gas

Rice husk (RH) co-combustion with natural gas in highly oxygen-enriched concentrations presents a net carbon-negative energy production opportunity while minimizing flue gas recycling. However, recent experiments have shown enhanced ash deposition rates in oxygen-enriched conditions, with deposition/shedding also being dependent on the particle size distribution (PSD) of the parent RH fuel. To uncover the causative mechanisms behind these observations, add-on models for ash deposition/shedding and radiative properties were employed in computational fluid dynamics simulations. The combustion scenarios investigated encompassed two types of RH (US RH, Chinese RH) with widely varying ash contents (by % mass) and inlet fuel PSD with air and O2/CO2 (70/30 vol %, OXY70) as oxidizers. Utilizing the measured fly-ash PSDs near the deposit surface and modeling the particle viscosity accurately, particle kinetic-energy (PKE)-based capture and shedding criteria were identified as the keys to accurate deposition/shedding rate predictions. The OXY70 scenarios showed higher ash-capturing propensities due to their lower PKE. Conversely, higher erosion rates were predicted in the AIR firing scenarios. In addition, the radiative characteristics across all the scenarios were dominated by the gases and were not sensitive to the fly-ash PSD. Therefore, the higher particle concentrations in the OXY70 conditions did not negatively impact the heat extraction.

Krishnamoorthy, Gautham (ORCID:0000000255205092)↗

Simultaneous Laser Raman-rayleigh-lif Measurements and Numerical Modeling Results of a Lifted Turbulent H2/N2 Jet Flame in a Vitiated Coflow

An experiment and numerical investigation is presented of a lifted turbulent H2/N2 jet flame in a coflow of hot, vitiated gases. The vitiated coflow burner emulates the coupling of turbulent mixing and chemical kinetics exemplary of the reacting flow in the recirculation region of advanced combustors. It also simplifies numerical investigation of this coupled problem by removing the complexity of recirculating flow. Scalar measurements are reported for a lifted turbulent jet flame of H2/N2 (Re = 23,600, H/d = 10) in a coflow of hot combustion products from a lean H2/Air flame ((empty set) = 0.25, T = 1,045 K). The combination of Rayleigh scattering, Raman scattering, and laser-induced fluorescence is used to obtain simultaneous measurements of temperature and concentrations of the major species, OH, and NO. The data attest to the success of the experimental design in providing a uniform vitiated coflow throughout the entire test region. Two combustion models (PDF: joint scalar Probability Density Function and EDC: Eddy Dissipation Concept) are used in conjunction with various turbulence models to predict the lift-off height (H(sub PDF)/d = 7,H(sub EDC)/d = 8.5). Kalghatgi's classic phenomenological theory, which is based on scaling arguments, yields a reasonably accurate prediction (H(sub K)/d = 11.4) of the lift-off height for the present flame. The vitiated coflow admits the possibility of auto-ignition of mixed fluid, and the success of the present parabolic implementation of the PDF model in predicting a stable lifted flame is attributable to such ignition. The measurements indicate a thickened turbulent reaction zone at the flame base. Experimental results and numerical investigations support the plausibility of turbulent premixed flame propagation by small scale (on the order of the flame thickness) recirculation and mixing of hot products into reactants and subsequent rapid ignition of the mixture.

Cabra, R.↗

Low Temperature Oxidation of Methylpropyl Ether

Ethers are potential high-performance fuels that can be produced from renewable carbon sources such as biomass, but their combustion properties are not clearly understood. The cetane numbers (CN) for ethers, a measure of autoignition tendency, can range from less than 10 to greater than 100, and predictive tools based upon the molecular structures are desirable in order to focus fuel development research. In this study, we used molecular modeling, chemical kinetic modeling, and experimental measurements to understand the important reactions that lead to the autoignition of methylpropyl ether (MPE), a model for ethers containing alkyl chains long enough for intramolecular hydrogen transfer. These reactions are known to be important in the autoignition of alkanes and ethers and we propose to investigate the relationship between alkyl and ether reaction pathways. A reaction mechanism for MPE was created using MIT's Reaction Mechanism Generator (RMG) and G4 quantum calculation results and models were tested using the ChemKin Pro suite of software. Results from atmospheric pressure, flow tube reactor experiments and rapid compression machine (RCM) experiments were used to test the kinetic model and suggest refinements.

33 ADVANCED PROPULSION SYSTEMS↗

The Reaction of o-Benzyne with Vinylacetylene: An Unexplored Way to Produce Naphthalene

Herein, the mechanisms and kinetics of the reaction of ortho-benzyne with vinylactylene have been studied by ab initio and density functional CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G(d,p) calculations of the pertinent potential energy surface combined with Rice-Ramsperger-Kassel-Marcus - Master Equation calculations of reaction rate constants at various temperatures and pressures. Under prevailing combustion conditions, the reaction has been shown to predominantly proceed by the biradical acetylenic mechanism initiated by the addition of C4H4 to one of the C atoms of the triple bond in ortho-benzyne by the acetylenic end, with a significant contribution of the concerted addition mechanism. Following the initial reaction steps, an extra six-membered ring is produced and the rearrangement of H atoms in this new ring leads to the formation of naphthalene, which can further dissociate to 1- or 2-naphthyl radicals. The o-C 6 H 4 + C 4 H 4 reaction is highly exothermic, by ~143 kcal/mol to form naphthalene and by 31-32 kcal/mol to produce naphthyl radicals + H, but features relatively high entrance barriers of 9-11 kcal/mol. Although the reaction is rather slow, much slower than the reaction of phenyl radical with vinylacetylene, it forms naphthalene and 1- and 2-naphthyl radicals directly, with their relative yields controlled by the temperature and pressure, and thus represents a viable source of the naphthalene core under conditions where ortho-benzyne and vinylacetylene are available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hybrid propulsion system with a gyro component for economic and dynamic operation

The design of a hybrid drive with gyro components is described and its drive components for a medium class private car are discussed. The gyro component affects the short-period output of the drive by accelerating and slowing down and -- because of the mechanical transfer of kinetic energy between the gyro and the vehicle -- it affects also the energy balance in the case of intermittent operation. Energy can be taken in as desired either in the form of fuel or as fuel and current. A high energy recovery efficiency as well as the favorable operating range of the interval combustion engine makes it possible to reduce the fuel consumption per unit distance travelled to almost half that for a private car with a traditional engine.

Giera, B.↗

High-Area-Ratio Rocket Nozzle at High Combustion Chamber Pressure: Experimental and Analytical Validation

Experimental data were obtained on an optimally contoured nozzle with an area ratio of 1025:1 and on a truncated version of this nozzle with an area ratio of 440:1. The nozzles were tested with gaseous hydrogen and liquid oxygen propellants at combustion chamber pressures of 1800 to 2400 psia and mixture ratios of 3.89 to 6.15. This report compares the experimental performance, heat transfer, and boundary layer total pressure measurements with theoretical predictions of the current Joint Army, Navy, NASA, Air Force (JANNAF) developed methodology. This methodology makes use of the Two-Dimensional Kinetics (TDK) nozzle performance code. Comparisons of the TDK-predicted performance to experimentally attained thrust performance indicated that both the vacuum thrust coefficient and the vacuum specific impulse values were approximately 2.0-percent higher than the turbulent prediction for the 1025:1 configurations, and approximately 0.25-percent higher than the turbulent prediction for the 440:1 configuration. Nozzle wall temperatures were measured on the outside of a thin-walled heat sink nozzle during the test fittings. Nozzle heat fluxes were calculated front the time histories of these temperatures and compared with predictions made with the TDK code. The heat flux values were overpredicted for all cases. The results range from nearly 100 percent at an area ratio of 50 to only approximately 3 percent at an area ratio of 975. Values of the integral of the heat flux as a function of nozzle surface area were also calculated. Comparisons of the experiment with analyses of the heat flux and the heat rate per axial length also show that the experimental values were lower than the predicted value. Three boundary layer rakes mounted on the nozzle exit were used for boundary layer measurements. This arrangement allowed total pressure measurements to be obtained at 14 different distances from the nozzle wall. A comparison of boundary layer total pressure profiles and analytical predictions show good agreement for the first 0.5 in. from the nozzle wall; but the further into the core flow that measurements were taken, the more that TDK overpredicted the boundary layer thickness.

Jankovsky, Robert S.↗