CN106503363A - A kind of transient state unstable airflow computational methods of roadway containing high concentration gas - Google Patents
A kind of transient state unstable airflow computational methods of roadway containing high concentration gas Download PDFInfo
- Publication number
- CN106503363A CN106503363A CN201610951446.1A CN201610951446A CN106503363A CN 106503363 A CN106503363 A CN 106503363A CN 201610951446 A CN201610951446 A CN 201610951446A CN 106503363 A CN106503363 A CN 106503363A
- Authority
- CN
- China
- Prior art keywords
- roadway
- branch
- gas
- equation
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 21
- 238000000205 computational method Methods 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 39
- 238000009423 ventilation Methods 0.000 claims abstract description 26
- 239000006185 dispersion Substances 0.000 claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims description 12
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 8
- ONUFESLQCSAYKA-UHFFFAOYSA-N iprodione Chemical compound O=C1N(C(=O)NC(C)C)CC(=O)N1C1=CC(Cl)=CC(Cl)=C1 ONUFESLQCSAYKA-UHFFFAOYSA-N 0.000 claims description 3
- 238000004445 quantitative analysis Methods 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 4
- 238000002485 combustion reaction Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 206010061217 Infestation Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Ventilation (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Feedback Control In General (AREA)
Abstract
Patent of the present invention discloses a kind of transient state unstable airflow computational methods of roadway containing high concentration gas, its computation model is made up of gas dispersion equation and roadway ventilation network equation, gas dispersion equation carries out finite difference using Crank Nicholson methods, obtains roadway point gas density using chasing method;Roadway ventilation network equation solves the air quantity (wind speed) of roadway branch using runge kutta method;Couple solution is carried out to two equations on this basis;Computation model of the present invention considers impact of the high concentration gas to roadway unstable airflow, it is achieved that the distinguished and admirable quantitative Analysis of the branch of roadway containing high concentration gas transient state.
Description
Technical field
The present invention relates to a kind of roadway transient state unstable airflow computational methods, more particularly to a kind of calculating is containing high concentration gas
The method of roadway transient state unstable airflow.
Background technology
The high concentration gas stability distinguished and admirable to roadway produces certain impact;By taking coal and gas prominent as an example, work as protrusion
After generation, high concentration gas accumulate in roadway, cause roadway wind current control, or even induce sending out for the secondary disasters such as gas explosion
Raw.Forefathers have carried out more in-depth study to fire period roadway transient state unstable airflow, and main consideration is acted in fire heating air pressure
Under, the distinguished and admirable infestation regularity of roadway.Zhou Yan, Wang Shengshen are respectively adopted static analysis method(Cloth Drake method)With dynamic-analysis method pair
The critical condition of Parallel Airway During Mine Five in Ascensional Ventilation wind current control has carried out relative analyses;Zhang Xingkai, Wang Zhencai etc. are according to mine
The distinguished and admirable flow performance of the combustion law and ventilating network of fire period, analyzes the distinguished and admirable transient of During Mine Fire Period
Draught principle, establishes mine fire MATHEMATICAL MODEL OF COMBUSTION and distinguished and admirable unsteady flow equation group, and gives calculation method, with
When by work out computer software achieve the simulation to mine fire combustion process and its distinguished and admirable change procedure.Gas air-flow exists
Migration in roadway is transient state, and its impact distinguished and admirable to roadway branch is also dynamic.Chinese scholars are to containing highly concentrated at present
Degree gas roadway transient state unstable airflow Journal of Sex Research is less, a kind of reliable computational methods of proposition of still needing.Result of study is anti-to mine
Calamity rescue, catastrophe are divulged information and effectively prevent secondary disaster accident to be respectively provided with important theoretical and realistic meaning.
Particular content
A kind of in view of deficiency of existing computational methods, the invention provides transient state unstable airflow of roadway containing high concentration gas
Computational methods, the method consideration high concentration gas impact distinguished and admirable to roadway branch transient state, can be realized using the computational methods
Any roadway branch of ventilation network is solved in air quantity Changing Pattern not in the same time.
For solving above-mentioned technical problem, patent of the present invention is adopted the following technical scheme that:Its computation model is by gas dispersion side
Journey
With roadway ventilation network equation
Constitute, wherein formula(1)Middle uiFor the wind speed of roadway branch i, xiIt is the space length along roadway branch i, t is the time, ciFor
The mash gas concentration of roadway branch i, ExFor gas dispersion coefficient, formula(2)Middle qk(k=1,2 ..., it is b) branch's air quantity, RiFor
Branch's i windages,For the ventilation blower blast of branch i,For the gas blast of branch i, cjiFor circuit matrix, cikFor loop square
Battle array cjiTransposed matrix, wherein gas dispersion equation(1)Method for solving is that had initially with Crank-Nicholson methods
Limit difference, then obtains roadway point gas density using chasing method, while roadway ventilation network equation(2)Method for solving is to adopt
Runge-Kutta method solves the air quantity of roadway branch(Wind speed), need to equation in solution procedure(1)And equation(2)Carry out coupling to ask
Solution.
Described gas dispersion coefficient Ex, it is characterised in that gas dispersion coefficientWherein r is pipe half
Footpath, α are coefficient of frictional resistance, and u is roadway branch airflow velocity.
Described branch's air quantity, it is characterised in that branch air quantity qk(k=1,2 ..., b) it is parameter to be asked, and the time
Function, i.e., required value is b unknown function, and b is the number of independent loop in ventilation network;qkIt is that judge roadway is distinguished and admirable steady
Qualitatively index, qk=0 represents distinguished and admirable stagnation, qk< 0 represents wind current control, qk> 0 represents that wind direction keeps constant.
Described ventilation blower blast, it is characterised in that fan pressurea0, a1..., anFor curve
Fitting coefficient;Given three or three blast with fan, air quantity data, solve fitting coefficient by method of least square.
Described gas blastCharacterized in that, For the average density of branch i, For tnThe gas density of j-th spatial mesh size of moment branch, For tnThe mash gas concentration of j-th spatial mesh size of moment branch, ρmAnd ρaThe density of respectively methane and air, h is space
Step-length, LiFor the length of roadway branch i, spatial mesh size numbers of the N (i) for branch i, ziAnd z (0)i(Li) it is respectively the beginning of branch i
End point absolute altitude.
Described gas dispersion equation(1)Method for solving, it is characterised in that entered using Crank-Nicholson methods
Row finite difference, obtains tri-diagonal system, and distinguished and admirable flowing in ventilating system is one-dimensional incompressible flow, each in ventilation network
Last spatial mesh size of tunnel is equal in the gas density of gas density moment back length with certain moment, and gas is released
Mix with distinguished and admirable after putting immediately, and tunnel intersection each distinguished and admirable in gas be mutually mixed immediately, using chasing method solve three pairs
Angle equation group, obtains roadway point gas density.
The method for solving of roadway ventilation network equation (2), it is characterised in that order Each variable is expressed as matrix form is:A=[ajk], C=
[cji], K=[kii], CT=[cik], D=(D1,D2,…Db)T, q=[q1,q2,…,qk]T(k=1,2 ..., b), wherein kii=
Ki, and A=CKCT, formula(2)Being expressed as matrix form isFrom improved Euler's formula(Runge-Kutta method)Solve
The air quantity of each roadway branch, on this basis, according to the basal area of roadway branch, draws the wind speed of roadway branch.
Need to equation in described solution procedure(1)And equation(2)Carry out couple solution, it is characterised in that to described
Formula(1)And formula(2)Carry out couple solution, t0The value of moment each parameter is, it is known that first solve formula(1)t1Moment gas density, then will
The value for solving substitutes into formula(2)Solve t1Moment air quantity, the air quantity for obtaining substitute into formula(1)Solve t2The gas density at moment, successively class
Push away concentration distribution and the air force that each time step tunnel is obtained.
Description of the drawings
Fig. 1 is the transient state unstable airflow computational methods overall schematic of roadway containing high concentration gas of the present invention;
Specific embodiment
The embodiment of the present invention will be described in further detail below;Due to for any roadway branch, gas density
Dispersion Equation is:Due to formula(1)Middle u and ExNot constant, but the function of time, the model is
Nonlinear partial differential equation, its analytic solutions are difficult to calculate, so numerical value of the present invention using the method calculating above formula of finite difference
Solution;In the differential equation, replace partial derivative with difference coefficient, obtain corresponding difference equation, differential side is obtained by solving difference equation
The approximation of journey solution.Before difference is carried out, each tunnel is temporally divided difference gridding with locus first.Assume tunnel
Length is L, divides N sections interval along air flow in laneway direction, then every section of siding-to-siding block length h=L/N, xjIt is that j-th spatial mesh size is saved away from beginning
The distance of point, is expressed as jh;Time same time step τ demarcation interval will be pressed simultaneously, then tnThe corresponding time at moment is n τ, adopts
Use Crank-Nicholson difference schemes.For in roadway, gas density Dispersion Equation presses Crank-Nicholson form difference
For
Abbreviation can be obtained:Wherein:
Gas dispersion model above is set up
On the basis of single tunnel, and any one mine ventilation system all has multi-mine roadway perhaps, distinguished and admirable in each associated tunnel
In migration process, branch has the continuous transmission of matter in front and back.It is thus desirable to setting up the network association of the gas density of each branch
Model.For this purpose, assuming:1)In ventilating system, distinguished and admirable flowing is one-dimensional incompressible flow;2)In ventilation network, each tunnel is last
One spatial mesh size is equal in the gas density of gas density moment back length with certain moment;3)Stand after gas release
Mix with distinguished and admirable, and tunnel intersection each distinguished and admirable in gas be mutually mixed immediately.If according to above-mentioned it is assumed that ventilation network is total
Branch's number is m, and the spatial mesh size number of i-th branch's division is N(i), then in ventilation network, the gas density of any branch i can be pressed
Following formula is calculated, i.e.,
In formula, i=1,2 ..., m;J=1,2 ..., N (i) -1.The arbitrarily gas density boundary condition of branch i is
In formula:Represent tn+1The gas density of last spatial mesh size of moment branch i;Represent tnMoment, branch i fell
The gas density of second spatial mesh size of number;And Qi' represent t respectivelyn+1Moment, the unit bodies that branch i is flowed into by contiguous branch
Product gas amount and unit volume air quantity.Obviously boundary condition will be obtained(5), it is necessary to first solveAnd Qi'.For this purpose, definition
MatrixThe distinguished and admirable relation of branch i and branch k in ventilation network is represented, and
When, remember δ=1, represent the distinguished and admirable inflow branch i of branch k;When, remember δ=0, represent branch k's
Distinguished and admirable outflow branch i or branch k and branch i is unconnected.Then have
In formula,Represent tn+1The air quantity of moment branch k.By formula(6)And formula(7)Substitution formula(5), obtain
As the present invention does not consider the changes of section in same tunnel, so for same tunnel its wind speed u and dispersion coefficient Ex
Do not change with locus, be only the function of time.So A [x can be madej]=a, B [xj]=b, C [xj]=c, D [xj]
=dj, and
In formulaT is represented respectivelynThe airflow velocity and dispersion coefficient of moment branch i.For branch i, if tnThe wind at moment
SpeedAnd dispersion coefficient, it is known that a, b, c can be calculated then with reference to spatial mesh size and time step, while according to each in the branch
The gas density of spatial mesh size, it can be deduced that dj(2≤j≤N-1), and the gas density boundary condition according to branch i, can obtain
Arrive tn+1Gas density x of two step-lengths in the first end in moment tunnel1And xN, the data of 2≤j in tunnel≤N-1 spatial mesh sizes are substituted into
Formula(3)And launch
Equation group is tri-diagonal system, can be solved with chasing method, so as to obtain tn+1The each space steps of moment branch i
Long gas density.Roadway ventilation network equation expression formula is:
Formula(2)Middle RiFor branch's i windages,For the ventilation blower blast of branch i, can be obtained by fan performance curve Equation for Calculating
Arrive;For the gas density in the gas blast of branch i, with roadway branch, the high difference correlation in tunnel.It is time and wind
The function of amount, air quantity qk(k=1,2 ..., b) it is parameter to be asked, and the function of time, i.e., required value is b unknown letter
Number.Formula(2)Complex, it is difficult to directly which to be solved, it is therefore desirable to which is converted according to correlation values solution, this
Invention uses improved Euler's formula(Runge-Kutta method)Method.Order
For convenience of solving, each variable is expressed as matrix form is:A=[ajk], C=[cji], K=[kii], CT=[cik], D=
(D1,D2,…Db)T, q=[q1,q2,…,qk]T(k=1,2 ..., b), wherein kii=Ki, and A=CKCT.So formula(2)It is expressed as
Matrix form is
I.e.
By formula(11)The b independent differential equation can be listed, equation can be solved.Using finite difference calculus, from improved Euler
Formula(Runge-Kutta method)Solve above formula and can reach enough precision.According to finite difference theory, first the time the same time will be pressed
Step-length τ demarcation interval, then tnThe corresponding time at moment is n τ.As the above analysis, formula(11)Actually air quantity is inclined
Lead with regard to time and the function of primary branch air quantity, so q' can be expressed ask=fk(t,q1,q2,…,qb) (k=1,2 ..., b),
Then according to improved Euler's formulaCan obtain
fk(t,q1,q2,…,qb) (k=1,2 ..., expression formula b) should be with regard to time and the letter of each primary branch air quantity
Number, so above formula correctly will be solved, needs its dependent variable to be converted into regard to time and the function of each primary branch air quantity.According to
Derivation above, also includes each branch's fan pressure in the expression formulaWith gas blastTwo variables.The blast of blower fan can
To be obtained by fan performance curve Equation for Calculating, what air pressure characteristics curve reflected is the relation of the blast with air quantity of ventilation blower.
Following fitting of a polynomial generally can use:hF=a0+a1q+a2q2+a3q3+...+anqn, a in formula0, a1..., anFor curve matching system
Number, q are blower fan place branch air quantity.The degree of polynomial of curve requires determination according to computational accuracy, gets 2 herein.I.e.
hF=a0+a1q+a2q2.In order to be fitted the active section of characteristic curve of fan, three point (q on its active section are given1,h1),(q2,
h2),(q3,h3), substitute into above formula, you can obtain the linear equation containing three unknown numbers.Solve this Simultaneous Equations and can be obtained by wind
Every fitting coefficient in pressure characteristic curve equation.
So for any branch i, if known blower fan place branch air quantity Qi, then branch's fan pressureCan be expressed as:
Additionally, gas blast to be solvedThe average density of each branch's gas should first be calculatedI.e.
Due to the mixed gas that the gas is methane and air, so
In formula, MiFor the air gross mass of branch i, kg;ViFor the total air amount of branch i, m3;AiFor the sectional area of branch i,
m2;H is spatial mesh size, m;For tnThe gas density of j-th spatial mesh size of moment branch, kg/m3;For tnMoment branch's jth
The gas density of individual spatial mesh size(Volume fraction);ρmFor methane density, kg/m3;ρaFor atmospheric density, kg/m3.In formula (1)
Speed u is calculated by air quantity q, and is solving formula using improved Euler's method(2)Formula is used during air quantity q again(1)In
Concentration c, so the present invention is by formula (1) and formula(2)Carry out couple solution.t0The value of moment each parameter is, it is known that first solve formula(1)t1
Then the value for solving is substituted into formula by moment gas density(2)Solve t1Moment air quantity, the air quantity for obtaining substitute into formula(1)Solve t2When
The gas density at quarter, the like i.e. be obtained each time step tunnel concentration distribution and air force.
Claims (8)
1. a kind of transient state unstable airflow computational methods of roadway containing high concentration gas, its computation model is by gas dispersion equation
With roadway ventilation network equation
Constitute, wherein u in formula (1)iFor the wind speed of roadway branch i, xiIt is the space length along roadway branch i, t is the time, ciFor
The mash gas concentration of roadway branch i, ExFor gas dispersion coefficient, in formula (2)For branch's air quantity, RiFor branch
I windages,For the ventilation blower blast of branch i,For the gas blast of branch i, cjiFor circuit matrix, cikFor circuit matrix cji
Transposed matrix, wherein gas dispersion equation (1) method for solving is to carry out finite difference initially with Crank-Nicholson methods
Point, then roadway point gas density is obtained using chasing method, while roadway ventilation network equation (2) method for solving is using dragon
Ge-Ku Tafa solve the air quantity (wind speed) of roadway branch, need to carry out couple solution to equation (1) and equation (2) in solution procedure.
2. one kind roadway containing high concentration gas transient state unstable airflow computational methods according to claim 1, its feature exist
In:Described gas dispersion coefficientWherein r is pipe radius, and α is coefficient of frictional resistance, and u is roadway branch
Wind speed.
3. one kind roadway containing high concentration gas transient state unstable airflow computational methods according to claim 1, its feature exist
In:Described branch air quantity qk(k=1,2 ..., b) it is parameter to be asked, and the function of time, i.e., required value is that b is individual unknown
Function, b are the number of independent loop in ventilation network;qkBe pass judgment on roadway unstable airflow index, qk=0 represents distinguished and admirable stopping
Stagnant, qk< 0 represents wind current control, qk> 0 represents that wind direction keeps constant.
4. one kind roadway containing high concentration gas transient state unstable airflow computational methods according to claim 1, its feature exist
In:Described ventilation blower blasta0, a1..., anFor curve matching coefficient;Given three or three with
The blast of fan, air quantity data, solve fitting coefficient by method of least square.
5. one kind roadway containing high concentration gas transient state unstable airflow computational methods according to claim 1, its feature exist
In:Described gas blast For the average density of branch i, For tnMoment point
The gas density of j-th spatial mesh size is propped up, For tnJ-th sky of moment branch
Between step-length mash gas concentration, ρmAnd ρaThe density of respectively methane and air, h is spatial mesh size, LiFor roadway branch i's
Length, spatial mesh size numbers of the N (i) for branch i, ziAnd z (0)i(Li) it is respectively the whole story node absolute altitude of branch i.
6. one kind roadway containing high concentration gas transient state unstable airflow computational methods according to claim 1, its feature exist
In:The method for solving of gas dispersion equation (1) carries out finite difference using Crank-Nicholson methods, obtains three pairs
Angle equation group, distinguished and admirable flowing in ventilating system is one-dimensional incompressible flow, and last is empty for ventilation network Zhong Ge roadways branch
Between step-length equal in the gas density of gas density moment back length with certain moment, gas release after immediately with distinguished and admirable
Mixing, and tunnel intersection each distinguished and admirable in gas be mutually mixed immediately, using chasing method solve tri-diagonal system, obtain well
Lane branch gas density.
7. one kind roadway containing high concentration gas transient state unstable airflow computational methods according to claim 1, its feature exist
In:The method for solving formula (2) of described roadway ventilation network equation (2), order Each variable is expressed as matrix form is:A=[ajk], C=
[cji], K=[kii], CT=[cik], D=(D1,D2,…Db)T, q=[q1,q2,…,qk]T(k=1,2 ..., b), wherein kii=
Ki, and A=CKCT, formula (2) is expressed as matrix form and isSolved from improved Euler's formula (Runge-Kutta method)
The air quantity of each roadway branch, on this basis, according to the basal area of roadway branch, draws the wind speed of roadway branch.
8. one kind roadway containing high concentration gas transient state unstable airflow computational methods according to claim 1, its feature exist
In:Need to carry out couple solution, t to equation (1) and equation (2) in described solution procedure0The value of moment each parameter is, it is known that first ask
Solution formula (1) t1Then the value for solving is substituted into formula (2) and solves t by moment gas density1Moment air quantity, the air quantity for obtaining substitute into formula
(1) t is solved2The gas density at moment, the like i.e. be obtained each time step tunnel concentration distribution and air force.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610951446.1A CN106503363A (en) | 2016-10-27 | 2016-10-27 | A kind of transient state unstable airflow computational methods of roadway containing high concentration gas |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610951446.1A CN106503363A (en) | 2016-10-27 | 2016-10-27 | A kind of transient state unstable airflow computational methods of roadway containing high concentration gas |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106503363A true CN106503363A (en) | 2017-03-15 |
Family
ID=58322249
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610951446.1A Pending CN106503363A (en) | 2016-10-27 | 2016-10-27 | A kind of transient state unstable airflow computational methods of roadway containing high concentration gas |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106503363A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108491663A (en) * | 2018-03-30 | 2018-09-04 | 中煤科工集团重庆研究院有限公司 | Coal mine underground space gas distribution rule calculation method based on mass conservation |
CN109446602A (en) * | 2018-10-15 | 2019-03-08 | 内蒙古科技大学 | A kind of Numerical Experimental Method of surface vertical borehole extraction super high seam gas |
CN109697318A (en) * | 2018-12-24 | 2019-04-30 | 山东蓝光软件有限公司 | The Gas explosion burns inversion method and device resolved based on gas drainage system simulation |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104564120A (en) * | 2014-11-11 | 2015-04-29 | 中国矿业大学 | Operation state control decision making method of mine ventilation system |
CN105697054A (en) * | 2016-04-12 | 2016-06-22 | 安徽理工大学 | Mode of arranging high-gas coal seam high dip tunnel for drawing out gas and gas extraction method |
-
2016
- 2016-10-27 CN CN201610951446.1A patent/CN106503363A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104564120A (en) * | 2014-11-11 | 2015-04-29 | 中国矿业大学 | Operation state control decision making method of mine ventilation system |
CN105697054A (en) * | 2016-04-12 | 2016-06-22 | 安徽理工大学 | Mode of arranging high-gas coal seam high dip tunnel for drawing out gas and gas extraction method |
Non-Patent Citations (4)
Title |
---|
周爱桃 等: "瓦斯风压诱导矿井风流灾变规律研究", 《中国矿业大学学报》 * |
周爱桃: "瓦斯突出冲击气流传播及诱导矿井风流灾变规律研究", 《中国博士学位论文全文数据库工程科技Ⅰ辑》 * |
屠世浩 等: "采空区渗透率演化规律及卸压开采效果研究", 《采矿与安全工程学报》 * |
程卫民: "《矿井通风与安全》", 28 February 2016 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108491663A (en) * | 2018-03-30 | 2018-09-04 | 中煤科工集团重庆研究院有限公司 | Coal mine underground space gas distribution rule calculation method based on mass conservation |
CN108491663B (en) * | 2018-03-30 | 2022-03-22 | 中煤科工集团重庆研究院有限公司 | Coal mine underground space gas distribution rule calculation method based on mass conservation |
CN109446602A (en) * | 2018-10-15 | 2019-03-08 | 内蒙古科技大学 | A kind of Numerical Experimental Method of surface vertical borehole extraction super high seam gas |
CN109697318A (en) * | 2018-12-24 | 2019-04-30 | 山东蓝光软件有限公司 | The Gas explosion burns inversion method and device resolved based on gas drainage system simulation |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Plate | Aerodynamic Characteristics of Atmospheric Boundary Layers. | |
CN106503363A (en) | A kind of transient state unstable airflow computational methods of roadway containing high concentration gas | |
CN106948853A (en) | A kind of global accurate sensors optimum placement method for surveying wind of mine | |
CN106092494A (en) | Drive aircraft pushes away resistance characteristic world conversion method | |
CN109377810A (en) | The monitoring of stope airflow parameter and automatic adjustment experiment porch and system | |
CN104657594A (en) | Determining method for exhaust mass flow of spacecraft under low pressure in cabin | |
Mary | Large eddy simulation of vortex breakdown behind a delta wing | |
Bauer et al. | Variational integrators for anelastic and pseudo-incompressible flows | |
CN104268126B (en) | Multi-modal automatic recognition method for air duct parameters of mine ventilation system | |
JP2005264671A (en) | Dynamic simulation device and program for road tunnel ventilation control | |
CN111914051B (en) | Flight path text description generation method based on geographic grid and neural network | |
Nguyen et al. | A hybridized multiscale discontinuous Galerkin method for compressible flows | |
Sun et al. | Analysis of numerical factors affecting large eddy simulation of pollutant diffusion around buildings | |
McPherson et al. | Ventilation network analysis | |
Ammar et al. | Coupling finite elements and reduced approximation bases | |
Rehm et al. | An efficient large eddy simulation algorithm for computational wind engineering: application to surface pressure computations on a single building | |
Abdelmeguid et al. | A comparison between the parabolic and partially-parabolic solution procedures for three-dimensional turbulent flows around ships' hulls | |
Reed et al. | Application of low speed preconditioning to the compressible Navier-Stokes equations | |
Kot | Numerical modelling of contaminant dispersion around buildings | |
Lombard | Conservative supra-characteristics method for splitting the hyperbolic systems of gasdynamics for real and perfect gases | |
Georgescu et al. | Numerical study of the secondary phase dispersion in a particle-laden flow | |
Sanchez | Investigating overturning high sided vehicles through modeling high Reynolds number incompressible flow around a rectangular cylinder near a plane wall boundary | |
CN103543024B (en) | Performance detecting method and device for porous medium burner | |
Balachandar | Signature of the transition zone in the tomographic results extracted through the Eigenfunctions of the two‐point correlation | |
COCKRELL, JR et al. | Generic hypersonic inlet module analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170315 |
|
WD01 | Invention patent application deemed withdrawn after publication |