# Research on the Heating of Deicing Fluid in a New Reshaped Coiled Tube

Authors: Mengli Wu; Chiyu Wang; Yunpeng Li; Qi Nie
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

Aircraft ground deicing operation is significant to ensure civil flight safety in winter. Helically coiled tube is the important heat exchanger in Chinese deicing fluid heating system. In order to improve the deicing efficiency, the research focuses on heat transfer enhancement of deicing fluid in the tube. Based on the field synergy principle, a new reshaped tube (TCHC) is designed by ring-rib convex on the inner wall. Deicing fluid is high viscosity ethylene-glycol-based mixture. Because of the power function relation between high viscosity and temperature, viscosity has a negative influence on heat transfer. The number of ring-ribs and inlet velocity are two key parameters to the heat transfer performance. For both water and ethylene glycol, the outlet temperature rises when the number of ring-ribs increases to a certain limit. However, the increasing of velocity reduces heating time, which results in lower outlet temperature. The heating experiment of the original tube is conducted. The error between experiment and simulation is less than 5%. The outlet temperature of TCHC increases by 3.76%. As a result, TCHC efficiently promotes the coordination of velocity and temperature fields by changing the velocity field. TCHC has enhanced heat transfer of high viscosity deicing fluid.


---

## Body

## 1. Introduction

In the snow and ice weather situation, ground deicing is significant to the safety of the civil flight. In recent years, aircraft ground centralized deicing has been developed in the world. The main structure of ground deicing fluid heating system [1] made by China is shown in Figure 1. Helically coiled tube is the important heat exchanger in the system. High viscosity deicing fluid needs to be quickly heated in the tube. In order to improve the deicing efficiency and ensure the flight punctuality, it is important to study the heat transfer enhancement of deicing fluid.Figure 1
The structure of aircraft ground deicing fluid heating system.The deicing fluid is complex ethylene-glycol-based mixture. Its viscosity changes nonlinearly with the increasing of temperature, which obviously affects heat transfer. Most researches on deicing fluid are mainly about material properties [2, 3]. However, few papers about heat transfer of have been openly reported. The structure deformation is a commonly used and effective method to heat transfer enhancement. Aiming at rapidly heating deicing fluid, our group has proposed a new tube based on field synergy principle, which is called transverse corrugation helically coiled tube (TCHC).Heat transfer enhancement has always been a hotspot. Field synergy principle is a novel theory in this aspect, which was recently proposed by GUO, a Chinese scholar. The principle indicates that the intensity of heat transfer is related to not only the velocity and the fluid properties, but also the synergy degree between the velocity and temperature fields [4–6]. The heat transfer enhancement can be achieved by improving the coordination between velocity and temperature fields. Based on this principle, many researches have been carried on the heat transfer enhancement of straight tubes by changing the flow velocity field [7–10]. The principle has been widely tested and applied by many investigators.Many scholars have studied the flow properties inside the helically coiled tube. Dean [11, 12] gave out the Navier-Stokes function of flow in the helically coiled tube. Resulting from the difference of pressure of flow inside, the secondary flow was proved to improve the heat transfer efficiency. Since late 20th century, the Computational Fluid Dynamics (CFD) technique has developed quickly [13, 14]. Palanichamy and Nagaraj [15] presented a numerical simulation of laminar heat transfer in circular tube flows with internal longitudinal fins using finite difference numerical code. Tian et al. [16] investigated numerically the effects of six geometric parameters and pin-fin arrangement in the flow direction on the thermohydraulic performanceNumerical methods have been employed to analyze the influence of geometrical configurations and flow parameters in heat transfer of helically coiled tube. Lin and Ebadian [17] and Di Piazza and Ciofalo [18] studied the effect of the pitch and the curvature of the tube on convective heat transfer. The k-ε model is a two-equation turbulent model proposed by Jones and Launder in 1972 [19]. The k-ε model is the most widely used complete turbulence model. It is incorporated in most CFD codes [20]. The realizable k-ε model, proposed by Shih et al., has been validated in simulation of a variety of flows including rotating homogeneous shear flows, channel flows, and boundary flows [21, 22]. Hüttl and Friedrich [23] indicated that curvature ratio had a strong impact on the secondary flow. Zachár [24] studied the natural convection induced heat transfer over the outer surface of the tube. Different geometrical configurations and flow parameters were considered. A fluid-to-fluid boundary condition was set to realize the actual flow configurations. The results of the inner side heat transfer rate were tested by existing empirical formulas and experimental results.Aiming at the heat transfer enhancement of deicing fluid, TCHC is designed based on field synergy principle. As ideal fluid, water is studied for contrast on the influence of viscosity. Due to complex properties, the deicing fluid is simplified by its main component, ethylene glycol. The number of ring-ribs and inlet velocity are two key parameters to the heat transfer performance. These parameters’ variation of ethylene glycol and water is simulated in TCHC and the original tube. The outlet temperature is chosen as the evaluation index. The heating experiment is conducted to verify the simulation.


## 2. Numerical Method

For heat transfer enhancement of deicing fluid, original helically coiled tube has been reshaped to TCHC based on field synergy principle.


### 2.1. Geometrical Configurations

The ring-rib is equal-pitch convex on the inner wall of TCHC. The structure of ring-rib inside TCHC has been presented in Figure2. The flow velocity changes when passing by the ring-rib. Thus, an additional disturbance is generated to enhance heat transfer.Figure 2
Structure of ring-rib inside TCHC.In order to study how ring-rib affects heat transfer, three types of TCHC are obtained by changing the number of ring-ribs. The original tube is set as NO. 1 tube for contrast. TCHC with different numbers of ring-ribs are set as No. 2 to No. 4 tube. The models of the original tube and TCHC are shown in Figure3. The detailed parameters of four types of tubes are shown in Table 1.Table 1
Parameters of the original tube and TCHC.
Tube type
NO. 1
NO. 2
NO. 3
NO. 4
Number of ring-ribs/per round
0
4
8
16
Inner diameter/d (mm)
14
14
14
14
Outer diameter/D (mm)
20
20
20
20
Ring-rib height/h (mm)
0
2
2
2Figure 3
Geometric models of the original tube and TCHC. (a) NO. 1 tube, (b) NO. 2 tube, (c) NO. 3 tube, and (d) NO. 4 tube.
(a)
(b)
(c)
(d)


### 2.2. Grid Generation

Structural grid technology is used to mesh three different grid models in hexahedral structure. The grid generations of TCHC and the original tube are displayed in Figure4. The thickness of wall is simplified, so the section by which the fluid passes is meshed. To show the mesh of tube clearly, solid simple display is a method by making the hollow tube look like solid tube.Figure 4
The gird generation (solid simple display). (a) TCHC and (b) the original tube.
(a)
(b)The grid quantities of TCHC with different numbers of ring-ribs are given in Table2. The grid distribution is uniform which satisfies the requirement of numerical simulation. To make the calculation more accurate, the mesh of ring-rib is more compact than the other parts.Table 2
Grid quantity.
Tube type
NO. 1
NO. 2
NO. 3
NO. 4
Grid quantity (104)
76
93
112
134For the original tube and TCHC, the skew values of grid are tested. More than 97% of grid is within 0.4. The grid meshing meets the requirement of the calculation.


## 2.1. Geometrical Configurations

The ring-rib is equal-pitch convex on the inner wall of TCHC. The structure of ring-rib inside TCHC has been presented in Figure2. The flow velocity changes when passing by the ring-rib. Thus, an additional disturbance is generated to enhance heat transfer.Figure 2
Structure of ring-rib inside TCHC.In order to study how ring-rib affects heat transfer, three types of TCHC are obtained by changing the number of ring-ribs. The original tube is set as NO. 1 tube for contrast. TCHC with different numbers of ring-ribs are set as No. 2 to No. 4 tube. The models of the original tube and TCHC are shown in Figure3. The detailed parameters of four types of tubes are shown in Table 1.Table 1
Parameters of the original tube and TCHC.
Tube type
NO. 1
NO. 2
NO. 3
NO. 4
Number of ring-ribs/per round
0
4
8
16
Inner diameter/d (mm)
14
14
14
14
Outer diameter/D (mm)
20
20
20
20
Ring-rib height/h (mm)
0
2
2
2Figure 3
Geometric models of the original tube and TCHC. (a) NO. 1 tube, (b) NO. 2 tube, (c) NO. 3 tube, and (d) NO. 4 tube.
(a)
(b)
(c)
(d)


## 2.2. Grid Generation

Structural grid technology is used to mesh three different grid models in hexahedral structure. The grid generations of TCHC and the original tube are displayed in Figure4. The thickness of wall is simplified, so the section by which the fluid passes is meshed. To show the mesh of tube clearly, solid simple display is a method by making the hollow tube look like solid tube.Figure 4
The gird generation (solid simple display). (a) TCHC and (b) the original tube.
(a)
(b)The grid quantities of TCHC with different numbers of ring-ribs are given in Table2. The grid distribution is uniform which satisfies the requirement of numerical simulation. To make the calculation more accurate, the mesh of ring-rib is more compact than the other parts.Table 2
Grid quantity.
Tube type
NO. 1
NO. 2
NO. 3
NO. 4
Grid quantity (104)
76
93
112
134For the original tube and TCHC, the skew values of grid are tested. More than 97% of grid is within 0.4. The grid meshing meets the requirement of the calculation.


## 3. Governing Equations and Boundary Conditions

### 3.1. Conservation Equations

The basic equations of fluid mechanics include continuity equation, momentum equation, and energy equation, presented as (1) and (3) according to White [25]. On the basis of a large number of experiments, the main idea of CFD is to construct different models according to different hydrodynamic problems. Different numerical methods are applied to solve the related differential equations.Continuity Equation
(1)
∂
ρ
∂
t
+
∇
·
ρ
U
=
0
,where ρ is the density, ∇ is the Hamilton operator, and U is the velocity vector.Momentum Equation
(2)
d
d
t
ρ
U
=
-
∇
P
+
∇
·
τ
+
F
,where P is the pressure, τ is the viscous-stress tensor, and F is the additional volume stress.Energy Equation
(3)
d
d
t
ρ
u
+
P
∇
·
U
=
∇
·
k
d
∇
T
+
Φ
,where u is the internal energy, kd is the thermal conductivity, T is the temperature, and Φ is the viscous dissipation term.


### 3.2. The Turbulent Model

The turbulentk-ε model is a semiempirical formula summed up from the experiment. It determines the viscosity coefficient by solving the k equation and the ε equation; then the turbulence stresses are determined.The functions of standardk-ε model [19] are calculated as(4)∂ρk∂t+∂ρkui∂xi=∂∂xjμ+μtσk∂k∂xj+Gk+Gb-ρε-YM+Sk,∂ρε∂t+∂ρεui∂xi=∂∂xjμ+μtσε∂ε∂xj+C1εεkGk+C3εGb-C2ερε2k+Sε,where k is turbulent kinetic energy, ε is turbulent dissipation rate. C1ε, C2ε, and C3ε are constants: C1ε=1.44,  C2ε=1.92, and  C3ε=0.09. σk and σε possess the significance of turbulent Prandtl numbers of k and ε for the transport processes: σk=1.0 and  σε=1.3. ui is the mean velocity component. Sk and Sε are user-defined source items. μt is the turbulent eddy viscosity: μt=ρCμk2/ε, where Cμ is suggested as a constant having the value 0.09. μ is the dynamic viscosity. Gb is the turbulent kinetic energy due to buoyancy. Gk is the stress source term due to velocity gradient. YM is the pulsation expansion term in compressible turbulent flow.The functions of realizablek-ε model [21] are defined as(5)∂ρk∂t+∂ρkuj∂xi=∂∂xju+utσk∂k∂xj+Gk+Gb-ρε-YM+Sk,∂ρε∂t+∂ρεui∂xi=∂∂xjμ+μtσε∂ε∂xj+ρC1Eε-ρC2ε2k+νε+C1εεkC3εGb+Sε,where E=SijSij, Sij is the velocity strain rate tensor expressed as Sij=1/2(∂ui/∂xj+∂uj/∂xi), ν is the kinematic viscosity, and C1 and C2 are constants in numerical simulation.Thek equation in the realizable k-ε turbulence model is similar to the standard k equation. The constant of the turbulence viscosity formula is changed in the former to satisfy the realistic conditions. Besides, the source term in the ε equation is no longer related to Gk. The Reynolds stress term does not appear in the equation, which is good for representing the transfer of the energy spectrum between different scales.


### 3.3. Boundary Conditions

The boundary conditions of simulation have been shown as follows:(1)
The solver choice is a three-dimensional form based on pressure, implicit difference scheme, and activation energy equation.(2)
The realizablek-ε turbulence model is selected to carry out the numerical calculations, which is suitable for helically coiled tube.(3)
The near wall adopts the Standard-Wall-Functions treatment.(4)
The SIMPLE algorithm is adopted for the coupling of pressure and velocity.(5)
Relaxation factor of pressure is set as 0.3. Residual energy is set as10-6.(6)
Momentum, energy, and turbulent flow are applied as second-order upwind scheme.The entry condition is applied as velocity inlet. According to the actual working conditions, the reference inlet velocity is 12.45 m/s, which is calculated by the maximum flux of deicing fluid. Then, four different inlet velocities are set into simulation. The heating boundary conditions of the tube wall are adopted for uniform wall temperature, which is constant, 373 K. No wall-slip condition is adapted to wall surface. The outlet is set as pressure outlet. The temperature of work medium at entrance is set as 278 K.To investigate the influence of high viscosity deicing fluid on heat transfer in TCHC, water is simulated firstly for contrast. As complex ethylene-glycol-based mixture, deicing fluid is simplified by its main component, ethylene glycol. The model of UDF is used to compile the viscosity-temperature function of ethylene glycol. The other physical parameters needed are set as constants, like densityρ, specific heat capacity cP, and thermal conductivity kd. The material of tube is steel 310. The physical parameters are shown in Table 3.Table 3
The physical parameters of work medium and material of tube.
Work medium
Water
Ethylene glycol
Steel 310
Densityρ (kg⋅m−3)
998.2
1111.4
8030.0
Thermal conductivitykd (W⋅m−1⋅K−1)
0.600
0.252
16.270
Specific heat capacitycp (J⋅kg⋅K−1)
4182.00
2415.00
502.48
Dynamic viscosityμ (kg⋅m−1⋅s−1)
0.001003


## 3.1. Conservation Equations

The basic equations of fluid mechanics include continuity equation, momentum equation, and energy equation, presented as (1) and (3) according to White [25]. On the basis of a large number of experiments, the main idea of CFD is to construct different models according to different hydrodynamic problems. Different numerical methods are applied to solve the related differential equations.Continuity Equation
(1)
∂
ρ
∂
t
+
∇
·
ρ
U
=
0
,where ρ is the density, ∇ is the Hamilton operator, and U is the velocity vector.Momentum Equation
(2)
d
d
t
ρ
U
=
-
∇
P
+
∇
·
τ
+
F
,where P is the pressure, τ is the viscous-stress tensor, and F is the additional volume stress.Energy Equation
(3)
d
d
t
ρ
u
+
P
∇
·
U
=
∇
·
k
d
∇
T
+
Φ
,where u is the internal energy, kd is the thermal conductivity, T is the temperature, and Φ is the viscous dissipation term.


## 3.2. The Turbulent Model

The turbulentk-ε model is a semiempirical formula summed up from the experiment. It determines the viscosity coefficient by solving the k equation and the ε equation; then the turbulence stresses are determined.The functions of standardk-ε model [19] are calculated as(4)∂ρk∂t+∂ρkui∂xi=∂∂xjμ+μtσk∂k∂xj+Gk+Gb-ρε-YM+Sk,∂ρε∂t+∂ρεui∂xi=∂∂xjμ+μtσε∂ε∂xj+C1εεkGk+C3εGb-C2ερε2k+Sε,where k is turbulent kinetic energy, ε is turbulent dissipation rate. C1ε, C2ε, and C3ε are constants: C1ε=1.44,  C2ε=1.92, and  C3ε=0.09. σk and σε possess the significance of turbulent Prandtl numbers of k and ε for the transport processes: σk=1.0 and  σε=1.3. ui is the mean velocity component. Sk and Sε are user-defined source items. μt is the turbulent eddy viscosity: μt=ρCμk2/ε, where Cμ is suggested as a constant having the value 0.09. μ is the dynamic viscosity. Gb is the turbulent kinetic energy due to buoyancy. Gk is the stress source term due to velocity gradient. YM is the pulsation expansion term in compressible turbulent flow.The functions of realizablek-ε model [21] are defined as(5)∂ρk∂t+∂ρkuj∂xi=∂∂xju+utσk∂k∂xj+Gk+Gb-ρε-YM+Sk,∂ρε∂t+∂ρεui∂xi=∂∂xjμ+μtσε∂ε∂xj+ρC1Eε-ρC2ε2k+νε+C1εεkC3εGb+Sε,where E=SijSij, Sij is the velocity strain rate tensor expressed as Sij=1/2(∂ui/∂xj+∂uj/∂xi), ν is the kinematic viscosity, and C1 and C2 are constants in numerical simulation.Thek equation in the realizable k-ε turbulence model is similar to the standard k equation. The constant of the turbulence viscosity formula is changed in the former to satisfy the realistic conditions. Besides, the source term in the ε equation is no longer related to Gk. The Reynolds stress term does not appear in the equation, which is good for representing the transfer of the energy spectrum between different scales.


## 3.3. Boundary Conditions

The boundary conditions of simulation have been shown as follows:(1)
The solver choice is a three-dimensional form based on pressure, implicit difference scheme, and activation energy equation.(2)
The realizablek-ε turbulence model is selected to carry out the numerical calculations, which is suitable for helically coiled tube.(3)
The near wall adopts the Standard-Wall-Functions treatment.(4)
The SIMPLE algorithm is adopted for the coupling of pressure and velocity.(5)
Relaxation factor of pressure is set as 0.3. Residual energy is set as10-6.(6)
Momentum, energy, and turbulent flow are applied as second-order upwind scheme.The entry condition is applied as velocity inlet. According to the actual working conditions, the reference inlet velocity is 12.45 m/s, which is calculated by the maximum flux of deicing fluid. Then, four different inlet velocities are set into simulation. The heating boundary conditions of the tube wall are adopted for uniform wall temperature, which is constant, 373 K. No wall-slip condition is adapted to wall surface. The outlet is set as pressure outlet. The temperature of work medium at entrance is set as 278 K.To investigate the influence of high viscosity deicing fluid on heat transfer in TCHC, water is simulated firstly for contrast. As complex ethylene-glycol-based mixture, deicing fluid is simplified by its main component, ethylene glycol. The model of UDF is used to compile the viscosity-temperature function of ethylene glycol. The other physical parameters needed are set as constants, like densityρ, specific heat capacity cP, and thermal conductivity kd. The material of tube is steel 310. The physical parameters are shown in Table 3.Table 3
The physical parameters of work medium and material of tube.
Work medium
Water
Ethylene glycol
Steel 310
Densityρ (kg⋅m−3)
998.2
1111.4
8030.0
Thermal conductivitykd (W⋅m−1⋅K−1)
0.600
0.252
16.270
Specific heat capacitycp (J⋅kg⋅K−1)
4182.00
2415.00
502.48
Dynamic viscosityμ (kg⋅m−1⋅s−1)
0.001003


## 4. Simulation with Water

On the basis of actual working conditions, the inlet velocities are set as 5 m/s, 7.5 m/s, 10 m/s, 12.5 m/s, and 15 m/s to calculate the outlet temperature of TCHC and the original tube. As ideal fluid, water is chosen as contrast to deicing fluid on the influence of viscosity. The outlet temperature is selected as the evaluation index. When outlet temperature is higher, the heat transfer performance is better.


### 4.1. Number of Ring-Ribs

Considering how the number of ring-ribs affects heat transfer, the outlet temperature of water(Tw) is analyzed in TCHC and the original tube under different inlet velocities. The results of heat transfer process are summarized in Figure 5.Figure 5
The outlet temperature of water in TCHC with different numbers of ring-ribs.As seen in Figure5, all the outlet temperatures of NO. 2, NO. 3, and NO. 4 tubes are higher than that of NO. 1 tube. When the number of ring-ribs increases, the temperature rise tends to flatten out. Consequently, outlet temperature cannot keep rising when the number of ring-ribs exceeds a certain limit. Ring-rib changes the flow state by increasing the radial velocity. The temperature field basically remains radial. Based on field synergy principle, ring-rib can enhance heat transfer by improving the coordination between the velocity and temperature fields.


### 4.2. Inlet Velocity

From the above discussion, the outlet temperature is affected by inlet velocity. To investigate the influence of different velocities, the outlet temperature differences of water(ΔTw) are displayed in Figure 6.(6)ΔTw=Twi-Twj,where i, j are the numbers of different inlet velocities.Figure 6
The outlet temperature difference trend of water at different inlet velocities.Seen from Figure6, ΔTw of the NO. 1 tube is higher than those of TCHC, especially at low inlet velocity (<10 m/s). To the same tube, when the inlet velocity increases, the outlet temperature gets lower. Water takes less time to pass through the tube. The shortened heating time directly affects the outlet temperature.Δ
T
w of all tubes tends to be flat within a high velocity range in Figure 6. ΔTw of TCHC changes not as greatly as that of the original tube. Because of ring-rib, TCHC is less affected by the inlet velocity compared with the original tube.


## 4.1. Number of Ring-Ribs

Considering how the number of ring-ribs affects heat transfer, the outlet temperature of water(Tw) is analyzed in TCHC and the original tube under different inlet velocities. The results of heat transfer process are summarized in Figure 5.Figure 5
The outlet temperature of water in TCHC with different numbers of ring-ribs.As seen in Figure5, all the outlet temperatures of NO. 2, NO. 3, and NO. 4 tubes are higher than that of NO. 1 tube. When the number of ring-ribs increases, the temperature rise tends to flatten out. Consequently, outlet temperature cannot keep rising when the number of ring-ribs exceeds a certain limit. Ring-rib changes the flow state by increasing the radial velocity. The temperature field basically remains radial. Based on field synergy principle, ring-rib can enhance heat transfer by improving the coordination between the velocity and temperature fields.


## 4.2. Inlet Velocity

From the above discussion, the outlet temperature is affected by inlet velocity. To investigate the influence of different velocities, the outlet temperature differences of water(ΔTw) are displayed in Figure 6.(6)ΔTw=Twi-Twj,where i, j are the numbers of different inlet velocities.Figure 6
The outlet temperature difference trend of water at different inlet velocities.Seen from Figure6, ΔTw of the NO. 1 tube is higher than those of TCHC, especially at low inlet velocity (<10 m/s). To the same tube, when the inlet velocity increases, the outlet temperature gets lower. Water takes less time to pass through the tube. The shortened heating time directly affects the outlet temperature.Δ
T
w of all tubes tends to be flat within a high velocity range in Figure 6. ΔTw of TCHC changes not as greatly as that of the original tube. Because of ring-rib, TCHC is less affected by the inlet velocity compared with the original tube.


## 5. Simulation with Deicing Fluid

### 5.1. Viscosity-Temperature Fitting of Deicing Fluid

Deicing fluid is complex ethylene-glycol-based mixture. The thermophysical properties change nonlinearly with temperature. Thus, deicing fluid is simplified by the main component, ethylene glycol. Referring to the experiment, the variation between viscosity and temperature is demonstrated in Table4.Table 4
The viscosity with the changing temperature of ethylene glycol.
Temperature (K)
267
273
289
293
313
333
373
Viscosity (mPa⋅s)
86.9
53.5
25.7
22.1
10.2
5.2
1.7Considering different function forms, the fitting of power function is most consistent with the actual point. The viscosity of ethylene glycol decreases nonlinearly with the increasing of temperature. The fitting curve is displayed in Figure7. Equation (7) is derived from the relationship between viscosity and temperature.Figure 7
The fitting of ethylene glycol’s viscosity-temperature.The power function of fitting is calculated as(7)μ=8×1029·T-11.56.


### 5.2. Number of Ring-Ribs

The outlet temperature of ethylene glycol(Te) is presented in Figure 8. Te gradually rises when the number of ring-ribs increases. Compared with Figure 5, Te rises less than Tw in general. High viscosity is inversely proportional to heat transfer because it thickens the boundary layer of fluid. When the inlet velocity increases, Te generally decreases because the heating time is shortened. Moreover, for both ethylene glycol and water, TCHC can achieve the heat transfer enhancement.Figure 8
The outlet temperature of ethylene glycol in TCHC with different numbers of ring-ribs.


### 5.3. Inlet Velocity

The temperature difference of ethylene glycol(ΔTe) changes with inlet velocity. The tendencies of ΔTe are shown in Figure 9.Figure 9
The outlet temperature difference trend of ethylene glycol at different inlet velocities.Similarly to water,ΔTe can be calculated as(8)ΔTe=Tei-Tej.Seen in Figure9, because of the power function relation between high viscosity and temperature, ΔTe varies greatly more than ΔTw compared with Figure 6 (especially 7.5 m/s–10 m/s). ΔTe of all tubes tends to be flat within a high velocity range.


## 5.1. Viscosity-Temperature Fitting of Deicing Fluid

Deicing fluid is complex ethylene-glycol-based mixture. The thermophysical properties change nonlinearly with temperature. Thus, deicing fluid is simplified by the main component, ethylene glycol. Referring to the experiment, the variation between viscosity and temperature is demonstrated in Table4.Table 4
The viscosity with the changing temperature of ethylene glycol.
Temperature (K)
267
273
289
293
313
333
373
Viscosity (mPa⋅s)
86.9
53.5
25.7
22.1
10.2
5.2
1.7Considering different function forms, the fitting of power function is most consistent with the actual point. The viscosity of ethylene glycol decreases nonlinearly with the increasing of temperature. The fitting curve is displayed in Figure7. Equation (7) is derived from the relationship between viscosity and temperature.Figure 7
The fitting of ethylene glycol’s viscosity-temperature.The power function of fitting is calculated as(7)μ=8×1029·T-11.56.


## 5.2. Number of Ring-Ribs

The outlet temperature of ethylene glycol(Te) is presented in Figure 8. Te gradually rises when the number of ring-ribs increases. Compared with Figure 5, Te rises less than Tw in general. High viscosity is inversely proportional to heat transfer because it thickens the boundary layer of fluid. When the inlet velocity increases, Te generally decreases because the heating time is shortened. Moreover, for both ethylene glycol and water, TCHC can achieve the heat transfer enhancement.Figure 8
The outlet temperature of ethylene glycol in TCHC with different numbers of ring-ribs.


## 5.3. Inlet Velocity

The temperature difference of ethylene glycol(ΔTe) changes with inlet velocity. The tendencies of ΔTe are shown in Figure 9.Figure 9
The outlet temperature difference trend of ethylene glycol at different inlet velocities.Similarly to water,ΔTe can be calculated as(8)ΔTe=Tei-Tej.Seen in Figure9, because of the power function relation between high viscosity and temperature, ΔTe varies greatly more than ΔTw compared with Figure 6 (especially 7.5 m/s–10 m/s). ΔTe of all tubes tends to be flat within a high velocity range.


## 6. Discussions by Field Synergy Principle

Field synergy principle is a novel theory for heat transfer enhancement. The coordination between velocity field and temperature field has a positive influence on the heat transfer enhancement [5, 6]. Equation (9) gives a more general insight on convective heat transfer. There are two ways to enhance heat transfer: (a) increasing Reynolds and/or Prandtl number; (b) increasing the value of the dimensionless integration [4].(9)Nux=RexPr⁡∫01U¯·∇T¯dy¯.U
¯
·
∇
T
¯ in the dimensionless integration in (9) can be expressed as(10)U¯·∇T¯=U¯T¯cos⁡β,where β is the included angle between the velocity vector and the temperature gradient (heat flow vector).Figure10 sketches the outlet velocity fields of NO. 1 and No. 4 tubes at the same inlet velocity. The movement in the NO. 4 tube is more obvious than that in NO. 1 tube. Because of the ring-rib, the fluid movement is strengthened from wall to center in TCHC. Ring-rib changes velocity field by breaking the thickness of the flow boundary layer. As seen in Figure 10, the flow in TCHC mixes well, so the temperature distribution becomes more uniform. TCHC has efficiently improved the coordination of velocity and temperature fields. For this reason, TCHC can achieve the heat transfer enhancement.Figure 10
The outlet velocity field. (a) NO. 1 tube and (b) NO. 4 tube.
(a)
(b)


## 7. Heating Experiment

The heating experiment of the original tube is conducted. The experimental results are demonstrated in Figure11. At steady working condition, the outlet temperature (Texp) is basically stable at 350 K.Figure 11
Heating experiment of the original tube.The error (δ) between experiment and simulation can be expressed as(11)δ=Tw-TexpTexp×100%.Through calculation,δ is less than 5%; thus the experiment verifies the correctness of numerical methods. At the same working conditions, the outlet temperature of TCHC increases by 3.76% compared with the experimental data. Therefore, TCHC has better heat transfer performance than the original tube.


## 8. Conclusion

The original helically coiled tube is reshaped into TCHC based on field synergy principle. The heat transfer of deicing fluid is analyzed in TCHC and the original tube. Water is selected as a contrast to deicing fluid on the influence of viscosity. The outlet temperature is chosen as the evaluation index. The number of ring-ribs and inlet velocity are two key parameters to the heat transfer performance. The results are summarized as follows:(1)
Both number of ring-ribs and inlet velocity can affect heat transfer greatly. For water and ethylene glycol, the outlet temperature rises when the number of ring-ribs increases to a certain limit. However, the increasing of velocity reduces heating time, which results in lower outlet temperature.(2)
High viscosity is inversely proportional to heat transfer because viscosity thickens boundary layer of fluid. The outlet temperature rise of ethylene glycol is not as obvious as that of water.(3)
Both water and ethylene glycol in TCHC are more efficiently heated compared with the original tube. Ring-rib changes the velocity field based on field synergy principle. In consequence, TCHC improves the coordination between velocity and temperature fields.(4)
The error between experiment and simulation is less than 5%. The outlet temperature of TCHC increases by 3.76% compared to that of the original tube.In general, TCHC has enhanced heat transfer of high viscosity deicing fluid. Therefore, TCHC is beneficial to deicing efficiency and flight punctuality.



---

*Source: 3254631-2017-12-31.xml*