Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 8521580, 12 pages
http://dx.doi.org/10.1155/2016/8521580

Research Article

Heat and Mass Transfer in a Thin Liquid Film over an Unsteady Stretching Surface in the Presence of Thermosolutal Capillarity and Variable Magnetic Field

Yan Zhang, Min Zhang, and Shujuan Qi

School of Science, Beijing University of Civil Engineering and Architecture, Beijing 100044, China

Received 20 January 2016; Accepted 30 June 2016

Academic Editor: Mohamed Abd El Aziz

Copyright © 2016 Yan Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The heat and mass transfer characteristics of a liquid film which contain thermosolutal capillarity and a variable magnetic field over an unsteady stretching sheet have been investigated. The governing equations for momentum, energy, and concentration are established and transformed to a set of coupled ordinary equations with the aid of similarity transformation. The analytical solutions are obtained using the double-parameter transformation perturbation expansion method. The effects of various relevant parameters such as unsteady parameter, Prandtl number, Schmidt number, thermocapillary number, and solutal capillary number on the velocity, temperature, and concentration fields are discussed and presented graphically. Results show that increasing values of thermocapillary number and solutal capillary number both lead to a decrease in the temperature and concentration fields. Furthermore, the influences of thermocapillary number on various fields are more remarkable in comparison to the solutal capillary number.

1. Introduction

In the recent years, researches on the flow and heat transfer of a liquid film on an unsteady stretching sheet have got more and more attentions for its wide applications. For example, during mechanical forming processes, such as polymer extrusion, melt spinning process, the process of shaping by forcing through a die, wire and fiber coating, and food stuff process, the flow of a liquid film on an unsteady stretching sheet will be met.

In 1970, Crane [1] studied the analytic solutions of two-dimensional boundary layer flow due to a stretching flat elastic sheet. Munawar et al. [2] considered time-dependent flow and heat transfer over a stretching cylinder. Shehzad et al. [3] investigated thermally radiative flow with internal heat generation and magnetic field. Furthermore, many scholars discussed other effects on the flow over a stretching sheet, such as three-dimensional flow [4, 5], heat and mass transfer [6–9], MHD [10], first-order chemical reaction [11], non-Newtonian fluids [12, 13], or different possible combinations of these above effects [14–16]. All of the above studies mainly focus on infinite fluid. In fact, the flow and heat transfer of a finite film are more suitable to describe industrial engineering. The hydrodynamics of the thin liquid film over a stretching sheet were first considered by Wang [17] who reduced the unsteady Navier-Stokes equations to the coupled nonlinear ordinary differential equations by similarity transformation and solved the problem using a kind of multiple shooting method (see Roberts and Shipman [18]). Subsequently, Wang [19] obtained the analytical solutions of a liquid thin film and confirmed the validity of the homotopy analysis method. On the basis of Wang’s work [17], several authors [20–28] explored finite fluid domain of both Newtonian and non-Newtonian fluids using various velocity and thermal boundary conditions. The combined effect of viscous dissipation and magnetic field on the flow and heat transfer in a liquid film over an unsteady stretching surface was presented by Abel et al. [29]. The thermocapillary effect in finite fluid domain was first discussed by Dandapat et al. [30]. Noor and Hashim [31] extended the flow problem to hydromagnetic case.

Marangoni convection is caused by surface-tension gradient at a free liquid-gas or liquid-liquid interface that occurs due to gradient of temperature or concentration in the course of heat or mass transfer. The surface-tension variation on the free liquid surface resulting from the temperature gradient or concentration gradient can induce motion within the fluid called thermocapillary flow or solutal capillary flow (thermal Marangoni convection or solutal Marangoni convection). Pop et al. [32] investigated thermosolutal Marangoni forced convection boundary layers. On the flow field of power-law fluid, Lin et al. [33] analyzed the effect of radiation on Marangoni convection flow and heat transfer in the fluids with variable thermal conductivity. Then, Lin et al. [34] dealt with thermosolutal Marangoni convection flow in the presence of internal heat generation. The surface tension plays an important role on the free liquid surface. Other studies about thermocapillary effect on a thin film can be found in [35–38].

The main objective of our study is to extend previous research to the solutal Marangoni effect and mass diffusion. By means of an exact similarity transformation, governing PDEs are reduced into coupled nonlinear ODEs. And the analytical solutions are obtained using the double-parameter transformation perturbation expansion method [11]. The influences of various relevant parameters such as unsteadiness parameter , Hartmann number Ma, the Prandtl number Pr, the Schmidt number Sc, the thermocapillary number , and the solutal capillary number on the flow field are elucidated through graphs and tables.

2. Mathematical Formulation

2.1. Governing Equations and Boundary Conditions

Consider the thin elastic sheet that emerges from a narrow slit at origin of the Cartesian coordinate system shown in Figure 1. A variable magnetic field normal to the stretching sheet is applied, where is a positive constant. And in the above model, we take concentration into consideration. A thin liquid film with uniform thickness rests on the horizontal sheet. By applying the boundary layer assumptions [39], the governing time-dependent equations for mass, momentum, energy, and concentration are given by

The boundary conditions are

where and are the velocity components of the fluid in the - and -directions, is the time, is the kinematic viscosity, is the electrical conductivity, is the magnetic field, is the density, is the temperature, is the thermal diffusivity, is the concentration, is the mass diffusivity, is dynamic viscosity, and is the uniform thickness of the liquid film. The dependence of surface tension on the temperature and concentration can be expressed as [40]

where and . The fluid motion within liquid film resulted from not only the viscous shear arising from the stretching of the elastic sheet but also the surface-tension gradient. The stretching velocity is assumed to be of the same form as that considered by Wang [17]:

where is a positive constant and denotes the initial stretching rate. With unsteady stretching (i.e., ), however, becomes the representative time scale of the resulting unsteady boundary layer problem. The adopted formulation of the velocity sheet in (5) is valid only for times unless . Also the temperature and the concentration of the surface of the elastic sheet are assumed to vary both along the sheet and with time, respectively, as [22, 30]

where and are the temperature and the concentration at the slit, respectively, and and are the reference temperature and the reference concentration in the case of , respectively. Equation (6) represents a situation in which the sheet temperature and concentration decrease from and at the slot in proportion to . The nonuniform distributions of temperature and concentration cause surface-tension gradient which leads to the fluid flow from lower surface tension to higher surface tension. According to the temperature and concentration boundary layer conditions (6), one can conclude that thermosolutal Marangoni convection flow is in accord with the flow direction of the thin film.

Figure 1: Physics model with interface condition and coordinate system.

2.2. Similarity Transformation

The special forms of the stretching velocity, surface temperature, and fluid concentration in (5) and (6), respectively, are chosen to allow (2) to be converted into a set of ordinary differential equations by means of similarity transformations:

where the similarity variable is given by Wang [19]

and is yet an unknown constant denoting the dimensionless film thickness, defined by

Substituting (7)–(9) into (1)–(3), (1)–(3) can be transformed into the following nonlinear ordinary differential equations:

subject to the boundary conditions

where is Hartmann number, is the dimensionless unsteadiness parameter, is the Prandtl number, is the Schmidt number, is the dimensionless film thickness, and the thermocapillary number and the solutal capillary number are defined as and , respectively.

The skin friction coefficient and the Nusselt number are defined as

where the wall skin friction and the sheet heat flux are

Substituting (14) into (13), then we get

where is the local Reynolds number.

2.3. Solution Approach

In order to solve (10), We used the double-parameter transformation perturbation expansion method that confirmed the validity by Zhang and Zheng [11]. We introduced an artificial small parameter ; the equations can be expanded to the form of power series. We can obtain analytical solutions in the form of series by comparing the coefficients of same power and solving the equation sets.

We transform the dependent variable and independent variable as follows:

where and is an artificial small parameter. We assume , , and ; here , , and are constants; the boundary conditions become

Equation (10) reduces to

And we can assume

with the boundary conditions

Let and , and we can obtain

Substituting (21) into (12), we can get four equations involving the four variables such as , , , and . Analytical solutions of (10)–(12) in the form of series can be obtained with this method.

3. Results and Discussion

The effects of physical parameters on velocity, temperature, and concentration fields are presented in Tables 1–7, respectively. Based on Table 1, the values of dimensionless film thickness and skin friction coefficient obtained are compared with the results by Wang [19] for the case , and good agreement can be seen. From Table 2, it can be observed that increasing the value of unsteadiness parameter will decrease the film thickness , the skin friction coefficient , the heat flux , and the diffusion flux but increase the free velocity , free temperature , and free concentration . In Table 3, as the magnetic parameter Ma increases, the hydrodynamics behavior is the same as that in Table 2 except for the skin friction coefficient .

Table 1: Comparison of values of film thickness β and skin friction coefficient with Ma = 0.

Wang [19]
Present results

0.8
2.151990
2.680940
2.152839
2.681920

1.0
1.543620
1.972385
1.543620
1.972385

1.2
1.127781
1.442631
1.127781
1.442625

1.4
0.821032
1.012784
0.821032
1.012780

Table 2: Values of β, , and , , , and with Ma = 1.2, , , Pr = 1, and Sc = 1.8 and S is varied.

0.8
1.706139
0.245858
0.190409
0.079435
2.791406
2.738673
3.835584

1.0
1.340872
0.387092
0.249442
0.118379
2.180379
2.252363
3.160992

1.2
1.088400
0.541510
0.309461
0.161276
1.714381
1.897617
2.678301

Table 3: Values of β, , and , , , and with , , , Pr = 1, and Sc = 1.8 and Ma is varied.

Ma

0.5
1.955885
0.228035
0.131483
0.047797
2.776933
3.21537
4.463988

1.5
1.629020
0.252600
0.211365
0.092719
2.802801
2.590747
3.640895

3.0
1.370781
0.283080
0.299672
0.154063
2.866680
2.089107
2.983500

Table 4: Values of β, , and , , , and with S = 0.8, Ma = 1, , Pr = 1, and Sc = 1.8 and is varied.

β

1
1.765648
0.242329
0.178412
0.070413
2.783885
2.850602
3.984435

3
1.895449
0.296754
0.143799
0.054085
2.978572
3.080157
4.285635

5
1.987330
0.335046
0.123865
0.044928
3.120442
3.24009
4.498705

8
2.089998
0.380759
0.106127
0.036553
3.279882
3.415187
4.733832

Table 5: Values of β, , and , , , and with , Ma = 1, , Pr = 1, and Sc = 1.8 and is varied.

1
1.765648
0.242329
0.178412
0.070413
2.783885
2.850602
3.984435

3
1.823807
0.265936
0.162696
0.062564
2.870952
2.95405
4.120261

5
1.870447
0.284213
0.148046
0.056928
2.942314
3.03796
4.229606

8
1.927131
0.30837
0.135444
0.050747
3.028503
3.136665
4.360637

Table 6: Values of β, , and , , , and with , Ma = 1, , , and Sc = 1.8 and Pr is varied.

Pr

1
1.765648
0.242329
0.178412
0.070413
2.783885
2.850602
3.984435

2
1.70801
0.216269
0.065849
0.079270
2.701304
4.083229
3.851739

5
1.675911
0.203976
0.009869
0.084604
2.654055
6.569892
3.775766

10
1.671722
0.20159
0.000697
0.085353
2.648593
9.447705
3.766435

Table 7: Values of β, , and , , , and with , Ma = 1, , , and Pr = 1 and Sc is varied.

Sc

0.5
1.884207
0.289866
0.144404
0.313734
2.963305
3.062276
2.011923

1
1.810615
0.259498
0.163349
0.160799
2.852006
2.932347
2.933503

1.5
1.77745
0.245319
0.170494
0.093966
2.80298
2.874491
3.62599

1.8
1.765648
0.242329
0.178412
0.070413
2.783885
2.850602
3.984435

From Tables 4 and 5, it can be seen that and have the same effects on velocity, temperature, and concentration fields. Increasing the value of the thermocapillary number or solutal capillary number can cause a descent in the free temperature and free concentration but a rise in the film thickness , the skin friction coefficient , the heat flux , and the free velocity . Table 6 shows that the values of the free concentration and the heat flux increase while the film thickness , the skin friction coefficient , the diffusion flux , the free velocity , and free temperature decrease with an increment in the Prandtl number Pr. Table 7 shows the effect of Schmidt number on velocity, temperature, and concentration fields.

The velocity, temperature, and concentration profiles for the hydromagnetic flow in a thin film over an unsteady stretching sheet with thermocapillary are presented graphically in Figures 2–7. Figures 2(a)–2(c) show the impacts of unsteadiness parameter on the velocity, temperature, and concentration, respectively. We can find the enhancement in the velocity is more remarkable compared with the temperature and concentration when increment of is the same. And unsteadiness parameter has the similar effects on the temperature and concentration fields. Figures 3(a)–3(c) show the effects of Hartmann number on the velocity, temperature, and concentration fields, which are the same as Table 3. From Figure 3(a), it can be seen that when the Hartmann number increases, at the beginning, the velocity decreases slightly due to the fact that magnetic field produces a drag, and then the velocity near the surface of the thin film rises up gradually under the Marangoni effect, which explains the emergence of the intersection. Meanwhile, the temperature and concentration rise significantly and monotonously as increases (see Figures 3(b) and 3(c)). It can be proved that the magnetic field promotes heat and mass transfer in a thin liquid film.

(a)

(b)

(c)

Figure 2: (a) Velocity profiles for different values of with , , , , and . (b) Temperature profiles for different values of with , , , , and . (c) Concentration profiles for different values of with , , , , and .

(a)

(b)

(c)

Figure 3: (a) Velocity profiles for different values of Ma with , , , , and . (b) Temperature profiles for different values of with , , , , and . (c) Concentration profiles for different values of with , , , , and .

(a)

(b)

(c)

Figure 4: (a) Velocity profiles for different values of with , , , , and . (b) Temperature profiles for different values of with , , , , and . (c) Concentration profiles for different values of with , , , , and .

(a)

(b)

(c)

Figure 5: (a) Velocity profiles for different values of with , , , , and . (b) Temperature profiles for different values of with , , , , and . (c) Concentration profiles for different values of with , , , , and .

Figure 6: Temperature profiles for different values of Pr with , , , , and .

Figure 7: Concentration profiles for different values of with , , , , and .

Figures 4(a) and 5(a) show the velocity profiles for different values of thermocapillary number and solutal capillary number . Similar to Figure 3(a), with increasing value of thermocapillary number or solutal capillary number, one can see that the velocity decreases initially mainly due to the magnetic field effect; then the velocity near the surface of the thin film obviously rises up mainly because of the Marangoni effect; furthermore, the two factors are balanced at the intersection. The effects of thermocapillary number on the temperature and concentration field are shown in Figures 4(b) and 4(c). Results demonstrate that the temperature and concentration both reduce when the thermocapillary number increases. Figures 5(b) and 5(c) illustrate that the fluid temperature and concentration are the decreasing functions of solutal capillary number. The lower value of thermocapillary number has higher temperature and concentration. Furthermore, decrement in the temperature and concentration profiles due to thermocapillary number is more remarkable in comparison to the solutal capillary number .

The effect of Pr on temperature is presented in Figure 6. It can be observed that heat transfer behaviors are strongly dependent on the value of the Prandtl number. Figure 6 shows that the temperature of the fluid decreases monotonously with the increasing of Pr. That is to say, Pr plays an opposite effect on the temperature field compared with and Ma. In addition, the effect of Schmidt number on the concentration distribution is displayed in Figure 7. Increasing the Schmidt number causes a decrease in the concentration profile, and it is clear that weaker mass transfer ratio is obtained with bigger Sc values.

4. Conclusions

In this study, a similarity analysis for thermosolutal capillarity and a magnetic field in a thin liquid film on an unsteady elastic stretching sheet has been studied. The effects of various relevant parameters on the velocity, temperature, and concentration fields have been discussed. The main points of presented analysis are listed below:

(1)The thermocapillary number and solutal capillary number both restrict the heat and mass transfer. Furthermore, the effects of the thermocapillary number on the temperature and concentration are more remarkable than solutal capillary number.(2)The velocity, temperature, and concentration all increase for the higher values of unsteadiness parameter.(3)Increasing Prandtl number can reduce the fluid temperature, and increasing Schmidt number causes a descent in the concentration.(4)Marangoni convection caused by surface-tension gradient creates a promotion in flow near the surface region of the film.

Nomenclature

b: | Positive constant
: | Time (s)
: | -component of velocity (ms⁻¹)
: | -component of velocity (ms⁻¹)
: | Horizontal coordinates (m)
: | Vertical coordinate (m)
: | Velocity similarity variable
: | Mass diffusivity (m²s⁻¹)
: | Magnetic field (T)
: | Initial magnetic field (T)
: | Hartmann number
: | Thermocapillary number
: | Solutal capillary number
: | Prandtl number
: | Dimensionless unsteadiness parameter
: | Stretching surface velocity
: | Temperature (K)
: | Temperature at slit (K)
: | Temperature of slit (K)
: | Reference temperature (K)
: | Concentration (kgm⁻³)
: | Concentration at slit (kgm⁻³)
: | Concentration of slit (kgm⁻³)
: | Reference concentration (kgm⁻³).

Greek Symbols

: | Constants
: | Dimensionless film thickness
: | Thermal diffusivity (m²s⁻¹)
: | Electrical conductivity (s (cm)⁻¹)
: | Similarity variable
: | Dynamic viscosity (Nsm⁻²)
: | Kinematic viscosity (m²s⁻¹)
: | Dimensionless film thickness
: | Temperature similarity variable
: | Concentration similarity variable
: | Artificial small parameter
: | Density (kgm⁻³)
: | Surface tension (Nm⁻¹).

Subscripts

: | Thermal quantity
: | Solutal quantity
: | At the slit
: | At the elastic sheet.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (nos. 21206009 and 21576023) and Funding Project for Academic Human Resources Development in Beijing University of Civil Engineering and Architecture (no. 21221214111).

References

L. J. Crane, “Flow past a stretching plate,” Zeitschrift für angewandte Mathematik und Physik ZAMP, vol. 21, no. 4, pp. 645–647, 1970.

S. Munawar, A. Mehmood, and A. Ali, “Time-dependent flow and heat transfer over a stretching cylinder,” Chinese Journal of Physics, vol. 50, no. 5, pp. 828–848, 2012.

S. A. Shehzad, Z. Abdullah, A. Alsaedi, F. M. Abbasi, and T. Hayat, “Thermally radiative three-dimensional flow of Jeffrey nanofluid with internal heat generation and magnetic field,” Journal of Magnetism & Magnetic Materials, vol. 397, pp. 108–114, 2016.

C. Y. Wang, “The three‐dimensional flow due to a stretching flat surface,” The Physics of Fluids, vol. 27, no. 8, pp. 1915–1917, 1984.

S. Munawar, A. Mehmood, and A. Ali, “Three-dimensional squeezing flow in a rotating channel of lower stretching porous wall,” Computers & Mathematics with Applications, vol. 64, no. 6, pp. 1575–1586, 2012.

M. E. Ali and E. Magyari, “Unsteady fluid and heat flow induced by a submerged stretching surface while its steady motion is slowed down gradually,” International Journal of Heat and Mass Transfer, vol. 50, no. 1-2, pp. 188–195, 2007.

M. M. Rashidi, M. Ali, B. Rostami, P. Rostami, and G.-N. Xie, “Heat and mass transfer for MHD viscoelastic fluid flow over a vertical stretching sheet with considering Soret and Dufour effects,” Mathematical Problems in Engineering, vol. 2015, Article ID 861065, 12 pages, 2015.

P. V. S. Narayana and D. H. Babu, “Numerical study of MHD heat and mass transfer of a Jeffrey fluid over a stretching sheet with chemical reaction and thermal radiation,” Journal of the Taiwan Institute of Chemical Engineers, vol. 59, pp. 18–25, 2016.

A. Mehmood, S. Munawar, and A. Ali, “Cooling of a hot stretching surface in the presence of across mass transfer phenomenon in a channel flow,” Zeitschrift fur Naturforschung A, vol. 69, no. 1-2, pp. 34–42, 2014.

S. A. Shehzad, T. Hayat, and A. Alsaedi, “Three-dimensional MHD flow of Casson fluid in porous medium with heat generation,” Journal of Applied Fluid Mechanics, vol. 9, no. 1, pp. 215–223, 2016.

Y. Zhang and L. C. Zheng, “Analysis of MHD thermosolutal Marangoni convection with the heat generation and a first-order chemical reaction,” Chemical Engineering Science, vol. 69, no. 1, pp. 449–455, 2012.

T. Hayat, A. Safdar, M. Awais, and S. Mesloub, “Soret and Dufour effects for three-dimensional flow in a viscoelastic fluid over a stretching surface,” International Journal of Heat and Mass Transfer, vol. 55, no. 7-8, pp. 2129–2136, 2012.

S. A. Shehzad, Z. Abdullah, F. M. Abbasi, T. Hayat, and A. Alsaedi, “Magnetic field effect in three-dimensional flow of an Oldroyd-B nanofluid over a radiative surface,” Journal of Magnetism & Magnetic Materials, vol. 399, pp. 97–108, 2016.

A. Mastroberardino, “Accurate solutions for viscoelastic boundary layer flow and heat transfer over stretching sheet,” Applied Mathematics and Mechanics, vol. 35, no. 2, pp. 133–142, 2014.

J. A. Khan, M. Mustafa, T. Hayat, and A. Alsaedi, “Three-dimensional flow of nanofluid over a non-linearly stretching sheet: an application to solar energy,” International Journal of Heat and Mass Transfer, vol. 86, pp. 158–164, 2015.

S. A. Shehzad, F. M. Abbasi, T. Hayat, F. Alsaadi, and G. Mousa, “Peristalsis in a curved channel with slip condition and radial magnetic field,” International Journal of Heat and Mass Transfer, vol. 91, pp. 562–569, 2015.

C. Y. Wang, “Liquid film on an unsteady stretching surface,” Quarterly of Applied Mathematics, vol. 48, no. 4, pp. 601–610, 1990.

S. M. Roberts and J. S. Shipman, Two Point Boundary Value Problems: Shooting Methods, Elsevier, New York, NY, USA, 1972.

C. Wang, “Analytic solutions for a liquid film on an unsteady stretching surface,” Heat and Mass Transfer, vol. 42, no. 8, pp. 759–766, 2006.

M. M. Nandeppanavar, K. Vajravelu, M. Subhas Abel, S. Ravi, and H. Jyoti, “Heat transfer in a liquid film over an unsteady stretching sheet,” International Journal of Heat and Mass Transfer, vol. 55, no. 4, pp. 1316–1324, 2012.

Y. Zhang, M. Zhang, and Y. Bai, “Flow and heat transfer of an Oldroyd-B nanofluid thin film over an unsteady stretching sheet,” Journal of Molecular Liquids, vol. 220, pp. 665–670, 2016.

H. I. Andersson, J. B. Aarseth, and B. S. Dandapat, “Heat transfer in a liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 43, no. 1, pp. 69–74, 2000.

C.-H. Chen, “Heat transfer in a power-law fluid film over a unsteady stretching sheet,” Heat and Mass Transfer, vol. 39, no. 8-9, pp. 791–796, 2003.

B. S. Dandapat, B. Santra, and K. Vajravelu, “The effects of variable fluid properties and thermocapillarity on the flow of a thin film on an unsteady stretching sheet,” International Journal of Heat and Mass Transfer, vol. 50, no. 5-6, pp. 991–996, 2007.

I.-C. Liu and H. I. Andersson, “Heat transfer in a liquid film on an unsteady stretching sheet,” International Journal of Thermal Sciences, vol. 47, no. 6, pp. 766–772, 2008.

C. Wang and I. Pop, “Analysis of the flow of a power-law fluid film on an unsteady stretching surface by means of homotopy analysis method,” Journal of Non-Newtonian Fluid Mechanics, vol. 138, no. 2-3, pp. 161–172, 2006.

C.-H. Chen, “Effect of viscous dissipation on heat transfer in a non-Newtonian liquid film over an unsteady stretching sheet,” Journal of Non-Newtonian Fluid Mechanics, vol. 135, no. 2-3, pp. 128–135, 2006.

N. F. Noor, O. Abdulaziz, and I. Hashim, “M{HD} flow and heat transfer in a thin liquid film on an unsteady stretching sheet by the homotopy analysis method,” International Journal for Numerical Methods in Fluids, vol. 63, no. 3, pp. 357–373, 2010.

M. S. Abel, N. Mahesha, and J. Tawade, “Heat transfer in a liquid film over an unsteady stretching surface with viscous dissipation in presence of external magnetic field,” Applied Mathematical Modelling, vol. 33, no. 8, pp. 3430–3441, 2009.

B. S. Dandapat, B. Santra, and H. I. Andersson, “Thermocapillarity in a liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 46, no. 16, pp. 3009–3015, 2003.

N. F. M. Noor and I. Hashim, “Thermocapillarity and magnetic field effects in a thin liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 53, no. 9-10, pp. 2044–2051, 2010.

I. Pop, A. Postelnicu, and T. Groşan, “Thermosolutal Marangoni forced convection boundary layers,” Meccanica, vol. 36, no. 5, pp. 555–571, 2001.

Y. H. Lin, L. C. Zheng, and X. X. Zhang, “Radiation effects on Marangoni convection flow and heat transfer in pseudo-plastic non-Newtonian nanofluids with variable thermal conductivity,” International Journal of Heat and Mass Transfer, vol. 77, no. 4, pp. 708–716, 2014.

Y. H. Lin, L. C. Zheng, X. X. Zhang, L. X. Ma, and G. Chen, “MHD pseudo-plastic nanofluid unsteady flow and heat transfer in a finite thin film over stretching surface with internal heat generation,” International Journal of Heat and Mass Transfer, vol. 84, pp. 903–911, 2015.

E. Lim and Y. M. Hung, “Thermophysical phenomena of working fluids of thermocapillary convection in evaporating thin liquid films,” International Communications in Heat and Mass Transfer, vol. 66, pp. 203–211, 2015.

D. M. Christopher and B.-X. Wang, “Similarity simulation for Marangoni convection around a vapor bubble during nucleation and growth,” International Journal of Heat and Mass Transfer, vol. 44, no. 4, pp. 799–810, 2001.

S. Maity, Y. Ghatani, and B. S. Dandapat, “Thermocapillary flow of a thin nanoliquid film over an unsteady stretching sheet,” Journal of Heat Transfer, vol. 138, no. 4, Article ID 042401, 2016.

E. A. Chinnov, “Wave-thermocapillary effects in heated liquid films at high Reynolds numbers,” International Journal of Heat and Mass Transfer, vol. 71, no. 1, pp. 106–116, 2014.

H. Schichting, Boundary Layer Theory, McGraw-Hill, New York, NY, USA, 6th edition, 1964.

K. Arafune and A. Hirata, “Thermal and solutal marangoni convection in In-Ga-Sb system,” Journal of Crystal Growth, vol. 197, no. 4, pp. 811–817, 1999.