Numerical Simulation of Hydraulic Transport of Sand-Water Mixtures in Pipelines ()
1. Introduction
Solid-liquid flow in pipelines has been a popular mode of transportation in chemical and waste-disposal industries. Circular pipes are normally used for long distance transportation of variety of materials in bulk quantities [1]. Both time and pumping power required for hydraulically transporting huge amounts of soils in a state of slurry (sand-water mixture) are crucial factors for successful dredging and reclaiming projects [2]. Thus, the economical hydraulic transportation of solids to a long distance has been attempted by using several methods [3,4].
During the several past decades, the complex computational scheme describing the microscopic processes in the solid-liquid flow and the computing costs has been a drawback to the development and use of commercial softwares. Ling et al. [5] investigated the double slurry flow in the pipe for the fully turbulent flow using Eulerian granular multiphase (EGM) model. Al Araby et al. [6] performed a numerical study on the single phase combined free and force convection in the entrance region of a horizontal pipe with its wall temperature fixed as a constant value. Lin and Ebadian [7] used a simplified three-dimensional algebraic slip mixture (ASM) model [8] and the RNC k-e turbulent model [9] while focusing on the developing process of volume of fraction and density distributions, mean velocity profiles.
Kim et al. [10] experimentally investigated the hydraulic transport characteristics of sand-water mixtures in circular and square pipelines by changing the Reynolds number and volumetric delivered concentration.
In the present paper, the numerical computation results are presented by first validating the numerical results with the measured data. The present study is also mainly focused on the heterogeneous and homogeneous flow regimes. The effect of the concentration and the geometrical changes on the relative solid effect will be discussed.
2. Numerical Method
2.1. Governing Equations
The continuity equation for the mixture is
(1)
where denotes the volume fraction of single material q among the materials that constitutes the mixture,
is the density of a material q, is the mass flow rate transferred from the pth material to the qth material.
The momentum equation for the single material q can be expressed as
(2)
where the stress-strain tensor of a material q, , is
(3)
where and is the shear and bulk viscosity of the material q, respectively. denotes the external force, is the lift force and is the virtual mass force. represents the interaction among the materials, is the static pressure.
(4)
(5)
where are the coefficients that represent the momentum transfer between the fluid and the soild, and can be denoted as follows
(5)
The shear stress of the solid includes the shear viscosity due to the collision and transfer among the particles and the bulk viscosity, and can be written as follows:
(6)
where is the component due to collision, is the component due to motion, and is the frictional component. In the present study, the k-e turbulence model is used for the eddy viscosity modeling.
2.2. Boundary Conditions
In the present study, the mean velocity inlet boundary condition and pressure outlet boundary condition are imposed on the inlet and the outlet of the pipeline, respectively. The turbulence intensity level is set to 4% for the intermediate velocities in the sand-water mixture flows. The pipe length L was set to 100D. The density of water is 998.2 kg/m3. The specific gravity (Ss) and diameter (d50) of Jumunsin sands are 2.65 and 0.54 mm, respectively [10]. Water temperature is set to 20˚C.
A multi-block unstructured, non-uniform grid system with hexahedral elements is used to discretize the computation domain. Using the symmetry boundary condition, the computational domain is set to the half of the physical domain. Heat transfer is neglected with the steady state slurry flow assumption. The second-order upwind and central difference are selected as the discretization schemes of both convection and diffusion terms, respectively. With the SIMPLEC algorithm for the coupling between the velocity and pressure, under-relaxation technique is applied in the dependent variables.
3. Results and Discussion
In order to validate the present method computed results are compared with measured data. Figure 1 shows the velocity profiles from the bottom of the circular pipe to the top. Present results are compared with the measured data in [11]. The diameters of the pipe and the sand in the experiment are 51.5 mm and 0.165 mm, respectively. The volumetric delivered concentration is 8.4%. Though the present computation slightly over-predicts at the center and under-predicts at the point where the flow velocity has the maximum value, the present results agree with the measured data.
Figure 2 shows the change of the hydraulic gradient due to the changes in the Reynolds number and the volumetric delivered concentration. The amount of solids delivered inside the pipeline is represented by the volumetric delivered concentration that is calculated as the ratio between solid and slurry flow rates. The Reynolds number is defined as, where Dh and μw represent the hydraulic diameter of the pipes and the viscosity of the water, respectively. Mean velocity is a bulk velocity of the mixture defined as the volumetric flow rate of a matter passing through a pipeline, Qm, divided by the cross-sectional area, A.
(7)
where Qm equals the sum of Qs and Qf that represent the
Figure 1. Comparison of velocity profile at Cvd = 8.4%, d = 0.165 mm and D = 51.5 mm.
Figure 2. Concentration (Cvd) effect to hydraulic gradient (Im) at d = 0.54 mm and D = 200 mm.
flow rates of the solid particles and the fluid, respectively. In this paper, the definition of the Reynolds number is also used for water and slurry flows together. Present results are compared with the measured data [2]. The diameters of the pipe and the sand in the experiment are 150 mm and 0.54 mm, respectively. It can be seen from the figure that present results agree well with the measured data when the Reynolds are within the range from 0.7 × 106 to 1.3 × 106.
Figure 3 shows the relation of the volumetric delivered concentration with the hydraulic gradient in the given Reynolds number. The energy loss due to friction can be assessed by using hydraulic gradient that means the pressure gradient for a mixture expressed in height of water per length of line.
(8)
As shown in Figure 3, the effect of the volumetric delivered concentration on hydraulic gradient increases as the Reynolds number decreases. When the Reynolds number is small, the increase in the volumetric delivered concentration has an effect of decreasing the hydraulic gradient.
Figure 4 shows the relation of the volumetric delivered concentration with solid effect in the given Reynolds number. In evaluating additional friction loss due to the solid particles, solid effect represented by is used. Iw means the friction gradient for water alone at the flow rate equal to the mixture flow rate. Figure 4 shows that the effect of the volumetric delivered concentration on solid effect increases as the Reynolds number decreases. When the Reynolds number is small, the solid effect increases with the volumetric delivered concentration stepping up. Thus, it can be said that, when the flow speed is low, the friction among the solid particles becomes a dominant factor, which results in the increase in the energy loss due to friction. When the flow speed is high, the solid particle distributions become uniform and the solid particles flow through the pipe with less probability of settling on the bed. Thus, the friction between the solid particles with the wall or deposited solid particles becomes small.
Figure 5 shows the effect of the solid particle size on the relative solid effect as a function of the Reynolds number. It can be seen from the figures that the large size solid particles has the large relative
Figure 3. Concentration (Cvd) effect on both hydraulic gradient (Im) hydraulic gradient (Im).
Figure 4. Concentration (Cvd) effect on both hydraulic gradient (Im) and solid effect (Im − Iw).
(a)(b)(c)
Figure 5. Correlation of data with relative solid effect. (a) D = 100 mm; (b) D = 400 mm; (c) D = 600 mm.
solid effect. However, when the solid particle size is very small (d < 0.1 mm), the particle size does not affect the value of the relative solid effect. Wilson et al. [12] found the following relationship.
(9)
where 0.22 represents the value of the relative solid effects when Vm = V50. The magnitude of M is a function of normal stress of solids (granular pressure). In Equation (9), Sw = 1 because the working fluid is water. Following Equation (9), we can represent the relative solid effect as a function of the Reynolds number.
(10)
According to the estimation of Clift et al. [13], the value of M is set to 1.7. It can be said from the figures that the relative solid effect is affected more by the particle size than the pipe diameter.
The effect of the pipe diameter on the relative solid effect can be represented as follows [14]:
(11)
Table 1 and Figure 6 show the effect of non-dimensionalized particle size effect on the relative solid effect as a function of the Reynolds number. As shown on Table 1, the non-dimensionalized particle size is set to 0.0025 for the investigated 5 cases. The value of M is also set to 1.7, then the value of A is obtained by computing the relative solid effect with the change in the mean velocity. It can be seen from the Table 1 and Figure 6, that, even though the non-dimensionalized particle size has the same value, the relative solid effect has different values. In case the diameters of both pipe and solid particle are large, the values of the relative solid effect are larger than the other cases. Thus, it can be said from the results in Figure 6 that the size of the pipe diameter compared to the solid particle size is not the dominant factor. The particle size itself is a critical parameter to the performance of the slurry transportation system.
4. Conclusions
In this paper, the effect of the concentration, and both particle and pipe sizes on the relative solid effect is numerically investigated using Eulerian granular multiphase (EGM) model and k-e turbulent model. Present method is validated by the computed values with the measured data.
It is found that the effect of the volumetric delivered concentration on both hydraulic gradient and solid effect increases as the Reynolds number decreases. When the Reynolds number is small, the increase in the volumetric delivered concentration has an effect of decreasing the
Table 1. Correlation coefficient at d/D = 0.0025.
Figure 6. Correlation of data with relative solid effect (Im − Iw)/(Sm − Sw).
hydraulic gradient whereas the solid effect increases with the volumetric delivered concentration stepping up.
The effect of the particle size and the pipe diameter on the relative solid effect is calculated and correlated using the previously developed simple model. It is found that the effect of the pipe diameter or the relative size of the particle size with regard to the pipe diameter is not the critical parameter for deciding the values of the relative solid effect in the sand-water mixture transportation. The particle size itself is shown to be the most important parameter when the concentration is fixed.
In future, the effect of the particle size on the deposition limit velocity will be discussed.
NOTES