CN104036075B - A kind of burner structure optimization method based on response phase method - Google Patents
A kind of burner structure optimization method based on response phase method Download PDFInfo
- Publication number
- CN104036075B CN104036075B CN201410246922.0A CN201410246922A CN104036075B CN 104036075 B CN104036075 B CN 104036075B CN 201410246922 A CN201410246922 A CN 201410246922A CN 104036075 B CN104036075 B CN 104036075B
- Authority
- CN
- China
- Prior art keywords
- point
- response
- ccd
- factor
- design
- 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.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of burner structure optimization method based on response phase method, concretely comprise the following steps：(1) determine that burner needs the structural factor considered, and transcoding, coding transform is carried out to it；(2) testing program is determined using Central Composite design CCD；(3) complete to test according to CCD testing programs；(4) response surface model is established, analyzes result of the test；(5) ridge analysis is carried out, determines that optimum structural parameter combines.The present invention designs CCD methods by Central Composite and determines testing program, and pass through response surface software analysis result of the test, it can not only determine that each factor influences on the linear term of burner combustion performance and quadratic term influences, the cross term of each factor, which can also be determined, to be influenceed, then the structural parameters for determining optimization by ridge analysis combine, so as to realize low nitrogen burning effect.
Description
【Technical field】
The present invention relates to burner structure optimisation technique field, is a kind of burner based on response phase method specifically
Structural optimization method.
【Background technology】
Environmental problem relation national security, more obtains the attention of government and enterprise level.On July 30th, 2013, middle stone
Change and announce to implement " clear water blue sky " Environmental Protection Plan, and one of nitrogen oxides pollutant that to be exactly it pay close attention to.Nitrogen oxides is big
Amount is present in the exiting flue gas of tubular heater, environmentally safe composition grave danger：Nitrogen oxides be photochemical fog and
One major reason of acid rain, the nitrogen oxides in vehicle exhaust react what is formed with hydrocarbon through ultraviolet irradiation
Toxic smog, referred to as photochemical fog.In addition, nitrogen oxides and the nitric acid and nitrous acid of the water reaction generation in air are acid rain
Composition.
The direction that low nitrogen burning is Petrochemical Enterprises all the time and heating furnace producer is pursued.At present, domestic manufacturer is to combustion
The structure of burner has carried out certain transformation, is realized by using combustion technologies such as combustion gas classifications to nitrogen oxidation in combustion process
The control of thing.Although these combustion technologies reduce NOx discharge to a certain extent, it can not still reach content and be less than
50ppm environmental requirement.How on the basis of existing burner structure using Rational structure size come the combustion of optimizing combustor
Performance is burnt, is be worthy of consideration the problem of.
Find that heating furnace producer and domestic scholars are to burner knot by the retrieval to prior art and Field Research
The optimization of structure mainly uses single_factor method.Li Yuanqing, Lu Xiaofeng exist《Nanjing University of Technology's journal (natural science edition)》2010,32
(5) burning has been investigated in the article " the CFD researchs of petrochemical industry tubular heater gas burner " delivered on using single_factor method
Influence of two factors of device spray angle and nozzle diameter to combustion process.Wang Juan, filoplume et al. exists《Oil Refining Technologies and engineering》
Also distinguish in the article " influence that burner structure is distributed to gas flame shape and in-furnace temperature " delivered on 2007,37 (8)
The influence of injection diameter, jet nozzle angle, eddy flow angle has been investigated using single_factor method.But burner structure is complex, Dan Yin
Element experiment can not determine the cross influence between each structural parameters, therefore it is not typically that object function (exports NO to optimize numerical valuex
Content) extremal region.
【The content of the invention】
It is an object of the invention to overcome the deficiencies of the prior art and provide a kind of burner structure based on response phase method is excellent
Change method.
The purpose of the present invention is achieved through the following technical solutions：
A kind of burner structure optimization method based on response phase method, it is concretely comprised the following steps：
(1) determine that burner needs the structural factor considered, and transcoding, coding transform is carried out to it；
Described structural factor is fire trough, and spray gun, stabilizer, burn road shape, diameter, Burners Positions, nozzle angle,
One or several kinds of combinations in nozzle diameter；
Described transcoding, coding transform is：The excursion of each structural factor is different in response surface design, some from change
The scope difference of amount is extremely greatly different；For the convenience being uniformly processed, the value of design parameter is made into transcoding, coding transform, it is established that structure because
The one-to-one relationship of plain level value and " coding "；Coding is to solve the trouble that dimension difference is brought to design and analysis；
(2) testing program is determined using Central Composite design CCD；
CCD=central composite design are writing a Chinese character in simplified form for Central Composite design；
Described Central Composite design CCD refers to：Central Composite design is a kind of classical second-order design, and testing site includes
Cube three kinds of point, central point and axial point；
For specially with three factors central composite designs, cube point is distributed in the cubical summit that the length of side is 1, center
Point be located at cube center, axial point is also known as initial point, and asterisk point is distributed in through in the axial direction of central point, except a coordinate be+
Outside α or-α, remaining coordinate is all 0；
(3) complete to test according to CCD testing programs；
Described refers to according to the completion experiment of CCD testing programs：According to CCD testing programs, using hot test or meter
Fluid operator mechanics CFD software completes the design experiment of all cubes of points, central point and axial point, passes through flue gas analyzer or meter
Fluid operator mechanics CFD post processings obtain the NO of testing sitexNumerical value.
(4) response surface model is established, analyzes result of the test；
Described establishes response surface model, analyzes result of the test, refers to：Test Point Data is inputted into response surface software, example
Such as, Statistica, Design-Expert or SAS；According to software analysis result, multinomial coefficient is extracted, obtains each factor
Polynomial relation between response；Output response surface graphics simultaneously judges influence situation of each factor to object function；Root
According to representative value analytical table, judge whether to need to carry out ridge analysis；
(5) ridge analysis is carried out, determines that optimum structural parameter combines.
Described ridge analysis, refers to：Using original design central point as the centre of sphere, each coordinate is constantly expanded since the centre of sphere
Greatly, optimal response value is found out in the footprint of the intersection point formation of hypersphere and response surface that R is radius, wherein, radius R is not
It can exceed that trial stretch.
Described intersection point also known as ridge point；
Compared with prior art, the positive effect of the present invention is：
The present invention is better than traditional single factor test optimization method, for example single factor test optimization method can not so determine that spray orifice is straight
Cross influence between footpath, spray angle, eddy flow angle.
The present invention designs CCD methods by Central Composite and determines testing program, and by response surface software analysis result of the test,
It can not only determine that each factor influences on the linear term of burner combustion performance and quadratic term influences, each factor can also be determined
Cross term influences, and the structural parameters for then determining optimization by ridge analysis combine, so as to realize low nitrogen burning effect.
【Brief description of the drawings】
The factor CCD testing sites of Fig. 1 tri- are distributed；
Fig. 2 response surface graphics；
(a) H, α influence on response；
(b) H, β influence on response；
(c) H, γ influence on response；
(d) α, β influence on response；
(e) α, γ influence on response；
(f) β, γ influence on response；
Each structure NO cloud atlas contrasts of Fig. 3；
(a) structure is optimized；
(b) 29th structures；
(c) 16th structures.
【Embodiment】
A kind of embodiment of the burner structure optimization method based on response phase method of the present invention presented below.
Embodiment 1
A kind of burner structure optimization method based on response phase method, it is concretely comprised the following steps：
(1) determine that burner needs the structural factor considered, and transcoding, coding transform is carried out to it；
Described structural factor is fire trough, and spray gun, stabilizer, burn road shape, diameter, Burners Positions, nozzle angle,
One or several kinds of combinations in nozzle diameter；
Described transcoding, coding transform is：The excursion of each structural factor is different in response surface design, some from change
The scope difference of amount is extremely greatly different；For the convenience being uniformly processed, the value of design parameter is made into transcoding, coding transform, it is established that structure because
The one-to-one relationship of plain level value and " coding "；Coding is to solve the trouble that dimension difference is brought to design and analysis；
(2) testing program is determined using Central Composite design CCD；
CCD=central composite design are writing a Chinese character in simplified form for Central Composite design；
Described Central Composite design CCD refers to：Central Composite design is a kind of classical second-order design, and testing site includes
Cube three kinds of point, central point and axial point；
For specially with three factors central composite designs, cube point is distributed in the cubical summit that the length of side is 1, center
Point be located at cube center, axial point is also known as initial point, and asterisk point is distributed in through in the axial direction of central point, except a coordinate be+
Outside α or-α, remaining coordinate is all 0；
(3) complete to test according to CCD testing programs；
Described refers to according to the completion experiment of CCD testing programs：According to CCD testing programs, using hot test or meter
Fluid operator mechanics CFD software completes the design experiment of all cubes of points, central point and axial point, passes through flue gas analyzer or meter
Fluid operator mechanics CFD post processings obtain the NO of testing sitexNumerical value.
(4) response surface model is established, analyzes result of the test；
Described establishes response surface model, analyzes result of the test, refers to：Test Point Data is inputted into response surface software, example
Such as, Statistica, Design-Expert or SAS；According to software analysis result, multinomial coefficient is extracted, obtains each factor
Polynomial relation between response；Output response surface graphics simultaneously judges influence situation of each factor to object function；Root
According to representative value analytical table, judge whether to need to carry out ridge analysis；
(5) ridge analysis is carried out, determines that optimum structural parameter combines.
Described ridge analysis, refers to：Using original design central point as the centre of sphere, each coordinate is constantly expanded since the centre of sphere
Greatly, optimal response value is found out in the footprint of the intersection point formation of hypersphere and response surface that R is radius, wherein, radius R is not
It can exceed that trial stretch.
It is specific as follows：
1st, determine to need the structural factor considered
The present embodiment mainly optimizes the burner structure of certain reboiler furnace, and the burner is combustion gas gradual combustor.By dividing
Analyse burner structure find, main and auxiliary spray gun relative altitude H, main jet rifle nozzle angle α, auxiliary spray tip angle β, γ etc. four because
Element is the principal element for influenceing burner combustion performance.The problems such as in view of calculating the time, other factors air channel length etc.
Do not account for for the time being.
2nd, transcoding, coding transform
The present embodiment illustrates the rule of transcoding, coding transform by taking H transcoding, coding transform as an example.In view of burner and the chi of burner hearth
Very little, H excursion is (Hmin=130, Hmax180), then H transcoding, coding transform is as follows：
Half length in section
Central point
Then H coding formula is
, can be with by section (H by this transcoding, coding transformmin, Hmax) be converted into (- 1,1).
Similarly, the transcoding, coding transform result of other factors can be obtained, it is seen that table 1.
Each factor transcoding, coding transform result of table 1
3rd, the design of CCD Central Composites and CFD simulations
At present, CCD be with most commonly used experimental design, testing site include three kinds of cube point, central point and axial point,
Accompanying drawing 1 show the distribution of three factor CCD testing sites.The present embodiment is designed using CCD Central Composites, and test sequence is random
's.After completing the design of CCD Central Composites, according to design table (schedule), whole simulated tests are completed using CFD software Fluent.CCD
Central Composite designs and result of the test is as shown in table 2.
Table 2CCD Central Composites design and result of the test
4th, response surface model is established, analyzes result of the test
Response surface analysis is carried out to result of the test using Statistica softwares, obtains returning using NO discharge capacitys as response
Return equation：
Y=9.39H+4.01H2-1.21α+0.13α2+4.68β+0.58β2+1.71γ+0.87γ2-1.75Hα+3.55Hβ
+1.22Hγ+1.65αβ+0.2αγ-1.27βγ
Analysis result can be shown in Table 3, from the table, during using NO discharge capacitys as response, model P ＜ 0.0001, show this
Quadratic equation model is notable, and fitting effect is good, and reliability is high.As P ＜ 0.05, that is, represent that this index is notable, P ＜
When 0.001, that is, show that this index is extremely notable.It can be seen that H linear terms, H quadratic terms and β linear terms are extremely significant model
, H, β product term are significant model terms, lose and intend item not significantly (P ＞ 0.1), and regression equation degree of fitting and confidence level are higher, examination
It is smaller to test error, illustrates that the model can be used for instructing burner structure optimization design.
The response phase method analysis result of table 3
The figure of Responds Surface Methodology is a three dimensions of each factors composition corresponding to specific response, Ke Yizhi
See ground and reflect the influence of each factor to response, can be analyzed from the response surface analysis figure of experiment gained between each factor
Interaction.Because response surface figure is graphics, the function containing two influence factors can only be expressed, therefore each response surface is three-dimensional
Figure is required for fixing two of which influence factor positioned at the position for being encoded to 0.The response surface graphics of each factor interaction
See accompanying drawing 2.
5th, ridge analysis
From representative value analytical table 4 as can be seen that the characteristic value of four factors has just has negative, it is saddle to show this Quadratic response
Face, without unique optimum value, it is therefore desirable to make ridge analysis, further determine that optimal response value.
The representative value analytical table of table 4
Ridge analysis be using the original design central point centre of sphere (this experiment go the centre of sphere as<- 1,0,0,0>), it is radius in R
Optimal response value is found out in the footprint that the intersection point of hypersphere and response surface (i.e. ridge point) is formed.The result of ridge analysis is pair
Each coordinate constantly expands since the centre of sphere, radius R no more than trial stretch, choose herein R 0.1,0.2 ..., 0.8 model
Enclose interior calculating ridge point.The present embodiment carries out ridge analysis using SAS softwares.
The ridge analysis of table 5
As shown in Table 5, as coding radius R increase, response NO are gradually reduced.As R ＞ 0.8, standard error compared with
Greatly, do not consider.The R=0.8 in this test level scope, minimum response value, now corresponding non-coding physical dimension is H=
140th, α=89, β=29, γ=92.Checking simulation is done using Fluent to calculate, obtain NO realities according to above-mentioned optimum structure parameter
Border analog result is 30.7ppm, basically identical with predicted value.Accompanying drawing 3 for optimization structure and CCD design in table the 16th group of structure and
The NO cloud atlas contrast of 29th group of structure, it can be found that the NO numerical value of optimization structure significantly reduces, significantly realizes NO emission reduction effects,
Project expection is reached.
Described above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
Member, without departing from the inventive concept of the premise, can also make some improvements and modifications, these improvements and modifications also should be regarded as
In protection scope of the present invention.
Claims (6)
1. a kind of burner structure optimization method based on response phase method, it is characterised in that concretely comprise the following steps：
(1) determine that burner needs the structural factor considered, and transcoding, coding transform is carried out to it；
(2) testing program is determined using Central Composite design CCD；
(3) complete to test according to CCD testing programs；
(4) response surface model is established, analyzes result of the test；
(5) ridge analysis is carried out, determines that optimum structural parameter combines；
Described ridge analysis, refers to：Using original design central point as the centre of sphere, each coordinate is constantly expanded since the centre of sphere,
Optimal response value is found out in the footprint of the intersection point formation of hypersphere and response surface that R is radius, wherein, radius R can not
More than trial stretch；
Described structural factor is specially：Main and auxiliary spray gun relative altitude H, main jet rifle nozzle angle α, auxiliary spray tip angle β,
The factors of γ tetra-；
Described establishes response surface model according to regression equation：
Y=9.39H+4.01H2-1.21α+0.13α2+4.68β+0.58β2+1.71γ+0.87γ2-1.75Hα+3.55Hβ+
α γ -1.27 β the γ of+1.65 α β of 1.22H γ+0.2 are established.
2. a kind of burner structure optimization method based on response phase method as claimed in claim 1, it is characterised in that described
The step of (1) in, described transcoding, coding transform is：The excursion of each structural factor is different in response surface design, some
The scope difference of independent variable is extremely greatly different；For the convenience being uniformly processed, the value of design parameter is made into transcoding, coding transform, it is established that knot
Structure factor level value and the one-to-one relationship of " coding ".
3. a kind of burner structure optimization method based on response phase method as claimed in claim 1, it is characterised in that described
The step of (2) in, described Central Composite design CCD refer to：Central Composite design is a kind of classical second-order design, testing site
Including a cube point, three kinds of central point and axial point.
4. a kind of burner structure optimization method based on response phase method as claimed in claim 1, it is characterised in that described
The step of (2) in, described Central Composite design CCD be specially with three factors central composite designs for, a cube point is distributed in
The length of side is 1 cubical summit, and central point is located at cube center, and axial point is also known as initial point, asterisk point, is distributed in
In the axial direction of heart point, in addition to a coordinate is+α or-α, remaining coordinate is all 0.
5. a kind of burner structure optimization method based on response phase method as claimed in claim 1, it is characterised in that described
The step of (3) in, it is described according to CCD testing programs complete experiment refer to：According to CCD testing programs, using hot test or
Person's computation fluid dynamics software completes the design experiment of all cubes of points, central point and axial point, by flue gas analyzer or
The post processing of person's computation fluid dynamics obtains the NO of testing sitexNumerical value.
6. a kind of burner structure optimization method based on response phase method as claimed in claim 1, it is characterised in that described
The step of (4) in, it is described to establish response surface model, analyze result of the test, refer to：Test Point Data input response is thin-skinned
Part, according to software analysis result, multinomial coefficient is extracted, obtains the polynomial relation between each factor and response；Output rings
Answer face graphics and judge influence situation of each factor to object function；According to representative value analytical table, judge whether to need to carry out
Ridge analysis.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410246922.0A CN104036075B (en) | 2014-06-05 | 2014-06-05 | A kind of burner structure optimization method based on response phase method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410246922.0A CN104036075B (en) | 2014-06-05 | 2014-06-05 | A kind of burner structure optimization method based on response phase method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104036075A CN104036075A (en) | 2014-09-10 |
| CN104036075B true CN104036075B (en) | 2018-02-09 |
Family
ID=51466846
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410246922.0A Expired - Fee Related CN104036075B (en) | 2014-06-05 | 2014-06-05 | A kind of burner structure optimization method based on response phase method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104036075B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104732084B (en) * | 2015-03-20 | 2017-07-11 | 上海交通大学 | The model building method of rapid digestion content of starch is reduced based on response phase method |
| CN107220400A (en) * | 2017-04-05 | 2017-09-29 | 中国石油化工股份有限公司 | A kind of high-temperature air burning organization optimization method based on response phase method |
| CN113569448B (en) * | 2021-07-13 | 2023-12-19 | 沈阳航空航天大学 | Optimization method for strain gauge sensitive gate structural parameters based on response surface method |
| CN114636169B (en) * | 2022-02-21 | 2023-01-20 | 北京航空航天大学 | Flame stabilizer perforation parameter determination method and device and radial flame stabilizer |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102360402A (en) * | 2011-10-26 | 2012-02-22 | 中冶南方工程技术有限公司 | Optimized design method for hydraulic cylinder supporting seat of cover carriage based on response surface |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103366043B (en) * | 2013-05-31 | 2016-01-27 | 天津大学 | A kind of bed piece supported at three point Optimization Design based on response surface model |
| CN103324798B (en) * | 2013-06-25 | 2016-01-27 | 福州大学 | Based on the stochastic response of interval response surface model |
-
2014
- 2014-06-05 CN CN201410246922.0A patent/CN104036075B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102360402A (en) * | 2011-10-26 | 2012-02-22 | 中冶南方工程技术有限公司 | Optimized design method for hydraulic cylinder supporting seat of cover carriage based on response surface |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104036075A (en) | 2014-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104036075B (en) | A kind of burner structure optimization method based on response phase method | |
| CN103235842B (en) | The acquisition methods of quadrangle tangential circle coal-fired boiler combustion characteristic and system | |
| CN103234219A (en) | Circumference air quantity adjustment method and system for pulverized coal boiler after fire coal kind changing | |
| Ribert et al. | Hybrid transported-tabulated strategy to downsize detailed chemistry for numerical simulation of premixed flames | |
| Tang et al. | Investigation of combustion structure and flame stabilization in an axisymmetric scramjet | |
| CN103216846B (en) | Smoke concentration adjustment method and smoke concentration adjustment system when fire coal variety of pulverized coal boiler is changed | |
| Bazooyar et al. | Numerical analysis of nitrogen oxides in turbulent lifted H2/N2 cabra jet flame issuing into a vitiated coflow | |
| Mazaheri et al. | Numerical optimization of laboratory combustor geometry for NO suppression | |
| CN103244953B (en) | Burner output adjustment method and system after changing of fire coal types of pulverized coal boiler | |
| Jamalabadi | Effect of fuel inject angle on non-premixed combustion of air/methane mixtures in vertical cylinder | |
| CN107220400A (en) | A kind of high-temperature air burning organization optimization method based on response phase method | |
| CN105498530A (en) | Method for instructing ammonia spray adjustment through introducing tracer gases in SCR system simulation | |
| Rahman et al. | Hydrogen Co-Firing Characteristics in a Single Swirl Burner: A Numerical Analysis | |
| Andreini et al. | Large Eddy Simulation Based Computational Fluid Dynamics Investigation of the Ignition Process in Lean Spray Burner | |
| Shrotriya et al. | Large eddy simulation of partially premixed flames with inhomogeneous inlets based on the DTF combustion model | |
| Aleksandrov et al. | Optimal design of a combustion chamber of gas turbine engine by a Combustion chamber 1D-2D computer program | |
| Rochette | Modeling and simulation of two-phase flow turbulent combustion in aeronautical engines | |
| Attili et al. | Mixing and entrainment of burned products in high Karlovitz number premixed jet flames | |
| Lan et al. | Effects of initial conditions on MILD combustion for diesel in a constant pressure combustion bomb | |
| CN105930597A (en) | Porous fuel injector and design method therefor | |
| CN117057238B (en) | Combustor stable combustion blunt body design method based on physical information operator network model | |
| Yakush et al. | Large Eddy Simulation of Hydrocarbon Fireballs | |
| Al-Khalidy | Determination of Flame Plume Characteristics utilising CFD and Experimental Approaches | |
| Fomenko et al. | Numerical study of the methane-air combustion in the direct-flow burner of the boiler 300 MW TGMP-314 boiler | |
| Liu | The combined effect of water mist system and longitudinal ventilation on fire and smoke behavior in tunnel fires |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180209 Termination date: 20190605 |