CN101900349A - Method for using fuse turbulence combustion information to simplify combustion reaction mechanism - Google Patents
Method for using fuse turbulence combustion information to simplify combustion reaction mechanism Download PDFInfo
- Publication number
- CN101900349A CN101900349A CN 201010189485 CN201010189485A CN101900349A CN 101900349 A CN101900349 A CN 101900349A CN 201010189485 CN201010189485 CN 201010189485 CN 201010189485 A CN201010189485 A CN 201010189485A CN 101900349 A CN101900349 A CN 101900349A
- Authority
- CN
- China
- Prior art keywords
- combustion
- information
- turbulence
- reaction mechanism
- odt
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
Abstract
The invention relates to a method for using fuse turbulence combustion information to simplify combustion reaction mechanism, wherein the obtained mechanism is applied to the analysis or computer numerical simulation of actual combustion process; and the method solves the present situation that the building of the reaction mechanism is all based on laminar flow combustion. The invention is characterized that the obtaining method of the turbulence flow combustion impulse information is proposed. When the simplify mechanism is built, multiple-point sampling is required in combustion field space, and simultaneously the combination data of speed and scalar is required; such large quantity of information is difficultly obtained just relaying on experiment, thus the construction of one-dimensional turbulence ODT numerical platform is put forwards; and the experiment measure data is used as the auxiliary and the verification of model. The reaction mechanism simplification application method of fuse turbulence impulse information is put forwards. The combustion key control parameter-mixed fraction and scalar dissipation rate type variable are selected as generalized components when combustion; and the simplified analysis method is put forwards, thereby obtaining the combustion simplification mechanism of vertex event information. In theory, the invention solves the two difficult problems of complete turbulence combustion information source and mechanism simplification method how to fuse turbulence impulse information; in technology, since the actual combustion nearly is turbulence combustion, the invention greatly expands the application range of combustion reaction mechanism and settles the basis for combustion analysis.
Description
Technical field
The present invention relates to merge turbulence combustion information, realize the method that fuel reaction mechanism is simplified, analysis or the Computer Numerical Simulation that simplification mechanism is applied to the Actual combustion process that obtain.
Background technology
Burning is one of general phenomenon the most in nature and the related industries device.Burning can adopt chemical equation to describe, and can be expressed as methyl hydride combustion: CH
4+ O
2=CO
2+ H
2O, such expression way is called one-step reaction.Combustion chemistry and theoretical chemistry or the like association area is through many decades research in the world, at present quite meticulous for the description of fuel reaction mechanism, almost considered all possible means of chemical reaction, such reaction mechanism is called the fuel detailed mechanism, and detailed mechanism comprises hundreds of extremely thousands of chemical equations, relates to tens of extremely hundreds of materials.For the hydrocarbon substance burning, many mechanism have been developed in the world, for instance: at low carbonizable substance (C
1-C
3Alkane or unsaturated hydrocarbons), developed GRI-Mech mechanism, Konnov mechanism, Leeds mechanism etc. at present in the world; U.S. Lawrence Livermore National Laboratory has developed from low-carbon (LC) C
1Up to the high-carbon substance C
16The fuel detailed mechanism.
Burning is the process of a complexity, in order to understand and to control combustion process, adopts computer at present both at home and abroad in a large number, and combustion process is analyzed or calculated.Theory analysis or numerical computations etc. need the chemical reaction mechanism of fuel, and promptly the reaction mechanism of fuel is one of of paramount importance data in combustion analysis and the prediction (or claiming file).Clearly, adopt one-step reaction too simple, it can not reflect the Actual combustion process, is bound to produce materials such as CO as methyl hydride combustion, but does not comprise CO etc. in the aforesaid single step chemical equation; But then, chemism is too loaded down with trivial details in detail.The detailed chemism of describing combustion process comprises too many reaction and component, and amount of calculation is too big when making Computer Analysis, and even more serious is that the reaction time yardstick of each chemical equation is widely different in the detailed mechanism, is 10
0Second is to 10
-10Second, the Mathematical Modeling that the greatest differences of time scale causes describing combustion process is strong rigidity, and finding the solution of stiff equation is very difficult, has more aggravated the consumption of computer resource simultaneously.For instance: Raman and Pitsch once adopted the methane mechanism of one 175 step chemical reaction, and a retention flame has been carried out analytical calculation, and reach more than 200 hours the running time on a 8-processor 600-MHZ computer.As seen, in order to carry out the analytical calculation of combustion process, need the mechanism of suitable scale, but then, the scale of reaction mechanism is increasing, as, 0.5 edition mechanism of Konnov comprises 1207 reactions, relates to 127 kinds of materials; The normal heptane mechanism of U.S. Lawrence Livermore National Laboratory development comprises 2540 reactions, relates to 557 kinds of materials; Methane that Gupta etc. work out and natural gas mechanism comprise 6874 elementary reactions, relate to 345 kinds of materials; Unite States Standard Technical Board NIST has compiled about 12000 pieces of documents about the hot phase chemical reaction, has set up the reaction mechanism database, and this database has comprised the parameter more than 38000 reactions.All detailed mechanisms of Gong Buing comprise redundant reaction and component in the world, and its reason mainly contains two: the one, and constructing of detailed mechanism tends to comprise additional component and reaction, even their importance can't be determined; The 2nd, the test specification of detailed mechanism is very wide, and the parameter area that mechanism is actual when using is comparatively narrow.When using detailed mechanism to simulate, except amount of calculation greatly, from numerical algorithm, the equation group rigidity has greatly retrained the use of numerical method, the problem that makes is difficult to find the solution.Simplify mechanism between one-step reaction and detailed mechanism, it can describe the of paramount importance chemistry route of combustion process simultaneously, and therefore, development detailed mechanism method for simplifying has very great practical meaning in engineering.
Research about the mechanism method for simplifying is continuing always in the world, and the method that has developed mainly contains sensitivity analysis method and time scale analytic approach.The sensitivity analysis method is the analytical method that a class early grows up, more present famous combustion analysis simulate commercial softwares, all has the sensitivity analysis module as CHEMKIN, RUN1DL and COSILAB etc., its use is comparatively simple, but the result of sensitivity analysis can not directly provide decoupling zero information, need carry out a large amount of post processing work, research and application person's experience are had higher requirements.The time scale analytic approach mainly comprises the intrinsic low-dimensional stream shape ILDM of Maas and Pope development and the calculating singular perturbation CSP of development such as Goussis and Lam etc.; the time scale analytic approach is carried out the time scale analysis to each chemical reaction or material in the reaction mechanism; thereby will react or material is divided into fast and slow two classes, thereby by to fast response or material is removed or modelling realizes that fuel mechanism simplifies.
As everyone knows, ignition quality in reality and the engineering all is the turbulent combustion process, turbulent flow has significant effects for burning, but, the domestic and international at present method of simplifying all about mechanism is all based on laminar combustion, and the simplification mechanism of Gou Jianing is not suitable for actual turbulent combustion analysis and calculating certainly like this.The energy, environment and safety etc. are closely related with turbulent combustion, the fuel simplification reaction mechanism of consideration turbulent flow pulsation is structured in theory with on using and is not achieved, thereby, can realize that the fuel mechanism method for simplifying that merges the turbulent flow pulsation information has huge demand at aspects such as scientific research, industrial production and security controls.
Summary of the invention
The invention provides a kind of method that the turbulent flow pulsation information realizes that fuel reaction mechanism is simplified that merges, it is characterized in that: (one)
Turbulent combustion is proposed The pulsation information acquisition methodsBased on the jet flame experimental provision, speed of uniting and scalar (temperature and concentration of component) are measured.The jet flame generator adopts reaction stream nozzle (coflow nozzle) structure in the same way, the temperature and the flow of control gas during air feed.Laser-Doppler, three dimensional particles kinetic analyzer LDA/APV and DPIV system are adopted in tachometric survey.Thermal imaging system and contactless ultrasonic wave CT technology are adopted in temperature survey.The analysis of combustion components adopts flame ion probe and online gas chromatograph-mass spectrometer GC-MS, in conjunction with Fourier infrared spectrum FTIR constituent analysis.Even more important, make up the numerical value platform of the detailed reaction mechanism of coupling fuel of whirlpool incident (turbulent flow pulsation) fusion, and select one dimension turbulent flow ODT model is suitable, its reason is: 1) construct when simplifying mechanism, need in coupling numerical value of detailed mechanism or experimental result space, sample, to obtain the information under the different burning conditions, the adaptability of the simplification mechanism that sampled point is many more, construct is strong more, it is difficult depending merely on experiment acquisition a large amount of information like this, therefore, need to make up corresponding numerical value platform; 2) ODT can catch the fine dimension amount, calculation cost is far below direct Numerical DNS simultaneously, ODT is a self-contained model, model does not rely on empirical parameter, it is for the complete instantaneous state according to the flow field of finding the solution of flow field, thereby it can catch typical turbulent flow phenomenon such as energy level string.Because ODT is an one-dimensional model, even calculating yardstick, it is divided into the Kolmogorov yardstick, its calculation cost still can bear.For the computing time of ODT, the applicant had once carried out the ODT simulation to the jet flame of a coupling methane complicated chemical (containing 279 reaction and 49 components), the Reynolds number is 10000, and the computing grid number is 1201, taking turns simulated time and be about 19 hours on the double-core high-performance computer.(2)
Propose to merge the reaction mechanism method for simplifying of turbulent flow pulsation informationControl turbulent combustion process has some critical quantity, as mixing mark and scalar dissipative shock wave etc., should consider these parameters when mechanism is simplified, with reflection turbulent combustion feature.Mix mark reflected fuel and oxidant mixability, controlled the fired state (is extent of reaction for premixed flame) of local, as when higher turbulent flow Reynolds counts Ret, heat and component molecular diffusion difference are generally second order to the influence of burning, when describing the key property of flame structure, can ignore this influence, at this moment mix the major control factor of mark for burning.We are once to ethene C
2H
4The stable state diffusion flame is analyzed and is calculated, and finds CO
2Can explain such flame fully with the mixing mark, and when increasing the unstable state factor, whole phase deviation takes place in result of calculation, and increases the scalar dissipative shock wave again, then can eliminate phase error.Simplification process when merging for more specifically describing whirlpool incident (turbulent flow pulsation), with the sensitivity analysis is example, the influence degree of the variation of certain component (or certain reaction) to other component only investigated in present sensitivity analysis, thereby determine redundant component or reaction, the application's method will be mixed mark and be included analysis in as a kind of what is called " broad sense component ", analogically, will be more comprehensive and appropriate as the analytical method that increases scalar dissipative shock wave class variable again.Owing to incorporated the whirlpool event information in analysis, the mechanism of building predicts the outcome approaching more actual combustion process.
The math equation that merges the reaction mechanism method for simplifying of turbulent flow pulsation information is described,
If Y
F, Y
OBe respectively fuel and oxidant mass fraction, β is oxidant-fuel equivalence ratio (being that the unit oxidant can consume β part fuel), introduces conversion X=Y
F-Y
O/ β mixes mark f and is defined as,
Wherein, X
1, X
2Represent X value respectively corresponding to fuel-side and oxygen side.Scalar dissipative shock wave χ is defined as,
Wherein, D is a diffusion coefficient, x
iExpression i direction x coordinate, f is for mixing mark.
Wherein, J
KiBe regular Jacobian matrix, the variation of another component j that perturbation caused of expression combustion components k concentration, its physical significance is the influence size of component k to component j.Of particular note: the component j here not only represents the component of burning and relating to, and also comprises mixing mark and scalar dissipative shock wave, and just as previously mentioned, we will mix mark and scalar dissipative shock wave conduct " broad sense component ".
Set a less real number ε (as ε=0.0001), if (J
Kj)
2Greater than ε, then think to should remain in k component outbalance in the former detailed reaction mechanism; If (J
Kj)
2Less than ε, think that then the k component is inessential, k component and its corresponding chemical reaction can be removed from former detailed reaction mechanism, thereby be realized detailed reaction mechanism simplification.
Compare with present existing method, technology, the present invention has two main outstanding features:
1 method that proposes has general research and using value, shows two aspects: choosing of (1) control turbulent flame critical quantity is pointed.Existing method for simplifying is all at laminar flame, and the simplification reaction mechanism of acquisition is applicable to laminar flame, is not suitable for turbulent flame, because the present invention has considered the critical quantity of control turbulent flame, institute obtains simplifies reaction mechanism applicable to turbulent flame.Can choose the critical quantity of controlling turbulent flame according to analyzing and simulate needs:, then only choose the mixing mark and get final product if only be concerned about the macroscopic properties of turbulent flame; If also be concerned about and stop working, restrike and the pollutant emission characteristic, then can choose and mix mark and scalar dissipative shock wave; If wish to analyze some special fire behaviors, also can add the extent of reaction critical quantity.Need to prove, choose after the critical quantity of controlling turbulent flame and the simplification reaction mechanism that obtains, go for the analysis and the simulation of all turbulent flames of the type; (2) choosing of simplified way has flexibility.The reaction mechanism method for simplifying of above-mentioned fusion turbulent flow pulsation information has adopted the Jacobian Matrix Technology, clearly, adopts component generating rate ratio to replace above-mentioned J
Kj, maybe will to mix mark and scalar dissipative shock wave also passable as modes such as stream shape dimensions, this makes method of the present invention have directive significance.
2. solved a turbulence combustion information acquisition methods difficult problem.Reaction mechanism is simplified needs complete combustion field data and information, it is to be difficult to realize that lot of data like this and information are depended merely on experiment, do not relate to the chemistry of fuel mechanism that merges turbulent flow information at present both at home and abroad and simplify, we think that a very big difficulty is to obtain detailed reliable turbulent combustion field information.Because the amount of calculation of direct Numerical DNS is too big, the DNS of the detailed chemism of coupling fuel can't realize at present, because the limiting it and also do not have a promotional value of hardware.The applicant realizes one dimension turbulent flow ODT method, by theory analysis and with a large amount of comparisons of international standard experimental data, confirm to adopt ODT to obtain turbulence combustion information and can solve this difficult problem.
Description of drawings
Fig. 1 is a schematic diagram of the present invention
The specific embodiment
Below embodiments of the present invention are further described:
The first step: based on the jet flame burner, the temperature and the flow of control gas during air feed, experiment measuring speed and scalar data.These experimental datas and not exclusively as turbulence combustion information source, and mainly as the Data Source of model and method validation.
Second step: make up one dimension turbulent flow ODT numerical value platform.Change the dimensionless number of quarterization combustion systems such as Reynolds number and Lewis number, calculate the turbulence reaction that obtains the coupling detailed mechanism and flow complete data.
The 3rd step: sampling.Select some characteristic points as the data analysis point, than burning, rich oil, lean burn etc., design an integral algorithm again, travel through whole zone as high temperature, low temperature, equivalent.
The 4th step: the simplification mechanism that merges the whirlpool incident is constructed.Choose the critical quantity (as mixing mark, scalar dissipative shock wave categorical variable) of control turbulent combustion process, they as " broad sense component ", are analyzed their variation.Select what type or what " broad sense component ", depend on the purposes of simplifying mechanism, for high Reynolds number burning, it is often enough to choose the mixing mark; When research is stopped working and restriked, can choose and mix mark and two parameters of scalar dissipative shock wave.The perturbation that generally needs when simplification mechanism is constructed to calculate a certain component can cause the great variation of mixing mark, adopt the matrix analysis may fall through, reason is to mix the linear combination that mark generally is some component of system, in addition, the increment of scalar dissipative shock wave categorical variable certainly obtains by its definition and transport equation, but the scalar dissipative shock wave model of development is mostly at perfect condition, as height at present
Number and even Lewis number etc.For solving these difficulties, these turbulent flow critical quantity can obtain from ODT result.
The 5th step: the fail-safe analysis and the checking of constructing reaction mechanism.The coupling reaction mechanism of constructing is analyzed and is calculated turbulent combustion, and its result and the result of the experiment and the detailed mechanism that is coupled are compared, and carries out fail-safe analysis and checking.
Claims (2)
1. one kind merges the method that turbulence combustion information realizes that fuel reaction mechanism is simplified.Solved a turbulence combustion information acquisition methods difficult problem, proposed to obtain the method for whirlpool event information by one dimension turbulent flow ODT numerical result.Owing to construct when simplifying mechanism, for making the adaptive surface of simplifying reaction mechanism wider, need sample at the multiple spot in the combustion field space, the associating data that need speed, temperature and concentration of component simultaneously, here the component that relates to reaches tens of kinds until hundreds of, and clearly, it is impossible depending merely on experiment acquisition a large amount of and complete data like this, we propose to make up the numerical value platform of ODT, and measured data of experiment is as the usefulness of auxiliary and modelling verification.Can the be coupled detailed chemical reaction mechanism of fuel of ODT, can catch the fine dimension amount, calculation cost is far below direct Numerical DNS simultaneously, ODT is a self-contained model, model does not rely on empirical parameter, it is for the complete instantaneous state according to the flow field of finding the solution of flow field, thereby it can catch typical turbulent flow phenomenon such as energy level string.Because ODT is an one-dimensional model, even calculating yardstick, it is divided into the Kolmogorov yardstick, its calculation cost still can bear.
2. the reaction mechanism that merges the turbulent flow pulsation information is simplified implementation method.Choose turbulent combustion and close key control parameter-mixing mark and scalar dissipative shock wave categorical variable, with they " broad sense components " as when burning, propose to simplify analytical method, i.e. the Jacobian analytic approach of concentration of component perturbation and will mix mark and mode that the scalar dissipative shock wave is tieed up as stream shape.The fuel simplification mechanism that obtains has like this comprised the whirlpool event information, has expanded the range of application of fuel reaction mechanism greatly.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010189485 CN101900349A (en) | 2010-05-27 | 2010-05-27 | Method for using fuse turbulence combustion information to simplify combustion reaction mechanism |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010189485 CN101900349A (en) | 2010-05-27 | 2010-05-27 | Method for using fuse turbulence combustion information to simplify combustion reaction mechanism |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN101900349A true CN101900349A (en) | 2010-12-01 |
Family
ID=43226120
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201010189485 Pending CN101900349A (en) | 2010-05-27 | 2010-05-27 | Method for using fuse turbulence combustion information to simplify combustion reaction mechanism |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101900349A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104461688A (en) * | 2014-11-27 | 2015-03-25 | 中国矿业大学 | Method for simplifying skeleton of detailed combustion chemical reaction mechanism |
| CN109767815A (en) * | 2018-12-25 | 2019-05-17 | 哈尔滨工业大学 | The method simplified based on the combustion reaction mechanisms under rate uncertainty |
-
2010
- 2010-05-27 CN CN 201010189485 patent/CN101900349A/en active Pending
Non-Patent Citations (2)
| Title |
|---|
| 《北京航空航天大学学报》 20020630 姜雪丽等 湍流燃烧反应机理的简化 291-293 1-2 第28卷, 第3期 2 * |
| 《燃烧科学与技术》 20061031 蒋勇等 基于一维全尺度湍流模型的氢气射流扩散火焰结构数值模拟 401-407 1-2 第12卷, 第5期 2 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104461688A (en) * | 2014-11-27 | 2015-03-25 | 中国矿业大学 | Method for simplifying skeleton of detailed combustion chemical reaction mechanism |
| CN104461688B (en) * | 2014-11-27 | 2017-09-01 | 中国矿业大学 | A kind of method that detailed chem-ical reaction mechanism skeleton simplifies |
| CN109767815A (en) * | 2018-12-25 | 2019-05-17 | 哈尔滨工业大学 | The method simplified based on the combustion reaction mechanisms under rate uncertainty |
| CN109767815B (en) * | 2018-12-25 | 2023-01-24 | 哈尔滨工业大学 | Method for simplifying combustion reaction mechanism based on rate uncertainty |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cellek et al. | Investigations on performance and emission characteristics of an industrial low swirl burner while burning natural gas, methane, hydrogen-enriched natural gas and hydrogen as fuels | |
| Hargis et al. | Methane ignition in a shock tube with high levels of CO2 dilution: consideration of the reflected-shock bifurcation | |
| Faravelli et al. | A new procedure for predicting NOx emissions from furnaces | |
| Yin et al. | Chemistry and radiation in oxy-fuel combustion: a computational fluid dynamics modeling study | |
| Cuoci et al. | Numerical modeling of NO x formation in turbulent flames using a kinetic post-processing technique | |
| Mei et al. | Dimensions of CH4-jet flame in hot O2/CO2 coflow | |
| Wang et al. | Combustion regimes of a jet diffusion flame in hot co-flow | |
| Funke et al. | Numerical and experimental evaluation of a dual-fuel dry-low-NOx micromix combustor for industrial gas turbine applications | |
| Le Cong et al. | Oxidation of natural gas, natural gas/syngas mixtures, and effect of burnt gas recirculation: Experimental and detailed kinetic modeling | |
| Fürst et al. | Optimization of chemical kinetics for methane and biomass pyrolysis products in moderate or intense low-Oxygen dilution combustion | |
| Monaghan et al. | Detailed emissions prediction for a turbulent swirling nonpremixed flame | |
| Funke et al. | Comparison of numerical combustion models for hydrogen and hydrogen-rich syngas applied for dry-low-nox-micromix-combustion | |
| Medwell et al. | Laminar flame calculations for analyzing trends in autoignitive jet flames in a hot and vitiated coflow | |
| Chen et al. | Numerical study of a cyclonic combustor under moderate or intense low-oxygen dilution conditions using non-adiabatic tabulated chemistry | |
| Castiñeira et al. | Computational fluid dynamics for simulation of wind-tunnel experiments on flare combustion systems | |
| Yadav et al. | Large eddy simulation of coal-ammonia flames with varied ammonia injection locations using a flamelet-based approach | |
| Resende et al. | Impacts of dilution on hydrogen combustion characteristics and nox emissions | |
| CN101900349A (en) | Method for using fuse turbulence combustion information to simplify combustion reaction mechanism | |
| Wang et al. | Auto-and forced-ignition temperatures of diffusion flames obtained through the steady RANS modeling | |
| Wang et al. | A modified two-dimensional numerical method for prediction of outer wall temperature distribution of rectangular micro-combustors | |
| Darbandi et al. | Numerical study of inlet turbulators effect on the thermal characteristics of a jet propulsion-fueled combustor and its hazardous pollutants emission | |
| Zhang et al. | Numerical and experimental assessment of a novel Multinozzle Burner with CO2 diluent to improve the emissions from a swirling flame in a combustion chamber | |
| Zhang et al. | Combustion mode and mixing characteristics of a reacting jet in crossflow | |
| Sayah et al. | Computational fluid dynamics modeling of a self-recuperative burner and development of a simplified equivalent radiative model | |
| Sadiki et al. | Thermodynamically consistent modelling of gas turbine combustion sprays |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20101201 |