Annular liquid jets at high Reynolds numbers

The flow of annular liquid jets at high Reynolds numbers is analysed by means of the finite element method and the full‐Newton iteration scheme. Results have been obtained for various values of the inner to the outer diameter ratio and for non‐zero surface tension, using extremely long meshes. The annular film moves far from the symmetry axis at low values of the Reynolds number. At higher Reynolds numbers, the film moves towards the axis of symmetry and appears close to very far downstream, forming a round jet. Asymptotic results for the radius of the resulting round jet are provided. Copyright © 2003 John Wiley & Sons, Ltd.     

雷列台阶-环槽端面密封机理与性能研究

基于质量守恒的JFO空化边界条件建立了雷列台阶-环槽机械密封端面润滑理论模型,采用有限单元法求解Reynolds控制方程,获得了端面膜压、密度比与液膜流线分布,分析了其密封机理与性能规律.结果表明:密封环端面内径侧的圆环浅槽和端面中部的圆环深槽组合结构可造成合理的空化现象,达到空化减漏的目的.其中,端面中部的圆环深槽是端面高压侧雷列台阶和端面低压侧圆环浅槽的隔离带,使得端面低压侧压力分布受端面高压侧压力分布的影响极小.密封流体进入低压侧圆环浅槽时,端面间隙突然发散,压力低于该工况温度下的饱和蒸汽压,整个圆环浅槽区液膜空化,达到零泄漏并出现大量的回流现象.外径侧雷列台阶提供良好的动压承载能力,实现了端面的非接触.雷列台阶和环槽织构的组合应用使得该机械密封具有优良的综合性能.     

Computer‐algebra multiple‐timescale method for spatially periodic thin‐film viscous‐flow problems

A flexible, fully automated, computer‐algebra algorithm is developed for solving a class of non‐linear partial‐differential evolution equations arising frequently in the modeling of two‐dimensional transient free‐surface viscous thin‐film flows. The method, which is formulated for solving spatially periodic problems, is based upon an explicit multiple‐timescale asymptotic approximation of the thin‐film thickness. It admits the resolution of diverse physical phenomena by employing a finite geometric progression of increasingly slow timescales. The method is implemented on a challenging test problem comprising the evolution of an annular film of viscous liquid, with a free surface, adhering to the exterior of a horizontal rotating circular cylinder; as a model for numerous industrially motivated coating flows, this benchmark problem has been analyzed in diverse numerical and theoretical studies, against whose results those of the present method are compared. The explicit algebraic form of the solution admits a study of large‐time evolutionary dynamics that lies beyond the reach of considerably more expensive conventional numerical solvers, thereby shedding new light on the hitherto‐undiscovered explicit dependence of large‐time evolutionary fluid dynamics in terms of independent parameters describing gravitational and capillary effects. The results obtained from the new computer‐algebra procedure are demonstrated to be in good agreement with those obtained from a bespoke efficient numerical integration method that is spectrally accurate in space and 8th‐order (Runge–Kutta) in time. Newly discovered mechanisms describing the decay of free‐surface wave modes, from arbitrary initial conditions to the steady state, are presented. Copyright © 2010 John Wiley & Sons, Ltd.     

Stability analysis of a polymer film casting problem

The polymer cast film process consists of stretching a molten polymer film between a flat die and a drawing roll. Drawing instabilities are often encountered and represent a drastic limitation to the process. Newtonian fluid film stretching stability is investigated using two numerical strategies. The first one is a ‘tracking’ method, which consists of solving Stokes equations in the whole fluid area (extrusion die and stretching path) by finite elements. The interface is determined to satisfy a kinematic equation. A domain decomposition meshing technique is used in order to account for a flow singularity resulting from the change in the boundary conditions between the die flow region and the stretching path region. A linear stability method is then applied to this transient kinematic equation in order to investigate the stability of the stationary solution. The second method is a direct finite element simulation in an extended area including the fluid and the surrounding air. The time‐dependent interface is captured by solving an appropriate level‐set function. The agreement between the two methods is fair. The influence of the stretching parameters (Draw ratio and drawing length) is investigated. For a long stretching distance, a critical Draw ratio around 20 delimitating stable and unstable drawing conditions is obtained, and this agrees well with the standard membrane models, which have been developed 40 years ago. When decreasing the stretching distance, the membrane model is no longer valid. The 2D models presented here point out a significant increase of the critical Draw ratio, and this is consistent with experimental results. Copyright © 2015 John Wiley & Sons, Ltd.     

Two-dimensional approach to the motion of a red blood cell in a plane couette flow of plasma

In relation to microrheology of blood, a theoretical approach to the motion of a red blood cell in a plane Couette flow between two parallel plates is made with emphasis on effects of wall. The red blood cell is assumed to be an elliptic cylindrical particle with a thin, inextensible membrane moving like a tank-tread along its perimeter and to contain a Newtonian fluid inside. Fluid motions are analysed numerically both inside and outside the particle on the basis of the Stokes equations, using the finite element method.A quasi-static equilibrium condition leads to the solution for the motion of the particle. It is shown that two types of motion exist (a stationary orientation motion and a flipping motion), depending on the viscosity ratio of inner to outer fluid, the axis ratio of the elliptic cylinder and the ratio of particle size to channel width. The results are applied to capillary blood flow.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid: 3D validation

A three‐dimensional finite element method for incompressible multiphase flows with capillary interfaces is developed based on a (formally) second‐order projection scheme. The discretization is on a fixed (Eulerian) reference grid with an edge‐based local h‐refinement in the neighbourhood of the interfaces. The fluid phases are identified and advected using the level‐set function. The reference grid is then temporarily reconnected around the interface to maintain optimal interpolations accounting for the singularities of the primary variables. Using a time splitting procedure, the convection substep is integrated with an explicit scheme. The remaining generalized Stokes problem is solved by means of a pressure‐stabilized projection. This method is simple and efficient, as demonstrated by a wide range of difficult free‐surface validation problems, considered in the paper. Copyright © 2005 John Wiley & Sons, Ltd.     

A substructure‐based iterative inner solver coupled with Uzawa's algorithm for the Stokes problem

A domain decomposition method with Lagrange multipliers for the Stokes problem is developed and analysed. A common approach to solve the Stokes problem, termed the Uzawa algorithm, is to decouple the velocity and the pressure. This approach yields the Schur complement system for the pressure Lagrange multiplier which is solved with an iterative solver. Each outer iteration of the Uzawa procedure involves the inversion of a Laplacian in each spatial direction. The objective of this paper is to effectively solve this inner system (the vector Laplacian system) by applying the finite‐element tearing and interconnecting (FETI) method. Previously calculated search directions for the FETI solver are reused in subsequent outer Uzawa iterations. The advantage of the approach proposed in this paper is that pressure is continuous across the entire computational domain. Numerical tests are performed by solving the driven cavity problem. An analysis of the number of outer Uzawa iterations and inner FETI iterations is reported. Results show that the total number of inner iterations is almost numerically scalable since it grows asymptotically with the mesh size and the number of subdomains. Copyright © 2003 John Wiley & Sons, Ltd.     

Integral vorticity transport method on an overset grid

We present an overset grid method for solution of the integro‐differential vorticity–velocity formulation of the Navier–Stokes equations for two‐dimensional, incompressible flow. The method uses a body‐fitted inner grid, on which vorticity is evolved semi‐implicitly, and a Cartesian outer grid with explicit vorticity evolution. The Biot–Savart integral is solved using an adaptive, optimized multipole acceleration method. The Biot–Savart integration is performed over all inner grid cells, over all ‘active cells’ of the outer grid that lie entirely outside of the inner grid, and over sub‐elements of a set of ‘overhanging’ cells of the outer grid that overlap part of the inner grid. A novel method is developed using a level‐set distance function to rapidly and easily partition the overhanging grid cells, which is essential for the Biot–Savart integration in order to avoid double‐counting vorticity in the overhanging region. A similar decomposition into outer, inner and overhanging cells is used in solving for pressure using a boundary‐element formulation, which requires evaluation of an integral over the vorticity field using a method similar to that used for the Biot–Savart integral. The new overset grid method is applied to flow past stationary and moving bodies in two dimensions and found to agree well with prior experimental and numerical results. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of surfactant on liquid film thickness in the Bretherton problem

The effects of insoluble and soluble surfactant on the motion of a long bubble propagating through a capillary tube are investigated computationally using a finite-difference/front-tracking method. Emphasis is placed on the effects of surfactant on the liquid film thickness between the bubble and the tube wall. The numerical method is designed to solve the evolution equations of the interfacial and bulk surfactant concentrations coupled with the incompressible Navier–Stokes equations. A non-linear equation of state is used to relate surface tension coefficient to surfactant concentration at the interface. Computations are first performed for soluble cases and then repeated for the corresponding clean and insoluble cases for a wide range of governing non-dimensional parameters in order to investigate the effects of surfactant and surfactant solubility. The computed film thickness for the clean case is found to be in a good agreement with Taylor’s law indicating the accuracy of the numerical method. We found that both the insoluble and soluble surfactant generally have a thickening effect on the film thickness, which is especially pronounced at low capillary numbers. This thickening effect strengthens with increasing sensitivity of surface tension to interfacial surfactant coverage mainly due to the enhanced Marangoni stresses along the liquid film. It is also observed that film thickening shows a non-monotonic behavior for variations in Peclet number. The validity of insoluble surfactant assumption is assessed for various non-dimensional numbers and it is demonstrated that insoluble assumption is valid only when capillary number is very low, i.e., Ca ≪ 1 and when surface tension is highly sensitive to interfacial surfactant coverage, i.e., the elasticity number is large.     

Finite element and implicit Runge–Kutta implementation of an acoustics–convection upstream resolution algorithm for the time‐dependent two‐dimensional Euler equations

The second of a two‐paper series, this paper details a solver for the characteristics‐bias system from the acoustics–convection upstream resolution algorithm for the Euler and Navier–Stokes equations. An integral formulation leads to several surface integrals that allow effective enforcement of boundary conditions. Also presented is a new multi‐dimensional procedure to enforce a pressure boundary condition at a subsonic outlet, a procedure that remains accurate and stable. A classical finite element Galerkin discretization of the integral formulation on any prescribed grid directly yields an optimal discretely conservative upstream approximation for the Euler and Navier–Stokes equations, an approximation that remains multi‐dimensional independently of the orientation of the reference axes and computational cells. The time‐dependent discrete equations are then integrated in time via an implicit Runge–Kutta procedure that in this paper is proven to remain absolutely non‐linearly stable for the spatially‐discrete Euler and Navier–Stokes equations and shown to converge rapidly to steady states, with maximum Courant number exceeding 100 for the linearized version. Even on relatively coarse grids, the acoustics–convection upstream resolution algorithm generates essentially non‐oscillatory solutions for subsonic, transonic and supersonic flows, encompassing oblique‐ and interacting‐shock fields that converge within 40 time steps and reflect reference exact solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

A 3D incompressible Navier–Stokes velocity–vorticity weak form finite element algorithm

The velocity–vorticity formulation is selected to develop a time‐accurate CFD finite element algorithm for the incompressible Navier–Stokes equations in three dimensions.The finite element implementation uses equal order trilinear finite elements on a non‐staggered hexahedral mesh. A second order vorticity kinematic boundary condition is derived for the no slip wall boundary condition which also enforces the incompressibility constraint. A biconjugate gradient stabilized (BiCGSTAB) sparse iterative solver is utilized to solve the fully coupled system of equations as a Newton algorithm. The solver yields an efficient parallel solution algorithm on distributed‐memory machines, such as the IBM SP2. Three dimensional laminar flow solutions for a square channel, a lid‐driven cavity, and a thermal cavity are established and compared with available benchmark solutions. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical implementation of the Crank–Nicolson/Adams–Bashforth scheme for the time‐dependent Navier–Stokes equations

This article considers numerical implementation of the Crank–Nicolson/Adams–Bashforth scheme for the two‐dimensional non‐stationary Navier–Stokes equations. A finite element method is applied for the spatial approximation of the velocity and pressure. The time discretization is based on the Crank–Nicolson scheme for the linear term and the explicit Adams–Bashforth scheme for the nonlinear term. Comparison with other methods, through a series of numerical experiments, shows that this method is almost unconditionally stable and convergent, i.e. stable and convergent when the time step is smaller than a given constant. Copyright © 2009 John Wiley & Sons, Ltd.     

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier–Stokes–Boussinesq equations governing the transport of momentum, mass and heat in a non-isothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a non-orthogonal co-ordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.     

Two‐level finite element method with a stabilizing subgrid for the incompressible Navier–Stokes equations

We consider the Galerkin finite element method for the incompressible Navier–Stokes equations in two dimensions. The domain is discretized into a set of regular triangular elements and the finite‐dimensional spaces employed consist of piecewise continuous linear interpolants enriched with the residual‐free bubble functions. To find the bubble part of the solution, a two‐level finite element method with a stabilizing subgrid of a single node is described, and its application to the Navier–Stokes equation is displayed. Numerical approximations employing the proposed algorithm are presented for three benchmark problems. The results show that the proper choice of the subgrid node is crucial in obtaining stable and accurate numerical approximations consistent with the physical configuration of the problem at a cheap computational cost. Copyright © 2008 John Wiley & Sons, Ltd.     

An ALE finite element method for a coupled Stefan problem and Navier–Stokes equations with free capillary surface

In this article, an ALE finite element method to simulate the partial melting of a workpiece of metal is presented. The model includes the heat transport in both the solid and liquid part, fluid flow in the liquid phase by the Navier–Stokes equations, tracking of the melt interface solid/liquid by the Stefan condition, treatment of the capillary boundary accounting for surface tension effects and a radiative boundary condition. We show that an accurate treatment of the moving boundaries is crucial to resolve their respective influences on the flow field and thus on the overall energy transport correctly. This is achieved by a mesh‐moving method, which explicitly tracks the phase boundary and makes it possible to use a sharp interface model without singularities in the boundary conditions at the triple junction. A numerical example describing the welding of a thin‐steel wire end by a laser, where all aforementioned effects have to be taken into account, proves the effectiveness of the approach.Copyright © 2012 John Wiley & Sons, Ltd.     

On the wrinkling of a pre-stressed annular thin film in tension

Asymptotic properties of the neutral stability curves for a linear boundary eigenvalue problem which models the wrinkling instability of an annular thin film in tension are considered. The film is subjected to imposed radial displacement fields on its inner and outer boundaries and, when these loads are sufficiently large, the film is susceptible to wrinkling. The critical values at which this onset occurs are dictated by the solution of a fourth-order ordinary differential eigensystem whose eigenvalue λ is a function of μ(?1), a quantity inversely proportional to the non-dimensional bending stiffness of the film, and n, the number of half-waves of the wrinkling pattern that sets in around the annular domain. Previously, Coman and Haughton [2006. Localised wrinkling instabilities in radially stretched annular thin films. Acta Mech. 185, 179-200] employed the compound matrix method together with a WKB technique to characterise the form of λ(μ,n) which essentially is related to a turning point in a reduced eigenproblem. The asymptotic analysis conducted therein pertained to the case when this turning point was not too close to the inner edge of the annulus. However, in the thin film limit μ→∞, the wrinkling load and the preferred instability mode are given by a modified eigenvalue problem that involves a turning point asymptotically close to the inner rim. Here WKB and boundary-layer asymptotic methods are used to examine these issues and comparisons with direct numerical simulations made.     

Numerical simulation of bubble and droplet deformation by a level set approach with surface tension in three dimensions

In this paper we present a three‐dimensional Navier–Stokes solver for incompressible two‐phase flow problems with surface tension and apply the proposed scheme to the simulation of bubble and droplet deformation. One of the main concerns of this study is the impact of surface tension and its discretization on the overall convergence behavior and conservation properties. Our approach employs a standard finite difference/finite volume discretization on uniform Cartesian staggered grids and uses Chorin's projection approach. The free surface between the two fluid phases is tracked with a level set (LS) technique. Here, the interface conditions are implicitly incorporated into the momentum equations by the continuum surface force method. Surface tension is evaluated using a smoothed delta function and a third‐order interpolation. The problem of mass conservation for the two phases is treated by a reinitialization of the LS function employing a regularized signum function and a global fixed point iteration. All convective terms are discretized by a WENO scheme of fifth order. Altogether, our approach exhibits a second‐order convergence away from the free surface. The discretization of surface tension requires a smoothing scheme near the free surface, which leads to a first‐order convergence in the smoothing region. We discuss the details of the proposed numerical scheme and present the results of several numerical experiments concerning mass conservation, convergence of curvature, and the application of our solver to the simulation of two rising bubble problems, one with small and one with large jumps in material parameters, and the simulation of a droplet deformation due to a shear flow in three space dimensions. Furthermore, we compare our three‐dimensional results with those of quasi‐two‐dimensional and two‐dimensional simulations. This comparison clearly shows the need for full three‐dimensional simulations of droplet and bubble deformation to capture the correct physical behavior. Copyright © 2009 John Wiley & Sons, Ltd.     

Flow simulation in the nip of a rigid forward roll coater

This paper shows the development of an efficient solution algorithm for the simulation of a forward roll coating flows with free surfaces. The technique is based on the method of successive approximation combined with a method for rapidly finding a good starting point, i.e. a good initial computational domain. Movement of the free surface from the initial to final positions is demonstrated using both kinematic and normal stress schemes. Adaptive domain decomposition at each pseudo time step is performed with no significant cost penalty. The flow field is computed using a finite element solution of the Navier–Stokes equation. The proposed scheme is shown to be flexible enough to accommodate different scenarios over the practical range of many applications, i.e. capillary numbers in the range of 0.01–300, and a wide range of gap settings and rotational speeds. A close match is found between simulation and experimental meniscus profiles as well as nip pressure profiles. Overall the technique is quite robust and is able to simulate a variety of coating flow situations without resorting to over‐simplifications. Copyright © 2005 John Wiley & Sons, Ltd.     

Asymptotic solutions for the Föppl – von Kármán equations governing deflections of thin axisymmetric annular plates

We are concerned with the deformation of thin, flat annular plates under a force applied orthogonally to the plane of the plate. This mechanical process can be described via a radial formulation of the Föppl – von Kármán equations, a set of nonlinear partial differential equations describing the deflections of thin flat plates. We are able to obtain analytical solutions for the radial Föppl – von Kármán equations with boundary conditions relevant for clamped, loosely clamped, and free inner and outer. This permits us to study the qualitative behavior of the out-of-plane deflections as well as the Airy stress function for a number of cases. Provided that an appropriate non-dimensionalization is taken, we find that the perturbation solutions are surprisingly valid for a wide variety of parameters, and compare favorably with numerical simulations in all cases (rather than just for small parameters). The results demonstrate that the ratio of the inner to outer radius of the annular plate will strongly influence the properties of the solutions, as will the specific boundary data considered. For instance, one may choose to fix the plate in place with a specific set of boundary conditions, in order to minimize the out-of-plane deflections. Other boundary conditions may result in undesirable behaviors.