Investigation on Pressure Drop of Fluid-Solid Mixture Flow through Pipes Using CFD and SK Model

The carrier fluid (air or water) is used to transport solid material from the source place to its destination point through pipeline. Using air as carrier fluid to transport solid material through pipeline is called Pneumo transport, whereas transporting material with water or any other liquid through pipeline is called as hydraulic transport. A large number of installations are now available globally to transport solid materials to short, medium, and long distances using water/air as carrier fluid. However, the design of such system of pipeline is still an empirical art. In the present investigation, one generalized mathematical model developed by Shrivastava and Kar (SK Model) and CFD models were used and compared with experimental results for pneumatic and hydraulic transport of granular solids. The motivation of present work is to find the accuracy of SK model based on analytical, empirical and semiempirical for the prediction of pressure drop and comparing the result with CFD based on mathematical equation for the mixture flow in the horizontal and vertical pipe lines. The comparison of pressure drop results obtained by using SK model and CFD model were validated with the experimental results for pneumatic and hydraulic transport of solids through. From the comparison results, it was observed that the results of pressure drop predicted by SK model are more accurate than the CFD models for all the cases considered.


Introduction
Fluid-solid mixture flow is form of multiphase flow which is simultaneous flow of two or several phases through a system such as pipe in the present investigation.From the literature survey [1] [2] [3], two-phase flow may be grouped according to the phases involved as gas-liquid mixture, gas-solid mixture, liquid-solid mixture, and two-immiscible-liquids mixture.The study of two phase flow is not only important from the viewpoint of fundamental phenomenon but also in practical applications (pneumatic transport and hydro-transport of particles in pipes) and natural phenomenon (e.g.sediment transport in water bodies, biological/ biomedical flow).Solid particles can be transported from one location to another using pipeline and carrier fluid like liquid, water, oil, and etc. and gases.Depending up on the carrier fluid the solid transport system may be called as hydraulic transport or Pneumo transport, when water or air is used as carrier fluid respectively.
Many investigators have carried out experimental on prediction of pressure drop and minimum conveying velocity on fluid-solid mixture flow through pipes.Titus et al. [4] simulated the effect of particle size (90 μm to 270 μm) liquid slurry using Eulerian-Eulerian two-fluid model in ANSYS CFX-15 in a horizontal pipeline to examine frictional pressure loss.They observed that decrease in frictional pressure loss as particle size increased at constant particle volume fraction.Furthermore, for a constant particle volume fraction, the radial distribution of particle concentration increased with increasing particle size, where high concentration of particles occurred at the bottom of the pipe.Kelessidis et al. [5] have performed experimental investigation on the flow of dilute solid-liquid mixtures in a horizontal concentric and 100% eccentric annulus.The experiments were performed in a 5 m long, 7 cm by 4 cm transparent annulus.A range of superficial liquid velocities between 0 and 2.32 m/s were studied, The authors observed that main parameters affecting the flow pattern transitions are the liquid superficial velocity, the pipe eccentricity and, to a lesser extent, the liquid viscosity.Besides frictional pressure drop with dilute slurries flowing in the concentric annulus is higher compared to liquid-only flowing for all three fluids An experimental investigation related to frictional pressure loss in a horizontal pipeline on the effect of various particle size was examined by Nabil et al. [6].
From the results, least pressure loss was recorded in fine particle size whereas the highest pressure loss was recorded in the coarse particle size for each constant particle volume fraction ranges from 5% to 30%.Besides, it was found that the particle velocity profiles were asymmetric in nature and dependent on particle size and fraction of volume.
Investigation is done on hydraulic transport of coarse solid particles to demonstrate as to how a generalized mathematical model (SK model), developed earlier by Shrivastava and Kar [7] can be utilized in designing the systems of pipeline for the hydraulic transport of coarse solid particles.The SK model has Journal of Applied Mathematics and Physics been validated thoroughly with a large number of experimental results for conveying coarse solid particles such as sand, gravels, pebbles, coal, etc. through pipes of different orientations in earlier investigations.By utilizing the results of the earlier investigations along with some new ones it has been clearly shown in this paper that the SK model estimates the values of head loss and critical velocity quite accurately as compared to the experimental results for conveying coarse solid particles hydraulically through the pipes of different diameters and orientations, irrespective of the particle-density, shape, size, size distribution, and the solid concentration in the mixture.The SK model has also been found useful in optimizing the various parameters of hydraulic transport from energy consumption point of view [8].
For any fluid-solid transport system the most important parameters required to design the system are the pressure drop and minimum transport velocity.Estimation of pressure drop is very complicated phenomena in multiphase flow.
The techniques for analyzing two phase flow include correlations, phenomenological models, simple analytical model, integral analysis, differential analysis, computational fluid dynamics (CFD) and artificial neural network (ANN) [9].
The physics of two-phase flow is more complex than for single flow due to the presence of dispersed (solid) phase so that this work deals with investigation on fluid (gas/liquid)-solid mixture flow through pipes.Such two-phase flow is widely encountered in pipe flows.One example of the gas-solid mixture through pipe in textile factory (air-cotton ball) that found nearby Mikelle city, Ethiopia in eastern Africa is presented in the application of this investigation.The data used for the study is based on both primary and secondary sources.The primary data are collected through MAA Garment and Textile Factory at Mekelle city, Tigray region, Ethiopia.The secondary sources of data that the researcher used are from different relevant books and Journals.In this study an attempt has been made to apply and compare the two methods of prediction pressure drop i.e.
Shrivastava and Kar (SK model) [2] [7] and CFD with experimental found in literatures for the fluid-solid mixture flow through pipe for different orientations.
The objective of the study is to determine the most accurate empirical/semi empirical correlations or computational fluid dynamics (CFD) models which can predict the most accurate results for pressure drop estimation and can be applying the method in industries.

Theory of SK Model
The total head loss ( t h ) per unit pipe length, for conveying heterogeneous mix- ture of fluid and solid through an horizontal, vertical, inclined and bending pipe under fully-accelerated and suspended conditions of flow, is the summation of head loss ( f h ) due to fluid-pipe friction and additional head loss s h due to the presence of solid particles in the flow, Shrivastava and Kar [2] were able to obtain an expression for t h in their mathematical SK model [7], which is as given below: The Darcy-Weisbach friction factor (f) of the above Equation ( 1) may be determined by using the Blasius equations as given below: where, Re-Reynolds numbers defined as . The values of terminal velocity ( T V ) and drag coefficient ( d C ) of the solid particles to be transported with a carrier fluid through a pipe may be determined either experimentally or analytically by using the properties of the carrier fluid and the solid particles and taking the mean equivalent spherical diameter (d) of the particles in the following Equations ( 3) to ( 6): ( ) 0.6 18.5 , 0.1 500 where, p Re -Reynolds number of particles is defined as ( ) And Mass flow rate of the fluid is defined as: The mathematical expression developed for critical velocity ( C V ) in the SK model by applying the well-known techniques of optimization to Equation ( 1) is as given below: ( ) where values of constants K and C may be determined experimentally from the plots of friction factor (f) against Reynolds number (Re) for the flow of fluid alone in pipe, which results into a relation of the type: The value of density of particle ( p ρ ) can be determined by using the rela- tionship between mass of particle ( p m ) which is directly measured and volume of the particle as shown below equation: Journal of Applied Mathematics and Physics The value of density of fluid is expressed by the equation below: 3. Computational Fluid Dynamics

Multiphase Modeling
Advances in computational fluid mechanics have provided the basis for further insight into the dynamics of multiphase flows.Currently there are two approaches for the numerical calculation of multiphase flows: the Euler-Lagrange approach and the Euler-Euler approach.

Euler-Lagrange Approach
The fluid phase is treated as a continuum by solving the Navier-Stokes equations, while the dispersed phase is solved by tracking a large number of particles, bubbles, or droplets through the calculated flow field.The dispersed phase can exchange momentum, mass, and energy with the fluid phase.This approach is made considerably simpler when particle-particle interactions can be neglected, and this requires that the dispersed second phase occupies a low volume fraction, even though high mass loading ( ) is acceptable as reported by Bartosik and Shook [10].The particle or droplet trajectories are computed individually at specified intervals during the fluid phase calculation.This makes the model appropriate for the modeling of spray dryers, coal and liquid fuel combustion, and some particle-laden flows, but inappropriate for the modeling of liquid-liquid mixtures, fluidized beds, or any application where the volume fraction of the second phase cannot be neglected.

Euler-Euler Approach
Manoj Kumar and Kaushal [9] were approached Euler method to investigate Modeling of sand-water slurry flow through horizontal pipe using CFD technique.In the Euler approach, the different phases are treated mathematically as interpenetrating continua.Since the volume of a phase cannot be occupied by the other phases, the concept of phasic volume fraction is introduced.These volume fractions are assumed to be continuous functions of space and time and their sum is equal to one.Conservation equations for each phase are derived to obtain a set of equations, which have similar structure for all phases.These equations are closed by providing constitutive relations that are obtained from empirical information, or, in the case of granular flows, by application of kinetic theory.

Results and Discussions
The input data used for prediction pressure drop using the SK model are collected from the Manual of MAA garment and Textile factory [11], Mekelle city, Tigray region, Ethiopia as shown in Table 1.

Comparisons of Pressure Drops Using SK Model and CFD Model
The equations given in the SK model predicts pressure drops more accurately than any other correlations available.However, many work in this area has been done utilizing the numerical methods and CFD techniques during last five decades to predict the values of pressure drop for conveying solids through pipes using water or air as carrier fluids.Hence it was decided to compare the performance of SK model with that of the CFD models by comparing the calculated results of pressure drop with the experimental values of the investigators [12] [13] [14].
In this section the objective is to examine the accuracy of SK model in predicting the values of pressure drops by comparing the results with the experimental values and those obtained utilizing the CFD model for hydraulic and pneumatic conveying of granular solids.

Hydraulic Transport
Experimental investigations have been conducted by Tamer Nabil et al.    Figure 4 and Figure 5 clearly depict that the effect of particles size on pressure drop in the plot of pressure drop against velocity of water for volume fraction of 5% and 25% respectively for the transportation of sand hydraulically through horizontal pipe made by using SK model.
It may be concluded from the curves obtained using SK and CFD model that the pressure drop is inversely proportional to particle size, i.e pressure drop is minimum for highest size of particle and as the size reduces pressure drop also increases throughout the velocity of water spectrum.However it is evident from these comparisons that, the results of pressure drop obtained by SK model match very closely with experimental values than the predictions made by CFD model.
The deviation from the experimental in the case of SK model was found to be ±7.54%whereas for CFD model it was observed to be ±27% as clearly reflected in the graph.The Properties of PVC granules and water along with calculated values of C d and V t given in Table 5 as described.
In another work made with by Bartosik and Shook et al. [16] for hydraulic transport, a CFD model was developed and the model was tasted with experimental result for conveying of PVC granules of 3.4 mm size through a vertical

Pneumatic Transport
An exhaustive numerical simulation using the CFD methodology and taking Euler-Euler approaching was performed for developing CFD model by Pandaba Patro et al. [17], in their research work for gas-solid flow, and the model was validated by comparing the predicted values of pressure drop of the investigators for conveying plastic particles of 200 µm and 3400 µm sizes pneumatically through 30 mm pipe.
In the present investigation the SK model was applied to calculate pressure drops for conveying solids of different mass loading ratio pneumatically at different volumetric concentrations and results of pressure drop have been compared with the predicted values using CFD model and also with the experimental  for mass loading 1, 2 and 3 respectively.The properties of air and calculated value of C d and V t for plastic solid material given in Table 6 as described.
Besides the above investigations numerical simulation was also conducted by Brundaban Patro [18] with Euler-Euler approach and considering four way      results than the predicted values using CFD model.The properties of air and calculated value of C d and V t for 200 µm given in Table 7 as described.

Conclusions
• From the above experimental literatures for prediction of pressure drop, the deviation between the experimental result and calculated value is found to be good agreement as similar with Shrivastava Kaushal.The validation of experimental results of pressure drop by using of empirical and semi-empirical correlation of SK model was concluded as more accurate.• The SK model is good agreement for prediction of pressure drop as the volume concentration is less than 25%.Because this model is not considered [15] to measure the pressure drops for transporting sands of three different sizes (d 50 = 0.2 mm, 0.7 mm and 1.4 mm) at different volumetric concentrations with water through a horizontal pipes of diameter 26.8 mm.These results of pressure drop have been compared with the predicted values using CFD model and also the experimental results.All these comparisons of pressure drops obtained SK model, CFD model and experimental results have been plotted in Figures1-3for sands of three different sizes d 50 = 0.2 mm, 0.7 mm and 1.4 mm respectively.It may be observed from Figures1-3, the SK model results are more closely and it was matching with experimental results than the predictions made by using CFD model for all the sizes of sand and all volumetric concentrations.The properties Journal of Applied Mathematics and Physics

Figure 1 .
Figure 1.Comparison of estimated pressure drop with for d 50 = 0.2 mm.

Figure 2 .
Figure 2. Comparison of estimated pressure drop for d 50 = 0.7 mm.

Figure 3 .
Figure 3.Comparison of estimated pressure drop for d 50 = 1.4 mm.

Table 3 .Table 4 .
Properties of sand and water, and calculated value of C d and V t for d 50 = 0.7 mm.Properties of sand and water, and calculated value of C d and V t for d 50 = 1.4 mm.Hydro-transport (Sand-water mixture flow through Horizontal pipe) et al.DOI: 10.4236/jamp.2019.71018226 Journal of Applied Mathematics and Physics

Figure 4 .
Figure 4. Effect of particles size pressure drop against velocity of water for (C v ) 5%.

Figure 5 .Figure 7
Figure 5.Effect of particles size pressure drop against velocity of water for (C v ) 25%.Table 5. Properties of PVC granules and water, and calculated value of C d and V t .Hydro-transport (liquid-solid mixture flow through vertical pipes) 40 mm D =

Figure 6 .
Figure 6.Comparison of estimated pressure drop using SK model, CFD model with C v = 10%.

Figure 7 .
Figure 7.Comparison of estimated pressure drop using SK model, CFD model with C v = 20%.

Figure 13 .Table 7 .
Figure 13.Comparison of CFD and SK for Pneumo-transport (β = 3) with 200 µm.Table 6.Properties of air, and calculated value of C d and V t for plastic solid material.Pneumatic transport (Gas-solid) 30 mm D = 200 m d = µ

Table 1 .
Specifications of the Materials from MAA Garment.

Table 2 .
Properties of sand and water and calculated value of C d and V t for d 50 = 0.2 mm.