Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic”

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 397 records · Page 22

Oxygen Diffusion and Reaction Kinetics in Continuous Fiber Ceramic Matrix Composites

Previous stressed oxidation tests of C/SiC composites at elevated temperatures (350 C to 1500 C) and sustained stresses (69 MPa and 172 MPa) have led to the development of a finite difference cracked matrix model. The times to failure in the samples suggest oxidation occurred in two kinetic regimes defined by the rate controlling mechanisms (i.e. diffusion controlled and reaction controlled kinetics). Microstructural analysis revealed preferential oxidation along as-fabricated, matrix microcracks and also suggested two regimes of oxidation kinetics dependent on the oxidation temperature. Based on experimental results, observation, and theory, a finite difference model was developed. The model simulates the diffusion of oxygen into a matrix crack bridged by carbon fibers. The model facilitates the study of the relative importance of temperature, the reaction rate constant, and the diffusion coefficient on the overall oxidation kinetics.

Halbig, Michael C.↗

Nonlinear Response of Layer Growth Dynamics in the Mixed Kinetics-Bulk-Transport Regime

In situ high-resolution interferometry on horizontal facets of the protein lysozyme reveal that the local growth rate R, vicinal slope p, and tangential (step) velocity v fluctuate by up to 80% of their average values. The time scale of these fluctuations, which occur under steady bulk transport conditions through the formation and decay of step bunches (macrosteps), is of the order of 10 min. The fluctuation amplitude of R increases with growth rate (supersaturation) and crystal size, while the amplitude of the v and p fluctuations changes relatively little. Based on a stability analysis for equidistant step trains in the mixed transport-interface-kinetics regime, we argue that the fluctuations originate from the coupling of bulk transport with nonlinear interface kinetics. Furthermore, step bunches moving across the interface in the direction of or opposite to the buoyancy-driven convective flow increase or decrease in height, respectively. This is in agreement with analytical treatments of the interaction of moving steps with solution flow. Major excursions in growth rate are associated with the formation of lattice defects (striations). We show that, in general, the system-dependent kinetic Peclet number, Pe(sub k) , i.e., the relative weight of bulk transport and interface kinetics in the control of the growth process, governs the step bunching dynamics. Since Pe(sub k) can be modified by either forced solution flow or suppression of buoyancy-driven convection under reduced gravity, this model provides a rationale for the choice of specific transport conditions to minimize the formation of compositional inhomogeneities under steady bulk nutrient crystallization conditions.

Vekilov, Peter G.↗

A Gas-Kinetic Scheme for Multimaterial Flows and Its Application in Chemical Reaction

This paper concerns the extension of the multicomponent gas-kinetic BGK-type scheme to multidimensional chemical reactive flow calculations. In the kinetic model, each component satisfies its individual gas-kinetic BGK equation and the equilibrium states of both components are coupled in space and time due to the momentum and energy exchange in the course of particle collisions. At the same time, according to the chemical reaction rule one component can be changed into another component with the release of energy, where the reactant and product could have different gamma. Many numerical test cases are included in this paper, which show the robustness and accuracy of kinetic approach in the description of multicomponent reactive flows.

Lian, Yongsheng↗

Measurement of the Turbulence Kinetic Energy Budget of a Turbulent Planar Wake Flow in Pressure Gradients

Turbulence kinetic energy (TKE) is a very important quantity for turbulence modeling and the budget of this quantity in its transport equation can provide insight into the flow physics. Turbulence kinetic energy budget measurements were conducted for a symmetric turbulent wake flow subjected to constant zero, favorable and adverse pressure gradients in year-three of research effort. The purpose of this study is to clarify the flow physics issues underlying the demonstrated influence of pressure gradient on wake development and provide experimental support for turbulence modeling. To ensure the reliability of these notoriously difficult measurements, the experimental procedure was carefully designed on the basis of an uncertainty analysis. Four different approaches, based on an isotropic turbulence assumption, a locally axisymmetric homogeneous turbulence assumption, a semi-isotropy assumption and a forced balance of the TKE equation, were applied for the estimate of the dissipation term. The pressure transport term is obtained from a forced balance of the turbulence kinetic energy equation. This report will present the results of the turbulence kinetic energy budget measurement and discuss their implication on the development of strained turbulent wakes.

Liu, Xiao-Feng↗

LSENS, A General Chemical Kinetics and Sensitivity Analysis Code for Homogeneous Gas-Phase Reactions: Code Description and Usage - Part 2

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part II of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part II describes the code, how to modify it, and its usage, including preparation of the problem data file required to execute LSENS. Code usage is illustrated by several example problems, which further explain preparation of the problem data file and show how to obtain desired accuracy in the computed results. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions. Part I (NASA RP-1328) derives the governing equations and describes the numerical solution procedures for the types of problems that can be solved by LSENS. Part III (NASA RP-1330) explains the kinetics and kinetics-plus-sensitivity-analysis problems supplied with LSENS and presents sample results.

Radhakrishnan, Krishnan↗

Kinetics of Organic Transformations Under Mild Aqueous Conditions: Implications for the Origin of Life and Its Metabolism

The rates of thermal transformation of organic molecules containing carbon, hydrogen, and oxygen were systematically examined in order to identify the kinetic constraints that governed origin-of-life organic chemistry under mild aqueous conditions. Arrhenius plots of the kinetic data were used to estimate the reaction half-life at 50 C, and to reveal the effect of functional groups on reactivity. This survey showed that hydrocarbons and organic substances containing a single oxygenated group were kinetically the most stable (i. e. acetate decarboxylation half-life was l0(exp 18) years at 50 C); whereas, organic substances containing two oxygenated groups in which one group was a beta-positioned carbonyl group were the most reactive (i. e. acetoacetate decarboxylation half-life was l0(exp-2) years at 50 C). Of all functional groups the beta-positioned carbonyl group (aldehyde or ketone) was the strongest activating group, giving rates of reaction that were up to 10(exp 24)-times faster than rates of similar molecules lacking the beta-carbonyl group. From this knowledge of organic reactivity and the inherent constraints of autocatalytic processes, we concluded that an origins-of-life process based on autocatalytic transformation of C,H,O-substrates was constrained to using the most reactive organic molecules that contain alpha- or beta-carbonyl groups, since small autocatalytic domains of plausible catalytic power that used less reactive substrates could not carry out chemical transformations fast enough to prevent catastrophic efflux (escape) of reaction intermediates. Knowledge of the kinetics of organic transformations is useful, not only in constraining the chemistry of the earliest autocatalytic process related to the origin of life, but also in establishing the relative reactivity of organic molecules on the early Earth and other planets that may or may not be related to the origin of life.

Weber, Arthur L.↗

Comparative kinetics and reciprocal inhibition of nitrate and nitrite uptake in roots of uninduced and induced barley (Hordeum vulgare L.) seedlings

Nitrate and NO2- transport by roots of 8-day-old uninduced and induced intact barley (Hordeum vulgare L. var CM 72) seedlings were compared to kinetic patterns, reciprocal inhibition of the transport systems, and the effect of the inhibitor, p-hydroxymercuribenzoate. Net uptake of NO3- and NO2- was measured by following the depletion of the ions from the uptake solutions. The roots of uninduced seedlings possessed a low concentration, saturable, low Km, possibly a constitutive uptake system, and a linear system for both NO3- and NO2-. The low Km system followed Michaelis-Menten kinetics and approached saturation between 40 and 100 micromolar, whereas the linear system was detected between 100 and 500 micromolar. In roots of induced seedlings, rates for both NO3- and NO2- uptake followed Michaelis-Menten kinetics and approached saturation at about 200 micromolar. In induced roots, two kinetically identifiable transport systems were resolved for each anion. At the lower substrate concentrations, less than 10 micromolar, the apparent low Kms of NO3- and NO2- uptake were 7 and 9 micromolar, respectively, and were similar to those of the low Km system in uninduced roots. At substrate concentrations between 10 and 200 micromolar, the apparent high Km values of NO3- uptake ranged from 34 to 36 micromolar and of NO2- uptake ranged from 41 to 49 micromolar. A linear system was also found in induced seedlings at concentrations above 500 micromolar. Double reciprocal plots indicated that NO3- and NO2- inhibited the uptake of each other competitively in both uninduced and induced seedlings; however, Ki values showed that NO3- was a more effective inhibitor than NO2-. Nitrate and NO2- transport by both the low and high Km systems were greatly inhibited by p-hydroxymercuribenzoate, whereas the linear system was only slightly inhibited.

NASA Program CELSS↗

Calcium kinetics with microgram stable isotope doses and saliva sampling

Studies of calcium kinetics require administration of tracer doses of calcium and subsequent repeated sampling of biological fluids. This study was designed to develop techniques that would allow estimation of calcium kinetics by using small (micrograms) doses of isotopes instead of the more common large (mg) doses to minimize tracer perturbation of the system and reduce cost, and to explore the use of saliva sampling as an alternative to blood sampling. Subjects received an oral dose (133 micrograms) of 43Ca and an i.v. dose (7.7 micrograms) of 46Ca. Isotopic enrichment in blood, urine, saliva and feces was well above thermal ionization mass spectrometry measurement precision up to 170 h after dosing. Fractional calcium absorptions determined from isotopic ratios in blood, urine and saliva were similar. Compartmental modeling revealed that kinetic parameters determined from serum or saliva data were similar, decreasing the necessity for blood samples. It is concluded from these results that calcium kinetics can be assessed with micrograms doses of stable isotopes, thereby reducing tracer costs and with saliva samples, thereby reducing the amount of blood needed.

Clinical Trial↗

Comparison of the Effects of Wave-Particle Interactions and the Kinetic Suprathermal Electron Population on the Acceleration of the Solar Wind

Kinetic effects due to wave-particle interactions and suprathermal electrons have been suggested in the literature as possible solar wind acceleration mechanisms. Ion cyclotron resonant heating, in particular, has been associated with some qualitative features observed in the solar wind. In terms of solar wind acceleration, however, it is interesting to compare the kinetic effects of suprathermal electrons with those due to the wave-particle interactions. The combined effects of the two acceleration mechanisms on the fast solar wind have been studied by Tam and Chang (1999a,b). In this study. we investigate the role of the suprathermal electron population in the acceleration of the solar wind. Our model follows the global kinetic evolution of the fast solar wind under the influence of ion cyclotron resonant heating, while taking into account Coulomb collisions, and the ambipolar electric field that is consistent with the particle distributions themselves. The kinetic effects due to the suprathermal electrons, which we define to be the tail of the electron distributions, can be included in the model as an option. By comparing the results with and without the inclusion of the suprathermal electron effects, we determine the relative importance of suprathermal electrons and wave-particle interactions in driving the solar wind. We find that although suprathermal electrons enhance the ambipolar electric potential in the solar wind considerably, their overall influence as an acceleration mechanism is relatively insignificant in a wave-driven solar wind.

Tam, S. W. Y.↗

Determining Kinetic Parameters for Isothermal Crystallization of Glasses

Non-isothermal crystallization techniques are frequently used to determine the kinetic parameters for crystallization in glasses. These techniques are experimentally simple and quick compared to the isothermal techniques. However, the analytical models used for non-isothermal data analysis, originally developed for describing isothermal transformation kinetics, are fundamentally flawed. The present paper describes a technique for determining the kinetic parameters for isothermal crystallization in glasses, which eliminates most of the common problems that generally make the studies of isothermal crystallization laborious and time consuming. In this technique, the volume fraction of glass that is crystallized as a function of time during an isothermal hold was determined using differential thermal analysis (DTA). The crystallization parameters for the lithium-disilicate (Li2O.2SiO2) model glass were first determined and compared to the same parameters determined by other techniques to establish the accuracy and usefulness of the present technique. This technique was then used to describe the crystallization kinetics of a complex Ca-Sr-Zn-silicate glass developed for sealing solid oxide fuel cells.

Ray, C. S.↗

Kinetic Temperature and Electron Density Measurement in an Inductively Coupled Plasma Torch using Degenerate Four-Wave Mixing

Laser wave mixing is presented as an effective technique for spatially resolved kinetic temperature measurements in an atmospheric-pressure radio-frequency inductively-coupled plasma. Measurements are performed in a 1 kW, 27 MHz RF plasma using a continuous-wave, tunable 811.5-nm diode laser to excite the 4s(sup 3)P2 approaches 4p(sup 3)D3 argon transition. Kinetic temperature measurements are made at five radial steps from the center of the torch and at four different torch heights. The kinetic temperature is determined by measuring simultaneously the line shape of the sub-Doppler backward phase-conjugate degenerate four-wave mixing and the Doppler-broadened forward-scattering degenerate four-wave mixing. The temperature measurements result in a range of 3,500 to 14,000 K+/-150 K. Electron densities measured range from 6.1 (+/-0.3) x 10(exp 15)/cu cm to 10.1 (+/-0.3) x 10(exp 15)/cu cm. The experimental spectra are analyzed using a perturbative treatment of the backward phase-conjugate and forward-geometry wave-mixing theory. Stark width is determined from the collisional broadening measured in the phase-conjugate geometry. Electron density measurements are made based on the Stark width. The kinetic temperature of the plasma was found to be more than halved by adding deionized water through the nebulizer.

Schafer, Julia↗

Evaluation of the Uncertainty in JP-7 Kinetics Models Applied to Scramjets

One of the challenges of designing and flying a scramjet-powered vehicle is the difficulty of preflight testing. Ground tests at realistic flight conditions introduce several sources of uncertainty to the flow that must be addressed. For example, the scales of the available facilities limit the size of vehicles that can be tested and so performance metrics for larger flight vehicles must be extrapolated from ground tests at smaller scales. To create the correct flow enthalpy for higher Mach number flows, most tunnels use a heater that introduces vitiates into the flow. At these conditions, the effects of the vitiates on the combustion process is of particular interest to the engine designer, where the ground test results must be extrapolated to flight conditions. In this paper, the uncertainty of the cracked JP-7 chemical kinetics used in the modeling of a hydrocarbon-fueled scramjet was investigated. The factors that were identified as contributing to uncertainty in the combustion process were the level of flow vitiation, the uncertainty of the kinetic model coefficients and the variation of flow properties between ground testing and flight. The method employed was to run simulations of small, unit problems and identify which variables were the principal sources of uncertainty for the mixture temperature. Then using this resulting subset of all the variables, the effects of the uncertainty caused by the chemical kinetics on a representative scramjet flow-path for both vitiated (ground) and nonvitiated (flight) flows were investigated. The simulations showed that only a few of the kinetic rate equations contribute to the uncertainty in the unit problem results, and when applied to the representative scramjet flowpath, the resulting temperature variability was on the order of 100 K. Both the vitiated and clean air results showed very similar levels of uncertainty, and the difference between the mean properties were generally within the range of uncertainty predicted.

Norris, A. T.↗

Constraints on Triple Oxygen Isotope Kinetics

Isotopic disequilibrium is not as well constrained as equilibrium, hindering interpretation of isotopic variations. Kinetic isotope effects, a subset of disequilibrium, are common in nature and have been assumed to be distinct from equilibrium and mass independent isotope effects based on underdeveloped criteria. Using basic physical principles, we provide needed mechanistic constraints on mass-dependent kinetic isotope effects for the triple oxygen isotope system. We find some kinetic isotope effects yield large isotopic variations, exceeding equilibrium, which could be mistaken for mass independent relationships. Meanwhile, other kinetic isotope effects are found to have triple oxygen isotope relationships that could be mistaken for equilibrium isotope effects. Comparison against prior case studies of thermal decomposition of calcite (CaCO3) and brucite (Mg(OH)2) further tests our results. Although oxygen is the focus here, our approach applies to any system with more than two isotopes.

Isotope Geochemistry↗

Zircon (U-Th)/He Impact Crater Thermochronometry and the Effects of Shock Microstructures on Helium Diffusion Kinetics

Accurate and precise age determination of impact cratering events remains challenging and often contentious; less than half of all known craters are regarded as accurately and precisely dated. Zircon (U-Th)/He (ZHe) dating of impactites can be employed to date medium to large impact structures as ZHe ages can be fully reset in minutes at T >1000°C, a plausible scenario in the central melt pool. In contrast, complete resetting of ZHe at 200-300°C, encountered near the crater margins or due to post-impact hydrothermal overprinting, may take >103-6 years. To test the reliability of ZHe impact dating, we have quantified the effects of shock-induced microstructures on helium diffusion kinetics in well-characterized variably shocked zircon. We investigated samples from two impact structures, the Chicxulub multi-ring crater and Ries complex crater, to compare diffusion kinetics from structures with different size, age, and hydrothermal system longevity. Shock microstructures were characterized by backscattered-electron imaging prior to determining the He diffusion kinetics by prograde and retrograde fractional-release experiments via light-bulb furnace with incremental step-heating (10°C) from 300°C to 600°C. Next, we examine the internal interconnectivity and sizes of the diffusion domains by EBSD. While we found that zircon with few shock microstructures exhibited no marked deviation from helium diffusion kinetics of undamaged zircon, zircon grains with planar microstructures and granular textures are characterized by a dramatic decrease in helium retentivity due to the reduction in the effective domain size and the introduction of interconnected fast diffusion pathways. A subset of grains were ZHe dated and showed that less deformed grains yielded a weighted mean age within error of the accepted impact ages, while the grains with planar microstructures or granular textures gave systematically younger ages. These new diffusion data and ZHe ages demonstrate that highly shocked grains are unsuitable for ZHe impact crater dating. Therefore, detailed characterization of impact-induced microstructures is critical for determining accurate ZHe impact ages and offers the possibility of investigating post-impact hydrothermal circulation.

Zircon↗

Constraints on triple oxygen isotope kinetics

Triple oxygen isotope measurements have been used as tracers for classifying meteorites, as indicators of atmospheric chemistry and as a proxy for primary productivity among other uses. These utilities rely on the mass dependency of most chemical reactions: δ17O is proportional to δ18O by roughly a factor of 0.5. However,in detail it is known that mass-dependency factorsvary. Constraints on mass dependency are well established for equilibrium isotope effects, but similar constraints for kinetic isotope effects are lacking.Here, we provide constraints for isotope effects generated by unidirectional chemical reactions (i.e. kinetic isotope effects, KIEs). A simplified Monte Carlo model is used to estimatethe range of kinetic isotope effects. It is found that the range of KIEs exceeds the range of equilibrium isotope effects such that some mass dependent kinetic isotope effects could be mistaken for being mass-independent. For the16O-17O-18O system,Δ′17Ovariations as large 0.8‰ can be generated from a single reaction step. Two case studies from published literature, on the thermal decomposition of calcite (CaCO3) and brucite (Mg(OH)2), are used to test the extreme limits of the modeled results. Isotope fractionationsgenerated fromdensity functional theory models of the surface reaction steps are shownto match the published data. The fractionation models are strictly mass dependent thus this result suggests thermal decomposition of calcite is also mass dependent. Although oxygen is the focus here, our approach applies to any system with more than two isotopes.

J. A. Hayles↗

Determining the Crosslinking and Degradation Reaction Kinetics in Photovoltaic Encapsulants Using Accelerated Aging: Preprint

Degradation of photovoltaic (PV) module encapsulant mechanical characteristics that lead to embrittlement and delamination remains a cause of failure in solar installations. A multiscale reliability model based on detailed molecular degradation reaction kinetics was previously published, connecting the encapsulant mechanical properties (elastic modulus, yield strength, and adhesion energy) to the degraded molecular structure and interfacial bond density to adjacent solar cell and glass substrates. The model, developed primarily for poly(ethylene-co-vinyl acetate) (EVA) encapsulants, remains to be experimentally validated. Determining the degradation and crosslinking kinetics of alternative encapsulants, such as polyolefin elastomer (POE) and EVA/POE/EVA composites (EPE), can further generalize the model. Importantly, the activation energy for crosslinking of fully cured PV encapsulant products from temperature or UV is presently unknown or unavailable. In this work, we subject EVA, POE, and EPE encapsulants to a set of accelerated aging conditions, varying the temperature, UV intensity, and relative humidity. We use DSC, FTIR, and Soxhlet extraction (gel content) to characterize the encapsulants' changing molecular structure. This allows for determination of the photochemical degradation and crosslinking kinetics of the encapsulants. Preliminary results show an increase in gel content (crosslinking) and a decrease in crystallinity of EVA, POE, and EPE encapsulants under hot-aerobic (90 deg C, 22% RH) and hot-anaerobic (90oC, sealed in N2 air) accelerated aging, even in the absence of UV and crosslinking initiating agents. For a least squares regression with an assumed first-order crosslinking kinetics, the crosslinking rate parameters were computed for the three encapsulants under hot-aerobic and hot-anaerobic aging conditions. FTIR analysis showed insignificant encapsulant degradation for hot-aerobic and hot-anaerobic aging, demonstrating the critical role of UV and moisture in accelerating degradation. Aging conditions with UV exposure and elevated humidity are presently in progress.

adhesion↗

Plasma-Coupled Flow Reactor Studies of Low-Temperature Plasma Assisted Kinetics of Methanol Blended with CO2

Ignition technologies based on low-temperature plasmas (LTP) have the potential to operate next generation engines at elevated pressures and increased dilution limits promoting higher efficiencies. From a practical standpoint, research on LTP igniters has shown to enhance combustion and ignition, improve flame stability, and extend the dilution limits of combustion. All of the aforementioned gains are complemented with higher ignition efficiencies. However, the biggest challenge with incorporating this technology into engines is the knowledge gap of how exactly plasma chemistry effects can enhance the basic combustion phenomena. This coupled with the lack of validated kinetic mechanisms for plasma-combustion chemistry is the biggest obstruction to recognize efficient ignition, especially for application relevant fuels and biofuels. In order to comprehensively evaluate the effects of LTP on an oxygenated fuel specific system, this present study examines the kinetics of methanol plasma-assisted pyrolysis and oxidation using a custom-built plasma flow reactor (PFR). Experimental regimes are also further extended to understand the effects of adding CO2 to the mixture and its consequence on reaction kinetics. The PFR is installed with a dielectric-barrier discharge (DBD) configuration to induce LTP into the fuel mixture. Non-equilibrium plasmas are generated by high-voltage pulses (nearing 20 KV) administered by a plasma pulser at high-pulse repetition rates (up to 10 kHz). In order to better understand and isolate plasma chemistry and its effect on neutral chemistry, the experiments were carried out by heavily diluting reactive mixtures in nitrogen at near isothermal conditions. This suppresses the effect of exothermic reactions on chemistry allowing stable intermediates and products to be detected and quantified using ex-situ GC/MS diagnostic methods. Experiments were carried out at 0.5 atm pressure and over a wide range of temperatures from 523 K to 1203 K. Experimental results depicted the enhancement in intermediates production as well as overall lower temperatures required for complete fuel consumption in the plasma specific cases as opposed to their pure thermal counterpart. Formation of oxygenated and nitrile compounds specific to plasma assisted pyrolysis cases illustrated the efficacy of LTP to introduce new reaction pathways accelerating fuel decomposition. Increase in reactivity at lower temperatures is sought to be an effect induced by plasma chemistry postulating the acceleration in intermediates production. The onset of thermal ignition in plasma assisted oxidation is seen 200 K earlier than thermal oxidation highlighting efficient fuel conversion to final byproducts. Enhanced collisional processes afforded by LTP is seen to perturb reaction pathways between oxygenated fuel radicals and N-atoms to alter overall chemical reactivity. New insights into kinetics are gained from this study highlighting governing plasma assisted combustion (PAC) pathways for oxygenated fuel reaction chemistry. The results from this study can be used to develop future mechanisms specific to plasma chemistry which can bridge the mechanistic knowledge gap for LTP ignition so that future engines can adopt LTPs in their design for efficient combustion.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗