CN104036075A - Combustor structure optimization method based on response surface method - Google Patents
Combustor structure optimization method based on response surface method Download PDFInfo
- Publication number
- CN104036075A CN104036075A CN201410246922.0A CN201410246922A CN104036075A CN 104036075 A CN104036075 A CN 104036075A CN 201410246922 A CN201410246922 A CN 201410246922A CN 104036075 A CN104036075 A CN 104036075A
- Authority
- CN
- China
- Prior art keywords
- response surface
- point
- ccd
- design
- structure optimization
- 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.)
- Granted
Links
Abstract
The invention relates to a combustor structure optimization method based on a response surface method. The combustor structure optimization method comprises the following specific steps: (1) determining structural factors, required to be considered, of a combustor, and performing coding conversion on the structural factors; (2) determining a test scheme by adopting CCD (Central Composite Design) method; (3) completing a test according to the CCD-based test scheme; (4) establishing a response surface model, and analyzing a test result; (5) performing ridge analysis, and determining an optimized structure parameter combination. According to the combustor structure optimization method, the test scheme is determined by using the CCD method, and the test result is analyzed by virtue of response surface software, so that the linear term influence and the quadratic term influence on the combustion performance of the combustor caused by all the factors can be determined, the cross term influence of all the factors can also be determined, the optimized structure parameter combination is further determined by the ridge analysis, and a low-nitrogen combustion effect is achieved.
Description
[technical field]
The present invention relates to burner structure optimisation technique field, specifically, is a kind of burner structure optimization method based on response surface method.
[background technology]
Environmental problem is related to national security, more obtains the attention of Government and enterprise aspect.On July 30th, 2013, China Petrochemical Industry announces to implement " clear water blue sky " Environmental Protection Plan, and one of oxides of nitrogen pollutant that to be exactly it pay close attention to.Oxides of nitrogen is present in the outlet flue gas of tubular heater in a large number, environmentally safe formation grave danger: oxides of nitrogen is a major reason of photo-chemical smog and acid rain, oxides of nitrogen in vehicle exhaust and hydrocarbon irradiate the toxic smog that reacts and form through ultraviolet ray, be called photo-chemical smog.In addition, oxides of nitrogen react with airborne water generate nitric acid and nitrous acid be the composition of acid rain.
Low nitrogen burning is the direction that Petrochemical Enterprises and heating furnace producer are pursued all the time.At present, domestic manufacturer has carried out certain transformation to the structure of burner, by adopting the combustion technologies such as combustion gas classification to realize the control to oxides of nitrogen in combustion process.Although these combustion technologies have reduced the discharge of NOx to a certain extent, still can not reach the environmental requirement of content lower than 50ppm.How on the basis of existing burner structure, adopting Rational structure size to carry out the burning performance of Optimizing Combustion device, is a considerable problem.
Find through the retrieval to prior art and Field Research, heating furnace producer and domestic scholars mainly adopt single factor method to the optimization of burner structure.Li Yuanqing, in the article of delivering " the CFD research of petrochemical industry tubular heater gas burner " of Lu Xiaofeng on " Nanjing University of Technology's journal (natural science edition) " 2010,32 (5), utilize single factor method to investigate injection angle of burners and the impact of two factors of nozzle diameter on combustion process.Wang Juan, in article that the people such as filoplume deliver in " oil Refining Technologies and engineering " 2007,37 (8) " burner structure on gas flame shape and stove in the impact of Temperature Distribution ", also adopt respectively single factor method to investigate the impact of injection diameter, jet nozzle angle, eddy flow angle.But burner structure is comparatively complicated, and single factor experiment cannot be determined the cross influence between each structural parameters, therefore optimizing numerical value is not generally that objective function (exports NO
xcontent) extremal region.
[summary of the invention]
The object of the invention is to overcome the deficiencies in the prior art, a kind of burner structure optimization method based on response surface method is provided.
The object of the invention is to be achieved through the following technical solutions:
A burner structure optimization method based on response surface method, its concrete steps are:
(1) determine the structural factor that burner need to be considered, and it is carried out to transcoding, coding transform;
Described structural factor is fire trough, spray gun, stabilizer, fire trough shape, diameter, Burners Positions, nozzle angle, one or several combinations in nozzle diameter;
Described transcoding, coding transform is: in response surface design, the variation range of each structural factor is different, and the scope difference of some independent variable is great disparity extremely; For the unified convenience of processing, the value of design parameter is made to transcoding, coding transform, set up the one-to-one relationship of structural factor level value and " coding "; Coding is the trouble of bringing to design and analysis in order to solve dimension difference;
(2) adopt Central Composite design CCD to determine testing program;
CCD=central composite design is writing a Chinese character in simplified form of Central Composite design;
Described Central Composite design CCD refers to: Central Composite design is a kind of second order design of classics, and testing site comprises a cube point, three kinds of central point and axial point;
Be specially with three factors central composite design, it is 1 cubical summit that cube point is distributed in the length of side, and central point is positioned at cube center, axial point claims again initial point, and asterisk point is distributed in the axial direction through central point, except a coordinate be+α or-α, all the other coordinates are all 0;
(3) complete test according to CCD testing program;
Described complete test according to CCD testing program and refer to: according to CCD testing program, adopt hot test or computation fluid dynamics software to complete all cubes of points, the design experiment of central point and axial point, by the NO of flue gas analyzer or computation fluid dynamics aftertreatment acquisition testing site
xnumerical value.
(4) set up response surface model, analytical test result;
The described response surface model of setting up, analytical test result, refers to: Test Point Data is inputted to response surface software, for example, Statistica, Design-Expert or SAS; According to software analysis result, extract multinomial coefficient, obtain the polynomial relation between each factor and response; Output response surface three-dimensional plot also judges the affect situation of each factor on objective function; According to representative value analytical table, judge whether to carry out ridge analysis;
(5) carry out ridge analysis, determine the parameter combinations of optimizing structure.
Described ridge analysis, refer to: taking original design central point as the centre of sphere, each coordinate is started to continuous expansion from the centre of sphere, within the scope of the track that the hypersphere that is radius and the intersection point of response surface form, find out optimal response value at R, wherein, radius R can not exceed trial stretch.
Described intersection point has another name called ridge point;
Compared with prior art, good effect of the present invention is:
The present invention is better than traditional single factors optimization method, such as single factors optimization method cannot be determined the cross influence between injection diameter, spray angle, eddy flow angle like this.
The present invention designs CCD method by Central Composite and determines testing program, and by response surface software analysis test findings, not only can determine linear term impact and the quadratic term impact of each factor on burner combustion performance, can also determine the cross term impact of each factor, then determine by ridge analysis the structural parameters combination of optimizing, thereby realize low nitrogen burning effect.
[brief description of the drawings]
Fig. 1 tri-factor CCD testing sites distribute;
Fig. 2 response surface three-dimensional plot;
(a) H, α affect response;
(b) H, β affect response;
(c) H, γ affect response;
(d) α, β affect response;
(e) α, γ affect response;
(f) β, γ affect response;
The each structure NO cloud atlas contrast of Fig. 3;
(a) optimize structure;
(b) 29th structure;
(c) 16th structure.
[embodiment]
The present invention is below provided a kind of embodiment of the burner structure optimization method based on response surface method.
Embodiment 1
A burner structure optimization method based on response surface method, its concrete steps are:
(1) determine the structural factor that burner need to be considered, and it is carried out to transcoding, coding transform;
Described structural factor is fire trough, spray gun, stabilizer, fire trough shape, diameter, Burners Positions, nozzle angle, one or several combinations in nozzle diameter;
Described transcoding, coding transform is: in response surface design, the variation range of each structural factor is different, and the scope difference of some independent variable is great disparity extremely; For the unified convenience of processing, the value of design parameter is made to transcoding, coding transform, set up the one-to-one relationship of structural factor level value and " coding "; Coding is the trouble of bringing to design and analysis in order to solve dimension difference;
(2) adopt Central Composite design CCD to determine testing program;
CCD=central composite design is writing a Chinese character in simplified form of Central Composite design;
Described Central Composite design CCD refers to: Central Composite design is a kind of second order design of classics, and testing site comprises a cube point, three kinds of central point and axial point;
Be specially with three factors central composite design, it is 1 cubical summit that cube point is distributed in the length of side, and central point is positioned at cube center, axial point claims again initial point, and asterisk point is distributed in the axial direction through central point, except a coordinate be+α or-α, all the other coordinates are all 0;
(3) complete test according to CCD testing program;
Described complete test according to CCD testing program and refer to: according to CCD testing program, adopt hot test or computation fluid dynamics software to complete all cubes of points, the design experiment of central point and axial point, by the NO of flue gas analyzer or computation fluid dynamics aftertreatment acquisition testing site
xnumerical value.
(4) set up response surface model, analytical test result;
The described response surface model of setting up, analytical test result, refers to: Test Point Data is inputted to response surface software, for example, Statistica, Design-Expert or SAS; According to software analysis result, extract multinomial coefficient, obtain the polynomial relation between each factor and response; Output response surface three-dimensional plot also judges the affect situation of each factor on objective function; According to representative value analytical table, judge whether to carry out ridge analysis;
(5) carry out ridge analysis, determine the parameter combinations of optimizing structure.
Described ridge analysis, refer to: taking original design central point as the centre of sphere, each coordinate is started to continuous expansion from the centre of sphere, within the scope of the track that the hypersphere that is radius and the intersection point of response surface form, find out optimal response value at R, wherein, radius R can not exceed trial stretch.
Specific as follows:
1, determine the structural factor that needs consideration
The present embodiment is mainly optimized the burner structure of certain reboiler furnace, and this burner is combustion gas gradual combustor.By the topology discovery of analytic combustion device, four factors such as main and auxiliary spray gun relative height H, main jet rifle nozzle angle α, auxiliary spray tip angle β, γ are the principal element that affects burner combustion performance.Consider the problems such as computing time, such as air channel length etc. of other factors do not considered for the time being.
2, transcoding, coding transform
The rule of the present embodiment explanation transcoding, coding transform as an example of the transcoding, coding transform of H example.Consider the size of burner and burner hearth, the variation range of H is (H
min=130, H
max180), the transcoding, coding transform of H is as follows:
Half interval length
Central point
The coding formula of H is
By this transcoding, coding transform, just can be by interval (H
min, H
max) be converted into (1,1).
In like manner, can obtain the transcoding, coding transform result of other factors, visible table 1.
The each factor transcoding, coding transform of table 1 result
3, the design of CCD Central Composite and CFD simulation
At present, CCD uses to obtain test design the most widely, and testing site comprises three kinds of a cube point, central point and axial point, the attached distribution that Figure 1 shows that three factor CCD testing sites.The present embodiment adopts the design of CCD Central Composite, and test sequence is random.After completing the design of CCD Central Composite, according to design table (schedule), utilize CFD software Fluent to complete whole simulation tests.The design of CCD Central Composite and test findings are as shown in table 2.
Table 2CCD Central Composite design and test findings
4, set up response surface model, analytical test result
Adopt Statistica software to carry out response surface analysis to test findings, obtain the regression equation taking NO discharge capacity as response:
Y＝9.39H+4.01H
2-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βγ
The visible table 3 of analysis result, from this table, during taking NO discharge capacity as response, model P < 0.0001, shows that this quadratic equation model is remarkable, and fitting effect is good, and reliability is high.In the time of P < 0.05, represent that this index is remarkable, when P < 0.001, show that this index is very remarkable.Visible, H linear term, H quadratic term and β linear term are extremely significant model terms, H, β product term are significant model terms, lose and intend not remarkable (P > 0.1) of item, regression equation degree of fitting and confidence level are higher, test error is less, illustrates that this model can be used for instructing burner structure optimal design.
Table 3 response surface method analysis result
The figure of response surface analysis method is the three dimensions that each factor corresponding to specific response forms, and can reflect intuitively the impact of each factor on response, from the response surface analysis figure of test gained, can analyze the interaction between each factor.Because response surface figure is three-dimensional plot, can only express the function containing two influence factors, therefore each response surface three-dimensional plot needs fixing wherein two influence factors to be positioned at the position that is encoded to 0.The response surface three-dimensional plot of each factor interaction is shown in accompanying drawing 2.
5, ridge analysis
Can find out from representative value analytical table 4, the eigenwert of four factors has just to be had negatively, shows that this quadratic response face is saddle face, there is no unique optimum value, therefore needs to do ridge analysis, further determines optimal response value.
Table 4 representative value analytical table
Ridge analysis is that (this test goes the centre of sphere as <-1 taking the original design central point centre of sphere, 0,0,0>), within the scope of the track that the hypersphere that is radius at R and the intersection point of response surface (being ridge point) form, find out optimal response value.The result of ridge analysis is that each coordinate is started to continuous expansion from the centre of sphere, and radius R can not exceed trial stretch, choose herein R 0.1,0.2 ..., calculate ridge point in 0.8 scope.The present embodiment adopts SAS software to carry out ridge analysis.
Table 5 ridge analysis
As shown in Table 5, along with the increase of coding radius R, response NO reduces gradually.In the time of R > 0.8, standard error is larger, does not consider.In this test level scope, R=0.8 when minimum response value, now corresponding non-coding physical dimension is H=140, α=89, β=29, γ=92.Utilize Fluent to do checking analog computation according to above-mentioned optimum structure parameter, obtaining NO realistic simulation result is 30.7ppm, basically identical with predicted value.Accompanying drawing 3 for optimize structure with CCD design table in the NO cloud atlas contrast of the 16th group of structure and the 29th group of structure, can find that the NO numerical value of optimizing structure significantly reduces, significantly realized NO emission reduction effect, reached project expection.
The above is only the preferred embodiment of the present invention; it should be pointed out that for those skilled in the art, without departing from the inventive concept of the premise; can also make some improvements and modifications, these improvements and modifications also should be considered within the scope of protection of the present invention.
Claims (7)
1. the burner structure optimization method based on response surface method, is characterized in that, concrete steps are:
(1) determine the structural factor that burner need to be considered, and it is carried out to transcoding, coding transform;
(2) adopt Central Composite design CCD to determine testing program;
(3) complete test according to CCD testing program;
(4) set up response surface model, analytical test result;
(5) carry out ridge analysis, determine the parameter combinations of optimizing structure;
Described ridge analysis, refer to: taking original design central point as the centre of sphere, each coordinate is started to continuous expansion from the centre of sphere, within the scope of the track that the hypersphere that is radius and the intersection point of response surface form, find out optimal response value at R, wherein, radius R can not exceed trial stretch.
2. a kind of burner structure optimization method based on response surface method as claimed in claim 1, it is characterized in that, in described step (1), described structural factor is fire trough, spray gun, stabilizer, fire trough shape, diameter, Burners Positions, nozzle angle, one or several combinations in nozzle diameter.
3. a kind of burner structure optimization method based on response surface method as claimed in claim 1, it is characterized in that, in described step (1), described transcoding, coding transform is: in response surface design, the variation range of each structural factor is different, and the scope difference of some independent variable is great disparity extremely; For the unified convenience of processing, the value of design parameter is made to transcoding, coding transform, set up the one-to-one relationship of structural factor level value and " coding ".
4. a kind of burner structure optimization method based on response surface method as claimed in claim 1, it is characterized in that, in described step (2), described Central Composite design CCD refers to: Central Composite design is a kind of second order design of classics, testing site comprises a cube point, three kinds of central point and axial point.
5. a kind of burner structure optimization method based on response surface method as claimed in claim 1, it is characterized in that, in described step (2), described Central Composite design CCD is specially with three factors central composite design, and it is 1 cubical summit that cube point is distributed in the length of side, central point is positioned at cube center, axial point claims again initial point, and asterisk point is distributed in the axial direction through central point, except a coordinate be+α or-α, all the other coordinates are all 0.
6. a kind of burner structure optimization method based on response surface method as claimed in claim 1, it is characterized in that, in described step (3), described complete test according to CCD testing program and refer to: according to CCD testing program, adopt hot test or computation fluid dynamics software to complete all cubes of points, the design experiment of central point and axial point, by the NO of flue gas analyzer or computation fluid dynamics aftertreatment acquisition testing site
xnumerical value.
7. a kind of burner structure optimization method based on response surface method as claimed in claim 1, it is characterized in that, in described step (4), the described response surface model of setting up, analytical test result, refers to: Test Point Data is inputted to response surface software, according to software analysis result, extract multinomial coefficient, obtain the polynomial relation between each factor and response; Output response surface three-dimensional plot also judges the affect situation of each factor on objective function; According to representative value analytical table, judge whether 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 true CN104036075A (en) | 2014-09-10 |
| CN104036075B 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) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104732084A (en) * | 2015-03-20 | 2015-06-24 | 上海交通大学 | Model construction method for reducing rapidly digestible starch content on basis of response surface method |
| CN107220400A (en) * | 2017-04-05 | 2017-09-29 | 中国石油化工股份有限公司 | A kind of high-temperature air burning organization optimization method based on response phase method |
| CN113569448A (en) * | 2021-07-13 | 2021-10-29 | 沈阳航空航天大学 | Response surface method-based strain gauge sensitive grid structure parameter optimization method |
| CN114636169A (en) * | 2022-02-21 | 2022-06-17 | 北京航空航天大学 | Flame stabilizer perforation parameter determination method and device and radial flame stabilizer |
Citations (3)
| 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 |
| CN103324798A (en) * | 2013-06-25 | 2013-09-25 | 福州大学 | Random model updating method based on interval response surface model |
| CN103366043A (en) * | 2013-05-31 | 2013-10-23 | 天津大学 | Three-point support optimization design method of machine tool bed based on response surface model |
-
2014
- 2014-06-05 CN CN201410246922.0A patent/CN104036075B/en not_active Expired - Fee Related
Patent Citations (3)
| 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 |
| CN103366043A (en) * | 2013-05-31 | 2013-10-23 | 天津大学 | Three-point support optimization design method of machine tool bed based on response surface model |
| CN103324798A (en) * | 2013-06-25 | 2013-09-25 | 福州大学 | Random model updating method based on interval response surface model |
Non-Patent Citations (1)
| Title |
|---|
| 刘波等: "管式加热炉空气分级燃烧器的CFD研究", 《石油学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104732084A (en) * | 2015-03-20 | 2015-06-24 | 上海交通大学 | Model construction method for reducing rapidly digestible starch content on basis of response surface method |
| 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 |
| CN113569448A (en) * | 2021-07-13 | 2021-10-29 | 沈阳航空航天大学 | Response surface method-based strain gauge sensitive grid structure parameter optimization method |
| CN113569448B (en) * | 2021-07-13 | 2023-12-19 | 沈阳航空航天大学 | Optimization method for strain gauge sensitive gate structural parameters based on response surface method |
| CN114636169A (en) * | 2022-02-21 | 2022-06-17 | 北京航空航天大学 | Flame stabilizer perforation parameter determination method and device and radial flame stabilizer |
| CN114636169B (en) * | 2022-02-21 | 2023-01-20 | 北京航空航天大学 | Flame stabilizer perforation parameter determination method and device and radial flame stabilizer |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104036075B (en) | 2018-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103235842B (en) | The acquisition methods of quadrangle tangential circle coal-fired boiler combustion characteristic and system | |
| Noor et al. | A preliminary study of control parameters for open furnace MILD combustion using CFD | |
| CN104036075A (en) | Combustor structure optimization method based on response surface method | |
| CN103225819B (en) | Pulverized-coal fired boiler changes the oxygen amount method of adjustment after coal-fired kind and system | |
| CN103234219B (en) | Pulverized-coal fired boiler changes the surrounding air amount adjustment method after coal-fired kind and system | |
| CN103278599A (en) | Low-density gas burner, burning test system and controlling method thereof | |
| CN109359268B (en) | Method for calculating penetration distance of heavy oil spray of marine diesel engine | |
| Elorf et al. | Effect of swirl strength on the flow and combustion characteristics of pulverized biomass flames | |
| CN103216846B (en) | Smoke concentration adjustment method and smoke concentration adjustment system when fire coal variety of pulverized coal boiler is changed | |
| 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 | |
| Klos et al. | Investigation of the effect of injection and control strategies on combustion instability in reactivity-controlled compression ignition engines | |
| Liu et al. | Numerical simulation of urea based SNCR process in a trinal-sprayed precalciner | |
| CN204003161U (en) | A kind ofly realize mixed uniformly gas engine mixer | |
| CN107220400A (en) | A kind of high-temperature air burning organization optimization method based on response phase method | |
| Chen et al. | The influence of fuel bias in the primary air duct on the gas/particle flow characteristics near the swirl burner region | |
| CN103324991B (en) | The method and system of pulverized-coal fired boiler NOx emission predictive | |
| Sun et al. | Experiment and Numerical Simulation Study of Low-nitrogen Combustion Technology inside Small Gas Boiler | |
| CN103528837B (en) | The method for testing performance of porous media combustor and device | |
| Bolegenova et al. | Staged supply of fuel and air to the combustion chamber to reduce emissions of harmful substances | |
| Wang et al. | The effect of in-cylinder temperature on the ignition initiation location of a pre-chamber generated hot turbulent jet | |
| Andreini et al. | Large Eddy Simulation Based Computational Fluid Dynamics Investigation of the Ignition Process in Lean Spray Burner | |
| Yang et al. | CFD modelling and numerical investigation of a large marine two-stroke dual fuel direct injection engine | |
| Fomenko et al. | Numerical study of the methane-air combustion in the direct-flow burner of the boiler 300 MW TGMP-314 boiler | |
| Wang et al. | Gas-Solid Flow and High-Efficiency and Low-NO x Combustion of a Novel Arch-Fired Boiler with Swirl Burners: Impact of Over-Fire Air |
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 |
Granted publication date: 20180209 Termination date: 20190605 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |