Optimal Homotopy Analysis of MHD Natural Convection Flow of Thixotropic Fluid under Subjection of Thermal Stratification: Boundary Layer Analysis ()
1. Introduction
The study of non-Newtonian fluid flow has attracted the attention of researchers over the past years owing to its numerous applications in industry. Some of the applications are encountered in paper production, food processing, glass-fibre and injection moulding, etc. The major difference between a Newtonian and non-Newtonian fluid is that Newtonian fluid, for instance water, does not resist much stress that is applied on them and does not exhibit any sign of strain or deformation unlike non-Newtonian fluid that will quickly exhibit signs of strain when stress is applied. Most of these non-Newtonian fluids are described as “time-dependent” due to the fact that after deformation it takes a fixed time to return to a more viscous state, and likewise the fluid flow is highly dependent on the viscosity. Fluids are referred to as shear thickening when the viscosity increases with the application of the shear stress over time we can categorize rheopetic fluids as shear thickening because they become more viscous over time when shaken or stressed. Fluids are also referred to as shear thinning when the viscosity decreases with the stress over time, thixotropic fluids exhibit this shear thinning property in the sense that, they become less viscous over time when shaken or stressed. In a more simple way, thixotropic fluids are described as fluids that are very thick (viscous) under their normal conditions, but become thin, less viscous and change their flow behavior when pressure is introduced. Many gels and colloids are thixotropic materials. Some natural examples of materials that may be thixotropic are honey, synovial fluid found in joints between some bones, sand, quicksand, clays and drilling muds Hendrickson [1] . In the practical wise, some of implications of the behavior of thixotropic materials are Liquefaction. Liquefaction is referred to as a phenomenon whereby saturated or partially saturated soil substantially loses strength and stiffness in response to an applied stress, usually earthquakes shaking or other sudden changes in stress conditions, causing it to behave like a liquid, that is, it temporarily turns firm ground into a liquid. Although, the impact of liquefaction has being felt some long time ago, but in recent times, a wide destruction made to residential properties in the eastern suburbs and satellite townships of Christchurch, New Zealand during the 2010 Canterbury earthquakes and also the Christchurch earthquakes that occurred in 2011 are caused by Liquefaction Davies et al. [2] . As a result of the complexities encountered working with this fluid, there are virtually very few published works that concentrate on boundary layer flows of the fluids. It was reported in the work of Sadeghy et al. [3] that Harris [4] [5] tried to address the boundary layer flow of thixotropic fluids. It was further reported that, he relied on a simple thixotropic fluid model (the so-called Harris model) to investigate the effects of a fluid’s thixotropic behaviour on the characteristics of the momentum boundary layer formed above a fixed plate and also, Hayat et al. [6] investigated MHD flow of Thixotropic fluid with variable thermal conductivity and thermal radiation, it was reported that the non-Newtonian parameters have quite opposite effects on the velocity and temperature distributions.
Convection on the other hand may be classified depending on how the fluid motion is initiated. Natural convection which is also termed buoyant or free convection is a mechanism of heat transportation in which the fluid motion is not initiated by an external source such as pump, fan, suction device, but only by fluid density differences in the fluid occurring due to temperature gradients and the driving force for natural convection is buoyancy which is as a result of differences in density. During this process of heating, the density changes in the boundary layer will cause the fluid to rise and be replaced by cooler fluid. Definitely, the cooler fluid will also be heated and rise. Recently, Animasaun [7] [8] discussed the convection which occurs in fluid flow along a vertical surface with a small magnitude of temperature (absolute zero) and an upper horizontal surface of a paraboloid of revolution, modified Boussinesq approximation models for both cases were presented.
Stratification of fluid is a formation of layers and occurs due to temperature variations and concentration differences. This natural process creates a transition zone of a temperature gradient between cold and hot fluid zones and hence in case of vertical natural convection it plays an important role in vertical temperature distribution. The dynamics of thermal stratification is essential for solar engineering because stratification may predict the possibilities of achieving higher energy efficiency. Oreyeni et al. [9] extensively discussed homotopy analysis of MHD free convective micropolar fluid flow along a vertical surface embedded in non-Darcian thermally-stratified medium. It was reported that stratification leads to decline in the velocity and temperature of the fluid respectively and the result is similar to what is observed in the work of Animasaun and Omowaye [10] on upper-convected Maxwell fluid flow with variable thermo-physical properties over a melting surface situated in hot environment subject to thermal stratification and the work of Hayat et al. [11] on thermal and concentration stratifications effects in radiative flow of a Jeffery fluid over a stretching sheet.
Internal energy generation is described as scientific method used in generating heat energy within a body by a chemical, electrical or nuclear process. Eldahab and Aziz [12] investigated the effect of non-uniform heat source with suction-blowing. Muthtamilselvan et al. [13] analyzed effect of non-uniform heat generation on unsteady MHD non-Darcian flow over a vertical stretching surface with variable properties. Crepeau and Clarksean [14] examined a similarity solution for a fluid with an exponentially decaying heat generation term and a constant temperature vertical plate under the assumption that the fluid has an internal volumetric heat generation. An exponentially form is used to account for the internal energy generation term. It was reported that the effect of internal heat generation is important in several applications i.e. reactor safety analysis, fire and combustion studies. In many situations, there may be appreciable temperature difference between the surface and the ambient fluid.
Physical changes of thixotropy which referred to the reversible changes from a flowable fluid to a solid-like elastic gel has been made to occur by changing the temperature. But due to the fact that thixotropic fluid exhibit time-dependent shear thinning property, viscosity of the fluid is assumed to be constant and thermal conductivity is assumed not to be constant and has been taken as function of temperature. A slight increase in the temperature enhances the transport phenomena by reducing the viscosity across the momentum boundary layer, hence the heat transfer rate at the wall is also affected. Sharma and Singh [15] presented effects of variable thermal conductivity and heat source/sink on mhd flow near a stagnation point on a linearly stretching sheet. Effects of variable thermal conductivity on heat and mass transfer with Jeffery fluid was investigated by Uwanta and Omokhuale [16] . Prasad et al. [17] discussed the effects of variable fluid properties on the hydromagnetic flow and heat transfer over a non-linearly stretching sheet. Omokhuale et al. [18] studied effect of Jeffery fluid on heat and mass transfer past a vertical porous plate with soret and variable thermal conductivity. Motivated by all the above works and applications, it is of interest to contribute to the body of knowledge by investigating the natural convection flow of thixotropic fluid in a thermally stratified medium and considering cases where the non-Newtonian parameters are properly accounted for. Furthermore, this study examines the characteristics behavior of thixotropic fluid on boundary layer.
2. Mathematical Formulation
A steady and incompressible magnetohydrodynamic (MHD) free convective boundary layer flow of thixotropic fluid over a vertical surface in a thermally stratified medium with temperature-dependent thermal conductivity is considered. The motion of an incompressible non-Newtonian fluid is induced by the stretching property of the vertical surface, buoyancy effect which is generated by gradients in the temperature field and space dependent internal heat generation. The flow is assumed to flow in x-direction which is along vertical surface and y-axis is normal to it. The temperature of the surface
is held uniform and embedded in a thermally stratified medium of variable ambient temperature
where
. A constant magnetic field of strength
is applied in the y-direction. Also the magnetic Reynolds number is assumed to be small so that the induced magnetic field is negligible in comparison to the applied magnetic field.
The equation of continuity can be written as;
(1)
For any incompressible fluid satisfying the thixotropic model with constant viscosity and temperature-dependent thermal conductivity, the momentum and energy equations can be simplified using the usual boundary layer theory approximations and obtain
(2)
(3)
Subject to boundary conditions
(4)
(5)
In order to justify the variation in the thermo-physical property of the thixotropic fluid as it flows past a vertical heated surface, classical Boussinesq's approximation is adopted such that the temperature at the surface is greater than temperature of the fluid at the free stream. It is valid to consider the mathematical model of temperature dependent thermal conductivity model of Charraudeau [19] as;
(6)
are the constant value of the coefficient thermal conductivity at the free stream. Where u and v are components of velocity in x and y directions respectively,
is the wall shrinking or stretching velocity, (
) for stretching, (
) for shrinking and (
) for static wall,
is the wall mass flux velocity,
is the fluid density,
is the kinematic viscosity,
is the electrical conductivity, T is the fluid temperature in the boundary layer,
is the free stream temperature,
is the thermal expansion coefficient,
is the thermal diffusivity,
and
are the non-Newtonian material constants,
is the volumetric coefficient of thermal expansion, A and B are the coefficient of exponentially decaying space and temperature dependent heat source respectively.
The momentum and energy equations can be transformed into the corresponding ordinary differential equations by the following transformations
(7)
where
is the independent dimensionless similarity variable. Thus u and v are given by
,
, substituting variables (7) into equations (2)-(5), we obtain the following ordinary differential equations:
(8)
(9)
The corresponding boundary conditions take the form;
(10)
(11)
In the above equations, primes denote differentiation with respect to
. The dimensionless velocity and temperature are represented as
and
,
is the Prandtl number,
is the magnetic parameter,
is the Modified Thermal Grashof number, A is the space dependent heat source parameter,
is the thermal stratification parameter,
is the temperature-dependent thermal conductivity parameter, B is the temperature dependent heat source parameter,
and
are the non-Newtonian parameters.
The dimensionless form of the skin friction coefficient is given as;
(12)
3. Optimal Homotopy Analysis Solutions
The analytic solutions of Equations (8) and (9) to the boundary conditions (10) and (11) have been computed by optimal homotopy analytic method (OHAM). Various researchers ( [20] [21] [22] ) have applied this method to compute flow problems. The suitable initial guesses and the linear operators satisfying the given boundary conditions (10) and (11) for
and
have been carefully selected. The initial guesses for this problem are of the form
(13)
Linear operators
and
are
(14)
(15)
The operators
and
have the following properties
(16)
In which
,
,
,
and
are constants.
(17)
(18)
Subject to boundary conditions
(19)
(20)
where q is the embedding parameters,
and
the non-zero auxilliary parameters and operators
and
are given by
(21)
(22)
Obviously, when
and
, zero order of deformation equations (17)-(18) leads to
(23)
(24)
Expanding
and
in Taylor series with respect to the embedding parameter q,
(25)
(26)
The auxilliary parameters
and
are properly chosen so that the series Equations (17)-(18) converge at
. Hence,
(27)
(28)
If we denote the special solutions
and
then the general solutions
and
are
(29)
(30)
Following the rule of solution expression, the rule of coefficient ergodicity and the rule of solution existence as discussed in [23] and [24] we choose auxiliary functions as
(31)
In order to analyze the salient characteristics of the problem, the results are presented in Figure 1(a)-Figure 8(b). The effect of non-Newtonian parameters
on velocity and temperature profiles are revealed in Figure 1(a), Figure 1(b) and Figure 2(a), Figure 2(b) for different cases of stratification
Figure 1. (a) Effect of non-Newtonian parameters
on velocity profile when
(b) Effect of non-Newtonian parameters
on velocity profile when
.
Figure 2. (a) Effect of non-Newtonian parameters
on temperature profile when
(b) Effect of non-Newtonian parameters
temperature profile when
.
, it is observed from Figure 1(a) for the case of very low stratification
, that an increase in the values of
from 0 through 0.3, 0.7 to 1.0 corresponds to a slight increase in the velocity distribution. Under normal circumstances, it is noticed that by increasing the values of
and
velocity and the boundary layer thickness increase. Also, from Figure 1(b) for a high value of stratification
, it is noticed that as the values of
increases, the velocity profile increases, but at the exact value of
all profiles merge and asymptotically tends to 0. From Figure 2(a) and Figure 2(b) it is observed that increase in
leads to decrease in the temperature distribution for both cases of
and
. It is noticed for the first case of stratification
that as
increases in magnitude, there is a notable decrease in the temperature profile for value of
which signifies the Newtonian fluid, and thereafter the profiles decreases slightly as
increases. Also, for the cases of
in Figure 2(b) it is observed that the temperature profile substantially decreases as the values of
increases. From Figure 3(a), it is seen that an increase in
when
corresponds to increase in the wall shear stress. Also, from Figure 3(b) it is noticed that as the
increases the wall shear stress decreases for the case of the non-Newtonian fluid i.e.
. The effect of thermal stratification parameter
on temperature distributions are presented in Figure 4(a) and Figure 4(b) for both cases of
(Newtonian) and
(non-Newtonian). It is depicted from the computations that for both cases of
the temperature profiles decrease substantially as the magnitude of
is increased from 0.1 through 0.4, 0.7 to 1.0 making each profile to satisfy the boundary condition associated with the model that is, (
). Apparently, when
reflecting in the boundary condition (
) which results in the decrease in the value of temperature
Figure 3. (a) Effect of non-Newtonian parameters
on the wall shear stress when
(b) Effect of Stratification parameter
on the wall shear stress when
.
Figure 4. (a) Effect of Stratification parameter
on temperature profile when
(b) Effect of Stratification parameter
on temperature profile when
.
as revealed in the temperature profiles in Figure 4(a) and Figure 4(b). In other words, this can be traced to the fact that as
increases, the surface temperature within the thermally stratified medium decreases. Physically, this decrease in the temperature distribution as the thixotropic fluid flows along a vertical surface may account for an increase in the viscosity of the fluid, and as a result of decrease in temperature the intermolecular forces binding the microstructure of thixotropic fluid becomes stronger and makes it possible for the fluid to be more viscous. Figure 5(a), Figure 5(b) and Figure 6(a), Figure 6(b) illustrate the influence of Magnetic parameter M on velocity and temperature distributions respectively for different cases of non-Newtonian parameters
. It is revealed that, as M increases at a fixed value of
, the velocity profiles decrease and the temperature profiles increase. Physically, this observation is due to the fact that, application of transverse magnetic field normal to the direction of the flow brings about a resistive force which has tendency to reduce the velocity of the fluid and hence makes the temperature of the fluid to rise. In order to further unravel the behaviour of the fluid, the effect of temperature-dependent thermal conductivity parameter
on temperature profiles are investigated in Figure 7(a) and Figure 7(b) for the different cases of
when
. It is observed from Figure 7(a) that increase in the values of
corresponds to increase in the temperature distribution for the case of
(Newtonian) and all profiles quickly merge at
and descend towards the free stream. Also, from Figure 7(b) it is observed that as
increases from 1 through 2.0, 3.0 to 4.0 for the case of
there is an obvious increase in temperature profile within the domain
and within the domain
all profiles merge together, thereafter at specific value of
the temperature decreases away from the surface and all profiles decay smoothly
Figure 5. (a) Effect of Magnetic parameter M on velocity profile when
(b) Effect of Magnetic parameter M on velocity profile when
.
Figure 6. (a) Effect of Magnetic parameter M on temperature profile when
(b) Effect of Magnetic parameter M on temperature profile when
.
to zero as
. Figure 8(a) and Figure 8(b) illustrates the influence of Grashof number
on the velocity and temperature profiles when
and
. From Figure 8(a) as
increases from 0.1 the velocity profile increases but as
increases through
all profiles conflate within the domain
and the velocity increases towards the free stream. Likewise, from Figure 8(b) it is observed that increase in
leads to considerable decrease in the temperature profile.
4. Convergence of the Homotopy Solutions
It is obvious that the series (27) and (28) consist of the non-zero auxiliary parameters
and
which can adjust and control the convergence. The
Figure 7. (a) Effect of Temperature-dependent thermal conductivity parameter
on temperature profile when
(b) Effect of Temperature-dependent thermal conductivity parameter
on temperature profile when
.
Figure 8. (a) Effect of Grashof number
on velocity profile when
(b) Effect of Grashof number
on temperature profile when
.
interval on
-axis for which the
-curve becomes parallel to the
-axis is recognized as the set of admissible values of
and
for which the solution series converges. These figures show that the ranges for the acceptable values of
and
are
and
. Figure 9(a) and Figure 9(b) depict the
-curves for this problem, we obtained the approximate optimal values of
and
at 10th-order of approximation as −0.797715 and −0.989828.
5. Results and Discussion
Computation has been carried out for various values of emerging parameters like; stratification parameter, magnetic parameter, temperature dependent
Figure 9. (a)
-curve of
obtained at 10th-order of approximation (b)
-curve of
obtained at 10th-order of approximation.
Table 1. Comparison of the results for the skin-friction coefficient
at
,
,
,
for various values of s.
Table 2. Numerical values of Skin-friction coefficient for various values s,
and
when
and
.
thermal conductivity parameter, non-Newtonian parameters and heat source parameters. Table 1 is being validated by direct comparison with the numerical results reported by Hayat et al. [6] when
. The numerical computations of skin-friction coefficient
for the values of s, M and
and
are presented in Table 2, it is shown that the values of skin-friction coefficient
increase by increasing s and decrease by increasing
and
.
6. Conclusions
Free convection MHD Thixotropic fluid flow past a vertical surface with temperature-dependent thermal conductivity, subjected to thermal stratification is investigated. Series solutions for the velocity and temperature fields were gotten and discussed. Effects of various parameters governing the flow fields were investigated. From this study, we conclude thus;
1) Increase in the values of
corresponds to increase in the velocity profile for both cases of
and
with opposite effect of the temperature profiles.
2) The wall shear stress increases as the values of
increase when
.
3) Temperature profile is a decreasing function of
for both cases of
and
.
4) The wall shear stress decreases with increase in values of
when
.