Next Article in Journal
The Effect of Dry Friction upon the Dynamics of a Short Eccentric Rotor: An Analytical and Experimental Study
Next Article in Special Issue
Surface Laser Micropatterning of Polyethylene (PE) to Increase the Shearing Strength of Adhesive Joints
Previous Article in Journal
Research on Internal Flow Field Characteristics of Straight-Groove Texture Using Three-Dimensional Modeling
Previous Article in Special Issue
Friction and Wear Characteristics of Bacterial Cellulose Modified by Microcellular Foaming Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Polymer Dispersion Effects on Drag, Heat Transfer, and Mass Transfer in Non-Newtonian Based Nanofluids

1
Department of Mathematics, COMSATS University Islamabad, Islamabad 44000, Pakistan
2
Department of Mathematics and Physics, University of Campania “L. Vanvitelli”, 81100 Caserta, Italy
3
Center for Applied Mathematics and Bioinformatics (CAMB), Gulf University for Science and Technology, Hawally 32093, Kuwait
*
Author to whom correspondence should be addressed.
Lubricants 2023, 11(8), 339; https://doi.org/10.3390/lubricants11080339
Submission received: 7 June 2023 / Revised: 21 July 2023 / Accepted: 6 August 2023 / Published: 9 August 2023
(This article belongs to the Special Issue Tribology of Polymer-Based Composites)

Abstract

:
This article explores the influence of polymers on the boundary layer flow, heat transfer, and mass transfer control of non-Newtonian-based nanofluids flowing past a stretching surface. The mathematical model incorporates the Oldroyd-B model to analyze the effects of polymers, while the Powell–Eyring and Reiner–Philippoff viscosity models are employed to study the behavior of non-Newtonian fluids. The dispersion model is adopted to account for nanofluid characteristics. Appropriate transformations yield governing equations with similar forms, which are solved numerically to investigate the impact of polymer inclusion on skin friction, Nusselt number, and Sherwood number. The study’s findings reveal that the addition of polymers to the non-Newtonian-based nanofluids leads to a reduction in heat and mass transport while enhancing skin drag. Detailed analysis of these effects sheds light on the underlying physical mechanisms.

Graphical Abstract

1. Introduction

In recent years, the study of rheology, or the science of the deformation and flow of matter, has become increasingly important in the field of polymers. As polymers are widely used in various scientific and technological fields, it is essential to understand their physical behavior and properties. At the same time, non-Newtonian fluids have also gained significant attention due to their occurrence in natural and industrial processes. The introduction of small amounts of polymers to non-Newtonian fluids can bring about a significant rheological change, leading to a wide range of remarkable phenomena. Nanofluids have shown immense potential in improving the heat transfer performance of fluids. However, the addition of nanoparticles often leads to an increase in viscosity, making the flow more challenging to control [1]. To address this issue, the use of polymers has been proposed as a means of reducing viscosity and controlling heat transfer [2]. This innovation has opened new avenues for research and development, leading to novel applications in the fields of thermal management, and energy storage.
In the realm of polymer fluid dynamics, two distinct approaches have been developed to account for the effects of polymer microstructure and concentration in order to describe the impact of viscoelasticity. The first approach is the molecular approach, which includes models such as the elastic dumbbell, rigid dumbbell, and chain. The second approach is the continuum approach, which includes models such as the convected Maxwell, corotational Maxwell, Oldroyd-8 constant, and many others [3]. Recent studies have explored the impact of polymers on fluid dynamics using both approaches. Celani et al. [4] used the Hookean dumbbell model to investigate the dynamics of polymers in a random flow, while Al-Yaari et al. [5] employed the same model to examine the impact of polymers as drag-reducing agents on horizontal oil–water flow. Another recent study focused on the turbulent Rayleigh–Benard (RB) convection of water in a cylindrical cell, where the addition of polymers resulted in reduced heat transport. This effect is believed to be due to the interplay between polymers and the flow of the boundary layer in turbulent RB convection [6]. Benzi et al. [7] utilized the Oldroyd-B model, a molecular approach, to investigate the impact of polymers on the heat transport of a boundary layer flow. Athar et al. [8] explored the effects of polymers on the flow of a Newtonian fluid over a heated and permeable stretching surface, utilizing the continuum approach. Marenduzzo et al. [9] used the same approach to study polymers with finite thickness, while Doufas et al. [10] investigated the effects of polymers in homogeneous flow fields under isothermal conditions. Zhang et al. [11] examined the use of superhydrophobic surfaces and their potential applications in drag and heat transfer, and Ahmad et al. [12] employed the FENE-P model to investigate the impact of magnetic effects on heat transfer and drag reduction.
Two distinct models have been proposed to study the characteristics of nanofluids. The first model, the two-phase model, was proposed by Tiwari and Das [13] and is based on the homogeneous principle. The second model, the dispersion model, was proposed by Buongiorno [14] and is based on the slip mechanism between nanoparticles and conducting fluid. This latter model takes into account the implications of Brownian and thermophoresis mechanisms, in addition to other slip mechanisms. Various studies have been conducted based on these models, with different features aimed at examining various characteristics of nanofluid flows (see, for example, [15,16,17,18,19,20,21] and the references cited therein). The nonlinear motions of axisymmetric ternary hybrid nanofluids in thermally radiated expanding or contracting permeable Darcy Walls with various forms and densities were studied by Raju et al. [22]. Kumar et al. [23] investigated the effects of hybrid nanoparticles and transpiration on linear and quadratic convection on 3D flow. Recently, Ahmed et. al. [24] studied the influence of polymers on drag reduction and heat transfer enhancement of a nanofluid past a uniformly heated permeable vertically stretching surface with a prime focus on analyzing the possible effects of polymer inclusion in the nanofluid on drag coefficient, Nusselt number, and Sherwood number. Reductions in the drag coefficient, Nusselt number, and Sherwood number are reported in this study because of polymer additives.
In this article, we analyze the flow of polymeric fluids in the presence of nanoparticles using the molecular approach. While there has been some research on these topics, no theoretical study has yet explored the interaction of polymeric fluids with non-Newtonian-based nanofluids. Therefore, we aim to contribute to this field by investigating the possible effects of polymer additives on the drag, heat, and mass transfer of the boundary layer flow of non-Newtonian-based nanofluids over a stretching sheet. To achieve our research objectives, we will use the Reiner–Philippoff and Powell–Eyring viscosity models, which are commonly employed in the study of non-Newtonian fluids. Our focus will be on understanding how the presence of polymers can affect the behavior of non-Newtonian-based nanofluids. Specifically, we aim to analyze how the addition of polymers can alter the drag force, as well as the heat and mass transfer characteristics of the boundary layer flow. Our study is primarily theoretical in nature and will be based on a careful analysis of the relevant mathematical models. We anticipate that our research will help to shed light on the fascinating ways in which polymers and nanoparticles interact in fluids, and will have important implications for a range of industrial applications. Ultimately, we hope that our work will serve as a springboard for future studies in this exciting area of research.

2. Mathematical Model

The aim of this study is to investigate the effects of polymers on the boundary layer flow and heat transport of a non-Newtonian-based nanofluid past a stretching sheet. The presence of polymers produces an extra stress term in the momentum equation of the flow. This stress tensor due to the presence of polymers, denoted by  i j , relies upon the amount of stretching of the polymers and is a characteristic of the dimensionless conformation tensor  β i j  of the polymers, which is derived from the ensemble average of the product of the polymer end-to-end distance vector. Using the Oldroyd-B model of polymers gives [7]:
i j = ν p κ ¯ ( β i j δ i j ) ,
where  κ ¯  is the largest relaxation time for polymers,  δ i j  is the Kronecker delta, and  ν p  is the kinematic viscosity due to the contribution of polymers in the solution.
We consider a horizontal surface that stretches with nonlinear velocity in the x-direction in a non-Newtonian-based nanofluid. The temperature  T  and the nanoparticle fraction  N  assume the constant values  T w  and  N w  at the surface (See Figure 1). Assuming that the nanoparticle concentration is dilute, the boundary layer equations for the conservation of mass, momentum, energy, and nanoparticle concentration in the presence of polymers are written as [7]:
u x + v y = 0 ,
u u x + v u y = 1 ρ f τ x y y + y ( ν p κ ¯ β x y ) ,
u T x + v T y = α ¯ 2 T y 2 + ϵ ( D B T y N y + D T T ( T y ) 2 ) ,
u N x + v N y = D B 2 N y 2 + D T T 2 T y 2 .
In Equations (2)–(5),  u ( x , y )  and  v ( x , y )  are the velocity components,  ρ f  is the density of the fluid,  α ¯  is the thermal diffusivity,  ϵ  is the volumetric volume expansion coefficient of the nanofluid, and  D T  and  D B  are the thermophoretic diffusion and Brownian diffusion coefficient, respectively. Also,  τ x y  is the shear stress that is correlated to the rate of strain nonlinearly in unique forms for different non-Newtonian fluids.
The appropriate boundary conditions of the problem are as follows [25]:
u ( x , y ) = U w , v ( x , y ) = V w , T ( x , y ) = T w , N ( x , y ) = N w at y = 0 , u ( x , y ) 0 , T T , N N as y .
Larson et al. [26] found out that the friction coefficient of the fluid depends on the shape of the polymer. The value of this coefficient is small when the polymers are in the equilibrium state and large when the polymers are stretched. Hence the relaxation time increases with the stretching of the polymers.  κ ¯  is a linear function of the end-to-end distance vector of polymers,  β ( β x x + β y y ) 1 / 2 .
κ ¯ = κ ¯ 0 ( 1 + a β 1 + a ) ,
where  κ ¯ 0  is the relaxation time of the polymers in the equilibrium shape and  a  is a positive parameter. In the equilibrium state,  β = 1  and  κ ¯  is equal to  κ ¯ 0  and when there is stretching,  β > 1  and  κ ¯ > κ ¯ 0 .
Understanding the behavior of a polymeric chain affected by a fluid flow is the first step in the study of polymer transport. When polymers are placed in a flow, they deform and stretch. Two effects determine the deformation of polymers; their stretching due to the velocity gradients and the elastic relaxation towards the equilibrium shape [4]. In the presence of polymers, even when they are not stretched, the viscosity of the solution increases. Now, the kinematic viscosity of the polymer solution without the stretching can be written as,  ν ¯ = ν + ν p . To proceed further, we transform our equations into ordinary differential equations using appropriate similarity transformations. For the stretching velocity of the form;  U ( x ) = U 0 x 1 3 , we have [25]:
ψ = U ( x ) x ν ¯ f ( η ) , τ = ρ U 0 3 ν ¯ g ( η ) , θ = T T T w T , ϕ = N N N w N ,
and the independent variable
η = U ( x ) ν ¯ x y ,
Equation (3) has unique forms for different non-Newtonian fluids models. Here, we will be using two models, Reiner–Philippoff and Powell–Eyring, to study the effects of polymers on the flow and heat transfer of the nanofluid and concentration of nanoparticles by using a non-Newtonian fluid as a base fluid.

2.1. Reiner–Philippoff Fluid

The stress deformation of a Reiner–Philippoff fluid is expressed as [25]:
u y = τ x y μ + μ 0 μ 1 + ( τ x y τ s ) 2 ,
where  τ s  is the reference shear stress,  μ 0  and  μ  is the zero and upper Newtonian limiting viscosity, respectively. The Reiner–Philippoff fluid model is one of the few models that exhibits all the pseudoplastic, Newtonian, and dilatant behaviors.
The transformed Equations (3)–(5) and (8) can be written as:
g + γ W i R e η ( 1 + a 1 + a β β x y ) = 1 3 f 2 2 3 f f ,
g = f λ γ 1 2 + ( 1 γ ) g 2 γ 1 2 + g 2 ,
θ + 2 3 P r f θ + N b θ ϕ + N t θ 2 = 0 ,
ϕ + 2 3 L e f ϕ + N t N b θ = 0 ,
where  W i  is the Weissenberg number, which is the ratio of elastic forces to the viscous forces, and  R e  is the Reynolds number, the ratio of inertial forces to the viscous forces. By substituting  x = L  such that  L  is the plate’s total length, we have:
W i τ 0 U 0 L 2 / 3 , R e U 0 L 4 / 3 ν ¯ ,
and  γ ν p ν + ν p  is a function of the concentration of the polymers and  λ  and  γ 1  are the Reiner–Philippoff model parameters,  P r  is the Prandtl number,  L e  is the Lewis number,  N b  is the Brownian motion parameter, and  N t  is the thermophoretic diffusion parameter, which are expressed as [25]:
γ 1 = τ s ρ U 0 3 ν ¯ , λ = μ 0 μ , N b = ϵ D B Δ N α ¯ , P r = ν ¯ α ¯ , L e = ν ¯ D B , N t = τ D T Δ T α ¯ T
The boundary conditions transform to:
f ( 0 ) = s , f ( 0 ) = 1 , θ ( 0 ) = 1 , ϕ ( 0 ) = 1 , f 0 , θ 0 , s 0 a s η .
For  λ = 1 γ , the above problem reduces for the flow of Newtonian fluid and when  γ = 0 , we get the flow of a Reiner–Philippoff fluid without polymers.

2.2. Powell–Eyring Fluid

For a Powell–Eyring fluid, we have the following stress-strain relationship [27]:
τ x y = μ u y + 1 ξ sinh 1 ( 1 C u y ) ,
where  μ  is the dynamic viscosity, and  ξ  and  C  are material constants. Considering the second-order approximation of the above function [27]:
sinh 1 ( 1 C u y ) 1 C u y 1 6 ( 1 C u y ) 3 ,
where  | 1 C u y | < < 1 .
The momentum equation for a Powell–Eyring fluid in the presence of polymers becomes,
u u x + v u y = ν 2 u y 2 + 1 ξ C ρ f 2 u y 2 1 2 ξ C 3 ρ f ( u y ) 2 2 u y 2 + y ( ν p κ ¯ β x y ) .
In the case of the Powell–Eyring model, we apply the same transformations as above, and the transformed Equation (17) attains the form:
1 3 f 2 2 3 f f = ( 1 γ ) f + K f K A f f 2 + γ W i R e η ( 1 + a 1 + a β β x y ) ,
where  K  and  A  are the fluid parameters and can be expressed as:
K = 1 ξ C μ , A = U 0 3 2 ν ¯ C 2
Equations (4) and (5) and the boundary conditions remain the same for the model.

2.3. Polymers End-to-End Distance Vector

To proceed further, we must provide the momentum equation with specific information on  β x y . Neglecting the thermal noise, we obtain the following equations for the components of a dimensionless polymer’s end-to-end distance vector [7]:
d q x d t = 1 2 κ ¯ ( q x q 0 x ) + q x u x + q y u y ,
d q y d t = 1 2 κ ¯ ( q y q 0 y ) + q x v x + q y v y ,
where  α  is the angle of equilibrium orientation and  q 0 x = cos α  and  q 0 y = sin α . We assume that the polymers follow streamlines of the flow so that:
d η d t = u η x + v η y .
Now, Equations (19) and (20) take the form:
f d q x d η = 3 4 ( 1 + a ) ( q x q 0 x ) W i ( 1 + a β ) 1 2 ( f η f ) q x 3 2 R e f q y ,
f d q y d η = 3 4 ( 1 + a ) ( q y q 0 y ) W i ( 1 + a β ) 1 6 R e ( 2 f η 2 f ) q x 1 2 ( f η f ) q y .
The quantity  β x y  is the average of  q x q y  over the polymers in a small region. This ensemble average is the same as an average over  α , which is uniformly distributed in  [ 0 , 2 π ] , thus  β x y = q x q y α . It is observed that the equation of motion remains unchanged on reflection and inversion, thus averaging  α  upon  [ 0 , 2 π ]  is the same as averaging upon  [ 0 , π / 2 ] . If  q ^ x  and  q ^ y  are values of  q x  and  q y  at  α = π / 4 , Equations (22) and (23) take the form:
f d q ^ x d η = 3 4 ( 1 + a ) ( q ^ x 1 / 2 ) W i ( 1 + a β ) 1 2 ( f η f ) q ^ x 3 2 R e f q ^ y ,
f d q ^ y d η = 3 4 ( 1 + a ) ( q ^ y 1 / 2 ) W i ( 1 + a β ) 1 6 R e ( 2 f η 2 f ) q ^ x 1 2 ( f η f ) q ^ y .
And thus we have:
β x y W i R e ( 1 + a β 1 + a ) f + ( q ^ x 1 / 2 ) ( q ^ y 1 / 2 ) = W i R e 1 + a β 1 + a f ( 1 + h ) ,
where
h ( η ) ( q ^ x ( η ) 1 / 2 ) ( q ^ y ( η ) 1 / 2 ) W i R e f ( η ) ( 1 + a β ( η ) ) / ( 1 + a )
with  β ( η ) = q ^ x 2 ( η ) + q ^ y 2 ( η ) .  Moreover, we have the viscosity as
ν ν p ( 1 + h ) ,
and the polymeric solution’s total viscosity is:
ν t o t = ν + ν p ν ¯ + ν p h ( η ) .
Substituting the expression of  β x y  from Equation (26) in (9) and (18), we have for the Reiner–Philippoff model:
g = 1 3 f 2 2 3 f f γ d d η ( f ( 1 + h ) ) ,
f = g ( γ 1 2 + g 2 λ γ 1 2 + g 2 ( 1 γ ) ) ,
and for the Powell–Eyring model:
1 3 f 2 2 3 f f = f ( 1 + K K A f 2 ) + γ d d η ( f h ) .
Solving (24), (25), and (27)–(29), and then with (30) gives  f ( η )  for the polymeric flow for Reiner–Philippoff and Powell–Eyring fluids, respectively.

3. Solution Methodology

In this study, the problem is solved using a numerical method iteratively. To obtain the solution to the problem without the presence of polymers, a nonlinear shooting method with a subroutine of the Runge–Kutta scheme is used, and the results are then utilized to solve Equations (24) and (25) for  q ^ x  and  q ^ y . The expression for  h ( η )  is obtained from (27), and this is used to solve (28), (29), and then with (30), again using the shooting method. The solution of  f ( η ) , the velocity profile for polymeric flow, attained from (28), (29), and (30), is then used in (24) and (25) to obtain a modified  h ( η ) . The convergence of the iterative scheme is depicted in Figure 2, which shows that convergence is achieved after just a few iterations. Finally, by using the solution of the momentum equation, the temperature and concentration profiles of the problem are obtained. To clarify the numerical method used to solve the set of equations for the given model, Figure 3 is presented, which shows the step-by-step approach.

4. Results and Discussions

This section investigates how polymers affect the drag coefficient, Nusselt number, and Sherwood number of non-Newtonian nanofluids. We present graphs for both the Reiner–Philippoff and Powell–Eyring models.

4.1. Impact of Polymers on Drag Coefficient

The skin friction for non-Newtonian fluids is influenced by various factors, including the type of fluid, its rheological properties, and the presence of polymers. Understanding these factors is crucial for predicting and optimizing the flow behavior of non-Newtonian fluids in various industrial applications. Skin friction is the shear stress exerted on the surface by the flowing viscous fluid and it can be mathematically expressed as:
C f = τ w ρ U w 2 ,
where  τ w  is shear stress at the surface. For a Reiner–Philippoff fluid,  τ w  can be written as [26]:
τ w = ( μ + μ 0 μ 1 + ( τ w τ s ) 2 ) u y | y = 0 ,
and for a Powell–Eyring fluid [27]:
τ w = μ u y + 1 β ( 1 C u y 1 6 ( 1 C u y ) 3 ) | y = 0
Using Expressions (31) and (32) in (30), skin friction for a Reiner–Philippoff fluid can be expressed as:
C f R e 1 / 2 = ( ( 1 γ ) + λ 1 + γ 1 + ( g ( 0 ) ) 2 γ 1 2 ) f ( 0 ) ,
and for a Powell–Eyring fluid, we have:
C f R e 1 / 2 = ( 1 γ + K ) f ( 0 ) K A 3 ( f ( 0 ) ) 3 .
Figure 4 depicts the behavior of skin friction ( S f ) of the Reiner–Philippoff-based nanofluid in the presence of polymers. Skin friction is plotted against the non-Newtonian fluid parameter,  λ ,  for different values of the polymer concentration parameter. The value  γ = 0  corresponds to the flow without the presence of polymers. In the present scenario, the Reiner–Phippoff model acts as dilatant for  λ < 1 γ , pseudoplastic for  λ > 1 γ ,  and becomes Newtonian for  λ = 1 γ . It is observed that skin friction for a dilatant fluid is higher than for a viscous fluid and a viscous fluid bears more skin friction than a pseudoplastic fluid. Also, one can note that skin friction increases with increasing  γ . This implies that the interaction between the polymer molecules and the sliding surface produces higher resistance, which increases skin drag. Figure 5 investigates the skin friction for a Powell–Eyring fluid in the presence of polymers. Skin friction is plotted against the non-Newtonian fluid parameter,  K , for different values of the polymer concentration parameter. The same behavior is seen as in the case of a Reiner–Philippoff fluid. As we increase the value of the fluid parameter, the skin friction decreases, and increasing the concentration of the polymers increases the skin friction of the fluid flow.
Figure 6 shows the influence of the presence of polymers on the percentage of drag enhancement ( % D E ) of the non-Newtonian-based nanofluid. Drag enhancement is defined as the difference between the skin friction with polymers and without polymers (when  γ = 0 ) present. It can be expressed as  D E = S f S f 0 S f 0  being the value without the presence of polymers. In Figure 6 % D E  is plotted against the polymer concentration parameter,  γ , for varying non-Newtonian parameter,  λ , and Weissenberg number. The Weissenberg number is a dimensionless number that compares the relaxation time to the characteristic stretching time. It is observed that the greater the concentration of polymers, the greater the resistance between the polymer molecules and the surface, which will increase the rate of enhancement in the skin drag; conversely, this increment decreases with increasing Reiner–Philippoff fluid parameter. Also, increasing  W i  implies that the relaxation time increases and the elastic force minimizes, allowing the velocity gradients to stretch the polymers. The fall in  % D E  with increasing Weissenberg number shows that increasing the relaxation time reduces the enhancement in the drag. The same effects on drag enhancement of varying  γ  and  W i  are observed for a Powell–Eyring fluid.

4.2. Impact of Polymers on Nusselt Number

Another important quantity of interest is the Nusselt number ( N u ) expressed as [26]:
N u = x q w k ( T w T )
where  q w  is the heat flux defined as:
q w = k T y | y = 0
Using Equation (36) in (35), we get:
N u r = R e 1 / 2 N u = θ ( 0 ) ,
which is known as the reduced Nusselt number.
In Figure 7, the behavior of the Nusselt number ( N u r ) of a Reiner–Philippoff-based nanofluid in the presence of polymers is discussed.  N u r  is plotted against the thermophoresis parameter,  N t ,  for nanoparticles with different values of the polymer concentration parameter. The zero value of  γ  corresponds to the flow without the presence of polymers. It is observed that increasing the value of  N t  reduces the  N u r . This implies that as the thermophoresis parameter increases, it will slow down the process of heat transfer. Also, the addition of polymers in the flow reduces the Nusselt number or the heat flux at the surface. Figure 8 shows  N u r  plotted against the Brownian diffusion parameter,  N b , for a Powell–Erying fluid. As can be seen, as the rate of diffusion between the particles increases, the heat flux or heat transfer rate decreases. Furthermore, it can be seen that this rate is also affected by increasing polymer concentration in the solution. The diffusion of the polymers in the solution also contributes to decreasing the process of heat transfer.
Figure 9 represents the effects of the presence of polymers on the percentage of heat reduction ( % H R ) of the non-Newtonian-based nanofluid flow. Heat reduction is defined as  H R = N u r N u r 0  ( N u r 0  being the value without the presence of polymers). In Figure 9 % H R  is plotted against the polymer concentration parameter for different values of Weissenberg number and non-Newtonian fluid parameter. This plot shows that if we increase the concentration of polymers in the solution, there is more heat reduction. Similarly, the reduction in  % H R  with increasing Weissenberg number shows that polymer stretching impedes the heat reduction process of the flow at the surface. The same results for heat reduction with varying  γ  and  W i  are observed for a Powell–Eyring fluid. Also, the heat reduction rate increases with the increasing value of the Reiner–Philippoff fluid parameter. Figure 10 indicates the influence of the presence of polymers on  % H R  as a function of  γ  for different  N b  for a Reiner–Philippoff fluid. It has already been observed that increasing the concentration of polymers in the solution increases heat reduction. One interesting thing to note is that the presence of polymers also affects the concentration of the nanoparticles in the flow. We see that as we increase the value of the Brownian diffusion parameter,  N b , heat reduction at the surface decreases. This implies that as the random movement of the particles increases due to collisions of the different molecules, the heat reduction process is impeded. Figure 11 shows the influence of polymers on  % H R  for different  N t  for a Powell–Eyring fluid. We find that when the value of the thermophoresis parameter increases, the heat transfer reduction at the surface decreases.

4.3. Impact of Polymers on Sherwood Number

Similarly, the Sherwood number can be defined as:
S h R e 1 / 2 = S h r = ϕ ( 0 ) ,
where  S h r  is the reduced Sherwood number. In Figure 12, the behavior of the Sherwood number ( S h r ) of a Reiner–Philippoff-based nanofluid in the presence of polymers is shown.  S h r  is plotted against the Brownian diffusion parameter,  N b , due to nanoparticles for different values of the polymer concentration parameter. A zero value of  γ  corresponds to the flow without the presence of polymers. It can be seen that increasing the value of  N b  increases the  S h r . This implies that as the Brownian diffusion parameter increases, convective mass transfer becomes dominant and the Sherwood number increases. Also, the addition of polymers in the flow reduces the Sherwood number or the mass flux at the surface. The same results are seen with the Powell–Eyring model.
Figure 13 presents the effect of the presence of polymers on the percentage of mass transfer reduction ( % M R ) of the flow, which is defined as  M R = S h r S h r 0  ( S h r 0  being the value without the presence of polymers). In Figure 13 % M R  is plotted against the polymer concentration parameter for different  W i  and  λ ,  which shows that if we increase the concentration of polymers in the solution, mass transfer at the surface decreases. Similarly, a decrease in  % M R  with increasing Weissenberg number shows that polymer stretching also impedes the mass reduction process of the flow at the surface. An increasing non-Newtonian flow parameter also causes a decrease in mass flux in the presence of polymers. In Figure 14, one interesting thing to note is that the polymers also affect the concentration of nanoparticles in the flow. We see that as we increase the value of the Brownian diffusion parameter,  N b , the percentage of mass reduction at the surface increases. This implies that as the random movement of the particles increases due to collisions of the different molecules, mass transfer due to convection becomes dominant and so the reduction process speeds up. Figure 15 shows the influence of the presence of polymers on  % M R  for different  N t  for a Powell–Eyring fluid. We find that when the value of the thermophoresis parameter increases, the heat transfer reduction at the surface decreases.
All the results presented in this article are in agreement with [7], in which Benzi et al. discussed the heat transport with polymers, and with [26], in which the authors studied the boundary layer flow of a non-Newtonian-based nanofluid.

5. Conclusions

In summary, this theoretical study investigated the influence of the presence of polymers on the boundary layer flow of non-Newtonian-based nanofluids. Various non-Newtonian fluid models and a molecular approach were employed to derive the governing equations. These equations were then solved numerically, and the results were analyzed. The key findings are as follows:
  • The addition of polymers to a non-Newtonian fluid generally increases the viscosity of the solution, although the viscosity reduction occurs to a lesser extent with polymer stretching.
  • The introduction of polymers in non-Newtonian-based nanofluids leads to an increase in the drag coefficient and a decrease in the Nusselt number and Sherwood number.
  • The concentration of polymers in the solution has a direct impact on skin friction, showing that an increase in polymer concentration results in higher skin friction.
  • Higher polymer concentrations correspond to greater drag enhancement, as well as increased heat and mass reduction.
  • The effects of thermophoresis and Brownian diffusion are not limited to reduced Nusselt number and Sherwood number; they also influence the behavior of the polymers.
This study primarily focused on the measurement of skin drag, heat transfer, and mass transfer in non-Newtonian-based nanofluids with polymer additives. Subjects that were not explored in this study present several avenues for further research, including investigating viscous dissipation effects and entropy generation and conducting experimental studies to support or challenge our findings.

Author Contributions

Software and writing—original draft preparation, A.S.; Conceptualization, and methodology, A.A.; validation and formal analysis, R.K.; writing—review and editing, R.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

  ν p viscosity due to polymers  ( m 2 s 1 )  
  κ ¯ relaxation time  ( s 1 )  
  μ dynamic viscosity  ( kg   m 1 s 1 )  
  W i Weissenberg number.
  γ polymer concentration parameter
  D T Thermophoretic diffusion coefficient  ( m 2 s 1 )  
  P r Prandtl number
  U w stretching velocity  ( ms 1 )  
  f ( η ) dimensionless stream function
  κ ¯ 0 κ ¯ 0  relaxation time in equilibrium  ( s 1 )  
  μ 0 ,     μ zero and upper Newtonian limiting viscosity
  T ,   N ambient temperature and concentration
  N u r reduced Nusselt number
  ξ material fluid constants  ( ms 2 kg 1 )  
  ρ density ( kgm 3 )  
  α ¯ thermal diffusivity  ( m 2 s 1 ) .
  λ ,   γ 1 Reiner–Philippoff fluid parameter
  R e Reynolds number
  D B Brownian diffusion coefficient ( ms 1 )  
  N t thermophoresis diffusion coefficient
  N b Brownian diffusion coefficient.
  T dimensional fluid temperature  ( K )  
  K ,   A Powell–Eyring fluid parameter
  L e Lewis number
  T w ,     N w temperature and concentration at the surface
  S f     skin friction coefficient
  u ,   v dimensional velocity components  ( ms 1 )  
  C material fluid constants  ( s 1 )  

References

  1. Apmann, K.; Fulmer, R.; Soto, A.; Vafaei, S. Thermal Conductivity and Viscosity: Review and Optimization of Effects of Nanoparticles. Materials 2021, 14, 1291. [Google Scholar] [CrossRef] [PubMed]
  2. Mackay, M.E.; Dao, T.T.; Tuteja, A.; Ho, D.L.; van Horn, B.; Kim, H.C.; Hawker, C.J. Nanoscale effects leading to non-Einstein-like decrease in viscosity. Nat. Mater. 2003, 2, 762–766. [Google Scholar] [CrossRef] [PubMed]
  3. Bird, R.B.; Giacomin, A.J. Polymer fluid dynamics: Continuum and molecular approaches. Annu. Rev. Chem. Biomol. Eng. 2016, 7, 479–507. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Celani, A.; Musacchio, S.; Vincenzi, D. Polymer transport in random flow. J.Stat. Phys. 2005, 118, 531–554. [Google Scholar]
  5. Al-Yaari, M.; Soleimani, A.; Abu-Sharkh, B.; Al-Mubaiyedh, U.; Al-Sarkhi, A. Effect of drag reducing polymers on oil–water flow in a horizontal pipe. Int. J. Multiph. Flow. 2009, 35, 516–524. [Google Scholar] [CrossRef]
  6. Ahlers, G.; Nikolaenko, A. Effect of a polymer additive on heat transport in turbulent Rayleigh-Bénard convection. Phys. Rev. Lett. 2010, 104, 034503. [Google Scholar]
  7. Benzi, R.; Ching, E.S.; Chu, V.W. Heat transport by laminar boundary layer flow with polymers. J. Fluid Mech. 2012, 696, 330–344. [Google Scholar] [CrossRef] [Green Version]
  8. Athar, M.; Ahmad, A. Behavior of fluid flow and heat transfer induced by a stretching surface in the presence of polymers. Phys. Scr. 2021, 96, 13. [Google Scholar] [CrossRef]
  9. Marenduzzo, D.; Micheletti, C.; Seyed-Allaei, H.; Trovato, A.; Maritan, A. Continuum model for polymers with finite thickness. J. Phys. A Math. Gen. 2005, 38, L277. [Google Scholar] [CrossRef]
  10. Doufas, A.K.; Dairanieh, I.S.; McHugh, A.J. A continuum model for flow-induced crystallization of polymer melts. J. Rheol. 1999, 43, 85–109. [Google Scholar] [CrossRef]
  11. Zhang, Y.; Zhang, Z.; Yang, J.; Yue, Y.; Zhang, H. A Review of Recent Advances in Superhydrophobic Surfaces and Their Applications in Drag Reduction and Heat Transfer. Nanomaterials 2022, 12, 44. [Google Scholar]
  12. Ahmad, A.; Ishaq, A.; Khan, Y. Influence of FENE-P fluid on drag reduction and heat transfer past a magnetized surface. Int. J. Mod. Phys. B 2022, 36, 2250145. [Google Scholar] [CrossRef]
  13. Tiwari, R.K.; Das, M.K. Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. Int. J. Heat Mass Transf. 2007, 50, 2002–2018. [Google Scholar]
  14. Buongiorno, J. Convective transport in nanofluids. J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  15. Eastman, J.A.; Choi SU, S.; Li, S.; Yu, W.; Thompson, L.J. Anomalously increased effective thermal conductivities of ethylene glycol-based nanofluids containing copper nanoparticles. Appl. Phys. Lett. 2001, 78, 718–720. [Google Scholar] [CrossRef]
  16. Sheikholeslami, M.; Rashidi, M.M.; Hayat, T.; Ganji, D.D. Free convection of magnetic nanofluid considering MFD viscosity effect. J. Mol. Liq. 2016, 219, 79–85. [Google Scholar] [CrossRef]
  17. Abbas, W.; Sayed, E.A. Hall current and joule heating effects on free convection flow of a nanofluid over a vertical cone in presence of thermal radiation. Therm. Sci. 2017, 21, 863–876. [Google Scholar] [CrossRef] [Green Version]
  18. Mahanthesh, B.; Gireesha, B.J.; Animasaun, I.L.; Muhammad, T.; Shashikumar, N.S. MHD flow of SWCNT and MWCNT nanoliquids past a rotating stretchable disk with thermal and exponential space dependent heat source. Phys. Scr. 2019, 94, 085214. [Google Scholar] [CrossRef]
  19. Shehzad, N.; Zeeshan, A.; Ellahi, R.; Vafai, K. Convective heat transfer of nanofluid in a wavy channel: Buongiorno’s mathematical model. J. Mol. Liq. 2016, 219, 970–977. [Google Scholar] [CrossRef]
  20. Zheng, L.; Zhang, C.; Zhang, X.; Zhang, J. Flow and radiation heat transfer of a nanofluid over a stretching sheet with velocity slip and temperature jump in porous medium. J. Frankl. Inst. 2013, 350, 2142–2156. [Google Scholar] [CrossRef]
  21. Jahan, S.; Sakidin, H.; Nazar, R.; Pop, I. Analysis of heat transfer in nanofluid past a convectively heated permeable stretching/shrinking sheet with regression and stability analyses. Results Phys. 2018, 10, 195–202. [Google Scholar] [CrossRef]
  22. Raju CS, K.; Ahammad, N.A.; Sajjan, K.; Shah, N.A.; Yook, S.J.; Kumar, M.D. Nonlinear movements of axisymmetric ternary hybrid nanofluids in a thermally radiated expanding or contracting permeable Darcy Walls with different shapes and densities: Simple linear regression. Int. Commun. Heat Mass Transf. 2022, 135, 106110. [Google Scholar]
  23. Kumar, M.D.; Raju CS, K.; Sajjan, K.; El-Zahar, E.R.; Shah, N.A. Linear and quadratic convection on 3D flow with transpiration and hybrid nanoparticles. Int. Commun. Heat Mass Transf. 2022, 134, 105995. [Google Scholar]
  24. Ahmad, A.; Athar, M.; Khan, Y. Influence of polymers on drag and heat transfer of nanofluid past stretching surface: A molecular approach. J. Cent. South Univ. 2022, 29, 3912–3924. [Google Scholar] [CrossRef]
  25. Ahmad, A. Flow of ReinerPhilippoff based nano-fluid past a stretching sheet. J. Mol. Liq. 2016, 219, 643–646. [Google Scholar]
  26. Larson, R.G.; Perkins, T.T.; Smith, D.E.; Chu, S. Hydrodynamics of a DNA molecule in a flow field. Phys. Rev. 1997, 55, 1794. [Google Scholar]
  27. Jalil, M.; Asghar, S.; Imran, S.M. Self similar solutions for the flow and heat transfer of Powell-Eyring fluid over a moving surface in a parallel free stream. Int. J. Heat Mass Transf. 2013, 65, 73–79. [Google Scholar]
Figure 1. Geometry of the problem.
Figure 1. Geometry of the problem.
Lubricants 11 00339 g001
Figure 2. Result of  h ( η )  after the first (dashed-dot), second (dashed), and sixth (solid) iteration with  R e = 4900 a = 0.01 L e = 5.0 N b = N t = 0.1 λ = 0.8 γ = 0.2 , and  W i = 2.8 .
Figure 2. Result of  h ( η )  after the first (dashed-dot), second (dashed), and sixth (solid) iteration with  R e = 4900 a = 0.01 L e = 5.0 N b = N t = 0.1 λ = 0.8 γ = 0.2 , and  W i = 2.8 .
Lubricants 11 00339 g002
Figure 3. Interpretation of the numerical solution of the model.
Figure 3. Interpretation of the numerical solution of the model.
Lubricants 11 00339 g003
Figure 4. S f  of a Reiner–Philippoff fluid as a function of the fluid parameter,  λ , for different values of  γ  at fixed  γ 1 = 1.0 W i = 2.8 a = 0.01 P r = 5.0 R e = 4900 N b = N t = 0.1 ,  and  L e = 5.0 .
Figure 4. S f  of a Reiner–Philippoff fluid as a function of the fluid parameter,  λ , for different values of  γ  at fixed  γ 1 = 1.0 W i = 2.8 a = 0.01 P r = 5.0 R e = 4900 N b = N t = 0.1 ,  and  L e = 5.0 .
Lubricants 11 00339 g004
Figure 5. S f  of a Powell–Eyring fluid as a function of the fluid parameter,  K ,  for different values of  γ  at fixed  A = 0.3 W i = 2.8 a = 0.01 P r = 5.0 R e = 4900 N b = N t = 0.1 ,  and  L e = 5.0 .
Figure 5. S f  of a Powell–Eyring fluid as a function of the fluid parameter,  K ,  for different values of  γ  at fixed  A = 0.3 W i = 2.8 a = 0.01 P r = 5.0 R e = 4900 N b = N t = 0.1 ,  and  L e = 5.0 .
Lubricants 11 00339 g005
Figure 6. Effects on the percentage of drag enhancement ( % D E ) due to the presence of polymers as a function of  γ  for fixed  γ 1 = 1.0 a = 0.01 P r = 5.0 R e = 4900 N b = 0.1 N t = 0.5 , and  L e = 5.0  for a Reiner–Philippoff fluid.
Figure 6. Effects on the percentage of drag enhancement ( % D E ) due to the presence of polymers as a function of  γ  for fixed  γ 1 = 1.0 a = 0.01 P r = 5.0 R e = 4900 N b = 0.1 N t = 0.5 , and  L e = 5.0  for a Reiner–Philippoff fluid.
Lubricants 11 00339 g006
Figure 7. N u r  as a function of  N t  for different  γ  with  γ 1 = 1.0 W i = 2.8 a = 0.01 P r = 5.0 R e = 4900 N b = 0.1 L e = 5.0 ,  and  λ = 0.5  of a Reiner–Philippoff fluid.
Figure 7. N u r  as a function of  N t  for different  γ  with  γ 1 = 1.0 W i = 2.8 a = 0.01 P r = 5.0 R e = 4900 N b = 0.1 L e = 5.0 ,  and  λ = 0.5  of a Reiner–Philippoff fluid.
Lubricants 11 00339 g007
Figure 8. N u r  as a function of  N b  for different  γ  with  A = 0.3 W i = 2.8 a = 0.01 P r = 5.0 R e = 4900 N t = 0.1 L e = 5.0 ,  and  K = 0.1  of a Powell–Eyring fluid.
Figure 8. N u r  as a function of  N b  for different  γ  with  A = 0.3 W i = 2.8 a = 0.01 P r = 5.0 R e = 4900 N t = 0.1 L e = 5.0 ,  and  K = 0.1  of a Powell–Eyring fluid.
Lubricants 11 00339 g008
Figure 9. Effects on the percentage of heat reduction ( % H R ) due to polymers as a function of  γ  for fixed  γ 1 = 1.0 a = 0.01 P r = 5.0 R e = 4900 N b = 0.1 N t = 0.5 , and  L e = 5.0  for a Reiner–Philippoff fluid.
Figure 9. Effects on the percentage of heat reduction ( % H R ) due to polymers as a function of  γ  for fixed  γ 1 = 1.0 a = 0.01 P r = 5.0 R e = 4900 N b = 0.1 N t = 0.5 , and  L e = 5.0  for a Reiner–Philippoff fluid.
Lubricants 11 00339 g009
Figure 10. Effects on the percentage of heat reduction ( % H R ) due to the presence of polymers as a function of  γ  for different  N b  with  γ 1 = 1.0 λ = 0.7 a = 0.01 P r = 5.0 R e = 4900 N t = 0.1 ,  and  L e = 5.0  for a Reiner–Philippoff fluid.
Figure 10. Effects on the percentage of heat reduction ( % H R ) due to the presence of polymers as a function of  γ  for different  N b  with  γ 1 = 1.0 λ = 0.7 a = 0.01 P r = 5.0 R e = 4900 N t = 0.1 ,  and  L e = 5.0  for a Reiner–Philippoff fluid.
Lubricants 11 00339 g010
Figure 11. Effects on the percentage of heat reduction ( % H R ) due to the presence of polymers as a function of  γ  for different  N t  with  K = 0.3 A = 0.3 a = 0.01 P r = 5.0 R e = 4900 N b = 0.1 ,  and  L e = 5.0  for a Powell–Eyring fluid.
Figure 11. Effects on the percentage of heat reduction ( % H R ) due to the presence of polymers as a function of  γ  for different  N t  with  K = 0.3 A = 0.3 a = 0.01 P r = 5.0 R e = 4900 N b = 0.1 ,  and  L e = 5.0  for a Powell–Eyring fluid.
Lubricants 11 00339 g011
Figure 12. S h r  as a function of  N b  for different  γ  with  γ 1 = 1.0 W i = 2.8 a = 0.01 P r = 7.0 R e = 4900 N t = 0.1 L e = 5.0 ,  and  λ = 0.5  for a Reiner–Philippoff fluid.
Figure 12. S h r  as a function of  N b  for different  γ  with  γ 1 = 1.0 W i = 2.8 a = 0.01 P r = 7.0 R e = 4900 N t = 0.1 L e = 5.0 ,  and  λ = 0.5  for a Reiner–Philippoff fluid.
Lubricants 11 00339 g012
Figure 13. Effect on  % M R  of the presence of polymers as a function of  γ  for fixed  γ 1 = 1.0 a = 0.01 P r = 7.0 R e = 4900 N b = N t = 0.1 ,  and  L e = 5.0  for a Reiner–Philippoff fluid.
Figure 13. Effect on  % M R  of the presence of polymers as a function of  γ  for fixed  γ 1 = 1.0 a = 0.01 P r = 7.0 R e = 4900 N b = N t = 0.1 ,  and  L e = 5.0  for a Reiner–Philippoff fluid.
Lubricants 11 00339 g013
Figure 14. Effect on % M R  of the presence of polymers as a function of  γ  for different  N b  with  γ 1 = 1.0 λ = 0.7 a = 0.01 P r = 5.0 R e = 4900 N t = 0.1 ,  and  L e = 5.0  for a Reiner–Philippoff fluid.
Figure 14. Effect on % M R  of the presence of polymers as a function of  γ  for different  N b  with  γ 1 = 1.0 λ = 0.7 a = 0.01 P r = 5.0 R e = 4900 N t = 0.1 ,  and  L e = 5.0  for a Reiner–Philippoff fluid.
Lubricants 11 00339 g014
Figure 15. Effect on  % M R  of the presence of polymers as a function of  γ  for different  N t  with  K = 0.3 A = 0.1 a = 0.01 P r = 5.0 R e = 4900 N b = 0.1 ,  and  L e = 5.0  for a Powell–Eyring fluid.
Figure 15. Effect on  % M R  of the presence of polymers as a function of  γ  for different  N t  with  K = 0.3 A = 0.1 a = 0.01 P r = 5.0 R e = 4900 N b = 0.1 ,  and  L e = 5.0  for a Powell–Eyring fluid.
Lubricants 11 00339 g015
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sahreen, A.; Ahmad, A.; Khan, R.; Nawaz, R. Polymer Dispersion Effects on Drag, Heat Transfer, and Mass Transfer in Non-Newtonian Based Nanofluids. Lubricants 2023, 11, 339. https://doi.org/10.3390/lubricants11080339

AMA Style

Sahreen A, Ahmad A, Khan R, Nawaz R. Polymer Dispersion Effects on Drag, Heat Transfer, and Mass Transfer in Non-Newtonian Based Nanofluids. Lubricants. 2023; 11(8):339. https://doi.org/10.3390/lubricants11080339

Chicago/Turabian Style

Sahreen, Ayesha, Adeel Ahmad, Razi Khan, and Rab Nawaz. 2023. "Polymer Dispersion Effects on Drag, Heat Transfer, and Mass Transfer in Non-Newtonian Based Nanofluids" Lubricants 11, no. 8: 339. https://doi.org/10.3390/lubricants11080339

APA Style

Sahreen, A., Ahmad, A., Khan, R., & Nawaz, R. (2023). Polymer Dispersion Effects on Drag, Heat Transfer, and Mass Transfer in Non-Newtonian Based Nanofluids. Lubricants, 11(8), 339. https://doi.org/10.3390/lubricants11080339

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