# Effect of Wind Stress on the Dynamics and Stability of  Nonisothermal  Power-Law Film  down an Inclined Plane

Authors: B. Uma
Journal: ISRN Mathematical Physics
 (2012)

---

## Abstract

Dynamics and stability of a nonisothermal power-law liquid film down an inclined plane is considered in the presence of interfacial shear. Linear stability characteristics of the power-law liquid film using normal mode approach reveal that isothermal and evaporating films are unstable for any value of power-law index while there exists a critical value of power-law index for the case of condensate film above which condensate film ow system is always stable. This critical value of power-law index increases with the increase in shear stress at the interface. Weakly nonlinear stability analysis using method of multiple scales divulges the existence of zones due to supercritical stability and subcritical instability. The nonlinear evolution equation is solved numerically in a periodic domain. The results reveal that (1) for an isothermal dilatant (pseudoplastic) liquids, the maximum wave amplitude is always smaller (larger) than that for a Newtonian liquid and the amplitude of permanent wave increases with the increase in interfacial shear; (2) condensation of pseudoplastic film happens for the earlier instant of time when the phase change parameter increases and the effect of interfacial shear makes the film more corrugated; (3) dilatant (pseudoplastic) evaporating liquid film attains rupture faster (slower) than that of Newtonian liquid film, and the interfacial shear does not influence the time at which rupture occurs.


---

## Body

## 1. Introduction

Gravity-driven flow of a thin film down a vertical or an inclined plane has attracted much attention due to its importance in many industrial applications such as film coating and interface heat and mass transfer processes in chemical technology and energetics.Yih [1] and Benjamin [2] first studied the linear theory for the isothermal falling film. Yih [3] formulated the problem in terms of long-wave asymptotics and determined the critical Reynolds number above which the instability would occur. Benney [4] extended the theory to nonlinear regime by deriving a nonlinear evolution equation. There have been number of extensions of this work as discussed by Chang and Demekhin [5] in their monograph.Although the theory of laminar film condensation flow due to gravity has been analyzed by Nusselt [6], the stability analysis of condensate or evaporating liquid film down a vertical or an inclined plane has been studied only after 1970s [7–9]. Burelbach et al. [10] and Bankoff [11] formulated a one-sided model and studied nonlinear stability and breakdown of evaporating/condensing horizontal static liquid film. Long-wave instabilities of heated falling films have been investigated by Bankoff [12] and Joo et al. [13]. Effect of interfacial phase change on the linear stability of thin films has been considered by Spindler et al. [14], Ünsal and Thomas [15], and Spindler [16]. They have showed that vapor recoil effects are destabilizing for a volatile liquid. Using perturbation methods, Ünsal and Thomas [17] have investigated the nonlinear stability of vertical condensate film flow. Hwang and Weng [18] have examined the finite-amplitude stability analysis of liquid film down a vertical plane with and without interfacial phase change and have shown that both supercritical stability and subcritical instability are possible for condensate film flow system.These results for a Newtonian film flows (isothermal, condensate, evaporating) cannot completely describe the dynamics of non-Newtonian film flows. For example, mud flows are pseudoplastic, saccharified honey is known to be a shear thickening fluid, whereas polymers exhibit generally a large elastic component. Non-Newtonian fluids generally exhibit a nonlinear relationship between shear stress and shear rate. These flows may be classified as inelastic and viscoelastic. The inelastic fluids may be subdivided as time-dependent fluids (thixotropic and rheopectic) and time-independent fluids (pseudoplastic, dilatant, Bingham plastic, pseudoplastic with yield stress). Apart from investigation on flow characteristics and stability analysis of viscoelastic fluid films along a vertical or on inclined plane [19–27], inelastic time-independent non-Newtonian fluids have received the greatest attention from rheologists [28–33] which has resulted in the development of a number of equations or models proposed to represent their flow behavior. The Ostwald de Waele power-law model represents several inelastic time-independent non-Newtonian fluids of practical interest and therefore has been used in this paper.The investigations on the stability characteristics of Newtonian fluid film down an inclined/a vertical plane show that in the linear theory, the film flow system is unstable for any Reynolds number. However, the finite-amplitude stability analysis of liquid films down a vertical wall by Hwang and Weng [18] with interfacial phase change reveals that the isothermal and evaporating Newtonian films are unstable for any Reynolds number while there exists a finite critical Reynolds number for the case of condensate film below which condensate Newtonian film flow system is always stable. This shows that the effect of mass transfer at the interface of a Newtonian fluid film strongly modifies the stability characteristics of the film flow when the phase change is considered. Usha and Uma [34] have extended their study to condensate/evaporating power-law liquid film down an inclined plane with or without interfacial phase change to show that isothermal and evaporating films are unstable for any value of power-law index “n” while there exists a critical value of power-law index “n” for the case of condensate film above which condensate film flow system is always stable. van der Walls interactions are taken into account by Gorla [35] to study the rupture of thin power-law liquid film on a cylinder. He has shown that rupture time for the dilatant fluids are higher than that of Newtonian and pseudoplastic fluids. This shows that a detailed numerical investigation on nonlinear evolution of nonisothermal power-law liquid film down an inclined plane is required to shed more light on understanding the effect of power-law index “n” on different kinds of permanent wave shapes (isothermal film), film condensation (condensate film), and film rupture (evaporating film).Recently, Pascal and D’Alessio [36] have studied the generation and structure of roll waves on the surface of a isothermal power-law liquid film down an inclined plane. To have a more realistic description of flows taking place in an environment, they have considered the effect of wind stress acting on the surface of a power-law film down an incline. Effect of interfacial-induced shear on thin films has been first attempted by Sheintuch and Dukler [37]. Later, Pascal [38] studied the effect of superficial shear stress in connection with the study of wind-aided spreading of oil on the sea.The chief motivation of the present study is to investigate the dynamics and stability of more realistic flows taking place in an environment. In this case, it becomes important to include the effects of superficial shear stress and phase change at the interface of a fluid film flowing down an inclined plane. The Ostwald de Waele power-law model is considered to represent the non-Newtonian fluid.In this paper, dynamics and stability of condensate/evaporating power-law liquid film flowing down an inclined plane with the effect of wind stress and phase change at the interface are considered. Linear stability analysis using normal mode approach reveals that cocurrent superficial wind stress at the interface destabilizes the film flow system. Weakly nonlinear stability analysis using method of multiple scales divulges the existence of zones due to supercritical stability and subcritical instability. The nonlinear evolution equation describing the shape of the free surface has been solved numerically in a periodic domain using Lee’s three-time level method. Applications of the results indicate that the amplitude of permanent wave increases with the increase in prescribed shear stress at the interface for isothermal power-law film and pseudoplastic liquid film has the largest amplitude compared to Newtonian and dilatant film. Also, evaporating liquid film attains rupture faster than that of Newtonian and condensate liquid film, and the prescribed shear stress at the interface does not influence the time at which rupture occurs.


## 2. Mathematical Formulation

A thin power-law liquid film flowing down an inclined planey¯=0 with the effect of prescribed superficial shear stress and phase change acting on the surface of a fluid layer y¯=h¯(x¯,t¯) (Figure 1) is considered.Figure 1
Schematic representation of a thin film down an inclined plane.The governing equations [18, 24, 34] are the two-dimensional mass, momentum, and energy balance equations for the power-law model given by(2.1)∂u¯∂x¯+∂v¯∂y¯=0,ρ(∂u¯∂t¯+u¯∂u¯∂x¯+v¯∂u¯∂y¯)=-∂p¯∂x¯+ρgsin⁡⁡θ(1-γ)+∂τ¯x¯x¯∂x¯+∂τ¯x¯y¯∂y¯,ρ(∂v¯∂t¯+u¯∂v¯∂x¯+v¯∂v¯∂y¯)=-∂p¯∂y¯-ρgcos⁡⁡θ(1-γ)+∂τ¯y¯x¯∂x¯+∂τ¯y¯y¯∂y¯,(∂T¯∂t¯+u¯∂T¯∂x¯+v¯∂T¯∂y¯)=Kρcp(∂2T¯∂x¯2+∂2T¯∂y¯2),
where K is the thermal conductivity, ρ is the density, cp is the liquid specific heat, g is the gravity, γ is the ratio of vapor density to liquid density, and T¯ is the temperature.The boundary conditions at the wall are the no-slip condition of velocity and a constant wall temperature (T¯w) given by(2.2)u¯=0,v¯=0,T¯=T¯won⁡y¯=0.
The effect of wind stress is implemented into the model by prescribing the shear stress τ¯w at the free surface of the fluid layer [36]. The boundary conditions at the liquid-vapor interface [17, 18, 34, 36] are the balance of normal and tangential stresses, the relation of interfacial energy balances, and the equality of liquid and saturated vapor temperatures (T¯s) and are given by(2.3)p¯+(2∂h¯∂x¯τ¯x¯y¯-(∂h¯∂x¯)2τ¯x¯x¯-τ¯y¯y¯)(1+(∂h¯∂x¯)2)-1+K2hfg-2(γ-1)ργ(∂T¯∂y¯-∂h¯∂x¯∂T¯∂x¯)2×(1+(∂h¯∂x¯)2)-1+σ∂2h¯∂x¯2[1+(∂h¯∂x¯)2]-3/2=p¯aon⁡y¯=h¯,∂h¯∂x¯(τ¯y¯y¯-τ¯x¯x¯)+(1-(∂h¯∂x¯)2)τ¯x¯y¯=τ¯won⁡y¯=h¯,K(∂T¯∂y¯-∂h¯∂x¯∂T¯∂x¯)-ρhfg(∂h¯∂t¯+u¯∂h¯∂x¯-v¯)=0on⁡y¯=h¯,T¯=T¯son⁡y¯=h¯,
where(2.4)τ¯x¯x¯=2μn∂u¯∂x¯(∂u¯∂y¯)n-1,τ¯x¯y¯=τ¯y¯x¯=μn(∂u¯∂y¯)n,τ¯y¯y¯=-2μn∂u¯∂x¯(∂u¯∂y¯)n-1,μn is the consistency coefficient, and n is the flow index. When the power-law exponent n is equal to 1, then the model describes the Newtonian fluid; If n<1, the fluid is said to be pseudoplastic or shear thinning, and, if n>1, the fluid is called dilatant or shear thickening. Using the dimensionless quantities defined by(2.5)u=u¯u¯0,v=v¯κu¯0(velocitiesinthexandydirectionsresp.),x=κx¯h¯0,y=y¯h¯0,t=κu¯0t¯h¯0(time),H=h¯h¯0(filmthickness),Fr=u¯02g(1-γ)h¯0(Froudenumber),u¯0=nn+1(ρg(1-γ)sin⁡⁡θμn)1/nh¯0(n+1)/n(referencevelocity),κ=2πh¯0λ(wavenumber),p=p¯-pa¯ρu¯02(pressure),Θ=T¯-T¯wT¯s-T¯w(temperature),Ren=ρu¯02-nh¯0nμn(Reynoldsnumber),Prn=μnu¯0n-1h¯01-ncpK(Prandtlnumber),Pen=Prn⋅Ren(Pecletnumber),ΔT¯=T¯s-T¯w,ξ=cpΔT¯hfg(hfg-latentheat),Fin=σρu¯02h¯0=WenRen-(3n+2)/(n+2)(1sin⁡⁡θ)(3n-2)/(n+2),Nd=ξ2γPrn2,τ=τ¯wρgh¯0sin⁡⁡θ(Shearstressattheinterface),Wen=σρ-1[(n+1n)n(3n-2)[(1-γ)g]-(3n-2)ν-4]1/(n+2)(Webernumber),
the nondimensional governing equations and the boundary conditions are obtained as(2.6)∂u∂x+∂v∂y=0,(2.7)κ(∂u∂t+u∂u∂x+v∂u∂y+∂p∂x)=sin⁡⁡θFr+κ2Ren∂τxx∂x+1Ren∂τxy∂y,(2.8)κ2(∂v∂t+u∂v∂x+v∂v∂y)=-∂p∂y-cos⁡⁡θFr+κRen(∂τyx∂x+∂τyy∂y),(2.9)∂2Θ∂y2=κPen(∂Θ∂t+u∂Θ∂x+v∂Θ∂y)-κ2∂2Θ∂x2,(2.10)u=0,v=0,Θ=0on⁡y=0,(2.11)p+1Ren[2κ∂H∂xτxy-κτyy-κ3(∂H∂x)2τxx](1+κ2(∂H∂x)2)-1+κ2Fin∂2H∂x2×(1+κ2(∂H∂x)2)-3/2+NdRen-2(γ-1)(∂Θ∂y-κ2∂H∂x∂Θ∂x)2(1+κ2(∂H∂x)2)-1=0on⁡y=H,(2.12)κ2∂H∂x(τyy-τxx)+(1-κ2(∂H∂x)2)τxy=Rensin⁡⁡θFrτon⁡y=H,(2.13)ξ(∂Θ∂y-κ2∂H∂x∂Θ∂x)-κPen(∂H∂t+u∂H∂x-v)=0on⁡y=H,(2.14)Θ=1on⁡y=H.
In (2.13), ξ>0 corresponds to condensate film flow (saturated vapor temperature is more than the wall temperature; (T¯s>T¯w)), ξ<0 corresponds to evaporating film flow (saturated vapor temperature is less than the wall temperature; (T¯s<T¯w)), and ξ=0 corresponds to isothermal film (saturated vapor temperature is equal to the wall temperature; (T¯s=T¯w)).It is to be noted that, whenn=1,τ=0,θ=90°, the above equations and boundary conditions reduce to evaporating or condensating Newtonian flow down a vertical wall investigated by Hwang and Weng [18] and, for τ=0, they reduce to the equations obtained for nonisothermal power-law fluid film flow investigated by Usha and Uma [34]. Since the long wavelength modes are the most unstable ones for the film flow, the physical quantities u,v,p, and Θ are expanded in powers of small wave number κ. Substituting these in (2.6)–(2.14) and collecting the coefficients of like powers of κ, the zeroth- and the first-order equations are obtained. Noting that Wen is large in practical applications, κ2Wen is taken to be of order one. Also, since the effect of Nd has been found to be negligible on stability [15], NdRen-2 is taken to be of order κ2. Further, in the analysis, Ren≃O(1) and Pen≃O(1). The first-order solutions, u1,v1, and Θ1, will have the time derivative Ht in the right hand side, which is replaced by the space derivative from the leading order representation [18] given by(2.15)Ht=ξκPenΘ0y-u0Hx+v0.
The zeroth- and the first-order (after replacing all the time derivatives with (2.15)) solutions are presented in the appendix.The solutions (A.1) and (A.2) are substituted into the expanded kinematic boundary condition (2.13) yields:(2.16)Ht=ξκPen[Θ0y+κΘ1y-κ2Hx(Θ0x+κΘ1x)]-{(u0+κu1)Hx-v0-κv1}.
The simplification of (2.16) now gives the generalized kinematic equation as(2.17)Ht+X(H)+A(H)Hx+B(H)Hxx+C(H)Hxxxx+D(H)Hx2+E(H)HxHxxx+O(κ2)=0,
where(2.18)X(H)=-ξκPen(1-ξ3)1H,A(H)=(Rensin⁡⁡θFr)1/n×{n2(2n+1)(3n+1)(H+τ)1/nH-ξnn+1×[16(H+τ)(n+1)/n-n2(2n+1)(3n+1)((H+τ)(3n+1)/n+τ(3n+1)/nH2)+2n3(2n+1)(3n+1)(4n+1)((H+τ)(4n+1)/n-τ(4n+1)/nH3)]}+Renξ(n+1)Pen(Rensin⁡⁡θFr)(2-n)/n×{-2n(H+τ)(2-n)/n(H-nτ)+((H+τ)(1-n)/n(n+1)H-n(H+τ)1/n(n+1)H2)×[((H+τ)(n+1)/n(H-nτ)+nτ(2n+1)/n)-n2n+1((H+τ)(2n+1)/n-τ(2n+1)/n)]+n22(n+1)(n+2)1H2[(H+τ)2(n+1)/n-τ2(n+1)/n]-n(n+2)1H(H+τ)(n+2)/n+(H+τ)2/n(H-nτ)H-n2n+1τ(n+1)/nH2[(H+τ)(n+1)/n-τ(n+1)/n]+τ(n+1)/nHn(H+τ)1/n},B(H)=κRen(Rensin⁡⁡θFr)(3-n)/n(H+τ)1/n×{3n2(n+1)(n+2)(2n+1)(H+τ)2(n+1)/nH-n(n+1)(n+2)(H+τ)(n+2)/nH2+(H+τ)2/nH3-3n3(n+1)(n+2)(2n+1)(3n+2)×[(H+τ)(3n+2)/n-τ(3n+2)/n]+2n2(n+1)(2n+1)τ(2n+1)/n×[nn+1((H+τ)(n+1)/n-τ(n+1)/n)-(H+τ)1/nH]+2n2(n+1)(2n+1)(H+τ)1/nH[(H+τ)(2n+1)/n-τ(2n+1)/n]-2n3(n+1)2(2n+1)(H+τ)(n+1)/n[(H+τ)(2n+1)/n-τ(2n+1)/n]+2n2(n+1)2(H+τ)2(n+1)/nH-3nn+1(H+τ)(n+2)/nH2}-κRen(n+1)cos⁡⁡θFr(Rensin⁡⁡θFr)(1-n)/n×{n2(n+1)(2n+1)(H+τ)1/n(H-nτ)H-nn+1[(H+τ)(n+1)/n(H-nτ)+nτ(2n+1)/n]+n2(n+1)(2n+1)[(H+τ)(2n+1)/n-τ(2n+1)/n]},C(H)=κ3Ren(n+1)(Rensin⁡⁡θFr)(1-n)/n×Fin{n2(n+1)(2n+1)(H+τ)1/n(H-nτ)H-nn+1[(H+τ)(n+1)/n(H-nτ)+nτ(2n+1)/n]+n2(n+1)(2n+1)[(H+τ)(2n+1)/n-τ(2n+1)/n]},D(H)=κRen(Rensin⁡⁡θFr)(3-n)/n×{n(6n2+19n+11)(n+1)2(n+2)(2n+1)(H+τ)(n+3)/nH-(7n+15)(n+1)(n+2)(H+τ)3/nH2+3n(H+τ)(3-n)/nH3-4n(n+1)(2n+1)τ(2n+1)/n(H+τ)(2-n)/nH+4n(n+1)(2n+1)(H+τ)(2-n)/nH[(H+τ)(2n+1)/n-τ(2n+1)/n]-3n2(H+τ)(1-n)/n(n+1)(n+2)(2n+1)(3n+2)[(H+τ)(3n+2)/n-τ(3n+2)/n]+2n2(n+1)2(2n+1)τ(2n+1)/n(H+τ)(1-n)/n[(H+τ)(n+1)/n-τ(n+1)/n]-2n2(n+1)2(2n+1)(H+τ)2/n[(H+τ)(2n+1)/n-τ(2n+1)/n]}-κRen(n+1)cos⁡⁡θFr(Rensin⁡⁡θFr)(1-n)/n{1n(H+τ)(1-n)/n(H-nτ)H+(H+τ)1/nH}-ξκPen(1-ξ3)1H,E(H)=κ3Ren(n+1)(Rensin⁡⁡θFr)(1-n)/nFin{1n(H+τ)(1-n)/n(H-nτ)H+(H+τ)1/nH}.


## 3. Stability Analysis

As the variation of the film thickness of the base flow is found to be very small for|κHx|≪1 using an analysis based on Nusselt assumption, the dimensionless film thickness is expressed as(3.1)H=1+η(x,t),
where η(x,t) is the perturbation of the stationary film thickness. The approximation |κHx|≪1 gives qualitative results for the constant film thickness assumption at the zeroth order. It is important to note that this constant film thickness approximation with long-wave perturbations are reasonable approximations only for certain segments of weakly condensing and evaporating flows. Substituting for H(x,t) in (2.18) and retaining terms up to the order of η3, the evolution equation for η is obtained as(3.2)ηt+X′η+Aηx+Bηxx+Cηxxxx=-[X′′η22+X′′′η36+(A′η+A′′η22)ηx+(B′η+B′′η22)ηxx+(C′η+C′′η22)ηxxxx+(D+D′η)ηx2+(E+E′η)ηxηxxx]+O(η4),
where the values of X,A,B,C,D,E, and their derivatives are evaluated at the dimensionless height of the film H=1. It is worth mentioning here that the analysis based on locally valid unsteady equation has been the subject of several investigations of laminar film condensation [7, 8, 15, 17, 18, 34, 39–42]. Equation (3.2) describes the behavior of finite-amplitude disturbances on the power-law film, and it is used to predict the time-wise behavior of an initially sinusoidal disturbance on the power-law film. It is important to remark that the simplified unsteady equation (3.2) is only locally valid.In order to understand the flow characteristics and the associated time-dependent properties of power-law liquid film down an inclined plane, typical values of the physical parameters considered by Lin and Hwang [24], Usha and Uma [34], and Pascal and D’Alessio [36] have been used for numerical evaluation and are given by(3.3)ρ=998Kg/m3,μn=1.002×10-3Pasn,h0=10-4m,σ=0.0727N/m,γ=0.000611,θ=60°.
The temperature at the interface is taken as T¯s=373K and the temperature difference between the wall and the interface as T¯s-T¯w=±47K. Under such temperature conditions, the phase change parameter ξ takes values |ξ|=0.0872 (with phase change) and |ξ|=0 (without phase change). Shear stress parameter τ varies between 0.0 and 2.0.


### 3.1. Linear Stability Analysis

For the linear stability analysis, the nonlinear terms of (3.2) are neglected and the linearized equation(3.4)ηt+X′η+Aηx+Bηxx+Cηxxxx=0
is obtained. Assuming the normal mode solution as(3.5)η=aei(x-dt)+a¯e-i(x-dt),
the complex wave celerity corresponding to linear stability problem is given by(3.6)d=dr+idi=A+i(B-C-X′),
where dr is the linear wave speed and di is the linear growth rate of the amplitudes. The flow is in a linearly unstable supercritical condition for di>0 and in a linearly stable subcritical condition for di<0. For di=0, the flow is neutrally stable. Zero linear growth rate of the perturbation (di=0) gives raise to a cut-off wave number κ=κc(n,Ren,Wen,θ,ξ,Prn,τ).Figure2 shows the neutral stability curve for isothermal, condensate, and evaporating power-law films with the effect of wind stress at the interface. It is observed from Figure 2 that the neutral stability curve plotted from the linear stability analysis separates κ-n plane into two regions depending on the value of phase change parameter ξ and shear stress parameter τ. Region of linear instability (di>0) increases with the increase in the shear stress τ at the interface and decreases with the increase in the phase change parameter ξ. The neutral stability curve for condensate film shows that there exists a critical value for n, above which condensate power-law fluid film is always stable, and this critical value increases with the increase in shear stress τ at the interface (Figure 2(d)). It is also observed that condensate power-law fluid film is more stable than the corresponding isothermal and evaporating power-law fluid film in the absence of wind stress.Neutral stability curves for isothermal (ξ=0), condensate (ξ=0.0872), and evaporating (ξ=-0.0872) films with the effect of wind stress.
(a)
τ=0(b)
τ=0.2(c)
τ=0.4(d)
ξ=0.0872Figure3 shows the temporal growth rate of power-law fluid given by (3.6). Temporal growth rate of the power-law liquid film increases with the increase in shear stress parameter τ and decrease in power-law index n. It is observed that temporal growth rate is lower for a dilatant (n=1.05) condensate (ξ>0) film when τ=0 and higher for a pseudoplastic (n=0.95) evaporating (ξ<0) film when τ>0. The results of the linear stability analysis are in good agreement with those of Lin and Hwang [24] (isothermal power-law liquid film for zero shear stress at the interface) and Usha and Uma [34] (power-law condensate/evaporating film for zero shear stress at the interface).Growth rate for different values of phase change parameter and wind stress parameter.
(a)
ξ=0,τ=0(b)
ξ=0.0872,τ=0(c)
ξ=-0.0872,τ=0(d)
ξ=0,τ=0.2(e)
ξ=0.0872,τ=0.2(f)
ξ=-0.0872,τ=0.2(g)
ξ=0,τ=0.4(h)
ξ=0.0872,τ=0.4(i)
ξ=-0.0872,τ=0.4


### 3.2. Weakly Nonlinear Stability Analysis

As the perturbed wave grows to a finite amplitude, linear stability theory cannot be used to predict the flow behavior accurately. Therefore, in order to examine whether the finite-amplitude disturbance in the linearly stable region causes instability (subcritical instability) and to investigate whether the subsequent nonlinear evolution of disturbances in the linearly unstable region develops into a new equilibrium state with a finite-amplitude (supercritical stability) or grows to be unstable, the nonlinear stability analysis is employed. The nonlinear stability analysis of (3.2) by the method of multiple scales [43] yields(3.7)(L0+ϵL1+ϵ2L2)(ϵη1+ϵ2η2+ϵ3η3)=-ϵ2N2-ϵ3N3+o(ϵ4),
where(3.8)η(ϵ,x,x1,t,t1,t2)=ϵη1+ϵ2η2+ϵ3η3,t1=ϵt,t2=ϵ2t,x1=ϵx,∂∂t⟶∂∂t+ϵ∂∂t1+ϵ2∂∂t2;∂∂x⟶∂∂x+ϵ∂∂x1,L0=∂∂t+X′+A∂∂x+B∂2∂x2+C∂4∂x4,L1=∂∂t1+A∂∂x1+2B∂2∂x∂x1+4C∂4∂x3∂x1,L2=∂∂t2+B∂2∂x12+6C∂4∂x2∂x12,N2=X′′2η12+A′η1η1x+B′η1η1xx+C′η1η1xxxx+Dη1x2+Eη1xη1xxx,N3=X′′η1η2+X′′′6η13+A′(η1η2x+η1xη2+η1η1x1)+B′(η1η2xx+2η1η1xx1+η1xxη2)+C′(η1η2xxxx+4η1η1xxxx1+η1xxxxη2)+D(2η1xη2x+2η1xη1x1)+E(η1xη2xxx+3η1xη1xxx1+η1xxxη2x+η1xxxη1x1)+A′′2η12η1x+B′′2η12η1xx+C′′2η12η1xxxx+D′η1η1x2+E′η1η1xη1xxx.
The solution of (3.7) at the order O(ϵ) is obtained by solving L0η1=0 and is in the form(3.9)η1=aei(x-drt)+a¯e-i(x-drt),
where a(x1,t1,t2) is the nonlinear amplitude function and a¯(x1,t1,t2) is its complex conjugate. The solution of the equation L0η2+L1η1=-N2 at the O(ϵ2) is in the form(3.10)η2=ea2e2i(x-drt)+e¯a¯2e-2i(x-drt).
Using the solutions for η1 and η2 in the O(ϵ3) equation given by L0η3+L1η2+L2η1=-N3, the equation for the perturbation amplitude a(x1,t1,t2) is obtained as(3.11)∂a∂t2+D1∂2a∂x12-ϵ-2dia+(E1+iF1)a2a¯=0
from the secular condition for O(ϵ3), where(3.12)D1=B-6C,F1=(X′′-5B′+17C′+4D-10E)ei+A′er+A′′2,E1=(X′′-5B′+17C′+4D-10E)er-A′ei+(X′′′2-32B′′+32C′′+D′-E′),er+iei=(-(X′′/2)+B′-C′+D-E)-iA′(16C-4B+X′).
The weakly nonlinear behavior of the fluid film can be investigated using (3.12). It is important to note that such an expansion is only valid for wave numbers close to neutral and not near critical when κ approaches zero. The solution of (3.12) for a filtered wave in which spatial modulation does not exist and the diffusion terms in (3.12) vanishes is obtained by taking a=a0e-ib(t2)t2. This leads to the Ginzburg-Landau equation given by(3.13)∂a0∂t2=(ϵ-2di-E1a02)a0,(3.14)∂(b(t2)t2)∂t2=F1a02.
The second term in (3.13) induced by the effect of nonlinearity can either accelerate or decelerate the exponential growth of the linear disturbance depending upon the signs of di and E1. The perturbed wave speed caused by the infinitesimal disturbances appearing in the nonlinear system can be modified using (3.13). The threshold amplitude ϵa0 is given by(3.15)ϵa0=diE1,
and the nonlinear wave speed is given as(3.16)Ncr=ϵ2b=dr+diF1E1.
It is observed from (3.15) that in the linearly unstable region (di>0), the condition for existence of a supercritical stable region is E1>0 and ϵa0 is the threshold amplitude. In the linearly stable region (di<0), if E1<0, then the flow has the behavior of subcritical instability and ϵa0 is the threshold amplitude. The condition for the existence of a subcritical stable region is E1>0, and E1=0 gives the condition of existence of a neutral stability curve.The neutral stability curves are obtained from (3.6) and (3.12) by equating to zero, the linear amplification rate di, and the nonlinear amplification rate E1. Figure 4 shows the regions of subcritical stability (di<0,E1>0) and subcritical instability (di<0,E1<0) in the linearly stable region and supercritical stability (di>0,E1>0) and supercritical explosive state (di>0,E1<0) in the linearly unstable region. It is observed that subcritical stable region (di<0,E1>0) decreases, and supercritical stability (di>0,E1>0) increases with the increase in shear stress parameter τ. However, for the condensate film (ξ>0), region of supercritical explosive state (di>0,E1<0) and supercritical stability (di>0,E1>0) exists only for the pseudoplastic fluids (n<1) which increases with the increase in τ.Neutral stability curves in theκ-n plane for different values of phase change parameter ξ (ξ=0—isothermal film flow; ξ=0.0872—condensate film flow; ξ=-0.0872—evaporating film flow) and wind stress parameters τ.
(a)
ξ=0,τ=0(b)
ξ=0.0872,τ=0(c)
ξ=-0.0872,τ=0(d)
ξ=0,τ=0.2(e)
ξ=0.0872,τ=0.2(f)
ξ=-0.0872,τ=0.2(g)
ξ=0,τ=0.4(h)
ξ=0.0872,τ=0.4(i)
ξ=-0.0872,τ=0.4Figures5 and 6 show the threshold amplitude and nonlinear wave speed in the supercritical stable region for various values of phase change parameter ξ and shear stress parameter τ. The flow system is stable if the finite amplitude of the disturbance is more than the threshold amplitude and explosively unstable otherwise. The results show that the threshold amplitude and nonlinear wave speed increases with increase in τ and decreases with increase in ξ and n.Threshold amplitude in the supercritical stable region.
(a)
ξ=0,τ=0(b)
ξ=0.0872,τ=0(c)
ξ=-0.0872,τ=0(d)
ξ=0,τ=0.2(e)
ξ=0.0872,τ=0.2(f)
ξ=-0.0872,τ=0.2(g)
ξ=0,τ=0.4(h)
ξ=0.0872,τ=0.4(i)
ξ=-0.0872,τ=0.4Nonlinear wave speed in the supercritical stable region.
(a)
ξ=0,τ=0(b)
ξ=0.0872,τ=0(c)
ξ=-0.0872,τ=0(d)
ξ=0,τ=0.2(e)
ξ=0.0872,τ=0.2(f)
ξ=-0.0872,τ=0.2(g)
ξ=0,τ=0.4(h)
ξ=0.0872,τ=0.4(i)
ξ=-0.0872,τ=0.4The above weakly nonlinear analysis shows the existence of both a supercritical stable region and a subcritical unstable region for isothermal and nonisothermal power-law liquid film down an inclined plane in the presence of shear stress at the interface. It is observed that the evolution of two-dimensional waves depends strongly on the initial disturbance wave number. There exists a value ofκs (obtained from equating denominator of E1=0) such that, when κs<κ<κc (obtained from di=0; cut-off wave number), the flow is supercritically stable and nonlinear equilibration occurs after the initial instability. The weakly nonlinear stability analysis shows that κs=κc/2, and the curve κ=κs separates the linearly unstable region into two portions where the nonlinear waves attain a finite equilibrium amplitude (di>0,E1>0) or reach an explosive state (di>0,E1<0), and the flow is supercritically equilibrated.


### 3.3. Nonlinear Analysis

In order to understand the mechanism responsible for the transfer of energy from the basic state to the disturbance, the evolution of finite-amplitude perturbations is considered. On the evolution of these perturbations, the influence of wind stress and mass transfer at the interface of a thin power-law liquid film down an inclined plane is examined by numerically solving the nonlinear evolution equation. The initial disturbance is taken to be a sinusoidal wave with small amplitude given by(3.17)H(x,0)=1-0.1cos⁡⁡(κx).
The evolution of waves with time is obtained by solving the evolution equation (2.18) by using Lee’s three-time level [44] finite-difference method in a periodic domain [-π/κ,π/κ]. The finite-difference scheme achieves linearity in the unknown Hi,j+1 by evaluating all coefficients of Hi,j+1 at a time level of known solution values, preserves stability by averaging Hi,j over three time levels, and maintains accuracy by using central difference approximations, where j denotes the number of time steps and i is a spatial grid point. Special attention is given to the discretization of the wave propagation term so as to overcome the solution blowup for very small initial wave numbers. The scheme possesses both conservative and transportive properties and allows to follow the solution over larger periods of time. The resulting algebraic equations are solved using MATLAB, and the iterations are continued until the local minimum layer thickness becomes smaller than the maximum error bound of 10-10. The computations are performed with Δt=0.01 and with the number of nodes along the spatial direction as N=100 and Δx=2π/Nκ. It is to be noted that the solution obtained with a decrease in Δt or an increase in N does not show any deviation from the values obtained for Δt=0.01, Δx=2π/Nκ, N=100. Further, it is confirmed that the wave that emerges on the film surface does not change its shape for a relatively long period of time.The weakly nonlinear analysis predicts the occurrence of a wave numberκs that separates the regimes of supercritical κs<κ<κc and subcritical κ<κs domains (Figure 4). For a wave number κ close to κc, the solution may evolve into a stable, almost sinusoidal wave of small finite amplitude. For κ close to κm=κc/2, the surface wave approaches the form of a solitary wave as observed in experimental results by Alekseenko et al. [45] and Liu et al. [46]. The evolution of the film in the supercritical stable region is analyzed by performing numerical simulations for initial disturbance wave number κ=κc,κm(=κc/2), and κs(=κc/2). The evolution of the film is examined for values of ξ in the range −0.0872 to 0.0872 and τ between 0 and 2.0.Figures7–9 show the evolution of the free surface configuration of isothermal (ξ=0) pseudoplastic (n=0.95), Newtonian (n=1.0), and dilatant (n=1.05) fluids at various instants of time for different values of shear stress parameter at the interface. Wave number is chosen in the range κs(=κc/2)<κ<κc (Figure 7—κ=0.06, κc=0.0868 when τ=0,κc=0.0971 when τ=0.2, and κc=0.1064 when τ=0.4; Figure 8—κ=0.05,κc=0.0569 when τ=0,κc=0.0635 when τ=0.2, and κc=0.0694 when τ=0.4; Figure 9—κ=0.03, κc=0.0365 when τ=0,κc=0.0409 when τ=0.2, and κc=0.0449 when τ=0.4). It is observed that the growth rate is much more important and the distortion of the free surface is conspicuous with increase in time. The wave amplitude decays after reaching a maximum, and the rear becomes longer and longer to attain a one-hump solitary-like wave at different instant of time for various values of τ. The large amplitude wave is then followed by a small amplitude capillary wave. In all the cases, it is observed that finite-amplitude permanent wave emerges. For a given n, amplitude of the wave increases with the increase in shear stress parameter τ. The corresponding evolution of maximum (Hmax) and minimum (Hmin) thickness is presented in Figure 10. The disturbance amplitude Hmax remains constant for initial times t~0–5 for pseudoplastic fluids (n=0.95), t~0–10 for Newtonian fluids (n=1.0), and t~0–20 for dilatant fluids (n=1.05) and then increases monotonically to attain a maximum amplitude. Isothermal pseudoplastic fluid attains the maximum amplitude in short time compared to Newtonian and dilatant fluid, and this amplitude further increases with the presence of shear stress at the interface. The distortion of the free surface is more significant when the effect of shear stress increases (Figure 11).Evolution of free surface for an isothermal (ξ=0) pseudoplastic (n=0.95) film when κ=0.06 and τ=0,0.2,0.4.
(a)(b)(c)Evolution of free surface for an isothermal (ξ=0) Newtonian (n=1.0) film when κ=0.05 and τ=0,0.2,0.4.
(a)(b)(c)Evolution of free surface for isothermal (ξ=0) dilatant (n=1.05) film when κ=0.03 and τ=0,0.2,0.4.
(a)(b)(c)Maximum and minimum amplitude of the waves down an inclined plane for an isothermal power-law film.
(a)
n=0.95,ξ=0.0(b)
n=1.0,ξ=0.0(c)
n=1.05,ξ=0.0Evolution of free surface for an isothermal film (ξ=0) when τ=2.0.
(a)(b)(c)Weakly nonlinear stability analysis predicted that for the condensate film (ξ>0), region of supercritical stability (di>0,E1>0) exists only for the pseudoplastic fluids (n=0.95). This is the case when the saturated vapor temperature is more than the wall temperature. When condensation is considered, the static layer grows thicker with time. Evolution of condensate pseudoplastic liquid film down an inclined plane is depicted in Figure 12 for ξ=0.01,0.02,0.04,and0.06 and τ=0. The result shows that pseudoplastic film down an inclined plane thickens faster as the phase change parameter ξ increases gradually in the absence of shear stress parameter τ. The effect of wind stress on the interface of the pseudoplastic film takes little longer time to thicken (Figure 13) and makes the film more corrugated as τ increases (Figures 13(a)(ii), 13(b)(ii), 13(c)(ii)).Evolution of free surface for condensating (ξ>0) pseudoplastic (n=0.95) film when κ=0.06 and τ=0.
(a)
n=0.95,ξ=0.01(b)
n=0.95,ξ=0.02(c)
n=0.95,ξ=0.04(d)
n=0.95,ξ=0.06Evolution of free surface for condensating (ξ=0.0872) pseudoplastic (n=0.95) film when κ=0.06 and τ=0,0.2,0.4.
(a)(b)(c)Figure14 shows the evolution of the free-surface shape for evaporating pseudoplastic, Newtonian, and dilatant fluids when ξ=-0.01 and −0.0872 with a zero shear stress at the interface. This is the case when the saturated vapor temperature is less than the wall temperature. It is observed that as n increases, depth of the liquid layer disappears in a finite time t=tr decreases leading to rupture instability. The rupture instability is thus augmented by the evaporative effect. Plot describes that the dilatant film evaporates faster than the pseudoplastic and Newtonian film. For a given n (given pseudoplastic, Newtonian, or dilatant fluid), evaporation of the film happens at earlier instant of time when the phase change parameter ξ decreases (Figure 14). Presence of wind stress on the interface does not influence the time of rupture tr for any n whereas it has little influence on the position where the rupture takes place (Figure 15).Evolution of free surface for an evaporating film (ξ<0) when κ=0.06 and τ=0.
(a)
n=0.95,ξ=-0.01(b)
n=1.0,ξ=-0.01(c)
n=1.05,ξ=-0.01(d)
n=0.95,ξ=-0.0872(e)
n=1.0,ξ=-0.0872(f)
n=1.05,ξ=-0.0872Evolution of free surface for an evaporating film (ξ=-0.00872) when κ=0.06 and τ=0.2,2.0.
(a)
n=0.95,τ=0.2(b)
n=1.0,τ=0.2(c)
n=1.05,τ=0.2(d)
n=0.95,τ=2.0(e)
n=1.0,τ=2.0(f)
n=1.05,τ=2.0


## 3.1. Linear Stability Analysis

For the linear stability analysis, the nonlinear terms of (3.2) are neglected and the linearized equation(3.4)ηt+X′η+Aηx+Bηxx+Cηxxxx=0
is obtained. Assuming the normal mode solution as(3.5)η=aei(x-dt)+a¯e-i(x-dt),
the complex wave celerity corresponding to linear stability problem is given by(3.6)d=dr+idi=A+i(B-C-X′),
where dr is the linear wave speed and di is the linear growth rate of the amplitudes. The flow is in a linearly unstable supercritical condition for di>0 and in a linearly stable subcritical condition for di<0. For di=0, the flow is neutrally stable. Zero linear growth rate of the perturbation (di=0) gives raise to a cut-off wave number κ=κc(n,Ren,Wen,θ,ξ,Prn,τ).Figure2 shows the neutral stability curve for isothermal, condensate, and evaporating power-law films with the effect of wind stress at the interface. It is observed from Figure 2 that the neutral stability curve plotted from the linear stability analysis separates κ-n plane into two regions depending on the value of phase change parameter ξ and shear stress parameter τ. Region of linear instability (di>0) increases with the increase in the shear stress τ at the interface and decreases with the increase in the phase change parameter ξ. The neutral stability curve for condensate film shows that there exists a critical value for n, above which condensate power-law fluid film is always stable, and this critical value increases with the increase in shear stress τ at the interface (Figure 2(d)). It is also observed that condensate power-law fluid film is more stable than the corresponding isothermal and evaporating power-law fluid film in the absence of wind stress.Neutral stability curves for isothermal (ξ=0), condensate (ξ=0.0872), and evaporating (ξ=-0.0872) films with the effect of wind stress.
(a)
τ=0(b)
τ=0.2(c)
τ=0.4(d)
ξ=0.0872Figure3 shows the temporal growth rate of power-law fluid given by (3.6). Temporal growth rate of the power-law liquid film increases with the increase in shear stress parameter τ and decrease in power-law index n. It is observed that temporal growth rate is lower for a dilatant (n=1.05) condensate (ξ>0) film when τ=0 and higher for a pseudoplastic (n=0.95) evaporating (ξ<0) film when τ>0. The results of the linear stability analysis are in good agreement with those of Lin and Hwang [24] (isothermal power-law liquid film for zero shear stress at the interface) and Usha and Uma [34] (power-law condensate/evaporating film for zero shear stress at the interface).Growth rate for different values of phase change parameter and wind stress parameter.
(a)
ξ=0,τ=0(b)
ξ=0.0872,τ=0(c)
ξ=-0.0872,τ=0(d)
ξ=0,τ=0.2(e)
ξ=0.0872,τ=0.2(f)
ξ=-0.0872,τ=0.2(g)
ξ=0,τ=0.4(h)
ξ=0.0872,τ=0.4(i)
ξ=-0.0872,τ=0.4


## 3.2. Weakly Nonlinear Stability Analysis

As the perturbed wave grows to a finite amplitude, linear stability theory cannot be used to predict the flow behavior accurately. Therefore, in order to examine whether the finite-amplitude disturbance in the linearly stable region causes instability (subcritical instability) and to investigate whether the subsequent nonlinear evolution of disturbances in the linearly unstable region develops into a new equilibrium state with a finite-amplitude (supercritical stability) or grows to be unstable, the nonlinear stability analysis is employed. The nonlinear stability analysis of (3.2) by the method of multiple scales [43] yields(3.7)(L0+ϵL1+ϵ2L2)(ϵη1+ϵ2η2+ϵ3η3)=-ϵ2N2-ϵ3N3+o(ϵ4),
where(3.8)η(ϵ,x,x1,t,t1,t2)=ϵη1+ϵ2η2+ϵ3η3,t1=ϵt,t2=ϵ2t,x1=ϵx,∂∂t⟶∂∂t+ϵ∂∂t1+ϵ2∂∂t2;∂∂x⟶∂∂x+ϵ∂∂x1,L0=∂∂t+X′+A∂∂x+B∂2∂x2+C∂4∂x4,L1=∂∂t1+A∂∂x1+2B∂2∂x∂x1+4C∂4∂x3∂x1,L2=∂∂t2+B∂2∂x12+6C∂4∂x2∂x12,N2=X′′2η12+A′η1η1x+B′η1η1xx+C′η1η1xxxx+Dη1x2+Eη1xη1xxx,N3=X′′η1η2+X′′′6η13+A′(η1η2x+η1xη2+η1η1x1)+B′(η1η2xx+2η1η1xx1+η1xxη2)+C′(η1η2xxxx+4η1η1xxxx1+η1xxxxη2)+D(2η1xη2x+2η1xη1x1)+E(η1xη2xxx+3η1xη1xxx1+η1xxxη2x+η1xxxη1x1)+A′′2η12η1x+B′′2η12η1xx+C′′2η12η1xxxx+D′η1η1x2+E′η1η1xη1xxx.
The solution of (3.7) at the order O(ϵ) is obtained by solving L0η1=0 and is in the form(3.9)η1=aei(x-drt)+a¯e-i(x-drt),
where a(x1,t1,t2) is the nonlinear amplitude function and a¯(x1,t1,t2) is its complex conjugate. The solution of the equation L0η2+L1η1=-N2 at the O(ϵ2) is in the form(3.10)η2=ea2e2i(x-drt)+e¯a¯2e-2i(x-drt).
Using the solutions for η1 and η2 in the O(ϵ3) equation given by L0η3+L1η2+L2η1=-N3, the equation for the perturbation amplitude a(x1,t1,t2) is obtained as(3.11)∂a∂t2+D1∂2a∂x12-ϵ-2dia+(E1+iF1)a2a¯=0
from the secular condition for O(ϵ3), where(3.12)D1=B-6C,F1=(X′′-5B′+17C′+4D-10E)ei+A′er+A′′2,E1=(X′′-5B′+17C′+4D-10E)er-A′ei+(X′′′2-32B′′+32C′′+D′-E′),er+iei=(-(X′′/2)+B′-C′+D-E)-iA′(16C-4B+X′).
The weakly nonlinear behavior of the fluid film can be investigated using (3.12). It is important to note that such an expansion is only valid for wave numbers close to neutral and not near critical when κ approaches zero. The solution of (3.12) for a filtered wave in which spatial modulation does not exist and the diffusion terms in (3.12) vanishes is obtained by taking a=a0e-ib(t2)t2. This leads to the Ginzburg-Landau equation given by(3.13)∂a0∂t2=(ϵ-2di-E1a02)a0,(3.14)∂(b(t2)t2)∂t2=F1a02.
The second term in (3.13) induced by the effect of nonlinearity can either accelerate or decelerate the exponential growth of the linear disturbance depending upon the signs of di and E1. The perturbed wave speed caused by the infinitesimal disturbances appearing in the nonlinear system can be modified using (3.13). The threshold amplitude ϵa0 is given by(3.15)ϵa0=diE1,
and the nonlinear wave speed is given as(3.16)Ncr=ϵ2b=dr+diF1E1.
It is observed from (3.15) that in the linearly unstable region (di>0), the condition for existence of a supercritical stable region is E1>0 and ϵa0 is the threshold amplitude. In the linearly stable region (di<0), if E1<0, then the flow has the behavior of subcritical instability and ϵa0 is the threshold amplitude. The condition for the existence of a subcritical stable region is E1>0, and E1=0 gives the condition of existence of a neutral stability curve.The neutral stability curves are obtained from (3.6) and (3.12) by equating to zero, the linear amplification rate di, and the nonlinear amplification rate E1. Figure 4 shows the regions of subcritical stability (di<0,E1>0) and subcritical instability (di<0,E1<0) in the linearly stable region and supercritical stability (di>0,E1>0) and supercritical explosive state (di>0,E1<0) in the linearly unstable region. It is observed that subcritical stable region (di<0,E1>0) decreases, and supercritical stability (di>0,E1>0) increases with the increase in shear stress parameter τ. However, for the condensate film (ξ>0), region of supercritical explosive state (di>0,E1<0) and supercritical stability (di>0,E1>0) exists only for the pseudoplastic fluids (n<1) which increases with the increase in τ.Neutral stability curves in theκ-n plane for different values of phase change parameter ξ (ξ=0—isothermal film flow; ξ=0.0872—condensate film flow; ξ=-0.0872—evaporating film flow) and wind stress parameters τ.
(a)
ξ=0,τ=0(b)
ξ=0.0872,τ=0(c)
ξ=-0.0872,τ=0(d)
ξ=0,τ=0.2(e)
ξ=0.0872,τ=0.2(f)
ξ=-0.0872,τ=0.2(g)
ξ=0,τ=0.4(h)
ξ=0.0872,τ=0.4(i)
ξ=-0.0872,τ=0.4Figures5 and 6 show the threshold amplitude and nonlinear wave speed in the supercritical stable region for various values of phase change parameter ξ and shear stress parameter τ. The flow system is stable if the finite amplitude of the disturbance is more than the threshold amplitude and explosively unstable otherwise. The results show that the threshold amplitude and nonlinear wave speed increases with increase in τ and decreases with increase in ξ and n.Threshold amplitude in the supercritical stable region.
(a)
ξ=0,τ=0(b)
ξ=0.0872,τ=0(c)
ξ=-0.0872,τ=0(d)
ξ=0,τ=0.2(e)
ξ=0.0872,τ=0.2(f)
ξ=-0.0872,τ=0.2(g)
ξ=0,τ=0.4(h)
ξ=0.0872,τ=0.4(i)
ξ=-0.0872,τ=0.4Nonlinear wave speed in the supercritical stable region.
(a)
ξ=0,τ=0(b)
ξ=0.0872,τ=0(c)
ξ=-0.0872,τ=0(d)
ξ=0,τ=0.2(e)
ξ=0.0872,τ=0.2(f)
ξ=-0.0872,τ=0.2(g)
ξ=0,τ=0.4(h)
ξ=0.0872,τ=0.4(i)
ξ=-0.0872,τ=0.4The above weakly nonlinear analysis shows the existence of both a supercritical stable region and a subcritical unstable region for isothermal and nonisothermal power-law liquid film down an inclined plane in the presence of shear stress at the interface. It is observed that the evolution of two-dimensional waves depends strongly on the initial disturbance wave number. There exists a value ofκs (obtained from equating denominator of E1=0) such that, when κs<κ<κc (obtained from di=0; cut-off wave number), the flow is supercritically stable and nonlinear equilibration occurs after the initial instability. The weakly nonlinear stability analysis shows that κs=κc/2, and the curve κ=κs separates the linearly unstable region into two portions where the nonlinear waves attain a finite equilibrium amplitude (di>0,E1>0) or reach an explosive state (di>0,E1<0), and the flow is supercritically equilibrated.


## 3.3. Nonlinear Analysis

In order to understand the mechanism responsible for the transfer of energy from the basic state to the disturbance, the evolution of finite-amplitude perturbations is considered. On the evolution of these perturbations, the influence of wind stress and mass transfer at the interface of a thin power-law liquid film down an inclined plane is examined by numerically solving the nonlinear evolution equation. The initial disturbance is taken to be a sinusoidal wave with small amplitude given by(3.17)H(x,0)=1-0.1cos⁡⁡(κx).
The evolution of waves with time is obtained by solving the evolution equation (2.18) by using Lee’s three-time level [44] finite-difference method in a periodic domain [-π/κ,π/κ]. The finite-difference scheme achieves linearity in the unknown Hi,j+1 by evaluating all coefficients of Hi,j+1 at a time level of known solution values, preserves stability by averaging Hi,j over three time levels, and maintains accuracy by using central difference approximations, where j denotes the number of time steps and i is a spatial grid point. Special attention is given to the discretization of the wave propagation term so as to overcome the solution blowup for very small initial wave numbers. The scheme possesses both conservative and transportive properties and allows to follow the solution over larger periods of time. The resulting algebraic equations are solved using MATLAB, and the iterations are continued until the local minimum layer thickness becomes smaller than the maximum error bound of 10-10. The computations are performed with Δt=0.01 and with the number of nodes along the spatial direction as N=100 and Δx=2π/Nκ. It is to be noted that the solution obtained with a decrease in Δt or an increase in N does not show any deviation from the values obtained for Δt=0.01, Δx=2π/Nκ, N=100. Further, it is confirmed that the wave that emerges on the film surface does not change its shape for a relatively long period of time.The weakly nonlinear analysis predicts the occurrence of a wave numberκs that separates the regimes of supercritical κs<κ<κc and subcritical κ<κs domains (Figure 4). For a wave number κ close to κc, the solution may evolve into a stable, almost sinusoidal wave of small finite amplitude. For κ close to κm=κc/2, the surface wave approaches the form of a solitary wave as observed in experimental results by Alekseenko et al. [45] and Liu et al. [46]. The evolution of the film in the supercritical stable region is analyzed by performing numerical simulations for initial disturbance wave number κ=κc,κm(=κc/2), and κs(=κc/2). The evolution of the film is examined for values of ξ in the range −0.0872 to 0.0872 and τ between 0 and 2.0.Figures7–9 show the evolution of the free surface configuration of isothermal (ξ=0) pseudoplastic (n=0.95), Newtonian (n=1.0), and dilatant (n=1.05) fluids at various instants of time for different values of shear stress parameter at the interface. Wave number is chosen in the range κs(=κc/2)<κ<κc (Figure 7—κ=0.06, κc=0.0868 when τ=0,κc=0.0971 when τ=0.2, and κc=0.1064 when τ=0.4; Figure 8—κ=0.05,κc=0.0569 when τ=0,κc=0.0635 when τ=0.2, and κc=0.0694 when τ=0.4; Figure 9—κ=0.03, κc=0.0365 when τ=0,κc=0.0409 when τ=0.2, and κc=0.0449 when τ=0.4). It is observed that the growth rate is much more important and the distortion of the free surface is conspicuous with increase in time. The wave amplitude decays after reaching a maximum, and the rear becomes longer and longer to attain a one-hump solitary-like wave at different instant of time for various values of τ. The large amplitude wave is then followed by a small amplitude capillary wave. In all the cases, it is observed that finite-amplitude permanent wave emerges. For a given n, amplitude of the wave increases with the increase in shear stress parameter τ. The corresponding evolution of maximum (Hmax) and minimum (Hmin) thickness is presented in Figure 10. The disturbance amplitude Hmax remains constant for initial times t~0–5 for pseudoplastic fluids (n=0.95), t~0–10 for Newtonian fluids (n=1.0), and t~0–20 for dilatant fluids (n=1.05) and then increases monotonically to attain a maximum amplitude. Isothermal pseudoplastic fluid attains the maximum amplitude in short time compared to Newtonian and dilatant fluid, and this amplitude further increases with the presence of shear stress at the interface. The distortion of the free surface is more significant when the effect of shear stress increases (Figure 11).Evolution of free surface for an isothermal (ξ=0) pseudoplastic (n=0.95) film when κ=0.06 and τ=0,0.2,0.4.
(a)(b)(c)Evolution of free surface for an isothermal (ξ=0) Newtonian (n=1.0) film when κ=0.05 and τ=0,0.2,0.4.
(a)(b)(c)Evolution of free surface for isothermal (ξ=0) dilatant (n=1.05) film when κ=0.03 and τ=0,0.2,0.4.
(a)(b)(c)Maximum and minimum amplitude of the waves down an inclined plane for an isothermal power-law film.
(a)
n=0.95,ξ=0.0(b)
n=1.0,ξ=0.0(c)
n=1.05,ξ=0.0Evolution of free surface for an isothermal film (ξ=0) when τ=2.0.
(a)(b)(c)Weakly nonlinear stability analysis predicted that for the condensate film (ξ>0), region of supercritical stability (di>0,E1>0) exists only for the pseudoplastic fluids (n=0.95). This is the case when the saturated vapor temperature is more than the wall temperature. When condensation is considered, the static layer grows thicker with time. Evolution of condensate pseudoplastic liquid film down an inclined plane is depicted in Figure 12 for ξ=0.01,0.02,0.04,and0.06 and τ=0. The result shows that pseudoplastic film down an inclined plane thickens faster as the phase change parameter ξ increases gradually in the absence of shear stress parameter τ. The effect of wind stress on the interface of the pseudoplastic film takes little longer time to thicken (Figure 13) and makes the film more corrugated as τ increases (Figures 13(a)(ii), 13(b)(ii), 13(c)(ii)).Evolution of free surface for condensating (ξ>0) pseudoplastic (n=0.95) film when κ=0.06 and τ=0.
(a)
n=0.95,ξ=0.01(b)
n=0.95,ξ=0.02(c)
n=0.95,ξ=0.04(d)
n=0.95,ξ=0.06Evolution of free surface for condensating (ξ=0.0872) pseudoplastic (n=0.95) film when κ=0.06 and τ=0,0.2,0.4.
(a)(b)(c)Figure14 shows the evolution of the free-surface shape for evaporating pseudoplastic, Newtonian, and dilatant fluids when ξ=-0.01 and −0.0872 with a zero shear stress at the interface. This is the case when the saturated vapor temperature is less than the wall temperature. It is observed that as n increases, depth of the liquid layer disappears in a finite time t=tr decreases leading to rupture instability. The rupture instability is thus augmented by the evaporative effect. Plot describes that the dilatant film evaporates faster than the pseudoplastic and Newtonian film. For a given n (given pseudoplastic, Newtonian, or dilatant fluid), evaporation of the film happens at earlier instant of time when the phase change parameter ξ decreases (Figure 14). Presence of wind stress on the interface does not influence the time of rupture tr for any n whereas it has little influence on the position where the rupture takes place (Figure 15).Evolution of free surface for an evaporating film (ξ<0) when κ=0.06 and τ=0.
(a)
n=0.95,ξ=-0.01(b)
n=1.0,ξ=-0.01(c)
n=1.05,ξ=-0.01(d)
n=0.95,ξ=-0.0872(e)
n=1.0,ξ=-0.0872(f)
n=1.05,ξ=-0.0872Evolution of free surface for an evaporating film (ξ=-0.00872) when κ=0.06 and τ=0.2,2.0.
(a)
n=0.95,τ=0.2(b)
n=1.0,τ=0.2(c)
n=1.05,τ=0.2(d)
n=0.95,τ=2.0(e)
n=1.0,τ=2.0(f)
n=1.05,τ=2.0


## 4. Conclusion

The influence of prescribed cocurrent superficial shear stress on the dynamics and stability of a condensate or evaporating power-law liquid film falling down an inclined plane has been analyzed by the method of long-wave perturbation. The interfacial boundary conditions include the effects of phase change across the interface. The evolution equation of the Benney type incorporating the effect of superficial shear stress on the surface of a nonisothermal power-law film has been derived. As the power-law exponent “n” decreases, the effective viscosity decreases, and, hence, “n” influences the Reynolds number Ren, Weber number Wen, and Prandtl number Prn. The results of the linear stability analysis reveal that the effect of decreasing the phase change parameter ξ or increasing the shear stress parameter τ is to destabilize the film flow system. Further, the dimensional quantities used to discuss the stability characteristics of the power-law model in terms of the power-law exponent “n” show that there exists a critical value of “n” for the condensate film, above which the film flow system is always stable (Figure 2), and this critical value of “n” increases with increase in shear stress parameter τ (Figure 2(d)), while isothermal and evaporating power-law fluid films are unstable for any value of power-law index “n”. The weakly nonlinear stability analysis of the power-law film flow system using long-wave theory is even more qualitative than the linearized stability results, and it reveals that both subcritical instability (di<0,E1<0) and supercritical stability (di>0,E1>0) are possible for isothermal and evaporating power-law fluid films and supercritical stability (di>0,E1>0) and supercritical explosive state (di>0,E1<0) are not possible for condensate dilatant fluid films (Figure 4).The nonlinear evolution of film thickness is found by numerically solving the evolution equation in a periodic domain. The numerical simulations confirm the results found on the basis of linear and weakly nonlinear stability analysis.  For an isothermal power-law liquid film, the finite amplitude waves are stable solutions of the evolution equation when the wave numbers are less than the cut-off wave numberκc.The permanent waves are nearly sinusoidal for initial wave numbers close to the cut-off wave number, and, for values much smaller than the cut-off value, the permanent waves are of solitary type. It is interesting to note that the range of supercritical stability broadens with an increase in the shear stress parameter τ for fixed values of other nondimensional parameters governing the film flow system. The nonlinear interactions in falling isothermal power-law films down an inclined plane exhibit a tendency toward permanent two-dimensional waves for very thin films at small wave numbers close to the cut-off wave number. The nonzero shear stress at the interface promotes the growth rate of the wave amplitude (Figure 10).The results reveal that (i) for an isothermal dilatant (pseudoplastic) liquids, the maximum wave amplitude is always smaller (larger) than that for a Newtonian liquid and the amplitude of permanent wave increases with the increase in prescribed shear stressτ at the interface in cocurrent direction; (ii) condensation of pseudoplastic film happens for the earlier instant of time when the phase change parameter increases and the effect of shear stress at the interface makes the film more corrugated; (iii) dilatant (pseudoplastic) evaporating liquid film attains rupture faster (slower) than that of Newtonian liquid film, and the prescribed shear stress τ at the interface does not influence the time of rupture and has little influence on the position where the rupture takes place.



---

*Source: 732675-2012-01-26.xml*