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.     

Finite amplitude long-wave instability of a thin viscoelastic fluid during spin coating

This paper investigates the stability of a thin incompressible viscoelastic fluid designated as Walters’ liquid B″ during spin coating. The long-wave perturbation method is proposed to derive a generalized kinematic model of the film flow. The method of normal mode is applied to study the linear stability. The amplitude growth rates and the threshold conditions are characterized subsequently and summarized as the by-products of the linear solutions. Using the multiple scales method, the weakly nonlinear stability analysis is studied for the evolution equation of a film flow. The Ginzburg–Landau equation is determined to discuss the threshold conditions of the various critical flow states. The study reveals that the rotation number and the radius of the rotating circular disk generate the destabilizing effects. Moreover, the viscoelastic parameter k indeed plays a more significant role in destabilizing the film flow than a thin Newtonian fluid during spin coating [27].     

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.     

Does Surface Tension Stabilize Liquid Films Inside a Rotating Horizontal Cylinder? Part 2: Multidimensional Disturbances

We examine the stability of a thin film of viscous fluid on the inside surface of a cylinder with horizontal axis, rotating about this axis. Depending on the parameters involved, the dynamics of the film can be described by several asymptotic models, one of which was examined by Benilov [ J. Fluid Mech. 501:105–124 (2004)]. It turned out that the linearized stability problem for this model admits infinitely many neutrally stable eigenmodes, which form a complete set. Despite that, the film is unstable with respect to exploding disturbances, which grow infinitely in a finite time. The present paper examines the effect of surface tension on the stability of the film. Two cases are considered: short-scale disturbances (such that the axial wavelength λ is much smaller than the radius R of the cylinder) and long-scale disturbances (for which  λ≳ R   ). In the former case, surface tension is a stabilizing influence, because it regularizes the exploding solutions and makes all eigenmodes asymptotically (not just neutrally) stable. The latter case was previously examined by Acrivos and Jin [ J. Eng. Math. 50:99–120 (2004)], who showed that surface tension destabilizes some of the eigenmodes. We argue, however, that the corresponding growth rate is much smaller than that of the so-called inertial instability.     

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.     

简支Mexwell模型粘弹性输流管道的稳定性分析

在弹性输流管道研究的基础上,采用递推格式的有限差分法,对简支Maxwell模型粘弹性输流管道（回转守恒系统）,探讨了其动力特性和稳定性问题,具体分析了材料的松弛时间对无量纲流速与前三阶模态的无量纲频率的实部及虚部之间的变化曲线的影响。发现发散临界流速随松弛时间的减小而降低,随后发生的耦合模态颤振临界流速随松弛时间的减小而增大;甚至在质量比较大时,随着松弛时间的减小,可推迟乃至不发生耦合模态颤振。当无量纲松弛时间达到10^3量级以上时,即可将其按弹性管道处理。甚至在H为10^2量级时,按弹性管道处理也不会带来太大的误差。     

On the Equation of a Rotating Film

We study positive periodic solutions to a nonautonomous nonlinear third-order ordinary differential equation of the theory of motion of a viscous incompressible fluid with free boundary. This equation describes the steady motions of a thin layer of a fluid film on the surface of a rotating horizontal cylinder in the gravity field. The linear operator on the left-hand side of the equation has a three-dimensional kernel. Moreover, the equation contains two nonnegative parameters proportional to the gravity acceleration and surface tension. Depending on these parameters the problem in question may have either two solutions or no solutions at all. We establish some qualitative properties of solutions to the problem: in particular, their asymptotic behavior at the extremal values of the parameters.     

Stability analysis of a single-walled carbon nanotube conveying pulsating and viscous fluid with nonlocal effect

For a single-walled carbon nanotube (CNT) conveying fluid, the internal flow is considered to be pulsating and viscous, and the resulting instability and parametric resonance of the CNT are investigated by the method of averaging. The partial differential equation of motion based on the nonlocal elasticity theory is discretized by the Galerkin method and the averaging equations for the first two modes are derived. The stability regions in frequency–amplitude plane are obtained and the influences of nonlocal effect, viscosity and some system parameters on the stability of CNT are discussed in detail. The results show that decrease of nonlocal parameter and increase of viscous parameter both increase the fundamental frequency and critical flow velocity; the dynamic stability of CNT can be enhanced due to nonlocal effect; the contributions of the fluid viscosity on the stability of CNT depend on flow velocity; the axial tensile force can only influence natural frequencies of the system however the viscoelastic characteristic of materials can enhance the dynamic stability of CNT. The conclusions drawn in this paper are thought to be helpful for the vibration analysis and structural design of nanofluidic devices.     

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.     

Analysis of viscoelastic fluid flow with temperature dependent properties in plane Couette flow and thin annuli

The steady state flow in very thin annuli has been studied analytically for the case where the annular gap is much smaller than the radius of the inner cylinder and for the outer cylinder rotating at constant angular speed and the inner cylinder at rest. The cylinders were subjected to two different thermal boundary conditions. The exponential effect of temperature on the relaxation time and the viscosity coefficient was accounted into the governing differential equations using Nahme’s law. Effects of viscous dissipation as well as εDe² (viscoelastic index for SPTT constitutive equation) on the dimensionless velocity and temperature profiles have been investigated. Results show that while the properties of the fluid depend on temperature, the velocity and temperature profiles are different compared to those obtained with constant physical properties. The Nahme–Griffith number increases whereas εDe² as a viscoelastic index decreases when temperature dependent physical properties are considered. In addition, the results indicate that the viscous dissipation has a sensible effect on heat transfer and the Nusselt number decreases with an increase in the Nahme–Griffith number.     

On the three-dimensional undulation instability of a rectangular viscoelastic composite plate in biaxial compression

Within the frame work of the second version of small precritical deformation in the three-dimensional linearized theory of stability of deformable bodies (TDLTSDB), the undulation instability problem for a simply supported rectangular plate made of a viscoelastic composite material is investigated in biaxial compression in the plate plane. The corresponding boundary-value problem is solved by employing the Laplace transformation and the principle of correspondence. For comparison and estimation of the accuracy of results given by the TDLTSDB, the same problem is also solved by using various approximate plate theories. The viscoelasticity properties of the plate material are described by the Rabotnov fractional-exponential operator. The numerical results and their discussion are presented for the case where the plate is made of a multilayered viscoelastic composite material. In particular, the variation range of problem parameters is established for which it is necessary to investigate the undulation instability of the viscoelastic composite plate by using the TDLTSDB. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 45, No. 1, pp. 93–108, January–February, 2009.     

Stability analysis of an interface between immiscible liquids in Hele-Shaw flow in the presence of a magnetic field

Hydrodynamic instabilities may occur when a viscous fluid is driven by a less viscous one through a porous medium. These penetrations are common in enhanced oil recovery, dendrite formation and aquifer flow. Recent studies have shown that the use of magnetic suspensions allow the external control of the instability. The problem is nonlinear and some further improvements of both theory and experimental observations are still needed and continue being a current source of investigation. In this paper we present a generalized Darcy law formulation in order to examine the growth of finger instabilities as a magnetic field is applied to the interface between the fluids in a Hele-Shaw cell. A new linear stability analysis is performed in the presence of magnetic effects and provides a stability criterion in terms of the non-dimensional physical parameters of the examined flow and the wavenumber of the finger disturbances. The interfacial tension inhibits small wavelength instabilities. The magnetic field contributes to the interface stability for moderate wavelength as it is applied parallel to the liquid-interface. In particular, we find an explicit expression, as a function of the susceptibility, for a critical angle between the interface and the magnetic field direction, in which its effect on the interface is neutral. We have developed a new asymptotic solution for the flow problem in a weak nonlinear regime. The first correction captures the second order nonlinear effects of the magnetic field, which tends to align the fingers with the field orientation and have a destabilizing effect. The analysis predicts that the non-linear effects at second order can counterbalance the first order stabilizing effect of a parallel magnetic field which results in a loss of effectiveness for controlling the investigated finger instabilities. The relevant physical parameters for controlling these finger instabilities are clearly identified by our non-dimensional analysis.     

Viscoelastic cylinder reinforced with thin elastic casing in a heterogeneous temperature field

The problem of a viscoelastic cylinder reinforced with a thin elastic casing with consideration given to a heterogeneous temperature field is presented. The method of numerical solution of the problem is discussed. The results of computation are presented in the form of curves.     

Numerical stability analysis of the Taylor-Couette flow in the two-dimensional case

The stability of the laminar flow between two rotating cylinders (Taylor-Couette flow) is numerically studied. The simulation is based on the equations of motion of an inviscid fluid (Euler equations). The influence exerted on the flow stability by physical parameters of the problem (such as the gap width between the cylinders, the initial perturbation, and the velocity difference between the cylinders) is analyzed. It is shown that the onset of turbulence is accompanied by the formation of large vortices. The results are analyzed and compared with those of similar studies.     

The free-parameter solution of the convection equations in a vertical cylinder with a volume heat source

An exact solution of the free-convection equations is constructed in the Oberbeck–Boussinesq approximation, describing the flow of a viscous heat-conducting fluid in a vertical cylinder of large radius when heated by radiation. The initial problem is reduced to an operator equation with an extremely non-linear operator, satisfying Schauder's theorem in C[0,1]. An iteration procedure is proposed for determining the independent parameter, that occurs in the solution, which enables three different values to be obtained and, correspondingly, three classes of solution of the initial problem. The linear stability of all the solutions obtained is investigated and it is shown that, for chosen values of the problem parameters, the most dangerous one is the plane wave mode and two instability mechanisms are present. The flow structure and the type of instability depend considerably on the values of the free parameter.     

UNSTEADY BOUNDARY LAYER FLOW OF FRACTIONAL OLDROYD-B VISCOELASTIC FLUID PAST A MOVING PLATE IN A POROUS MEDIUM

This paper considers the unsteady boundary layer flow over a moving flat plate embedded in a porous medium with fractional Oldroyd-B viscoelastic fluid. The governing equations with mixed time-space fractional derivatives are solved numerically by using the finite difference method combined with an L1-algorithm. The effect of various physical parameters on the velocity and average skin friction are discussed and graphically illustrated in detail.Results show that the porosity € and fractional derivative α enhance the flow of Oldroyd-B viscoelastic fluid within porous medium, but fractional derivative βweakens the flow. Moreover, it is found that the average skin friction coefficient rises with the increase of fractional derivative β.     

3D Analyses of the symmetric local stability loss of the circular hollow cylinder made from viscoelastic composite material

The 3D approach was employed for investigations of the symmetric local stability loss of the circular hollow cylinder made from the viscoelastic composite materials. This approach is based on investigations of the development of the initial rotationally symmetric infinitesimal local imperfections of the circular hollow cylinder within the scope of 3D geometrically nonlinear field equations of the theory of viscoelasticity for anisotropic bodies. The numerical results of the critical force and critical time are presented and discussed. For comparison and estimation of the accuracy of the results given by the 3D approach, the same problem is also solved by using various approximate shell theories. The viscoelasticity properties of the plate material are described by the fractional–exponential operator. The numerical results and their discussion are presented for the case where the cylinder is made of a uni-directional fibrous viscoelastic composite material. In particular, it is established that the difference between the critical times obtained by employing 3D and third order refined shell theories becomes more non-negligible if the values of the external compressive force are close to the critical compressive force which is obtained at t = ∞ (t denotes a time).     

Dynamics of elastic bodies connected by a thin soft viscoelastic layer

A dynamic study was performed on a structure consisting of two three-dimensional linearly elastic bodies connected by a thin soft nonlinear Kelvin–Voigt viscoelastic adhesive layer. The adhesive is assumed to be viscoelastic of Kelvin–Voigt generalized type, which makes it possible to deal with a relatively wide range of physical behavior by choosing suitable dissipation potentials. In the static and purely elastic case, convergence results when geometrical and mechanical parameters tend to zero have already been obtained using variational convergence methods. To obtain convergence results in the dynamic case, the main tool, as in the quasistatic case, is a nonlinear version of Trotter?s theory of approximation of semigroups acting on variable Hilbert spaces. The limit problem involves a mechanical constraint imposed along the surface to which the layer shrinks. The meaning of this limit with respect to the relative behavior of the parameters is discussed. The problem applies in particular to wave phenomena in bonded domains.     

Coexistence and stability of a competition—diffusion system in population dynamics

The coexistence and stability of the population densities of two competing species in a bounded habitat are investigated in the present paper, where the effect of dispersion (transportation) is taken into consideration. The mathematical problem involves a coupled system of Lotka-Volterra-type reaction-diffusion equations together with some initial and boundary conditions, including the Dirichlet, Neumann and third type. Necessary and sufficient conditions for the coexistence and competitive exclusion are established and the effect of diffusion is explicitly given. For the stability problem, general criteria for the stability and instability of a steady-state solution are established and then applied to various situations depending on the relative magnitude among the physical parameters. Also given are necessary and sufficient conditions for the existence of multiple steady-state solutions and the stability or instability of each of these solutions. Special attention is given to the Neumann boundary condition with respect to which some threshold results for the coexistence and stability or instability of the four uniform steady states are characterized. It is shown in this situation that only one of the four constant steady states is asymptotically stable while the remaining three are unstable. The stability or instability of these states depends solely on the relative magnitude among the various rate constants and is independent of the diffusion coefficients.     

Pattern formation in microwave heated ceramics: cylinders and slabs

Analyses of microwave heating of a thin ceramic cylinder and a thin ceramic slab in a single mode, highly resonant cavityare presented. Realistic assumptions regarding the effectiveelectrical conductivity, thermal parameters, and physical dimensionsare adhered to throughout. Consequently, the models developedherein incorporate most of the features of actual experiments.They incorporate both the effects of cavity detuning and alocal electric field perturbation on the heating process. The models presented take the form of one- and two-dimensionalreaction–diffusion equations which contain a functionaland an inhomogeneous source term for the cylinder and slab,respectively. The development of these equations is the product of a systematic modelling process that involves S-matrix theory,a small Biot number asymptotic analysis, and a matched asymptotic analysis of a non-standard electromagnetic scattering problem.The one-dimensional equation for the cylinder reveals boththe mathematical structure and physical mechanism for the formationof hot-spots. The two-dimensional equation supports a hot stripe pattern, due to preferential electromagnetic heating, whichbecomes unstable and evolves into an oval-like spot. Accuratenumerical methods which approximate the solutions of theseequations and their stability are presented and these agreequalitatively with experiments and predict observed trends. Received 1 June, 1999. ⁺ antoine@mip.ups-tlse.fr