CN109767815A - The method simplified based on the combustion reaction mechanisms under rate uncertainty - Google Patents

The method simplified based on the combustion reaction mechanisms under rate uncertainty Download PDF

Info

Publication number
CN109767815A
CN109767815A CN201811585839.0A CN201811585839A CN109767815A CN 109767815 A CN109767815 A CN 109767815A CN 201811585839 A CN201811585839 A CN 201811585839A CN 109767815 A CN109767815 A CN 109767815A
Authority
CN
China
Prior art keywords
reaction
combustion
combustion reaction
model
under
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
Application number
CN201811585839.0A
Other languages
Chinese (zh)
Other versions
CN109767815B (en
Inventor
周伟星
林圣强
韩之雄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201811585839.0A priority Critical patent/CN109767815B/en
Publication of CN109767815A publication Critical patent/CN109767815A/en
Application granted granted Critical
Publication of CN109767815B publication Critical patent/CN109767815B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Feedback Control In General (AREA)

Abstract

The present invention is the method simplified based on the combustion reaction mechanisms under rate uncertainty.PSR zero dimension homogeneous combustion device fuel combustion process is simulated according to detailed reaction mechanism model, the characteristic parameter of production concentration, temperature variations and ignition delay time under acquisition setting condition, using the component of GDRG method choice close coupling, pass through FAST method, based under rate condition of uncertainty, global sensitivity analysis is carried out to combustion reaction mechanisms, combustion reaction mechanisms velocity coefficient is optimized using least square method, so that Reduced mechanisms precisely predict that burning performance parameter, last optimizing application algorithm optimize model.The time needed for numerical value calculates is greatly reduced under the premise of guaranteeing that reaction mechanism is complete reliable in the present invention.

Description

The method simplified based on the combustion reaction mechanisms under rate uncertainty
Technical field
It is the burning under a kind of uncertainty based on rate the present invention relates to the computer modeling technique field of combustion process The method that reaction mechanism simplifies.
Background technique
The burning of fuel provides power for most of equipment in production and life, and the burning of hydrocarbon fuel is shared wherein Ratio is more than 60%.With the gradual depletion of fossil energy with the increasingly raising of environmental requirement, the burning effect of fuel how is improved Rate and the pollution for reducing burning generation have become the task of top priority of combustion field.The burning of fuel is related to flowing, and conducts heat and changes The coupling for learning more of reaction etc., is an extremely complex process.Computer technology rises to the simulation for solving fuel combustion Possibility is provided, especially in terms of the design of burner, the control of the tissue and pollutant of combustion process, numerical value calculating is played the part of The role drilled is more and more important.It is related to hundreds and thousands of intermediate products and elementary reaction in the combustion process of fuel, quantity is such Huge process makes the convergence rate of entire calculating process extremely slow, or even the case where not restraining occurs.Therefore, it is calculating In the case that precision is met the requirements, intermediate product and elementary reaction for certain contents lower than given threshold then be can be ignored not Meter, so that the calculating process of reaction is simplified, the calculating time is substantially reduced.Therefore, the simplification of the combustion mechanism of hydrocarbon fuel Important research direction as combustion field.
At this stage, the mechanism simplification method of burning mainly has two major classes: the first seed type is then the number for reducing elementary reaction Amount carries out, such as SA and PCA method, and the second class method is then the component reduced in reaction Kinetics Model, such as DRG, PFA and DRGEP etc..In the biggish situation of reaction Kinetics Model, two methods are generally used in conjunction with each other, i.e. the first stage first leads to The method for crossing screen fraction simplifies model, and second stage is then by deleting the chemical reaction of redundancy with other means Subtract, finally reaches the purpose of mechanism simplification.Here it is noted that either elementary reaction delete or intermediate product is deleted Subtract, is required to follow certain reaction rule, the reaction mechanism of finally obtained simplification is also in accordance with reaction rule.
Although the method that above-mentioned combustion mechanism simplifies can substantially combustion reaction mechanism, numerical value is calculated and Speech simplifies obtained result and still has two o'clock deficiency: first, Reduced mechanisms obtained by the above method are for special circumstances Under, the applicability of combustion case especially within the engine is unsatisfactory, especially under the mechanism application conditions of big component, Its effect is poor.Second it is that the method that above-mentioned reaction mechanism simplifies does not consider the uncertainty for reflecting velocity coefficient to mechanism Influence, may cause unessential component and elementary reaction be reinforced, and important component and elementary reaction is weakened.This kind Situation, which can seriously affect, characterizes important elementary reaction and main response path in fuel combustion process, and mechanism simplification Amplitude is bigger, then the uncertainty of mechanism is higher
Summary of the invention
The present invention is to solve existing issue, provides what the combustion reaction mechanisms under a kind of uncertainty based on rate simplified Method, the present invention provides following technical schemes:
A kind of combustion reaction mechanisms under the uncertainty based on rate simplify method, comprising the following steps:
Step 1: simulating PSR zero dimension homogeneous combustion device fuel combustion process according to detailed reaction mechanism model, The characteristic parameter of production concentration, temperature variations and ignition delay time under acquisition setting condition;
Step 2: 1000 combustion mechanism models are simplified by GDRG method, to the detailed reaction mechanism model In the combination of non-close coupling deleted;
Step 3: by FAST method, under rate condition of uncertainty, global sensibility is carried out to combustion reaction mechanisms Analysis, according to the sensitivity coefficient of the obtained each elementary reaction of FAST method, set sensitivity coefficient threshold value to elementary reaction into Row is deleted;
Step 4: combustion reaction mechanisms velocity coefficient is optimized using least square method, so that Reduced mechanisms are accurate Predict burning performance parameter;
Step 5: by being compared with burning detailed reaction mechanism model, after mechanism obtained is carried out numerical value calculating The ignition delay time that arrives to ignition delay time and burning detailed reaction mechanism compares, when error is in claimed range, then It is not required to adjust, requires, then the sensitivity coefficient threshold value of setting is adjusted when error is greater than.
Preferably, the step 2 specifically:
Step 1: counting the probability of non-coupled component in the way of probability by GDRG method;
Step 2: 1000 sample points of uniform sampling in the parameter space of uncertain velocity coefficient;
Step 3: simplifying using DRG method to 1000 combustion mechanism models, to the non-coupled combination in model It is deleted.
Preferably, the step 3 specifically:
Step 1: measuring the global sensibility system of each elementary reaction in a manner of disturbing elementary reaction by FAST method Number,;
Step 2: setting sensitivity coefficient threshold value, threshold value elementary reaction below is deleted;
Step 3: obtaining the global sensitivity coefficient of elementary reaction by following formula:
Wherein, yiFor the parameter (such as concentration, temperature) of the i-th outlet;NsSegmentation to be sampled to curve is counted, Middle sampled point presses s on curvek=2k π/Ns, k=1,2,3 ... NsIt is sampled;It is the Fourier of the i-th outlet parameter Coefficient.
Preferably, sensitivity coefficient is represented using Fourier coefficient.
Preferably, the step 4 specifically:
Step 1: being optimized using method of least squares to objective function, so that the function f of velocity coefficientt(Aii,Ei) become In 0;
Step 2: being calculate by the following formula to obtain corresponding suitable velocity coefficient parameter Aii,Ei:
Wherein, Ai, βi,EiFor the pre-exponential factor during combustion reaction, dimensionless group and reaction activity;Yk,DAnd TD Respectively indicate the outgoing quality score obtained by detailed mechanism and temperature;WithIt respectively indicates and passes through PSR equation calculation obtain in parameters and temperature such as the concentration in t=τ moment reaction device exit;
Step 3: establishing PSR equation by following formula:
Wherein, WkFor the molal weight (g/mol) of k-th of component, ωkFor the molar yield of k-th of component, ρ is density, hk,iAnd hk,oRespectively indicate the enthalpy of low k component import and export, k=1,2...K.
Preferably, during combustion reaction mechanisms model simplification, sensitivity coefficient threshold size is according to simplified mould The width of the type scope of application is adjusted.
The invention has the following advantages:
The present invention is based on simplifying under rate condition of uncertainty to the combustion mechanism of hydrocarbon fuel, tradition can be avoided There is the phenomenon that inaccurately deleting component during simplification in method.By GDRG method to the component in detailed mechanism model It is screened, obtains the component of close coupling, so that the component of non-close coupling can be deleted in wide scope threshold value, make mould The simplification that type obtains.The sensitivity analysis of overall importance that disturbance indirectly obtains elementary reaction, and and GDRG are carried out by FAST method Method combines, and can cut down the quantity of component in reaction to the greatest extent.Elementary reaction is deleted by FAST method Later, so that the quantity of the elementary reaction in simplified model is reduced, model is further simplified.By least square method to filtering out Elementary reaction in reflection velocity coefficient reasonably optimized so that the knot of the result of simplified model and detailed mechanism model Error between fruit minimizes.The present invention can be adjusted during the simplification of mechanism model, be guaranteed to the greatest extent The reliability and accuracy of model.When the parameter relevant to ignition quality that numerical simulation result is characterized is mainly firing delay Between, production concentration distribution and variation in reactor etc..Sampled data in detailed mechanism model is that above-mentioned calculating process mentions For foundation, in allowable range of error, also can sampled data be increased or be reduced, when providing more suitable firing delay Between, product distribution with variation etc. typical combustions characteristic parameter.
Detailed description of the invention
Fig. 1 is the flow chart of the method simplified based on the combustion reaction mechanisms under rate uncertainty.
Fig. 2 is the schematic diagram that component is coupled by elementary reaction in reaction.
Fig. 3 is that the wide scope down-firing being calculated using the obtained methyl hydride combustion mechanism model of the present invention is postponed The comparing result of time and document.
Fig. 4 is that the close limit down-firing being calculated using the obtained methyl hydride combustion mechanism model of the present invention is postponed The comparing result of time and document.
Specific embodiment
Below in conjunction with specific embodiment, describe the invention in detail.
Specific embodiment one:
According to Fig. 1, the present invention provides the method that the combustion reaction mechanisms under a kind of uncertainty based on rate simplify, The following steps are included:
Step 1: simulating PSR zero dimension homogeneous combustion device fuel combustion process according to detailed reaction mechanism model, The characteristic parameter of production concentration, temperature variations and ignition delay time under acquisition setting condition;
Step 2: 1000 combustion mechanism models are simplified by GDRG method, to the detailed reaction mechanism model In the combination of non-close coupling deleted;
Step 3: by FAST method, under rate condition of uncertainty, global sensibility is carried out to combustion reaction mechanisms Analysis, according to the sensitivity coefficient of the obtained each elementary reaction of FAST method, set sensitivity coefficient threshold value to elementary reaction into Row is deleted;
Step 4: combustion reaction mechanisms velocity coefficient is optimized using least square method, so that Reduced mechanisms are accurate Predict burning performance parameter;
Step 5: by being compared with burning detailed reaction mechanism model, after mechanism obtained is carried out numerical value calculating The ignition delay time that arrives to ignition delay time and burning detailed reaction mechanism compares, when error is in claimed range, then It is not required to adjust, requires, then the sensitivity coefficient threshold value of setting is adjusted when error is greater than.
Specific embodiment two:
The present embodiment is the GRI3.0 combustion mechanism for not considering NO_x formation, which includes 35 components and 217 A reaction.It under certain condition, can the Accurate Prediction spark delay time using the simplification reaction mechanism that DRG method obtains.
The document A criterion based on for being Lu et al. is compared with calculated result of the invention computational singular perturbation for the identification of quasi steady The meter of state species:A reduced mechanism for methane oxidation with NO chemistry Calculate result.
Step 1, detailed mechanism the model calculation: the calculating primary condition of consideration are as follows: pressure P=1~30atm, equivalent Than Phi=0.5~1.5, temperature range is T=1000~1600K.Multiple temperature spots during this are calculated, point is obtained The typical combustions characteristic parameter such as ignition delay time.
Step 2 passes through GDRG (global DRG) method, and the probability of non-coupled component is counted in the way of probability, grasps It is as follows to make process: 1000 sample points of uniform sampling first in the parameter space of uncertain velocity coefficient;Then DRG is applied Method simplifies 1000 combustion mechanism models, to obtain the probability that each component is non-close coupling component.Herein it Preceding to need to set a threshold value between 0 and 1, the component for being greater than the threshold value retains, for being less than the threshold value Component further delete.According to the non-close coupling probability under condition of uncertainty, which is less than to the coupling probability in model The component of value is deleted.Methane and oxygen are critical reactants in the present embodiment.
Step 3, component global sensitivity analysis: by FAST Fourier's overall situation susceptibility process, to disturb elementary reaction Mode measure the global sensitivity coefficient of each elementary reaction, according to the sensitivity coefficient threshold value of setting, below the threshold value Elementary reaction deleted.The global sensitivity coefficient of elementary reaction is obtained in the present invention by following formula:
In above formula, yiFor the parameter (such as concentration, temperature) of the i-th outlet;NsSegmentation to be sampled to curve is counted, Wherein sampled point presses s on curvek=2k π/Ns, k=1,2,3 ... NsIt is sampled;Be the i-th outlet parameter Fu in Leaf system number.In the present invention, sensitivity coefficient is represented using Fourier coefficient.
Finally on the basis of above-mentioned documents, mechanism model is further simplified, the group further deleted Divide as shown in the table.
Step 4, mechanism optimization: ignition delay time etc. is joined to what is chosen by step 3 using least square method There is number the elementary reaction of great influence to optimize, so that the analog result of Reduced mechanisms is forced to the analog result of detailed mechanism Closely.Following objective function is optimized using least square method in the present invention so that velocity coefficient function ft (Ai, β i, Ei tend to 0) to find most suitable velocity coefficient parameter Ai, β i, Ei.
In above formula, Ai, βi,EiFor the pre-exponential factor during combustion reaction, dimensionless group and reaction activity;Yk,DWith TDRespectively indicate the outgoing quality score obtained by detailed mechanism and temperature;WithIt respectively indicates and passes through PSR equation calculation obtain in parameters and temperature such as the concentration in t=τ moment reaction device exit.
Step 5 is adjusted by different threshold values, obtains the simplification machine of different component quantity Yu elementary reaction quantity Model is managed,
PSR equation are as follows:
In above formula, WkFor the molal weight (g/mol) of k-th of component, ωkFor the molar yield of k-th of component, ρ is close Degree, hk,iAnd hk,oRespectively indicate the enthalpy of low k component import and export.
Process is deleted by above-mentioned, so that have 19 components in the model simplified, 99 reactions, by selection, for it In 33 elementary reactions optimize.
Methyl hydride combustion is calculated using above-mentioned model, calculated result and bibliography are compared, result such as Fig. 3 It is shown, under wide scope the ignition delay time of (P=1~30atm, Phi=0.5~1.5, T=1000~1600K) and in detail The result worst error of mechanism model is 4%, thus under the conditions of the present embodiment, the calculated result of the reaction model of methyl hydride combustion Precision with higher.
Specific embodiment three:
Compared with calculated result of the invention be GRI3.0 combustion mechanism calculated result, this example is for methane Combustion model under specific operation is simplified.
In the present embodiment, primary condition is calculated are as follows: pressure P=1atm, equivalent proportion Phi=1, temperature range T=1000 ~1600K.Multiple temperature spots during this are calculated, the typical combustions characteristic parameters such as ignition delay time are obtained.
Similar, former mechanism model carries out further component according to step 2 in above-mentioned specific embodiment two to step 3 It deletes, final to delete wherein 6 components, last mechanism model is 16 components, and result is as shown in the table.
It repeats step 4 in above-mentioned specific embodiment two and arrives step 5, so that there is 16 components in the model simplified, 58 anti- It answers, by selection, is optimized for 33 reactions therein.Methyl hydride combustion is calculated using above-mentioned model, will be calculated As a result compare with bibliography, result as shown in figure 4, under close limit (P=1atm, Phi=1, T=1000~ The result of ignition delay time and detailed mechanism model 1600K) is very close, and precision is higher.The thus mechanism in close limit The simplification of model can more have specific aim, and computational accuracy is also higher.
The above is only the preferred embodiment of the method simplified based on the combustion reaction mechanisms under rate uncertainty, Protection scope based on the simplified method of the combustion reaction mechanisms under rate uncertainty is not limited merely to above-described embodiment, all The technical solution belonged under thinking all belongs to the scope of protection of the present invention.It should be pointed out that for those skilled in the art For, several improvements and changes without departing from the principles of the present invention, such modifications and variations also should be regarded as of the invention Protection scope.

Claims (6)

1. the method that the combustion reaction mechanisms under a kind of uncertainty based on rate simplify, it is characterized in that: the following steps are included:
Step 1: simulating PSR zero dimension homogeneous combustion device fuel combustion process according to detailed reaction mechanism model, obtains The characteristic parameter of production concentration, temperature variations and ignition delay time under imposing a condition;
Step 2: 1000 combustion mechanism models are simplified by GDRG method, in the detailed reaction mechanism model The combination of non-close coupling is deleted;
Step 3: by FAST method, under rate condition of uncertainty, global sensibility point is carried out to combustion reaction mechanisms Analysis sets sensitivity coefficient threshold value and carries out to elementary reaction according to the sensitivity coefficient of the obtained each elementary reaction of FAST method It deletes;
Step 4: combustion reaction mechanisms velocity coefficient is optimized using least square method, so that Reduced mechanisms are precisely predicted Burning performance parameter;
Step 5: by with burning detailed reaction mechanism model comparison, will mechanism obtained carry out numerical value calculating after arrive point The ignition delay time that arrives of ignition delay time and burning detailed reaction mechanism compares, and when error is in claimed range, is then not required to Adjustment is required when error is greater than, is then adjusted to the sensitivity coefficient threshold value of setting.
2. the method that the combustion reaction mechanisms under a kind of uncertainty based on rate according to claim 1 simplify, special Sign is: the step 2 specifically:
Step 1: counting the probability of non-coupled component in the way of probability by GDRG method;
Step 2: 1000 sample points of uniform sampling in the parameter space of uncertain velocity coefficient;
Step 3: simplifying using DRG method to 1000 combustion mechanism models, the non-coupled combination in model is carried out It deletes.
3. the method that the combustion reaction mechanisms under a kind of uncertainty based on rate according to claim 1 simplify, special Sign is: the step 3 specifically:
Step 1: the global sensitivity coefficient of each elementary reaction is measured in a manner of disturbing elementary reaction by FAST method,;
Step 2: setting sensitivity coefficient threshold value, threshold value elementary reaction below is deleted;
Step 3: obtaining the global sensitivity coefficient of elementary reaction by following formula:
Wherein, yiFor the parameter (such as concentration, temperature) of the i-th outlet;NsSegmentation to be sampled to curve is counted, wherein adopting Sampling point presses s on curvek=2k π/Ns, k=1,2,3 ... NsIt is sampled;It is the Fourier system of the i-th outlet parameter Number.
4. the method that the combustion reaction mechanisms under a kind of uncertainty based on rate according to claim 3 simplify, special Sign is: representing sensitivity coefficient using Fourier coefficient.
5. the method that the combustion reaction mechanisms under a kind of uncertainty based on rate according to claim 1 simplify, special Sign is: the step 4 specifically:
Step 1: being optimized using method of least squares to objective function, so that the function f of velocity coefficientt(Aii,Ei) tend to 0;
Step 2: being calculate by the following formula to obtain corresponding suitable velocity coefficient parameter Aii,Ei:
Wherein, Ai, βi,EiFor the pre-exponential factor during combustion reaction, dimensionless group and reaction activity;Yk,DAnd TDRespectively Indicate the outgoing quality score obtained by detailed mechanism and temperature;WithIt respectively indicates through the side PSR Journey be calculated in parameters and temperature such as the concentration in t=τ moment reaction device exit;
Step 3: establishing PSR equation by following formula:
Wherein, WkFor the molal weight (g/mol) of k-th of component, ωkFor the molar yield of k-th of component, ρ is density, hk,iWith hk,oRespectively indicate the enthalpy of low k component import and export, k=1,2...K.
6. the method that the combustion reaction mechanisms under a kind of uncertainty based on rate according to claim 1 simplify, special Sign is: during combustion reaction mechanisms model simplification, sensitivity coefficient threshold size is according to the simplified model scope of application Width be adjusted.
CN201811585839.0A 2018-12-25 2018-12-25 Method for simplifying combustion reaction mechanism based on rate uncertainty Active CN109767815B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811585839.0A CN109767815B (en) 2018-12-25 2018-12-25 Method for simplifying combustion reaction mechanism based on rate uncertainty

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811585839.0A CN109767815B (en) 2018-12-25 2018-12-25 Method for simplifying combustion reaction mechanism based on rate uncertainty

Publications (2)

Publication Number Publication Date
CN109767815A true CN109767815A (en) 2019-05-17
CN109767815B CN109767815B (en) 2023-01-24

Family

ID=66450930

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811585839.0A Active CN109767815B (en) 2018-12-25 2018-12-25 Method for simplifying combustion reaction mechanism based on rate uncertainty

Country Status (1)

Country Link
CN (1) CN109767815B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114171128A (en) * 2021-12-14 2022-03-11 西安交通大学 Method for simplifying deep cracking model of hydrocarbon fuel
CN114613443A (en) * 2022-03-10 2022-06-10 哈尔滨工业大学(深圳) Method and device for constructing combustion mechanism model of ammonia-containing fuel
CN115248978A (en) * 2022-07-01 2022-10-28 北京理工大学 Method for simplifying combustion dynamics mechanism of long-chain normal alcohol fuel
CN115292883A (en) * 2022-06-26 2022-11-04 哈尔滨工程大学 Combustion chamber performance online monitoring and predicting method and system based on chemical reactor network method
CN116798531A (en) * 2023-06-29 2023-09-22 西南科技大学 Fusion algorithm-based hydrocarbon fuel mechanism simplification method for supersonic combustion chamber

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0773408A1 (en) * 1995-11-07 1997-05-14 Hitachi, Ltd. Furnace inside state estimation control apparatus of pulverized coal combustion furnace
JP2001143731A (en) * 1999-11-16 2001-05-25 Daikin Ind Ltd Fuel cell system
CN101900349A (en) * 2010-05-27 2010-12-01 中国科学技术大学 Method for using fuse turbulence combustion information to simplify combustion reaction mechanism
CN102799775A (en) * 2012-07-06 2012-11-28 广东电网公司电力科学研究院 Method for calculating temperature of tube wall of platen superheater of ultra supercritical boiler by assistance of numerical simulation
CN104461688A (en) * 2014-11-27 2015-03-25 中国矿业大学 Method for simplifying skeleton of detailed combustion chemical reaction mechanism
CN104880479A (en) * 2015-05-12 2015-09-02 华中科技大学 Two-section type solid fuel staged reaction dynamic analysis equipment
CN106407594A (en) * 2016-10-10 2017-02-15 中国神华能源股份有限公司 Method and device for simulating oxygen-enriched combustion
CN106960092A (en) * 2017-03-22 2017-07-18 哈尔滨工程大学 The regular empirical parameter automatic calibrating method of weber burning
CN107273693A (en) * 2017-06-22 2017-10-20 武汉科技大学 A kind of hydrocarbon fuel mechanism simplification method
CN107832573A (en) * 2017-07-28 2018-03-23 中冶华天工程技术有限公司 Predict the numerical computation method of flue gas recirculation sintering matter thermal coupling process
CN107944207A (en) * 2017-10-31 2018-04-20 南京航空航天大学 3 aviation kerosine alternative fuel of RP simplifies reaction model computational methods
CN108362819A (en) * 2018-03-06 2018-08-03 江苏中烟工业有限责任公司 A method of cigarette burn rate is predicted based on detection Flow Field Distribution

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0773408A1 (en) * 1995-11-07 1997-05-14 Hitachi, Ltd. Furnace inside state estimation control apparatus of pulverized coal combustion furnace
JP2001143731A (en) * 1999-11-16 2001-05-25 Daikin Ind Ltd Fuel cell system
CN101900349A (en) * 2010-05-27 2010-12-01 中国科学技术大学 Method for using fuse turbulence combustion information to simplify combustion reaction mechanism
CN102799775A (en) * 2012-07-06 2012-11-28 广东电网公司电力科学研究院 Method for calculating temperature of tube wall of platen superheater of ultra supercritical boiler by assistance of numerical simulation
CN104461688A (en) * 2014-11-27 2015-03-25 中国矿业大学 Method for simplifying skeleton of detailed combustion chemical reaction mechanism
CN104880479A (en) * 2015-05-12 2015-09-02 华中科技大学 Two-section type solid fuel staged reaction dynamic analysis equipment
CN106407594A (en) * 2016-10-10 2017-02-15 中国神华能源股份有限公司 Method and device for simulating oxygen-enriched combustion
CN106960092A (en) * 2017-03-22 2017-07-18 哈尔滨工程大学 The regular empirical parameter automatic calibrating method of weber burning
CN107273693A (en) * 2017-06-22 2017-10-20 武汉科技大学 A kind of hydrocarbon fuel mechanism simplification method
CN107832573A (en) * 2017-07-28 2018-03-23 中冶华天工程技术有限公司 Predict the numerical computation method of flue gas recirculation sintering matter thermal coupling process
CN107944207A (en) * 2017-10-31 2018-04-20 南京航空航天大学 3 aviation kerosine alternative fuel of RP simplifies reaction model computational methods
CN108362819A (en) * 2018-03-06 2018-08-03 江苏中烟工业有限责任公司 A method of cigarette burn rate is predicted based on detection Flow Field Distribution

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JING WANG等: "Kinetic Characteristics of Coal Char Combustion in Oxygen-Enriched Environment", 《IEEE》 *
侯金丽等: "基于敏感性分析的氧/甲烷燃烧反应简化机理", 《航空动力学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114171128A (en) * 2021-12-14 2022-03-11 西安交通大学 Method for simplifying deep cracking model of hydrocarbon fuel
CN114171128B (en) * 2021-12-14 2024-02-02 西安交通大学 Method for simplifying hydrocarbon fuel deep cracking model
CN114613443A (en) * 2022-03-10 2022-06-10 哈尔滨工业大学(深圳) Method and device for constructing combustion mechanism model of ammonia-containing fuel
CN114613443B (en) * 2022-03-10 2023-04-21 哈尔滨工业大学(深圳) Construction method and device of ammonia-containing fuel combustion mechanism model
CN115292883A (en) * 2022-06-26 2022-11-04 哈尔滨工程大学 Combustion chamber performance online monitoring and predicting method and system based on chemical reactor network method
CN115248978A (en) * 2022-07-01 2022-10-28 北京理工大学 Method for simplifying combustion dynamics mechanism of long-chain normal alcohol fuel
CN116798531A (en) * 2023-06-29 2023-09-22 西南科技大学 Fusion algorithm-based hydrocarbon fuel mechanism simplification method for supersonic combustion chamber
CN116798531B (en) * 2023-06-29 2024-03-22 西南科技大学 Fusion algorithm-based hydrocarbon fuel mechanism simplification method for supersonic combustion chamber

Also Published As

Publication number Publication date
CN109767815B (en) 2023-01-24

Similar Documents

Publication Publication Date Title
CN109767815A (en) The method simplified based on the combustion reaction mechanisms under rate uncertainty
CN107133460A (en) A kind of online dynamic prediction method of boiler flyash carbon content
CN105069537A (en) Constructing method of combined air quality forecasting model
CN111352408B (en) Multi-working-condition process industrial process fault detection method based on evidence K nearest neighbor
CN109508818A (en) A kind of online NOx prediction technique based on LSSVM
CN112163380A (en) System and method for predicting furnace oxygen concentration based on numerical simulation neural network
CN107067100A (en) Wind power anomalous data identification method and device for identifying
CN110675920A (en) MI-LSTM-based boiler NOxPrediction method
CN102982390A (en) Boiler combustion efficiency predicting method based on support vector machine incremental algorithm
US11979697B2 (en) Methods and internet of things systems for obtaining natural gas energy metering component
CN116432123A (en) Electric energy meter fault early warning method based on CART decision tree algorithm
Yan et al. A robust probabilistic wind power forecasting method considering wind scenarios
Ye et al. Modeling and optimization of the NOX generation characteristics of the coal-fired boiler based on interpretable machine learning algorithm
Dashti et al. Molecular descriptors-based models for estimating net heat of combustion of chemical compounds
Li et al. Product tri‐section based crude distillation unit model for refinery production planning and refinery optimization
CN113533622A (en) Coal quality prediction method for coal mill based on neural network
CN106950826A (en) Gasification device PID parameter setting method based on ant colony algorithm
CN114861969A (en) LSTM-based electric power system oscillation risk prediction method
Liu et al. Research on optimization of aircraft climb trajectory considering environmental impact
CN114638169A (en) Method and device for calculating time-varying fault probability of transformer and computer readable storage medium
Vlah et al. Solving bilevel optimal bidding problems using deep convolutional neural networks
Yang et al. ERP and DTW-based Transformer-customer Identification
Kim et al. Feasibility study on AI-based prediction for CRUD induced power shift in PWRs
Deng et al. Research on Soft‐Sensing Methods for Measuring Diene Yields Using Deep Belief Networks
CN109598643A (en) A kind of annealing furnace subregion gas flow distribution method, equipment and storage equipment

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant