MHD dissipative Casson nanofluid liquid film flow due to an unsteady stretching sheet with radiation influence and slip velocity phenomenon

Elham Alali; Ahmed M. Megahed

doi:10.1515/ntrev-2022-0031

Article Open Access

MHD dissipative Casson nanofluid liquid film flow due to an unsteady stretching sheet with radiation influence and slip velocity phenomenon

Elham Alali and Ahmed M. Megahed

Published/Copyright: January 15, 2022

Published by

Become an author with De Gruyter Brill

Submit Manuscript Author Information

From the journal Nanotechnology Reviews Volume 11 Issue 1

Abstract

The problem of non-Newtonian Casson thin film flow of an electrically conducting fluid on a horizontal elastic sheet was studied using suitable dimensionless transformations on equations representing the problem. The thin film flow and heat mechanism coupled with mass transfer characteristics are basically governed by the slip velocity, magnetic field, and the dissipation phenomenon. The present numerical analysis by the shooting method was carried out to study the detailed, fully developed heat and mass transfer techniques in the laminar thin film layer by solving the competent controlling equations with eight dominant parameters for the thin liquid film. Additionally, the predicted drag force via skin-friction coefficient and Nusselt and Sherwood numbers were correlated. In view of the present study, a smaller magnetic parameter or a smaller slip velocity parameter exerts very good influence on the development of the liquid film thickness for the non-Newtonian Casson model. Furthermore, a boost in the parameter of unsteadiness causes an increase in both velocity distribution and concentration distribution in thin film layer while an increase in the same parameter causes a reduction in the film thickness. Likewise, the present results are observed to be in an excellent agreement with those offered previously by other authors. Finally, some of the physical parameters in this study, which can serve as improvement factors for heat mass transfer and thermophysical characteristics, make nanofluids premium candidates for important future engineering applications.

Keywords: nanofluid Casson thin film; unsteady stretching sheet; viscous dissipation; thermal radiation; slip velocity; Joule heating

1 Introduction

The mechanism of heat transfer occurring through a thin liquid film past an elastic sheet is given further attention, because it has enormous applications in industry, many technology fields and engineering disciplines. Typical applications include foodstuff processing reactor fluidization, polymer processing, fiber coating, and solar cell encapsulation. In recent years, due to the immense importance of this field, the rheological physical properties of thin films have attracted the attention of many researchers. Wang [1] pioneered a rare exact solution for unsteady fluid film flow problems and heat transfer properties regarding an accelerating elastic sheet. His novel similarity transformations were found to fulfill the full equations of Navier–Stokes together with the equations of the boundary layer. Thus, the study of the topic of flow of liquid thin film together with the mechanism of heat transfer past a flexible sheet has been the leading topic of a huge number of scientific projects in the past [2,3, 4,5]. Later on, Chen [6] took into account the Marangoni impact and forced and mixed convection phenomenon on the problem which investigate the flow of the fluid and heat transfer within a non-Newtonian power-law thin film due to stretched surface. Liu and Andersson [7] examined the fluid flow problem together with the heat transfer process which occur in the thin liquid film layer which is forced by a flexible sheet surface with variable stretching rate and variable surface temperature. Recently, Abel et al. [8], Noor et al. [9], Noor and Hashim [10], Nandeppanavar et al. [11], Liu and Megahed [12], and Khader and Megahed [13] have examined the heat transfer mechanism through a liquid thin film due to an unsteady stretching sheet for different physical assumptions such as viscous dissipation, thermocapillarity phenomenon, thermal radiation, internal heat generation, and external magnetic field. Among the previous non-Newtonian studied problems, the non-Newtonian Casson model has received more attention. Casson fluid has distinct characteristics in which it can treat fluids which have an infinite viscosity at the absence of the shear rate. This type of fluid is usually essential in bio-engineering activities, food manufacturing, and some of the important metallurgical processes [14,15,16].

Nowadays, nanotechnology is an attractive area of research due to its wide use in traditional industry and technological applications such as polymers, insulating materials, and heat exchange media. Nanoparticles are distinguished by nanometer cereal sizes which possess unrivaled electrical and chemical characteristics. A formidable advantage of this novel scientific field is that nanomaterials can be scattered in fluids, especially the traditional heat transfer fluids [17]. For more detailed information on the topic of nanofluid science and historical aspects of the enhancement mechanism of convective heat transfer, we refer to the paper by Kakaç and Pramuanjaroenkij [18]. As pointed out by Khan and Pop [19], the thermal conductivity of some fluids can be improved by hanging micro (nano) material particles in fluids. Also, they reported that the heat transfer mechanism can be improved by adding microscale particles to the principle fluid. The potential mechanism of heat generation or absorption properties on the nanofluid flow and heat transfer through saturated porous media was examined by Hamad and Pop [20]. On the other hand, they considered that the field of nanofluid technology is one of the serious sciences that may control the next senior technological revolution of this century. The highest precision of the manufacturing product in the microfluid flow due to a heated porous elastic sheet process can be obtained by accurate control of flow and heat transfer [21,22]. Very recent studies on the topic of nanofluid flow and heat transfer with different important physical assumptions can be consulted in refs [23,24, 25,26].

Based on the research above and motivated by the possible technological applications regarding the nanofluid issues of non-Newtonian Casson fluids due to the flexible sheet, the purpose of the present work is to present a numerical solution for the physical problem which describes the simultaneous fluid flow mechanism and heat mass transfer properties of a laminar thin film layer of an electrically conducting non-Newtonian Casson nanofluid subject to an unsteady elastic sheet with radiation mechanism, viscous dissipation phenomenon, and slip velocity. The findings of this theoretical scientific project are presumed to be helpful to the thermo community in the field of nanofluid flow, especially in the cooling process efficiency through some future industrial techniques.

2 Formulation of the problem

First, we must recall the definition of nanofluids from the literature review. The heat transfer liquids that are newly used in more technological applications are nanofluids. This important type of fluid is advanced by hanging micro-sized solid particles in the base fluid. Due to the chaos movement of the hanging particles in a liquid, this motion is described by a Brownian motion with a Brownian diffusion coefficient D B . Also, in nanofluid movement, there is a distinct response for the suspended particles to the temperature gradient, which can be observed in the form of a thermophoretic phenomenon with a thermophoretic diffusion coefficient D T . Here, we assume an electrically conducting, radiative non-Newtonian Casson nanofluid flow, which is characterized by β parameter due to an unsteady elastic sheet simultaneously with the impact of the magnetic field and including the influence of chemical reaction. Additionally, the proposed magnetic field is assumed to depend nonlinearly on the time according to the following relation:

(1) B ( t ) = B 0 1 − a t ,

clearly that at t = 0 the magnetic field strength B ( t ) becomes equal to B 0 and a is a constant with dimension 1 t . Likewise, the influences of Brownian motion, slip velocity phenomenon, and thermophoresis through a sheet are also considered. Further, all physical properties are assumed to be constants. By recalling all of the Boussinesq and Rosseland diffusion approximations, the governing equations for our physical investigation can be written as [27]:

(2) ∂ u ∂ x + ∂ v ∂ y = 0 ,

(3) ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y = ν 1 + 1 β ∂ 2 u ∂ y 2 − σ B 2 ρ f u ,

(4) ∂ T ∂ t + u ∂ T ∂ x + v ∂ T ∂ y = α ∂ 2 T ∂ y 2 + ν c p 1 + 1 β ∂ u ∂ y 2 + σ B 2 ρ f c p u 2 − 1 ρ f c p ∂ q r ∂ y + τ D B ∂ C ∂ y ∂ T ∂ y + D T T 0 ∂ T ∂ y 2 ,

(5) ∂ C ∂ t + u ∂ C ∂ x + v ∂ C ∂ y = D B ∂ 2 C ∂ y 2 + D T T 0 ∂ 2 T ∂ y 2 .

In the overhead equations, u and v are the nanofluid velocity components. Moreover, the symbols α and ρ f represent the fundamental fluid thermal diffusivity and the density of the Casson nanofluid, respectively. Another property for the present physical problem is the ratio of the microparticle heat capacity to the fundamental fluid heat capacity, which is mathematically symbolized by τ . Also, the nanofluid electric conductivity is measured by σ and the nanofluid kinematic viscosity is denoted by ν . Furthermore, from the aforementioned equations we observe a robust coupling between the nanofluid temperature T and the nanoparticle volume fraction C . Here, we must refer to the fact that the present non-Newtonian Casson model is valid for all values of β greater than zero. Also, we must mention that when β → 0 our model reduces to the Newtonian model. Finally, q r which is basically associated with the energy equation denotes the radiative heat flux and it is further defined and simplified linearly previously by everyone who studied the thermal radiation phenomenon, see for example, Raptis [28,29]. The proposed nanofluid is incompressible and hence the boundary conditions can undergo the following equations:

(6) u = U + N 1 1 + 1 β ∂ u ∂ y , v = 0 , T w = T 0 − T r b x 2 2 ν ( 1 − a t ) − 3 2 , C w = C 0 − C r b x 2 2 ν ( 1 − a t ) − 3 2 , at y = 0 ,

(7) ∂ u ∂ y y = h = 0 , ∂ T ∂ y y = h = 0 , ∂ C ∂ y y = h = 0 , ( v ) y = h = d h d t ,

where T w is the sheet temperature, T 0 is the fluid temperature at the slot, T r is the constant reference temperature, C 0 is the fluid concentration at the slot, C w is the nanoparticle volume fraction, and C r is the constant reference nanoparticle volume fraction. The physical meaning which can be concluded here from the first portion of equation (6) is that our model is basically based on the fact of the slip velocity phenomenon which can not be ignored especially for the non-Newtonian models. N 1 is the velocity slip factor which can be presumed to proportion to ( 1 − a t ) 1 2 for the similarity solution, so it can take the form N 1 = N 0 ( 1 − a t ) 1 2 , where N 0 is the initial value for the velocity slip factor ( N 1 = N 0 at t = 0 ). Also, h ( t ) is the uniform thickness of the thin elastic liquid film, which is located on the horizontal stretching sheet (Figure 1).

Figure 1

Physical configuration.

Non-dimensional transformations are normally offered as follows [27]:

(8) η = b ν ( 1 − a t ) y , ψ = ν b ( 1 − a t ) x f ( η ) ,

(9) T = T 0 − T r b x 2 2 ν ( 1 − a t ) 1 − a t θ ( η ) , C = C 0 − C r b x 2 2 ν ( 1 − a t ) 1 − a t ϕ ( η ) ,

where η is the non-dimensional variable, f is the dimensionless stream function, θ is the dimensionless temperature, and ϕ is the dimensionless concentration. It is interesting to note here that both the components of the fluid velocity u and v can be represented as the derivatives of the stream function ψ as follows:

(10) u = ∂ ψ ∂ y , v = − ∂ ψ ∂ x .

Applying the previous non-dimensional transformations help us to reach the following governing ordinary differential equations:

(11) 1 + 1 β f ‴ + f f ″ − f ′ 2 − S η 2 f ″ + f ′ − M f ′ = 0 ,

(12) 1 + R Pr θ ″ + f θ ′ − 2 f ′ θ − S 2 ( 3 θ + η θ ′ ) + Ec 1 + 1 β f ″ 2 + Nt θ ′ 2 + Nb θ ′ ϕ ′ + M Ec f ′ 2 = 0 ,

(13) 1 Sc ϕ ″ + f ϕ ′ − 2 f ′ ϕ − S 2 ( 3 ϕ + η ϕ ′ ) + Nt ScNb θ ″ = 0 ,

together with the following associated boundary conditions [27]:

(14) f ( 0 ) = 0 , f ′ ( 0 ) = 1 + λ 1 + 1 β f ″ , θ ( 0 ) = 1 , ϕ ( 0 ) = 1 ,

(15) f ( γ ) = γ 2 S , f ″ ( γ ) = 0 , θ ′ ( γ ) = 0 , ϕ ′ ( γ ) = 0 .

Now, we observe that the momentum equation (11) contains three parameters, β which is previously defined as a Casson parameter, S = a b which is defined as the unsteadiness parameter, and M = σ B 0 2 ρ b which represents the magnetic number. On the other hand, in the energy equation there are five basic parameters, the radiation parameter R = 16 σ ∗ T 0 3 3 k ∗ κ , the Prandtl number Pr = ν α , the Eckert number Ec = U 2 c p ( T w − T 0 ) , the thermophoresis parameter Nt = τ D T ( T w − T 0 ) ν T 0 , and the Brownian motion parameter Nb = τ D B ( C w − C 0 ) ν . Furthermore, the remaining parameters can be defined as the Schmidt number Sc = ν D B , the slip velocity parameter λ = N 0 b ν , and the dimensionless film thickness γ . Additionally, the value of the film thickness γ can be calculated via equation (8) as follows:

(16) γ = b ν ( 1 − a t ) 1 2 h .

Now, our main aim is to determine the unknown γ . Through the flow of fluid in the thin film layer, the film thickness varies with time and it can be easily obtained after differentiating the previous equation with respect to t , to give:

(17) d h d t = − a γ 2 ν b ( 1 − a t ) 1 2 .

Three important physical quantities, the local skin-friction coefficient Cf x , the wall temperature gradient Nu x , and the wall concentration gradient Sh x , are given by

(18) Cf x Re 1 2 = − 1 + 1 β f ″ ( 0 ) , Nu x Re − 1 2 = − θ ′ ( 0 ) , Sh x Re − 1 2 = − ϕ ′ ( 0 ) ,

where f ″ ( 0 ) , θ ′ ( 0 ) , and ϕ ′ ( 0 ) are dimensionless shear stress, the temperature gradient at the stretching sheet, and the concentration gradient at the stretching sheet, respectively. Also, Re = U w x ν which appears in the last equation is the local Reynolds number.

3 Results and discussion

The obtained systems of equations which appear in equations (11)–(13) are non-linear, coupled, ordinary differential equations, which have no exact (analytical) solution. So, they ought to solve numerically together with the physical boundary condition equations (14)–(15). The shooting method discussed by Conte and de Boor [30] has been confirmed to be efficient and appropriate for the solution of these kinds of equations. So, because of the accuracy of the shooting method, it is employed in this article. Prior to investigating nanofluid flow and heat mass transfer techniques, the accuracy and validity of the obtained results, which resulted from employing the numerical shooting method, are first validated through the following comparison. The numerical data with previously published ones are introduced below in the following table for the thickness film γ and the skin-friction coefficient − f ″ ( 0 ) for assorted values of S . Apparently, an excellent agreement is achieved. After that, results are obtained in Table 1 for various combinations of the physical parameters mentioned above. The curves in Figure 2 show that at lower unsteadiness parameter values S , larger dimensionless film thicknesses are obtained together with a diminishing for both the dimensionless velocity and the velocity values f ′ ( γ ) at free surface. Also, the same figure clearly shows that the velocity of the non-Newtonian nanofluid over the sheet f ′ ( 0 ) is an increasing function of unsteadiness parameter S , which physically means that nanofluid flow with great S grows faster than the flow with small S . It is also found that the dimensionless velocity f ′ ( η ) through the liquid thin film layer enhances with augmenting the unsteadiness parameter. Likewise, at lower unsteadiness parameter values S , smaller temperature θ ( γ ) at the ambient and smaller free surface concentration values ϕ ( γ ) are obtained as observed in Figures 3 and 4, respectively. Indeed, the nanofluid concentration chooses the high values of unsteadiness parameter to pass to its improvement.

$Figure 2 f ′ ( η ) f^{\prime} \left(\eta ) for assorted quantities of S S .$

Figure 2

f ′ ( η ) for assorted quantities of S .

$Figure 3 θ ( η ) \theta \left(\eta ) for assorted quantities of S S .$

Figure 3

θ ( η ) for assorted quantities of S .

$Figure 4 ϕ ( η ) \phi \left(\eta ) for assorted quantities of S S .$

Figure 4

ϕ ( η ) for assorted quantities of S .

Table 1

Thickness film γ and the values of the drag force − f ″ ( 0 ) for assorted value of S with β → ∞ and M = λ = 0

S	Qasim et al. [27]		Present results
	γ	− γ f ″ ( 0 )	γ	− γ f ″ ( 0 )
0.4	4.9814540	5.6494474	4.98145434	5.649439985
0.6	3.1317100	3.7427863	3.13170973	3.742785976
0.8	2.1519940	2.6809656	2.15199391	2.680964459
1.0	1.5436160	1.9723849	1.54361609	1.972384593
1.2	1.1277810	1.4426252	1.12777901	1.442624708
1.4	0.8210320	1.0127802	0.82103192	1.012779957

The velocity profiles as elucidated in Figure 5 indicate that the film thicknesses are larger for the minimal values of magnetic number M in the thin film layer. Also, the same behavior is observed for the fluid velocity along the sheet f ′ ( 0 ) .

$Figure 5 f ′ ( η ) f^{\prime} \left(\eta ) for assorted quantities of M M .$

Figure 5

f ′ ( η ) for assorted quantities of M .

Variations of temperature θ ( η ) as a function of η and concentration profiles ϕ ( η ) for assorted values of the same magnetic number M are introduced in Figures 6 and 7, respectively. It is motivating to observe that, a boost in the value of M is a reason for the enhancement of the distribution of temperature, the concentration distribution, the temperature θ ( γ ) at the free surface, and the free surface concentration ϕ ( γ ) . After carefully having studied the effects of magnetic parameter on the temperature field, it could be observed that the magnetic parameter seems to act as a heating factor in the nanofluid flow.

$Figure 6 θ ( η ) \theta \left(\eta ) for assorted quantities of M M .$

Figure 6

θ ( η ) for assorted quantities of M .

$Figure 7 ϕ ( η ) \phi \left(\eta ) for assorted quantities of M M .$

Figure 7

ϕ ( η ) for assorted quantities of M .

Figure 8 shows the physical model predictions for the thickness of the liquid film and the dimensionless velocity distribution through the thin film region according to the variation in β . It is evident that the Newtonian fluid considerations ( β → ∞ ) have a more worthy impact on the film thickness in which the thickness becomes minimum when the fluid is closer to the Newtonian type ( β → ∞ ). Additionally, a paramount feature of the present results in the same figure is the prediction of both the velocity values f ′ ( γ ) at the free surface and the velocity of the non-Newtonian nanofluid along the sheet f ′ ( 0 ) in which they reach the maximum as β → ∞ .

$Figure 8 f ′ ( η ) f^{\prime} \left(\eta ) for assorted quantities of β \beta .$

Figure 8

f ′ ( η ) for assorted quantities of β .

The increased Casson parameter β resulted in maximum temperature θ ( γ ) at the free surface and extreme free surface concentration ϕ ( γ ) as well as highly distribution for both temperature and concentration inside the liquid film layer as clearly observed in Figures 9 and 10. Physically, the Casson parameter should be sufficiently small, to achieve a satisfactory nanofluid cooling which is very much required, especially in some of the important technological applications.

$Figure 9 θ ( η ) \theta \left(\eta ) for assorted quantities of β \beta .$

Figure 9

θ ( η ) for assorted quantities of β .

$Figure 10 ϕ ( η ) \phi \left(\eta ) for assorted quantities of β \beta .$

Figure 10

ϕ ( η ) for assorted quantities of β .

Figure 11 depicts that by enhancing the values of the slip velocity parameter λ , both the velocity f ′ ( η ) inside thin film layer and the film thickness γ decreases, while the completely contrary behavior is observed for the velocity values f ′ ( γ ) at the free surface.

$Figure 11 f ′ ( η ) f^{\prime} \left(\eta ) for assorted quantities of λ \lambda .$

Figure 11

f ′ ( η ) for assorted quantities of λ .

The profiles of temperature θ ( η ) and the profiles of concentration ϕ ( η ) for various values of the slip velocity parameter λ are demonstrated in Figures 12 and 13, respectively. It is noted that both the temperature and the concentration of the fluid flow increase with the increasing trend in λ through the film layer. Besides, we observe that the slip velocity parameter λ has considerable effects on both the temperature θ ( γ ) at the free surface and the free surface concentration ϕ ( γ ) . When λ increases, both θ ( γ ) and ϕ ( γ ) enhance. Physically, within the film layer, both the nanofluid warming and the nanofluid concentration depend strongly on the great slip velocity parameter.

$Figure 12 θ ( η ) \theta \left(\eta ) for assorted quantities of λ \lambda .$

Figure 12

θ ( η ) for assorted quantities of λ .

$Figure 13 ϕ ( η ) \phi \left(\eta ) for assorted quantities of λ \lambda .$

Figure 13

ϕ ( η ) for assorted quantities of λ .

From Figure 14, it is clear that for big values of the parameter of radiation R , the temperature curve inside the thin film layer differs little from the small values of R ; however, when the parameter of radiation R enhances, these profiles show enhancement behavior for both the temperature throughout the film layer and the temperature θ ( γ ) at the free surface as well as a fixed thin film thickness. The physical interpretation for the following behavior is that the thermal radiation phenomenon can also be utilized to predict high warming rates of purely viscous nanofluids.

$Figure 14 θ ( η ) \theta \left(\eta ) for assorted quantities of R R .$

Figure 14

θ ( η ) for assorted quantities of R .

The impact of Eckert number Ec on the temperature curves is depicted in Figure 15. From these plots, we observe that the effect of viscous dissipation phenomenon results in promoting the thermal energy throughout the film layer as well as the temperature θ ( γ ) at the free surface. These observed results hold good for assorted values of the Eckert number. Clearly that, the Eckert number appears to act as a stimulating factor for warming the film region in which it has a very significant influence on the thermal film layer.

$Figure 15 θ ( η ) \theta \left(\eta ) for assorted quantities of Ec.$

Figure 15

θ ( η ) for assorted quantities of Ec.

Temperature distributions θ ( η ) and concentration distributions ϕ ( η ) for various values of thermophoresis parameter Nt are shown in Figures 16 and 17, respectively. Clearly that a larger thermophoresis parameter results in a higher dimensionless free surface temperature θ ( γ ) and dimensionless free surface concentration ϕ ( γ ) . In addition, θ ( η ) reaches its minimum value at Nt = 0 .

$Figure 16 θ ( η ) \theta \left(\eta ) for assorted quantities of Nt.$

Figure 16

θ ( η ) for assorted quantities of Nt.

$Figure 17 ϕ ( η ) \phi \left(\eta ) for assorted quantities of Nt.$

Figure 17

ϕ ( η ) for assorted quantities of Nt.

Figures 18 and 19 display, respectively, the variation of θ ( η ) and ϕ ( η ) with the Brownian parameter Nb. It is observed that at a given location η , θ ( η ) enhances with augment in the Brownian parameter Nb. Furthermore, it is clear from Figure 18 that the increasing Brownian parameter would reason the free surface temperature θ ( γ ) to enhance. Moreover, the impact of the same parameter on the dimensionless concentration ϕ ( η ) has quite opposite trend as illustrated in Figure 19.

$Figure 18 θ ( η ) \theta \left(\eta ) for assorted quantities of Nb.$

Figure 18

θ ( η ) for assorted quantities of Nb.

$Figure 19 ϕ ( η ) \phi \left(\eta ) for assorted quantities of Nb.$

Figure 19

ϕ ( η ) for assorted quantities of Nb.

Table 2 gives the values of Cf x Re 1 2 , Nu x Re − 1 2 , and Sh x Re − 1 2 , which represent, respectively, the values of skin-friction coefficient, the values of the local Nusselt number, and the values of the local Sherwood number, which is defined in equation (18) for several values of the governing parameters. Tabular data elucidate that the values of Nu x Re − 1 2 decreased for all values of λ , R , Ec, Nt, and Nb. Thus, we find that there is a drastic change in the heat transfer characteristics of the Casson thin film flow when the aforementioned several effects are taken into account. Both the unsteadiness parameter S and Casson parameter β tend to decrease the values of Cf x Re 1 2 and thereby also the local Sherwood number. Furthermore, the impact of magnetic number on Cf x Re 1 2 is found to be more noticeable at large values, whereas a totally opposite trend is noted for the velocity slip parameter. Likewise, both the values of Nu x Re − 1 2 and the values of Sh x Re − 1 2 decrease when λ enhances. The cause of this physical behavior is that the thin film thickness becomes thin as the slip velocity parameter increases, which results in smaller values for both Nu x Re − 1 2 and Sh x Re − 1 2 at the surface. Finally, it is found that the local Sherwood number reduces with the growth of the magnetic number. However, it slightly increases with increase in both the radiation parameter and the Eckert number.

Table 2

Values for Cf x Re 1 2 , Nu x Re − 1 2 , and Sh x Re − 1 2 for various values of S , M , β , λ , R , Ec, Nt, and Nb with Pr = 3.97 and Sc = 0.6

S	M	β	λ	R	Ec	Nt	Nb	Cf x ( Re x ) 1 2	Nu x Re − 1 2	Sh x Re − 1 2
0.8	1.0	0.5	0.2	0.5	0.1	0.01	0.5	1.626601	1.941910	1.116211
1.0	1.0	0.5	0.2	0.5	0.1	0.01	0.5	1.544731	2.105241	1.130440
1.2	1.0	0.5	0.2	0.5	0.1	0.01	0.5	1.392272	2.287381	1.077141
1.5	1.0	0.5	0.2	0.5	0.1	0.01	0.5	1.017640	2.405210	0.818659
1.2	0.0	0.5	0.2	0.5	0.1	0.01	0.5	1.252251	2.301941	1.238081
1.2	1.0	0.5	0.2	0.5	0.1	0.01	0.5	1.392272	2.287381	1.077141
1.2	2.0	0.5	0.2	0.5	0.1	0.01	0.5	1.484881	2.242960	0.935370
1.2	0.5	0.2	0.2	0.5	0.1	0.01	0.5	1.541242	2.260850	1.216691
1.2	0.5	0.5	0.2	0.5	0.1	0.01	0.5	1.330440	2.296791	1.156580
1.2	0.5	1.0	0.2	0.5	0.1	0.01	0.5	1.197443	2.317281	1.103971
1.2	0.5	0.5	0.0	0.5	0.1	0.01	0.5	2.478701	2.340581	1.405230
1.2	0.5	0.5	0.2	0.5	0.1	0.01	0.5	1.330440	2.296791	1.156580
1.2	0.5	0.5	0.5	0.5	0.1	0.01	0.5	0.718359	2.225960	0.796689
1.2	0.5	0.5	0.2	0.0	0.1	0.01	0.5	1.330440	2.673711	1.154312
1.2	0.5	0.5	0.2	0.5	0.1	0.01	0.5	1.330440	2.296791	1.156580
1.2	0.5	0.5	0.2	1.0	0.1	0.01	0.5	1.330440	2.037012	1.158091
1.2	0.5	0.5	0.2	0.5	0.0	0.01	0.5	1.330440	2.335921	1.156342
1.2	0.5	0.5	0.2	0.5	0.5	0.01	0.5	1.330440	2.140260	1.157521
1.2	0.5	0.5	0.2	0.5	1.0	0.01	0.5	1.330440	1.944522	1.158710
1.2	0.5	0.5	0.2	0.5	0.1	0.0	0.5	1.330440	2.309911	1.168690
1.2	0.5	0.5	0.2	0.5	0.1	0.1	0.5	1.330440	2.185721	1.054761
1.2	0.5	0.5	0.2	0.5	0.1	0.2	0.5	1.330440	2.113840	1.002950
1.2	0.5	0.5	0.2	0.5	0.1	0.01	0.1	1.330440	2.700031	1.094920
1.2	0.5	0.5	0.2	0.5	0.1	0.01	0.2	1.330440	2.592491	1.133412
1.2	0.5	0.5	0.2	0.5	0.1	0.01	0.3	1.330440	2.489610	1.146401

4 Conclusion

The theoretical analysis presented in this work accounts for viscous dissipation, thermophoresis, and the slip velocity effects on the thin film flow and heat mass transfer along the film layer for non-Newtonian Casson fluid which exposed to thermal radiation and magnetic field. This analysis represents an improvement over previous studies. Although strictly no exact solution can be obtained for the constitutive equations used here, the numerical shooting method that is employed in the present analysis has been demonstrated to yield useful results. The numerical procedure via the shooting method is given so that estimations of the film thickness, skin friction coefficient, local Nusselt number, local Sherwood number, and the velocity profile can be conveniently carried out. Comparison of the film thickness and skin-friction coefficient for the Casson thin film flow considered indicates that the results obtained exhibit good accuracy. The numerical calculations via shooting method have shown that:

The unsteadiness parameter, magnetic parameter, slip velocity parameter, and the Casson parameter have a very significant influence on the development of the film thickness.
The magnetic parameter, the unsteadiness parameter, the slip velocity parameter, and the Casson parameter have very significant effects on decreasing the rate of mass transfer.
The radiation parameter, Brownian parameter, Eckert number, and thermophoresis parameter do have a significant effect on the results, but the magnitude of the effect appears to be directly linked to the rate of heat transfer.
As a future work, we can apply an analytic approach such as homotopy analysis method to obtain convergent series solutions, especially for the non-Newtonian nanofluid problems under different physical assumptions.

Acknowledgements

The authors are thankful to the reviewers for their useful suggestions and comments which improved the quality of this article.

Funding information: The authors state no funding involved.
Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.
Conflict of interest: The authors state no conflict of interest.

References

[1] Wang CY. Liquid film on an unsteady stretching surface. Quart Appl Math. 1990;48:601–10. 10.1090/qam/1079908Search in Google Scholar

[2] Andersson HI, Aarseth JB, Braud N, Dandapat BS. Flow of a power-law fluid on an unsteady stretching surface. J Non-Newtonian Fluid Mech. 1996;62:1–8. 10.1016/0377-0257(95)01392-XSearch in Google Scholar

[3] Andersson HI, Aarseth JB, Dandapat BS. Heat transfer in a liquid film on an unsteady stretching surface. Int J Heat Mass Transfer. 2000;43:69–74. 10.1016/S0017-9310(99)00123-4Search in Google Scholar

[4] Chen CH. Heat transfer in a power-law Fluid film over a unsteady stretching sheet. Heat Mass Transfer. 2003;39:791–6. 10.1007/s00231-002-0363-2Search in Google Scholar

[5] Dandapat BS, Santra B, Andersson HI. Thermocapillarity in a liquid film on an unsteady stretching surface. Int J Heat Mass Transfer. 2003;46:3009–15. 10.1016/S0017-9310(03)00078-4Search in Google Scholar

[6] Chen CH. Marangoni effects on forced convection of power-law liquids in a thin film over a stretching surface. Phys Lett A. 2007;370:51–7. 10.1016/j.physleta.2007.05.024Search in Google Scholar

[7] Liu IC, Andersson HI. Heat transfer in a liquid film on an unsteady stretching sheet. Int J Therm Sci. 2008;47:766–72. 10.1016/j.ijthermalsci.2007.06.001Search in Google Scholar

[8] Abel MS, Mahesha N, Tawade J. Heat transfer in a liquid film over an unsteady stretching surface with viscous dissipation in presence of external magnetic field. Appl Math Model. 2009;33:3430–41. 10.1016/j.apm.2008.11.021Search in Google Scholar

[9] Noor NFM, Abdulaziz O, Hashim I. MHD flow and heat transfer in a thin liquid film on an unsteady stretching sheet by the homotopy analysis method. Int J Numer Meth Fluids. 2010;63:357–73. 10.1002/fld.2078Search in Google Scholar

[10] Noor NFM, Hashim I. Thermocapillarity and magnetic field effects in a thin liquid film on an unsteady stretching surface. Int J Heat Mass Transfer. 2010;53:2044–51. 10.1016/j.ijheatmasstransfer.2009.12.052Search in Google Scholar

[11] Nandeppanavar MM, Vajravelu K, Abel MS, Ravi S, Jyoti, H. Heat transfer in a liquid film over an unsteady stretching sheet. Int J Heat Mass Transfer. 2012;55:1316–24. 10.1016/j.ijheatmasstransfer.2011.09.007Search in Google Scholar

[12] Liu IC, Megahed AM. Numerical study for the flow and heat transfer in a thin liquid film over an unsteady stretching sheet with variable fluid properties in the presence of thermal radiation. J Mech. 2012;28:291–7. 10.1017/jmech.2012.32Search in Google Scholar

[13] Khader MM, Megahed AM. Numerical studies for flow and heat transfer of the Powell-Eyring fluid thin film over an unsteady stretching sheet with internal heat generation using the Chebyshev finite difference method. J Appl Mech Tech Phys. 2013;54:440–50. 10.1134/S0021894413030139Search in Google Scholar

[14] Hussain SM, Jamshed W, Kumar V, Kumar V, Sooppy Nisar K, Eid MR, et al. Computational analysis of thermal energy distribution of electromagnetic Casson nanofluid across stretched sheet: Shape factor effectiveness of solid-particles. Energy Reports. 2021;7:7460–77. 10.1016/j.egyr.2021.10.083Search in Google Scholar

[15] Jamshed W, Devi SSU, Safdar R, Redouane F, Nisar KS, Eid MR. Comprehensive analysis on copper-iron (II, III)/oxide-engine oil Casson nanofluid flowing and thermal features in parabolic trough solar collector. J Taibah Univ Sci. 2021;15:619–36. 10.1080/16583655.2021.1996114Search in Google Scholar

[16] Alotaibi H, Althubiti S, Eid MR, Mahny KL. Numerical treatment of MHD flow of Casson nanofluid via convectively heated non-linear extending surface with viscous dissipation and suction/injection effects. Comput Mater Continua. 2021;66:229–45. 10.32604/cmc.2020.012234Search in Google Scholar

[17] Vasu V, Rama Krishna K, Kumar ACS. Analytical prediction of forced convective heat transfer of fluids embedded with nanostructured materials (nanofluids). Pramana J Phys. 2007;69:411–21. 10.1007/s12043-007-0142-1Search in Google Scholar

[18] Kakaç S, Pramuanjaroenkij A. Review of convective heat transfer enhancement with nanofluids. Int J Heat Mass Transf. 2009;52:3187–96. 10.1016/j.ijheatmasstransfer.2009.02.006Search in Google Scholar

[19] Khan WA, Pop I. Boundary-layer flow of a nanofluid past a stretching sheet. Int J Heat Mass Transf. 2010;53:2477–83. 10.1016/j.ijheatmasstransfer.2010.01.032Search in Google Scholar

[20] Hamad MA, Pop I. Scaling transformations for boundary layer flow near the stagnation-point on a heated permeable stretching surface in a porous medium saturated with a nanofluid and heat generation/absorption effects. Transport Porous Media. 2011;87:25–39. 10.1007/s11242-010-9683-8Search in Google Scholar

[21] Hamad MAA, Ferdows M. Similarity solution of boundary layer stagnation-point flow towards a heated porous stretching sheet saturated with a nanofluid with heat absorption/generation and suction/blowing: A Lie group analysis. Commun Nonlinear Sci Numer Simulat. 2012;17:132–40. 10.1016/j.cnsns.2011.02.024Search in Google Scholar

[22] Makinde OD, Khan WA, Khan ZH. Buoyancy effects on MHD stagnation point flow and heat transfer of a nanofluid past a convectively heated stretching/shrinking sheet. Int J Heat Mass Transf. 2013;62:526–33. 10.1016/j.ijheatmasstransfer.2013.03.049Search in Google Scholar

[23] Shamshuddin MD, Abderrahmane A, Koulali A, Eid MR, Shahzad F, Jamshed W. Thermal and solutal performance of Cu/CuO nanoparticles on a non-linear radially stretching surface with heat source/sink and varying chemical reaction effects. Int Commun Heat Mass Transf. 2021;129:105710. 10.1016/j.icheatmasstransfer.2021.105710Search in Google Scholar

[24] Amine BM, Redouane F, Mourad L, Jamshed W, Eid MR, Al-Kouz W. Magnetohydrodynamics natural convection of a triangular cavity involving Ag-MgO/water hybrid nanofluid and provided with rotating circular barrier and a quarter circular porous medium at its right-angled corner. Arabian J Sci Eng. 2021;46:12573–97. 10.1007/s13369-021-06015-6Search in Google Scholar

[25] Jamshed W, Nasir NAAM, Mohamed Isa SSP, Safdar R, Shahzad F, Nisar KS, et al. Thermal growth in solar water pump using Prandtl-Eyring hybrid nanofluid: a solar energy application. Scientific Reports. 2021;11:18704. 10.1038/s41598-021-98103-8Search in Google Scholar

[26] Jamshed W, Eid MR, Nisar KS, Nasir NAAM, Edacherian A, Saleel CA, et al. A numerical frame work of magnetically driven Powell-Eyring nanofluid using single phase model. Scientific Reports. 2021;11:16500. 10.1038/s41598-021-96040-0Search in Google Scholar

[27] Qasim M, Khan ZH, Lopez RJ, Khan WA. Heat and mass transfer in nanofluid over an unsteady stretching sheet using Buongiorno’s model. Eur Phys J Plus. 2016;131:1–16. 10.1140/epjp/i2016-16016-8Search in Google Scholar

[28] Raptis A. Flow of a micropolar fluid past a continuously moving plate by the presence of radiation. Int J Heat Mass Transf. 1998;41:2865–6. 10.1016/S0017-9310(98)00006-4Search in Google Scholar

[29] Raptis A. Radiation and viscoelastic flow. Int Comm Heat Mass Transfer. 1999;26:889–95. 10.1016/S0735-1933(99)00077-9Search in Google Scholar

[30] Conte SD, De Boor C. Elementary numerical analysis: An algorithmic approach. New York: McGraw-Hill; 1972. Search in Google Scholar

Received: 2021-10-30

Revised: 2021-12-04

Accepted: 2022-01-03

Published Online: 2022-01-15

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

https://doi.org/10.1515/ntrev-2022-0031

Keywords for this article

nanofluid Casson thin film; unsteady stretching sheet; viscous dissipation; thermal radiation; slip velocity; Joule heating

Creative Commons

BY 4.0