Next Article in Journal
A Shack-Hartmann Sensor for Single-Shot Multi-Contrast Imaging with Hard X-rays
Next Article in Special Issue
Stability Analysis of Mixed Convection Flow towards a Moving Thin Needle in Nanofluid
Previous Article in Journal
Control of Carrier Concentration by Ag Doping in N-Type Bi2Te3 Based Compounds
Previous Article in Special Issue
Stability Analysis of Stagnation-Point Flow in a Nanofluid over a Stretching/Shrinking Sheet with Second-Order Slip, Soret and Dufour Effects: A Revised Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MHD Steady/Unsteady Porous Boundary Layer of Cu–Water Nanofluid with Micropolar Effect over a Permeable Surface

by
Kamal Raslan
1,
Selim Mohamadain
2,
Mohamed Abdel-wahed
3,* and
Elsayed M. Abedel-aal
1,*
1
Mathematics Department, Faculty of Science, Al-azhar University, Cairo 002, Egypt
2
Mathematics Department, Faculty of Science, Tanta University, Tanta 002, Egypt
3
Basic Science Department, Faculty of Engineering, Benha University, Cairo 002, Egypt
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2018, 8(5), 736; https://doi.org/10.3390/app8050736
Submission received: 5 March 2018 / Revised: 26 March 2018 / Accepted: 3 April 2018 / Published: 7 May 2018
(This article belongs to the Special Issue Nanofluids and Their Applications)

Abstract

:
This work provides a mathematical model for the cooling process of a moving surface, in the presence of a uniform external magnetic field and thermal radiation, through a porous medium by using a weak concentration micropolar nanofluid. The model—based on the conservation equations of the unsteady case in the momentum and thermal boundary layer—takes into consideration the effect of the suction process. The conservation equations were transformed into ordinary differential equations using similar transformation techniques. The equations were solved numerically for the general case and analytically for the steady case. The rate of heat transfer, couple shear stress, and surface shear stress are deduced. We discuss the impact of these physical characteristics on the mechanical properties of the surface that will be cooled.

1. Introduction

We enhanced the thermal conductivity of the coolant based on its essential effects on many engineering applications, such as the heat treatment of metals [1]. The aim of the cooling process is to improve the metals’ mechanical properties, such as ductility and strength. Creating more useful metals is one of the most important aims of researchers. In this process, the metal is heated to a certain temperature and then cooled in a series of specific operations. Nanofluids have properties that make them useful in many heat transfer applications, such as microelectronics, fuel cells, pharmaceutical processes, hybrid powered engines, engine cooling/vehicle thermal management, domestic refrigerators, heat exchangers, and the reduction of flue gas temperature in boilers [2]. Compared to the base fluid, thermal conductivity and the convective heat transfer coefficient were enhanced. In 1995, the authors of [3] proposed the idea of enhancing the thermal conductivity of the fluids using these nanoparticles. This was the first push for researchers to conduct studies that investigated the effects of nanoparticle types, sizes, and concentrations on the heat transfer characteristics and took into consideration some of the effects such as the thermal radiation, variable heat flux, magnetic field, and suction/injection processes [4]. We studied the magnetohydrodynamics (MHD) free convection of Al2o3–water nanofluid by considering the thermal radiation problem [5]. This study examined the flow and heat transfer over an unsteady shrinking sheet with suction in nanofluids [6]. The natural convection flow of a nanofluid over a linear stretching sheet, in the presence of a magnetic field, was discussed [7]. The mixed convection boundary layer flow from a horizontal circular cylinder was embedded in a porous medium, which was filled with a nanofluid [8]. We observed the boundary layer flow over a stretching/shrinking surface under an external uniform shear flow with a convective surface boundary condition in a nanofluid [9]. The effect of a magnetic field on the flow and heat transfer of a nanofluid over an unsteady continuous moving surface—and in the presence of suction/injection—was observed [10]. Buongiorno’s model was used to study the heat and mass transfer flow in a thin nanofluid film over an unsteady stretching sheet [11]. The study focused on the effects of thermophoresis and Brownian motion on the boundary layer flow of a nanofluid in the presence of thermal stratification—due to solar energy [12]. The flow and heat elements are transferred, in a nanofluid, over a moving surface with a non-linear velocity and variable thickness and in the presence of thermal radiation [13]. The boundary layer flow was maneuvered over a moving surface, which was embedded into a nanofluid, in the presence of a magnetic field and suction/injection [14]. We analyzed the effects of thermal radiation and heat generation on the mechanical properties of an unsteady, continuously moving cylinder in a nanofluid in the presence of suction/injection [15]. The flow and heat transfer was observed over a moving surface in a nanofluid with a non-linear velocity and variable thickness and in the presence of Brownian motion.
Micropolar fluids may represent fluids that consist of rigid, randomly oriented (or spherical) particles that are suspended in a viscous medium, where the deformation of fluid particles is ignored. The microrotation of the suspended microparticles within the fluid is discussed through the theory of micropolar fluids. In this theory, the fluids can be subjected to the surface and body couples, in which the material points in the volume elements can both move around the centers of mass and deform. Many researchers, such as the authors of [16], have investigated the micropolar fluid theory by studying the behavior of the unsteady boundary layer flow of a micropolar fluid near the rear stagnation point of a plane surface [17]. The unsteady mixed convection boundary flow of a micropolar fluid with uniform suction/injection over a stretching/shrinking sheet was analyzed [18]. We investigated the flow and heat transfer in a micropolar fluid along a stretching surface. According to the authors of [19], the effect of the thermal radiation interaction of the boundary layer flow of a micropolar fluid along a heated plate was embedded in a porous medium with variable heat flux [20]. The effect of variable heat flux on the heat transfer characteristics of a stretching surface in micropolar fluids was studied [21]. This study examined the influence of suction/injection and heat generation on the MHD flow of a micropolar fluid along a stretched permeable surface [22]. The heat transfer characteristic of a micropolar fluid along a heated vertical porous plate with variable suction and heat flux was compared to the Newtonian fluid [23]. The effects of the thermal radiation and heat generation were observed with respect to the flow and heat transfer characteristics over an unsteady moving surface in a micropolar fluid in the presence of suction/injection [24]. The suction/injection effects were studied on the micropolar fluid flow along a continuously moving plate in the presence of radiation [25]. The effects of microrotation and nanoparticle volumetric fraction on the fluid flow and heat transfer behavior in the boundary layer were observed [26]. We studied the effect of the microrotation of nanoparticles on the momentum boundary layer of a permeable stretching sheet in micropolar nanofluids with suction [27]. The flow and heat transfer of a micropolar fluid on a porous stretching sheet was observed [28]. The study recorded the effects of Soret and the non-uniform heat source/sink on the MHD non-Darcian convective flow along a stretching sheet in a micropolar fluid with the radiation effect [29]. We examined the effects of radiative and Joule heating on the MHD flow of a micropolar fluid with a partial slip and convective boundary condition problem [30], as well as on the micropolar nanofluid flow with MHD and viscous dissipation effects towards a stretching sheet with a multimedia feature [31]. We used the Chebyshev Spectral Newton Iterative Scheme to analyze the hydromagnetic Hiemenz flow of a micropolar fluid over a non-linear stretching/shrinking sheet [32], and investigated the MHD orthogonal stagnation-point flow of a micropolar fluid in the presence of a magnetic field that is parallel to the velocity at infinity [33]. The Lorentz force effect was observed on a mixed convection micropolar flow in a vertical conduit [34], along with the flow and heat transfer of a weak concentration micropolar nanofluid over a steady/unsteady moving surface [35]. The Molecular Dynamics Study of Thermodynamic Properties of Nanoclusters for Additive Manufacturing was considered [36]. We observed the nanoscale liquid film that was sheared between strong wetting surfaces and studied the effects of the interface region on the flow [37], as well as the molecular dynamics simulations of oscillatory Couette flows with slip boundary conditions [38]. Viscous heating in the nanoscale shear driven liquid flows, as well as [39] the molecular dynamics of bimetallic nanoparticles—including the case of AuxCuy alloy clusters [40]—was studied. The thermal resistance at a liquid–solid interface, which is dependent on the ratio of thermal oscillation frequencies, along with [41] the molecular simulation of thermal transport across hydrophilic interfaces, was examined.
This work presents an investigation of the effect of microrotation of the nanoparticles along a moving plate, with time-dependent velocity and temperature, through a boundary layer that is subjected to thermal radiation and a suction process, in the presence of a uniform magnetic field through a porous medium. The study considered two cases, namely the steady weak concentration—which was solved analytically using the hypergeometric function—and the unsteady motion—which was solved numerically.

2. Formulation of the Problem

The study analyzed a steady/unsteady flow of an incompressible micropolar nanofluid over a permeable moving surface under the effect of thermal radiation with heat flux q r and uniform suction S w . This was analyzed in the presence of thermal radiation and a uniform magnetic field of strength B0 through a porous medium. The temperature of the outer surface is T w , while T is the temperature of the free stream that is far away from the surface (see Figure 1). It is assumed that the nanoparticles have a uniform size and shape. In addition, both the fluid phase and the nanoparticles are assumed to be in thermal equilibrium state. The nanoparticles are also assumed to be very small in size, which means that the slip velocity between the phases is presumed to be negligible. The conservation equations for the unsteady micropolar nanofluid flow are presented below [23,26].
u x + v y = 0 ,
u t + u u x + v u y = ( μ n f + k ρ n f ) 2 u y 2 + k ρ n f ω y ν K p u σ B ο 2 ρ u ,
ω t + u ω x + v ω y =   γ n f ( ρ n f )   j 2 ω y 2 k ( ρ n f )   j ( 2 ω + u y ) ,   and
T t + u T x + v T y = α n f 2 T y 2 1 ( ρ C p ) n f q y .
These equations are subjected to the following boundary conditions.
u = U w , v = S w ,   T = T w ,   ω = n u y a t y = 0 , u 0 ,   ω 0 ,   T = T   a s   y ,
where n is the boundary parameter, with 0     n     1 . The case corresponds to n results in the vanishing of the anti-symmetric part of the stress tensor and represents a weak concentration. In this paper, this was the only case that was considered.
The spin gradient viscosity is defined below [34].
γ n f = ( μ n f + k 2 ) j = μ f ( 1 ( 1 ϕ ) 2.5 + K 2 ) j .
The properties of the nanofluid are defined below [34] (see Table 1).
μ n f = ( μ f ( 1 ϕ ) 2.5 ) , ρ n f = ( 1 ϕ ) ρ f + ϕ ρ s ,     α n f = k n f ( ρ C p ) n f ,
( ρ C p ) n f = ( 1 ϕ ) ( ρ C p ) f + ϕ ( ρ C p ) s , k n f k f = ( k s + 2 k f ) 2 ϕ ( k f k s ) ( k s + 2 k f ) + ϕ ( k f k s ) .
Using Rosseland’s approximation for radiation, the radiative heat flux is simplified in Equation (8).
q r = 4 σ 3   α T 4 y .
Hence, by expanding T 4 in a Taylor series about T and neglecting the higher order terms, we get Equation (9).
T 4 4 T T 3 3 T 4 .
Using Equations (8) and (9) in the energy Equation (4), we get Equation (10).
T t + u T x + v T y = α n f 2 T y 2 + ( 16 σ   T 3 3 α ( ρ C p ) n f ) 2 T y 2 .
Assuming that the velocity and the fluid temperature are in the form below.
U w ( x , t ) = b x 1 α t ,   T w ( x , t ) = T + c x 1 α t ,
where b , c , and α are constants with b , c , and α > 0 .
By introducing the dimensionless functions f ,   g , and θ , and the similarity variable η ,
η = b υ ( 1 α t ) y ,   ψ = b υ ( 1 α t ) x f ( η ) ,   ω = b 3 υ ( 1 α t ) 3 x g ( η ) , θ ( η ) = T T T w T ,
where ψ ( x , y ) is a stream function that is defined as u = ψ / y and v = ψ x , which both satisfies the mass conservation in Equation (1) and uses the dimensionless functions in Equation (12), we get the following equations.
( ( 1 ϕ ) 2.5 + K β ) f + f f f 2 A ( f + η 2 f ) + K   B β g ( M + kp )   f = 0 ,
( ( 1 ϕ ) 2.5 + 1 2 K β ) g + f g f g A ( 3 2 g + η 2 g ) K β ( 2 g + f ) = 0 ,   and
( L + 4 R 3 H ) θ + P r ( f θ f θ ) A   P r ( θ + η 2 θ ) = 0 .
The boundary condition (5) is calculated in Equation (16).
f ( 0 ) = S w ,   f ( 0 ) = 1 ,   g ( 0 ) = n f ( 0 ) ,   θ ( 0 ) = 1 ,   f ( ) = 0 ,   g ( ) = 0 ,   θ ( ) = 0 ,
where
β = ( 1 ϕ + ϕ ρ s ρ f ) ,   H = ( 1 ϕ + ϕ ( ρ C p ) s ( ρ C p ) f ) ,   B = ν ( 1 α t ) b   j ,   L = ( k n f k f ) H , A = b a , M = σ B ο 2 ( 1 a t ) b ρ f ,   k p = ν ( 1 a t ) b K p ,   K = k μ ,   R = 4 σ T 3 α k f ,   P r = ν f ( ρ C p ) f k f .

3. Solutions and Results

3.1. Steady Weak Concentration Boundary Layer

For the steady weak concentration boundary layer, we considered the unsteadiness parameter ( A ) = 0 and n = 0.5 . Thus, Equations (13)–(15) became the equations below.
( ( 1 ϕ ) 2.5 + K β ) f + f f f 2 + K   B β g ( M + kp ) f = 0 ,
( ( 1 ϕ ) 2.5 + 1 2 K β ) g + f g f g K β ( 2 g + f ) = 0 ,   and
( L + 4 R 3 H ) θ + P r ( f θ f θ ) = 0 .
The boundary condition (16) became Equation (20).
f ( 0 ) = S w , f ( 0 ) = 1 , g ( 0 ) = 1 2 f ( 0 ) , θ ( 0 ) = 1 ,
f ( ) = 0 ,   g ( ) = 0 ,   θ ( ) = 0 .

3.1.1. Exact Solution of the Momentum Boundary Layer Problem

In this section, general solutions for Equations (17) and (18) are suggested in the form below.
f ( η ) = S w + 1 m ( 1 e m η )   a n d   g ( η ) = m 2 e m η ,
which are satisfied by the following conditions.
f ( 0 ) = S w ,   f ( 0 ) = 1 ,   g ( 0 ) = 1 2 f ( 0 ) ,   θ ( 0 ) = 1 ,
f ( ) = 0 ,   g ( ) = 0 ,   θ ( ) = 0 .
Upon substituting Equation (21) into Equations (17) and (18), and using the boundary conditions from Equation (22), we get the following equation.
m = β ( 1 ϕ ) 2.5 2 ( 1 + 1 2 K ( 1 ϕ ) 2.5 ) { S w + ( S w ) 2 + 4 ( 1 + 1 2 K ( 1 ϕ ) 2.5 ) ( 1 + M + k p ) β ( 1 ϕ ) 2.5   } .
Therefore, the linear velocity and angular velocity components take the form below.
u = ( b x   e m η ) ,   v = b υ f   [ S w + 1 m ( 1 e m η ) ] .

3.1.2. Exact Solution of the Momentum Boundary Layer Problem

By substituting Equation (19) into Equation (21), one can obtain the following ordinary differential equation.
θ + ( 3 H P r 4 R + 3 H L ) [ ( S w + 1 m ( 1 e m η ) ) θ ( e m η ) θ ] = 0 ,
with boundary conditions that are given by the following equation.
θ ( 0 ) = 1 ,   θ ( ) = 0 .
Equation (25) is a homogenous linear ordinary differential of the second order with a variable coefficient, which can be solved by introducing the change variable.
ε = λ e m 2 m η ,
where λ = ( 3 H P r 4 R + 3 H L ) . Equation (25) can then be written as below.
ε d 2 θ d ε 2 + [ 1 ξ ε ] d θ d ε + θ = 0 ,
where ξ = λ m 2 ( 1 + m   S w ) , with the following boundary conditions.
θ ( λ m 2 ) = 1   and   θ   ( ) = 0 .
Equation (28) is similar to Kummer’s differential equation, which has a Kummer confluent hypergeometric function, 1 F 1 , as a solution.
θ ( ε ) = ( ε m 2 λ ) ξ 1 F 1 ( ξ 1 ; ξ + 1 ; ε ) 1 F 1 ( ξ 1 ; ξ + 1 ; λ m 2 ) .
The solution of Equation (30) in terms of the similarity variable ( η ) can be found as described below.
θ ( η ) = ( e m η ) ξ 1 F 1 ( ξ 1 ; ξ + 1 ; ε ) 1 F 1 ( ξ 1 ; ξ + 1 ; λ m 2 ) .
In addition, the surface temperature gradient is described below.
θ ( 0 ) = ( m ξ ) + ( λ ( ξ 1 ) m ( ξ + 1 ) ) 1 F 1 ( ξ ; ξ + 2 ; λ m 2 ) 1 F 1 ( ξ 1 ; ξ + 1 ; λ m 2 ) .

3.2. Unsteady Weak Concentration Boundary Layer

The equations for the unsteady case at n = 0.5 , Equations (13)–(15), and the boundary conditions, Equation (16), are solved numerically. In order to solve this system, we first need to find the missing initial conditions f ( 0 ) ,   g ( 0 ) , and   θ ( 0 ) by constructing a program for the system as an initial value problem with the missing initial conditions. In order to solve this system, we have to suggest a proper value of η η . Following this, the missing and necessary values are calculated, using the root technique, by the code. Once the values of f ( 0 ) ,   g ( 0 ) , and   θ ( 0 ) have no change with an increase of η (i.e., the boundary conditions f ( η ) = 0 ,   g ( η ) = 0 ,   θ ( η ) = 0 are satisfied at the suggested value of η ), the system is ready to be solved, using the classical Runge–Kutta method from the fourth order with a step size of Δ η = 0.01 . To support the steady case results that were obtained by this method, the values in Table 2 are compared with the results for C f R x and g ( 0 ) . The numerical solution was reported in Ezzah et al. [26] and the exact solution was obtained in this study.

4. Physical Aspect of the Problem

For certain important scientific purposes, it is useful to investigate the flow characteristics. The functions f ( 0 ) ,   g ( 0 ) ,   and   θ ( 0 ) allow us to determine the skin friction coefficient, the dimensionless couple stress, and the Nusselt number, which indicate the surface shear stress, the couple stress, and the rate of heat transfer, respectively. It is useful to mention the final quality of the surface that is controlled by the surface shear stress and couple stress. In addition, hardness, flexibility, and ductility are important mechanical properties of the surface and depend on the rate of heat transfer from the hot surface (the cooling rate).

4.1. Surface Shear Stress

The shear stress is defined below.
τ w = [ ( μ n f + k ) u y + k ω ] .
The local skin friction is defined below.
C f = 2 τ w ρ n f U w 2 = 2 R e ( 1 β ( 1 ϕ ) 2.5 + K ( 1 n ) β ) f ( 0 ) .

4.2. Surface Couple Stress

M w = γ n f [ ω y ] .
The dimensionless couple stress is defined below.
M x = 2 M w ρ n f υ n f U w = ( 2 + K ( 1 ϕ ) 2.5 2 ) g ( 0 ) .

4.3. Surface Heat Flux

q w = k n f [ T y ] .
The local Nusselt number is defined below.
N u = x q w k f ( T w T ) = k n f k f [ R e θ ( 0 ) ] .

5. Discussion

The following discussion is based upon a numerical solution—supported by a solution form that was obtained for the nonlinear system—describing the micropolar nanofluid boundary layer over a moving flat surface, which is an emulation of the cooling process. The model considers the steady and unsteady boundary layer, which contains Cu (Copper) and is subjected to an external magnetic field, thermal radiation, and suction through porous media. The linear and angular velocities, along with the temperature profiles, are plotted and presented in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8. In addition, the values of the local skin friction, couple stress, and local Nusselt number for Cu particles are presented in Table 2 and Table 3.

5.1. Material Parameter

Newtonian or non-Newtonian fluids depend on the ratio between the vortex viscosity and the dynamics viscosity, which is called the material parameter (K), so that K = 0 represents Newtonian fluids when there is no vortex viscosity. Figure 2a,b characterize the linear and angular velocity profiles under the variation of the material parameter. We can see that the linear velocity increases when this parameter rises but the angular velocity decreases near the surface—because of microrotation—before reversing and increasing farther from the surface. The changes in the surface skin friction, couple stress, and Nusselt number are presented in Table 3 and Table 4, and physically indicate the surface shear stress, the couple stress, and the rate of heat transfer. In general, it is observed that the presence of the viscosity vortex ( K 0 ), due to the microstructures microrotation, increases near the surface shear stress and the rate of heat transfer. Physically, the weak angular velocity near the surface resulted from a weak concentration of the microstructures.

5.2. Steady/Unsteady Boundary Layer

The value of the unsteadiness parameter (A) indicates the type of motion such that for A = 0, the problem is reduced to the steady case, and A ≠ 0 indicates the unsteady motion. Moreover, the positive sign of this unsteadiness parameter indicates accelerated motion and the negative sign indicates deceleration. Figure 3a,b indicate the effect of the unsteadiness parameter on the linear and angular velocities. It is observed that increasing this parameter leads to a decrease in the linear and angular velocities when it is close to the surface. This effect is reversed when it is far away from the surface. Figure 3c indicates that the unsteadiness parameter has a finite effect on the boundary layer temperature. It is important to mention that only the steady motion, indicated by a black line, is plotted in Figure 3a–c. On the other hand, Table 3 and Table 4 present a comparison of the physical characteristics of the steady and unsteady boundary layers through the values of the surface skin friction, couple stress, and Nusselt number. In the case of unsteady motion, the surface skin friction value is larger for Newtonian fluids (K = 0) than it is for non-Newtonian fluids (K ≠ 0). This means that the skin friction is more affected by Newtonian fluids than by non-Newtonian fluids. Additionally, in the case of unsteady motion, we observe the same effect for the couple stress and the local Nusselt number. It is observed that increasing the surface heat flux, caused by the increase of the Nusselt number, will accelerate the rate of cooling.

5.3. Nanoparticle Concentration

The ratio ϕ represents the nanoparticle concentration within the fluid. It is clear that for ϕ = 0 , the problem is converted to the regular micropolar fluid. This implies that this ratio controls the dynamic viscosity and the thermal conductivity of the fluid, such that the linear velocity of the flow increases by increasing the nanoparticle concentration (see Figure 4a). On the other hand, Figure 4b indicates that the presence of the nanoparticles increases the angular velocity when it is close to the surface. As a result of presence of the nanoparticles, Figure 4c indicates that the boundary layer temperature increases when the thermal conductivity grows. Table 3 illustrates the influence of the nanoparticle concentration on the surface skin friction, couple stress, and Nusselt number. It is observed that the nanoparticle concentration is more effective on the skin friction when there is no suction effect than it is when there is a suction effect. Additionally, it is clear that the nanoparticle concentration is more effective on the couple shear stress when there is a suction effect than it is when there is no suction effect. Moreover, according to the steadiness parameter and the material parameter, the rate of heat transfer increases when the local Nusselt number rises. In general, increasing the skin friction or the couple stress has a negative effect on the quality of the outer surface of the plate. When the Nusselt number increases, the rate of heat transfer from the surface rises, which has a positive effect on certain of the mechanical properties, such as hardness and strength.

5.4. Suction Parameter Effect

Additionally, the effect of the suction parameter (Sw) on the velocities and the temperature are presented in Figure 5a–c, respectively. Figure 5a,b indicate that the linear and angular velocities decrease with increasing values of Sw. The fluid velocities are found to decrease with the increase in Sw, which means that suction causes a dragging force that decreases the fluid velocity in the boundary layer. Figure 5c indicates that the temperature in the boundary layer decreases with increasing values of Sw, which causes a decrease in the rate of heat transfer from the fluid to the surface. The thermal boundary layer is thinner in the case of suction.

5.5. Permeability Effect

Figure 6a depicts the effect of the permeability of the porous medium parameter (kp) on the velocity distribution profile. It is clear that as kp increases, the linear velocity decreases along the boundary layer thickness. This is expected, since the holes of the porous medium become larger. Figure 6b indicates the influence of kp on the angular velocity distribution. It is obvious that as kp increases, the angular velocity decreases.

5.6. Magnetic–Radiation Effect

Figure 7a,b indicate the effect of the magnetic field parameter (M) on the velocities. With the increase of M, the linear velocity and the angular velocity of the flow decrease. Figure 8 illustrates the temperature profiles for different values of the radiation parameter (R). The figure indicates that as the value of R increases, the temperature profiles increase.

6. Conclusions

The heat transfer rate decreases in the case of the unsteady boundary layer when compared to the case of the steady boundary layer. The rate of heat transfer increases more for non-Newtonian fluids than it does for Newtonian fluids. Furthermore, Newtonian fluids are more effective for increasing the surface shear stress and couple stress than non-Newtonian fluids. Moreover, the presence of nanoparticles increases the rate of heat transfer for the non-Newtonian boundary layer when compared to Newtonian fluids. Lastly, the rate of heat transfer, the surface shear stress, and the couple shear stress decrease during the suction process.

Acknowledgments

The authors would like to express their very sincere thanks to the editor and the reviewers for their valuable comments and suggestions.

Author Contributions

Kamal Raslan and Selim Mohamadain suggested the problem and Mohamed Abdel-wahed with Elsayed M. Abedel-aal analyzed the data and Elsayed M. Abedel-aal wrote it.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

x ,   y Cartesian coordinate system
b   and   c Constants
u ,   v Velocity components along the x-, y-directions
ω Angular velocity
t Time
ψ ( x , y ) The stream function
U w ,   V w Velocity of the stretching sheet along x- and y-directions
T Fluid temperature
n boundary parameter
k vortex viscosity
K p Porous medium parameter
S w Suction parameter
K Material parameter
B 0 Magnetic field
A Unsteadiness parameter
q r Radiative heat flux
α Mean absorption coefficient
j Micro-inertia per unit mass
P r Prandtl number
R Radiation parameter
f, gDimensionless components of velocities
τ w Shear stress
C f Skin friction
M x Dimensionless couple stress
N u Nusselt number
ρ Density of the fluid
γ Spin gradient viscosity
σ Stefen–Boltzmann constant
α n f thermal diffusivity of the nanofluid
k Thermal conductivity
μ Dynamic viscosity
υ kinematic viscosity
η Similarity variable
ϕ Nanoparticles concentration
θ Dimensionless temperature.
Subscripts
f Fluid phase
n f Nanofluid
s Solid particles
w Condition of the wall
Ambient condition

References

  1. Elbashbeshy, E.M.A.; Emam, T.G.; Abdel-Wahed, M.S. Effect of heat treatment process with a new cooling medium (nanofluid) on the mechanical properties of an unsteady continuous moving cylinder. J. Mech. Sci. Technol. 2013, 27, 3843–3850. [Google Scholar] [CrossRef]
  2. Das, S.K.; Choi, S.U.S.; Yu, W.; Pradeep, T. Nanofluids: Science and Technology; John Wiley & Sons: Hoboken, NJ, USA, 2007; ISBN 978-0-470-07473-2. [Google Scholar]
  3. Choi, S.U.S.; Eastman, J.A. Enhancing conductivity of fluids with nanoparticles. ASME Fluid Eng. Div. 1995, 231, 99–105. [Google Scholar]
  4. Sheikholeslami, M.; Hayat, T.; Alsaedi, A. MHD free convection of Al2O3–water nanofluid considering thermal radiation: A numerical study. Int. J. Heat Mass Transf. 2016, 96, 513–524. [Google Scholar] [CrossRef]
  5. Azizah, M.; Syakila, A.; Pop, I. Flow and heat transfer over an unsteady shrinking sheet with suction in nanofluids. Int. J. Heat Mass Transf. 2011, 55, 1888–1897. [Google Scholar]
  6. Hamad, M. Analytical solution of natural convection flow of a nanofluid over a linearly stretching sheet in the presence of magnetic field. Int. Commun. Heat Mass Transf. 2011, 38, 487–492. [Google Scholar] [CrossRef]
  7. Nazar, R.; Tham, L.; Pop, I. Mixed convection boundary layer flow from a horizontal circular cylinder embedded in a porous medium filled with a nanofluid. Transp. Porous Med. 2011, 86, 517–536. [Google Scholar] [CrossRef]
  8. Yacob, N.A.; Ishak, A.; Nazar, R.; Pop, I. Boundary layer flow past a stretching/shrinking surface beneath an external uniform shear flow with a convective surface boundary condition in a nanofluid. Nanoscale Res. Lett. 2011, 6, 1–7. [Google Scholar] [CrossRef] [PubMed]
  9. Elbashbeshy, E.M.A.; Emam, T.G.; Abdel-Wahed, M.S. Effect of magnetic field on flow and heat transfer of a nanofluid over an unsteady continuous moving surface in the presence of suction/injection. Int. J. Appl. Math. 2012, 14, 436–442. [Google Scholar]
  10. Qasim, M.; Khan, Z.; Lopez, R.; Khan, W. Heat and mass transfer in nanofluid thin film over an unsteady stretching sheet using Buongiorno’s model. Eur. Phys. J. Plus 2016, 131, 16. [Google Scholar] [CrossRef]
  11. Anbuchezhian, N.; Srinivasan, K.; Chandrasekaran, K.; Kandasamy, R. Thermophoresis and Brownian motion effects on boundary layer flow of nanofluid in presence of thermal stratification due to solar energy. Appl. Math. Mech. 2012, 33, 765–780. [Google Scholar] [CrossRef]
  12. Elbashbeshy, E.M.A.; Emam, T.G.; Abdel-Wahed, M.S. Flow and heat transfer over a moving surface with non-linear velocity and variable thickness in a nanofluid in the presence of thermal radiation. Can. J. Phys. 2013, 91, 699–708. [Google Scholar]
  13. Elbashbeshy, E.M.A.; Emam, T.G.; Abdel-Wahed, M.S. An exact solution of boundary layer flow over a moving surface embedded into a nanofluid in the presence of magnetic field and suction/injection. Heat Mass Transf. 2014, 50, 57–64. [Google Scholar] [CrossRef]
  14. EElbashbeshy, M.A.; Emam, T.G.; Abdel-Wahed, M.S. The effect of thermal radiation and heat generation on the mechanical properties of unsteady continuous moving cylinder in a nanofluid in the presence of suction or injection. Therm. Sci. 2015, 19, 1591–1601. [Google Scholar] [CrossRef]
  15. Abdel-wahed, M.S.; Elbashbeshy, E.M.A.; Emam, T.G. Flow and heat transfer over a moving surface with non-linear velocity and variable thickness in a nanofluid in the presence of Brownian motion. Appl. Math. Comput. 2015, 254, 1591–1601. [Google Scholar] [CrossRef]
  16. Yain, Y.; Amin, N.; Pop, I. Unsteady boundary layer flow of a micropolar fluid near rear stagnation point of a plane surface. Int. J. Therm. Sci. 2003, 42, 995–1001. [Google Scholar]
  17. Sandeep, N.; Sulochana, C. Dual solution for unsteady mixed convection flow of MHD micropolar fluid over a stretching/shrinking sheet with non-uniform heat source/sink. Eng. Sci. Technol. 2015, 18, 738–745. [Google Scholar] [CrossRef]
  18. Abo-Eldahab, E.M.; El Aziz, M.A. Flow and heat transfer in a micropolar fluid past a stretching surface embedded in a non-Darcian porous medium with uniform free stream. Appl. Math. Comput. 2005, 162, 881–899. [Google Scholar] [CrossRef]
  19. Rahman, M.; Sultana, T. Radiative Heat transfer flow of micropolar fluid with variable heat flux in a porous medium. Nonlinear Anal. Model. Control 2008, 13, 17–87. [Google Scholar]
  20. Ishak, A.; Nazar, R.; Pop, I. Heat transfer over a stretching surface with variable heat flux in micropolar fluids. Phys. Lett. A 2007, 372, 559–561. [Google Scholar] [CrossRef]
  21. Khedr, M.E.M.; Chamkha, A.J.; Bayomi, M. MHD flow of a micropolar fluid past a stretched permeable surface with heat generation or absorption. Nonlinear Anal. Model. Control 2009, 14, 27–40. [Google Scholar]
  22. Rahman, M.; El-tayeb, A.; Rahman, S.M.M. Thermo-micropolar fluid flow along a vertical permeable plate with uniform surface heat flux in the presence of heat generation. Therm. Sci. 2009, 13, 23–36. [Google Scholar] [CrossRef]
  23. Saad, A.A.; Abdel-Wahed, M.S. Heat transfer over unsteady moving surface with heat generation and thermal radiation in micropolar fluid in the presence of suction/injection. Int. J. Energy Technol. 2011, 3, 1–7. [Google Scholar]
  24. L-Arabawy, H.A.M.E. Effects of suction/injection on the flow of micropolar fluid past a continuously moving plate in the presence of radiation. Int. J. Heat Mass Transf. 2003, 46, 1471–1477. [Google Scholar] [CrossRef]
  25. Hussain, S.T.; Nadeem, S.; Haq, R.U. Model-based analysis of micropolar nanofluid flow over a stretching surface. Eur. Phys. J. Plus 2014, 129, 161. [Google Scholar] [CrossRef]
  26. Ezzah, L.A.; Syakila, A.; Pop, I. Flow over a permeable stretching sheet in micropolar nanofluids with suction. AIP Conf. Proc. 2014, 1605, 428–433. [Google Scholar]
  27. Turkyilmazoglu, M. Flow of a micropolar fluid due to a porous stretching sheet and heat transfer. Int. J. Non-Linear Mech. 2016, 83, 59–64. [Google Scholar] [CrossRef]
  28. Mabood, F.; Ibrahim, S.M.; Rashidi, M.M.; Shadloo, M.S. Giulio Lorenzini, Nonuniform heat source/sink and Soret effects on MHD non-Darcian convective flow past a stretching sheet in a micropolar fluid with radiation. Int. J. Heat Mass Transf. 2016, 93, 674–682. [Google Scholar] [CrossRef]
  29. Ramzan, M.; Farooq, M.; Hayat, T.; Chung, J.D. Radiative and Joule heating effects in the MHD flow of a micropolar fluid with partial slip and convective boundary condition. J. Mol. Liq. 2016, 221, 394–400. [Google Scholar] [CrossRef]
  30. Hsiao, K. Micropolar nanofluid flow with MHD and viscous dissipation effects towards a stretching sheet with multimedia feature. Int. J. Heat Mass Transf. 2017, 112, 983–990. [Google Scholar] [CrossRef]
  31. Mahmood, A.; Chen, B.; Ghaffari, A. Hydromagnetic Hiemenz flow of micropolar fluid over a nonlinearly stretching/shrinking sheet: Dual solutions by using Chebyshev Spectral Newton Iterative Scheme. J. Magn. Magn. Mater. 2016, 416, 329–334. [Google Scholar] [CrossRef]
  32. Borrelli, A.; Giantesio, G.; Patria, M.C. MHD orthogonal stagnation-point flow of a micropolar fluid with the magnetic field parallel to the velocity at infinity. Appl. Math. Comput. 2015, 264, 44–60. [Google Scholar] [CrossRef]
  33. Abdel-Wahed, M.S. Lorentz force effect on mixed convection micropolar flow in a vertical conduit. Eur. Phys. J. Plus 2017, 132, 195. [Google Scholar] [CrossRef]
  34. Abdel-wahed, M.S. Flow and heat transfer of a weak concentration micropolar nanofluid over steady/unsteady-moving surface. Appl. Phys. A 2017, 123, 195. [Google Scholar] [CrossRef]
  35. Vo, T.Q.; Kim, B.H. Molecular Dynamics Study of Thermodynamic Properties of Nanoclusters for Additive Manufacturing. Int. J. Precis. Eng. Manuf. Green Technol. 2017, 4, 301–306. [Google Scholar] [CrossRef]
  36. Vo, T.Q.; Park, B.; Park, C.; Kim, B. Nano-scale liquid film sheared between strong wetting surfaces: Effects of interface region on the flow. J. Mech. Sci. Technol. 2015, 29, 1681–1688. [Google Scholar] [CrossRef]
  37. Priezjev, N.V. Molecular dynamics simulations of oscillatory Couette flows with slip boundary conditions. Microfluid Nanofluid 2013, 14, 225–233. [Google Scholar] [CrossRef]
  38. Kim, B.H.; Beskok, A.; Cagin, T. Viscous heating in nanoscale shear driven liquid flows. Microfluid Nanofluid 2010, 9, 31–40. [Google Scholar] [CrossRef]
  39. Rodrıguez-Lopez, J.L.; Montejano-Carrizales, J.M.; Jose-Yacaman, M. Molecular dynamics study of bimetallic nanoparticles: The case of AuxCuy alloy clusters. Appl. Surf. Sci. 2003, 219, 56–63. [Google Scholar] [CrossRef]
  40. Kim, B. Thermal resistance at a liquid–solid interface dependent on the ratio of thermal oscillation frequencies. Chem. Phys. Lett. 2012, 554, 77–81. [Google Scholar] [CrossRef]
  41. Murad, S.; Puri, I.K. Molecular simulation of thermal transport across hydrophilic interfaces. Chem. Phys. Lett. 2008, 467, 110–113. [Google Scholar] [CrossRef]
Figure 1. Physical model and coordinate system.
Figure 1. Physical model and coordinate system.
Applsci 08 00736 g001
Figure 2. Influence of the material parameter (K). (a) Velocity profile and (b) angular velocity profile.
Figure 2. Influence of the material parameter (K). (a) Velocity profile and (b) angular velocity profile.
Applsci 08 00736 g002
Figure 3. Influence of the unsteadiness parameter (A). (a) Velocity profile, (b) angular velocity profile, and (c) temperature profile.
Figure 3. Influence of the unsteadiness parameter (A). (a) Velocity profile, (b) angular velocity profile, and (c) temperature profile.
Applsci 08 00736 g003
Figure 4. Influence of the nanoparticle concentration. (a) Velocity profile, (b) angular velocity profile, and (c) temperature profile.
Figure 4. Influence of the nanoparticle concentration. (a) Velocity profile, (b) angular velocity profile, and (c) temperature profile.
Applsci 08 00736 g004
Figure 5. Influence of suction parameter (Sw). (a) Velocity profile, (b) angular velocity profile, and (c) temperature profile.
Figure 5. Influence of suction parameter (Sw). (a) Velocity profile, (b) angular velocity profile, and (c) temperature profile.
Applsci 08 00736 g005
Figure 6. Influence of the porous media parameter (kp). (a) Velocity profile and (b) Angular velocity profile.
Figure 6. Influence of the porous media parameter (kp). (a) Velocity profile and (b) Angular velocity profile.
Applsci 08 00736 g006
Figure 7. Influence of the magnetic parameter (M). (a) Velocity profile and (b) angular velocity profile.
Figure 7. Influence of the magnetic parameter (M). (a) Velocity profile and (b) angular velocity profile.
Applsci 08 00736 g007aApplsci 08 00736 g007b
Figure 8. Temperature profile with increasing radiation parameter (R).
Figure 8. Temperature profile with increasing radiation parameter (R).
Applsci 08 00736 g008
Table 1. Thermophysical properties of water and copper (Cu) [12].
Table 1. Thermophysical properties of water and copper (Cu) [12].
PropertiesFluid (Water)Cu
C p ( J / k g   K ) 4179385
D e n s i t y   ( ρ )   ( k g / m 3 ) 997.18933
T h e r m a l   c o n d u c t i v i t y   ( k )   ( W / m   K ) 0.613400
Table 2. Comparison of C f R x and g ( 0 ) for several values of the suction parameter ( S w ) at the unsteadiness parameter (A) = radiation parameter (R) = 0, n = 0.5, nanoparticle concentration ( ϕ ) = 0.1 and the material parameter (K) = 2, magnetic field parameter ( M ) = 0.1 ,   permeability   of   the   porous   medium   parameter   ( k p ) = 0.1 .
Table 2. Comparison of C f R x and g ( 0 ) for several values of the suction parameter ( S w ) at the unsteadiness parameter (A) = radiation parameter (R) = 0, n = 0.5, nanoparticle concentration ( ϕ ) = 0.1 and the material parameter (K) = 2, magnetic field parameter ( M ) = 0.1 ,   permeability   of   the   porous   medium   parameter   ( k p ) = 0.1 .
ParameterPresent ResultsEzzah [26]
Exact SolutionNumerical Solution
S w C f R e g ( 0 ) C f R e g ( 0 ) C f R e g ( 0 )
23.91821.91903.91821.91903.91831.9191
35.274123.47705.274123.47705.27413.4770
Table 3. Values of f ( 0 ) , g ( 0 ) ,   and   θ ( 0 ) for the steady case [A = 0] at R = 1, Prandtl number (Pr) = 6.2, M = 0.1 ,   k p = 0.1 .
Table 3. Values of f ( 0 ) , g ( 0 ) ,   and   θ ( 0 ) for the steady case [A = 0] at R = 1, Prandtl number (Pr) = 6.2, M = 0.1 ,   k p = 0.1 .
K S w ϕ f ( 0 ) g ( 0 ) θ ( 0 ) C f R e N u / R e M x
0−0.50.14.0383460.9794470.531071−8.0766920.531071−1.958894
0.13.5840310.7128130.518383−7.1898170.690297−1.699500
0.23.1307940.5762390.564942−4.2204490.986223−1.152478
00.04.4648031.7286920.924024−13.394400.924024−4.321730
0.13.9138611.4216751.020822−5.6721461.359363−2.843350
0.23.3780400.9902401.344222−5.8571062.346619−2.263903
0.50.04.9362953.6250161.828576−9.8725901.828576−7.250032
0.14.2740452.5932311.895725−8.5740342.524416−6.182824
0.23.6448101.7692112.522203−4.9133654.403029−3.538422
0.5−0.50.03.4055830.8039440.522100−10.216740.522100−2.009860
0.13.1241870.6866800.551381−4.5277140.734239−1.373360
0.22.8142920.5742420.594486−4.8796361.037798−1.312841
00.03.7434951.4505750.992598−7.4869900.992598−2.901150
0.13.4004681.1895471.054966−6.8215771.404831−2.836137
0.23.0301690.9267401.126033−4.0848021.965724−1.853480
0.50.04.1411002.6613021.888458−12.423301.888458−6.653255
0.13.7152032.0445831.944352−5.3842422.589170−4.089166
Table 4. Values of f ( 0 ) , g ( 0 ) ,   and   θ ( 0 ) for the unsteady case [A = 1] at R = 1, Pr = 6.2, M = 0.1 ,   k p = 0.1 .
Table 4. Values of f ( 0 ) , g ( 0 ) ,   and   θ ( 0 ) for the unsteady case [A = 1] at R = 1, Pr = 6.2, M = 0.1 ,   k p = 0.1 .
K S w ϕ f ( 0 ) g ( 0 ) θ ( 0 ) C f R e N u / R e M x
0−0.50.04.1905533.1020111.339903−8.3811061.339903−6.204022
0.13.7170992.4766671.365301−7.4567611.818084−5.904910
0.23.2453271.9154561.393946−4.3748442.433422−3.830912
00.04.6130034.2631781.777243−13.839001.777243−10.65794
0.14.0437743.2760031.810996−5.8604222.411588−6.552006
0.23.4901672.4403922.177493−6.0515203.801266−5.579264
0.50.05.0790285.9190602.410378−10.158052.410378−11.83812
0.14.3998174.3677752.450397−8.8263413.263038−10.41372
0.23.7538983.1282623.065966−5.0604205.352281−6.256524
0.5−0.50.03.62466382.4812601.366782−10.873991.366782−6.203150
0.13.3086562.0725411.386910−4.7950551.846860−4.145082
0.22.9654901.6707361.410602−5.1417952.462499−3.819664
00.03.9583213.3015331.810306−7.9166421.810306−6.603066
0.13.57987032.6795051.837633−7.1814712.447059−6.388520
0.23.1779862.0960671.869445−4.2840663.263504−4.192134
0.50.04.3307444.4344992.446977−12.992232.446977−11.08624
0.13.8782163.4898792.479820−5.6204873.302219−6.979758
0.23.4084462.64438082.517496−5.9098264.394812−6.045627

Share and Cite

MDPI and ACS Style

Raslan, K.; Mohamadain, S.; Abdel-wahed, M.; Abedel-aal, E.M. MHD Steady/Unsteady Porous Boundary Layer of Cu–Water Nanofluid with Micropolar Effect over a Permeable Surface. Appl. Sci. 2018, 8, 736. https://doi.org/10.3390/app8050736

AMA Style

Raslan K, Mohamadain S, Abdel-wahed M, Abedel-aal EM. MHD Steady/Unsteady Porous Boundary Layer of Cu–Water Nanofluid with Micropolar Effect over a Permeable Surface. Applied Sciences. 2018; 8(5):736. https://doi.org/10.3390/app8050736

Chicago/Turabian Style

Raslan, Kamal, Selim Mohamadain, Mohamed Abdel-wahed, and Elsayed M. Abedel-aal. 2018. "MHD Steady/Unsteady Porous Boundary Layer of Cu–Water Nanofluid with Micropolar Effect over a Permeable Surface" Applied Sciences 8, no. 5: 736. https://doi.org/10.3390/app8050736

APA Style

Raslan, K., Mohamadain, S., Abdel-wahed, M., & Abedel-aal, E. M. (2018). MHD Steady/Unsteady Porous Boundary Layer of Cu–Water Nanofluid with Micropolar Effect over a Permeable Surface. Applied Sciences, 8(5), 736. https://doi.org/10.3390/app8050736

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop