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.     

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 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.     

Solution of the laminar viscous flow in a semi-porous channel in the presence of a uniform magnetic field by using the homotopy analysis method

In this paper, the problem of laminar viscous flow in a semi-porous channel in the presence of a transverse magnetic field is presented and the homotopy analysis method (HAM) is employed to compute an approximation to the solution of the system of nonlinear differential equations governing the problem. It has been attempted to show the capabilities and wide-range applications of the homotopy analysis method in comparison with the numerical method in solving this problem. The obtained solutions, in comparison with the numeric solutions admit a remarkable accuracy. A clear conclusion can be drawn from the numerical method’s (NM) results that the HAM provides highly accurate solutions for nonlinear differential equations.     

Analytical solution of magnetohydrodynamic stagnation-point flow of a power-law fluid towards a stretching surface

A new kind of analytic technique, namely the homotopy analysis method (HAM), is employed to give an explicit analytical solution of the steady two-dimensional stagnation-point flow of an electrically conducting power-law fluid over a stretching surface when the surface is stretched in its own plane with a velocity proportional to the distance from the stagnation-point. A uniform transverse magnetic field is applied normal to the surface. An explicit analytical solution is given by recursive formulae for the first-order power-law (Newtonian) fluid when the ratio of free stream velocity and stretching velocity is not equal to unity. For second and real order power-law fluids, an analytical approach is proposed for magnetic field parameter in a quite large range. All of our analytical results agree well with numerical results. The results obtained by HAM suggest that the solution of the problem under consideration converges.     

Thermo-magneto-dynamic stresses and perturbation of magnetic field vector in a non-homogeneous hollow cylinder

An analytical method is presented to investigate thermo-magneto-elastic stresses and perturbation of the magnetic field vector in a conducting non-homogeneous hollow cylinder under thermal shock. The interaction between the deformation and the magnetic field vector in a non-homogeneous hollow cylinder is considered by adding a Lorentz’s electro-magneto-force into the equation of thermo-elastic motion of the non-homogeneous hollow cylinder in an axial magnetic field. The exact solution for magneto-thermo-dynamic stresses and perturbation responses of an axial magnetic field vector in a conducting non-homogeneous hollow cylinder was obtained by using finite integral transforms. From numerical calculations, the dynamic characteristics on both thermo-magneto-stresses and perturbation of the axial magnetic field vector in the conducting non-homogeneous hollow cylinder is revealed and discussed.     

Nonlinear Kelvin-Helmholtz instability of two miscible ferrofluids in porous media

The nonlinear theory of the Kelvin-Helmholtz instability is employed to analyze the instability phenomenon of two ferrofluids through porous media. The effect of both magnetic field and mass and heat transfer is taken into account. The method of multiple scale expansion is employed in order to obtain a dispersion relation for the first-order problem and a Ginzburg–Landau equation, for the higher-order problem, describing the behavior of the system in a nonlinear approach. The stability criterion is expressed in terms of various competing parameters representing the mass and heat transfer, gravity, surface tension, fluid density, magnetic permeability, streaming, fluid thickness and Darcy coefficient. The stability of the system is discussed in both theoretically and computationally, and stability diagrams are drawn.     

Electromagnetic interaction of a conducting cylinder with a magnetic dipole caused by steady rotation

The motion of a conductor in a magnetic field induces eddy currents whose interaction with the field produces Lorentz forces opposing the motion. One can determine the velocity of the conductor from the force on the magnet system since the latter is equal but opposite to the Lorentz force on the conductor. This contactless method is known as Lorentz force velocimetry (LFV). We study an idealized configuration of LFV, i.e. a rotating solid cylinder interacting with a point dipole. The understanding of parameter influences in this setup can be helpful for more realistic configurations. We use a purely kinematic approach appropriate for low magnetic Reynolds numbers. Numerical results for small and large distances between dipole and cylinder have been obtained with the commercial software COMSOL Multiphysics. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetic Suppression of Thermocapillary Convection in a Cylindrical Specimen

Convective flows may adversely affect the homogeneity of semiconductor crystals grown by the floating-zone method. The intensity of the convective flow can be reduced by an external magnetic field. The study simulates the convective flow in a cylindrical specimen. The mathematical model consists of the system of quasihydrodynamic (QHD) equations for a quasineutral electrically conducting fluid in an external electromagnetic field. The numerical results show how the convective flow velocities and configurations depend on the magnetic field strength.     

MHD flow of a Newtonian fluid over a stretching sheet: an approximate solution

An approximate solution to the problem of steady laminar flow of a viscous incompressible electrically conducting fluid over a stretching sheet is presented. The approach is based on the idea of stretching the variables of the flow problem and then using least squares method to minimize the residual of a differential equation. The effects of the magnetic field on the flow characteristics are demonstrated through numerical computations with different values of the Hartman number.     

Nonlinear Kelvin-Helmholtz instability of two miscible ferrofluids in porous media

The nonlinear theory of the Kelvin-Helmholtz instability is employed to analyze the instability phenomenon of two ferrofluids through porous media. The effect of both magnetic field and mass and heat transfer is taken into account. The method of multiple scale expansion is employed in order to obtain a dispersion relation for the first-order problem and a Ginzburg–Landau equation, for the higher-order problem, describing the behavior of the system in a nonlinear approach. The stability criterion is expressed in terms of various competing parameters representing the mass and heat transfer, gravity, surface tension, fluid density, magnetic permeability, streaming, fluid thickness and Darcy coefficient. The stability of the system is discussed in both theoretically and computationally, and stability diagrams are drawn. Received: July 25, 2002; revised: April 16, 2003     

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.     

Mathematical modeling of coupled turbulent flow and solidification in a single belt caster with electromagnetic brake

A mathematical model has been developed to simulate turbulent fluid flow and solidification in the presence of a DC magnetic field in an extended nozzle for metal delivery to a single belt caster. This paper reports on predicted effects of DC magnetic field conditions in modifying flows and solidification behavior in the metal delivery system. It is shown that the application of a DC magnetic brake to the proposed system can result in a reasonably uniform feeding of melt onto the cooled moving belt. This, in turn, optimises the rate of even shell growth along the chilled substrate. In order to account for the effects of turbulence, a revised low-Reynolds k–ε turbulent model was employed. A Darcy-porosity approach was used to simulate fluid flow within the mushy solidification region. Simulations were carried out for plain carbon steel strip casting. The fully coupled transport equations were numerically solved using the finite volume method. The computed flow patterns were compared with those reported in the literature. The performance of the magnetic flow control device proposed in this work is evaluated and compared with flow modifications obtained by inserting a ceramic filter within the reservoir.     

Hydromagnetic fluctuating flow of a couple stress fluid through a porous medium

The equations of a polar fluid of hydromagnetic fluctuating through a porous medium are cast into matrix form using the state space and Laplace transform techniques the resulting formlation is applied to a variety of problems. The solution to a problem of an electrically conducting polar fluid in the presence of a transverse magnetic field and to a probem for the flow between two parallel fixed plates is obtained. The inversion of the Laplace transforms, is carried out using a numerical approach. Numerical results for the velocity, angular velocity distribution and the induced magnetic field are given and illustrated graphically for each problems.     

Rotation implementation of a circular cylinder in incompressible flow via staggered grid approach

In this paper, we present a finite difference method for the implementation of the rotation of a circular cylinder in the incompressible flow field by solving the two-dimensional unsteady Navier-Stokes equations. The approach is to use staggered grid method so that the accuracy and order of convergence of the associated algorithms can be maintained. The proposed method is easy to be implemented and is effective. A set of simulations for the flow dynamics is provided to show the computational results.     

Numerical simulation of a flexible fiber deformation in a viscous flow by the immersed boundary-lattice Boltzmann method

In this paper, deformation of a mass-less elastic fiber with a fixed end, immersed in a two-dimensional viscous channel flow, is simulated numerically. The lattice-Boltzmann method (LBM) is used to solve the Newtonian flow field and the immersed-boundary method (IBM) is employed to simulate the deformation of the flexible fiber interacting with the flow. The results of this unsteady simulation including fiber deformation, fluid velocity field, and variations of the fiber length are depicted in different time-steps through the simulation time. Similar trends are observed in plots representing length change of fibers with different values of stretching constant. Also, the numerical solution reaches a steady state equivalent to the fluid channel flow over a flat plate.     

A stable boundary elements method for magnetohydrodynamic channel flows at high Hartmann numbers

The article is devoted to extension of boundary element method (BEM) for solving coupled equations in velocity and induced magnetic field for time dependent magnetohydrodynamic (MHD) flows through a rectangular pipe. The BEM is equipped with finite difference approach to solve MHD problem at high Hartmann numbers up to 10⁶. In fact, the finite difference approach is used to approximate partial derivatives of unknown functions at boundary points respect to outward normal vector. It yields a numerical method with no singular boundary integrals. Besides, a new approach is suggested in this article where transforms 2D singular BEM's integrals to 1D nonsingular ones. The new approach reduces computational cost, significantly. Note that the stability of the numerical scheme is proved mathematically when computational domain is discretized uniformly and Hartmann number is 40 times bigger than length of boundary elements. Numerical examples show behavior of velocity and induced magnetic field across the sections.     

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).     

Locomotion based on the control of the shape of magnetic fluid surfaces and of magnetizable media

The realization of locomotion based on the deformation of a free surface of a magnetic fluid layer in a traveling magnetic field is studied. A plane flow of an incompressible viscous magnetic fluid layer on a horizontal surface in a nonuniform magnetic field and a plane two-layers flow of incompressible viscous magnetic fluids between two parallel solid planes in a magnetic field is considered. Also the flow of an incompressible viscous magnetic fluid layer on a cylinder in a nonuniform magnetic field is an object of investigation. The deformation and the motion of a body made by a magnetizable polymer in an alternating magnetic field are experimentally studied. The cylindrical body (worm) which is located in a cylindrical tube is analyzed. These effects can be used in designing autonomous mobile robots without a hard cover. Such robots can be employed in clinical practice and biological investigations. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical study on viscoelastic fluid flow past a rigid body

Viscoelastic non-Newtonian fluids can be achieved by adding a small amount of polymer additives to a Newtonian fluid. In this paper, numerical simulations are used to investigate the influence of such polymer additives on the behavior of flow past a circular cylinder. A numerical method is proposed that discretizes the non-linear viscoelastic system on a uniform Cartesian grid, with a penalization method to model the presence of the cylinder. The drag of the cylinder and the flow behavior under the effect of different Reynolds numbers ($\mathrm{Re}$), Weissenberg numbers (Wi) and polymer viscosity ratios (ε) are studied. Numerical results show that different flow characteristics are exhibited in different parameter zones. The polymer viscosity ratio plays an important role at low Weissenberg and Reynolds numbers, but as the Reynolds and Weissenberg numbers increase, the influence of ε weakens. The drag force of the cylinder is mostly affected by the Reynolds and Weissenberg numbers. At low Reynolds numbers, the drag of the cylinder and the flow fields are only affected by a large value of Wi when the elastic forces are strong. Non-trivial drag reduction occurs only when there is vortex shedding in the wake flow, whereas drag enhancement happens when the vortex shedding is inhibited.