Inversion transformation of Stokes flows bounded by parallel planes

The positive inversion transformation applied to a two-dimensional Stokes flow around bodies leads alike to a Stokes flow. This fact can be exploited to find new two-dimensional Stokes flow solutions around inverse bodies. Some features of this method, such as the relations between the reference and inverse fluid velocity fields, are presented followed by an application to examples of cellular flow between two parallel plates induced by rotating or translating cylinder. Thus hydrodynamic characteristics of flow around circular bodies obtained by inversion of the plates are straightforward deduced. Typical fluid flow patterns around two circular cylinders in contact placed in the centre of a rotating or a translating circular cylinder are thus illustrated.     

Circle theorems for steady Stokes flow

Two circle theorems for two-dimensional steady Stokes flow are presented. The first theorem gives an expression for the stream function for a Stokes flow past a circular cylinder in terms of the stream function for a slow and steady irrotational flow in an unbounded incompressible viscous fluid. The second theorem gives a more general expression for the stream function for another Stokes flow past the circular cylinder in terms of the stream function for a slow and steady rotational flow in the same fluid.     

Controlling chaos in a fluid flow past a movable cylinder

The model of a two-dimensional fluid flow past a cylinder is a relatively simple problem with a strong impact in many applied fields, such as aerodynamics or chemical sciences, although most of the involved physical mechanisms are not yet well known. This paper analyzes the fluid flow past a cylinder in a laminar regime with Reynolds number, Re, around 200, where two vortices appear behind the cylinder, by using an appropriate time-dependent stream function and applying non-linear dynamics techniques. The goal of the paper is to analyze under which circumstances the chaoticity in the wake of the cylinder might be modified, or even suppressed. And this has been achieved with the help of some indicators of the complexity of the trajectories for the cases of a rotating cylinder and an oscillating cylinder.     

A note on plane Stokes flow past a shear free impermeable cylinder

Summary The slow steady two-dimensional motion of a viscous incompressible fluid in the unbounded region exterior to a shear free circular cylinder which is impermeable is examined. It is shown that the above problem requires a certain consistency condition for the existence of a solution. In addition, a circle theorem for the biharmonic equation is presented, for the above plane Stokes flow. Some examples are also given.     

Numerical solution of the free‐surface viscous flow on a horizontal rotating elliptical cylinder

The numerical solution of the free‐surface fluid flow on a rotating elliptical cylinder is presented. Up to the present, research has concentrated on the circular cylinder for which steady solutions are the main interest. However, for noncircular cylinders, such as the ellipse, steady solutions are no longer possible, but there will be periodic solutions in which the solution is repeated after one full revolution of the cylinder. It is this new aspect that makes the investigation of noncircular cylinders novel. Here we consider both the time‐dependent and periodic solutions for zero Reynolds number fluid flow. The numerical solution is expedited by first mapping the fluid film domain onto a rectangle such that the position of the free‐surface is determined as part of the solution. For the time‐dependent case a simple time‐marching method of lines approach is adopted. For the periodic solution the discretised nonlinear equations have to be solved simultaneously over a time period. The resulting large system of equations is solved using Newton's method in which the form of the Jacobian enables a straightforward decomposition to be implemented, which makes matrix inversion manageable. In the periodic case all derivatives have been approximated pseudospectrally with the time derivative approximated by a differentiation matrix which has been specially derived so that the weight of fluid is algebraically conserved. Of interest is the solution for which the weight of fluid is at its maximum possible value, and this has been obtained by increasing the weight until a consistency break‐down occurs. Time‐dependent solutions do not produce the periodic solution after a long time‐scale but have protuberances which are constantly appearing and disappearing. Periodic solutions exhibit spectral accuracy solutions and maximum supportable weight solutions have been obtained for ranges of eccentricity and angular velocity. The maximum weights are less than and approximately proportional to those obtained for the circular case. The shapes of maximum weight solutions is distinctly different from sub‐maximum weight solutions. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Near trapped wave phenomena in an array of truncated cylinders in a perpendicular arrangement in front of a vertical breakwater

The present paper deals with the investigation of the interaction phenomena of regular waves with arrays of circular cylinders either bottom-mounted or floating in front of a vertical breakwater of infinite and finite length. The acting loads on the array have been modeled using an analytical solution to the diffraction problem. The method of images is applied to define the fluid flow around the bodies of the array in front of a breakwater of infinite length, whereas the finite-length breakwater is represented by an elliptical cylinder with zero semi–minor axis. The multiple scattering approach and the direct matrix inversion method are compared to describe the hydrodynamic interaction phenomena among the members of the array. Focus is placed on the so called “near trapped” mode phenomena, created by an array of finite length, associated to the presence of the breakwater supplemented by representative numerical results concerning the exciting loads on the cylindrical bodies.     

Simulation of the flow past a circular cylinder using an unsteady panel method

In the present work, an in-house UnSteady Double Wake Model (USDWM) is developed for simulating general flow problems behind bodies. The model is presented and used to simulate flows past a circular cylinder at subcritical, supercritical, and transcritical flows. The flow model is a two-dimensional panel method which uses the unsteady double wake technique to model flow separation and its dynamics. In the present work the separation location is obtained from experimental data and fixed in time. The highly unsteady flow field behind the cylinder is analyzed in detail. The results are compared with experiments and Unsteady Reynolds-Averaged Navier Stokes (URANS) simulations and show good agreement in terms of the vortex shedding characteristics, drag, and pressure coefficients for the different flow regimes.     

Expansions at small Reynolds numbers for the flow past a porous circular cylinder

The purpose of this article is to use the method of matched asymptotic expansions (MMAE) in order to study the two-dimensional steady low Reynolds number flow of a viscous incompressible fluid past a porous circular cylinder. We assume that the flow inside the porous body is described by the continuity and Brinkman equations, and the velocity and boundary traction fields are continuous across the interface between the fluid and porous media. Formal expansions for the corresponding stream functions are used. We show that the force exerted by the exterior flow on the porous cylinder admits an asymptotic expansion with respect to low Reynolds numbers, whose terms depend on the characteristics of the porous cylinder. In addition, by considering Darcy's law for the flow inside the porous circular cylinder, an asymptotic formula for the force on the cylinder is obtained. Also, a porous circular cylinder with a rigid core inside is considered with Brinkman equation inside the porous region. Stress jump condition is used at the porous–liquid interface together with the continuity of velocity components and continuity of normal stress. Some particular cases, which refer to the low Reynolds number flow past a solid circular cylinder, have also been investigated.     

Nonlinear characteristics of vortex-induced vibration at low Reynolds number

Numerical simulations of vortex-induced vibration of a two-dimensional elastic circular cylinder under the uniform flow are calculated when Reynolds number is 200. In order to achieve the vortex-induced vibration, two-dimensional incompressible Navier–Stokes equations are solved with the space–time finite element method, the equations of the cylinder motion are solved with the new explicit integral method and the emeshing is achieved by the spring analogy technology. Considering vortex-induced vibration with the low reduced damping parameters, the variety trends of the lift coefficient, the drag coefficient, the displacement of cylinder are analyzed under different oscillating frequencies of cylinder. The nonlinear phenomena of locked-in, beat and phaseswith are captured successfully. The limit cycle and bifurcation of lift coefficient and displacement are analyzed. Besides, the Poincare sections of the lift coefficient are used for discussing the bifurcation of periodic solution. There are some differences in nonlinear characteristics between the results of the one degree of freedom cylinder model and those of the two degrees of freedom cylinder model. The streamwise vibration has a certain effect on the lateral vibration.     

A note on the orientation of symmetric rigid bodies sedimenting in a power-law fluid

We study the terminal orientation of symmetric bodies translating in a quiescent liquid modeled by the power-law fluid. We are able to show by invoking the symmetries of the sedimenting body and the Stokes flow field that at small Reynolds numbers, the competition of inertial and shear-thinning (or shear-thickening) contributions to the torque does not cause the tilt angle that is observed in experiments performed on viscoelastic liquids with shear-thinning properties.     

Axial eigenfunctions and inhomogeneous solutions for the three-dimensional Stokes system with an application to natural convection in a cylinder

This paper presents two contributions to the analysis of three-dimensional slow viscous flows in cylinders of circular section. First the vector axial eigenfunctions for this geometry, namely those that satisfy homogeneous boundary conditions on the flat end walls, are derived. Secondly a method is presented to find particular solutions to the inhomogeneous Stokes equations in this geometry. These new results, together with some results obtained earlier, are used to analyse slow natural convection in a vertical cylinder completely filled with a viscous liquid. The fluid motion is generated by the differential heating of the walls of the cylinder. The natural convection flow field is shown to be a superposition of an inhomogeneous field, the fields generated by the vector eigenfunctions and a Stokes flow field. A by-product of this work has been the identification of constraints on the boundary data that have to be satisfied in order for the eigenfunction expansions to work; this knowledge will be useful when attempts are made to prove the completeness of these Stokes flow eigenfunctions.Received: June 30, 2003; revised: February 26, 2004     

Axial eigenfunctions and inhomogeneous solutions for the three-dimensional Stokes system with an application to natural convection in a cylinder

This paper presents two contributions to the analysis of three-dimensional slow viscous flows in cylinders of circular section. First the vector axial eigenfunctions for this geometry, namely those that satisfy homogeneous boundary conditions on the flat end walls, are derived. Secondly a method is presented to find particular solutions to the inhomogeneous Stokes equations in this geometry. These new results, together with some results obtained earlier, are used to analyse slow natural convection in a vertical cylinder completely filled with a viscous liquid. The fluid motion is generated by the differential heating of the walls of the cylinder. The natural convection flow field is shown to be a superposition of an inhomogeneous field, the fields generated by the vector eigenfunctions and a Stokes flow field. A by-product of this work has been the identification of constraints on the boundary data that have to be satisfied in order for the eigenfunction expansions to work; this knowledge will be useful when attempts are made to prove the completeness of these Stokes flow eigenfunctions.     

Exact solutions of problems of the three-dimensional flow of ideal viscous incompressible fluids near cylindrical surfaces

Three-dimensional flows of an incompressible fluid, the parameters of which depend on two coordinates and time, are considered. The stream surfaces of such flows are cylindrical. The equations of continuity and the Navier-Stokes equations can be transformed to relations, one of which is the equation for the stream function the other is the integral of the equations relating the pressure and the stream function, and the third is a linear equation for the projection of the velocity vector onto the axis parallel to the generatrix of the cylindrical surfaces. The problems of modelling the flows are considered on the basis of the exact solutions of the Navier-Stokes equations and Euler's equations using examples. Relations for the distribution of the flow parameters in the channel created by hyperbolical cylinders are derived for the case of unsteady inviscid flow. The streamlines of these flows are situated on the side surfaces of the hyperbolical cylinders and intercept the generatrices of the cylinders at certain indirect angles. The flow around a circular cylinder and the flow of fluid inside an elliptic cylinder are considered in the case of steady inviscid flow. The streamlines on the circular cylinder are arranged transverse to the cylinder (the projection of the velocity vector onto the coordinate axis, parallel to the generatrix of the cylinder, is equal to zero). Far from the cylinder the streamlines are also situated on a cylindrical surfaces, but not transverse to the cylinder, making certain indirect angles with the generatrix. Viscous three-dimensional flows, possessing a certain symmetry, are considered. In the case of radial symmetry the streamlines are helical lines. The non-planar Couette flow between parallel moving planes is characterized by the fact that the velocity vectors, being situated in the same plane, are collinear, while the velocity vectors in parallel planes are not collinear. Relations for viscous steady three-dimensional flows, using well-known relations, obtained for the stream function of two-dimensional flows, are given.     

MILU-CG方法和绕旋转圆柱流动的数值研究

本文采用以修正的不完全LU分解作预处理器的共轭梯度法(MILU-CG),结合高阶隐式差分格式,改进了作者(1992)提出的基于区域分解、有限差分法与涡法杂交的数值方法(HDV).系统地研究了雷诺数Re=1000,200,旋转速度比α∈(0.5,3.25)范围内,绕旋转圆柱从突然起动到充分发展,长时间内尾流旋涡结构和阻力、升力系数的变化规律.计算所得流线与实验流场显示相比,完全吻合.首次揭示了临界状态时的旋涡结构特性,并指出最佳升阻比就在该状态附近得到.     

Drag force and virtual mass of a cylindrical porous shell in potential flow

We study the two-dimensional potential flow due to a circular cylinder in motion relative to an unbounded fluid. The cylinder consists of a thin, circular porous shell with fluid inside. The full nonlinear hydrodynamic problem is solved by Fourier expansion of Green's theorem. The truncated series is determined numerically by sampling points around the circle. A dimensionless shell parameter is introduced. For homogeneous porous shells, a maximal drag force occurs at the value 0.433 for the shell parameter, but the virtual mass is a monotonous function of the shell parameter. For an inhomogeneous shell, we have found a maximal value for the virtual mass which is 5% above the value for a rigid cylinder. Some of the results may be relevant to offshore engineering, especially in connection with porous coating of platform legs to reduce the total force.     

Modelling two-dimensional flow past arbitrary cylindrical bodies using boundary element formulations

This paper presents an efficient method of solving Queen's linearized equations for steady plane flow of an incompressible, viscous Newtonian fluid past a cylindrical body of arbitrary cross-section. The numerical solution technique is the well known direct boundary element method. Use of a fundamental solution of Oseen's equations, the ‘Oseenlet’, allows the problem to be reduced to boundary integrals and numerical solution then only requires boundary discretization. The formulation and solution method are validated by computing the net forces acting on a single circular cylinder, two equal but separated circular cylinders and a single elliptic cylinder, and comparing these with other published results. A boundary element representation of the full Navier-Stokes equations is also used to evaluate the drag acting on a single circular cylinder by matching with the numerical Oseen solution in the far field.     

速度-涡量法数值求解具有表面吹吸圆柱的绕流问题

用速度-涡量法数值求解了具有表面吹吸圆柱的绕流问题.所得高阶隐式差分方程,采用以修正的不完全LU分解作预处理器的共轭梯度法(MILU-CG),高效解出.研究了雷诺数Re=100时,各种吹吸位置、吹吸强度对圆柱尾流涡旋结构和阻力、升力系数的影响规律.指出,在圆柱肩部的吸气和在圆柱尾部的吹气,可有效地抑制尾流涡旋结构在垂直来流方向上的非对称性,达到减小升力的目的.对在圆柱肩部吸气的情形,合适选择吸气强度,还可有效减小圆柱在来流方向上所受的阻力.     

Couette flow between corrugated cylinders

This note reports some analytic Stokes solutions to the flows outside a rotating corrugated cylinder and between two rotating corrugated cylinders of either epitrochoidal or hypotrochoidal cross-sections. The results are obtained by using conformal mappings and general biharmonic solutions inside and outside of a unit circle.     

Effect of aligned magnetic field on the steady viscous flow past a circular cylinder

The steady viscous incompressible and slightly conducting fluid flow around a circular cylinder with an aligned magnetic field is simulated for the range of Reynolds numbers 100 ? Re ? 500 using the Hartmann number, M. The multigrid method with defect correction technique is used to achieve the second order accurate solution of complete non-linear Navier–Stokes equations. The magnetic Reynolds number is assumed to be small. It is observed that volume of the separation bubble decreases and drag coefficient increases as M is increased. We noticed that the upstream base pressure increases slightly with increase of M whereas downstream base pressure decreases with increase of M. The effect of the magnetic field on the flow is discussed with contours of streamlines, vorticity, plots of surface pressure and surface vorticity.     

Separation of streamlines for spatially periodic flow at non-zero Reynolds numbers

The flow generated by a small rotating circular cylinder at the center of a corrugated outer cylinder is considered. By using a Stokes expansion, the first order correction in the Reynolds numberR is found for the creeping flow solution. An approximate critical Reynolds numberR _c is found at which separation appears, and it is expressed in terms of the boundary parameters. Separation is found to occur in the concave regions of the boundary skewed opposite to the direction of rotation of the inner cylinder. By partially solving for the second order correction in the Stokes expansion, it is found that an increase inR causes an increase in the torque exerted on the outer boundary.This work was supported in part by a grant from NSERC.