At present, aero-engines face a major need to widen the ignition envelope. In order to provide a technical support to expand the high altitude ignition envelope of aero-engines, in this article we propose a novel igni...At present, aero-engines face a major need to widen the ignition envelope. In order to provide a technical support to expand the high altitude ignition envelope of aero-engines, in this article we propose a novel ignition technology, i.e., "precombustion plasma jet ignition technology". In this paper, we also design a pre-combustion plasma jet igniter. Its discharge characteristics, jet characteristics, and ignition effects are studied. The results show that increasing the equivalent ratio of jet gas can enhance the discharge stability and increase the duty cycle. At the same time, it can reduce working power and energy consumption. The increase of equivalent ratio in jet gas can enhance the length and ignition area of plasma jet.In the process of ignition, the pre-combustion plasma jet igniter has obvious advantages, suchn as shortening the ignition delay time and enlarging the ignition boundary. When the airflow velocity is 39.11 m/s and the inlet air temperature is80℃, compared with the spark igniter and the air plasma jet igniter, the pre-combustion plasma jet igniter has an ignition boundary that is expanded by 319.8% and 55.7% respectively.展开更多
In this paper,the methodology of the directed relation graph with error propagation and sensitivity analysis（DRGEPSA）,proposed by Niemeyer et al.（Combust Flame 157：1760-1770.2010）.and its differences to the origi...In this paper,the methodology of the directed relation graph with error propagation and sensitivity analysis（DRGEPSA）,proposed by Niemeyer et al.（Combust Flame 157：1760-1770.2010）.and its differences to the original directed relation graph method are described.Using DRGEPSA,the detailed mechanism of ethylene containing 71 species and 395 reaction steps is reduced to several skeletal mechanisms with different error thresholds.The 25-species and 131-step mechanism and the 24-species and115-step mechanism are found to be accurate for the predictions of ignition delay time and laminar flame speed.Although further reduction leads to a smaller skeletal mechanism with 19 species and 68 steps,it is no longer able to represent the correct reaction processes.With the DRGEPSA method,a detailed mechanism for n-dodecane considering low-temperature chemistry and containing 2115 species and8157 steps is reduced to a much smaller mechanism with249 species and 910 steps while retaining good accuracy.If considering only high-temperature（higher than 1000 K）applications,the detailed mechanism can be simplified to even smaller mechanisms with 65 species and 340 steps or48 species and 220 steps.Furthermore,a detailed mechanism for a kerosene surrogate having 207 species and 1592 steps is reduced with various error thresholds and the results show that the 72-species and 429-step mechanism and the66-species and 392-step mechanism are capable of predicting correct combustion properties compared to those of the detailed mechanism.It is well recognized that kinetic mechanisms can be effectively used in computations only after they are reduced to an acceptable size level for computation capacity and at the same time retaining accuracy.Thus,the skeletal mechanisms generated from the present work are expected to be useful for the application of kinetic mechanisms of hydrocarbons to numerical simulations of turbulent or supersonic combustion.展开更多
Ignition delay times for n-decane/O2/Ar mixtures were measured behind reflected shock waves using endwall pressure and CH＊ emission measurements in a heated shock tube. The initial postshock conditions cover pressure...Ignition delay times for n-decane/O2/Ar mixtures were measured behind reflected shock waves using endwall pressure and CH＊ emission measurements in a heated shock tube. The initial postshock conditions cover pressures of 0.09-0.26 MPa, temperatures of 1 227-1 536 K, and oxygen mole fractions of 3.9%-20.7% with an equivalence ratio of 1.0. The correlation formula of ignition delay dependence on pressure, temperature, and oxygen mole fraction was obtained. The current data are in good agreement with available low-pressure experimental data, and they are then compared with the prediction of a kinetic mechanism. The current measurements extend the kinetic modeling targets for the n-decane combustion at low pressures.展开更多
A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentration...A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentrations less than 10-6 mol·L-1 do not activate any chemical reaction. So, within each step, the time steplength Δt of the algorithm is determined from the fastest reaction rate maxR by the formula Δt = 10-6mol·L-1/max R. All the reversible elementary reactions occur simultaneously;however, by a simple book-keeping technique, the updating of species concentrations, within each step of the algorithm, is performed within each elementary reaction separately. The above proposed simple algorithm for Chemical Kinetics is applied to a simple model for hydrogen combustion with only five reversible elementary reactions (Initiation, Propagation, First and Second Branching, Termination by wall destruction) with six species (H2, O2, H, O, HO, H2O). These five reversible reactions are recommended in the literature as the most significant elementary reactions of hydrogen combustion [1] [2]. Based on the proposed here simple algorithm for Chemical Kinetics, applied to the global mechanism of proposed five reversible elementary reactions for hydrogen combustion, a simple and short computer program has been developed with only about 120 Fortran instructions. By this proposed program, the following are obtained: 1) The total species concentration of hydrogen combustion, starting from the sum of initial reactants concentrations [H2] + [O2], gradually diminishes, due to termination reaction by wall destruction, and tends to the final concentration of the product [H2O], that is to the 2/3 of its initial value, in accordance to the established overall stoichiometric reaction of hydrogen combustion 2H2 + O2 → 2H2O. 2) Time-histories for concentrations of main species H2, O2, H, H2O of hydrogen combustion, in explosion and equilibrium regions, obtained by the proposed program, are compared to corresponding ones obtained by a展开更多
As it is very difficult to release boron energy completely, kinetic mechanism of boron is not clear, which leads to the lack of theoretical guidance for studying how to accelerate boron combustion. A new semi-empirica...As it is very difficult to release boron energy completely, kinetic mechanism of boron is not clear, which leads to the lack of theoretical guidance for studying how to accelerate boron combustion. A new semi-empirical boron combustion model is built on the King combustion model, which contains a chemical reaction path; two new methods of plasma-assisted boron combustion based on kinetic and thermal effects respectively are built on the ZDPLASKIN zero-dimensional plasma model. A plasma-supporting system is constructed based on the planar flame, discharge characteristics and the spectral characteristics of plasma and boron combustion are analyzed. The results show that discharge power does not change the sorts of excited-particles, but which can change the concentration of excited-particles. Under this experimental condition,plasma kinetic effect will become the strongest at the discharge power of 40 W; when the discharge power is less than 40 W,plasma mainly has kinetic effect, otherwise plasma has thermal effect. Numerical simulation result based on plasma kinetic effect is consistent with the experimental result at the discharge power of 40 W, and boron ignition delay time is shortened by 53.8% at the discharge power of 40 W, which indicates that plasma accelerates boron combustion has reaction kinetic paths, while the ability to accelerate boron combustion based on thermal effect is limited.展开更多
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51776223 and 91741112)
文摘At present, aero-engines face a major need to widen the ignition envelope. In order to provide a technical support to expand the high altitude ignition envelope of aero-engines, in this article we propose a novel ignition technology, i.e., "precombustion plasma jet ignition technology". In this paper, we also design a pre-combustion plasma jet igniter. Its discharge characteristics, jet characteristics, and ignition effects are studied. The results show that increasing the equivalent ratio of jet gas can enhance the discharge stability and increase the duty cycle. At the same time, it can reduce working power and energy consumption. The increase of equivalent ratio in jet gas can enhance the length and ignition area of plasma jet.In the process of ignition, the pre-combustion plasma jet igniter has obvious advantages, suchn as shortening the ignition delay time and enlarging the ignition boundary. When the airflow velocity is 39.11 m/s and the inlet air temperature is80℃, compared with the spark igniter and the air plasma jet igniter, the pre-combustion plasma jet igniter has an ignition boundary that is expanded by 319.8% and 55.7% respectively.
基金supported by the National Natural Science Foundation of China(Grant 11172309)supported by the China's Programme of Introducing Talents of Discipline to Universities-111 Project under Grant No.B08009the Thousand Talents Program
文摘In this paper,the methodology of the directed relation graph with error propagation and sensitivity analysis（DRGEPSA）,proposed by Niemeyer et al.（Combust Flame 157：1760-1770.2010）.and its differences to the original directed relation graph method are described.Using DRGEPSA,the detailed mechanism of ethylene containing 71 species and 395 reaction steps is reduced to several skeletal mechanisms with different error thresholds.The 25-species and 131-step mechanism and the 24-species and115-step mechanism are found to be accurate for the predictions of ignition delay time and laminar flame speed.Although further reduction leads to a smaller skeletal mechanism with 19 species and 68 steps,it is no longer able to represent the correct reaction processes.With the DRGEPSA method,a detailed mechanism for n-dodecane considering low-temperature chemistry and containing 2115 species and8157 steps is reduced to a much smaller mechanism with249 species and 910 steps while retaining good accuracy.If considering only high-temperature（higher than 1000 K）applications,the detailed mechanism can be simplified to even smaller mechanisms with 65 species and 340 steps or48 species and 220 steps.Furthermore,a detailed mechanism for a kerosene surrogate having 207 species and 1592 steps is reduced with various error thresholds and the results show that the 72-species and 429-step mechanism and the66-species and 392-step mechanism are capable of predicting correct combustion properties compared to those of the detailed mechanism.It is well recognized that kinetic mechanisms can be effectively used in computations only after they are reduced to an acceptable size level for computation capacity and at the same time retaining accuracy.Thus,the skeletal mechanisms generated from the present work are expected to be useful for the application of kinetic mechanisms of hydrocarbons to numerical simulations of turbulent or supersonic combustion.
基金supported by the National Natural Science Foundation of China (91016002)
文摘Ignition delay times for n-decane/O2/Ar mixtures were measured behind reflected shock waves using endwall pressure and CH＊ emission measurements in a heated shock tube. The initial postshock conditions cover pressures of 0.09-0.26 MPa, temperatures of 1 227-1 536 K, and oxygen mole fractions of 3.9%-20.7% with an equivalence ratio of 1.0. The correlation formula of ignition delay dependence on pressure, temperature, and oxygen mole fraction was obtained. The current data are in good agreement with available low-pressure experimental data, and they are then compared with the prediction of a kinetic mechanism. The current measurements extend the kinetic modeling targets for the n-decane combustion at low pressures.
文摘A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentrations less than 10-6 mol·L-1 do not activate any chemical reaction. So, within each step, the time steplength Δt of the algorithm is determined from the fastest reaction rate maxR by the formula Δt = 10-6mol·L-1/max R. All the reversible elementary reactions occur simultaneously;however, by a simple book-keeping technique, the updating of species concentrations, within each step of the algorithm, is performed within each elementary reaction separately. The above proposed simple algorithm for Chemical Kinetics is applied to a simple model for hydrogen combustion with only five reversible elementary reactions (Initiation, Propagation, First and Second Branching, Termination by wall destruction) with six species (H2, O2, H, O, HO, H2O). These five reversible reactions are recommended in the literature as the most significant elementary reactions of hydrogen combustion [1] [2]. Based on the proposed here simple algorithm for Chemical Kinetics, applied to the global mechanism of proposed five reversible elementary reactions for hydrogen combustion, a simple and short computer program has been developed with only about 120 Fortran instructions. By this proposed program, the following are obtained: 1) The total species concentration of hydrogen combustion, starting from the sum of initial reactants concentrations [H2] + [O2], gradually diminishes, due to termination reaction by wall destruction, and tends to the final concentration of the product [H2O], that is to the 2/3 of its initial value, in accordance to the established overall stoichiometric reaction of hydrogen combustion 2H2 + O2 → 2H2O. 2) Time-histories for concentrations of main species H2, O2, H, H2O of hydrogen combustion, in explosion and equilibrium regions, obtained by the proposed program, are compared to corresponding ones obtained by a
基金Project supported by the National Natural Science Foundation of China(Grant No.11372356)
文摘As it is very difficult to release boron energy completely, kinetic mechanism of boron is not clear, which leads to the lack of theoretical guidance for studying how to accelerate boron combustion. A new semi-empirical boron combustion model is built on the King combustion model, which contains a chemical reaction path; two new methods of plasma-assisted boron combustion based on kinetic and thermal effects respectively are built on the ZDPLASKIN zero-dimensional plasma model. A plasma-supporting system is constructed based on the planar flame, discharge characteristics and the spectral characteristics of plasma and boron combustion are analyzed. The results show that discharge power does not change the sorts of excited-particles, but which can change the concentration of excited-particles. Under this experimental condition,plasma kinetic effect will become the strongest at the discharge power of 40 W; when the discharge power is less than 40 W,plasma mainly has kinetic effect, otherwise plasma has thermal effect. Numerical simulation result based on plasma kinetic effect is consistent with the experimental result at the discharge power of 40 W, and boron ignition delay time is shortened by 53.8% at the discharge power of 40 W, which indicates that plasma accelerates boron combustion has reaction kinetic paths, while the ability to accelerate boron combustion based on thermal effect is limited.
