The Thin Film Flow of Walter’s B Fluid over the Surface of a Stretching Cylinder with Heat and Mass Transfer Analysis
Received Date: Jul 23, 2017 / Accepted Date: Jul 31, 2018 / Published Date: Aug 08, 2018
Keywords: Heat and mass transfer flow of the Walter’s B fluids; Liquid film; Unstable stretching cylinder; Skin friction; Local nusselt and sherwood numbers; HAM and numerical method
Thin film flow is one of important natural phenomena. Its enormous usages and application in the expanse of industries, engineering and technology in a last few years clarify its importance in nature. The analysis of thin film flow of practical applications is a challenging interplay between fluid mechanics and fluid dynamic. Wire and fibre coating is one of its important applications. Extrusion of polymer and metal, processing of food stuff, continuous casting, plastic sheets drawing, and fluidization of reactor, exchanges, and chemical processing equipment are some of its common applications. Keeping in view such a countless applications researcher did a lot of work on it. Wang  was the main researcher to investigate liquid film on a time dependent stretching sheet. The similar work was further improved by Ushah and Sridharan . The heat transmission effect on liquid film analysis on horizontal sheet was studied by Liu and Andersson . In their work they used numerical methods to obtain the solution and also discuss the confined flow parameters. Aziz et al.  perceived the consequence of inner heat production of a thin liquid film flow on a time dependent stretched sheet. Recently Tawade et al.  examined thin liquid flow over an un-steady stretched sheet with thermal radiation, in the presence of magnetic field using RK-Fehlberg and Newton-Raphson schemes for solution of non-linear equations.
Thin fluid film flow of Non-Newtonian liquid increases in many life geographies which is used mostly in cylindrical shapes. Several researcher [6-9] investigated Power law fluid by using different cases in unsteady stretching surface. Megahe  and Abolbashari et al.  scrutinized thin film flow of Caisson fluid flows and heat transmission in the existence of inconsistence heat flux and viscid dissipation with slip velocity. Recently Qasim et al.  deliberates the Nano-fluid thin film on a time dependent stretched surface taking the Buongiorno’s model.
Walter’s B fluid is one of the important subclasses of non-Newtonian fluids. Several researchers investigated Walter’s B fluid with effect of heat and MHD. Very less study is available on it in the form of thin films flow in cylinder. Hussain and Ullah  examined boundary layer flow using Walter’s B fluid on a stretched cylinder with temperature reliant viscosity. Manjunatha et al.  studied thermal dusty fluid flow over a permeable stretched cylinder in the occurrence of non-uniform source or sink. Hayat et al.  investigated axisymmetric third order liquid flow on a stretching cylinder within the effect of MHD. Qasim et al.  considered the MHD Boundary Layer having Slip condition Flow under Heat Transition of the Ferrofluid on a Stretching Cylinder with a suggested Heat Flux. Sheikholeslami  studied the consequence of constant value of suction of Nanofluid flows and heat transition over a stretched pipe. Manjunatha et al.  studied the radiation influence and heat transmission of MHD Dusty Fluid flow over a Stretched Cylinder Surrounded in a Porous Medium in the occurrence of Heat Sources. Abdulhameed et al.  considered the oscillating and heat transmission and association of different pressure waveforms in a circular cylinder. Mahdy  examined heat transmission of a Casson fluid flow over a stretched cylinder while considering dufour and soret effects. Abdul Hakeem et al.  observed the influence of heat radioactivity of Walter’s B fluid on a stretched surface with non-uniform heat source/sink and flexible deformation. Pandey et al.  deliberate the Features of Walter’s B Visco-Elastic Nano-fluid Layer and Heat transition.
In the field of engineering most of mathematical problems are complex in nature and the exact solution is almost very difficult or even not possible. So Numerical and Analytical methods are used to find the approximate solution. In these methods Homotopy Analysis Method is one of the important and popular techniques for the solution of such type of problems. It is an auxiliary method and its foremost advantage is to apply on the differential equations which are nonlinear without discretization and linearization. Liao [23-30] investigated this technique for the answer of non-linear equations and successfully proved that this method is rapidly convergent to the approximated solutions. Apart from this, this method provides series solutions in the form of single variable functions. Finding solution with this method is important because it involves all the physical parameters of the problem and we can easily discuss its behaviour. The authors [31-34] used this procedure to answer highly non-linear and coupled equations due to it fast convergence.
As it is already mentioned that most of the published work related to thin films are available about viscous fluids in Cartesian systems, where very less literature is available on the non-Newtonian fluids in Cylindrical system. The advancement of research in the field of thin films and its application in industry and technology is the need of the existing era. Therefore, the determination of existing research is based on the study of non-Newtonian thin film flow in cylindrical system. So, for this aim Walter’s B fluid flow is considered . HAM is used for the solution modelled equations which are nonlinear and coupled. The effect of modelled parameters has been studied graphically .
Consider the steady axi-symmetric and incompressible flow of Walter’s B fluid in a stretched cylinder having radius a alongside axial direction, where z is alongside cylinder and r is along to radial direction. The temperature at the cylinder surface is represented by Tb and the concentration is symbolized by Cb respectively. The stretching velocity of the cylinder is Ww = 2cz, in which c>0 represented the stretched constant, the size of the cylinder condenses but the outside radius of the pipe are taken fixed. is the temperature and concentration fields, in which Tb, Cb is the initial temperature and concentration and Tγ,Cγ is the reference temperature and concentration respectively. Viscous dissipation and the body forces are neglected in the flow field. Also, supposition is made that the base fluid and the nanometer size particle exist in thermal equilibrium, dominate the deficiency of slip.
Considering these suppositions in cylindrical coordinates (γ,z), the mass of conservation equation, the conservation of momentum equation, energy and concentration equations can be identified as:
The applicable Boundary Constrains for the flow pattern are taken as:
Taking u(γ,z)and w(γ,z) are the velocities components, ρ indicates the density of the fluid, k0 indicates the thermal conductivity, p is the pressure, α represents the thermal diffusivity, T denotes the local temperature, σ indicates thickness of the fluid film, Cp is mentioned for specific heat, Cs is used for concentration susceptibility, D indicates the mass diffusivity, C symbolizes concentration field, where is Tm represents the mean temperature of the fluid.
Familiarising the succeeding similarity transformations:
The prime identifies the derivative w.r.t the similarity variable ζ, c>0 indicates stretching constraint, a specifies exterior radius of cylinder where the exterior radius is represented by b. f (ς ),Θ(ς),Φ(ς ) Indicate the non-dimensional velocity profile, temperature gradient and concentration field.
Using the similarity transformation i.e., eqn. (7) in eqns. (1-6), justifies the mass of conservation equation automatically and the left over equations are transformed to non-linear ordinary differential scheme as:
The non-dimensional boundary condition of the problem are:
The physical non-dimensional parameter after simplification:
Where, represents the Reynolds quantities, is Walter’s B fluid constraint, designates the Prandtl number, represents the Dufour number,denotes the Schmidt number and is used for the Soret number.
The Skin friction value, local Nusselt number and Sherwood number can be expressed for the given flow problem as:
The non-dimensional forms for the above mentioned physical properties are:
Where is the local Reynolds number.
Solution by HAM
The boundary condition in eqn. (11), are answered approximately by Homotopy Analysis Method (HAM). The solutions enclose the auxiliary parameters which is used to stabilise and adjust the convergence of solutions eqns. (8-10).
The initial solutions are chosen as follow:
The linear operator for the problem are denoted by
Which has given the following general solution:
Where are constant:
The subsequent non-linear operators for the problem are selected in the below given form:
The basic idea for HAM solution is defined in eqns. (23-30), where is the Zeroth-order system for eqn. (8-10) as:
The equivalent boundary condition for the system of equation are:
Where ρ ∈[0,1] is the embedding constraint, are the axillary constraints use for the convergence of the results. For ρ=0 and ρ=1 we get:
Expand the non-dimensional velocity, temperature and concentration fields by Taylor’s series when.
The axillary parameter are designated in a fashion that the series in eqn. (25) converges at ρ =1 injecting ρ =1 in eqn. (25), we achieve:
The ϖ th - order system for eqn. (8-10) are the following:
The consistent boundary conditions are:
HAM Solution Convergence
The convergence of eqn. (27), exclusively can be influenced by the auxiliary constraints . It is a selection in a manner which controls and converge the series answer Figure 1. The possibility segments of are plotted -curves of and for 20th order approximated HAM solution. The effective region of is . The convergence of the HAM method through -curves for velocity profile, temperature gradient and concentration field has been presented in Figures 2-4 respectively.
Results And Discussion
The performance of various physical constraints, enclosed in the flow field, like Reynolds quantities, Walter’s B fluid parameter, Dufour number, Schmidt numbers, Soret number and Prandtl number on a non-dimensional velocity profile, temperature gradient and concentration field have been presented through graphs. The association of HAM approximation with Numerical method (NDSolve Method) has been shown through graphs and tables. Figures 2-4 show the -curves of f (ς ),Θ(ς ),Φ(ς ) for the 20th-order probable HAM solution. Figure 5 and 6 shows that the convergent region for the problem shrinks when the value of Schmidt number and Soret number are increased. Figure 7 displays the velocity drop for the Walter’s B fluid parameter A. The transport rate is decreased by increasing the value of A, and as a result the boundary-layer thickness is decreased. It is observed that thickness of the boundary-layer and velocity field reduce by rising of Walter’s B fluid parameter A. Whereas the elastic stress constraint rises with the thermal boundary layer thickness. While momentum boundary-layer thickness is decreased by increasing the values of A, consequence that the velocity field at the surface of the cylinder increases. Figure 8 illustrates that Walter’s B fluid constraint has the same effect on the velocity rate of change f '(ς ) . Figure 9 presents the effect of temperature gradient for several values of Walter’s B fluid constraint A and it is detected that temperature drops when Walter’s B fluid parameter A is increased. The velocity distribution and rate of change of velocity field for several numbers of the Reynolds are specified in Figures 10 and 11. It is mentioned that the Reynolds quantity characterises the comparative consequence of the inertia force against the viscous force and due to inertia the velocity is decreased. Actually the Reynolds quantity is the proportion of inertial and viscous forces. So in the Reynolds quantity the inertial force plays the central role as compared to viscous force. Therefore for by increasing larger Reynolds number the velocity decreases to zero and the flow degeneration is very slow. The same behaviour can be observed in Figures 12 and 13 for the temperature gradient and concentration profile i.e., temperature and concentration profile increase when Re increases. The reason for that bigger Reynolds number implies that inertial force exists as an overpowering agent in opposing the viscous force. Thin film constraint shows a special role in the fluid flow. From plot 14, it appears that the velocity profile reduces with larger values of film thickness parameter Figure 14. The more thick is the fluid film the more resistance offers in the flow and leads to retard the movement of the liquid. By increasing values of film thickness parameter the mass of fluid growths and as a result boundary layer thickness is decreased which causes retardation in the flow motion which is inversely related to the mass of the fluid. The thin film parameter plays a vital role in temperature distribution. Temperature reduces with larger magnitudes of thin film parameter which is observable in plot. 15. Heat is absorbed in the fluid film size, therefore, temperature of the fluid falls down as a result cooling effect is produced Figure 15. It is a fact that in non- Newtonian liquids the concentration is comparatively greater than that of the Newtonian liquids. The concentration is strongly dependent on film size parameter. The Concentration profile rises up with increasing values of thin film parameter as exhibited by Figure 16. The reason is that the film thickness has direct relation with viscosity and thermal conductivity which increase the concentration. The ratio of the momentum and thermal diffusivities represented the Prandtl quantity. In Figure 17 it is deliberated that the temperature gradient is a reducing function of Pr quantity. The temperature filed at the face of the cylinder attained some negative values for various numbers of Pr, i.e., heat is continuously transported from the surface of the cylinder to the fluid. Fluids that’s having inferior Prandtl numbers holding upper thermal conductivities, so by this phenomenon heat are transferred from the sheet surface quicker with-in fluids with greater Prandtl quantity Pr. While in Figure 18 the concentration field displays a poor performance beside Prandtl values Pr because of the weakening of boundary layer. The consequence of the Dufour number of the heat and mass transference rate at the surface of the cylinder is demonstrated in Figures 19 and 20. Figure 19 shows the influence of Dufour quantity on the temperature field which reveals that Heat transfer is directly proportional to Dufour number. Actually, the Dufour Number is the proportion of the concentrations difference to temperature difference. When the Dufour Number Du increases then concentration automatically increases consequently temperature falls down. Figure 20 is prepared for the effect of Du number on Concentration. Since Du defines the ratio of the concentrations difference to temperature difference. By increasing Du, the viscosity improves and as a result the concentration increases and plotted that heat transmission rate reduces while the mass transmission rate rises with the increase of Dufour number. The same effects is observed for heat and mass transmission rates over Soret number. i.e., Figure 21 deliberates the effect of Soret number Sγ over temperature gradient. It is seen that the heat transference rate increases with the rise of Soret number Sγ. The cause is that Soret number shows the proportion of the temperatures difference to concentrations difference. This parameter exhibits that when Sγ increases the temperature also increases. It means that the Soret number reductions the rate of heat transmission. Figure 22 displays the contribution of Sγ in concentration. The Schmidt number is defined as the proportion of the viscid boundary films thickness and concentration boundary films thickness, which plays a vital role in heat and mass transmission procedures. When the Schmidt quantity is improved alternatively decreases the classes of concentration field and boundary-layer thickness and a small growth in the fluid temperature. Figure 23 represents the effect of appropriate parameter Sc. The significance of the Schmidt number on the mass transmission rate is demonstrated in Figure 24. It is also shown that the mass transmission rate falls with the increase in Schmidt number. It is due to the reason that the concentration boundary layer reductions with growth of Schmidt number Sc. The convergence of the solutions is shown in the Table 1 and it is observed from the table that momentum equation converges to 13th order of approximations whereas temperature distribution equation converges to 27th order of approximations and concentration equation converge up to 32nd order of approximations. Table 2 displays the performance of several implanted physical constraints on skin friction quantity. It is shown that skin friction quantity rises with the increase in fluid film thickness parameter and Reynolds number while decreases the performance for Walter’s B fluid constant. Table 3 is designed for the influence of various dynamic limitations on local Nusselt number. It is detected that the increase of Nusselt number is directly related to fluid film thickness parameter, Reynolds number and Prandtl number, while it is inversely related to Dufour number. Table 4 indicates the effect of included physical parameters on local Sherwood number and it is investigated that Sherwood number is increased for fluid film thickness parameter, Reynolds number and Schmidt number while decreases for Soret number. The association of HAM and numerical methods for the velocity profile, temperature gradient and concentration field are displayed in Figures 25-27 and Tables 5-7 and shown that these two methods has adjacent solutions.
|Order of approximations||− f '' (1)||−Θ' (1)||−Φ' (1)|
Table 1: Displays the HAM approximated solutions convergent, when , Pr=0.7, Sr=0.2, β=1.1, Sc=0.2, Re=0.02, A=0.01, Du=0.2.
|β||Re||A||f '' (1)|
Table 2: Demonstrates the Numerical results for skin friction factor for several physical constraints, when , Pr=0.7, Sr=0.2, =1.1, Sc=0.2, Re=0.02, A=0.01, Du=0.2.
Table 3: Demonstrates the Numerical values of local Nusselt number for various physical parameters, when , Pr=0.7, Sr=0.2, =1.1, Sc=0.2, Re=0.02, A=0.01, Du=0.2.
Table 4: Demonstrates of the Sherwood quantity for various physical constraints, when , Pr=0.7, Sr=0.2, =1.1, Sc=0.2, Re=0.02, A=0.01, Du=0.2.
|ζ||HAM Result f (ζ)||Numerical result||Absolute error|
Table 5: The association between HAM and ND-Solve for non-dimensional velocity f (ζ) when , Pr=1.0, Sr=0.6, =1.4, Sc=0.6, Re=1.0, A=0.7, Du=0.5.
|ζ||HAM Result Θ(ζ)||Numerical result||Absolute error|
Table 6: The association between HAM and ND-Solve for temperature gradient Θ(ζ) when , Pr=1.0, Sr=0.6, =1.4, Sc=0.6, Re=1.0, A=0.7, Du=0.5.
|ζ||HAM Result Θ(ζ)||Numerical result||Absolute error|
Table 7: The association between HAM and ND-Solve for concentration field Θ(ζ) when , Pr=1.0, Sr=0.6, =1.4, Sc=0.6, Re=1.0, A=0.7, Du=0.5.
In the present work heat transmission distribution is examined analytically during thin film flow of a Walter’s B fluid on a stretched cylinder. Applying an appropriate similarity transformations on the major equations of the flow problem and has been converted to nonlinear ODEs. The obtained nonlinear equations are solved by HAM and numerical (ND solve) methods. The Consequence of the several parameters are deliberated. The following assumptions can be observed from the analytical solution:
1. The velocity profile and rate of change of velocity profile both decreases by increasing the value of Walter’s B fluid constant and Reynolds number.
2. The same effect is concluded in temperature filed for Walter’s B fluid constant and Reynolds number.
3. The effect of Prandtl number is directly related to concentration field and inversely related to temperature filed.
4. The temperature distribution increases by increasing Soret, Dufour and Schmidt numbers during fluid motion.
5. While the concentration field is decreased by increasing the values of Soret, Dufour and Schmidt numbers during fluid motion.
According to the prevailing literature this is the first effort concerning the heat distribution of a thin film flow of a Walter’s B fluid on a stretching cylinder.
All the writers state that they have no challenging securities.
TG and WK modelled the problem and solved. MI participated in the physical discussion of the problem. All authors recited and acknowledged the last manuscript.
- Wang CY (1990) Liquid film on an unsteady stretching surface. Q Appl Math 48: 601-610.
- Usha R, Sridharan R (1995) The Axisymmetric Motion of a Liquid Film on an Unsteady Stretching Surface. J Fluids Eng 117: 81-85.
- Liu IC, Andersson IH (2008) Heat transfer in a liquid film on an unsteady stretching sheet. Int J Thermal Sci47: 766-772.
- Aziz RC, Hashim I, Alomari AK (2011) Thin film flow and heat transfer on an unsteady stretching sheet with internal heating. Meccanica 46: 349-357.
- Tawade J, Abel MS, Metri PG (2016) Thin film flow and heat transfer over an unsteady stretching sheet with thermal radiation, internal heating in presence of external magnetic field. Physics Flu 3: 29-40.
- Andersson HI, Aarseth JB, Braud N, Dandapat BS (1996) Flow of a power-law fluid film on an unsteady stretching surface. J Non-Newtonian Fluid Mech 62: 1-8.
- Chen CH (2003) Heat transfer in a power-law liquid film over an unsteady stretching sheet. Heat and Mass Transfer 39: 791-796.
- Wang C, Pop I (2006) Analysis of the flow of a power-law liquid film on an unsteady stretching surface by means of homotopy analysis method. J Non-Newtonian Fluid Mech 138: 161-172.
- Chen CH (2006) Effect of viscous dissipation on heat transfer in a non-Newtonian liquid film over an unsteady stretching sheet. J Non-Newtonian Fluid Mech 135: 128-135.
- Megahed AM (2015) Effect of slip velocity on Casson thin film flow and heat transfer due to unsteady stretching sheet in presence of variable heat flux and viscous dissipation. Appl Math Mech Engl 36: 1273-1284.
- Abolbashari MH, Freidoonimehr N, Rashidi MM, Abolbashari MH, Nazari F, et al. (2015) Analytical modelling of entropy generation for Casson nano-fluid flow induced by a stretching surface. Adv Powder Technol 26: 542-552.
- Qasim M, Khan ZH, Lopez RJ, Khan WA (2016) Heat and mass transfer in nanofluid thin film over an unsteady stretching sheet using Buongiorno’s model. Eur Phys J Plus 131: 1-16.
- Hussain A, Ullah A (2016) Boundary layer flow of a Walter’s B fluid due to a stretching cylinder with temperature dependent viscosity. Alexandria Eng J 55: 3073-3080.
- Manjunatha PT, Gireesha BJ, Prasannakumara BC (2014) Thermal analysis of conducting dusty fluid flow in a porous medium over a stretching cylinder in the presence of non-uniform source/sink. Int J Mech Mater Eng 9: 1-10.
- Hayat T, Shafiq A, Alsaedi A (2015) MHD axisymmetric flow of third grade fluid by a stretching cylinder. Alexandria Eng J 54: 205-212.
- Qasim M, Khan ZH, Khan WA, Shah IA (2014) MHD Boundary Layer Slip Flow and Heat Transfer of Ferrofluid along a Stretching Cylinder with Prescribed Heat Flux. PLoS ONE 9:1-6.
- Sheikholeslami M (2014) Effect of uniform suction on nanofluid flow and heat transfer over a cylinder. J Braz Soc Mech Sci Eng 37: 1623-1633.
- Manjunatha PT, Gireesha BJ, Prasannakumara BC (2017) Effect of Radiation on Flow and Heat Transfer of MHD Dusty Fluid Over a Stretching Cylinder Embedded in a Porous Medium in Presence of Heat Source. Int J Appl Comput Math 3: 293-310.
- Abdulhameed M, Vieru D, Sharidan S (2016) Comparison of different pressure waveforms for heat transfer performance of oscillating flow in a circular cylinder. Eng Sci Technol Int J 19: 1040-1049.
- Mahdy A (2015) Heat Transfer and Flow of a Casson fluid due to A Stretching Cylinder with the Soret and Dufour Effects. J Eng Phys Thermophys 88: 928-936.
- Hakeem AKA, Ganesh NV, Ganga B (2014) Effect of heat radiation in a Walter’s liquid B fluid over a stretching sheet with non-uniform heat source/sink and elastic deformation. J King Saud Univer- Eng Sci 26: 168-175.
- Pandey SD, Nema KV, Sanddep T (2016) Characteristic of Walter’s B Visco-Elastic Nanofluid Layer Heated from Below. Int J Energy Eng 6: 7-13.
- Liao SJ (1992) The proposed homotopy analysis method for the solution of nonlinear problems. Shangai Jiao Tong University, Shanghai.
- Liao SJ (1999) An Explicit, totally Analytic Approximate Solution for Blasius Viscous Flow Problems. Int J Non-Linear Mech. 34: 759-778.
- Liao SJ (2003) Beyond Perturbation: Introduction to the Homotopy Analysis Method. Chapman and Hall, CRC, Boca Raton, p: 336.
- Liao SJ (2003) On the Analytic Solution of Magnetohydrodynamic Flows of Non-Newtonian Fluids over a Stretching Sheet. J Fluid Mech 488: 189-212.
- Liao SJ (2004) On Homotopy Analysis Method for Nonlinear Problems. Appl Math Comput 147: 499-513.
- Liao SJ (2006) An Analytic Solution of Unsteady Boundary Layer Flows Caused by Impulsively Stretching Plate. Commun Nonlinear Sci Numer Simu 11: 326-339.
- Liao SJ (2012) Homotopy analysis method in nonlinear differential equations. Springer and Higher Education Press Heidelberg.
- Liao SJ (2010) An optimal homotopy-analysis approach for strongly nonlinear differential equations. Comm Nonlinear Sc Num Simulation. 15: 2003-2016.
- Abbasbandy S (2006) The application of homotopy analysis method to nonlinear equations arising in heat transfer. Phys Lett A 360: 109-113.
- Abbasbandy S (2007) Homotopy analysis method for heat radiation equations. Int Commun Heat Mass Trans 34: 380-388.
- Abbasbandy S (2007) The application of homotopy analysis method to solve a generalized Hirota-Satsuma coupled KdV equation. Physics Letters A 361: 478-483.
- Abbasbandy S, Shirzadi A (2011) A new application of the homotopy analysis method: Solving the Sturm-Liouville problems. Communications in Nonlinear Science and Numerical Simulation 16: 112-126.
- Sheikholeslami M (2014) Effect of uniform suction of nanofluid flow and heat transfer over a cylinder. Braz Soc Mech Sci Eng 37: 1623-1633.
- Khan NS, Gul T, Islam S, Khan I, Alqahtani AM, et al. (2017) Magneto hydrodynamic Nano liquid Thin Film Sprayed on a Stretching Cylinder with Heat Transfer. Appl Sci 7: 1-25.
Citation: Khan W, Gul T (2018) The Thin Film Flow of Walter’s B Fluid over the Surface of a Stretching Cylinder with Heat and Mass Transfer Analysis. J Appl Computat Math 7: 411. DOI: 10.4172/2168-9679.1000411
Copyright: © 2018 Khan W, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Select your language of interest to view the total content in your interested language
Share This Article
- Total views: 511
- [From(publication date): 0-0 - Aug 18, 2019]
- Breakdown by view type
- HTML page views: 472
- PDF downloads: 39