RETRACTED ARTICLE: Effect of electromagnetic field on the stability of viscoelastic fluid film flowing down an inclined plane

The stability of thin electrically conducting viscoelastic fluid film flowing down on a non-conducting inclined plane in the presence of electromagnetic field is investigated under induction-free approximation. Surface evolution equation is derived by long-wave expansion method. The stabilizing role of Hartman number M (magnetic field) and the destabilizing role of the viscoelastic property \({\varGamma}\) and the electric parameter E on such fluid film are established through the linear stability analysis of the surface evolution equation. Investigation shows that at small values of Hartman number (M), the influence of electric parameter (E) on the viscoelastic parameter \({(\varGamma)}\) is insignificant, while for large values of M, E introduces more destabilizing effect at low values of \({\varGamma}\) than that at high values of \({\varGamma }\). An interesting result also perceived from our analysis is that the stabilizing effect of Hartman number (M) is decreasing with the increase of the values of \({\varGamma}\) and E, even it gives destabilizing effect after a certain high value of the electric field depending on the high value of \({\varGamma}\). The weakly nonlinear study reveals that the increase of \({\varGamma}\) decreases the explosive and subcritical unstable zones but increases the supercritical stable zone keeping the unconditional zone almost constant.     

Nonlinear evolution of the travelling waves at the surface of a thin viscoelastic falling film

The problem of hydrodynamic instability of a thin condensate viscoelastic liquid film flowing down on the outer surface of an axially moving vertical cylinder is investigated. In order to improve the accuracy of numerical results, the viscoelastic and heat transfer parameters have been included into the governing equations. Also, the analytical solutions are obtained by utilizing the long-wave perturbation method. The influence of some physical parameters is discussed in both linear and nonlinear steps of the problem. It has been revealed that the stability of the film flow is weakened when the radius of cylinder and the temperature difference are reduced. Moreover, it is found that the increment of down-moving motion of the cylinder can enhance the flow stability. Further, the thin film flow can be destabilized by the viscoelastic property. The results show that both supercritical stability and subcritical instability can take place within the film flow system given appropriate conditions. Moreover, the absence of Reynolds number leads to an obvious difference in the behavior of some physical parameters.     

Hydromagnetic linear instability analysis of Giesekus fluids in plane Poiseuille flow

The effects of a fluid’s elasticity are investigated on the instability of plane Poiseuille flow on the presence of a transverse magnetic field. To determine the critical Reynolds number as a function of the Weissenberg number, a two-dimensional linear temporal stability analysis will be used assuming that the viscoelastic fluid obeys Giesekus model as its constitutive equation. Neglecting terms nonlinear in the perturbation quantities, an eigenvalue problem is obtained which is solved numerically by using the Chebyshev collocation method. Based on the results obtained in this work, fluid’s elasticity is predicted to have a stabilizing or destabilizing effect depending on the Weissenberg number being smaller or larger than one. Similarly, solvent viscosity and also the mobility factor are both found to have a stabilizing or destabilizing effect depending on their magnitude being smaller or larger than a critical value. In contrast, the effect of the magnetic field is predicted to be always stabilizing.     

Finite-amplitude long-wave instability of Bingham liquid films

The long-wave perturbation method is employed to investigate the weakly nonlinear hydrodynamic stability of a thin Bingham liquid film flowing down a vertical wall. The normal mode approach is first used to compute the linear stability solution for the film flow. The method of multiple scales is then used to obtain the weak nonlinear dynamics of the film flow for stability analysis. It is shown that the necessary condition for the existence of such a solution is governed by the Ginzburg–Landau equation. The modeling results indicate that both the subcritical instability and supercritical stability conditions can possibly occur in a Bingham liquid film flow system. For the film flow in stable states, the larger the value of the yield stress, the higher the stability of the liquid film. However, the flow becomes somewhat unstable in unstable states as the value of the yield stress increases.     

Stability characteristics of periodic streaming fluids in porous media

We study the linear stability of a three-layer flow of immiscible liquids located in a periodic normal electric field. We consider certain porous media assumed to be uniform, homogeneous, and isotropic. We analytically and numerically simulate the system of linear evolution equations of such a medium. The linearized problem leads to a system of two Mathieu equations with complex coefficients of the damping terms. We study the effects of the streaming velocity, permeability of the porous medium, and the electrical properties of the flow of a thin layer (film) of liquid on the flow instability. We consider several special cases of such systems. As a special case, we consider a uniform electric field and solve the transition curve equations up to the second order in a small dimensionless parameter. We show that the dielectric constant ratio and also the electric field play a destabilizing role in the stability criteria, while the porosity has a dual effect on the wave motion. In the case of an alternating electric field and a periodic velocity, we use the method of multiple time scales to calculate approximate solutions and analyze the stability criteria in the nonresonance and resonance cases; we also obtain transition curves in these cases. We show that an increase in the velocity and the electric field promote oscillations and hence have a destabilizing effect.     

Hydromagnetic instability of a power-law liquid film flowing down a vertical cylinder using numerical approximation approach techniques

The long-wave perturbation method is employed to investigate the hydromagnetic stability of a thin electrically-conductive power-law liquid film flowing down the external surface of a vertical cylinder in a magnetic field. The validity of the numerical results is improved through the introduction of the flow index and the magnetic force into the governing equation. In contrast to most previous studies presented in the literature, the solution scheme employed in this study is based on a numerical approximation approach rather than an analytical method. The normal mode approach is used to analyze the stability of the film flow. The modeling results reveal that the stability of the film flow system is weakened as the radius of the cylinder is reduced. However, the flow stability can be enhanced by increasing the intensity of the magnetic field and the flow index, respectively. In general, the optimum conditions can be found through the use of a system to alter stability of the film flow by controlling the applied magnetic field.     

Nonlinear stability characterization of thin Newtonian film flows traveling down on a vertical moving plate

The paper presents both the linear and nonlinear stability theories for the characterization of thin Newtonian film flows traveling down along a vertical moving plate. The linear model is first developed to characterize the flow behavior. After showing the inadequacy of the linear model in representing certain flow characteristics, the nonlinear kinematics model is then developed to represent the system. The long-wave perturbation method is employed to derive the generalized kinematic equations with free film surface condition. The linear model is solved by using the normal mode method for three different, namely, the quiescent, up-moving and down-moving, moving conditions. Subsequently, the elaborated nonlinear film flow model is solved by the method of multiple scales. The modeling results clearly indicate that both subcritical instability and supercritical stability conditions are possible to occur in the film flow system. The effect of the down-moving motion of the vertical plate tends to enhance the stability of the film flow.     

Defect correction method for viscoelastic fluid flows at high Weissenberg number

We study a defect correction method for the approximation of viscoelastic fluid flow. In the defect step, the constitutive equation is computed with an artificially reduced Weissenberg parameter for stability, and the resulting residual is corrected in the correction step. We prove the convergence of the defect correction method and derive an error estimate for the Oseen‐viscoelastic model problem. The derived theoretical results are supported by numerical tests for both the Oseen‐viscoelastic problem and the Johnson‐Segalman model problem. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Elastoviscoplastic fluid flow in non-circular tubes: Transversal field and interplay of elasticity and plasticity

An analytical study of elastoviscoplastic fluid flow in tubes of non-circular cross section is presented. The constitutive structure of the fluid is described by a linear frame invariant combination of the Phan-Thien−Tanner model of viscoelastic fluids and the Bingham model of plastic fluids. Non-circular tube cross sections are modeled by the shape factor method a one-to-one mapping of the circular base contour into a wide spectrum family of arbitrary tube contours. Field variables are expanded into asymptotic series in terms of the elasticity measure, the Weissenberg number We, coupled with an asymptotic expansion in terms of the geometrical mapping parameter ε leading to a set of hierarchical momentum balance equations which are solved successively up to and including the third order in We when the secondary field appears for the first time. The computational algorithm developed is applied to the study of the non-rectilinear flow in tubes with triangular and square cross sections. We find that the presence of the yield stress dampens the intensity of the purely viscoelastic vortices, the higher the yield stress the lower the intensity of the vortices in the cross-section, and the further away the vortices are from the center of the cross section as compared to the purely viscoelastic vortices. The results also evidence that viscoelasticity increases the axial flow for given viscoplastic conditions and pressure drop, and consequently increases the rate of flow, a phenomenon that may find applications in optimizing material transportation.     

Stability of distributed-parameter systems with retarded argument

Lyapunov functions are used to investigate the stability of processes described by a system of linear partial differential equations with retarded argument (for example: magnetohydrodynamic processes, elastic vibrations in aircraft, etc.). Some equations of the system may not involve time derivatives (for example, the equation of continuity in incompressible fluid flow, and the equation for the magnetic induction vector in the theory of electromagnetic phenomena). Such equations also arise when the order of a partial differential equation is reduced by introducing new notation for the space derivatives. A method is developed for investigating the stability of processes described by a system of this kind, some of whose equations do not contain time derivatives. Two constructions of the Lyapunov functions, as different integral quadratic forms, are proposed. Sufficient conditions for stability, in the form of inequalities relating the coefficients of the system, are established. As an example, the stability of the vibrations of a stretched string in a viscoelastic medium due to a distributed control force is considered.     

Explicit stability conditions for integro-differential equations with periodic coefficients

New explicit stability conditions are derived for a linear integro-differential equation with periodic operator coefficients. The equation under consideration describes oscillations of thin-walled viscoelastic structural members driven by periodic loads. To develop stability conditions two approaches are combined. The first is based on the direct Lyapunov method of constructing stability functionals. It allows stability conditions to be derived for unbounded operator coefficients, but fails to correctly predict the critical loads for high-frequency excitations. The other approach is based on transforming the equation under consideration in such a way that an appropriate ‘differential’ part of the new equation would possess some reserve of stability. Stability conditions for the transformed equation are obtained by using a technique of integral estimates. This method provides acceptable estimates of the critical forces for periodic loads, but can be applied to equations with bounded coefficients only. Combining these two approaches, we derive explicit stability conditions which are close to the Floquet criterion when the integral term vanishes. These conditions are applied to the stability problem for a viscoelastic bar compressed by periodic forces. The effect of material and structural parameters on the critical load is studied numerically. © 1998 B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Mathematical Modeling and Simulation of Antibubble Dynamics

In this study, we propose a mathematical model and perform numerical simulations for the antibubble dynamics. An antibubble is a droplet of liquid surrounded by a thin film of a lighter liquid, which is also in a heavier surrounding fluid. The model is based on a phase-field method using a conservative Allen-Cahn equation with a space-time dependent Lagrange multiplier and a modified Navier-Stokes equation. In this model, the inner fluid, middle fluid and outer fluid locate in specific diffusive layer regions according to specific phase filed (order parameter) values. If we represent the antibubble with conventional binary or ternary phase-field models, then it is difficult to have stable thin film. However, the proposed approach can prevent nonphysical breakup of fluid film during the simulation. Various numerical tests are performed to verify the efficiency of the proposed model.     

Couette flow of viscoelastic fluid with constitutive relation involving general Caputo-type fractional derivative

In this paper, the general Caputo-type fractional differential operator introduced by Pr. Anatoly N. Kochubei is applied to the linear theory of viscoelasticity. Firstly, using the general Caputo-type derivative, a generalized linear viscoelastic constitutive equation is proposed for the first time. Secondly, the momentum equation for the plane Couette flow of viscoelastic fluid with the constitutive relation is given as an integrodifferential equation and the analytical solution of the equation is established by employing the separation of variables method. Lastly, for special cases of the general constitutive relation, the analytical solutions are obtained in terms of the Mittag-Leffler functions.     

Rimming flow of a power-law fluid: Qualitative analysis of the mathematical model and analytical solutions

Rimming flow of a non-Newtonian fluid on the inner surface of a horizontal rotating cylinder is investigated. Simple lubrication theory is applied since the Reynolds number is small and liquid film is thin. For the steady-state flow of a power-law fluid the mathematical model reduces to a simple algebraic equation regarding the thickness of the liquid film. The qualitative analysis of this equation is carried out and the existence of two possible solutions is rigorously proved. Based on this qualitative analysis, different regimes of the rimming flow are defined and analyzed analytically. For the particular case, when the flow index in a power-law constitutive equation is equal to 1/2, the problem reduces to the fourth order algebraic equation which is solved analytically by Ferrari method.     

Stability of a thin liquid layer,spreading on a quasi–rotating disk

In many industrial processes solids are coated to obtain specific surface properties, as e.g. corrosion resistance, mechanical (wear) resistance, optical, or electrical properties. Even today many coating processes are not fully understood and the choice of parameters is largely based on experience. Hence, a prediction of the complete hydrodynamic process and the appearance of instabilities in its dependency on the parameters appears highly desirable. This would serve to optimize the quality of the coating. A common coating technique is the so-called spin coating. The coating agent is dissolved or suspended in a liquid, brought onto the solid, spread by rotation, and the carrier liquid is finally removed by evaporation or by chemical reactions. In this article an evolution equation is derived from lubrication theory, valid for thin liquid layers. The model involves a dynamic contact angle, centrifugal, capillary, and gravitational forces. The evolution equation can be solved analytically, provided the capillary number is small. Then a coupled linear stability analysis of the contact line and the free interface is performed. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Multiscale analysis and numerical simulation for stability of the incompressible flow of a Maxwell fluid

How to predict the stability of a small-scale flow subject to perturbations is a significant multiscale problem. It is difficult to directly study the stability by the theoretical analysis for the incompressible flow of a Maxwell fluid because of its analytical complexity. Here, we develop the multiscale analysis method based on the mathematical homogenization theory in the stress–stream function formulation. This method is used to derive the homogenized equation which governs the transport of the large-scale perturbations. The linear stabilities of the large-scale perturbations are analyzed theoretically based on the linearized homogenized equation, while the effect of the nonlinear terms on the linear stability results is discussed numerically based on the nonlinear homogenized equation. The agreements between the multiscale predictions and the direct numerical simulations demonstrate the multiscale analysis method is effective and credible to predict stabilities of flows.     

Finite element method for viscoelastic fluids flows: approximation of the Maxwell model

We discuss about the finite element approximation of viscoelastic fluid flow. We consider a fluid obeying the Oldroyd model and particularly we study the purely viscoelastic case, the so-called Maxwell model, important in practice for the applications. We examine two kinds of methods used for the approximation of the Maxwell model: method using a splitting technique and finite element method satisfying inf-sup conditions relating tensor and velocity. We present numerical results for these methods and we discuss about their stability. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-linear stability characterization of the thin micropolar liquid film flowing down the inner surface of a rotating vertical cylinder

Linear and non-linear stability analysis for characterization of micropolar film flowing down the inner surface of a rotating infinite vertical cylinder is given. A generalized non-linear kinematic model is derived to represent the physical system and is solved by the long wave perturbation method in the following procedure. First, the normal mode method is used to characterize the linear behaviors. Then, an elaborated non-linear film flow model is solved by using the method of multiple scales to characterize flow behaviors at various states of sub-critical stability, sub-critical instability, supercritical stability, and supercritical explosion. The modeling results indicate that by increasing the rotation speed, Ω, and the radius of cylinder, R, the film flow will generally stabilize the flow system.     

Breakup of finite fluid films

We study the dewetting process of thin fluid films that partially wet a solid surface. Using long wave (lubrication) approximation, we formulate a nonlinear partial differential equation governing the evolution of the film thickness, h. This equation includes the effects of capillarity, gravity, and additional conjoining/disjoining pressure term to account for intermolecular forces. We perform standard linear stability analysis of an infinite flat film, and identify the corresponding stable, unstable and metastable regions. Within this framework, we analyze the evolution of a semi-infinite film of length L in one direction. The numerical simulations show that for long and thin films, the dewetting fronts of the film generate a pearling process involving successive formation of ridges at the film ends and consecutive pinch-off behind these ridges. On the other hand, for shorter and thicker films, the evolution ends up by forming a single drop. The time evolution as well as the final drops pattern shows a competition between the dewetting mechanisms caused by nucleation and by free surface instability. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)