# Analytical study of Bingham fluid flow through a conical tube/Analitine binghamo skyscio srauto tekejimo kuginiu vamzdziu studija.

1. IntroductionThe blood vessel can be idealized as a tube with tissue space as circular porous bed (Vide Guha and Chuadhury, 1985).The rotating viscometer data of Rand et al. [1], Bugliarello et al. [2] and Chien et al. [3] suggest the non-Newtonian behaviour of blood. Lew et al. [4] suggested chyme as a non-Newtonian material having plastic-like properties. In view of this the biofluid flow in a living body can possess the non-Newtonian behaviour in general. In order to have a better understanding of blood flow in arteries and veins and chyme flow in stomach. It is necessary to consider the biofluid to be an yield stress fluid. One of the models for yield stress fluids is Bingham model. Hence the study of Bingham fluid flow through a conical tube with permeable wall is of considerable importance in medicine.

Bird et al [5] investigated the Bingham fluid flow in a rigid circular tube. Rathy [6] studied the flow of a Bingham fluid in a channel and in an annulus with impermeable walls. Vajravelu et al. [7] made a study on the Bingham fluid in a circular tube with permeable wall. The velocity field is obtained using Beavers and Joseph [8] slip condition at the permeable wall. Buckingham-Reiner equation for the flow is obtained.

The Bingham fluid flow between two permeable beds is discussed by Goverdhan et al. [9].The flow in the channel is assumed to be governed by Bingham model. The flow in the permeable beds is governed by Darcy's law. The velocity distribution is obtained. Some results are deduced and discussed. Comparini [10] discussed a one-dimensional model for the time dependent flow of a Bingham fluid between two parallel plates. The global existence and uniqueness of classical solution to the problem is proved.

Ravana et al. [11] studied the free surface flow of a Bingham fluid in an inclined channel over a permeable bed. The flow in the channel is described by Bingham model, whereas the flow in the permeable bed is according to Darcy's law. The velocity field, the shear stress, the mass flow rate and its fractional increase are obtained.

The problem of rotational motion of a Bingham fluid in the gap between two coaxial cylinders, the outer one being at rest and the inner one moving at given angular velocity is solved by Comparini [12].

Narahari [13] discussed unsteady flow of a Bingham fluid between two permeable beds having different permeabilities. The velocity distribution in the porous and non-porous regions are obtained. Some deductions are made and the results are discussed.

Sankara Reddy et al. [14] made a detailed study on the Bingham fluid flow in an inclined channel bounded by two permeable beds. The velocity distribution in the porous and non-porous regions are obtained. The temperature variation with Ec.Pr is discussed.

Viswanatha Reddy et al. [15] made a study on Bingham fluid flow in an annulus. The velocity field, the mass flow rate and its fractional increase are obtained. The results are deduced and discussed. Helical flow of a powerlaw fluid in a thin annulus with permeable walls is investigated by Vajravelu et al. [16]. It is observed that velocity increases due to permeable nature of the annulus.

In this paper, Bingham fluid flow through a conical tube with permeable wall is investigated. The velocity distribution, the volume rate of flow and its fractional increase are obtained. The results are deduced and discussed through graphs.

2. Mathematical formulation

Consider the flow of a Bingham fluid through a conical tube of length L with permeable wall. The flow takes place due to pressure gradient, and the porous medium is homogeneous with permeability k. The flow surrounded by the porous medium is governed by the Bingham model, and the flow in the porous medium is governed by the Darcy's law.

The tube has a radius [R.sub.0] at the entrance and radius [R.sub.L] at the exit (Fig. 1). The tube radius at any distance z from the inlet is given by

R(z) = [R.sub.0] + ([R.sub.L] - [R.sub.0]/L) z. (1)

The flow is axi-symmetric. Cylindrical polar coordinate system is used. The following assumptions are made in deriving the governing equations:

--the flow is steady and incompressible;

--the flow is in axial direction;

--all physical quantities except the pressure are function of z only;

--the body forces are negligible.

In view of the above assumptions, the governing equations and boundary conditions of the flow take the following form:

2.1. Governing equation

1/r d/dr (r[[tau].sub.rz]) = -(dp/dz), (2)

where

[[tau].sub.rz] = -[[mu].sub.0] [dv.sub.z]/dz + [[tau].sub.0]. (3)

2.2. Boundary conditions

[[tau].sub.rz] = finite at r = 0; (4)

[dv.sub.z]/dr = 0 at r = [R.sub.p] (5)

[v.sub.z] = [u.sub.B] at r = R; (6)

[dv.sub.z]/dr = [alpha]/[square root of k] ([u.sub.B] - Q) at r = R, (7)

where

Q = -[k/[mu][dp/dz]. (8)

By integrating equation Eq. 2 and using the boundary condition (4), we obtain:

[[tau].sub.rz] = - (dp/dz) r/2 a. (9)

Substituting (9) in (3), we get:

[[tau].sub.0] - [[mu].sub.0] [dv.sub.z]/dz = - (dp/dz) r/2, (10)

where [[tau].sub.0], [[tau].sub.rz], [alpha], k , [u.sub.B], Q, [[mu].sub.0], [u.sub.av], [v.sub.z], p and [R.sub.p] are the yield stress, Shear stress, slip parameter, permeability, slip velocity, Darcy velocity, coefficient of viscosity at the interface r = 0, average velocity, velocity component in z-direction, pressure and radius of the plug flow region.

2.3. Dimensionless equations

The flow physical parameters are made dimensionless as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

2.3.1. Radius equation

Eq. (1) takes the form (* are omitted here after):

R = 1 + ([R.sub.L] - l)z. (11)

2.3.2. Governing equation

Eq. (10) takes the form (* are omitted here after):

[[tau].sub.0]/r - [1/r][[dv.sub.z]/dr] = P/2, (12)

where P = - dp/dz.

2.3.3. Dimensionless boundary conditions

Eqs. (4)-(8) take the form (* are omitted here after):

[dv.sub.z]/dr = 0 at r = [R.sub.p]; (13)

[v.sub.z] = [u.sub.B] at r = R; (14)

[dv.sub.z] = [alpha]/[square root of k] ([u.sub.B] - [D.sub.a]P) at r = R, (15)

where [D.sub.a] = k/[R.sup.2.sub.0].

3. Solution of problem

Solving (12) and using the boundary condition (14) we obtain the velocity field around the plug flow region as:

[v.sub.1] = - P/4 ([r.sup.2] - [R.sup.2]) + [[tau].sub.0] (r - R) + [u.sub.B], (16)

when [R.sub.p] [less than or equal to] r [less than or equal to] - R.

Using (15) in (16), we get the slip velocity at the porous wall as:

[u.sub.B] = P[square root of ([D.sub.a]/[alpha])] [[alpha] [square root of [D.sub.a]] - R/2 + [[tau].sub.0]/P]. (17)

Using (13) in (16), we obtain the relation between [[tau].sub.p] and [R.sub.p] as:

[R.sub.p] = 2[[tau].sub.0]/P. (18)

Taking = r = [R.sub.p] in Eq. (16) and using the relation (18), we get the velocity field in the plug flow region as:

[v.sub.2] = -P/4 [(R - [R.sub.p]).sup.2] + [u.sub.B], (19)

when 0 [less than or equal to] r [less than or equal to] [R.sub.p].

The flow in the porous region (r [greater than or equal to] [R.sub.p]) is governed by the Darcy law which is given by:

Q = P[D.sub.a]. (20)

4. Deductions

Case (a): when R (z) = 1 and the permeability parameter k [right arrow] 1 in the Eqs. (16), (17) and (19) we get the velocity field outside the plug flow region r [greater than or equal to] [R.sub.p]:

[v.sub.1] = - P/4 ([r.sup.2] - l) + [[tau].sub.0] (r - 1) + [u.sub.B], (21)

and that in the plug flow region r [less than or equal to] [R.sub.p].

[v.sub.2] = - P/4 [(1 - [R.sub.p]).sup.2], (22)

which are in agreement with those of Bird et al (1960).

Case (b): when R (z) = 1, the Eqs. (16), (17) and

(19) we get the velocity field for the flow of a Bingham fluid flow through a circular pipe with permeable wall as follows:

[v.sub.1] = - P/4 [r.sup.2] + [[tau].sub.0]r + [u.sub.B] - [[tau].sub.0], (23)

when [R.sub.p] [less than or equal to] r [less than or equal to] 1;

[u.sub.B] = P[square root of [D.sub.a]]/[alpha] [[alpha] [square root of [D.sub.a]] - 1/2 + [[tau].sub.0]/P]; (24)

[v.sub.2] = P/4 [(1 - [R.sub.p]).sup.2] + [u.sub.B], (25)

when r [greater than or equal to] [R.sub.p]. These equations agree with those of Vajravelu and Sreenadh (1987).

Case (c): when R (z) = 1 as [[tau].sub.0] [right arrow] 0, the velocity field for the flow with permeable wall is given by:

[v.sub.z] = - P/4 [r.sup.2] + [u.sub.B] + P/4, (26)

when [R.sub.p] [less than or equal to] r [less than or equal to] 1.

[u.sub.B] = P [square root of [D.sub.a]]/[[alpha] [[alpha][square root of [D.sub.a]] [[alpha] [square root of [D.sub.a]] - 1/2]. (27)

These results are in a good agreement with those of Sreenadh and Arunachalam (1983).

5. Volume rate flow

The Volume rate of flow for the flow of a Bingham fluid flow through a conical tube is:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (28)

where [B.sub.1] = [R.sup.4] + 3[R.sup.4.sub.P] - 4R[R.sup.3.sub.p] and [B.sub.2] = [R.sup.3] - 2[R.sup.3.sup.p] + 3R[R.sup.2.sub.p].

When the permeability parameter k [right arrow] 0 (i.e., [D.sub.a] [right arrow] 0) (28) reduces to:

[Q.sup.1] = p[pi]/8 ([B.sub.1]) - [pi][[tau].sub.0]/3 ([B.sub.2]), (29)

where [B.sub.1] = [R.sup.4] + 3[R.sup.4.sub.p] - 4R[R.sup.3.sub.p] and [B.sub.2] = [R.sup.3] - 2[R.sup.3.sub.p] + 3R[R.sup.2.sub.p].

6. Fractional increase

The fractional increase in the volume rate of flow of the Bingham fluid through a conical tube with permeable wall over what it would be if the wall of the tube were impermeable is:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (30)

where [B.sub.1] = [R.sup.4] + 3[R.sup.4.sub.p] - 4R[R.sup.3.sub.p] and [B.sub.2] = [R.sup.3] - 2[R.sup.3.sub.p] + 3R[R.sup.2.sub.p].

7. Discussion of the results

The velocity profiles are shown in Figs. 2 - 12 for different values of [[tau].sub.0], Darcy number [D.sub.a], P = 10 and [alpha] = 1. It is observed that the velocity attains maximum value at r = 0 and decreases with the increment in r.

For fixed r and [D.sub.a], the velocity decreases due to increase in z. The slip velocity at the permeable wall decreases with the increment in the value of [D.sub.a]. On comparing the velocity profiles for permeable and impermeable conical tubes, it is found that the velocity is enhanced due to the permeability of the wall of the conical tube. It is also observed that the velocity remains constant from the axis (i.e., r = 0) up to some value of r (which is the Plug flow region) and then decreases to a value at the permeable wall (which is the non Plug flow region).

For fixed r and [[tau].sub.0], the velocity increases with decrement in [D.sub.a]. As [D.sub.a] decreases the gap between the velocity curves becomes smaller for any fixed [[tau].sub.0]. For larger [[tau].sub.0], there is an increase in plug flow region. From Fig. 13 it is observed that, for the fixed Darcy number [D.sub.a], the velocity decreases along the axis with the increase in [[tau].sub.0]. It is the other way at the permeable wall.

The fractional increase in volume rate of flow is numerically evaluated and is depicted in Fig. 14. It is observed that for a fixed Darcy number [D.sub.a], the fractional increases with yield stress [[tau].sub.0].

6. Conclusion

In this work, we have investigated a Bingham fluid flow through a conical tube with permeable wall. The velocity distribution, the volume rate of flow and its fractional increase are obtained. The Velocity profiles obtained for different values of the yield stress, and a fixed Darcy number, pressure gradient along the flow equal to 10 and slip parameter equal to 1, show that the velocity attains maximum value on the tube axis and decreases with the increment of the rayon.

For a fixed rayon and fixed pressure gradient along the flow, the velocity decreases along the z-axis (when z increase). The slip velocity at the permeable wall decreases with the increment in the value of the Darcy number. The comparison between the velocity profiles for permeable and impermeable conical tubes shows that the velocity is enhanced due to the permeability of the wall of the conical tube. The velocity remains constant in the Plug flow region and decreases in the non Plug flow region.

For fixed radius and fixed yield stress, the velocity increases with decrement in the Darcy number. As Darcy number decreases the gap between the velocity curves becomes smaller for any fixed yield stress. For larger yield stress, there is an increase in plug flow region. For fixed Darcy number, the velocity decreases along the axis flow with the increase in yield stress. It is shown that for a fixed Darcy number, the fractional increases with yield stress.

Nomenclature

[D.sub.a] - Darcy number; L - length of conical tube, mm; Q - Darcy velocity, [ms.sup.-1] ; [R.sub.0] - entrance radius, mm; [R.sub.L] - exit radius, mm; [R.sub.p] - radius of the plug flow region, mm; p - pressure, N/[m.sup.2] ; r - radius, mm; [u.sub.B] - slip velocity, [ms.sup.-1]; [u.sub.av] - average velocity, [ms.sup.-1] ; [v.sub.z] - velocity component in z - direction, [ms.sup.-1]; k - permeability; z - flow axis; [alpha] - slip parameter; [[mu].sub.0] - coefficient of viscosity at the interface, [kgm.sup.-1][s.sup.-1]; [[tau].sub.0] - yield stress, N/[m.sup.2]; [[tau].sub.rz] - shear stress, N/[m.sup.2].

cross ref http://dx.doi.org/10.5755/j01.mech.19.6.5988

Acknowledgements

This work has been supported by Research Laboratory of Industrial Technologies, University Ibn Khaldoun of Tiaret-Algeria.

Received September 19, 2012

Accepted October 10, 2013

References

[1.] Rand, P.W.; Lacombe, E.; Hunt, H.E.; Austin, W.H. 1964. Viscosity of normal blood under Normothermic and Hypothermic conditions, J. Appl. Physiol. 19: 117122.

[2.] Bugliarello, G.; Kapur, C.; Hsiao, G. 1965. Proceedings of the fourth International Congress on Rheology, 4, Symposium of Biorheology (Edited by Copley, A.L.), Inter Sciences New York, 351(3).

[3.] Chien, S.; Usami, S.; Taylor, H.M.; Lundberg, J.L.; Gregre-Sons, M.T. 1965. Effects of haemotocrit and plasma proteins on human blood Rheology at low shear rates, J. Appl. Phys. 21: 81-87.

[4.] Lew, H.S.; Fung, Y.C.; Lowenstein, C.B. 1971. Peristaltic carrying and mixing of chyme in small intestine, J. Biomechanics 4: 297-315. http://dx.doi.org/10.1016/0021-9290(71)90036-4.

[5.] Bird, R.B.; Stewart, W.E.; LightFoot, E.N. 1960. Transport Phenomena, John Wiley & Sons, Inc., New York, 780 p.

[6.] Rathy, R.K. 1976. An Introduction to Fluid Dynamics, Oxford and IBH Publishing Co., New Delhi, Bombay, Calcutta.

[7.] Vajravelu, K.; Sreenadh, S.; Ramakrishna, S.; Arunachalam, P.V. 1987. Bingham fluid flow through a circular pipe with permeable wall, ZAMM 67: 568-569.

[8.] Bervers, G.S.; Joseph, D.D. 1967. Boundary conditions at a naturally permeable wall, J. Fluid Mech. 30: 197-207. http://dx.doi.org/10.1017/S0022112067001375.

[9.] Goverdhan, S.; Sreenadh, S.; Arunachalam, P.V. 1991. Bingham fluid flow between two permeable, Proc. Nat. Acad. Sci. India. 61A: 19-27.

[10.] E. Comparini, E. 1992. A one-dimensional Bingham flow, Journal of Mathematical Analysis and Applications 169: 127-139. http://dx.doi.org/10.1016/0022-247X(92)90107-O.

[11.] Ravana, G.; Ramakrishna, S.; Sreenadh, S.; Arunachalam, P.V. 1996. Bingham fluid flow through an inclined channel over a permeable bed, Bulletin of Pure and Applied Sciences 15E(1): 47-60.

[12.] Comparini, E. 1997. Regularization procedures of singular free boundaries problems in rotational Bingham flows, ZAMM 77(3): 1-12.

[13.] Narahari, M.; Sreenadh, S.; Arunachalam, P.V. 1999. Unsteady flow of a Bingham fluid between two permeable beds, J. Indian Acad. Math. 22(2): 225-242.

[14.] Sankara Reddy, P.; Sreenadh, S.; Arunachalam, P.V. 2000. Bingham fluid flow in an inclined channel bounded by two permeable beds, J. Indian Acad. Math. 22(1): 143-162.

[15.] Viswanatha Reddy, G.; Sreenadh, S.; Arunachalam, P.V. 2001. Bingham fluid flow in an an nulus, J. Indian Acad. Math. 23(2).

[16.] Vajravelu, K.; Ramesh, K.; Sreenadh S.; Arunachalam, P.V. 2003. Pulsatile flow between per meable beds, Int. J. of Non-linear Mechanics 38: 999-1005. http://dx.doi.org/10.1016/S0020-7462(02)00045-8.

N. Sad Chemloul, Department of Mechanical Engineering, University Ibn Khaldoun BP 14000, Tiaret - Algeria, E-mail: sad 1260@yahoo.fr

Printer friendly Cite/link Email Feedback | |

Author: | Chemloul, N. Sad |
---|---|

Publication: | Mechanika |

Article Type: | Report |

Geographic Code: | 4EUFR |

Date: | Nov 1, 2013 |

Words: | 2965 |

Previous Article: | Simulation of turbulent piloted methane non-premixed flame based on combination of finite-rate/eddy-dissipation model/Turbulentines valdomos... |

Next Article: | Research of vibrations in induction machines during transient processes using a ferraris sensor/Indukciniu masinu vibraciju pereinamuju procesu metu... |

Topics: |