Method for calculating head loss coefficient and branch pipe flow distribution of multi-point water intake

Abstract

The present invention relates to a method for calculating a head loss coefficient and branch pipe flow distribution of a multi-point water intake. The method includes the following steps: inputting and determining a structure and size information of the multi-point water intake; determining a head loss coefficient; determining a head loss coefficient form of a confluence port; pre-estimating or adjusting a trial flow; performing trial calculation on the head loss coefficient of each branch pipe; comparing the head losses calculated by two flow paths; and carrying out overall calculation. According to the present invention, the flow of each branch pipe is adjusted and the head loss coefficient is repeatedly calculated by utilizing iterative calculation of the head loss coefficients of the flows in two directions, so that the head loss and the flow distribution proportion can be finally and accurately calculated.

Claims

1. A method for calculating a head loss coefficient and branch pipe flow distribution of a multi-point water intake, for designing a structure of a multi-point water intake facility which comprises a main pipe, provided with a plurality of confluence ports, buried in a water bottom seabed, branch pipes each vertically connected to a corresponding one of the confluence ports, and water-absorbing heads each immersed in water and arranged at a top end of each branch pipe, the method comprising the following steps: step 1: inputting and determining a structure and size information of the main pipe, the plurality of confluence ports and the branch pipes, wherein, the structure and size information comprise cross sectional area A of the main pipe of the multi-point water intake, a shape of each branch pipe, a length l.sub.i of the branch pipes and cross sectional area A.sub.i of the branch pipes, a smooth radius r of a confluence port of the branch pipe and the main pipe, a spacing L between the branch pipes, a flow Reynolds number Re and a rough height of a pipe material, wherein i is a serial number of the branch pipes arranged from upstream to downstream and is equal to 1, 2 . . . , n; step 2: determining local pipeline resistance coefficients of the water-absorbing heads and on-way friction resistance coefficients of the main pipe and the branch pipes, wherein, determining a local resistance coefficient K.sub.in of the water-absorbing head and an on-way friction resistance coefficient K.sub.f of the main pipe and the branch pipes; step 3: determining a head loss coefficient at the confluence port according to the structure and size information and the local pipeline resistance coefficients of the water-absorbing heads and on-way friction resistance coefficients of the main pipe and the branch pipes, wherein, determining a calculation form of a head loss coefficient K.sub..fwdarw.c of a confluence pipe according to shapes of the branch pipes and the main pipe, a flow area ratio, a confluence port angle and the smooth radius and according to a Gardel formula, wherein a loss coefficient of an i-th confluence port comprises a flowing loss coefficient of the branch pipes K.sub.i, b.fwdarw.c of an i-th branch pipe and a penetrating flow loss coefficient K.sub.i, t.fwdarw.c of flowing through the i-th confluence port from the upstream main pipe, the form of the Gardel formula is as follows: $K_{t .fwdarw. c} = 0.03 {(1 - R_{Q})}^{2} - R_{Q}^{2} {1 + (1.62 - r^{0.5}) [\frac{\cos (-)}{R_{A}} - 1] - 0.38 (1 - R_{A})} + (2 - R_{A}) (1 - R_{Q}) R_{Q}$ $K_{b .fwdarw. c} = - 0.92 {(1 -)}^{2} - R_{Q}^{2} {(1.2 - r^{0.5}) [\frac{\cos (-)}{R_{A}} - 1] + 0.8 [1 - R_{A}^{- 2}] - (R_{A}^{- 1} - 1) \cos (-)} + (2 - R_{A}) (1 - R_{Q}) R_{Q}$ wherein K.sub.t.fwdarw.c is a penetrating flow loss coefficient of the main pipe, K.sub.b.fwdarw.c is a flowing loss coefficient of the branch pipes, the subscript t represents upstream main pipe flowing, the subscript c represents main pipe downstream, is an included angle between the branch pipe and the main pipe upstream, =/2 at the T confluence port, a dimensionless number R.sub.Q is a flow ratio of a branch flow to a downstream main flow, R.sub.A is an area ratio of the branch pipe to the main pipe, and r is a smooth radius r of a confluence port of the branch pipe and the main pipe when considering a smooth angle confluence port situation; step 4: pre-estimating or adjusting trial flows and upstream accumulative flows of the branch pipes, wherein, from i=2, pre-estimating a trial dimensionless flow Q.sub.1 of the branch pipe i and an upstream accumulative flow ${.Math.}_{j = 1}^{i - 1} j_{}$ of the branch pipe i according to experience or existing data, or adjusting a trial flow Q.sub.i of the branch pipe i and an upstream accumulative flow ${.Math.}_{j = 1}^{i - 1} j_{}$ of the branch pipe i according to the trial result, wherein i is a value from 2 to n; step 5: calculating a head loss K.sub.i at a downstream confluence port of the branch pipe i according to the head loss coefficient at the confluence port, and the trial flows and upstream accumulative flows of the branch pipes, wherein, calculating a K.sub.i of the confluence port from two flow paths respectively according to the branch pipe flow Q.sub.i and the upstream accumulative flow ${.Math.}_{j = 1}^{i - 1} j_{}$ obtained by trial calculation: a first flow path is an upstream penetrating flow path p.sub.i-1: ${K_{1}^{} .Math.}_{p_{.Math. i - 1}} = K_{i - 1}^{} + K_{i - 1 ~ i, f} {({.Math.}_{j = 1}^{i - 1} j_{})}^{2} + K_{i, t .fwdarw. c} {({.Math.}_{j = 1}^{i} Q_{j}^{})}^{2}$ wherein K.sub.i.sup.|p.sup.i-1 is a trial-calculation head loss of a total confluence flow at the i point of the i-th T-shaped confluence port downstream in the main pipe from the first branch pipe p.sub.1 to the (i1)-th branch pipe p.sub.i-1 from the main pipe upstream penetrating flow to the downstream i point, K.sub.i-1 is a trial head loss of the position i1, and K.sub.i-1i,f is an on-way friction resistance loss of the main pipe from the position i1 to the position i section, a second flow path is a branch pipe flowing flow path p.sub.i: ${K_{i}^{} .Math.}_{p_{i}} = K_{i, in} {(\frac{A}{A_{i}})}^{2} i_{2} + K_{i, f} {(\frac{A}{A_{i}})}^{2} Q_{i}^{2} + K_{i, b .fwdarw. c} {({.Math.}_{j = 1}^{i} Q_{j}^{})}^{2}$ wherein K.sub.i.sup.|p.sup.i is a trial-calculation head loss at the i point of the i-th T-shaped confluence port downstream in the main pipe from the i-th branch pipe to the water flow path p.sub.i; step 6: comparing the head losses calculated for the first flow path and the second flow path, judging whether the head losses are equal, if yes, proceeding to a next branch pipe, and if no, returning to the step 4 and adjusting the trial-calculation flow, so as to obtain the head losses and flow distribution proportions of the branch pipes, wherein, every time performing one-time adjustment on the branch pipe flow Q.sub.i and trial calculation on the head loss K.sub.i of the branch pipe with a serial number i, that is, comparing the head losses calculated by the two flow paths p.sub.i-1 and p.sub.i for one time, judging whether they are equal, setting a threshold error which is used to judge whether the head losses are equal according to an expected flow precision requirement:
|K.sub.i.sup.|p.sup.iK.sub.i.sup.|p.sup.i-1|<error, if yes, proceeding to a next branch pipe, and if no, returning to the step 4 and adjusting the trial-calculation flow, so that the cyclic iterative calculation is performed from i=2 to i=n, that is, the calculation of the most downstream confluence port is completed; and step 7: calculating an total head loss coefficient of the multi-point water intake facility according to the head losses and flow distribution proportions of the branch pipes, so as to obtain a total head loss of the multi-point water intake facility for designing the structure of the multi-point water intake facility, when calculating the most downstream main pipe position n, that is, i=n, obtaining a loss coefficient K.sub.n of a composite pipe here, then an overall head loss coefficient K and an overall head loss H of a water intake system are respectively: $K = K_{n} / {(.Math. i)}^{2}$ $H = K \frac{U^{2}}{2 g}$ wherein in the formula, $U = A$ is an average flow velocity of the most downstream main pipe, and g is a gravitational acceleration.

2. The method according to claim 1, wherein when on-way diameters of the main pipe are different, the head loss K in the above steps is corrected: $K_{correct} = {K_{origin} (\frac{A_{d}}{A_{u}})}^{2}$ A.sub.d being area of a downstream main pipe, A.sub.u being area of an upstream main pipe, K.sub.origin being an uncorrected head loss, and K.sub.correct being a corrected head loss.

Description

DETAILED BRIEF DESCRIPTION OF THE DRAWINGS

(1) The present invention is further described with reference to the accompanying drawings and embodiments.

(2) FIG. 1 is a distribution schematic diagram of a structure and serial number of a series multi-point water intake, a flow path and a head loss of the method of the present invention applied in an embodiment I of the present invention;

(3) FIG. 2 is a calculation schematic diagram of the method according to an embodiment I of the present invention;

(4) FIG. 3 is a schematic diagram of an application example: a 3-point water intake of a certain power plant, in an embodiment I of the present invention;

(5) FIG. 4 is a schematic diagram of an application example: a calculation process of three flow paths p and three positions, in an embodiment I of the present invention; and

(6) FIG. 5 is a structural schematic diagram of a model size of a 6-point position model test in an embodiment I of the present invention.

DESCRIPTION OF THE EMBODIMENTS

Embodiment I

(7) This embodiment is a method for calculating a head loss coefficient and branch pipe flow distribution of a multi-point water intake. A multi-point water intake facility calculated by the method includes: a main pipe 001 buried in a water bottom seabed, the main pipe is provided with a plurality of confluence ports 002 (generally are T-shaped confluence ports), branch pipes each vertically connected to a corresponding one of the confluence ports 003, and a water-absorbing head 004 immersed in water is arranged at a top end of each branch pipe, as shown in FIG. 1. In FIG. 1, the number of the multi-point water intakes is generalized as four, and the case where the cross sectional area A of the main pipe on way is equal is taken as an example. In practice, there may be more water intakes, such as six and seven. Of course, there may be fewer water intakes, such as two and three. The purpose of setting a plurality of water intakes is to disperse water intake.

(8) The specific operation process and principle of the method in the embodiment are as follows, and the operation flow is shown in FIG. 2: firstly, structural information of the multi-point water intake is determined: inputting cross sectional area A of the multi-point water intake, a shape, length l.sub.i and cross sectional area of a branch pipe, a spacing L between the branch pipes, a flow Reynolds number Re and a rough height of a pipe material, wherein the A A.sub.i may be on-way changeable or on-way equal. If the on-way area of the main pipes is not equal, the area A.sub.d of the downstream main pipe and the area A.sub.u of the upstream main pipe of each section are determined, the square of their ratio is calculated, and later the head loss coefficient is standardized, that is,

(9) 0 $K = {K^{} (\frac{A_{d}}{A_{u}})}^{2} .$

(10) A local resistance and an on-way resistance coefficient of a single structure are determined: determining a local resistance coefficient K.sub.in of the water-absorbing head and an on-way friction resistance coefficient K.sub.f of the main pipe and the branch pipes.

(11) The forms of the head loss coefficients K.sub.i,b.fwdarw.c and K.sub.i,t.fwdarw.c of the T-shaped confluence port are determined according to the cross section area ratio and the shape type of the branch pipes and the main pipe. K.sub.i,b.fwdarw.c and K.sub.i,t.fwdarw.c, are generally determined by an area ratio of branch/main pipes, a branch/main flow ratio and a critical distance. In application practice, for the determined multi-point water intake system, parameters such as the area ratio of the branch/main pipelines and the distance are determined constants; therefore, values of K.sub.i,b.fwdarw.c and K.sub.i,t.fwdarw.c are only functions of the flow ratio. The loss coefficients of the T-shaped confluence port are K.sub.i,b.fwdarw.c and K.sub.i,t.fwdarw.c. Cyclic iterative calculation: assuming the dimensionless inflow Q.sub.1 and Q.sub.2 of the branch pipes, the head loss coefficient K.sub.{circle around (2)} of an arrival position {circle around (2)} of a flow path p.sub.1 and a flow path p.sub.2 is calculated respectively based on this assumption. Generally speaking, p.sub.1 calculation includes the following parts: K.sub.{circle around (2)}|.sub.p1=head loss coefficient K.sub.{circle around (1)}|.sub.p1 at the position {circle around (1)} +on-way friction resistance loss coefficient K.sub.12,fQ.sub.1.sup.2 of the section 1-2 of the main pipe +overflow head coefficient K.sub.2,t.fwdarw.c(Q.sub.1+Q.sub.2).sup.2 of confluence port 2
wherein K.sub.{circle around (1)}|.sub.p1=local resistance coefficient

(12) $K_{1, in} {(\frac{A}{A_{1}})}^{2} 1_{2}$
of water-absorbing head 1 +on-way friction resistance loss coefficient

(13) $K_{1, f} {(\frac{A}{A_{1}})}^{2} 1_{2}$ of the vertical pipe 1 +penetrating flow head loss coefficient K.sub.1,t-tQ.sub.1.sup.2 of confluence port 1
p.sub.2 calculation includes the following parts: K.sub.{circle around (2)}|.sub.p2=local resistance coefficient

(14) $K_{2, in} {(\frac{A}{A_{2}})}^{2} 2_{2}$ of inlet 2 +on-way friction resistance coefficient

(15) $K_{2, f} {(\frac{A}{A_{2}})}^{2} 2_{2}$ of vertical pipe 2 +branch flow confluence head loss coefficient K.sub.2,b.fwdarw.t(Q.sub.1+Q.sub.2).sup.2 of confluence port 2

(16) Then whether the loss coefficients calculated by the two paths are equal (K.sub.{circle around (2)}.sup.|p1=K.sub.{circle around (2)}.sup.|p2) is judged, or whether a difference value between the two meets a set convergence threshold (|K.sub.{circle around (2)}.sup.|p1K.sub.{circle around (2)}.sup.|p2|<error) is judged, wherein error is the convergence threshold with a value of 0.01 Q. If the loss coefficients are not equal or the difference between the two is greater than the set threshold, flow distribution is adjusted, K custom character .sup.|p.sup.i-1 and K.sup.|p.sup.i are calculated again until the loss coefficients are equal, and the loss coefficient K.sub.2 at the position {circle around (2)} is calculated. flow distribution may be adjusted by a gradient descent method.

(17) Incremental calculation is continued along the pipeline serial number according to the above calculation process, the head loss coefficient of a custom character position is calculated by the loss coefficient K.sub.i-1 of the i1 position calculated in the previous step and the assumed flow of the branch pipe i: calculating the head loss coefficient K.sup.|p.sup.i-1 of an arrival position of a flow path K.sup.p.sup.i-1 (after confluence of the branch pipes 1, 2, . . . , i1), and calculating a loss coefficient K custom character .sup.|p.sup.i of an arrival position of a flow path p.sub.i (entering from the branch pipe i), wherein K.sup.|p.sup.i-1=a head loss coefficient K.sub.i-1 of the position i1, +on-way friction resistance head loss

(18) $K_{i - 1 i, f} {({.Math.}_{j = 1}^{i - 1} j_{})}^{2}$ in the section i1 to i of the main pipe +overflow head loss coefficient

(19) $K_{i, t .fwdarw. c} {({.Math.}_{j = 1}^{i - 1} j_{})}^{2}$ of confluence port i K custom character .sup.|p.sup.i=local resistance coefficient

(20) $K_{i, in} {(\frac{A}{A_{1}})}^{2} i_{2}$ of inlet i +on-way friction

(21) $K_{i, f} {(\frac{A}{A_{i}})}^{2} Q_{i}^{2}$ resistance coefficient in vertical pipe i +branch flow confluence head loss coefficient

(22) $K_{i, b .fwdarw. c} {({.Math.}_{j = 1}^{i} Q_{i}^{})}^{2}$ of confluence port i

(23) If comparison between K custom character .sup.|p.sup.i and K.sup.|p.sup.i-1 does not meet a set convergence criterion, a flow Q.sub.i of the pipe i is assumed again until the convergence criterion is met, and the loss coefficient K.sub.i at the position is obtained. 1) When the most downstream main pipe position is calculated, that is, i=n, a loss coefficient K.sub.n of a composite pipe here may be obtained. 2) Then an overall head loss coefficient and an overall head loss of a water intake system are respectively:

(24) 0 $K = K_{n} / {(.Math. Q_{i})}^{2}$ $H = K \frac{U^{2}}{2 g}$
In the formula,

(25) $U = \frac{Q}{A} =$
is an average flow velocity of the most downstream main pipe.

(26) The multi-point water intakes arranged in series may be regarded as a continuous confluence pipeline. For a continuous confluence pipe with a constant main pipe on-way section, with the gradual confluence of the water flow, the total flow velocity of the main pipe gradually increases, but static pressure gradually decreases. For a short composite pipe, the static pressure decreasing towards the downstream end along the way will increase the flow of the branch pipe at the downstream end. Finally, the confluence quantity of the branch pipe at the near end is large and the confluence quantity of the branch pipe at the tail end is minimal.

(27) The core of solving the flow distribution and the total loss coefficient of the series multi-point water intake is to clarify the form of a constraint function K ( )the expression form of the local head loss and the on-way loss. A method for determining a head loss coefficient in a constraint formula is as follows:

(28) K.sub.in is a local head loss coefficient of the water flow entering from a water-absorbing head to the branch pipe, and the general expression is:
K.sub.i,in=.sub.i
Wherein, is a local resistance coefficient of an inlet of the water-absorbing head. In general, the water-absorbing heads of the multi-point water intakes have the consistent shape and size and have the same local head loss coefficient which may be selected as a constant value by inquiring a local resistance coefficient table (a water inlet with a top cap+ a gradient shrinkage section) according to the structural size form of the water-absorbing heads.

(29) K.sub.f represents an on-way head loss coefficient, and the general expression is: Branch pipe:

(30) $K_{i, f} =_{i} \frac{l_{i}}{d_{i}}$ Main pipe:

(31) $K_{i - 1 ~ i, f} = \frac{L_{i - 1 ~ i}}{D}$

(32) .sub.i is an on-way friction resistance coefficient of the pipeline, which is related to a Reynolds number Re and relative roughness of a vertical pipe. It may be calculated by an empirical formula according to the flow state:

(33) $K_{p} = \frac{Q_{t} K_{t .fwdarw. c} + Q_{b} K_{b .fwdarw. c}}{Q_{c}}$

(34) In the actual engineering, a spacing L.sub.i-1i between the branch pipes is generally about 3-5 D. At this time, the local head loss of the confluence is absolutely dominant compared with the main pipe on-way loss, and an on-way head loss of a main pipe section between the adjacent vertical pipes may be completely covered by h.sub.i,.fwdarw.c. Therefore, the main pipe on-way head loss h.sub.i-1i,f0 is ignored.

(35) K.sub..fwdarw.c represents a branch flow confluence loss at the confluence port and a main pipe penetrating flow loss coefficient. The local head loss coefficient at the confluence port is closed related to a flow ratio, an area ratio and a smooth radius r, that is, K.sub.i.fwdarw.c=function(R.sub.Ai, R.sub.Qi, r); therefore, the constraint formula function is a non-linear form. In the actual application, R.sub.Ai and r have been determined, so f is only a function of R.sub.Qi.

(36) for the first branch pipe of the multi-point water intake of an upstream seal (Q.sub.0=0) a local loss coefficient at the confluence port may be expressed as K.sub.1,b.fwdarw.c=function(R.sub.A1,R.sub.Q1=1,r). However, in most cases, the head loss coefficient of the bending confluence pipeline needs to be necessarily corrected.

APPLICATION EXAMPLE

(37) a certain power plant will adopt a 3-point water intake in the form of a buried culvert under a seabed to take water and cool. The structure of a water intake system is as shown in FIG. 3. FIG. 4 shows the meaning of each operational symbol of the structure of the water intake system, wherein the confluence port custom character may be simplified as a special T-shaped confluence port. At the beginning of the design of the water intake, it is necessary to estimate inflow distribution of three branch pipes and the overall head loss coefficient. The local resistance coefficient of the water-absorbing head (water inlet) of the water intake vertical pipe is assumed as 0.25; and the friction resistance coefficient of the pipeline is estimated as 0.015. The multi-point water intake pipeline belongs to a short composite pipe. In most cases, the friction resistance head loss may be ignored. In order to make the calculation example general, this example is considered. It may be seen from FIG. 3 that the branch pipes are equally spaced, the spacing is about 3 times the diameter D of the main pipe, and a length of each branch pipe is 7 times the pipe diameter. According to the existing pipeline confluence and research conclusion, this structure may ignore the neighborhood interference between the branch pipes. A head loss coefficient of a single confluence port may be calculated by using a previous test fitted curve or a Gardel formula. In FIG. 3, the size marking unit is centimeter.

(38) in view of the constant area A of the main pipe in this example, the area A.sub.i, the length l.sub.i and the spacing L.sub.i of each branch pipe are equal (i=1,2,3), so

(39) $R_{A i} = {(\frac{d}{D})}^{2} = 0 .25$ $K_{i, i n_{= 0.25}}$ $K_{i, f} = f \frac{l_{i}}{d_{i}} = 0.0 1 5 \frac{7}{1} = 0.1 05$ $K_{i - 1 ~ i, f} = f \frac{L_{i - 1 ~ i}}{D} = 0.0 1 5 \frac{6}{2} = 0.0 4 5$

(40) The specific calculating process is as follows: step 1: referring to the flow distribution in FIG. 4, the flow distribution proportions are preliminarily pre-estimated as Q.sub.1=0.28, Q.sub.2=0.32 and Q.sub.3=0.4. According to the Gardel formula of the T-shaped confluence head loss coefficient, the head loss coefficient of the confluence port is a function of a branch/main area ratio and a branch/main flow ratio:

(41) $K_{1, b .fwdarw. c} = {function}_{1} (\frac{A_{1}}{A}, \frac{Q_{1}^{}}{Q_{1}^{}}) = {function}_{1} (0.2 5, 1) 8.$ $K_{2, b .fwdarw. c} {function}_{1} (\frac{A_{2}}{A}, \frac{Q_{2}^{}}{(Q_{1}^{} + Q_{2}^{})}) = {function}_{1} (0.2 5, 0.5 3) 2.4$

(42) According to the Gardel formula of the T-shaped confluence head loss coefficient, the head loss coefficient of the confluence port is also a function of a branch/main area ratio and a branch/main flow ratio:

(43) $K_{2, t .fwdarw. c} {function}_{2} (\frac{A_{2}}{A}, \frac{Q_{2}^{}}{Q_{3}^{}}) = {function}_{2} (0.2 5, \frac{0.3 2}{0.2 8 + 0.3 2}) 0.69$ Step 2: a total loss coefficient K.sub.{circle around (2)} (corresponding to the flow velocity at a position {circle around (2)}) from a flow path p.sub.1 to the position {circle around (2)}: K.sub.{circle around (2)}|.sub.p1=local resistance coefficient

(44) $K_{1, i n} {(\frac{A}{A_{1}})}^{2} Q_{1}^{2}$ of inlet 1 +on-way friction resistance coefficient

(45) $K_{1, f} {(\frac{A}{A_{1}})}^{2} Q_{1}^{2}$ of vertical pipe 1 +branch flow confluence head loss coefficient K.sub.1,b.fwdarw.cQ.sub.1.sup.2 of confluence port 1 +on-way friction resistance coefficient K.sub.12, fQ.sub.1.sup.2 in the section 1-2 of the main pipe +overflow head loss coefficient K.sub.2,t.fwdarw.c(Q.sub.1+Q.sub.2).sup.2 of confluence port 2

(46) 0 $= (0.25 + 0.105) {(\frac{1}{0.25})}^{2} {0.28}^{2} + (8. + 0.045) {0.28}^{2} + 0.69 {(0.28 + 0.32)}^{2} = 1.1$ Step 3: a total head loss coefficient K.sub.{circle around (2)} from a flow path p.sub.2 to the position {circle around (2)}: K.sub.{circle around (2)}|.sub.p2=local resistance coefficient

(47) $K_{2, i n} {(\frac{A}{A_{2}})}^{2} Q_{2}^{2}$ of inlet 2 +on-way friction resistance coefficient

(48) $K_{2, f} {(\frac{A}{A_{2}})}^{2} Q_{2}^{2}$ of vertical pipe 2 +branch flow confluence head loss coefficient K.sub.2,b.fwdarw.c(Q.sub.1+Q.sub.2).sup.2 of confluence port 2

(49) $= (0.25 + 0.105) {(\frac{1}{0.25})}^{2} {0.32}^{2} + 2.4 {(0.28 + 0.32)}^{2} = 1.16$ Step 4: a difference value between values of |K.sub.{circle around (2)}|.sub.p1K.sub.{circle around (2)}|.sub.p2| obtained by two calculation methods is compared, |K.sub.{circle around (2)}|.sub.p1K.sub.{circle around (2)}|.sub.p2|=0.06, and there is a certain error.

(50) Flow distribution proportions are re-adjusted as Q.sub.1=0.285, Q.sub.2=0.315 and Q.sub.3=0.4. Then a flow path p.sub.1:

(51) $K_{2, t .fwdarw. c} {function}_{2} (0.25, \frac{0.315}{0.6}) = 0.68,$
so that K.sub.{circle around (2)}|.sub.p1==1.13 is obtained; and a flow path p.sub.2:

(52) $K_{2, b .fwdarw. c} {function}_{1} (0.25, \frac{0.315}{0.6}) = 2.3,$
so that K.sub.{circle around (2)}|.sub.p2=1.11 is obtained.

(53) At this time, it is considered that the loss coefficients calculated by the two flow paths are basically equal, so Q.sub.1=0.285, Q.sub.2=0.315 and

(54) $K \frac{1.1 3 + 1.1 1}{2} = 1.1 2$
are determined. Step 5: a head loss coefficient of a position {circle around (3)} is continuously calculated. Firstly, a head loss coefficient from a flow path p.sub.1+2 to the position {circle around (3)} is calculated. According to the above calculation, K.sub.{circle around (2)}=1.12; and

(55) $K_{3, t .fwdarw. c} {function}_{2} (\frac{A_{3}}{A}, \frac{Q_{2}^{}}{Q_{3}^{}}) = {function}_{2} (0.2 5, \frac{0.4}{1}) 0.52$
is deduced through the pre-estimated Q.sub.3=0.4 and Gardel formula of the T-shaped confluence head loss coefficient.

(56) Therefore, the total head loss coefficient from the flow path p.sub.1+2 to the position {circle around (3)} is:
K.sub.{circle around (3)}|.sub.p1+2=K.sub.{circle around (2)} +on-way friction resistance loss coefficient K.sub.23, f(Q.sub.1+Q.sub.2).sup.2 of the section 2-3 of the main pipe +overflow head loss coefficient K.sub.3,t.fwdarw.c(Q.sub.1+Q.sub.2+Q.sub.3).sup.2 of confluence port 3=1.12+0.0450.6.sup.2+0.521.0.sup.2=1.66 Step 6: a total head loss coefficient K.sub.{circle around (3)} from a flow path p.sub.3 to the position {circle around (3)} is calculated. According to the Gardel formula:

(57) $K_{3, b .fwdarw. c} = {function}_{1} (\frac{A_{3}}{A}, \frac{Q_{3}^{}}{(Q_{1}^{} + Q_{2}^{} + Q_{3}^{})}) = {function}_{1} (0.2 5, 0 .4) 1.4$

(58) Therefore, the total head loss coefficient from the flow path p.sub.3 to the position {circle around (3)} is: K.sub.{circle around (3)}|.sub.p3=local resistance coefficient

(59) $K_{3, i n} {(\frac{A}{A_{3}})}^{2} Q_{3}^{2}$ of inlet 3 +on-way friction resistance coefficient

(60) 0 $K_{3, f} {(\frac{A}{A_{3}})}^{2} Q_{3}^{2}$ of vertical pipe 3 +branch flow confluence head loss coefficient K.sub.3,b.fwdarw.t(Q.sub.1+Q.sub.2+Q.sub.3).sup.3 of confluence port 3

(61) $= (0.25 + 0.105) {(\frac{1}{0.25})}^{2} {0.4}^{2} + 1.4 {1.}^{2} = 1.86$ Step 7: K.sub.{circle around (3)} values obtained by the two calculation methods are compared. |K.sub.{circle around (3)}|.sub.p1K.sub.{circle around (3)}|.sub.p2|=0.2, that is, the head loss coefficient calculated by the flow path p.sub.3 is too large, so it is re-adjusted as a small flow Q.sub.3=0.38 (without adjusting Q.sub.1 and Q.sub.2), and the step 5 and the step 6 are repeated. The K.sub.{circle around (3)} of the two paths may be calculated to be 1.65 and 1.66 respectively. It may be considered that the two are equal, and accurate K.sub.{circle around (3)}=1.655 is obtained. Step 8: a next position custom character is continuously calculated, and the method is as same as the above. In this example, there are only three vertical pipes, that is, the position {circle around (3)} is the most downstream main pipe position. Then an overall head loss coefficient of a water intake system is:

(62) $K = \frac{K}{{(.Math. Q_{i})}^{2}} = \frac{1.655}{{(0.2 8 5 + 0.315 + 0.3 8)}^{2}} = 1.7 8$ Step 9: the above calculation has obtained the total head loss coefficient of the water intake system and the flow distribution proportion of each branch pipe. According to the design input data, the total head loss and the absolute flow of the branch pipe may be calculated:

(63) $U = \frac{Q}{A} = \frac{1 0}{2 . 0^{2} / 4} = 3.18 m / s$ $H = K \frac{U^{2}}{2 g} = 1.7 3 \frac{9.5 5^{2}}{2 9.8 1} = 0.894 m$ $Q_{1} = 3 0 \frac{0.2 8 5}{(0.2 8 5 + 0.3 1 5 + 0.3 8)} = 8.72 m^{3} / s$ $Q_{2} = 3 0 \frac{0.3 1 5}{(0.2 8 5 + 0.3 1 5 + 0.3 8)} = 9.64 m^{3} / s$ $Q_{3} = 3 0 \frac{0.3 8}{(0.2 8 5 + 0.3 1 5 + 0.3 8)} = 11.63 m^{3} / s$

(64) At this point, a flow velocity U of the most downstream main pipe, the overall head loss H and flows of the three branch pipes Q.sub.1, Q.sub.2 and Q.sub.3 are calculated, and the calculation process ends.

Simulation Experiment

(65) According to this embodiment, the flow distribution and the head loss coefficient of the method are calculated by an experimental model of a round branch pipe 6-point water intake shown in FIG. 5 so as to verify the correctness of the method according to this embodiment. In FIG. 5, the marking size unit is centimeter.

(66) (1) 6 round water-absorbing vertical pipes are arranged at the top of the main pipe at equal intervals, the upper end of each of the water-absorbing vertical pipes is an immersed mushroom head type water-absorbing head, environmental water flows into the water-absorbing head from the side, enters the vertical pipe and then sequentially flows into the main pipe, and the specific size of the model is shown in FIG. 5.

(67) (2) A head loss coefficient of the water-absorbing head inlet: since the six water-absorbing heads have the consistent shape and size, the local head loss coefficients are the same. According to the test measurement, the local resistance coefficient .sub.i of the used water-absorbing head 0.574 (lateral water inlet 0.5+round reducing pipe 0.074).
K.sub.i,in=.sub.i.sub.0.574

(68) (3) An on-way friction resistance loss coefficient of the vertical pipe: the on-way friction resistance loss coefficient is calculated by using a Colobrook-White formula and according to .sub.i and Re.sub.i, an absolute roughness of a model pipe is 2.0 m, and the relative roughness height of the branch pipe

(69) $_{i} = \frac{}{d} = 2 1 0^{- 5} .$
A transition section with a 1d length at the downstream of a gradually changed section of the water-absorbing head is ignored when the on-way loss is calculated, that is, an effective friction resistance distance is li=ld=0.3 m.

(70) $K_{i, f} =_{i} (_{i}, {Re}_{i}) \frac{l_{i}}{d} =_{i} (2 1 0^{- 5}, \frac{4 Q_{i}}{d}) \frac{l_{i}}{d}$

(71) A spacing L between the branch pipes is equal to 3D, so the on-way loss h.sub.i-1i,f0 of the main pipe is ignored.

(72) (4) A confluence port loss coefficient: the spacing among six continuous confluence ports are equal and are six times the diameter of the branch pipe. According to the previous research, most local water losses caused by confluence occur within this spacing, that is, the mutual influence of the confluence ports is small.

(73) A confluence energy loss coefficient K.sub..fwdarw.cR.sub.Q adopts a formula of an optimized quadratic function type:
K.sub.t.fwdarw.c(R.sub.Q)|.sub.R.sub.A.sub.=0.25,radii=0=0.233R.sub.Q.sup.2+1.598R.sub.Q
K.sub.b.fwdarw.c(R.sub.Q)|.sub.R.sub.A.sub.=0.25,radii=0=1.06R.sub.Q.sup.2+5.421R.sub.Q1

(74) (5) In this test, the upstream of the confluence port 1 is closed (that is, Q.sub.0=0), and there is a section of convex stagnant water area (part of construction space is generally remained at the end of a pipe shield in engineering practice). The local loss coefficient of the confluence port is similar to the confluence of R.sub.Q=1, but due to the influence of the reverse reflux in the stagnant water area, the local head loss coefficient needs to be calculated specially.

(75) (6)

(76) $K_{1, i n} = {_{1 0.574,} (\frac{A}{A_{1}})}^{2} = 16 and K_{1, f} =_{1} (_{1}, {Re}_{1}) \frac{l}{d} = 0.0178 \frac{l}{d} = 0.1065$
are substituted into an equation to obtain K.sub.1,b.fwdarw.c=5.542.

(77) (7) A relationship among the number of the branch pipes, the flow distribution and the total loss coefficient is calculated:

(78) According to an iterative calculation method provided by this embodiment, the flow distribution and the overall head loss coefficient of the multi-point water intake of different numbers of branch pipes may be solved. The iteration steps are completed by computer programming optimization iteration.

(79) When Q=25.11 L/s, the flow and the head loss coefficient of the n=1-6 branch pipes are calculated according to the measured K.sub.1=16.43. It may be seen from Table 1 and Table 2 that by the method according to the embodiment, the head loss distribution and flow distribution of different numbers of pipelines may be calculated accurately.

(80) TABLE-US-00001 TABLE 1 Comparison between the flow distribution of each branch pipe and the head loss coefficient of the downstream of the confluence port, and the test result Branch pipe Branch pipe Branch pipe Branch pipe Branch pipe Branch pipe Type Branch pipe 1# 2# 3# 4# 5# 6# Actually K 16.01/16.43 3.51 2.30 1.56 1.21 1.01 measured Flow 10.6% 11.9% 14.0% 16.4% 20.7% 26.4% in the distribution experiment proportion Calculated K.sub.n 16.43 4.64 2.43 1.60 1.16 0.87 by this Flow 11.4% 12.7% 14.4% 17.1% 20.3% 24.1% method distribution proportion

Embodiment II

(81) This embodiment is an improvement of the embodiment I and a refinement of the main pipe in the embodiment I. In this embodiment, the on-way diameters of the main pipes are equal, but when the on-way diameters of the main pipes change, the loss coefficient K may be corrected:

(82) $K_{correct} = {K_{origin} (\frac{A_{d}}{A_{u}})}^{2}$

(83) wherein A.sub.d is area of a downstream main pipe, A.sub.u is area of an upstream main pipe, K.sub.origin is an uncorrected head loss, and K.sub.correct is a corrected head loss.

(84) Finally, it should be noted that the above is only used to describe the technical solution of the present invention, but not to limit it. Although the present invention is described in detail with reference to the form of the 3-point water intake, the mushroom head type water-absorbing head and the Gardel formula of the confluence flow head loss coefficient, those of ordinary skill in the art should understand that the technical solutions (such as different forms of water-absorbing heads, different numbers of water-absorbing vertical pipes and different confluence flow head loss coefficient formulas) of the present invention may be corrected or equivalently replaced without departing from the spirit and scope of the technical solutions of the present invention.

Method for calculating head loss coefficient and branch pipe flow distribution of multi-point water intake

Assignee

Inventors

Cpc classification

Classification Explorer

G06F30/18

PHYSICS

Classification Explorer

G06Q10/0631

PHYSICS

Classification Explorer

G06F2113/14

PHYSICS

Classification Explorer

G06F30/28

PHYSICS

International classification

Classification Explorer

G06F30/28

PHYSICS

Classification Explorer

G06F113/14

PHYSICS

Classification Explorer

G06F30/18

PHYSICS

Classification Explorer

G06Q10/0631

PHYSICS

Abstract

Claims

Description