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

Examining experimental nitrogen-based emissions trends from ammonia/diesel and ammonia/hydrogen/diesel combustion

Ammonia has garnered interest as an alternative fuel for power sectors with heavy payload and distance requirements, such as shipping. In this study, ammonia was used in a dual-fuel compression-ignition combustion strategy to overcome some of its technical barriers, using a diesel pilot to ignite a premixed mixture of ammonia and air. Mixtures of premixed ammonia and hydrogen were also explored to evaluate whether the inclusion of hydrogen improves nitrogen-based emissions from the combustion process. A single-cylinder version of a Cummins ISB 6.7 L engine platform was used to experimentally study these effects at various global air/fuel ratios and hydrogen energy fractions. Hydrogen inclusion produced pronounced NOx and N2O emissions, while inclusion of trapped residuals increased N2O but reduced NOx. The two most recent and relevant mechanisms available in the literature—those from Xu and Zhang-Ren-Kokjohn—were used in a chemical kinetics analysis to examine the observed differences in the NOx trends from two dual-fuel ammonia/diesel datasets: (1) in which a portion of the premixed ammonia was substituted with hydrogen and (2) in which the effect of hot trapped residuals was evaluated with only ammonia/air premixed mixtures. The analysis with both mechanisms showed agreement with experimental trends; however, contributions from thermal vs. fuel-borne NOx pathways showed disagreement. A reaction pathway analysis showed that the HNO to NO pathway was the key to NO formation in the mixture.

Tyrewala, Daanish [ORNL] (ORCID:0000000208599324)↗

Isomer-Dependent Reaction Mechanisms of Cyclic Ether Intermediates: cis -23-Dimethyloxirane and trans -23-Dimethyloxirane.

Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β-hydroperoxyalkyl radicals (β-˙QOOH) or from reactions of HÒO with alkenes. The cis- and trans-isomers of 2,3-dimethyloxirane are intermediates of n-butane oxidation, and while rate coefficients for β-˙QOOH → 2,3-dimethyloxirane + OH are reported extensively, subsequent reaction mechanisms of the cyclic ethers are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe the consumption of 2,3-dimethyloxirane by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. The present research examines the isomerdependence of 2,3-dimethyloxirane reaction mechanisms in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred via the detection of products from Cl-initiated oxidation of both cis-2,3-dimethyloxirane and trans-2,3-dimethyloxirane using multiplexed photoionization mass spectrometry (MPIMS). The experiments were conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, the enthalpies of stationary points on the ˙R + O 2 surfaces were computed at the ccCA-PS3 level of theory. In total, 28 barrier heights were computed on the 2,3-dimethyloxiranylperoxy surfaces. Two notable aspects are low-lying pathways that form resonance-stabilized ketohydroperoxide-type radicals caused by ˙QOOH ring-opening when the unpaired electron is localized adjacent to the ether group, and cis-trans isomerization of ˙R and ˙QOOH radicals, via inversion, which enable reaction pathways otherwise restricted by stereochemistry. Several species were identified in the MPIMS experiments from ring opening of 2,3-dimethyloxiranyl radicals. Neither of the two conjugate alkene isomers prototypical of ˙R + O 2 reactions were detected. Products were also identified from decomposition of ketohydroperoxide-type radicals. The present work provides the first analysis of 2,3-dimethyloxirane oxidation chemistry and reveals that consumption pathways are complex and require the expansion of submechanisms in chemical kinetics mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shock tube investigation of high-temperature, extremely-rich oxidation of several co-optima biofuels for spark-ignition engines

To reduce the reliance on fossil fuels in the transportation sector and increase combustion efficiency, the Co-Optima initiative from US Department of Energy identified the top 10 biofuels for downsized, boosted, spark-ignition engines. Most of these biofuels have detailed reaction mechanisms available in literature developed based on studies at temperatures lower than 1700 K and an equivalence ratio of less than five. As such, the performance of these detailed mechanisms at high temperature and extremely rich conditions are unknown. It is important to validate kinetic mechanisms at these conditions because they are conducive to soot formation. Prediction of soot by chemical kinetic models relies on the prediction of underlying benchmark species like carbon monoxide and ethylene. In this work, we conduct high temperature (1700–2050 K) and high equivalence ratio(Φ=8.6) oxidation of these biofuels, namely 2,4,4-trimethyl-1-pentene (α-diisobutylene), ethanol, cyclopentanone, methyl acetate, and 2-methylfuran, blended in ethylene behind reflected shock waves at 4–4.7 atm pressure. Carbon monoxide and ethylene time histories are measured simultaneously using a continuous feedback quantum cascade laser near 4.9 µm and a tunable CO 2 gas laser at 10.532 nm, respectively. Results show that ethanol blend forms more carbon monoxide than other biofuels and consumes ethylene faster than the biofuel blends in the temperature range considered. Overall, the performance of different mechanisms in literature are evaluated against the experimental results. The novel reaction mechanism ‘the Co-Optima model’, which includes the sub mechanisms for all the biofuels in this study, was found to be the best mechanism for the experimental conditions studied.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An Experimental Kinetics Study of Isopropanol Pyrolysis and Oxidation behind Reflected Shock Waves

Isopropanol has potential as a future bio-derived fuel and is a promising substitute for ethanol in gasoline blends. Even so, little has been done in terms of high-temperature chemical kinetic speciation studies of this molecule. To this end, experiments were conducted in a shock tube using simultaneous CO and H2O laser absorption measurements. Water and CO formation during isopropanol pyrolysis was also examined at temperatures between 1127 and 2162 K at an average pressure of 1.42 atm. Species profiles were collected at temperatures between 1332 and 1728 K and at an average pressure of 1.26 atm for equivalence ratios of 0.5, 1.0, and 2.0 in highly diluted mixtures of 20% helium and 79.5% argon. Species profiles were also compared to four modern C3 alcohol mechanisms, including the impact of recent rate constant measurements. The Li et al. (2019) and Saggese et al. (2021) models both best predict CO and water production under pyrolysis conditions, while the AramcoMech 3.0 and Capriolo and Konnov models better predict the oxidation experimental profiles. Additionally, previous studies have collected ignition delay time (τign) data for isopropanol but are limited to low pressures in highly dilute mixtures. Therefore, real fuel–air experiments were conducted in a heated shock tube with isopropanol for stoichiometric and lean conditions at 10 and 25 atm between 942 and 1428 K. Comparisons to previous experimental results highlight the need for real fuel–air experiments and proper interpretation of shock-tube data. The AramcoMech 3.0 model over predicts τign values, while the Li et al. model severely under predicts τign. The models by Capriolo and Konnov and Saggese et al. show good agreement with experimental τign values. A sensitivity analysis using these two models highlights the underlying chemistry for isopropanol combustion at 25 atm. Additionally, modifying the Li et al. model with a recently measured reaction rate shows improvement in the model’s ability to predict CO and water profiles during dilute oxidation. Finally, a regression analysis was performed to quantify τign results from this study.

09 BIOMASS FUELS↗

Infrared Fingerprint and Unimolecular Decay Dynamics of the Hydroperoxyalkyl Intermediate (•QOOH) in Cyclopentane Oxidation

A transient carbon-centered hydroperoxyalkyl intermediate (•QOOH) in the oxidation of cyclopentane is identified by IR action spectroscopy with time-resolved unimolecular decay to hydroxyl (OH) radical products that are detected by UV laser-induced fluorescence. Two nearly degenerate •QOOH isomers, β- and γ-QOOH, are generated by H atom abstraction of the cyclopentyl hydroperoxide precursor. Fundamental and first overtone OH stretch transitions and combination bands of •QOOH are observed and compared with anharmonic frequencies computed by second-order vibrational perturbation theory. An OH stretch transition is also observed for a conformer arising from torsion about a low-energy CCOO barrier. Definitive identification of the β-QOOH isomer relies on its significantly lower transition state (TS) barrier to OH products, which results in rapid unimolecular decay and near unity branching to OH products. A benchmarking approach is utilized to compute high-accuracy stationary point energies, most importantly TS barriers, for cyclopentane oxidation (C 5 H 9 O 2 ), building on higher level reference calculations for ethane oxidation (C 2 H 5 O 2 ). The experimental OH product appearance rates are compared with computed statistical microcanonical rates using RRKM theory, including heavy-atom tunneling, thereby validating the computed TS barrier. The results are extended to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclopentane combustion via master-equation modeling. Furthermore, the various torsional and ring puckering states of the wells and transition states are explicitly considered in these calculations.

Chemical calculations↗

Electrocatalytic water oxidation performance in an extended porous organic framework with a covalent alliance of distinct Ru sites

The rational synthesis of durable, earth-abundant efficient electrocatalysts for the oxygen evolution reaction (OER) from water is one of the most important routes for storing renewable energy and minimizing fossil fuel combustion. The prime hurdles for effectively utilizing commercial RuO 2 as (OER) electrocatalysts are its very low stability, catalyst deactivation, and high cost. In this work, we explored a Ru-integrated porous organic polymer (Ru@Bpy-POP) by a facile one-pot Friedel–Crafts alkylation strategy between redox-active (Ru(demob) 3 Cl 2 ) and a carbazole unit, which is composed of unique features including an extended framework unit, isolated active sites, and tunable electrode kinetics. Ru@Bpy-POP can serve as a bridge between a Metal–Organic Framework (MOF) and POP-based catalytic systems with a balanced combination of covalent bonds (structural stability) and open metal sites (single site catalysis). Ru@Bpy-POP, deposited on a three-dimensional nickel foam electrode support, exhibits a promising electrocatalytic OER activity with an ultra-low ruthenium loading compared to a benchmark RuO 2 catalyst, providing an overpotential of about 270 mV to reach 10 mA cm –2 in an alkaline medium. Moreover, a high current density of 248 mA cm –2 was achieved for the Ru@Bpy-POP catalyst at only 1.6 V (vs. RHE), which is much higher than 91 mA cm –2 for commercial RuO 2 . The robust, albeit highly conjugated, POP framework not only triggered facile electro-kinetics but also suppressed aggregation and metallic corrosion during electrolysis. In particular, the benefits of covalent integration of distinct Ru sites into the framework can modulate intermediate adsorption and charge density, which contributes to its exceptional OER activity. All of the critical steps involved in OER are complemented by Density Functional Theory (DFT) calculations, which suggest that electrocatalytic water oxidation proceeds from a closed-shell configuration to open-shell electronic configurations with high-spin states. Finally, these open-shell configurations are more stable than their closed-shell counterparts by 1 eV, improving the overall catalytic activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗

An Experimental and Numerical Investigation to Characterize the Low-Temperature Heat Release in Stoichiometric and Lean Combustion

This work reports on an experimental and modeling study on the low-temperature heat release (LTHR) characteristics for three RON 90 binary blends (n-heptane blended with isooctane, toluene and ethanol) in a Cooperative Fuel Research (CFR) engine at lean and stoichiometric conditions that are representative of homogeneous charge compression ignition (HCCI) and spark-ignition (SI) end-gas combustion conditions, respectively. An analysis of the end-gas temperature-pressure (T-P) trajectories was performed to identify the intake conditions leading to similar T-P trajectories between the two lambdas for each fuel blend. A heat release analysis was then conducted for the identified cases, where fuel-to-fuel differences in LTHR were identified and found to be sensitive to the operating condition. Simulations were conducted for these cases using a recently updated chemical kinetic model and a 0-D engine model, where good qualitative and reasonable quantitative agreements in LTHR were obtained. Sensitivity analysis was also performed directly on the rates of LTHR, to understand the controlling chemical reactions of LTHR, providing further insights into the fuel-to-fuel differences. The results demonstrate the significant promoting effect of n-heptane on LTHR rates, while inhibiting effects were seen for ethanol and toluene. Also highlighted was the importance of H atom abstraction reactions from the chemistry of each fuel component, which could lead to contradictory fuel behavior depending on the locations of the H site of the abstraction reaction due to the different ensuing pathways for the primary fuel radicals.

LTHR, RON 90 binary blends, sensitivity analysis↗

Oxidation and pyrolysis of methyl propyl ether

The ignition, oxidation, and pyrolysis chemistry of methyl propyl ether (MPE) was probed experimentally at several different conditions, and a comprehensive chemical kinetic model was constructed to help understand the observations, with many of the key parameters computed using quantum chemistry and transition state theory. Experiments were carried out in a shock tube measuring time variation of CO concentrations, in a flow tube measuring product concentrations, and in a rapid compression machine (RCM) measuring ignition delay times. The detailed reaction mechanism was constructed using the Reaction Mechanism Generator software. Sensitivity and flux analyses were used to identify key rate and thermochemical parameters, which were then computed using quantum chemistry to improve the mechanism. Validation of the final model against the 1-20 bar 600-1500 K experimental data is presented with a discussion of the kinetics. The model is in excellent agreement with most of the shock tube and RCM data. Strong non-monotonic variation in conversion and product distribution is observed in the flow-tube experiments as the temperature is increased, and unusually strong pressure dependence and significant heat release during the compression stroke is observed in the RCM experiments. These observations are largely explained by a close competition between radical decomposition and addition to O2 at different sites in MPE; this causes small shifts in conditions to lead to big shifts in the dominant reaction pathways. The validated mechanism was used to study the chemistry occurring during ignition in a diesel engine, simulated using Ignition Quality Test (IQT) conditions. At the IQT conditions, where the MPE concentration is higher, bimolecular reactions of peroxy radicals are much more important than in the RCM.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of Gasoline Composition and Operating Parameters on the Response of End-Gas Autoignition to Nitric Oxide in a Lean-Burn DI-SI Engine

Lean operation of spark-ignition engines can lead to engine thermal efficiency gains and lower NOx emissions due to reduced combustion temperatures. Yet, lean operation could still face challenges in end-gas autoignition and knock generation due to higher intake pressures and trapped NO in the residual gas. Here, this study evaluates the impact of NO on end-gas autoignition for two gasoline fuels with similar octane rating but different composition: high cycloalkane fuel (HCA) and high olefin fuel (HO). Experiments were performed at stoichiometric and lean (λ = 2) conditions and at two engine speeds of 1400 rpm and 2000 rpm. Accompanying chemical kinetics simulations in CHEMKIN revealed that the mechanisms controlling the effect of NO on autoignition are similar λ = 2 and λ = 1, with NO + HO 2 = NO 2 + OH being the main pathway for enhancing reactivity by promoting low-temperature heat release (LTHR). The compositionally different fuels reacted differently to NO seeding and engine speed, and differences were augmented at λ = 2 compared to λ = 1 as the end-gas autoignition shifted to the low temperature regime. HO, which has inherent low temperature chemistry, was strongly impacted by engine speed at low NO seeding levels, with no noticeable peak of LTHR detected at 2000 rpm. On the other hand, LTHR of HCA was marginally affected by shortened residence time at higher engine speed as NO + HO 2 reaction was not greatly affected by shorter time scales, since HO 2 production was sustained even at 2000 rpm to support OH generation from NO + HO 2 . Contrary to HO, HCA exhibited greater sensitivity to NO seeding, as the increased OH production at higher NO concentrations offset the OH-quenching effect of cyclopentane, which accounts for 28.6% of HCA’s composition. Consequently, a sensitivity analysis revealed that fuels with weak inherent low-temperature chemistry, like HCA, are likely to be more sensitive to variations in NO concentration and charge temperature, whereas fuels with strong low temperature chemistry are more sensitive to variations in end-gas λ and intake pressure.

Energy↗

Implementation of new mixture rules has a substantial impact on combustion predictions for H 2 and NH 3

Complex-forming reactions comprise a substantial fraction of all important combustion reactions and are central to combustion behavior. Despite being often called “pressure-dependent” reactions, their rate constants depend on not only the pressure but also the composition. While modern combustion codes allow arbitrarily high accuracy in treating pressure dependence, recent work has consistently demonstrated dramatic failures of essentially all available treatments of mixture dependence. In situations where mixture dependence is treated at all, it is inevitably treated through specification of pressure-dependent rate constants for a set of pure bath gases, which are then combined to estimate the rate constant in a mixture via a “mixture rule.” While there had been a generally unquestioning confidence in these mixture rules, they had, in reality, been scarcely tested until the last decade, when comparisons against master equation calculations revealed order-of-magnitude errors for important pressure-dependent reactions. New mixture rules, based on the reduced pressure, have recently been proposed and shown to reproduce master equation calculations for broad classes of complex-forming reactions very accurately. Here, in this work, we present an implementation of one such new mixture rule (“LMR-R”) in Cantera and then use it to enable simulations that use new high-accuracy ab initio data for individual bath gases (for the first time, since codes previously could not accommodate the complex bath gas dependence). Demonstrations focus on combustion of H 2 and NH 3 , where (1) high-accuracy ab initio data are available and (2) the impact is expected to be large due to the high fractions of efficient colliders (e.g., H 2 O and NH 3 ) in the burned and unburned gases. Indeed, we find the impact of this treatment to be substantial and may explain previous modeling difficulties for these important carbon-free fuels, particularly for NH 3 , whose extraordinarily high third-body efficiency (~20) is often omitted from kinetic models.

Ammonia↗

Coupling of Lipid Peroxidation and Criegee Intermediate Mediated Autoxidation in the Heterogeneous Oxidation of Linoleic Acid Aerosols

Autoxidation is an established mechanism for the degradation of organic molecules, which is relevant in the atmosphere, combustion processes, the environment, and the rancidification of lipids (commonly known as lipid peroxidation). Autoxidation proceeds via radical chain reactions involving hydroxyl (•OH), peroxy (RO2•), and alkoxy radicals, which are also prominent oxidants in the atmosphere. Recent reports have provided evidence for an alternative autoxidation mechanism driven instead by Criegee intermediates (CIs), which are produced from the reaction of β-hydroxy peroxy radicals (β-OH-RO2•). This work evaluates the contributions of these two mechanisms in the •OH initiated heterogeneous oxidation of linoleic acid (LA) aerosols. Reaction kinetics and product distributions are monitored using a vacuum ultraviolet photoionization aerosol time-of-flight mass spectrometer. To explain the observed kinetics, a kinetic model is developed that incorporates both the conventional peroxidation and alternative CI-mediated autoxidation mechanisms. We observe that the CI-mediated autoxidation pathways enhance the heterogeneous autoxidation rate, while the peroxidation reactions, although present, contributes less to the overall oxidation rate. α-Acyloxyalkyl hydroperoxides (AAHPs) are identified as key indicators for bimolecular reactions of CI with LA, highlighting the role of LA as a CI scavenger. Moreover, the measured functionalized LA products with hydroxyl or carbonyl group(s), serve as markers for the peroxidation reactions. In summary, this work presents a quantitative framework to understand the coupled reaction network of •OH, RO2•, β-OH-RO2• radicals, and CI in driving heterogeneous autoxidation, which is crucial for understanding degradation mechanisms of organic molecules in the environment and atmosphere.

Criegee intermediate↗

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Conjugate heat transfer↗

Advanced sorbents for modular oxygen production for REMS gasifiers

Under a DOE funded effort, Thermosolv LLC has been developing a sorbent-based oxygen production technology for gasification and oxy-combustion applications. The sorbents utilize oxygen-storage properties of certain perovskites to (1) selectively adsorb oxygen at moderate temperature from compressed air and (2) release the adsorbed oxygen into a vacuum or a sweep gas such as CO2 and/or steam. Through cyclic operations of multiple sorbent beds such a process can be made continuous. Pressure drop across the sorbent bed and other similar design considerations dictate that the sorbent be in the form of a high surface-area-to-volume pellet and yet still possess sufficient crush strength and integrity. Process efficiency is determined by the sorption and desorption kinetics and sorbent capacity. Typical process involves cycle times in the order of 100 seconds, a time too short to fully utilize the full volume of the sorbent and the sorption capacity of the sorbent pellets. As a part of this project, Thermosolv LLC undertook the development of oxygen sorbents as a high-surface area supported sorbent on a light-weight inert support to (1) reduce the cost, (2) increase the productivity, and (3) reduce the overall weight of the reactor. Using a previously developed perovskite (LSCF 1991) with its well-characterized performance, composite sorbent pellets each consisting of an inert core coated by a thin layer of the functional perovskite material were produced. Several low-density, inexpensive inert solids supports in the 1/8”-1/4” size were selected to keep the pressure drop across the sorbent bed in the acceptable range. A number of commercial technologies including spray coating, dip coating, incipient impregnation of porous substrate followed by thermal annealing at various temperatures and duration were employed and tested to produce robust composite supported sorbent pellets. A sense of optimum coating thickness was developed by testing sorbent pellets of various increasing diameter pellets. LSCF-1991extrudates ranging in size from 1/32” to 3/16” were tested in a TGA and in a fixed-bed pressure swing test set-up for sorption/desorption cycles of interest. The data show that for the operational conditions of interest the coating thickness for a supported sorbent pellet needs to be approximately 1/32” (about 0.8 mm). Subsequent work thereby concentrated on developing composite pellets of various substrates, shapes and sizes with coating layer of about this size range. Candidate support materials were chosen based on cost, inertness, mechanical strength, thermal expansion and chemical stability with sorbent material, and in a size range to give an acceptable pressure drop in fixed-bad reactor configurations. Coating application methods used for application of the sorbent onto the support included precipitation, spray coating and dip coating from sorbent slurries and sol gels. For all coated supports where we could successfully apply a uniform coating of desired thickness with an acceptable handling performance in terms of exfoliation, TGA-based cyclic sorption/desorption testing was performed to determine cycling capacity of oxygen. Among nearly twenty different support materials tested, best adhesion performance was obtained from stainless saddle supports. In the bench-scale fixed-bad tests stainless steel saddle support composite pellets with a coating thickness in the 0.6 mm or so range, the composite sorbent pellet performance approached up to 95% of that of the 2 mm parent material pellets. In a parallel approach to reducing the cost of the sorbent, alternate sorbent formulations replacing/reducing the amount of cobalt in the LSCF family were also investigated. Successful formulations that could match the performance of LSCF 1991 were identified based on TGA cyclic tests as LSCF 1919 and LSF 1910. In the range of operational envelop of cycle times and other relevant process operating conditions, the reduced Co formulations showed comparable performance in the bench-scale fixed-bad cyclic operations. As a part of this project, we also attempted substitution of La and Sr with Ca and Ba, and Co with Mn, Ni and Cu with little success, but the overall project goal of reducing sorbent cost in terms of raw material and manufacturing expenses was successful.

01 COAL, LIGNITE, AND PEAT↗

An energy model of droplet impingement on an inclined wall under isothermal and non-isothermal environments

The observation of spray-wall interaction is of great importance to understand the dynamics that occur during fuel impingement onto the chamber wall or piston surfaces in internal combustion engines. It is found that the maximum spreading length of an impinged droplet can provide a quantitative estimation of heat transfer and energy transformation for spray-wall interaction. Furthermore, it influences the air-fuel mixing and hydrocarbon and particle emissions at combusting conditions. In this paper, an analytical model of different droplet-wall impingement conditions is developed in terms of βm (dimensionless maximum spreading length, the ratio of maximum spreading length to initial droplet diameter) to understand the detailed impinging dynamic process. These conditions are grouped as: a single diesel droplet impinging on the wall with different inclined angles (α); cold wall - heated droplet and heated wall - cold droplet impingement when inclined angle of the wall is 0°, respectively. The analytical model is built up based on the energy conservation that considers kinetic energy, gravitation energy, and surface energy before impingement, as well as viscous dissipation, gravitation energy, adhesion energy, deformation energy, and heat energy after impingement. The experimental work of diesel droplet impinging on an inclined wall is performed at a certain range of the Weber number (We of 33 to 420) with various inclined angles (α of 0° to 45°), while for inclined angle is 0°, droplet and wall temperature are varied from 25°C to 150°C to study the effects of the inclined angle and temperature on the temporal evolution of the post-impingement characteristics (i.e. droplet spreading length, dynamic contact angle). The analytical model is validated and evaluated at the aforementioned experimental operating points. The validated model can be employed to predict maximum spreading length of the droplet impinged on the wall. It is further utilized to determine the transition from capillary regime to kinetic regime, then to viscous regime at different inclined angle of the wall.

42 ENGINEERING↗

Plasma-assisted chemical-looping combustion: Mechanistic insights into low temperature methane oxidation with CuO

The low-temperature oxidation of CH 4 by CuO in a coaxial, fixed bed, double dielectric barrier discharge (DBD) reactor was investigated with time-dependent species measurements by an electron-ionization molecular beam mass spectrometer (EI-MBMS). In the experiment, 10% methane carried by noble gasses was flown at 50 sccm through 1 g CuO dispersed in quartz wool both under plasma and non-plasma conditions, while time-dependent gas-phase species profiles were collected. Plasma conditions were explored from 300 to 600 °C while the non-plasma conditions were set from 300 to 900 °C. Mechanistic insights into the oxidation of CH 4 by CuO with plasma discharge at lower temperatures (≤ 600 °C) were obtained by quantifying the fuel oxidation, intermediate species, and CO 2 production in comparison to the non-plasma conditions. Here, we observed significant enhancement of fuel oxidation from the plasma discharge between 400 and 500 °C. The CO 2 production at 500 °C with plasma was greater than that at 700 °C without plasma, reducing fuel oxidation temperature by 200+ °C. During tests, three distinct reaction stages were observed: a gas-phase transport limited stage, a surface reaction limited kinetic stage, and an oxygen ion diffusion limited stage. It was observed that plasma greatly improved the reactivity of the second stage at low temperature. In addition, no carbon deposits were observed on the resultant particles, even under the presence of plasma. M. species such as C 4 H 2 and C 6 H 6 not previously observed or predicted in CuO/CH 4 chemical looping were observed, with some species such as CH 3 OH only becoming detectable as total flowrate was increased from 50 to 1500sccm. A non-plasma reaction pathway for CH 4 based the observed species from the MBMS spectrum and previous predictions from reactive molecular dynamics simulations was created, providing a framework from which more complex plasma CuO mechanisms can be crafted from.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗

CRN Modeling of Ammonia RQL Combustion using a Partially-Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However low flammability and a propensity for high nitrogen oxide (NOx) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NOx emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel, while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NOx concentrations below that of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing non-adiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and post-flame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors’ past work, this arrangement has demonstrated the significant sensitivity of exit NOx to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and non-premixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NOx emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have shown to create conditions which are conducive to increased N2O formation in the lean stage. The results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia..

advanced gas turbines↗