CN113312771A - Calculation method and application of limited wind speed of side key smoke exhaust of tunnel - Google Patents

Calculation method and application of limited wind speed of side key smoke exhaust of tunnel Download PDF

Info

Publication number
CN113312771A
CN113312771A CN202110598123.XA CN202110598123A CN113312771A CN 113312771 A CN113312771 A CN 113312771A CN 202110598123 A CN202110598123 A CN 202110598123A CN 113312771 A CN113312771 A CN 113312771A
Authority
CN
China
Prior art keywords
wind speed
tunnel
smoke
fire
limited wind
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
CN202110598123.XA
Other languages
Chinese (zh)
Other versions
CN113312771B (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.)
Wuhan University of Science and Engineering WUSE
Original Assignee
Wuhan University of Science and Engineering WUSE
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 Wuhan University of Science and Engineering WUSE filed Critical Wuhan University of Science and Engineering WUSE
Priority to CN202110598123.XA priority Critical patent/CN113312771B/en
Publication of CN113312771A publication Critical patent/CN113312771A/en
Application granted granted Critical
Publication of CN113312771B publication Critical patent/CN113312771B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Ventilation (AREA)

Abstract

The invention relates to a method for calculating limited wind speed of side key smoke exhaust of a tunnel and application thereof, comprising the following specific steps of: determining a limited wind speed ucDefining, determining influencing factors, establishing a limited wind speed ucA relationship to an influencing factor; determining basic dimension, and establishing limited wind speed ucA dimensional relation with the influencing factor; determining basic physical quantity, obtaining dimensionless items of influence factors according to the pi theorem, and deducing to obtain the limited wind speed ucThe dimensionless calculation formula of (1); obtaining the numerical value of the smoke backflow length L of the fire under different working conditions through numerical simulation, obtaining the limited wind speed through the smoke backflow length, drawing the simulation result into a scatter diagram, and determining the influence factors on the limited wind speed ucThe influence of (a); performing data fitting on the result to obtain the value of each unknown coefficient in the dimensionless relational expression, and further establishing the dimensionless limited wind speed ucThe calculation formula of (2). Has beneficial effectsThe method is simple, is suitable for different tunnels, has great influence on the efficient utilization of the smoke exhaust capacity of the tunnels, and provides guidance for fire smoke control and fire rescue of the side key smoke exhaust tunnel.

Description

Calculation method and application of limited wind speed of side key smoke exhaust of tunnel
Technical Field
The invention relates to the technical field of tunnel fire smoke control, in particular to a calculation method and application of a limited wind speed of side key smoke exhaust of a tunnel.
Background
Fire is one of the greatest risks for tunnel operation. The tunnel has a long, narrow and closed characteristic, so that the tunnel has a unique space structure and fire characteristics compared with a common ground building. Once a fire disaster occurs in the tunnel, smoke cannot be discharged in time, people are very easy to suffocate or be poisoned, vehicles in the tunnel cannot drive away in time to form secondary damage, and serious disastrous accidents of group death and group injury are caused. The limited wind speed is controlled, so that the smoke exhaust capacity of the tunnel can be efficiently utilized, and the smoke is controlled within a certain range, so that people can be evacuated and rescued.
Limited wind speed u in tunnel fire scenecDefinition of (1): when the flue gas countercurrent length is insensitive to the flow velocity of the changed smoke outlet, the corresponding minimum longitudinal induced wind speed is called as the limited wind speed ucThe flue gas countercurrent length refers to the distance between the fire source and the smoke outlet in the tunnel traveling direction, as shown in FIG. 2, and the limited wind speed ucThe proper wind speed of the longitudinal induced wind for ensuring that the flue gas countercurrent length is not increased due to the reduction of the flow speed of the smoke outlet is a control index for controlling the non-diffusion and non-spread of the flue gas. The current research mainly explores a method for calculating the critical wind speed, the critical wind speed is the longitudinal wind speed which just enables the flue gas countercurrent length to be zero, the upstream of a fire source can be ensured to have no flue gas, but the critical wind speed is adopted to control the flue gas to spread, on one hand, the critical wind speed is theoretically larger than the limited wind speed, so the flue gas cannot be discharged to the nearest upstream smoke outlet; on the other hand, the prediction model of the critical wind speed of the fire in the side point type smoke discharging mode is related to the lateral dimensionless smoke discharging amount, the dimensionless fire source heat release rate and the distance between the dimensionless smoke discharging port and the fire source, so that the countercurrent of the smoke is controlled by the longitudinal wind speed predicted by the critical wind speed model and still influenced by the change of the smoke discharging amount.
At present, most of existing limited wind speed calculation methods are used for determining limited wind speed according to the ratio of the smoke backflow length to the tunnel height under the condition that a smoke exhaust port is located at the top of a tunnel, the smoke exhaust port in the tunnel in a side key smoke exhaust mode is located on the side wall of the tunnel, and compared with top smoke exhaust, the stress state and the motion state of smoke are changed at the moment, and a calculation method for limited wind speed of fire disasters of the tunnel in the side key smoke exhaust mode is not available at home and abroad.
Disclosure of Invention
The invention aims to solve the technical problem of providing a calculation method and application of limited wind speed of key smoke discharge fire on the side part of a tunnel, and aims to solve the technical problem.
The technical scheme for solving the technical problems is as follows:
a method for calculating limited wind speed of side key smoke discharge fire of a tunnel comprises the following specific steps:
s1: determining limited wind speed u in tunnel fire scenecEstablishing said limited wind speed ucThe relationship to the influencing factor: f (v, Q, ρ)0,CP,T0,g,HD)＝0；
The influencing factors comprise the heat release rate Q and the air density rho of the fire source0Constant pressure specific heat of air CpAir temperature T0Acceleration of gravity g, tunnel height HDThe number n of physical parameters is 7, and the unit of the heat release rate Q of the fire source is kg.m2/s3Air density ρ0Unit of (b) is kg/m3Constant pressure specific heat of air CpHas the unit of m2/s2K, air temperature T0Has the unit of K and the unit of the gravity acceleration g is m/s2Height H of tunnelDThe unit of (a) is m;
s2: determining basic dimension according to the unit of the influence factor, representing the influence factor by the basic dimension, and establishing the limited wind speed ucA dimensional relationship to the influencing factor;
s3: determining basic dimension according to the unit of the influence factor, representing the influence factor by the basic dimension, and establishing the limited wind speed ucA dimensional relationship to the influencing factor;
in the step S3, the basic dimensions include mass M, time T, length L, and temperature T, and the basic dimension number η is 4;
in the step S3, the limited wind speed ucThe relationship with the influencing factor is as follows:
f(Lt-1,ML2t-3,ML-3,L2t-2T-1,T,Lt-2,L)＝0；
s4: determining the basic physical quantity of the influence factors, obtaining the dimensionless parameter number n-eta-3 of the influence factors according to the pi theorem, and then obtaining the circulation quantity according to the pi theoremSelection principle, selecting variable H directly related to length LDSelecting a variable T directly related to the temperature T0Selecting a variable g directly related to the time t, and selecting a variable ρ directly related to the mass M0The 4 cyclic quantities are combined with all other n-eta physical parameters in turn to form a dimensionless parameter pi as cyclic quantity1、π2、π3Converting the relational expression in the step S3 into a dimensionless relational expression, and then obtaining the limited wind speed ucThe dimensionless calculation formula of (1);
in the step S4, the dimensionless parameter π of the influencing factor1、π2、π3Comprises the following steps:
Figure BDA0003091937170000031
the limited wind speed ucThe dimensionless formula of calculation is:
Figure BDA0003091937170000032
namely, it is
Figure BDA0003091937170000033
uc *Is dimensionless wind speed, Q**Is a dimensionless power;
s5: FDS numerical simulation part: establishing a scaled geometric model of the tunnel, setting a virtual fire source as propane combustion, simulating different fire source heat release rates, setting the smoke outlet flow velocity v to gradually increase from small to large to obtain different smoke countercurrent lengths L fitted by FDS software, and when the trend that the smoke countercurrent length L decreases along with the increase of the smoke outlet flow velocity v becomes slow, the induced wind speed generated by the smoke outlet flow velocity v in the tunnel at the moment is the actual limited wind speed u under the fire conditioncAnd measuring and recording the limited wind speed ucThe number of the influencing factors is that two smoke outlets are arranged on one side wall in the tunnel along the longitudinal direction, and fire is arranged between the two smoke outletsThe fire source is positioned on the central line of the middle lane of the unidirectional three lanes, steady fire or t square fire is adopted, a plurality of temperature measuring points are uniformly arranged in the reduced scale geometric model at intervals along the lengthwise direction of the top wall, and the limited wind speed u fitted by FDS software is extractedcThe cross section of the reduced scale geometric model is provided with a plurality of flow velocity measuring points, and the limited wind speed u fitted by FDS software is extractedc；
S6: and (3) solid simulation: establishing a tunnel reduced scale geometric model, setting different fire conditions in the reduced scale geometric model, including using a burner as a fire source, simulating different fire source heat release rates, setting different exhaust port flow velocities v, controlling different flue gas countercurrent lengths L, changing the longitudinal wind speed under the fire condition, observing the time when the flue gas countercurrent length is not influenced by the exhaust port flow velocity, and measuring the actual limited wind speed u at the timecAnd measuring and recording the limited wind speed ucThe tunnel is internally provided with two smoke outlets along the longitudinal direction, a fire source is arranged between the two smoke outlets and is positioned on the central line of a middle lane of the unidirectional three lanes, a plurality of temperature sensors are uniformly arranged in the reduced scale geometric model along the longitudinal direction of the top wall at intervals, and the limited wind speed u is measuredcThe wind speed meter is arranged on one section of the reduced scale geometric model to measure the limited wind speed uc；
S7: the limited wind speed u under different fire conditions is obtained according to the steps S5 and S6cAnd the numerical values of the influencing factors, respectively drawing the simulation results into scatter diagrams, and if the fitted line graphs are similar, selecting the scatter diagram obtained in the step 5 to fit data;
s8: performing data fitting on the resulting scattergrams of the scattergram obtained in step S5 to obtain unknown coefficients k in the dimensionless calculation formula in step S42Taking the value of (a), and taking the obtained uc *＝k1Q**Drawing a scatter diagram and performing data fitting to obtain the unknown coefficient k in the dimensionless calculation formula in the step S41To establish said limited wind speed u in a dimensionless mannercThe calculation formula of (2):
Figure BDA0003091937170000041
on the basis of the technical scheme, the invention can be further improved as follows.
Further, the application of the method for calculating the limited wind speed of the side key smoke exhaust of the tunnel comprises the following steps:
step 1, according to a calculation formula L of the reverse flow length of the smoke of the tunnel side-direction key smoke discharge fire*＝76(v*)-1.67(Q*)0.56(l*)-0.23And the acquired fire information: q, rho0、Cp、T0、g、HD、L、l、WDCalculating the smoke countercurrent length L, L of the tunnel in which the fire occurs*＝L/HD，l*＝l/HD，Q*＝Q/(HD 3/2g1/2WDT0Cpρ0) Wherein W isDThe tunnel width, v the smoke outlet flow speed and l the distance from the smoke outlet to the fire source;
step 2, if L is larger than L, calculating the limited wind speed u of the side key smoke discharge of the tunnelcAnd the fans of the station halls at the two ends of the accident section of the tunnel simultaneously exhaust air to the accident section to form longitudinal wind, and the flow velocity of the longitudinal wind is greater than ucUntil the temperature sensor or CO alarm in the tunnel detects that the length L of the flue gas counter flow is equal to L, then the flow speed of the longitudinal wind is reduced to uc；
And 3, if the L of the wind discharged by the fans of the station halls at the two ends of the accident section of the tunnel is less than L, reducing the flow velocity of longitudinal wind to be less than ucUntil the fume counterflow length L equals to L, then the flow speed of longitudinal wind is raised to uc。
Compared with the prior art, the invention has the beneficial effects that:
the method is simple, parameters can be set according to the actual conditions of the tunnel, and the method is suitable for different tunnels adopting side key smoke exhaust modes. The method is scientific and effective, applies pi theorem and dimension analysis, and has more theoretical basis. There is no calculation method for the limited wind speed in the side key smoke exhaust mode. The method considers the main influence factors of the heat release rate of the fire source in the process of dimension relation derivation, and the obtained result has innovation and practical engineering significance. The limited wind speed of the tunnel fire has great influence on the efficient utilization of the tunnel smoke discharge capacity, and guidance can be provided for the fire rescue of the tunnel fire under the condition of side key smoke discharge through predicting and calculating the limited wind speed of the tunnel fire.
Drawings
FIG. 1 is a flowchart of a method for calculating a limited wind speed of side-focused smoke evacuation of a tunnel according to the present invention;
FIG. 2 is a schematic view of a limited wind speed in a side-focused smoke evacuation mode according to the present invention;
FIG. 3 is a schematic view of a reduced-size model tunnel built in the present invention;
FIG. 4 is a graph of the limited wind speed versus the heat release rate of the fire source in the present invention;
FIG. 5 is a diagram of the result of a numerical simulation of a limited wind speed according to the present invention;
FIG. 6 is a cross-sectional view of the FDS software at the Velocity measurement point to directly measure induced wind Velocity.
In the drawings, the components represented by the respective reference numerals are listed below:
1. the tunnel comprises a one-way three-lane tunnel, 2, a smoke exhaust port, 3, a temperature measuring point, 4 and a flow velocity measuring point.
Detailed Description
The principles and features of this invention are described in connection with the drawings and the detailed description of the invention, which are set forth below as examples to illustrate the invention and not to limit the scope of the invention.
As shown in fig. 1 to 5, the invention provides a method for calculating limited wind speed of side-focused smoke discharge fire of a tunnel, which comprises the following specific steps:
s1: determining limited wind speed u in tunnel fire scenecDefinition of (1): when the flue gas countercurrent length is insensitive to the flow velocity of the changed smoke outlet, the corresponding minimum longitudinal induced wind speed is called as the limited wind speed ucThe flue gas countercurrent length refers to the distance between the fire source and the smoke outlet in the tunnel advancing direction.
S2: determining limited wind speed u in tunnel fire scenecEstablishing said limited wind speed ucThe relationship to the influencing factor: f (v, Q, ρ)0,CP,T0,g,HD)＝0；
The influencing factors comprise the heat release rate Q and the air density rho of the fire source0Constant pressure specific heat of air CpAir temperature T0Acceleration of gravity g, tunnel height HDThe number n of physical parameters is 7, and the unit of the heat release rate Q of the fire source is kg.m2/s3Air density ρ0Unit of (b) is kg/m3Constant pressure specific heat of air CpHas the unit of m2/s2K, air temperature T0Has the unit of K and the unit of the gravity acceleration g is m/s2Height H of tunnelDThe unit of (d) is m.
S3: determining basic dimension according to the unit of the influence factor, representing the influence factor by the basic dimension, and establishing the limited wind speed ucA dimensional relationship to the influencing factor;
in the step S3, the basic dimensions include mass M, time T, length L, and temperature T, and the basic dimension number η is 4;
in the step S3, the limited wind speed ucThe relationship with the influencing factor is as follows:
f(Lt-1,ML2t-3,ML-3,L2t-2T-1,T,Lt-2,L)＝0；
s4: determining the basic physical quantity of the influence factors, obtaining the dimensionless parameter number n-eta-3 of the influence factors according to the pi theorem, and selecting the variable H directly related to the length L according to the selection principle of the circulation quantity in the pi theoremDSelecting a variable T directly related to the temperature T0Selecting a variable g directly related to the time t, and selecting a variable ρ directly related to the mass M0Using the 4 circulation quantities and all other parameters in other n-eta physical parameters to form dimensionless parameter pi as circulation quantity1、Π2、Π3Converting the relational expression in the step S3 into a dimensionless relational expression, and obtaining the dimensionless relational expressionSaid limited wind speed ucThe dimensionless calculation formula of (1);
in the step S4, non-dimensional parameter Π of the influencing factor1、Π2、Π3Comprises the following steps:
Figure BDA0003091937170000071
the limited wind speed ucThe dimensionless formula of calculation is:
Figure BDA0003091937170000072
namely, it is
Figure BDA0003091937170000081
uc *Is dimensionless wind speed, Q**Is a dimensionless power;
s5: FDS numerical simulation part: and establishing a scaled geometric model of the tunnel, setting a virtual fire source for propane combustion, and simulating different heat release rates of the fire source. The size of the smoke outlet is 0.3m multiplied by 0.1m, the flow velocity v of the smoke outlet is gradually increased from small to large to obtain different smoke countercurrent lengths L fitted by FDS software, and when the trend that the smoke countercurrent length L is reduced along with the increasing of the flow velocity v of the smoke outlet becomes slow, the induced wind speed generated by the flow velocity v of the smoke outlet in the tunnel at the moment is the actual limited wind speed u under the fire conditioncAnd measuring and recording the limited wind speed ucThe scaled geometric model of the tunnel is as follows: 10 reduced size, length x width x height 50m x 1.1m x 0.45 m; two smoke outlets 2 are arranged on one side wall in the unidirectional three-lane tunnel 1 along the longitudinal direction, a fire source is arranged between the two smoke outlets, the fire source is positioned on the central line of a middle lane of the unidirectional three lanes, steady fire or t square fire is adopted, and the size is 0.6m multiplied by 0.2 m; the power of the fire source is 9.5-94.9 kW, the corresponding full-size fire source heat release rate is 3-30 MW, a plurality of temperature measuring points are uniformly arranged in the reduced-scale geometric model at intervals along the longitudinal direction of the top wall, and FDS software is extracted and fitted to obtainLimited wind speed ucThe cross section of the reduced scale geometric model is provided with a plurality of flow velocity measuring points, and the limited wind speed u fitted by FDS software is extractedc；
S5: entity experiment part: the method comprises the steps of establishing a reduced scale geometric model of the tunnel, setting different fire conditions in the reduced scale geometric model, and simulating different heat release rates of the fire source by using a burner as the fire source. Setting different flow velocities v of the smoke outlet with the size of 0.3m multiplied by 0.1m, controlling the counter flow length L of different smoke gases, changing the longitudinal wind speed under the fire condition, observing the moment when the counter flow length of the smoke gases is not influenced by the flow velocity of the smoke outlet, and measuring the actual limited wind speed u at the momentcAnd measuring and recording the limited wind speed ucThe scaled geometric model of the tunnel is as follows: 10 reduced size, length x width x height 10m x 1.1m x 0.45 m; two smoke outlets 2 are arranged on one side wall in the unidirectional three-lane tunnel 1 along the longitudinal direction, a fire source is arranged between the two smoke outlets, the fire source is positioned on the central line of a middle lane of the unidirectional three lanes, and the length, the width and the size of the fire source are multiplied by 15cm and 15 cm; the power of the fire source is 9.5-94.9 kW, and the corresponding full-size heat release rate of the fire source is 3-30 MW. A plurality of temperature sensors are uniformly installed at intervals in the length-wise direction of the top wall in the reduced scale geometric model to measure the limited wind speed ucThe wind speed meter is arranged on one section of the reduced scale geometric model to measure the limited wind speed uc；
S6: the limited wind speed u under different fire conditions is obtained according to the step S5cAnd the numerical values of the influencing factors, and drawing the simulation result into a scatter diagram;
s7: performing non-linear fitting on the resulting scattergrams of all scattergrams to obtain each unknown coefficient k in the dimensionless calculation formula in the step S41、k2To establish said limited wind speed u in a dimensionless mannercThe calculation formula of (2):
Figure BDA0003091937170000091
the above fire dynamics simulation software (FDS) is prior art.
In step S5, the FDS numerical simulation section, the unidirectional three-lane tunnel 1 is 1:10 reduced in size length × width × height 50m × 1.1m × 0.45 m; two smoke outlets 2 are arranged on one side wall in the unidirectional three-lane tunnel 1, a plurality of temperature measuring points 3 are uniformly arranged on the inner top wall of the unidirectional three-lane tunnel at intervals, and a plurality of flow velocity measuring points 4 are arranged on the section of a certain position of the tunnel.
In step S5, in the physical experiment part, the unidirectional three-lane tunnel 1 is 1:10 reduced in size, length × width × height, 10m × 1.1m × 0.45 m; two smoke outlets 2 are arranged on one side wall in a unidirectional three-lane tunnel 1, a plurality of temperature measuring points 3 are uniformly arranged on the inner top wall of the unidirectional three-lane tunnel at intervals, a section of a certain position of the tunnel is provided with a plurality of flow velocity measuring points 4, each temperature measuring point 3 is provided with a k-type thermocouple with the diameter of 1mm and the probe of 0.1m in a manner that can be thought by a person skilled in the art, each thermocouple sensor is connected with a 32-channel paperless recorder through a circuit, the temperature sensor detects the temperature of the corresponding temperature measuring point 3 and sends a corresponding temperature signal to a controller, and the controller receives and stores the corresponding temperature signal. Each flow velocity measuring point 4 is provided with a flow velocity measuring point in a manner that can be thought by those skilled in the art, a hot wire anemometer with the model number (Anemomaster KA23) is adopted to detect the wind velocity at the corresponding flow velocity measuring point 4, and a corresponding wind velocity signal is sent to the controller, and the controller receives and stores the corresponding flow velocity signal. The controller, the temperature sensor and the wind speed sensor all adopt the prior art, and control circuits among the temperature sensors, the wind speed sensors and the controller are also in the prior art. On the basis of the above conditions, values of the heat release rate of the fire source and the flow rate of the smoke outlet are sequentially changed to form different smoke backflow lengths, under the fire condition of each heat release rate of the fire source, the flow rate of the smoke outlet and the smoke backflow length, the longitudinal wind speed is changed to observe and measure to obtain the limited wind speed values under different fire conditions, and the specific simulation result is shown in table 1:
table 1 FDS numerical simulation results:
Figure BDA0003091937170000101
the simulation procedure is exemplified by the first group: when the first group Q is 9.5kW, the changed flow speed of the smoke outlet is 3m/s to 11m/s, and finally the smoke countercurrent length is reduced from 5.95m to 0.75m, and the reduction process finds that when the smoke countercurrent length reaches 1.05m/s, the flow speed of the smoke outlet is increased, and the smoke countercurrent length L is reduced slowly, so that the smoke countercurrent length L becomes insensitive. The flow Velocity of the smoke outlet corresponding to the 1.05m is 8m/s, induced wind speeds are generated at two ends of the accident tunnel due to the fact that the flow Velocity of the smoke outlet is 8m/s, the value of the induced wind speed is the limited wind speed corresponding to Q being 9.5kW, the induced wind speed value of 0.51m/s is directly measured through a Velocity measuring point of the FDS, and the measuring point is shown in figure 6. The simulation process for the remaining groups is similar.
Table 2 results of the physical experiments:
Figure BDA0003091937170000111
because the experiment condition is limited, the actual tunnel is not long enough, only a plurality of groups of experiments with small heat release rate of the fire source can be carried out, and the maximum fan power of the experiment can only lead the flow velocity of the smoke outlet to reach 11 m/s. So the whole experimental data is less. And in winter, the difference between the numerical simulation and the environmental temperature of a laboratory is large, and the final experiment result has errors with the FDS simulation result. The error of the first group of experimental simulation is small, and the experimental trend of the other groups is the same. This experiment verifies the correctness of the results of the FDS simulation. Drawing a scatter diagram according to the numerical simulation results of the tables 1 and 2 to obtain the limited wind speed ucWith respect to the variation relationship of each influencing factor, specifically as shown in fig. 4, the FDS simulation result and the entity experiment result have the same trend, and thus the obtained functions are similar:
FIG. 4 shows a limited wind speed ucThe relation graph of the heat release rate Q of the fire source can be obtained from the graph of FIG. 4cThe heat release rate Q of the fire source is increased and then tends to be stable;
data fitting of FIG. 4 with 0rigin software yields uc *And Q**Has the function relationship of
Figure BDA0003091937170000112
I.e. 0 < Q**When k is less than or equal to 0.52, k21 is ═ 1; when Q is**When > 0.52, k2＝0；
To determine the coefficient k in the above formula1And k1Value of "" A, B, D, Ec *＝k1Q**The calculated values of (A) are plotted in FIG. 5. it can be seen from FIG. 5 that the numerical simulation results all fluctuate around a straight line, k10.34, correlation coefficient 0.99, k1' is 0.4. Shows the dimensionless limited wind speed uc *The numerical calculation formula of (a) and the numerical simulation result of (b) are consistent. Will k1Substituting 0.34 into the above formula can obtain the dimensionless limited wind speed uc *The calculation formula of (2) is as follows:
Figure BDA0003091937170000121
in the formula: u. ofc *: dimensionless restricted wind speed, Q**: dimensionless fire source heat release rate.
By the method, the numerical value of the limited wind speed under different heat release rates of the fire source can be quickly obtained.
The invention has the beneficial effects that: parameters can be set according to the actual conditions of the tunnel, and the tunnel smoke exhausting device is suitable for different tunnels adopting side key smoke exhausting modes. The method is scientific and effective, applies theorems and dimension analysis and has more theoretical basis. Dimensionless limited wind speed u for fire without tunnelc *The method takes the most important influence factors into consideration in the process of dimension relation derivation, and the obtained result has innovation and practical engineering significance. The limited wind speed of the tunnel fire has great influence on the efficient utilization of the tunnel smoke discharge capacity, and guidance can be provided for the fire rescue of the tunnel fire under the condition of side key smoke discharge through predicting and calculating the limited wind speed of the tunnel fire.
An application of a calculation method of a limited wind speed of side key smoke exhaust of a tunnel,
step 1. Tunnel laterology according to patent CN111027176AFormula L for calculating smoke countercurrent length of key smoke discharge fire*＝76(v*)-1.67(Q*)0.56(l*)-0.23And the acquired fire information: q, rho0、Cp、T0、g、HD、L、l、WDCalculating the smoke countercurrent length L, L of the tunnel in which the fire occurs*＝L/HD，l*＝l/HD，Q*＝Q/(HD 3/2g1/2WDT0Cpρ0) Wherein W isDThe tunnel width, v the smoke outlet flow speed and l the distance from the smoke outlet to the fire source;
step 2, if L is larger than L, calculating the limited wind speed u of the side key smoke discharge of the tunnelcAnd the fans of the station halls at the two ends of the accident section of the tunnel simultaneously exhaust air to the accident section to form longitudinal wind, and the flow velocity of the longitudinal wind is greater than ucUntil the temperature sensor or CO alarm in the tunnel detects that the length L of the flue gas counter flow is equal to L, then the flow speed of the longitudinal wind is reduced to uc；
And 3, if the L of the wind discharged by the fans of the station halls at the two ends of the accident section of the tunnel is less than L, reducing the flow velocity of longitudinal wind to be less than ucUntil the fume counterflow length L equals to L, then the flow speed of longitudinal wind is raised to uc。
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims (2)

1. A method for calculating limited wind speed of side key smoke discharge fire of a tunnel comprises the following specific steps:
s1: determining limited wind speed u in tunnel fire scenecEstablishing said limited wind speed ucThe relationship to the influencing factor: f (v, Q, ρ)0,CP,T0,g,HD)＝0；
The influencing factors comprise the heat release rate Q and the air density rho of the fire source0Air-fixingSpecific heat under pressure CpAir temperature T0Acceleration of gravity g, tunnel height HDThe number n of physical parameters is 7, and the unit of the heat release rate Q of the fire source is kg.m2/s3Air density ρ0Unit of (b) is kg/m3Constant pressure specific heat of air CpHas the unit of m2/s2K, air temperature T0Has the unit of K and the unit of the gravity acceleration g is m/s2Height H of tunnelDThe unit of (a) is m;
s2: determining basic dimension according to the unit of the influence factor, representing the influence factor by the basic dimension, and establishing the limited wind speed ucA dimensional relationship to the influencing factor;
s3: determining basic dimension according to the unit of the influence factor, representing the influence factor by the basic dimension, and establishing the limited wind speed ucA dimensional relationship to the influencing factor;
in step S3, the basic dimensions include mass M, time T, length L, temperature T, and basic dimension number
Figure FDA0003091937160000011
Is 4;
in the step S3, the limited wind speed ucThe relationship with the influencing factor is as follows:
f(Lt-1,ML2t-3,ML-3,L2t-2T-1,T,Lt-2,L)＝0；
s4: determining the basic physical quantity of the influence factors, and obtaining the dimensionless parameter number of the influence factors according to the pi theorem
Figure FDA0003091937160000012
Then according to the selection principle of circulation quantity in pi theorem, selecting variable H directly related to length LDSelecting a variable T directly related to the temperature T0Selecting a variable g directly related to the time t, and selecting a variable ρ directly related to the mass M0Using the 4 circulation amounts and other n-All other parameters in eta physical parameters are combined in sequence to form dimensionless parameter pi1、Π2、Π3Converting the relational expression in the step S3 into a dimensionless relational expression, and then obtaining the limited wind speed ucThe dimensionless calculation formula of (1);
in the step S4, non-dimensional parameter Π of the influencing factor1、Π2、Π3Comprises the following steps:
Figure FDA0003091937160000021
the limited wind speed ucThe dimensionless formula of calculation is:
Figure FDA0003091937160000022
namely, it is
Figure FDA0003091937160000023
uc *Is dimensionless wind speed, Q**Is a dimensionless power;
s5: FDS numerical simulation part: establishing a scaled geometric model of the tunnel, setting a virtual fire source as propane combustion, simulating different fire source heat release rates, setting the smoke outlet flow velocity v to gradually increase from small to large to obtain different smoke countercurrent lengths L fitted by FDS software, and when the trend that the smoke countercurrent length L decreases along with the increase of the smoke outlet flow velocity v becomes slow, the induced wind speed generated by the smoke outlet flow velocity v in the tunnel at the moment is the actual limited wind speed u under the fire conditioncAnd measuring and recording the limited wind speed ucThe method comprises the following steps of (1) extracting values of influence factors, arranging two smoke outlets along the longitudinal direction on one side wall in a tunnel, arranging a fire source between the two smoke outlets, enabling the fire source to be located on a central line of a middle lane of a unidirectional three-lane, adopting steady fire or t square fire, uniformly arranging a plurality of temperature measuring points at intervals along the longitudinal direction of a top wall in a reduced scale geometric model, and extracting limited wind speed fitted by FDS softwareucThe cross section of the reduced scale geometric model is provided with a plurality of flow velocity measuring points, and the limited wind speed u fitted by FDS software is extractedc；
S6: and (3) solid simulation: establishing a tunnel reduced scale geometric model, setting different fire conditions in the reduced scale geometric model, including using a burner as a fire source, simulating different fire source heat release rates, setting different exhaust port flow velocities v, controlling different flue gas countercurrent lengths L, changing the longitudinal wind speed under the fire condition, observing the time when the flue gas countercurrent length is not influenced by the exhaust port flow velocity, and measuring the actual limited wind speed u at the timecAnd measuring and recording the limited wind speed ucThe tunnel is internally provided with two smoke outlets along the longitudinal direction, a fire source is arranged between the two smoke outlets and is positioned on the central line of a middle lane of the unidirectional three lanes, a plurality of temperature sensors are uniformly arranged in the reduced scale geometric model along the longitudinal direction of the top wall at intervals, and the limited wind speed u is measuredcThe wind speed meter is arranged on one section of the reduced scale geometric model to measure the limited wind speed uc；
S7: the limited wind speed u under different fire conditions is obtained according to the steps S5 and S6cAnd the numerical values of the influencing factors, respectively drawing the simulation results into scatter diagrams, and if the fitted line graphs are similar, selecting the scatter diagram obtained in the step 5 to fit data;
s8: performing data fitting on the result scatter points of all the scatter diagrams obtained in the step (5) to obtain an unknown coefficient k in the dimensionless calculation formula in the step (S4)2Taking the value of (a), and taking the obtained uc *＝k1Q**Drawing a scatter diagram and performing data fitting to obtain the unknown coefficient k in the dimensionless calculation formula in the step S41To establish said limited wind speed u in a dimensionless mannercThe calculation formula of (2):
Figure FDA0003091937160000031
2. the application of the method for calculating the limited wind speed of side-focused smoke exhaust of the tunnel according to claim 1 comprises the following steps:
step 1, according to a calculation formula L of the reverse flow length of the smoke of the tunnel side-direction key smoke discharge fire*＝76(v*)-1.67(Q*)0.56(l*)-0.23And the acquired fire information: q, rho0、Cp、T0、g、HD、L、l、WDCalculating the smoke countercurrent length L, L of the tunnel in which the fire occurs*＝L/HD，l*＝l/HD，Q*＝Q/(HD 3/2g1/2WDT0Cpρ0) Wherein W isDThe tunnel width, v the smoke outlet flow speed and l the distance from the smoke outlet to the fire source;
step 2, if L is larger than L, calculating the limited wind speed u of the side key smoke discharge of the tunnelcAnd the fans of the station halls at the two ends of the accident section of the tunnel simultaneously exhaust air to the accident section to form longitudinal wind, and the flow velocity of the longitudinal wind is greater than ucUntil the temperature sensor or CO alarm in the tunnel detects that the length L of the flue gas counter flow is equal to L, then the flow speed of the longitudinal wind is reduced to uc；
And 3, if the L of the wind discharged by the fans of the station halls at the two ends of the accident section of the tunnel is less than L, reducing the flow velocity of longitudinal wind to be less than ucUntil the fume counterflow length L equals to L, then the flow speed of longitudinal wind is raised to uc。
CN202110598123.XA 2021-05-31 2021-05-31 Calculation method and application of limited wind speed of side key smoke exhaust of tunnel Active CN113312771B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110598123.XA CN113312771B (en) 2021-05-31 2021-05-31 Calculation method and application of limited wind speed of side key smoke exhaust of tunnel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110598123.XA CN113312771B (en) 2021-05-31 2021-05-31 Calculation method and application of limited wind speed of side key smoke exhaust of tunnel

Publications (2)

Publication Number Publication Date
CN113312771A true CN113312771A (en) 2021-08-27
CN113312771B CN113312771B (en) 2022-05-20

Family

ID=77376276

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110598123.XA Active CN113312771B (en) 2021-05-31 2021-05-31 Calculation method and application of limited wind speed of side key smoke exhaust of tunnel

Country Status (1)

Country Link
CN (1) CN113312771B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110275302A1 (en) * 2008-10-24 2011-11-10 Mosen Limited Improved tunnel ventilation device
US20130137356A1 (en) * 2010-07-27 2013-05-30 Ivor Pavetic Method And System For Tunnel Ventilation In Normal Conditions And In Conditions Of Fire
CN108022496A (en) * 2018-01-05 2018-05-11 清华大学 Tunnel fire test device and test method of a kind of top containing crossbeam
CN110909474A (en) * 2019-11-27 2020-03-24 武汉科技大学 Method for determining effective air volume of air supply of vertical shaft of tunnel between subway sections
CN111027176A (en) * 2019-11-07 2020-04-17 武汉科技大学 Method for calculating flue gas countercurrent length of tunnel side-direction key smoke discharge fire
CN111075493A (en) * 2019-12-05 2020-04-28 西南交通大学 Design method for controlling wind speed by fire smoke of plateau tunnel
CN111462609A (en) * 2020-04-21 2020-07-28 合肥工业大学 Fire burning and ventilation control system of spiral tunnel group

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110275302A1 (en) * 2008-10-24 2011-11-10 Mosen Limited Improved tunnel ventilation device
US20130137356A1 (en) * 2010-07-27 2013-05-30 Ivor Pavetic Method And System For Tunnel Ventilation In Normal Conditions And In Conditions Of Fire
CN108022496A (en) * 2018-01-05 2018-05-11 清华大学 Tunnel fire test device and test method of a kind of top containing crossbeam
CN111027176A (en) * 2019-11-07 2020-04-17 武汉科技大学 Method for calculating flue gas countercurrent length of tunnel side-direction key smoke discharge fire
CN110909474A (en) * 2019-11-27 2020-03-24 武汉科技大学 Method for determining effective air volume of air supply of vertical shaft of tunnel between subway sections
CN111075493A (en) * 2019-12-05 2020-04-28 西南交通大学 Design method for controlling wind speed by fire smoke of plateau tunnel
CN111462609A (en) * 2020-04-21 2020-07-28 合肥工业大学 Fire burning and ventilation control system of spiral tunnel group

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
XUEPENG JIANG等: "Study on induced airflow velocity of point smoke extraction in road tunnel fires", 《TUNNELLING AND UNDERGROUND SPACE TECHNOLOGY》 *
姜学鹏等: "侧部点式排烟隧道火灾临界风速研究", 《中国安全科学学报》 *
程婧园: "隧道火灾临界风速计算公式的比较分析", 《消防科学与技术》 *

Also Published As

Publication number Publication date
CN113312771B (en) 2022-05-20

Similar Documents

Publication Publication Date Title
Chen Computational fluid dynamics for HVAC: successes and failures
CN111027176B (en) Calculation method for smoke countercurrent length of tunnel lateral key smoke discharging fire
Dullenkopf et al. The effect of incident wake conditions on the mean heat transfer of an airfoil
CN109583585B (en) Construction method of power station boiler wall temperature prediction neural network model
Zhang et al. An inverse method based on CFD to quantify the temporal release rate of a continuously released pollutant source
Heiselberg et al. Experimental and CFD evidence of multiple solutions in a naturally ventilated building.
CN103728339B (en) A kind of real-time identification method for average heat resistance of heat-exchange equipment on thermal power boiler side
CN109259291A (en) For predicting the numerical method of pipe tobacco Mass and heat transfer rule in cylinder cut tobacco dryer
CN113312771B (en) Calculation method and application of limited wind speed of side key smoke exhaust of tunnel
Pekař et al. Circuit heating plant model with internal delays
Zhang et al. An online model for predicting dynamic response of a packed liquid desiccant dehumidifier
JP3746700B2 (en) Control method of concentrate burner
Zhao et al. Wind load of low-rise building based on fluent equilibrium atmospheric boundary layer
Bottcher et al. Experimental analysis of wind ventilation of poultry buildings
CN105973319A (en) Method for computing hydraulic characteristic of sewage system of control rod driving mechanism
CN113868952A (en) Online detection method for temperature field and ring thickness distribution in rotary kiln
Lu et al. Indoor aerosol particle deposition and distribution: numerical analysis for a one-zone ventilation system
JPH06288390A (en) Operation control method of turbo blower or compressor piping system
CN113482693B (en) Tunnel smoke-proof partition division simulation device and method, calculation model and application
CN114519288A (en) Icing power transmission line flow field finite element analysis method based on S-A turbulence model
Silva et al. Evatuation system of exhaust fans used on ventilation system in commercial broiler house
CN217717587U (en) Tunnel flue gas mass flow rate proportion measuring device
CN115789854B (en) Method and system for measuring and adjusting effective coefficient of air flow organization of air-conditioning room
CN115901855A (en) Tunnel smoke mass flow rate proportion testing method
JPH03260798A (en) Smoke behavior analyzing device

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