Slow motions of a circular cylinder experiencing slip near a plane wall

A theoretical study is presented for the two-dimensional creeping flow caused by a long circular cylindrical particle translating and rotating in a viscous fluid near a large plane wall parallel to its axis. The fluid is allowed to slip at the surface of the particle. The Stokes equations for the fluid velocity field are solved in the quasi-steady limit using cylindrical bipolar coordinates. Semi-analytical solutions for the drag force and torque acting on the particle by the fluid are obtained for various values of the slip coefficient associated with the particle surface and of the relative separation distance between the particle and the wall. The results indicate that the translation and rotation of the confined cylinder are not coupled with each other. For the motion of a no-slip cylinder near a plane wall, our hydrodynamic drag force and torque results reduce to the closed-form solutions available in the literature. The boundary-corrected drag force and torque acting on the particle decrease with an increase in the slip coefficient for an otherwise specified condition. The plane wall exerts the greatest drag on the particle when its migration occurs normal to it, and the least in the case of motion parallel to it. The enhancement in the hydrodynamic drag force and torque on a translating and rotating particle caused by a nearby plane wall is much more significant for a cylinder than for a sphere.     

Slow motion of a slip spheroid along its axis of revolution

A combined analytical and numerical study of the Stokes flow caused by a rigid spheroidal particle translating along its axis of revolution in a viscous fluid is presented. The fluid is allowed to slip at the surface of the particle. The general solution for the stream function in prolate and oblate spheroidal coordinates can be expressed in an infinite-series form of semi-separation of variables. The slip boundary condition incorporating the shear stress at the particle surface is applied to this general solution to determine its unknown coefficients of the leading orders. The solution of these coefficients can be either numerical results obtained from a boundary-collocation method or explicit formulas derived analytically. The drag force exerted on the spheroidal particle by the fluid is evaluated with good convergence behavior for various values of the slip parameter and aspect ratio of the particle. The agreement between our hydrodynamic drag results and the relevant numerical solutions obtained previously using a singularity method is excellent. Although the drag force acting on the translating spheroid normalized by that on a corresponding sphere with equal equatorial radius increases monotonically with an increase in the axial-to-radial aspect ratio for a no-slip spheroid, it decreases monotonically as this aspect ratio increases for a perfect-slip spheroid. The normalized drag force exerted on a spheroid with a given surface slip coefficient in between the no-slip and perfect-slip limits is not a monotonic function of its aspect ratio. For a spheroid with a fixed aspect ratio, its drag force is a monotonically decreasing function of the slip coefficient of the particle.     

Creeping-flow rotation of a slip spheroid about its axis of revolution

The problem of rotation of a rigid spheroidal particle about its axis of revolution in a viscous fluid is studied analytically and numerically in the steady limit of negligible Reynolds number. The fluid is allowed to slip at the surface of the particle. The general solution for the fluid velocity in prolate and oblate spheroidal coordinates can be expressed in an infinite-series form of separation of variables. The slip boundary condition on the surface of the rotating particle is applied to this general solution to determine the unknown coefficients of the leading orders, which can be numerical results obtained from a boundary collocation method or explicit formulas derived analytically. The torque exerted on the spheroidal particle by the fluid is evaluated for various values of the slip parameter and aspect ratio of the particle. The agreement between our hydrodynamic torque results and the available analytical solutions in the limiting cases is good. It is found that the torque exerted on the rotating spheroid normalized by that on a sphere with radius equal to the equatorial radius of the spheroid increases monotonically with an increase in the axial-to-radial aspect ratio for a no-slip or finite-slip spheroid and vanishes for a perfectly slip spheroid. For a spheroid with a specified aspect ratio, the torque is a monotonically decreasing function of the slip capability of the particle.     

Particle and bubble dynamics in a creeping flow

Dynamics of a solid particle and non-deformable gaseous bubble in viscous fluid are studied analytically and numerically within the framework of creeping flow regime (flow at vanishingly small Reynolds numbers). Equations of motion for the particle and bubble include the consideration of the buoyancy force, Stokes drag force and memory-integral drag force. Exact analytical solutions are obtained and categorised in terms of inclusion (particle or bubble) density with respect to the density of a surrounding fluid. Through the analytical and numerical solutions, the dynamics of solid particle and air bubble in water have been found to behave differently especially at the early stages of motion, whereas some qualitative similarities exist in the long-term asymptotic.     

Hydrodynamic force and torque models for a particle moving near a wall at finite particle Reynolds numbers

This research work is aimed at proposing models for the hydrodynamic force and torque experienced by a spherical particle moving near a solid wall in a viscous fluid at finite particle Reynolds numbers. Conventional lubrication theory was developed based on the theory of Stokes flow around the particle at vanishing particle Reynolds number. In order to account for the effects of finite particle Reynolds number on the models for hydrodynamic force and torque near a wall, we use four types of simple motions at different particle Reynolds numbers. Using the lattice Boltzmann method and considering the moving boundary conditions, we fully resolve the flow field near the particle and obtain the models for hydrodynamic force and torque as functions of particle Reynolds number and the dimensionless gap between the particle and the wall. The resolution is up to 50 grids per particle diameter. After comparing numerical results of the coefficients with conventional results based on Stokes flow, we propose new models for hydrodynamic force and torque at different particle Reynolds numbers. It is shown that the particle Reynolds number has a significant impact on the models for hydrodynamic force and torque. Furthermore, the models are validated against general motions of a particle and available modeling results from literature. The proposed models could be used as sub-grid scale models where the flows between particle and wall can not be fully resolved, or be used in Lagrangian simulations of particle-laden flows when particles are close to a wall instead of the currently used models for an isolated particle.     

Axisymmetric creeping motion of a slip spherical particle in a nonconcentric spherical cavity

A semianalytical study of the creeping flow caused by a spherical fluid or solid particle with a slip surface translating in a viscous fluid within a spherical cavity along the line connecting their centers is presented in the quasisteady limit of small Reynolds number. In order to solve the Stokes equations for the flow field, a general solution is constructed from the superposition of the fundamental solutions in the two spherical coordinate systems based on both the particle and cavity. The boundary conditions on the particle surface and cavity wall are satisfied by a collocation technique. Numerical results for the hydrodynamic drag force exerted on the particle are obtained with good convergence for various values of the ratio of particle-to-cavity radii, the relative distance between the centers of the particle and cavity, the relative viscosity or slip coefficient of the particle, and the slip coefficient of the cavity wall. In the limits of the motions of a spherical particle in a concentric cavity and near a cavity wall with a small curvature, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force exerted on the particle for all cases is a monotonic increasing function of the ratio of particle-to-cavity radii, and becomes infinite in the touching limit. For a specified ratio of particle-to-cavity radii, the drag force is minimal when the particle is situated at the cavity center and increases monotonically with its relative distance from the cavity center to infinity in the limit as it is located extremely away from the cavity center. The drag force acting on the particle, in general, increases with an increase in its relative viscosity or with a decrease in its slip coefficient for a given configuration, but surprisingly, there are exceptions when the ratio of particle-to-cavity radii is large.     

Some Aspects of the Lift Force on a Spherical Bubble

Several situations in which a spherical bubble experiences a lift force are examined, especially through the use of computational results obtained by solving the full Navier–Stokes equations. The lift force is computed over a wide range of Reynolds number for the case of pure shear flow, pure strain and solid body rotation. Using these results, the validity of asymptotic solutions derived in the limit of low Reynolds number or inviscid flow is discussed. A general expression of the lift force valid for low to moderate shears is proposed. It is shown that for such shears, the lift force in a complex flow can be predicted by superposing the results obtained in pure strain flow and solid body rotation flow. Finally, the interaction force experienced by two bubbles rising side-by-side is studied. The computational results reveal that, at variance with the predictions of potential theory, the sign of this force changes when the Reynolds number or the separation distance between the bubbles decreases below a critical value. All these results are discussed in terms of vorticity. The respective role played by the vorticity generated at the bubble surface and by the one that is eventually present in the unperturbed flow is emphasized.     

Interception of two spheroidal particles in shear flow

A numerical method is developed for simulating the motion of particles with arbitrary shape in an effectively infinite or bounded viscous flow. The particle translational and angular velocities are computed directly by solving an integral equation of the second kind arising from the double-layer representation for Stokes flow, subject to a specified force and torque. The deflated integral equation is solved by the method successive substitutions using a spectral boundary-element method. In the case of force- and torque-free particles, the contribution of the particles to the effective viscosity of the suspension is expressed in terms of the distribution density of the double-layer potential over the particle surfaces. The method is applied to investigate the interception of two spheroidal particles in simple shear flow, with emphasis on the net particle displacement and shift in the phase of rotation after separation, and on the transient signature of the interception on the effective viscosity of the suspension. The net particle displacement normal to the streamlines of the shear flow is found to have both positive and negative values depending on the relative particle configuration before interception. For moderate particle separations, the phase of rotation is shifted only slightly with respect to the Jeffery orbit executed by a particle in isolation.     

Rotary oscillations of axi‐symmetric bodies in an axi‐symmetric viscous flow with slip: Numerical solutions for sphere and spheroids

Rotary oscillations of several axi‐symmetric bodies in axi‐symmetric viscous flows with slip are investigated. A numerical method based on the Green's function technique is used wherein the relevant Helmholtz equation, as obtained from the unsteady Stokes equation, is converted into a surface integral equation. The technique is benchmarked against a known analytical solution, and accurate numerical results for local stress and torque on spheres and spheroids as function of the frequency parameter and the slip coefficients are obtained. It is found that in all cases, slip reduces stress and torque, and increasingly so with the increasing frequency parameter. The method discussed here can be potentially extended to the realistic case of an oscillating disk viscometer. Copyright © 2003 John Wiley & Sons, Ltd.     

Creeping motion of a slip spherical particle in a circular cylindrical pore

A combined analytical–numerical study for the creeping flow caused by a spherical fluid or solid particle with a slip-flow surface translating in a viscous fluid along the centerline of a circular cylindrical pore is presented. To solve the axisymmetric Stokes equations for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the pore wall by the Fourier transforms and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force acting on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle, the slip parameter of the pore wall, and the ratio of radii of the particle and pore. For the motion of a fluid sphere along the axis of a cylindrical pore, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force for all cases is a monotonic increasing function of the ratio of particle-to-pore radii, and approaches infinity in the limit. Except for the case that the cylindrical pore is hardly slip and the value of the ratio of particle-to-pore radii is close to unity, the drag force exerted on the particle increases monotonically with an increase in its relative viscosity or with a decrease in its slip coefficient for a constant ratio of radii. In a comparison for the pore shape effect on the axial translation of a slip sphere, it is found that the particle in a circular cylindrical pore in general acquires a lower hydrodynamic drag than in a spherical cavity, but this trend can be reversed for the case of highly slippery particles and pore walls.     

Resistance of rough plates in a rotating fluid

Generalizing Navier’s partial slip condition, the flow due to a rough or striated plate moving in a rotating fluid is studied. It is found that the motion of the plate, the fluid surface velocity, and the shear stress are in general not in the same direction. The solution is extended to the case of finite depth, or Couette slip flow in a rotating system. In this case an optimum depth for minimum drag is found. The solutions are also closed form exact solutions of the Navier–Stokes equations. The results are fundamental to flows with Coriolis effects.     

The motion of rigid rod-like particles suspended in non-homogeneous flow fields

Equations are written for the velocities of rotation and translation of rigid rod-like particles suspended in arbitrary Stokes flows. These make use of the first approximation from slender body theory for the evaluation of drag forces parallel and transverse to the particle axis, and neglect couples induced by shear stress at the particle surface. They are therefore asymptotically valid as the particle axis ratio becomes large. Simple forms of the equations, applying in constant viscosity flows, are solved, where possible analytically and otherwise numerically, and results obtained for particle motion in planar Poiseuille and sink flows. These are discussed and displayed in terms of appropriate dimensionless groups in a comprehensive set of plots, that can conveniently be used to provide information on translational and rotational velocities, and orientation and displacement as a function of time, including particle slip along and across streamlines, for a wide range of cases. In this way the effects of non-homogeneity in the flow fields are quantified.     

Separation of screw-sensed particles in a homogeneous shear field

The Stokes motions of three-dimensional screw-sensed slender particles in a homogeneous shear field are investigated, including the effects of buoyancy. Conclusions are drawn about the possibility of achieving a separation of mixtures of right- and left-handed particles. The linearity of the Stokes equations allows complex flows to be solved by adding the effects of the several terms which describe the flow in which the particle is immersed. The homogeneous shear flow considered here consists of three such terms; solutions for a series of 12 unit motions are sufficient to determine the hydrodynamic resistance tensors. The forces and torques experienced by screw-sensed particles are calculated from these 51 resistance tensors, using slender-filament theory. The results allow an estimate of the range of buoyancy parameters for which gravitational sedimentation can be neglected. The fundamental component of the particle motion is a rotation, at approximately the same angular velocity as that of the fluid. Superimposed on this are variations, of large period, in the particle orientation. A phase plane analysis is used to find the terminal orientations. Very long calculation times are required for the phase portrait. An approximate method based on azimuthally-averaged equations is developed to avoid the requirements for long time integration.     

Interaction force and residual stress models for volume-averaged momentum equation for flow laden with particles of comparable diameter to computational grid width

An interaction model considering the local stress on a particle surface is developed based on a volume averaging technique. With a scope to apply to turbulence modulation caused by particles of diameter comparable to the Kolmogorov length scale, grid width for resolving vortical structures outside the boundary layer of the particles is set to be close to the particle diameter. The interaction force in the volume-averaged momentum equation is modelled based on analytical solutions of two fundamental flows. For the uniform flow case, the nonlinear effect of the first-order term in a series expansion with respect to the particle Reynolds number is found to be essential for the anisotropic Eulerian distribution of the interaction force. For the shear flow case, the anisotropic distribution of the interaction force is also essential, and it is modelled based on the Stokes’s solution. Considering that the length scale of the averaging volume is determined to be comparable to the grid width and the particle diameter, the residual stress, which originates from the volume averaging of the nonlinear term in the momentum equation, is also modelled based on an undisturbed linear shear flow. According to the test simulation using the interaction force and residual stress models of the fundamental flows, the anisotropic interaction force model essentially improves the representation of the flow field and the mechanical work, and the effect of the residual stress is found to be reasonably reproduced by the present model.     

Shear flow over a porous particle

Linear axisymmetric Stokes flow over a porous spherical particle is investigated. An exact analytic solution for the fluid velocity components and the pressure inside and outside the porous particle is obtained. The solution is generalized to include the cases of arbitrary three-dimensional linear shear flow as well as translational-shear Stokes flow. As the permeability of the particle tends to zero, the solutions obtained go over into the corresponding solutions for an impermeable particle. The problem of translational Stokes flow around a spherical drop (in the limit a gas bubble or an impermeable sphere) was considered, for example, in [1,2]. A solution of the problem of translational Stokes flow over a porous spherical particle was given in [3]. Linear shear-strain Stokes flow over a spherical drop was investigated in [2].Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 113–120, May–June, 1995.     

Drag forces of interacting spheres in power-law fluids

Drag forces of interacting particles suspended in power-law fluid flows were investigated in this study. The drag forces of interacting spheres were directly measured by using a micro-force measuring system. The tested particles include a pair of interacting spheres in tandem and individual spheres in a cubic matrix of multi-sphere in flows with the particle Reynolds number from 0.7 to 23. Aqueous carboxymethycellulose (CMC) solutions and glycerin solutions were used as the fluid media in which the interacting spheres were suspended. The range of power-law index varied from 0.6 to 1.0. In conjunction to the drag force measurements, the flow patterns and velocity fields of power-law flows over a pair of interacting spheres were also obtained from the laser assisted flow visualization and numerical simulation.

Both experimental and computational results suggest that, while the drag force of an isolated sphere depends on the power-index, the drag coefficient ratio of an interacting sphere is independent from the power-law index but strongly depends on the separation distance and the particle Reynolds number. Our study also shows that the drag force of a particle in an assemblage is strongly positions dependent, with a maximum difference up to 38%.     

Boundary conditions for stokes flows near a porous membrane

A theoretical development is carried out to model the boundary conditions for Stokes flows near a porous membrane, which, in general, allows non-zero slip as well as normal flow at the surface. Two types of models are treated: an infinitesimally thin plate with a periodic array of circular apertures and a series of parallel slits. For Stokes flows, the mean normal flux and slip velocity are proportional to the pressure difference across the membrane and the average shear stress at the membrane, respectively. The appropriate proportionality constants which depend on the membrane geometry are calculated as functions of the porosity. An interesting feature of the results is that the slip at the membrane has, in general, a direction different from that of the applied shear for these models.     

An SPH model for free surface flows with moving rigid objects

This paper presents a computational model for free surface flows interacting with moving rigid bodies. The model is based on the SPH method, which is a popular meshfree, Lagrangian particle method and can naturally treat large flow deformation and moving features without any interface/surface capture or tracking algorithm. Fluid particles are used to model the free surface flows which are governed by Navier–Stokes equations, and solid particles are used to model the dynamic movement (translation and rotation) of moving rigid objects. The interaction of the neighboring fluid and solid particles renders the fluid–solid interaction and the non‐slip solid boundary conditions. The SPH method is improved with corrections on the SPH kernel and kernel gradients, enhancement of solid boundary condition, and implementation of Reynolds‐averaged Navier–Stokes turbulence model. Three numerical examples including the water exit of a cylinder, the sinking of a submerged cylinder and the complicated motion of an elliptical cylinder near free surface are provided. The obtained numerical results show good agreement with results from other sources and clearly demonstrate the effectiveness of the presented meshfree particle model in modeling free surface flows with moving objects. Copyright © 2013 John Wiley & Sons, Ltd.     

Dynamics of a cavitation bubble near a solid surface and the induced damage

Numerical and experimental studies of the dynamics of a cavitating bubble near a resilient metal surface were performed. To augment the experimental flow visualizations of a collapsing bubble, numerical simulations were conducted to more thoroughly identify the collapse dynamics and analyze the flow. A bubble collapse was captured using a high-speed camera and back illumination. The metal sample was made of pure aluminum placed near a collapsing cavitation bubble at various distances from the metal surface. Width, depth, and volume of the induced material deformations were measured using an optical microscope and a three-dimensional profilometer and then compared against existing experimental data from the literature. The cavitating bubble’s dynamics and the related flow were simulated numerically using the open source finite volume based flow solver CavitatingFOAM. This code solved the Navier–Stokes equations for compressible two-phase flows using an Euler–Euler approach, including the barotropic equations of state. Bubble shapes, collapse times, and obtained damage parameters were compared to experimental observations. Impact velocities, pressures, shear rates, and various flow phenomena were discussed, providing broad insight into bubble dynamics and the induced damage.     

Slow Motion of a Porous Eccentric Spherical Particle-in-Cell Models

A combined analytical?Cnumerical method is presented for the quasisteady axisymmetrical flow of an incompressible viscous fluid past an assemblage of porous eccentric spherical particle-in-cell models. The flow inside the porous particle is governed by the Brinkman model and the flow in the fictitious envelope region is governed by Stokes equations. In order to solve the Stokes equations for the flow field, a general solution is constructed from the superposition of the basic solutions in the two spherical coordinate systems based on both the porous particle and fictitious spherical envelope. Boundary conditions on the particle??s surface and fictitious spherical envelope that correspond to the Happel, Kuwabara, Kvashnin, and Cunningham/Mehta-Morse models are satisfied by a collocation technique. The drag of these eccentric porous particles relative to the drag experienced by a centered porous particle are investigated as functions of the effective distance between the center of the porous particle and the fictitious envelope, the volume ratio of the porous particle over the surrounding sphere and a coefficient that is proportional to the inverse of the permeability. In the limits of the motions of the porous particle in the concentric position with cell surface and near the cell surface with a small curvature, the numerical values of the normalized drag force are in good agreement with the available values in the literature.