SciELO - Scientific Electronic Library Online

 
vol.28 número3Critical flow regions in tissue artificial heart valve assessed by laser doppler anemometer in continuous flowSome observations on wear and damages in cemented carbide tools índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

Compartilhar


Journal of the Brazilian Society of Mechanical Sciences and Engineering

versão impressa ISSN 1678-5878versão On-line ISSN 1806-3691

J. Braz. Soc. Mech. Sci. & Eng. v.28 n.3 Rio de Janeiro jul./set. 2006

http://dx.doi.org/10.1590/S1678-58782006000300003 

TECHNICAL PAPERS

 

Unsteady MHD flow of a dusty non-Newtonian Bingham fluid through a circular pipe

 

 

Hazem A. Attia

Department of Mathematics, College of Science; Al-Qasseem University; P.O. Box 237, Buraidah 81999; Kingdom of Saudi Arabia

 

 


ABSTRACT

In this paper, the transient magnetohydrodynamic (MHD) flow of a dusty incompressible electrically conducting non-Newtonian Bingham fluid through a circular pipe is studied taking the Hall effect into consideration. A constant pressure gradient in the axial direction and an uniform magnetic field directed perpendicular to the flow direction are applied. The particle-phase is assumed to behave as a viscous fluid. A numerical solution is obtained for the governing nonlinear equations using finite differences.

Keywords: Fluid mechanics, magneto-fluid mechanics, circular pipe flow, non-Newtonian fluid, Bingham fluid


 

 

Introduction

The flow of a dusty and electrically conducting fluid through a circular pipe in the presence of a transverse magnetic field has important applications such as MHD generators, pumps, accelerators, and flowmeters. The performance and efficiency of these devices are influenced by the presence of suspended solid particles in the form of ash or soot as a result of the corrosion and wear activities and/or the combustion processes in MHD generators and plasma MHD accelerators. When the particle concentration becomes high, mutual particle interaction leads to higher particle-phase viscous stresses and can be accounted for by endowing the particle phase by the so-called particle-phase viscosity. There have been many articles dealing with theoretical modelling and experimental measurements of the particle-phase viscosity in a dusty fluid (Soo 1969, Gidaspow et al. 1986, Grace 1982, and Sinclair et al. 1989).

The flow of a conducting fluid in a circular pipe has been investigated by many authors (Gadiraju et al. 1992, Dube et al. 1975, Ritter et al. 1977, and Chamkha 1994). Gadiraju et al. (1992) investigated steady two-phase vertical flow in a pipe. Dube et al. (1975) and Ritter et al. (1977) reported solutions for unsteady dusty-gas flow in a circular pipe in the absence of a magnetic field and particle-phase viscous stresses. Chamkha (1994) obtained exact solutions which generalize the results reported in Dube et al. 1975 and Ritter et al. 1977 by the inclusion of the magnetic and particle-phase viscous effects. It should be noted that in the above studies the Hall effect is ignored.

A number of industrially important fluids such as molten plastics, polymers, pulps and foods exhibit non-Newtonian fluid behavior (Nakayama et al. 1988). Due to the growing use of these non-Newtonian materials, in various manufacturing and processing industries, considerable efforts have been directed towards understanding their flow characteristics. Many of the inelastic non-Newtonian fluids, encountered in chemical engineering processes, are known to follow the so-called "power-law model" in which the shear stress varies according to a power function of the strain rate (Metzner et al. 1965). It is of interest in this paper to study the influence of the magnetic field as well as the non-Newtonian fluid characteristics on the dusty fluid flow properties in situations where the particle-phase is considered dense enough to include the particulate viscous stresses.1

In the present study, a new element is added to the problem studied by Attia (2003) by taking the Hall effect into consideration. Therefore, the unsteady flow of a dusty non-Newtonian Bingham fluid through a circular pipe is investigated considering the Hall effect. The carrier fluid is assumed viscous, incompressible and electrically conducting. The particle phase is assumed to be incompressible pressureless and electrically non-conducting. The flow in the pipe starts from rest through the application of a constant axial pressure gradient. The governing nonlinear momentum equations for both the fluid and particle-phases are solved numerically using the finite difference approximations. The effect of the Hall current, the non-Newtonian fluid characteristics and the particle-phase viscosity on the velocity of the fluid and particle-phases are reported.

 

Governing Equations

Consider the unsteady, laminar, and axisymmetric horizontal flow of a dusty conducting non-Newtonian Bingham fluid through an infinitely long pipe of radius "d" driven by a constant pressure gradient. A uniform magnetic field is applied perpendicular to the flow direction. The Hall current is taken into consideration and the magnetic Reynolds number is assumed to be very small, consequently the induced magnetic field is neglected (Sutton et al. 1965). We assume that both phases behave as viscous fluids and that the volume fraction of suspended particles is finite and constant (Chamkha 1994). Taking into account these and the previously mentioned assumptions, the governing momentum equations can be written as

where t is the time, r is the distance in the radial direction, V is the fluid-phase velocity, Vp is the particle-phase velocity, r is the fluid-phase density, rp is the particle-phase density, P/z is the fluid pressure gradient, j is the particle-phase volume fraction, N is a momentum transfer coefficient (the reciprocal of the relaxation time, the time needed for the relative velocity between the phases to reduce e-1 of its original value (Chamkha 1994), s is the fluid electrical conductivity, m=sgBo is the Hall parameter, g is the Hall factor (Sutton et al. 1965), Bo is the magnetic induction, µp is the particle-phase viscosity which is assumed constant, and µ is the apparent viscosity of the fluid which is given by,

where µo is the plastic viscosity of a Bingham fluid and to is the yield stress. In this work, r, rp, µp, j, and Bo are all constant. It should be pointed out that the particle-phase pressure is assumed negligible and that the particles are being dragged along with the fluid-phase.

The initial and boundary conditions of the problem are given as

where "d" is the pipe radius.

Equations (1)-(3) constitute a nonlinear initial-value problem which can be made dimensionless by introducing the following dimensionless variables and parameters

By introducing the above dimensionless variables and parameters as well as the expression of the fluid viscosity defined above, Eqs. (1)-(3) can be written as (the bars are dropped),

The volumetric flow rates and skin-friction coefficients for both the fluid and particle phases are defined, respectively, as (Chamkha 1994)

 

Results and Discussion

Equations (6) and (7) represent coupled system of non-linear partial differential equations which are solved numerically under the initial and boundary conditions (8) using the finite difference approximations. A linearization technique is first applied to replace the nonlinear terms at a linear stage, with the corrections incorporated in subsequent iterative steps until convergence is reached. The computational domain is divided into meshes each of dimension Dt and Dr in time and space, respectively. Then the Crank-Nicolson implicit method is used at two successive time levels (Mitchell et al., 1980; Evans et al., 2000). An iterative scheme is used to solve the linearized system of difference equations. The solution at a certain time step is chosen as an initial guess for next time step and the iterations are continued till convergence, within a prescribed accuracy. Finally, the resulting block tri-diagonal system is solved using the generalized Thomas-algorithm (Mitchell et al., 1980; Evans et al., 2000). Computations have been made for a=1 and k=10. Grid-independence studies show that the computational domain 0<t<¥ and 0<r<1 can be divided into intervals with step sizes Dt=0.0001 and Dr=0.005 for time and space respectively. Smaller step sizes do not show any significant change in the results. Convergence of the scheme is assumed when all of the unknowns V, Vp, V/r, and Vp/r for the last two approximations differ from unity by less than 10-6 for all values of r in 0<r<1 at every time step. It should be mentioned that the results obtained herein reduce to those reported by Dube et al. (1975) and Chamkha (1994) for the cases of non-magnetic, inviscid particle-phase (B=0), and Newtonian fluid. These comparisons lend confidence in the accuracy and correctness of the solutions.

Imposing of a magnetic field normal to the flow direction gives rise to a drag-like or resistive force and it has the tendency to slow down or suppress the movement of the fluid in the pipe, which in turn, reduces the motion of the suspended particle-phase. This is translated into reductions in the average velocities of both the fluid- and the particle-phases and, consequently, in their flow rates. In addition, the reduced motion of the particulate suspension in the pipe as a result of increasing the strength of the magnetic field causes lower velocity gradients at the wall. This has the direct effect of reducing the skin-friction coefficients of both phases. Including the Hall parameter decreases the resistive force imposed by the magnetic field due to its effect in reducing the effective conductivity. Therefore, the Hall parameter leads to an increase in the average velocities of both the fluid- and the particle-phases and, consequently, in their flow rates and the velocity gradients at the wall.

Figures 1 and 2 present the time evolution of the profiles of the velocity of the fluid V and dust particles Vp, respectively, for various values of the Bingham number tD and for m=0, Ha=0.5 and B=0.5. Both V and Vp increase with time and V reaches the steady-state faster than Vp for all values of tD. It is clear from Figs. 1 and 2 that increasing tD, which increases the driving force for V, increases V and, consequently, increases Vp while its effect on their steady-state times can be neglected.

 

 

 

 

Figures 3 and 4 present the time evolution of the profiles of the velocity of the fluid V and dust particles Vp, respectively, for various values of the Bingham number tD and for m=1, Ha=0.5 and B=0.5. It is indicated in the figures that increasing m increases V and, in turn, Vp due to the decrease in the effective conductivity (s/(1+m2)) which reduces the damping magnetic force on V. It is shown that the influence of the Hall parameter m on V is more apparent for higher values of tD.

 

 

 

 

Table 1 presents the steady state values of the fluid-phase volumetric flow rate Q, the particle-phase volumetric flow rate Qp, the fluid-phase skin friction coefficient C, and the particle-phase skin friction coefficient Cp for various values of the parameters tD and m and for Ha=0.5 and B=0.5. It is clear that increasing the parameter m increases Q, Qp, C, and Cp for all values of tD. This comes from the fact that increasing m increases the velocities and their gradients which increases the average velocities of both the fluid- and the particle-phases and, consequently, increases their flow rates and skin-friction coefficients of both phases. It is also shown that increasing tD increases Q, Qp, C, and Cp for all values of m as a result of increasing the velocities of both phases.

 

 

Table 2 presents the steady state values of the fluid-phase volumetric flow rate Q, the particle-phase volumetric flow rate Qp, the fluid-phase skin friction coefficient C, and the particle-phase skin friction coefficient Cp for various values of the parameters m and

B and for Ha=0.5 and tD=0. It is clear that, increasing m increases Q, Qp, C, and Cp for all values of B and its effect becomes more pronounced for smaller values of B. Increasing the parameter B decreases the quantities Q, Qp, and C, but increases Cp for all values of m. This can be attributed to the fact that increasing B increases viscosity and therefore the flow rates of both phases as well as the fluid-phase wall friction decreases considerably. However, since Cp in defined as directly proportional to B, it increases as B increases at all times.

 

 

Conclusion

The unsteady MHD flow of a particulate suspension in an electrically conducting non-Newtonian Bingham fluid in a circular pipe is studied considering the Hall effect. The governing nonlinear partial differential equations are solved numerically using finite differences. The effect of the magnetic field parameter Ha, the Hall parameter, the non-Newtonian fluid characteristics (Bingham number tD), and the particle-phase viscosity b on the transient behavior of the velocity, volumetric flow rates, and skin friction coefficients of both fluid and particle-phases is studied. It is shown that increasing the magnetic field decreases the fluid and particle velocities, while increasing the Hall parameter increases both velocities. It is found that increasing the parameter m increases Q, Qp, C, and Cp for all values of tD. The effect of the Hall parameter on the quantities Q, Qp, C, and Cp becomes more pronounced for smaller values of b.

 

References

Attia, H.A., 2003, "Unsteady flow of a dusty conducting non-Newtonian fluid through a pipe", Can. J. Phys., Vol. 81(3), pp. 789-795.        [ Links ]

Chamkha, A.J., 1994, "Unsteady flow of a dusty conducting fluid through a pipe", Mechanics Research Communications, Vol. 21(3), pp. 281-286.        [ Links ]

Dube, S.N. and Sharma, C.L., 1975, "A note on unsteady flow of a dusty viscous liquid in a circular pipe", J. Phys. Soc. Japan, Vol. 38(1), pp. 298-310.        [ Links ]

Evans, G.A., Blackledge, J.M. and Yardley, P.D., 2000, "Numerical methods for partial differential equations", Springer Verlag, New York.        [ Links ]

Gadiraju, M., Peddieson, J. and Munukutla, S., 1992, "Exact solutions for two-phase vertical pipe flow", Mechanics Research Communications, Vol. 19(1), pp. 7-13.        [ Links ]

Gidaspow, D., 1986, "Hydrodynamics of fluidization and heat transfer: super computer modeling", Appl. Mech. Rev., Vol. 39, pp. 1-23.        [ Links ]

Grace, J.R., 1982, "Fluidized-Bed Hydrodynamic, Handbook of Multiphase Systems", McGraw-Hill, New York.        [ Links ]

Metzner, A.B., 1965, "Heat Transfer in non-Newtonian fluid", Adv. Heat Transfer, Vol. 2, pp. 357-397.        [ Links ]

Mitchell, A.R. and Griffiths, D.F., 1980, "The finite difference method in partial differential equations", John Wiley & Sons, New York.        [ Links ]

Nakayama, A. and Koyama, H., 1988, "An analysis for friction and heat transfer characteristics of power-law non-Newtonian fluid flows past bodies of arbitrary geometrical configuration", Warme-und Stoffubertragung, Vol. 22, pp. 29-37.        [ Links ]

Ritter, J.M. and Peddieson, J., 1977, "Transient two-phase flows in channels and circular pipes", Proceedings of the 6th Canadian Congress of Applied Mechanics, 1977.        [ Links ]

Sinclair, J.L. and Jackson, R., 1989, "Gas-particle flow in a vertical pipe with particle-particle interactions", AICHE J., Vol. 35, pp. 1473-1486.        [ Links ]

Soo, S.L., 1969, "Pipe flow of suspensions", Appl. Sci. Res., Vol. 21, pp. 68-84.        [ Links ]

Sutton, G.W. and Sherman, A., 1965, "Engineering Magnetohydrodynamics", McGraw-Hill, New York.        [ Links ]

 

 

Paper accepted May 2006.

 

 

Technical Editor: Aristeu da Silveira Neto.

Creative Commons License Todo o conteúdo deste periódico, exceto onde está identificado, está licenciado sob uma Licença Creative Commons