# CFD Study of Liquid Sodium inside a Wavy Tube for Laminar Convectors: Effect of Reynolds Number, Wave Pitch, and Wave Amplitude

Authors: Syed Murtuza Mehdi; Maaz Akhtar; Ahmad Hussain; Dheya Shuja Alothmany; Shahid Aziz
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

Metallic tubes have been widely used as primary heat transfer elements in laminar convectors for domestic and aerospace heating purpose. This paper uses CFD tool to investigate the heat output and pressure drop of liquid sodium flowing inside a circular tube having a wavy profile throughout its length. The wavy tube can be utilized in laminar liquid metal convectors as basic heat transfer element. The effect of Reynolds number (500≤Re≤2000) wave pitch (25mm≤λ≤100mm) and wave amplitude (2mm≤a≤6mm) on the heat output and pressure drop has been numerically studied. Based on the CFD results important controlling parameters have been identified and it is concluded that the heat output from the wavy tube is affected by the wave pitch and the wave amplitude while the pressure drop is mostly affected by the Reynolds number and wave amplitude.


---

## Body

## 1. Introduction

Laminar flows involve several attractive features like low thermodynamic and hydrodynamic irreversibility, silent and vibration-free operation and longer life and greater reliability, Doty et al. [1]. Another important feature of laminar flows is that the equations of motion (i.e., the Navier-Stokes equations for fluid dynamics) do not involve any time averaging technique like the turbulent flows, which ultimately makes the numerical results pertaining to this class of fluid motion directly and more reliably applicable to engineering designs. For instance, Kays [2] numerically studied the heat transfer in laminar flows inside circular tube, Ko [3] numerically studied the entropy generation and optimum Reynolds number for laminar developing flow inside double sine ducts, Chen and Chiou [4] studied the heat transfer in the entrance region of a laminar liquid metal pipe flow, Sahin and Mansour [5] numerically studied entropy generation in laminar fluid flow inside a circular pipe, Mehdi and Choi [6] performed CFD study of liquid sodium within thermally developing region under laminar flow, Gedik et al. [7] utilized CFD to study the laminar magnetohydrodynamics flow of liquid metal inside a circular tube, and Makinde and Eegunjobi [8] numerically studied the effects of convective heating on entropy generation for laminar flow inside a channel with permeable walls. It is important to mention that in laminar flows the convective effects are low which places limit on heat and mass transport; however, the transport of heat in laminar flows can be augmented by utilizing several alternatives which includes the use of low Prandtl number fluids such as liquid metals [6], the use of nanofluids as indicated by the work of Yang et al. [9], and increase in heat transfer area as indicated by the work of Kim et al. [10].Consequently laminar flows can become useful for domestic and aerospace heating applications. In fact, there exist many interesting heat exchanger design concepts that are based on laminar fluid flow; for instance, Stignor [11] presented the concept of laminar flat tube heat exchangers; Miner and Ghoshal [12] studied liquid metals as coolant for laminar microcoolers; similarly Tawk et al. [13] investigated different design parameters for microcoolers to be exploited in electronic cooling; Lin et al. [14] investigated the performance of a laminar heat pipe heat exchanger; Iqbal and Syed [15] presented the concept of laminar finned double-pipe heat exchanger; Saji et al. [16] developed a compact laminar microtube heat exchanger; Sagnik and Saha [17] presented the concept of utilizing simultaneously transverse corrugations and twisted tape inserts in a circular tube for augmenting the convective laminar heat transfer, whereas Sui et al. [18] investigated laminar flow and heat transfer in wavy microchannels.Based on the extensive use of laminar flows in different types of heat exchangers, here we present the concept of wavy tube as basic heat transfer element for laminar convectors which can be used in domestic and aerospace heating applications. In this paper computational fluid dynamics approach has been used to investigate the performance (heat output and pressure drop) of the wavy tube operating under fixed set of convective boundary condition for ambient air. The length of the tube is such that the flow always remains hydrodynamically developing. Also in this study liquid sodium is used as a heat transfer fluid since liquid metals due to their high molecular heat conduction can enhance the heat transport in laminar flows.


## 2. Fabrication of Wavy Tube

A possible fabrication technique for the wavy tube is depicted schematically in Figures1(a)–1(d). The process starts with splitting a tube of desired diameter in longitudinal direction. One half of the tube is then clamped within a set of die that is designed on the basis of desired wave parameters (wave pitch and amplitude). With the application of external force the desired profile can be impregnated on the tube. The other half of the tube is then processed in similar manner and finally the two halves can be joined together by using a suitable joining process such as fusion welding. The geometry of the tube with different waviness profile (i.e., different values of wave pitch and amplitude) was built using the readymade geometry features available in the graphic user interface (GUI) of the COMSOL Multiphysics software. Thus the different wavy profiles were generated in a way that the flow radius increases and decreases by an amount equal to the wave amplitude on either side of the base radius. The value of the base radius was kept constant at 13 mm. The shape of the waves remains like a portion of a smooth circular arc as seen in Figure 1(d).Figure 1
Schematic diagram depicting the conceptual steps for the fabrication of wavy tube. (a) Axially split tube of desired diameter. (b) Impregnating the desired wave length (λ) and wave amplitude (a) using a die. (c) Joining of the two parts and (d) final tube with desired wave length (λ) and wave amplitude (a).
(a) (b) (c) (d)


## 3. Problem Formulation

Steady hydrodynamically developing laminar flow of liquid sodium was studied inside a wavy tube for different values ofλ and a (see Figure 1(d)). The diameter of the tube was kept constant at 26 mm with an axial length of 0.55 m. The simulated Reynolds number range was 500 to 2000. It is well known that for steady laminar flows the hydrodynamic developing length can be estimated by the equation (Lh≈20ReD), where Lh is the hydrodynamic developing length of the flow, D is the diameter of the tube, and Re is the Reynolds number [19]. According to the above equation the developing length of the problem under consideration remains larger than the length of the tube considered in this study which means that the flow always remains hydrodynamically developing. Since the geometry under consideration is axisymmetric, the governing equations, that is, (1)–(3), were numerically solved for half of the computational domain using a finite element commercial code COMSOL Multiphysics.(1)ρV→·∇VV→=∇·-pI+μ∇V→+∇V→t-23μ∇·V→I+F,(2)∇·V→=0,(3)ρCV→·∇T-∇·k∇T=0.Equation (1) is the famous Navier-Stokes equation for fluid motion which solves for the velocity and pressure field within the flow domain, (2) is the continuity equation (conservation of mass) for an incompressible fluid, and lastly (3) is the energy equation used to calculate the temperature field. In (1) ρ is the fluid density, V→ is the velocity vector (composed of the axial and radial velocity components), ∇ is the del operator of vector calculus, p is the pressure of the fluid, I is the identity matrix, μ is the dynamic viscosity of the fluid, t indicates transpose of the matrix, and F is the body force such as gravity. Also in (3) C, T, and k represent the specific heat, temperature, and thermal conductivity of the working fluid, respectively. The different thermophysical properties of liquid sodium were estimated at the inlet temperature (T) in Kelvin from the following [20]:(4)ρ=1014-0.235T,C=-3.001×106T2+1658-0.8479T+4.454×10-4T2,k=104-0.047T,ln⁡μ=556.835T-0.3958ln⁡T-6.4406.After obtaining the complete solution (i.e., the velocity, pressure, and temperature field) the heat output and the pressure drop were calculated by using (5) and (6), respectively.(5)Qw=∫s-kdTdrwall,(6)Δpw=∫0R2πrpdr∫0R2πrdr.In (5) Qw represents the heat output from the wavy tube and is calculated as the surface integral of the radial conductive heat flux at the wall where r is the radial coordinate. Also in (6) Δpw represents the pressure drop along the wavy tube and p is the pressure distribution at the inlet of the tube.


## 4. Computational Details

COMSOL Multiphysics was used to solve the governing partial differential equations, that is, (1)–(3), within the computational domain. It is to be noted that the present problem is a coupled multiphysics fluid flow problem; therefore, direct PARDISO solver was used which treats the problem as coupled and highly nonlinear. The convergence criterion for all solution variables was set to 10−5. Both the stream line and cross wind stabilization techniques were used as these stabilization techniques are important for finite element based fluid dynamics simulations. The cross wind stabilization parameter was set to 0.1. Second-order Lagrange elements were used for the velocity and temperature variables while linear Lagrange elements were used for the pressure variable. The Lagrange elements approximate the solution variables with a piecewise polynomial function on each element of the grid. Further information on the finite element methods can be retrieved from [21, 22].For the present study, theR/38 (where R is the radius of the tube) criterion for laminar flows was used to design the grid [23]. The R/38 criterion states that the height of the wall adjacent element must be less than or equal to the ratio R/38. However, in order to properly resolve the boundary layer, approximately 5 elements were embedded within a small height of R/38 adjacent to the tube wall. The part of the grid near the wall is shown in Figure 2(a) which shows that presence of 5 elements within a height of 0.342 mm, while Figure 2(b) shows the overall and the magnified view of the grid used for the simulations. It can be noticed from Figure 2(b) that structured quadrilateral grid was used with grid density being larger near the solid wall. The grid independency of the final results (i.e., the heat output and the pressure drop) was confirmed by running the simulations at three different grid densities (4800, 9600, and 19200 elements). Figures 3(a) and 3(b) indicate that no substantial deviation exists between the three grids which highlights the grid independency of the computed results.Figure 2
CFD grid used for the simulations. (a) Number of elements and their height near the wall within a distance ofR/38. (b) Overall mesh and the enlarged view of the grid showing the distribution of the elements in the domain.
(a) (b)Figure 3
Grid independent results for (a) the heat output and (b) pressure drop.
(a) (b)For the fluid dynamics part of the problem, a uniform velocity based on the flow Reynolds number was defined at the inlet (V→=U=constant) of the tube. No slip condition (∂V→/∂r=0) was applied at the solid, nonporous, and stationary wall of the tube and zero gage pressure (i.e., pout = patm) was applied at the outlet. For the heat transfer part, at the inlet, a uniform temperature (T = 105°C) was defined. At the wall a convection condition (qconv′′=h∞(Tw-T∞)), where qconv′′ is the convective heat flux and Tw is the temperature of the wall, was assumed with a fixed heat transfer coefficient of h∞ = 10 W/m2K for still air and an ambient air temperature of T∞ = 25°C. At the outlet, convective flow was assumed which mathematically states zero temperature gradient at the exit boundary (∂T/∂z=0). The initial condition for the fluid flow in the entire domain was set to a uniform velocity that was equal to the inlet velocity and zero gage pressure, while for the heat transfer it was uniform temperature equal to the inlet temperature.


## 5. Results and Discussions

Figure4 indicates that at any given value of λ and a the increase in the heat output with the Re is not significant; therefore, at particular ambient conditions the wave pitch λ and amplitude a are the only important parameters that affect the heat output from the wavy tube. Consequently, Figure 5 shows that the normalized heat output increases if λ decreases and a increases. Decreasing the pitch and increasing the amplitude increase the effective length of the tube which means increase in the overall heat transfer area. This increase in the heat transfer area is the main cause of the enhanced heat output from the wavy tube. For the case of λ = 25 mm and a = 6 mm the heat transfer area increases by 48% which causes an increase of approximately 46% in the heat output as compared to the plane circular tube of same dimensions (see Figure 5). As an example of space heating, 50 wavy-peristaltic tubes having λ = 25 mm and a = 6 mm combined together in form of a laminar convector can increase the temperature of air inside 10×10×10m3 room from 10°C to 22°C in just 60 minutes. On the other hand decreasing λ and increasing a cause the pressure drop across the tube to increase. This effect is shown in Figure 6. From Figure 6 it can also be seen that, unlike the heat output, this time Re also has significant effect on the pressure drop. Figure 7 represents the typical streamline patterns within the diverging or bulge part of the wavy tube. The streamline pattern confirms the presence of circulation within the diverging part of the wavy tube. It can be observed that the size of the circulation zone increases as the wave amplitude a increases. The presence of circulation also indicates adverse pressure gradient effects which is also the reason of high overall pressure drop for the wavy tube. It is worth mentioning that streamline patterns indicating circulation within the diverging parts of the wavy tube were also identified in the works of Kim et al. [10] and Yang et al. [9].Figure 4
Heat output from the wavy tube as a function of Re,λ, and a.Figure 5
Normalized heat output from the wavy tube as a function ofλ and a.Figure 6
Pressure drop across the wavy tube as a function of Re,λ, and a.Figure 7
Typical streamline pattern within the diverging part of the wavy tube (a) atλ=25 mm and a=2 mm, (b) at λ=25 mm and a=4 mm, and (c) at λ=25 mm and a=6 mm.
(a) (b) (c)It is to be noted that high thermal conductivity of liquid sodium is capable of counterbalancing the low convective effects of a laminar flow; therefore, it is possible to operate wavy tubes at low Re in order to reduce the pressure drop across the tube and at the same time increase the heat output by reducingλ and increasing a. For instance, at λ = 25 mm and a = 6 mm, the pressure drop across the tube at Re = 2000 is 35 Pa; however, if Re is reduced to 500 the pressure drop reduces to 3.7 Pa (i.e., 89% reduction in the pressure drop) with only 3.4% reduction in the heat output of the tube. From Figure 6 it can also be noticed that compared to the wave pitch λ the wave amplitude a has a more prominent effect on the pressure drop.Based on the above results, it can be concluded that Re anda are the controlling parameters for optimizing the pressure drop across the wavy tube, whereas the heat output can be optimized by controlling λ and a independent of Re.Since, from Figure5, it is evident that the heat output increases if λ decreases and a increases, further simulations were performed to study the effect of heat transfer coefficient h∞ and the temperature of the ambient air T∞ on the exit temperature of sodium at Re = 500, λ = 25 mm, and a = 6 mm. It was found that with an inlet temperature of 105°C the exit temperature of liquid sodium always remains above its melting point of 97.8°C for all values of h∞≤10 W/m2K and T∞≥ 0°C. Therefore, Figure 8 depicts the variation of exit sodium temperature and the heat output as a function of T∞. From Figure 8 it can be concluded that the exit temperature of sodium remains higher and almost constant with the variation in T∞ whereas the heat output decreases as T∞ increases; however, this decrease is only less than 10% and therefore remains insignificant for the type of thermal system studied in this work.Figure 8
Variation of exit sodium temperature and heat output as a function of ambient air temperature ath∞=10 W/m2K, Re=500, λ=25 mm, and a=6 mm.


## 6. Conclusions

In summary, effect of Reynolds number Re, wave amplitudea, and wave pitch λ is studied numerically for flow of liquid sodium under laminar developing flow regime inside a wavy tube. It was found that the Reynolds number has insignificant effect on the heat output and can be ignored making the wave pitch and wave amplitude the controlling parameters for the optimization of the heat output. On the other hand, the pressure drop across the wavy tube is mostly affected by the Reynolds number and the wave amplitude. It is demonstrated that the wavy tube under the studied flow regime can be used as efficient heat transfer elements in laminar convectors for domestic and aerospace heating applications.



---

*Source: 6146195-2016-11-07.xml*