Mathematical Problems in Engineering
Volume 2018 (2018), Article ID 4852361, 16 pages
https://doi.org/10.1155/2018/4852361

Research Article

Determination and Application of Comprehensive Specific Frictional Resistance in Heating Engineering

Yanan Tian,¹ Haosen Wang,² Yongjiang Shi,¹ Jiayu Rong,¹ Zhenhua Ma,² and Cong Xie²

¹College of Energy and Environmental Engineering, Hebei Institute of Architecture and Civil Engineering, Zhangjiakou, Hebei 075000, China
²College of Mathematics and Physics, Hebei Institute of Architecture and Civil Engineering, Zhangjiakou, Hebei 075000, China

Correspondence should be addressed to Yongjiang Shi; yy1114916122@163.com

Received 9 July 2017; Revised 12 November 2017; Accepted 22 November 2017; Published 4 January 2018

Academic Editor: Andrea L. Facci

Copyright © 2018 Yanan Tian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In this study, we analyze the deficiencies of specific frictional resistance in heating engineering. Based on economic specific frictional resistance, we put forward the concept of comprehensive specific frictional resistance, which considers the multiple factors of technology, economy, regulation modes, pipe segment differences, and medium pressure. Then, we establish a mathematical model of a heating network across its lifespan in order to develop a method for determining the comprehensive specific frictional resistance. Relevant conclusions can be drawn from the results. As an application, we have planned the heating engineering for Yangyuan County in China, which demonstrates the feasibility and superiority of the method.

1. Introduction

With the development of modern society, improvements in heating engineering have attracted attention. In the field of heating engineering, the specific frictional resistance refers to the resistance loss of pipe per meter [1, 2], which is very important in China. In this paper, we discuss specific frictional resistance during the designing or planning stage, when specific frictional resistance acts as an intermediary to select a suitable pipe diameter.

At present, studies mainly involve two aspects of specific frictional resistance: one is the weighted average specific frictional resistance; the other is economic specific frictional resistance. For a heating network with regulated change in flow rate at different stages, the specific frictional resistance undergoes a step change due to the change in flow rate. In order to measure and reflect this characteristic of the specific frictional resistance, weighted average specific frictional resistance is recommended for selecting pipe diameters. Hong [3] discussed the differences between weighted average frictional resistance and traditional specific frictional resistance and analyzed the power-saving effect by taking a heating network with flow rate regulation as an example. The economic specific frictional resistance is an ideal value relative to the specific frictional resistance from the economic perspective and is widely used in China at present. Ge et al. [4] established a mathematical model, which took the minimum cost of the heating network as the objective function, and discussed the relationship between the economic specific frictional resistance and the economic heating radius. Liu et al. [5] offered the concept of the new economic specific frictional resistance by establishing the transportation distance limit model and the transport energy consumption ratio model. Wang et al. [6] took the indoor heating horizontal single tube and vertical single tube system as examples to establish a mathematical model with minimum cost as the objective function, by which they obtained the economic specific frictional resistance under different conditions. Li et al. [7] studied the hydraulic calculation methods and the existing problems on the condensate pipe net, and they used the differential method to calculate the economic specific frictional resistance of the condensate pipe net.

Through analysis of the various research [3–11] on specific frictional resistance, we identified problems related to the value of specific frictional resistance, as follows:

(i) In China, the value of specific frictional resistance used in engineering still follows the traditional economic specific frictional resistance value of 30~70 Pa/m, which was introduced by the former Soviet Union in the last century. With the rise in prices and changes in materials, this value is not applicable to the current state of practice in heating engineering.

(ii) The value of specific frictional resistance should be related to the regulation modes, which can bring about different changes in flow rate and medium temperature. However, the traditional model of economic specific frictional resistance treats the different regulation modes in the same way [3–11].

(iii) With the development of the central heating network, pipe networks have become more massive and complex, which has led to increasing differences among pipe segments. Thus, it is not suitable to use the same specific frictional resistance standard of 30~70 Pa/m to determine the pipe diameter in a network system. Instead, each pipe segment should be differentiated. That is to say, each pipe segment should have a specific frictional resistance of its own [3–11].

(iv) The distribution of pressure in the pipes is now more complicated because of the mass and complexity of the pipe networks. It has a direct relationship to the value of specific frictional resistance whether the pressure meets the requirements. However, the value of traditional economic specific frictional resistance of 30~70 Pa/m does not take into account the constraints of the medium pressure [3–11].

In order to solve the above problems, we propose the concept of comprehensive specific frictional resistance. It refers to a measure of resistance loss of pipe per meter that incorporates a number of related factors. By means of it, heating engineering applications can achieve better economic effects during their lifespan [12–14]. The related factors considered in this paper include technology, economy [15, 16], operation regulation modes [17, 18], pipe segment differences, and the medium pressure [19–22].

2. The Mathematical Model

In this section, we give the method for obtaining the value of comprehensive specific friction resistance as follows.

2.1. Factors Related to Comprehensive Specific Frictional Resistance Theory

Through the analysis of technology, economy, operation regulation modes, pipe segment differences, and the medium pressure in a heating system, we have drawn several conclusions about the comprehensive specific friction resistance model theory:

(i) Pipe Segment Differences. With pipe networks becoming massive and complex, the difference among pipe segments becomes more apparent. For example, it is common for the pipe diameter close to the heat source to be DN1400, while the pipe diameter can be reduced to DN200 at the user end point. This difference should not be ignored. In order to account for the differences among pipe segments, a single pipe segment is taken as an object in establishing the mathematical model in this paper. In this way, we get different values corresponding to different pipe segments.

(ii) Operation Regulation Modes. We take into account the influence of operation regulation modes on comprehensive specific frictional resistance due to the fact that different operation regulation modes incur different electricity consumption and different heat loss.

(iii) Medium Pressure. The value of comprehensive specific frictional resistance takes into account the medium pressure constraint, which can prevent heating accidents on the operational stage, such as pipe cracking or medium vaporization. In the following study, the medium pressure constraint is mainly embodied in the constraint conditions of the mathematical model.

(iv) Technology and Economy. The Present Value Cost is chosen as the objective function [23, 24] by analyzing the advantages and disadvantages of each economic effect evaluation index.

(v) Technology and Economy. The value of specific frictional resistance is sensitive to the factors of primary network construction fee (including civil engineering cost and pipe installation fee), electricity fee, and fuel cost across the lifespan of heating engineering application. These factors must be embodied in the mathematical model. Heat source construction fee, heat exchange station construction fee, thermal user fuel cost, water fee, maintenance cost, and labor wages are irrelevant factors that can be ignored in the mathematical model.

2.2. Establishment of the Objective Function

In this section, we establish the objective function taking the single segment as the object. The objective function is the Present Value Cost of pipe segment across the lifespan of years, and the cost includes civil engineering cost, pipe installation fee, electricity fee, and fuel cost.

For the different regulations, the civil engineering cost and the pipe installation fee are the same, but the electricity fee and fuel cost are different, because the heat load varies continuously during heating season. The following is the establishment process of the objective function.

2.2.1. The Civil Engineering Cost

Taking the directly buried laying method as an example, a set of data about pipe diameter and the corresponding civil engineering cost from Chinese markets is shown in Figure 1.

Figure 1: Fitting curve of heating pipe diameter and civil engineering cost.

The fitted curve is , whose -squared value is as high as 99.95%. Hence, the relationship between the civil engineering cost and the pipe diameter can be expressed by a quadratic polynomial, as follows:

where is the civil engineering cost of pipe with a diameter (Ұ/m), is the pipe diameter (m), and , , and are the fitting polynomial coefficients of civil engineering cost, which are determined by the specific situation.

As we all know, the Darcy-Weisbach Formula is

where is the head loss (m), is the coefficient of frictional resistance [25–27], is the length of pipe (m), is the velocity of medium in the pipe (m/s), and is the gravitational acceleration (m/s²).

The specific frictional resistance refers to the resistance loss of pipe per meter, so we get formula (3) from formula (2):

where ρ is the density of medium (kg/m³) and is the specific frictional resistance (Pa/m).

Therefore, the specific frictional resistance, pipe diameter, and flow rate have the following relationship:

where is the flow rate (kg/s).

For a heating network, the determination of the comprehensive specific frictional resistance needs to account for the change in flow rate during operation. However, the selection of the diameter through the flow rate and the specific frictional resistance still relates to the most unfavorable situation. Therefore, the flow rate in formula (4) is the design flow rate . That is,

where is the pipe diameter of pipe segment (m), is the design flow rate of pipe segment (kg/s), and is the specific frictional resistance of pipe segment (Pa/m).

Applying the basic formula of specific heat capacity ( is the specific heat capacity, which for water is 4.2 × 10³ J/(kg·°C), m is the quality (kg), is the temperature difference (°C), and is the heat quantity (J)) to heating engineering, the following formula can be obtained:

where is the thermal load (W), is the temperature difference (°C), and is the flow rate of medium (kg/s).

Substituting formula (6) into formula (5), we get

where , is the design temperature difference between supply and return water of pipe segment (°C), and is the design thermal load of pipe segment (w).

Substituting formula (7) into formula (1), the civil engineering cost of pipe segment can be given by

where is the civil engineering cost of pipe segment (Ұ) and is the length of pipe segment (m).

2.2.2. The Pipe Installation Fee

The data from Chinese markets about pipe diameter and the corresponding pipe installation fee is shown in Figure 2.

Figure 2: Fitting curve of heating pipe diameter and pipe installation fee.

The fitted curve is , whose -squared value is as high as 99.73%. Hence, the relationship between the pipe installation fee and the pipe diameter can also be expressed by a quadratic polynomial, as follows:

where is the installation fee of pipe with a diameter (Ұ/m) and , , and are the fitting polynomial coefficients of the installation fee, which are determined by the specific situation.

Substituting formula (7) into formula (9), the installation fee of the pipe segment can be given by

where is the installation fee of the pipe segment (Ұ).

2.2.3. The Electricity Fee

(i) For a Heating System with Quality Regulation. A heating system with quality regulation is regulated simply by changing the water-supply temperature, so the flow rate of pipe segment is always the design flow rate . Because the specific frictional resistance R refers to the resistance loss of pipe per meter, the annual electricity consumption of pipe segment is

where is the annual electricity consumption of pipe segment (J) and is the run time of pipe segment per year (s).

So, the annual electricity fee for the pipe segment can be given by

where is the annual electricity fee for pipe segment (Ұ), is the electrovalence (Ұ/J), and is the pump efficiency.

Substituting formula (6) into formula (12), we get

(ii) For a Heating System with Flow Rate Regulation. A heating system with quality regulation is regulated by changing the flow rate. The large variation on the thermal load during the entire heating period results in a wide range of flow rate changes, so calculating the electricity fee according to the design flow rate does not adequately describe the actual situation. In order to estimate the electricity consumption accurately, we need to know the local thermal load changes and the corresponding hours to draw a thermal load duration diagram [28]. Under common conditions, the thermal load duration of the pipe segment in heating engineering can be described as Figure 3.

Figure 3: Thermal load duration diagram of the pipe segment .

According to the definition of specific frictional resistance , the annual electricity consumption of pipe segment is

where is the actual flow rate of pipe segment (kg/s), is the flow rate (kg/s)-duration (s) function, and is the duration corresponding to the minimum thermal load during a heating season (s).

The annual electricity fee for the pipe segment is

Substituting formula (6) into formula (15), the annual electricity fee for the pipe segment can be given by

where , is the actual thermal load of pipe segment (w) and is the thermal load (w) – duration (s) function, as shown in Figure 3.

2.2.4. The Fuel Cost

The fuel cost includes two parts: one is heat loss fuel cost; the other is thermal user fuel cost. Due to the irrelevance of the specific frictional resistance to thermal user fuel cost, there is no need to embody it in the mathematical model. However, massive, complicated pipeline networks increase heat loss along the route, so heat loss fuel cost should be included in the total cost rather than ignored. The heat loss fuel cost is different for different regulation modes, because it is associated with the medium temperature.

(i) For a Heating System with Flow Rate Regulation. The heat transfer coefficient refers to the amount of heat transfer per unit time per unit surface area for a 1°C temperature difference, so for a heating network adopting the flow rate regulation mode, the consumption of heat loss fuel is

where is the consumption of heat loss fuel (J), is the heat transfer coefficient of the pipe surface (w/m²·°C), is the ambient temperature of pipe segment (°C), and is the design temperature of medium in pipe segment (°C). For the water-supply pipe, , and for the water-return pipe, , , where is the design water-supply temperature of pipe segment (°C), is the design water-return temperature of pipe segment (°C).

So the annual heat loss fuel cost is

where is the annual heat loss fuel cost of the pipe segment (Ұ) and is the fuel price (Ұ/J).

Substituting formula (7) into formula (18), the annual heat loss fuel cost of the pipe segment can be given by

(ii) For a Heating System with Quality Regulation. For heat-supply networks adopting the quality regulation mode, calculating the heat loss fuel cost by the design temperature is unreasonable due to the variance of the medium temperature in the pipe. As we know, compared with the design water-supply temperature and design water-return temperature, a reduction in actual water-supply temperature of Δ1 is equal to an increase of Δ2 of the actual water-return temperature in the quality regulation networks [8]. This is described in Figure 4.

Figure 4: Duration of actual water-supply or water-return temperature.

So, the following formula can be given:

where is the actual temperature of medium in pipe segment (°C) and is the actual temperature difference between supply and return water of pipe segment (°C), as shown in Figure 4.

According to the definition of heat transfer coefficient , the consumption of heat loss fuel is

Substituting formula (20) into formula (21), we get the consumption of heat loss fuel:

So, the annual heat loss fuel cost of the pipe segment is

Substituting formula (7) into formula (23), the annual heat loss fuel cost of the pipe segment can be given by

For a heating system with quality regulation, the actual flow rate is always the design flow rate , so we know

Substituting formula (25) into formula (24), we get

2.2.5. The Objective Function

For pipe segment , the objective function is the Present Value Cost in its lifespan of years. Supposing the discount rate is , it can be described as follows.

(i) For a Heating System with Quality Regulation

(ii) For a Heating System with Flow Rate Regulation

where PC is the Present Value Cost (Ұ), is the discount rate, and is the life period (years), , , , .

2.3. Establishment of the Constraint Conditions

2.3.1. The Constraint Conditions of Medium Pressure

The pressure in the pipe must meet certain requirements in order to ensure normal operation of the system and prevent pipe burst or medium vaporization. As the research object is the primary network system, the lower pressure limit must guarantee that high-temperature water does not vaporize in the pipeline, while the upper limit is the pressure that the pipe and pipe fittings can bear. Additionally, the water-return temperature is generally less than 100°C, and for any water-supply pipe segment , there must be a corresponding water-return pipe segment whose pressure is lower than the pressure in the pipe segment . Therefore, only the water-supply pipes need to obey the pressure constraints. In addition, in order to ensure that the pressure in the pipe segment meets the requirements, the starting and ending points of the pipe must also meet the requirements.

A branched network diagram without distributed pressure pumps is shown in Figure 5. Pipe segments constitute the trunk line. Under normal circumstances, the constant pressure point is located at the entrance of the circulating pumps, that is, the end of pipe segment 1′. Assuming that the pressure at the constant pressure point is , the pressure at any section of any pipe segment can be calculated.

Figure 5: A branched network diagram.

According to the Bernoulli equation (Bernoulli’s principle), the pressure at the start of pipe segment can be given by

So the pressure constraint condition of pipe segment can be given by

where is the pressure at the start of pipe segment (Pa), is the flow velocity of the medium in pipe segment (m/s), is the topographic height at the start of pipe segment (m), is the pressure at the end of pipe segment (the constant pressure point) (Pa), is the flow velocity of the medium in pipe segment (m/s), is the topographic height at the end of pipe segment (m), is the resistance loss at heat exchange station (Pa), is the ratio of the local resistance to the resistance along the path, is the saturation steam pressure corresponding to the design water-supply temperature (Pa), is the maximum pressure that the pipe and pipe fittings can bear (Pa), and is the volumetric weight (N/m³), .

In order to simplify formula (30), the item can be moved as follows:

There is little difference among different pipelines in terms of flow velocity, for which the range is 0~4 m/s [29, 30]; that is, the value of is −8000~8000 Pa. Compared with and , its order of magnitude is low. Therefore, formula (31) can be simplified to the following form:

In this way, not only has the formula been simplified, but also a certain margin has been set for the upper and lower limits.

From the Antoine equation [31], we know that  Pa, so the pressure constraint condition at the start of pipe segment can be given by

where .

As described above, the pressure at the end of pipe segment is equal to the pressure at the start of pipe segment because the velocity difference has been set as a margin for the upper and lower limits, which means the pressure constraint at the end of pipe segment has been validated in the previous step by the pressure constraint of pipe segment . Therefore, the pressure constraint at the end of pipe segment can be ignored.

2.3.2. The Constraint Condition of Design Temperature Difference between Supply and Return Water

The design temperature difference between supply and return water should be limited to a certain range. Too small a temperature difference can result in reduced and inadequate heat exchange, while too large a temperature difference may need a very high water-supply temperature, which may be difficult to achieve or may incur a very high cost in practice. Therefore, the design temperature difference between supply and return water should be limited by

2.3.3. The Constraint Condition of the Flow Velocity

The flow velocity in the pipes should be restricted to a certain range, as follows:

Substituting formulas (6) and (7) into formula (35), the constraint conditions of the flow velocity can be given by

2.4. Establishment of a Mathematical Model

Based on the above analysis, a mathematical model of Present Value Cost to determine the comprehensive specific frictional resistance of pipe segment can be given as follows.

(i) For a Heating System with Quality Regulation

(ii) For a Heating System with Flow Rate Regulation

For the water-supply pipes, the constraint conditions are the inequalities (39), (40), and (41), while, for the water-return pipes, the constraint conditions are the inequalities (39) and (40).

2.5. The Solvability Analysis of the Mathematical Model

The objective function PC(, ) we established is a continuous real valued function map, and the value ranges of the decision variable and are compact sets. Therefore, we can know by the Weierstrass Theorem [32] that there must be a solution vector (, ) that makes the objective function a minimum.

2.6. The Analysis of Algorithm

The mathematical model established above belongs to the two-dimensional extreme value problem under inequality constraints. The Complex method is selected as the optimization algorithm and is used to solve the models, whose function is to solve the -dimensional extremum problem under the equality and inequality constraints.

2.7. The Value of Comprehensive Specific Frictional Resistance

According to the actual engineering, the basic parameters in the objective function and the constraint conditions can be determined, while the total heat load of the heat consumers that the pipe segment provides heat for during the heating period can be calculated. Substitute all of them into the mathematical models (37)~(41) and solve the mathematical model. We can obtain the comprehensive specific frictional resistance of the pipe segment and the design temperature difference between supply and return water .

3. Application

3.1. Case

We have planned the heating engineering for Yangyuan County in China. The heating system is based on flow rate regulation as the main regulation mode, considering its energy-saving advantage. In the planning process of the network, the two schemes, comprehensive specific frictional resistance and economic specific frictional resistance, are used to optimize the pipe diameters.

According to market research, the calculation parameters have been selected as follows.

Coefficient of frictional resistance , discount rate , electrovalence  Ұ/J, fuel price  Ұ/J, the lifespan years, run time per year  s (equivalent to 168 days), pump efficiency , ambient temperature of pipe segments °C, heat transfer coefficient of pipe surface  w/(m²·°C), design water-return temperature °C, fitting polynomial coefficients of civil engineering cost , , and , fitting polynomial coefficients of pipe installation fee , , and , the ratio of local resistance to resistance along the path , the resistance loss at heat exchange station  MPa,  MPa, upper bound of flow velocity  m/s, lower bound of flow velocity  m/s, upper bound of design temperature difference between supply and return water °C, lower bound of design temperature difference between supply and return water °C, upper bound of pressure  MPa,  kg/m³, and  N/m³, and the duration of outdoor temperature in Zhangjiakou are shown in Table 1.

Table 1: Duration of outdoor temperature in Zhangjiakou.

Number
Outdoor temperature (°C)
Duration of outdoor temperature (hours)

>5
4728

5
4032

4
3850

3
3667

2
3483

1
3299

0
3114

2927

2740

2552

2363

2173

1981

1788

1593

1397

1197

996

790

579

361

120

According to the site investigations, and to follow the trends of the existing pipeline network, the whole heating area was divided into three independent heating areas as shown in Figure 6: the planning floor area of heating area A is 9.61 × 10⁶ m², the planning floor area of heating area B is 6.57 × 10⁶ m², and the planning floor area of heating area C is 4.75 × 10⁶ m².

Figure 6: The map of heating areas A, B, and C.

In this study, we take heating area C as an example. In Figure 7, only water-supply pipes are marked, while the position and trend of water-return pipes are in agreement with the corresponding water-supply pipes. In the following description, the numbers with “′” represent the water-return pipes, while the numbers without “′” represent the water-supply pipes.

Figure 7: Heating network in heating area C.

The entire calculation process is as follows:

① Choose the most distant loop from the heat source (the trunk line), that is, ;

② Calculate the and of the water-return pipes in the most distant loop with the models (38), (39), and (40);

③ Follow the order to calculate the and of the water-supply pipes in the most distant loop with the models (38), (39), (40), and (41);

④ For the other branch loops, calculate the and of the water-return pipes that have not been calculated with the models (38), (39), and (40), and then calculate the and of the water-supply pipes that have not been calculated, following the order from far to near with the models (38), (39), (40), and (41).

The calculation results are shown in Table 2. As expected, the range of values for traditional economic specific frictional resistance is 30~70 Pa/m.

Table 2: Diameter optimization result.

Number
Thermal load
(MW)
Design water-supply temperature
(°C)
Design temperature difference
(°C)
Design flow rate
(kg/s)
Economic specific frictional resistance
Comprehensive specific frictional resistance

Pipe diameter
(m)
Specific frictional resistance
(Pa/m)
Flow rate
(m/s)
Specific frictional resistance
(Pa/m)
Pipe diameter
(m)
Corrected pipe diameter
(m)
Flow rate
(m/s)

100
184.00
130
60
732.59
0.7
51
1.90
65
0.666
0.7
1.9

101
168.61
130
60
671.32
0.7
43
1.74
67
0.39
0.6
2.4

102
160.05
130
60
637.24
0.7
38
1.66
68
0.624
0.6
2.3

103
108.22
130
60
430.88
0.6
38
1.52
77
0.521
0.5
2.2

104
45.50
130
60
181.16
0.4
51
1.44
105
0.346
0.35
1.9

105
30.61
130
60
121.87
0.35
45
1.27
122
0.287
0.3
1.7

106
20.74
130
60
82.58
0.3
45
1.17
142
0.238
0.25
1.7

107
10.62
120
50
50.74
0.25
42
1.03
197
0.183
0.20
1.6

108
15.39
130
60
61.28
0.25
61
1.25
161
0.206
0.2
2.0

109
9.18
130
60
36.55
0.2
66
1.16
199
0.161
0.15
2.1

110
6.21
130
60
24.73
0.2
30
0.79
235
0.133
0.15
1.4

111
8.56
130
60
34.08
0.2
58
1.08
206
0.155
0.15
1.9

112
35.79
130
60
142.50
0.4
31
1.13
115
0.309
0.3
2.0

113
8.66
130
60
34.48
0.2
59
1.10
206
0.156
0.15
2.0

114
9.05
130
60
36.03
0.2
64
1.15
201
0.159
0.15
2.0

115
18.08
130
60
71.99
0.3
34
1.02
151
0.223
0.2
2.3

116
9.66
130
60
38.46
0.2
73
1.22
197
0.164
0.15
2.2

117
8.42
130
60
33.52
0.2
56
1.07
206
0.154
0.15
1.9

118
16.04
130
60
63.86
0.25
66
1.30
158
0.210
0.2
2.0

119
7.82
130
60
31.14
0.2
48
0.99
215
0.148
0.15
1.8

120
8.22
130
60
32.73
0.2
53
1.04
209
0.152
0.15
1.9

121
62.72
130
60
249.72
0.5
32
1.27
94
0.402
0.4
2.0

122
17.95
130
60
71.47
0.3
33
1.01
152
0.222
0.2
2.3

123
10.53
130
60
41.93
0.25
29
0.85
188
0.172
0.15
2.4

124
34.24
130
60
136.33
0.35
56
1.42
117
0.302
0.3
1.9

125
16.86
130
60
67.13
0.25
73
1.37
155
0.215
0.2
2.1

126
8.33
130
60
33.17
0.2
55
1.06
209
0.153
0.15
1.9

127
9.05
130
60
36.03
0.2
64
1.15
201
0.159
0.15
2.0

128
7.68
130
60
30.58
0.2
46
0.97
215
0.147
0.15
1.7

129
7.21
130
60
28.71
0.2
41
0.91
221
0.143
0.15
1.6

130
9.87
130
60
39.30
0.25
25
0.80
194
0.166
0.15
2.2

131
10.12
130
60
40.29
0.25
26
0.82
192
0.168
0.15
2.3

100′
184.00
130
60
732.59
0.7
51
1.90
64
0.645
0.7
2.37

101′
168.61
130
60
671.32
0.7
43
1.74
65
0.630
0.6
2.25

102′
160.05
130
60
637.24
0.7
38
1.66
74
0.525
0.6
2.19

103′
108.22
130
60
430.88
0.6
38
1.52
100
0.350
0.5
1.88

104′
45.50
130
60
181.16
0.4
51
1.44
115
0.290
0.35
1.72

105′
30.61
130
60
121.87
0.35
45
1.27
134
0.241
0.3
1.68

106′
20.74
130
60
82.58
0.3
45
1.17
176
0.174
0.25
2.39

107′
10.62
130
60
42.28
0.25
29
0.86
151
0.209
0.15
1.95

108′
15.39
130
60
61.28
0.25
61
1.25
185
0.163
0.2
2.07

109′
9.18
130
60
36.55
0.2
66
1.16
219
0.135
0.15
1.40

110′
6.21
130
60
24.73
0.2
30
0.79
192
0.157
0.15
1.93

111′
8.56
130
60
34.08
0.2
58
1.08
109
0.312
0.15
2.02

112′
35.79
130
60
142.50
0.4
31
1.13
192
0.158
0.3
1.95

113′
8.66
130
60
34.48
0.2
59
1.10
188
0.161
0.15
2.04

114′
9.05
130
60
36.03
0.2
64
1.15
141
0.226
0.15
1.47

115′
18.08
130
60
71.99
0.3
34
1.02
183
0.167
0.25
2.18

116′
9.66
130
60
38.46
0.2
73
1.22
192
0.156
0.15
1.90

117′
8.42
130
60
33.52
0.2
56
1.07
148
0.213
0.15
2.03

118′
16.04
130
60
63.86
0.25
66
1.30
200
0.150
0.2
1.76

119′
7.82
130
60
31.14
0.2
48
0.99
195
0.154
0.15
1.85

120′
8.22
130
60
32.73
0.2
53
1.04
89
0.407
0.15
1.99

121′
62.72
130
60
249.72
0.5
32
1.27
142
0.225
0.4
1.46

122′
17.95
130
60
71.47
0.3
33
1.01
176
0.174
0.25
2.37

123′
10.53
130
60
41.93
0.25
29
0.85
111
0.306
0.15
1.93

124′
34.24
130
60
136.33
0.35
56
1.42
145
0.218
0.3
2.14

125′
16.86
130
60
67.13
0.25
73
1.37
195
0.155
0.2
1.88

126′
8.33
130
60
33.17
0.2
55
1.06
188
0.161
0.15
2.04

127′
9.05
130
60
36.03
0.2
64
1.15
200
0.149
0.15
1.73

128′
7.68
130
60
30.58
0.2
46
0.97
206
0.145
0.15
1.62

129′
7.21
130
60
28.71
0.2
41
0.91
181
0.168
0.15
2.22

130′
9.87
130
60
39.30
0.25
25
0.80
180
0.170
0.15
2.28

131′
10.12
130
60
40.29
0.25
51
1.90
64
0.645
0.15
2.37

The optimization results in Table 2 show that the design water-supply temperature for pipe segment 107 is 120°C, which is different from all the others. The reasons for this are the long distance from the heat source to the pipe and the high elevation of the pipe, which lead to a lot of pressure reduction in the pipe segment. Therefore, the only way to ensure the pipe pressure is higher than the vaporization pressure is to decrease the water-supply temperature. If the pressure constraint formula (41) is removed, the result for pipe segment 107 will be , °C, at which temperature the medium will vaporize. Because there must be one water-supply temperature for the entire engineering, so this optimization result is not desirable. There are two solutions to this problem: ① the entire heating system is recalculated according to a lower water-supply temperature in all the results until the optimization results of all the segments are same; ② the heat consumers after the pipe segments whose optimization results of water-supply temperature are different from the others take the indirect connection form. In this case, taking into account the fact that only one segment pipe is different from the others, it is obvious that the second solution is more reasonable. That is, pipe segment 107 and the heat user 57 should connect with the trunk line indirectly.

From the perspective of the entire engineering application, based on the calculation of the various costs and incomes, including heat source construction fee, heat exchange station construction fee, primary network construction fee, electricity fee, heat loss fuel cost, thermal user fuel cost, water fee, maintenance cost, labor wage, heating charges, and residual value of fixed assets, the economic effect evaluation indexes of the two schemes are shown in Table 3.

Table 3: Comparison table of economic effect evaluation indexes.

Schemes
Index
Payback Period
(years)
Net Present Value
(10000)
Rate of Return on Investment
(%)
Internal Rate of Return
(%)
Present Value Cost
(10000)
Initial Investment
(10000)

Scheme one
(Comprehensive specific frictional resistance)
11.39
40872.73
39.13
11.7
32116.94
104463.1

Scheme two
(Economic specific frictional resistance)
12.04
37460.97
34.15
11.1
34477.33
109679.8

By comparison, we can see that under the premise of ensuring the same heating quality, the economic effect evaluation indexes of scheme one are better. Compared to scheme two, the Payback Period is lowered, the Initial Investment is reduced by 5%, the Present Value Cost is decreased by Ұ23.6039 million, the Net Present Value has an increase of Ұ34.1176 million, the Rate of Return on Investment is increased by 4.98%, and the Internal Rate of Return is increased by 0.6%, which demonstrate the feasibility and superiority of the comprehensive specific frictional resistance.

3.2. Discussion

Through the analysis of the case and the mathematical models, we show there are many factors that influence the value of comprehensive specific frictional resistance, including the coefficient of frictional resistance , heat transfer coefficient of pipe surface , electrovalence , fuel price , discount rate , lifespan , pump efficiency , and primary network construction fee. The various values of the above parameters in different engineering applications will lead to the different values of comprehensive specific frictional resistance. In addition, the comprehensive specific frictional resistance has different sensitivities for different parameters.

Using the calculation parameters in Section 3.1, and setting the design water-supply temperature °C, we get the sensitivities as shown in Figure 8.

Figure 8: Effects of various factors on comprehensive specific frictional resistance.

Similarly, using the calculation parameters in Section 3.1, and setting the design water-supply temperature °C, the comprehensive specific frictional resistance corresponding to each flow rate can be solved, by which we can draw the curves of flow rate-Comprehensive specific frictional resistance as in Figure 9.

Figure 9: Supply and return water flow rate-comprehensive specific frictional resistance.

Some conclusions can be drawn from the Figure 9:

Different from the traditional method, the concrete value of the comprehensive specific frictional resistance can be calculated for each pipe segment, rather than resulting from an arbitrary choice in the range of 30~70 Pa/m, so the method of comprehensive specific frictional resistance has a higher degree of refinement. In addition, the comprehensive specific frictional resistances are different for different pipe segments in a heating system, which reflects the pipe segment differences. Comprehensive specific frictional resistance decreases with an increase of design flow rate, and this decreasing trend slows down with the progressive increase in flow rate. Therefore, when the design flow rate of the pipe segment is in a certain range, the comprehensive specific frictional resistance is between 30~70 Pa/m, which is consistent with the traditional economic specific frictional resistance, but for the pipe segments beyond this range, it is unreasonable to determine the comprehensive specific frictional resistance according to the traditional method.

Comparing the two images in Figure 9, there is little difference in the comprehensive specific frictional resistance between the water-supply pipes and the corresponding water-return pipes, but we still find that a small difference may lead to different diameters, such as in pipe segments 115 and 115′ or 122 and 122′ in Table 2.

We know from Figure 9 that the comprehensive specific frictional resistance with flow rate regulation is bigger than it is with quality regulation for the same flow rate, and this difference increases with decreases in the design flow rate. Therefore, it is necessary to distinguish between different regulation modes regarding the issue of comprehensive specific frictional resistance.

For the pipe segments with universal flow rate, comparing with the range of traditional economic specific frictional resistance of 30~70 Pa/m, the value of comprehensive specific frictional resistance is higher and can reach up to 200 Pa/m. There are three reasons for this: ① the value of 30~70 Pa/m was obtained from experience with directly connected heating networks in the past [8]. For indirectly connected heating networks, according to the experience of Nordic countries, the value of specific frictional resistance will be higher than the above specified value; ② in recent years, the growth trend of electrovalence has been relatively flat, but the growth in pipe materials and pipe fittings has been large, resulting in an increase in specific frictional resistance; ③ for flow rate regulation, the calculation of comprehensive specific frictional resistance considers the change of the actual flow rate during operation. By this calculation, the computational electricity consumption is less than the electricity consumption calculated by design flow rate, so the value of comprehensive specific frictional resistance is greater than the traditional value.

4. Conclusions

In this paper, we put forward the concept of comprehensive specific frictional resistance, which aims to improve the precision degree of specific frictional resistance so as to promote the development of heating engineering. We provide a method to obtain the value of comprehensive specific frictional resistance by establishing a mathematical model, which is applicable to the optimization of pipe diameter in primary heating networks. Although the comprehensive specific frictional resistance has the disadvantage of a complex calculation process, it reflects the pipe segment differences, differentiates the different operation regulation modes, and takes into account the safety problems caused by pipe pressure. Therefore, it is of high practical value.

The ideal value of comprehensive specific frictional resistance is affected by the factors mentioned in this paper and also needs to be enriched and developed through further research.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

References

H. V. Larsen, B. Bøhm, and M. Wigbels, “A comparison of aggregated models for simulation and operational optimisation of district heating networks,” Energy Conversion and Management, vol. 45, no. 7-8, pp. 1119–1139, 2004.

H. Shi, The theory research of economic friction of heating pipe network [M.S. thesis], Xi'an University of Architecture and Technology, Xi'an, China, 2013.

C. Hong, “A brief analysis of the weighted mean specific frictional resistance,” District Heating, vol. 41, pp. 49–53, 1992.

A. Ge, C. Qiao, and Y. Zhao, “The relationship between range of heat-supply service and optimal specific frictional head loss,” Energy Conservation, vol. 35, no. 12, pp. 18–21, 2015.

M. Liu, R. Wu, and Q. Yang, “Transportation energy consumption and economical specific frictional resistance design of sewage heat pump system,” Energy Conservation Technology, vol. 32, no. 188, pp. 487–491, 2014.

Z. Wang, Z. Dong, and R. Wang, “Study of economic specific frictional resistance on hot water heating system for civil buildings,” Architectural Science, vol. 26, no. 10, pp. 31–34, 2010.

X. Li, P. Zhou, and M. Duan, “Calculation of optimal specific frictional head loss in condensate network,” Gas and Heat, vol. 26, no. 3, pp. 59–62, 2006.

P. He and G. Sun, Heating Project, China Architecture and Building Press, Beijing, China, 2000.

Ministry of Housing and Urban-Rural Construction of the People's Republic of China, State Technical Method for Civil Building Project Design, China Planning Press, China.

X. Sun, Research of the economic specific frictional resistance of primary network of hot water heating[M.S. thesis], Taiyuan University of Technology, Shanxi, China, 2010.

J. Liu, “Analysis and calculation of economic specific frictional resistance in hot water heating network,” Journal of Nanjing Institute of Architectural Engineering, vol. 15, no. 4, pp. 38–43, 1998.

M. Manfren, P. Caputo, and G. Costa, “Paradigm shift in urban energy systems through distributed generation: Methods and models,” Applied Energy, vol. 88, no. 4, pp. 1032–1048, 2011.

C. Y. Xu and Q. Chen, “An entransy dissipation-based method for global optimization of district heating networks,” Energy & Buildings, vol. 48, no. 1, pp. 50–60, 2011.

M. Bojic and N. Trifunovic, “Linear programming optimization of heat distribution in a district-heating system by valve adjustments and substation retrofit,” Building and Environment, vol. 35, no. 2, pp. 151–159, 2000.

J. Fu, Industrial Technology Economics, Tsinghua University Press, Beijing, China, 1996.

F. Wu and F. Ye, Engineering Economics, China Machine Press, Beijing, China, 2007.

S. Xu and Y. Jia, “The analysis of the operational regulation mode in the system of district heating,” Journal of Hebei University of Engineering, vol. 25, no. 2, pp. 36–39, 2008.

J. Chen, Research on energy saving of the operation adjustment mode of heating system [M.S. thesis], North China Electric Power University, Beijing, China, 2014.

B. Babiarz, “An introduction to the assessment of reliability of the heat supply systems,” International Journal of Pressure Vessels and Piping, vol. 83, no. 4, pp. 230–235, 2006.

H. Boudali and J. B. Dugan, “Erratum: "A continuous-time Bayesian network reliability modeling and analysis framework" (IEEE Transactions on Realibility),” IEEE Transactions on Reliability, vol. 57, no. 3, p. 532, 2008.

O. Fujiwara and A. U. De Silva, “Algorithm for reliability-based optimal design of water networks,” Journal of Environmental Engineering (United States), vol. 116, no. 3, pp. 575–587, 1990.

I. C. Goulter and A. V. Coals, “Quantitative approaches to reliability assessment in pipe networks,” Journal of Transportation Engineering, vol. 112, no. 3, pp. 287–301, 1986.

D. A. Savic and G. A. Walters, “Genetic algorithms for least-cost design of water distribution networks,” Journal of Water Resources Planning and Management, vol. 123, no. 2, pp. 67–77, 1997.

K. Vairavamoorthy and M. Ali, “Optimal design of water distribution systems using genetic algorithms,” Computer-Aided Civil and Infrastructure Engineering, vol. 15, no. 5, pp. 374–382, 2000.

X. Xia and L. Feng, “Analysis of the influence of absolute equivalent roughness,” Shanghai Gas, vol. 54, no. 2, pp. 10–12, 2010.

F. Chen, H. Yu, Z. Sun, and Y. Guo, “Determination of inner wall roughness of pipeline,” Nature Gas and Oil, vol. 25, no. 6, pp. 8–10, 2007.

Y. Dou and X. He, “Calculation method for hydraulic resistance of pipe,” Journal of Hydraulic Engineering, vol. 40, no. 7, pp. 28–31, 1995.

M. Xu and P. He, “Determination method of heating load continuity diagram,” District Heating, vol. 4, no. 1, pp. 3–10, 1986.

L. Yang and H. Tao, “Engineering experimental research on energy-saving of economical flow velocity and pipe diameter optimization in water supply system,” Water & Wastewater Engineering, vol. 35, pp. 334–336, 2009.

H.-B. Jiang, “Cost function, optimal pipe diameter and economical flow velocity of water supply system,” Harbin Gongye Daxue Xuebao/Journal of Harbin Institute of Technology, vol. 36, no. 1, pp. 84–90, 2004.

Z. Chen, Z. Cai, and C. Qian, Chemical Engineering Thermodynamics, Chemical Industry Press, China, 2015.

J. Chen, Real Function Theory, Science Press, China, 1978.