# Mixed Convection Boundary Layer Flow from a Solid Sphere with Newtonian Heating in a Micropolar Fluid

Authors: M. Z. Salleh; R. Nazar; I. Pop
Journal: SRX Physics
 (2010)

---

## Abstract

The steady mixed convection boundary layer flow from a solid sphere in a micropolar fluid, generated by Newtonian heating in which the heat transfer from the surface is proportional to the local surface temperature, is considered. The governing boundary layer equations are first transformed into a system of nondimensional equations via the non-dimensional variables, and then into nonsimilar equations before they are solved numerically using an implicit finite-difference scheme known as the Keller-box method. Numerical solutions are obtained for the skin friction coefficient, wall temperature and heat transfer coefficient, as well as the velocity and temperature profiles with several parameters considered, namely the mixed convection parameter, the material or micropolar parameter, and the Prandtl number.


---

## Body

## 1. Introduction

The essence of the theory of micropolar fluid flow lies in the extension of the constitutive equation for Newtonian fluid, so that more complex fluids such as particle suspensions, liquid crystal, animal blood, lubrication, and turbulent shear flows can be described by this theory. The theory of micropolar fluid was first proposed by Eringen [1]. Extensive review of the theory and applications can be found in the review article by Ariman et al. [2] and quite recent papers by Rees and Bassom [3], Pop et al. [4], Nazar et al. [5–8], and so forth. This clearly shows the fast development of the theory of micropolar fluid. In summary, all the papers above considered the boundary condition of either a constant wall temperature or constant heat flux.Generally, in modeling the boundary layer flow and heat transfer of these problems, the boundary conditions that are usually applied are either a constant wall temperature (CWT) or a constant wall heat flux (CHF). However, there is a class of boundary layer flow and heat transfer problems in which the surface heat transfer depends on the surface temperature. Perhaps the simplest case of which is when there is a linear relation between the surface heat transfer and surface temperature. This situation arises in conjugate heat transfer problems (see, e.g., Merkin and Pop [9]), and also when there is Newtonian heating of the convective fluid from the surface. The latter case has been first discussed in detail by Merkin [10]. The situation with Newtonian heating arises in what are usually termed conjugate convective flows, where the heat is supplied to the convecting fluid through a bounding surface with a finite heat capacity. This results in the heat transfer rate through the surface being proportional to the local difference in the temperature with the ambient conditions. However, the Newtonian heating conditions have been used only recently by Lesnic et al. [11–13] and Pop et al. [14] to study the free convection boundary layer over vertical and horizontal surfaces as well as over a small inclined flat plate from the horizontal surface embedded in a porous medium. The asymptotic solution near the leading edge and the full numerical solution along the whole plate domain have been obtained numerically, whilst the asymptotic solution far downstream along the plate has been obtained analytically. Chaudhary and Jain [15, 16] studied the unsteady free convection boundary layer flow past an impulsively started vertical infinite flat plate with Newtonian heating.This configuration occurs in many important engineering devices, for example, in heat exchanger where the conduction in solid tube wall is greatly influenced by the convection in the fluid flowing over it. Furthermore, for conjugate heat transfer around fins where the conduction within the fin and the convection in the fluid surrounding it must be simultaneously analyzed in order to obtain the vital design information and also in convection flow setup when the bounding surfaces absorb heat by solar radiation [15, 16]. This results in the heat transfer rate through the surface being proportional to the local difference in the temperature with the ambient conditions.Recently, Salleh et al. [17, 18] employed an implicit finite-difference scheme, namely the Keller-box method to obtain numerical solutions for mixed convection boundary layer flow near the lower stagnation point of a solid sphere and forced convection boundary layer flow at a forward stagnation point with Newtonian heating in viscous fluids, respectively.Therefore, the aim of the present paper is to study the problem of mixed convection boundary layer flow from a solid sphere with Newtonian heating in a micropolar fluid. The governing boundary layer equations are first transformed into a system of non-dimensional equations via the non-dimensional variables, and then into non-similar equations before they are solved numerically by the Keller-box method, as described in the books by Na [19] and Cebeci and Bradshaw [20].  To the best of our knowledge, this present problem (for the case of Newtonian heating) has not been considered before, so that the reported results are new.


## 2. Analysis

Consider a heated sphere of radiusa, which is immersed in a viscous and incompressible micropolar fluid of ambient temperature T∞, which is subjected to a Newtonian heating (NH) as shown in Figure 1. The convective forced flow is assumed to be moving upward, while the gravity vector g acts downward in the opposite direction, where the coordinates x̅ and y̅ are chosen such that x̅ measures the distance along the surface of the sphere from the lower stagnation point and y̅ measures the distance normal to the surface of the sphere. We assume that the equations are subjected to a Newtonian heating of the form proposed by Merkin [10]. Under the Boussinesq and boundary layer approximations, the basic equations are (Nazar et al. [6, 8])Figure 1
Physical model and coordinate system.(1)∂∂x̅(r̅u̅)+∂∂y̅(r̅v̅)=0,u̅∂u̅∂x̅+v̅∂u̅∂y̅=u̅edu̅edx̅+(μ+κρ)∂2u̅∂y̅2±gβ(T-T∞)sin(x̅a)+κρ∂H̅∂y̅,ρj(u̅∂H̅∂x̅+v̅∂H̅∂y̅)=-κ(2H̅+∂u̅∂y̅)+γ∂2H̅∂y̅2,u̅∂T∂x̅+v̅∂T∂y̅=νPr∂2T∂y̅2
subject to the boundary conditions(2)u̅=v̅=0,∂T∂y̅=-hsT(NH)orT=Tw(CWT)or∂T∂y̅=-qwk(CHF),H̅=-n∂u̅∂y̅ony̅=0u̅→u̅e(x̅),T→T∞,H̅→0asy̅→∞,
where u̅ and v̅ are the velocity components along the x̅ and y̅ directions, respectively, H̅ is the microrotation component normal to the x̅-y̅ plane, γ is the spin gradient viscosity, κ is the vortex viscosity, T is the local temperature, Tw is the wall temperature, g is the gravity acceleration, α=ν/Pr is the thermal diffusivity, β is the thermal expansion coefficient, ν=μ/ρ is the kinematic viscosity, μ is the dynamic viscosity, ρ is the density, Pr is the Prandtl number, k is the thermal conductivity, qw is the constant heat flux from the wall, and hs is the heat transfer parameter. It is worth mentioning that in boundary conditions (2), n is a constant and 0≤n≤1. The value n=0, which indicates H̅=0 at the wall, represents concentrated particle flows in which the particle density is sufficiently great that microelements close to the wall are unable to rotate or called as “strong” concentration of microelements [21, 22]. The case corresponding to n=1/2 results in the vanishing of antisymmetric part of the stress tensor and represents “weak” concentration of microelements [23]. In this case, the particle rotation is equal to fluid vorticity at the boundary for fine particle suspension. When n=1, we have flows which are representative of turbulent boundary layer [24]. The case of  n=1/2 is considered in the present study.Letr̅(x̅) be the radial distance from the symmetrical axis to the surface of the sphere and u̅e the local free stream velocity, which are given by(3)r̅(x̅)=asin(x̅a),u̅e(x̅)=32U∞sin(x̅a),γ=(μ+κ2)j,
where j is the micro-inertia per unit mass.We introduce now the following non-dimensional variables:(4)x=x̅a,y=Re1/2(y̅a),u=u̅U∞,v=Re1/2(v̅U∞),H=(aU∞)Re-1/2H̅,ue(x)=u̅e(x̅)U∞,hs=Re1/2a,r(x)=r̅(x̅)a,θ=T-T∞T∞(NH)orθ=T-T∞Tw-T∞(CWT)orθ=(kaqw)Re1/2(T-T∞)(CHF),
where Re=U∞a/ν is the Reynolds number. Using (4), the system of (1) then becomes(5)∂∂x(ru)+∂∂y(rv)=0,(6)u∂u∂x+v∂u∂y=ueduedx+(1+K)∂2u∂y2+K∂H∂y+λθsinx,(7)u∂H∂x+v∂H∂y=-K(2H+∂u∂y)+(1+K2)∂2H∂y2,(8)u∂θ∂x+v∂θ∂y=1Pr∂2θ∂y2
subject to  the boundary conditions(9)u=v=0,∂θ∂y=-(1+θ)(NH),orθ=1(CWT)or∂θ∂y=-1(CHF),H=-12∂u∂yony=0,u→ue(x),θ→0,H→0asy→∞,
where λ is the mixed convection parameter and K is the material parameter which are given by(10)λ=GrRe2(NHandCWT)orλ=GrRe5/2(CHF),K=κμ
with Gr=gβT∞a3/ν2(NH) or Gr=gβ(Tw-T∞)a3/ν2 (CWT) or Gr=gβ(aqw/k)a3/ν2 (CHF)  being the Grashof number. It is worth mentioning that in both cases of CWT and CHF, λ>0 is for the aiding or assisting flow (heated sphere) and λ<0  is for the opposing flow (cooled sphere), while for the present case of NH, the value of λ considered is only for λ>0.Further, the stream functionψ is introduced and it is defined as(11)u=1r∂ψ∂y,v=-1r∂ψ∂x.
Thus, since sinx/x→1 as x→0, it is appropriate to introduce the transformation(12)ψ=xr(x)f(x,y),θ=θ(x,y),H=xh(x,y).Equations (6) to (8) then become(13)(1+K)∂3f∂y3+(1+xsinxcosx)f∂2f∂y2-(∂f∂y)2+94sinxcosxx+K∂h∂y+λsinxxθ=x(∂f∂y∂2f∂x∂y-∂f∂x∂2f∂y2),(1+K2)∂2h∂y2+(1+xsinxcosx)f∂f∂y-∂f∂yh-K(2h+∂2f∂y2)=x(∂f∂y∂h∂x-∂f∂x∂h∂y),1Pr∂2θ∂y2+(1+xsinxcosx)f∂θ∂y=x(∂f∂y∂θ∂x-∂f∂x∂θ∂y)
subject to the boundary conditions(14)f=∂f∂y=0,∂θ∂y=-(1+θ)(NH)orθ=1(CWT)or∂θ∂y=-1(CHF),h=-12∂2f∂y2ony=0∂f∂y→32sinxx,θ→0,h→0asy→∞.At the lower stagnation point of the sphere,x≈0, (13) reduces to the following ordinary differential equations:(15)(1+K)f′′′+2ff′′-f′2+Kh′+λθ+94=0(1+K2)h′′+2fh′-f′h-K(2h+f′′)=01Prθ′′+2fθ′=0
and the boundary conditions (14) become(16)f(0)=f'(0)=0,θ′(0)=-(1+θ(0))(NH)orθ(0)=1(CWT)orθ′(0)=-1(CHF),h=-12f′′(0)ony=0,f′→32,θ→0,h→0asy→∞,
where primes denote differentiation with respect to y.The quantities of physical interest are the skin friction coefficient,Cf and the local wall temperature distribution, θw(x), which are given by(17)Re1/2Cf=x(1+K2)(∂2f∂y2)y=0,θw(x)=-1-(∂θ∂y)y=0(NH),
where  Cf=τw/ρU∞2 and  τw=μ(∂u̅/∂y̅)y̅=0 is the wall skin friction or wall shear stress.


## 3. Solution Procedure

Equations in (13) subject to the boundary conditions (14) are solved numerically using the Keller- box method as described in the books by Na [19] and Cebeci and Bradshaw [20]. The solution is obtained in the following four steps.(i)
Reduce (13) to a first-order system.(ii)
Write the difference equations using central differences.(iii)
Linearize the resulting algebraic equations by Newton’s method, and write them in the matrix-vector form.(iv)
Solve the linear system by the block tridiagonal elimination technique (see Salleh et al. [25] for the details of this method).


## 4. Results and Discussion

Equations in (13) subject to the boundary conditions (14) are solved numerically using an efficient implicit finite-difference scheme with various parameters, namely the mixed convection or buoyancy parameter λ, the material or micropolar parameter K and the Prandtl number Pr. In this paper, the numerical solutions start at the lower stagnation point of the sphere, x≈0, as given in (15) subjects to the boundary conditions (16) and proceed round the sphere up to the separation point xs. The values of the skin friction coefficient, Cf and the local wall temperature distribution, θw(x) have been obtained at different positions 0≤x≤1200 for following values of the material parameter K=0 (Newtonian fluid), 1, 2, and 3 (micropolar fluid), with various values of the mixed convection parameter λ.  It should be noticed that the present results were obtained up to the value of x=1200, as given by Nazar et al. [6, 8]. Values of Pr considered are Pr = 0.7, 1, and 7 which correspond to air, electrolyte solution and water, respectively.Due to the decoupled boundary layer equation (15) when the mixed convection parameter λ=0andK=0, there is a unique value of the reduced skin friction coefficient, f′′(0)=2.4104 for all Prandtl number Pr, which is in good agreement with the value f′′(0)=2.4151 found by Nazar et al. [6–8]. The values of f′′(0), -θ′(0) and θ(0) for the cases of CWT, CHF, and NH are shown in Tables 1, 2, and 3, respectively. Some results reported by Nazar et al. [6, 8] for the cases of CWT and CHF and Salleh et al. [17] for the case of NH in viscous fluid (K=0) are also included in these tables. It is found that the agreement between the previously published results with the present ones is very good. We can conclude that this method works efficiently and we are also confident that the results presented here are accurate.Table 1
Values off′′(0)and-θ′(0) for various values of λ with Pr=7 and K=1 (CWT).
λf′′(0)-θ′(0)Nazar et al. [8]PresentNazar et al. [8]Present-5.01.05491.04621.58341.5681-2.01.44651.44301.67261.65810.01.70421.70201.72751.71273.02.10182.10891.80811.794410.02.95682.95541.96171.9434Table 2
Values off′′(0)andθ(0) for various values of λ with Pr=0.7 and K=1 (CHF).
λf′′(0)θ(0)Nazar et al. [6]PresentNazar et al. [6]Present-3.00.33830.30351.65921.6781-2.00.97900.97251.46741.47080.01.71351.71121.32551.32742.02.29412.29231.24551.247110.04.00303.99111.08891.0909Table 3
Values off′′(0)andθ(0) for various values of λ with fixed Pr=1 and K=0 (viscous fluid) (NH).
λf′′(0)θ(0)Salleh et al. [17]PresentSalleh et al. [17]Present0.053.54853.548565.426865.42681.05.21615.21628.73648.73675.08.12738.12733.94953.94957.09.11519.11513.40473.404710.010.373010.37302.92582.9258Tables4(a) and 4(b) show comparisons on the values of  θ(0) and -θ′(0) as well as f′′(0), respectively, for different boundary conditions (CWT, CHF, and NH cases) with fixed Pr=7 and K=1 (micropolar fluid). Table 4(a) shows the values of the heat transfer coefficient -θ′(0) for the case of CWT, the wall temperature θ(0) for the case of CHF, and the wall temperature θ(0) or the heat transfer coefficient -θ′(0) for the case of NH. It is noticed from Table 4(a) that for the case of CWT, as λ increases, the heat transfer coefficient -θ′(0) also increases. This table also shows that for the case of CHF, as λ increases, the wall temperature θ(0) decreases. However, for the case of NH, as λ increases, both the heat transfer coefficient and the wall temperature decrease. The trend for NH case is similar to the CHF case but different from the CWT case. On the other hand, Table 4(b) shows the values of the skin friction coefficient f′′(0) for the cases of CWT, CHF, and NH. It shows that for all the three cases, as λ increases, the skin friction coefficient f′′(0) also increases. The trend for NH case is similar to the CWT and CHF cases, but the values of f′′(0) are closer to the CWT case than the CHF case.Table 4(a)
Comparison on the values ofθ(0) and -θ′(0) from the present solution of (15) for different boundary conditions with Pr=7 and K=1 (micropolar fluid).
λCWTCHFNH-θ′(0)θ(0)θ(0)-θ′(0)-2.01.66700.5878——-1.01.69800.5847——0.01.72720.57901.27482.27482.01.78050.56881.24102.24104.01.82860.55981.17502.17506.01.87250.55171.11902.11908.01.91590.54451.07452.074510.01.95090.53791.03782.037820.02.11010.51150.91651.9165(b)
Comparison on the values off′′(0) from the present solution of (15) for different boundary conditions with Pr=7 and K=1 (micropolar fluid).
λCWTCHFNHf′′(0)f′′(0)f′′(0)-2.01.47481.6244—-1.01.61791.6761—0.01.75681.75681.75682.02.02361.91022.11764.02.27852.05472.39376.02.52352.19182.63678.02.76032.32272.856810.02.99012.44833.059720.04.05903.01603.9136The values off′′(0) and θ(0) for various values of λ, K and Pr are shown in Tables 5 and 6, respectively. From Table 5, it is found that as  K increases, the values of f′′(0) are higher and θ(0) are lower for K=1 than those for K=2 and 3. On the other hand, from Table 6, it can be seen from this table that the values of f′′(0) and θ(0) are higher for Pr = 0.7 than those for Pr = 1 and 7.Table 5
Values off′′(0)andθ(0) for various values of λ and K with fixed Pr=7 (NH).
λK=1K=2K=3f′′(0)θ(0)f′′(0)θ(0)f′′(0)θ(0)0.051.76621.37361.45971.56771.27771.73821.01.93321.29881.59801.49701.39901.65515.02.49501.15272.05831.30551.79921.43337.02.72681.10192.24691.24501.96231.364310.03.03861.03782.50011.17552.18071.2853Table 6
Values off′′(0)andθ(0) for various values of λ and Pr with fixed K=2 (NH).
λPr=0.7Pr=1Pr=7f′′(0)θ(0)f′′(0)θ(0)f′′(0)θ(0)0.055.3893582.99663.3508269.98921.45971.56771.05.934533.97674.010018.97001.59801.49705.07.51889.77015.54756.58722.05831.30557.08.12217.83126.09225.46272.24691.245010.08.90996.28426.79084.52842.50011.1755Tables7 and 8 show the values of Cf and θw(x) for Pr=0.7, K=1 (micropolar fluid) and various values of λ, respectively. It can be seen from these tables that as the parameter λ increases, Cf increases while θw(x) decreases.  The actual value of λ=λk, which first gives no separation is difficult to determine exactly as it has to be found from the equations. However, the numerical solution indicates that the value of  λk which first gives no separation lies between 0.05 and 1 for fixed Pr=0.7 and K=1.Table 7
Values of the local skin friction coefficientCf for Pr=0.7, K=1 and various values of λ(NH).
xλ0.0050.020.05157000.00000.00000.00000.00000.00000.00001002.22612.22802.23192.29952.51902.58622004.67694.68034.68704.80015.19155.31113007.07507.07977.08897.24197.79207.96024009.46159.46739.47899.667610.365710.579750011.690711.697611.711311.933812.771813.196660013.755313.763113.778814.031215.080715.812070015.753415.762215.779616.060017.450518.288680017.467917.477517.496517.809619.486920.370490019.026019.036219.056619.433521.268522.1885100020.215520.237020.693422.575823.5257110021.182321.696323.555724.5200120022.153223.941824.8979Table 8
Values of the wall temperature distributionθw(x) for Pr=0.7, K=1 and various values of λ (NH).
xλ0.0050.020.05157003972.1601997.6703402.734424.93857.78926.340410014946.04233741.49301500.588778.259618.048213.354520015962.96453995.30361601.776382.918818.818513.867830016578.98764149.11611663.146185.800719.325314.214640017114.56424282.91421716.587588.370519.809614.554250017663.74274420.16911771.457591.059220.342414.866760018250.94454566.97091830.179393.973320.901915.061070018940.41944739.37551899.169897.425721.448415.322280019744.13354940.37231979.6231101.465722.051515.727690020692.22735177.49642074.5533106.177222.809516.295510005460.62042187.9111111.704523.787717.040211002324.9537118.337025.044417.98911200126.516126.659919.1931Tables9 and 10 present the values of the local skin friction coefficient Cf and the wall temperature distribution θw(x) for λ=1, K=1 and various values of Pr, respectively. It is found that, for fixed λ=1 and K=1, as Pr increases, both the Cf and θw(x) decreases. On the other hand, for fixed Pr, as x increases, that is, from the lower stagnation point of the sphere, x≈0, and proceeds round the sphere up to the point, x=1200, both the values of Cf and θw(x) increase.Table 9
Values of the local skin friction coefficientCf for λ=1, K=1 and various values of Pr(NH).
xPr0.717000.00000.00000.00001002.29951.50530.14902004.80013.16090.38063007.24194.75980.52114009.66796.37230.766850011.93387.85490.888760014.03129.28171.065770016.060010.77711.320780017.809612.04951.413590019.433513.20691.6429100020.693414.02391.6813110021.696314.67901.8646120022.153214.91311.9266Table 10
Values of the wall temperature distributionθw(x) for λ=1, K=1 and various values of Pr(NH).
xPr0.7170024.938513.78101.298810078.259647.57983.491320082.918850.55833.842330085.800752.39054.060340088.370554.01044.216650091.059255.69984.381460093.973357.45984.534870097.425759.25234.6891800101.465761.17984.8683900106.177263.49155.07581000111.704566.39025.33831100118.337070.06175.66351200126.516174.75206.0838The velocity and temperature profiles near the lower stagnation point of the sphere,x≈0,  for various values of  K with fixed Pr=7 and λ=1 are plotted in Figures 2 and 3, respectively. From these figures, it is shown that as K increases, the velocity profiles decrease and the temperature profiles increase.Figure 2
Velocity profiles for various values ofK with fixed Pr=7andλ=1 (NH).Figure 3
Temperature profiles for various values ofK with fixed Pr=7andλ=1 (NH).Figures4 and 5 illustrate the velocity and temperature profiles near the lower stagnation point  of the sphere, x≈0, for various values of λ with fixed Pr=1 and K=1, respectively. We found that the velocity profiles increase while the temperature profiles decrease when the mixed convection parameter λ increases. Also, it is noticed that there are overshoots of the velocity profiles when λ>1 where these overshoots take place higher for λ=10 than for λ=1.Figure 4
Velocity profiles for various values ofλ with fixed Pr=1andK=1 (NH).Figure 5
Temperature profiles for various values ofλ with fixed Pr=1andK=1 (NH).Figures6 and 7 illustrate the velocity and temperature profiles near the lower stagnation point of the sphere, x≈0, for various values of  Pr with fixed K=2 and λ=1, respectively.  It is seen that, as the Prandtl number Pr increases, both the velocity and temperature profiles decrease. At large Pr(≫1), the thermal boundary layer is thinner than at smaller Pr. This is because for small values of Pr(≪1), the fluid is highly conductive. Physically, if Pr increases, the thermal diffusivity decreases and this phenomena leads to the decreasing of energy ability that reduces the thermal boundary layer. It is also noticed that there are overshoots of the velocity profiles when Pr≤1 where these overshoots are higher for Pr = 0.7 than for Pr = 1.Figure 6
Velocity profiles for various values ofPr with fixed K=2andλ=1 (NH).Figure 7
Temperature profiles for various values ofPr with fixed K=2andλ=1 (NH).


## 5. Conclusions

In this paper, we have numerically studied the problem of mixed convection boundary layer flow from a solid sphere in a micropolar fluid, generated by Newtonian heating. It is shown that the mixed convection or buoyancy parameterλ, the material or micropolar parameter K, and the Prandtl number Pr affect the flow and heat transfer characteristics. We can conclude that (for the case of NH):(i)
whenPr and K are fixed, an increase in parameter λ leads to the decrease of both the heat transfer coefficient -θ'(0) and the wall temperature θ(0);(ii)
whenPr and λ are fixed, an increase in parameter K leads to the decrease of the skin friction coefficient f′′(0) and the increase of the wall temperature θ(0);(iii)
whenPr and K are fixed, an increase in parameter λ leads to the increase of the skin friction coefficient f′′(0) and the decrease of the wall temperature θ(0);(iv)
whenK and λ are fixed, an increase in parameter Pr leads to the decrease of the skin friction coefficient f′′(0) and the wall temperature θ(0);(v)
near the lower stagnation point of the sphere, the velocity profiles decrease while the temperature profiles increase when the parameterK increases;(vi)
near the lower stagnation point of the sphere, the velocity profiles increase while the temperature profiles decrease when the parameterλ increases;(vii)
near the lower stagnation point of the sphere, both the velocity and temperature profiles decrease when the Prandtl numberPr increases.



---

*Source: 736039-2010-03-02.xml*