# Dual Solutions on Thermosolutal Marangoni Forced Convection Boundary Layer with Suction and Injection

Authors: R. A. Hamid; N. M. Arifin; R. Nazar; F. M. Ali; I. Pop
Journal: Mathematical Problems in Engineering
 (2011)

---

## Abstract

This paper considers the extended problem of the thermosolutal Marangoni forced convection boundary layer by Pop et al. (2001) when the wall is permeable, namely, there is a suction or injection effect. The governing system of partial differential equations is transformed into a system of ordinary differential equations, and the transformed equations are solved numerically using the shooting method. The effects of suction or injection parameterf0 on the velocity, temperature, and concentration profiles are illustrated and presented in tables and figures. It is shown that dual solutions exist for the similarity parameter β less than 0.5.


---

## Body

## 1. Introduction

Thermosolutal Marangoni forced convection boundary layer flow refers to the thermal and solutal concentration in Marangoni forced convection boundary layer flow due to the surface tension gradients. The study of Marangoni convection has attracted the interest of many researchers in recent years. This is mainly because of its vast contributions in the industrial field especially in the art work of dyeing on the ground (Kuroda [1]) and in the field of crystal growth (Arafune and Hirata [2]). Some of the relevant works are done by Christopher and Wang [3] who studied the effects of Prandtl number on the Marangoni convection over a flat plate. They have also presented a similarity solution for Marangoni flow for both the momentum and the energy equations assuming a developing boundary layer along a surface. Besides, Chamkha et al. [4] have dealt with a steady coupled dissipative layer, called Marangoni mixed convection boundary layer. The mixed convection boundary layer is generated when besides the Marangoni effects there are also buoyancy effects due to the gravitational and external pressure gradient effects. Furthermore, in the following papers, numerical solutions on Marangoni boundary layers in various geometries were discovered, analyzed, and discussed, for instance, papers by Golia and Viviani [5], Dressler and Sivakumaran [6], Al-Mudhaf and Chamkha [7], and Magyari and Chamkha [8]. An excellent paper on Marangoni boundary layer flow along the interface of the immiscible nanofluid by Arifin et al. [9] has been recently published.  They show that the results  indicate that dual solutions exist when β<0.5. The paper complements also the work by Golia and Viviani [5] concerning the dual solutions in the case of adverse pressure gradient.Investigations into the effects of suction and injection on the boundary layer flow have generated the interest of many researchers nowadays. Pop and Watanabe [10] analyzed the effects of uniform suction or injection on the boundary layer flow and heat transfer on a continuous moving permeable surface.  Meanwhile, Hamza [11] obtained the similarity solution for a flow between two parallel plates (rectangular or circular) approaching or receding from each other with suction or injection at the porous plate. Shojaefard et al. [12] investigated the numerical study concerning flow control by suction and injection on a subsonic airfoil. They discovered that the surface suction can significantly increase the lift coefficient; meanwhile, the injection decreases the skin friction.The existence of dual solutions or second solutions has aroused a lot of interest of many authors. We mention here some of the papers that discovered the dual solution in their problems. For example, papers by de Hoog et al. [13], Ingham [14], Ramachandran et al. [15], Xu and Liao [16] and Ishak et al. [17–19]. Ridha [20] considered the dual solutions of two coupled third-degree nonlinear ordinary differential equations associated with the incompressible viscous laminar flow along a corner. Further, Xu and Liao [16], obtained the dual solutions of boundary layer flow over an upstream moving plate using the homotopy analysis method (HAM). Ishak et al. [19] discovered the dual solution of the classical Blasius problem. In their study, they considered the boundary layer flow over a static flat plate and introduced a new parameter, namely, the velocity ratio parameter to analyze the case when both the flat plate and the free stream are in moving situations. Dual solutions are found to exist when the plate and the free stream move in the opposite directions.In this paper, we extend the problem by Pop et al. [21] by taking the effects of suction or injection for the boundary layer flow. We would also like to investigate whether the dual solutions exist in the above problem. To validate our findings, we have compared the present results with those of Pop et al. [21] for the case of impermeable surface. It is found that the results are in very good agreement. The effects of suction and injection parameter on the surface velocity, temperature, and concentration as well as on velocity, temperature, and concentration profiles are discussed and illustrated in tables and figures.


## 2. Mathematical Formulation

We consider the steady two-dimensional flow along the interfaceS of two Newtonian immiscible fluids, where x and y are the axes of a Cartesian coordinate system. We also assume that the temperature and concentration at the interface are Ts(x) and Cs(x), respectively. Besides, the surface is assumed to be permeable so as to allow for possible suction or injection at the wall (Al-Mudhaf and Chamkha [7]). Under the usual boundary-layer approximations, the basic governing equations are (see Pop et al. [21])(2.1)∂u∂x+∂v∂y=0,u∂u∂x+v∂u∂y=ueduedx+ν∂2u∂y2,u∂T∂x+v∂T∂y=α∂2T∂y2,u∂C∂x+v∂C∂y=D∂2C∂y2,with the boundary conditions,(2.2)v=vw,T=Ts(x),C=Cs(x),μ∂u∂y=σT∂T∂x+σC∂C∂xony=0,(2.3)u⟶ue(x),T⟶TmC⟶Cmasy⟶∞,
where u and v are the velocity components along x and y axes,  T is the fluid temperature, C is the solutal concentration, ue(x) is the external velocity, α is the thermal diffusivity, ν is the kinematic viscosity, D is the mass diffusivity, and vw is the constant suction (vw>0) or injection (vw<0) velocity, μ is the dynamic viscosity, and σT and σC are the rates of change of surface tension with temperature and solute concentration, respectively. As mentioned by Pop et al. [21], the forth condition of (2.2) represents the Marangoni coupling conditions at the interface (balance of the surface tangential momentum), having considered for the surface tension the linear relation given by(2.4)σ=σm-σT(T-Tm)-σC(C-Cm),σT=-∂σ∂T,σC=-∂σ∂C.
Following Pop et al. [21], we define the following nondimensional variables: (2.5)x=L0+XL,y=δLY,u=UcU,v=δUcV,T=Tm+tΔT,C=Cm+ϕΔC,ue(x)=UcUe(x),vw=UcLδ2Vw,
where  L0 locates the origin of the curvilinear abscissa x, L is the extension of the relevant interface S,andΔT and ΔC are positive increments of temperature and solute concentration linked to the temperature and solute concentration gradients imposed on the interface, respectively. Further, δ is a scale factor in the direction normal to the interface, and Uc is the reference velocity which is defined as δ=Re-1/3 and Uc=ν/(Lδ2) with Re=σTΔTL/νμ being the Reynolds number. Substituting (2.5) into (2.1), we obtain the following nondimensional equations: (2.6)∂U∂X+∂V∂Y=0,U∂U∂X+V∂U∂Y=UedUedX+∂2U∂Y2,U∂t∂X+V∂t∂Y=1Pr∂2t∂Y2,U∂ϕ∂X+V∂ϕ∂Y=1Sc∂2ϕ∂Y2,
and the boundary conditions (2.2) and (2.3) reduce to(2.7)V=Vw,t=ts(X),ϕ=ϕs(X),∂U∂Y=∂t∂X+ε∂ϕ∂XonY=0,U⟶Ue(X),t⟶0,ϕ⟶0asY⟶∞,
where Pr is the Prandtl number, Sc is the Schmidt number, and ε is the Marangoni parameter and is defined as(2.8)ε=MaCMaT=σCΔCσTΔT,
with MaC=σCΔCL/μα and MaT=σTΔTL/μα being the solutal Marangoni number and thermal Marangoni number, respectively. It is also noted that positive Vw is for fluid suction and negative for fluid injection at the wall. Equations (2.6) can be transformed into the corresponding ordinary differential equations by the following transformations (see Pop et al. [21]):
(2.9)U=u0X(2β-1)/3f′(η),V=13u0l0x(β-2)/3((2-β)ηf′(η)-(1+β)f(η)),t=-t0Xβθ(η),ϕ=-c0Xβh(η),Ue=u0X(2β-1)/3,η=Yl0x(2-β)/3,
where f′(η),θ(η), and h(η) represent the velocity, temperature, and concentration profiles in the similarity plane and η being the similarity variable. The constant scale factors u0, t0,c0, and l0 are chosen in order to simplify the equations, and these must satisfy the following conditions:(2.10)u0l02=3β+1,t0l0u0=1β,c0t0=1.
If we take t0=1 or c0=1 (due to arbitrariness of the temperature difference ΔT or the solutal concentration difference ΔC), then l0 and u0 are uniquely determined as(2.11)l0=(31+β)1/3β-1/3,u0=(31+β)1/3β2/3.
The transformed ordinary differential equations are(2.12)f′′′+ff′′+2β-1β+1(1-f′2)=0,1Prθ′′+fθ′-3β1+βf′θ=0,1Sch′′+fh′-3β1+βf′h=0,
along with the boundary conditions(2.13)f(0)=f0,f′′(0)=-1-ε,θ(0)=1,h(0)=1,f′(∞)=1,θ(∞)=0,h(∞)=0,
where f0(>0) is the constant suction parameter and f0(<0) is the constant injection parameter.


## 3. Results and Discussion

The system of transformed governing equations (2.12) along with the boundary equations (2.13) is solved numerically using the shooting method. Tables 1, 2, and 3 illustrate the influence of the suction and injection parameter f0=1,0, and −1 on the surface velocity, f′(0), surface temperature, -θ′(0), and surface concentration, -h′(0) for different values of the similarity parameter, β and different values of the Schmidt number, Sc in the case of Marangoni parameter, ε=0,1, and −1. Further, the Prandtl number, Pr, is taken to be Pr=0.7 corresponding to air and the Schmidt number; Sc has the following values: Sc=0.22 (hydrogen), 0.6 (water), 0.75 (oxygen), and 0.78 (ammonia). It should be noticed that the results given in the parentheses ( ) are the second (dual) solutions. Results obtained by Pop et al. [21] for the case of impermeable surface are also included in these tables. It is clearly seen that the results agree well in all the three tables.Table 1
Values off′(0), -θ′(0), and -h′(0) for ε=0 (absent Marangoni effect) and f0=0 (impermeable surface) for different values of β when Pr=0.7.
Pop et al. [21]Presentβf0f′(0)-θ′(0)-h′(0)f′(0)-θ′(0)-h′(0)0.1−120.42222 (7.52155)2.73049 (1.39352)1.33457 (0.55567)2.50450 (1.25965)2.83556 (1.45391)2.89636 (1.48832)05.63217 (1.13296)1.66616 (0.63301)0.81816 (0.34081)1.52112 (0.57999)1.73479 (0.65856)1.77485 (0.67361)12.18587 (0.32671)1.51789 (0.81825)0.69038 (0.36368)1.36228 (0.72258)1.59367 (0.86650)1.63858 (0.89559)0.2−16.88979 (3.22734)1.70167 (0.89556)0.91512 (0.42147)1.57757 (0.81507)1.75892 (0.93265)1.79192 (0.95395)03.14160 (0.96668)1.40729 (0.58243)0.72749 (0.30936)1.29013 (0.53105)1.46279 (0.60744)1.49521 (0.62224)11.87556 (0.23194)1.54931 (0.65695)0.72307 (0.29068)1.39337 (0.57422)1.62523 (0.69967)1.67021 (0.72573)0.3−14.09080 (2.49526)1.38372 (0.82414)0.78472 (0.39756)1.29037 (0.75041)1.42662 (0.85849)1.45130 (0.87834)02.44523 (0.97452)1.34800 (0.54113)0.71088 (0.26932)1.23812 (0.48650)1.40007 (0.56816)1.43051 (0.58430)11.73727 (0.24512)1.59029 (0.43442)0.75025 (0.20410)1.43243 (0.37519)1.66711 (0.46669)1.71261 (0.48694)0.4−13.09833 (2.2474)1.26054 (0.14478)0.73410 (0.36768)1.17886 (0.56070)1.29804 (7.84348)1.31959 (1.83203)02.14605 (1.08908)1.34240 (0.35840)0.71485 (0.22042)1.23425 (0.35627)1.39369 (0.34761)1.42366 (0.33547)11.65874 (0.34764)1.62937 (−1.28801)0.77466 (0.07600)1.46934 (−0.48622)1.70719 (−2.35420)1.75326 (−3.77805)0.5−12.639241.212010.715101.134971.247371.2677001.984261.354160.724961.245831.405531.435561.984271.354170.725061.245831.405531.4355611.607811.664440.796021.502331.743221.789861.0−11.980271.211940.722701.136091.246831.2669201.696871.442330.778941.328691.496241.527741.696871.442330.777861.328661.496231.5277411.494251.788010.867681.618111.870421.919152.0−11.742131.291450.770201.210091.328981.3506201.565951.550620.842321.429441.608071.641661.565951.550610.840971.429421.608071.6416611.428711.906940.937621.729181.993032.043905.0−11.616211.386580.824981.298401.427341.4508701.488561.659440.904151.530471.720581.756321.488571.659430.903221.530451.720571.7563111.384782.019461.003091.834072.109122.16207Sc0.220.600.750.780.220.600.750.78Table 2
Values off′(0), -θ′(0), and -h′(0) for ε=1 and f0=0 (impermeable surface) for different values of β when Pr=0.7.
Pop et al. [21]Presentβf0f′(0)-θ′(0)-h′(0)f′(0)-θ′(0)-h′(0)0.1−121.97653 (10.27675)2.83813 (1.70953)1.38014 (0.66173)2.60161 (1.54850)2.94816 (1.78238)3.01186 (1.82402)07.19144 (2.24508)1.85863 (0.80540)0.89765 (0.36166)1.69425 (0.71946)1.93639 (0.84714)1.98180 (0.87179)13.12331 (0.85654)1.67897 (0.87403)0.75490 (0.35695)1.50685 (0.76213)1.76256 (0.93064)1.81201 (0.96479)0.2−18.28421 (4.98050)1.87820 (1.20322)0.98806 (0.50472)1.73707 (1.08332)1.94347 (1.25792)1.98114 (1.28916)04.29569 (2.03295)1.61240 (0.81721)0.80832 (0.35273)1.47475 (0.72457)1.67759 (0.86250)1.71568 (0.88933)12.61520 (0.83925)1.69789 (0.78169)0.77746 (0.30148)1.52666 (0.67010)1.78105 (0.83928)1.83026 (0.87434)0.3−15.32651 (3.83896)1.59751 (1.08646)0.87461 (0.46868)1.48397 (0.97392)1.64991 (1.13936)1.68012 (1.16998)03.39676 (1.97392)1.55082 (0.80880)0.79147 (0.32899)1.42074 (0.70805)1.61248 (0.85897)1.64851 (0.88895)12.37773 (0.89485)1.73187 (0.60914)0.80281 (0.22770)1.55938 (0.50839)1.81562 (0.66417)1.86517 (0.69865)0.4−14.18935 (3.36425)1.48697 (1.27934)0.83076 (0.40355)1.38419 (1.31977)1.53439 (1.29417)1.56172 (1.30698)02.97529 (1.98120)1.53913 (2.00757)0.79355 (0.26050)1.41135 (−2.95626)1.59972 (1.64605)1.63513 (1.55295)12.23941 (1.00247)1.76664 (2.39461)0.82624 (0.08520)1.59238 (10.44973)1.85121 (2.02516)1.90123 (1.90673)0.5−13.618181.439790.817241.341681.485061.5111602.733621.545070.806161.417691.605491.640812.733631.545080.806241.417691.605491.6408212.148371.798880.851031.622831.884301.934811.0−12.701041.425200.821171.329571.469431.4949602.274601.617110.854361.485851.679391.715812.274601.617110.854421.485851.679391.7158111.941891.916520.922211.733212.005332.057812.0−12.329531.491490.861591.391281.537941.5647902.051961.715970.912721.577991.781431.819722.051971.715960.912761.577991.781431.8197211.820822.032780.990581.841822.125172.179745.0−12.123081.578020.911431.471551.627471.6560601.916541.819170.972011.673871.888121.928431.916541.819170.972031.673871.888111.9284311.739022.144101.055231.945592.240032.29666Sc0.220.600.750.780.220.600.750.78Table 3
Values off′(0), -θ′(0), and -h′(0) for ε=-1 and f0=0 (impermeable surface) for different values of β when Pr=0.7.
Pop et al. [21]Presentβf0f′(0)-θ′(0)-h′(0)f′(0)-θ′(0)-h′(0)0.1−11 (18.69400)0.42585 (2.60605)0.31791 (1.28181)0.41446 (2.39223)0.43041 (2.70538)0.43284 (2.76284)01 (3.43343)0.78611 (1.34348)0.44073 (0.68254)0.72780 (1.23088)0.81370 (1.39675)0.82982 (1.42786)11 (0.02861)1.26652 (0.88734)0.58489 (0.41760)1.13660 (0.79323)1.33015 (0.93423)1.36796 (0.96233)0.2−11 (4.97600)0.52173 (1.42764)0.37065 (0.79313)0.50312 (1.32976)0.52969 (1.47255)0.53409 (1.49836)01 (1.00014)0.87518 (0.87524)0.49066 (0.49068)0.81026 (0.81033)0.90590 (0.90596)0.92384 (0.92390)11 (−0.28079)1.34283 (0.63887)0.63095 (0.31836)1.20818 (0.57116)1.40866 (0.67331)1.44774 (0.69417)0.3−11 (1.91347)0.59576 (0.90172)0.41163 (0.56117)0.57163 (0.85234)0.60631 (0.92399)0.61223 (0.93668)01 (0.24145)0.94510 (0.56101)0.52985 (0.35365)0.87499 (0.52701)0.97827 (0.57702)0.99765 (0.58635)11 (−0.40926)1.40391 (0.34253)0.66752 (0.20927)1.26538 (0.31198)1.47153 (0.35879)1.51165 (0.36889)0.4−11 (0.69584)0.65507 (0.51595)0.44459 (0.37337)0.62654 (0.49813)0.66768 (0.52364)0.67480 (0.52792)01 (−0.08927)1.00170 (0.25272)0.56157 (0.21647)0.92739 (0.24983)1.03686 (0.25371)1.05739 (0.25419)11 (−0.43782)1.45403 (−0.21212)0.69734 (0.07165)1.31228 (−0.12659)1.52315 (−0.26250)1.56413 (−0.29559)0.5−110.704050.480980.672130.718300.72637011.048660.595090.970981.085441.1069311.048660.595220.970981.085441.10692111.495980.727521.351541.566361.608071.0−110.859500.564150.816030.879160.89038011.199770.676911.110811.241871.2664611.199760.677001.110811.241861.26645111.633390.806291.479841.708001.752172.0−110.997280.639360.943671.021701.03569011.335130.751171.236111.381981.4093511.335120.751221.236111.381981.40935111.758850.878241.596821.837431.883895.0−111.122210.708341.059451.150931.16742011.458650.819451.350461.509851.5397511.458650.819491.350451.509841.53974111.874890.944731.704891.957192.00582Sc0.220.600.750.780.220.600.750.78Figures1, 2, and 3 present the variations of the surface velocity, temperature, and concentration with β, respectively, when Pr=0.7,Sc=0.78 in the case of ε=0 (Marangoni effect neglected) with the effect of the parameter f0. The dashed line refers to the second solution of the problem. From these figures, we can see that the second solutions exist for β<0.5. This discovery is consistent with the statement given in Pop et al. [21], that for β<0.5, the solutions are not unique. Not only that it can be observed from the figures but the effects of suction or injection parameter can be clearly seen when β>0.4. The imposition of suction (f0>0) at the surface has the tendency to reduce the velocity but increase the temperature and concentration gradients along the interface. Meanwhile, the opposite results are observed for the case of surface injection (f0<0), an increase in the surface velocity and decrease in the surface temperature and concentration gradients.  Hence, consequently, the interface heat transfer and the interface mass transfer will decrease as well.Figure 1
Variation of the surface velocity withβ when Pr=0.7, Sc=0.78, and ε=0.Figure 2
Variation of the surface temperature withβ when Pr=0.7, Sc=0.78, and ε=0.Figure 3
Variation of the surface concentration withβ when Pr=0.7, Sc=0.78, and ε=0.Figures4 to 12 illustrate the effects of parameter f0 on the velocity f′(η), temperature θ(η), and concentration h(η) profiles when Pr=0.7 (air), Sc=0.78 (ammonia), and β=0.35 in the case of ε=0,1, and −1. The influence of f0 on the velocity profiles is depicted in Figures 4, 5, and 6. The figures show that the suction parameter decreases the velocity profiles, and the injection parameter increases the velocity profiles for the case of the Marangoni parameter ε=0 and 1. Further, in Figures 7, 8, and 9, we can see the effects of f0 on the temperature profiles. The figures indicate that for ε=0,1, and −1, imposition of the suction parameter tends to reduce the temperature profiles, whereas the injection parameter increases the profiles. It should be noticed that the second solutions in these cases have positive values. Besides, it can also be seen that reduction in Marangoni parameter ε causes the dual solutions to be more stable. On the other hand, the effects of suction or injection parameter on concentration profiles h(η) are shown in Figures 10, 11, and 12 and are similar to the concentration profiles.Figure 4
Effects off0 on the velocity profiles when Pr = 0.7, Sc = 0.78, β = 0.35, and ε=0.Figure 5
Effects off0 on the velocity profiles when Pr = 0.7, Sc = 0.78, β = 0.35, and ε=1.Figure 6
Effects off0 on the velocity profiles when Pr = 0.7, Sc = 0.78, β = 0.35, and ε=-1.Figure 7
Effects off0 on the temperature profiles when Pr = 0.7, Sc = 0.78, β = 0.35, and ε=0.Figure 8
Effects off0 on the temperature profiles when Pr = 0.7, Sc = 0.78, β = 0.35, and ε=1.Figure 9
Effects off0 on the temperature profiles when Pr = 0.7, Sc = 0.78, β = 0.35, and ε=-1.Figure 10
Effects off0 on the concentration profiles when Pr = 0.7, Sc = 0.78, β = 0.35, and ε=0.Figure 11
Effects off0 on the concentration profiles when Pr = 0.7, Sc = 0.78, β = 0.35, and ε=1.Figure 12
Effects off0 on the concentration profiles when Pr = 0.7, Sc = 0.78, β = 0.35, and ε=-1.


## 4. Conclusions

The present work deals with the thermosolutal Marangoni forced convection boundary layer flow as considered by Pop et al. [21]. We extended the paper by taking into consideration the effects of suction or injection. On the other hand, we also studied the existence of dual similarity solutions in the present problem. Further, the governing equations are transformed into ordinary differential equations and are then solved numerically using the shooting method. The effects of suction or injection parameter on the flow and heat transfer characteristics are studied. In general, imposition of suction is to decrease the velocity, temperature, and concentration profiles, whereas injection shows the opposite effects. On the other hand, the second (dual) solutions are discovered to exist when similarity parameter β<0.5.



---

*Source: 875754-2011-03-09.xml*