Sliding of fine particles on the slip surface of rising gas bubbles: Resistance of liquid shear flows

In this paper a model was developed to describe the shear flow resistance force and torque acting on a fine particle as it slides on the slip surface of a rising gas bubble. The shear flow close to the bubble surface was predicted using a Taylor series and the numerical data obtained from the Navier–Stokes equations as a function of the polar coordinates at the bubble surface, the bubble Reynolds number, and the gas hold-up. The particle size was considered to be sufficiently small relative to the bubble size that the bubble surface could be locally approximated to a planar interface. The Stokes equation for the disturbance shear flows was solved for the velocity components and pressure using series of bispherical coordinates and the boundary conditions at the no-slip particle surface and the slip bubble surface. The solutions for the disturbance flows were then used to calculate the flow resistance force and torque on the particle as a function of the separation distance between the bubble and particle surfaces. The resistance functions were determined by dividing the actual force and torque by the corresponding (Stokes) force and torque in the bulk phase. Finally, numerical and simplified analytical rational approximate solutions for force correction factors for sliding particles as a function of the (whole range of the) separation distance are presented, which are in good agreement with the exact numerical result and can be readily applied to more general modelling of the bubble–particle interactions.     

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.     

Asymptotic theory of the Boltzmann system,for a steady flow of a slightly rarefied gas with a finite Mach number: General theory

A steady rarefied gas flow with Mach number of the order of unity around a body or bodies is considered. The general behaviour of the gas for small Knudsen numbers is studied by asymptotic analysis of the boundary-value problem of the Boltzmann equation for a general domain. The effect of gas rarefaction (or Knudsen number) is expressed as a power series of the square root of the Knudsen number of the system. A series of fluid-dynamic type equations and their associated boundary conditions that determine the component functions of the expansion of the density, flow velocity, and temperature of the gas is obtained by the analysis. The equations up to the order of the square root of the Knudsen number do not contain non-Navier–Stokes stress and heat flow, which differs from the claim by Darrozes (in Rarefied Gas Dynamics, Academic Press, New York, 1969). The contributions up to this order, except in the Knudsen layer, are included in the system of the Navier–Stokes equations and the slip boundary conditions consisting of tangential velocity slip due to the shear of flow and temperature jump due to the temperature gradient normal to the boundary.     

Simulation of coating flows with slip effects

This work is concerned with the numerical prediction of wire coating flows. Both annular tube‐tooling and pressure‐tooling type extrusion–drag flows are investigated for viscous fluids. The effects of slip at die walls are analysed and free surfaces are computed. Flow conditions around the die exit are considered, contrasting imposition of no‐slip and various instances of slip models for die wall conditions. Numerical solutions are computed by means of a time marching Taylor–Galerkin/pressure–correction finite element scheme, that demonstrates how slip conditions on die walls mitigate stress singularities at the die exit. For pressure‐tooling and with appropriate handling of slip, reduction in shear rate at the die exit may be achieved. Maximum shear rates for tube‐tooling are about one quarter of those encountered in pressure‐tooling. Equivalently, extension rates peak at land entry, and tube‐tooling values are one third of those observed for pressure‐tooling. With slip and tube‐tooling, peak shear values at die exit may be almost completely eliminated. Nevertheless, in contrast to the pressure‐tooling scenario, this produces larger peak shear rates upstream within the land region than would otherwise be the case for no‐slip. Copyright © 2000 John Wiley & Sons, Ltd.     

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.     

The effect of the slip condition on Stokes and Couette flows due to an oscillating wall: exact solutions

Stokes and Couette flows produced by an oscillatory motion of a wall are analyzed under conditions where the no-slip assumption between the wall and the fluid is no longer valid. The motion of the wall is assumed to have a generic sinusoidal behavior. The exact solutions include both steady periodic and transient velocity profiles. It is found that slip conditions between the wall and the fluid produces lower amplitudes of oscillations in the flow near the oscillating wall than when no-slip assumption is utilized. Further, the relative velocity between the fluid layer at the wall and the speed of the wall is found to overshoot at a specific oscillating slip parameter or vibrational Reynolds number at certain times. In addition, it is found that wall slip reduces the transient velocity for Stokes flow while minimum transient effects for Couette flow is achieved only for large and small values of the wall slip coefficient and the gap thickness, respectively. The time needed to reach to steady periodic Stokes flow due to sine oscillations is greater than that for cosine oscillations with both wall slip and no-slip conditions.     

Tractions, Balances, and Boundary Conditions for Nonsimple Materials with Application to Liquid Flow at Small-Length Scales

Using a nonstandard version of the principle of virtual power, we develop general balance equations and boundary conditions for second-grade materials. Our results apply to both solids and fluids as they are independent of constitutive equations. As an application of our results, we discuss flows of incompressible fluids at small-length scales. In addition to giving a generalization of the Navier–Stokes equations involving higher-order spatial derivatives, our theory provides conditions on free and fixed boundaries. The free boundary conditions involve the curvature of the free surface; among the conditions for a fixed boundary are generalized adherence and slip conditions, each of which involves a material length scale. We reconsider the classical problem of plane Poiseuille flow for generalized adherence and slip conditions.     

BY-PASS MECHANISM OF TRANSITION TO TURBULENCE

The by-pass mechanism of transition for a wall-bounded shear layer is explained for the case when an infinite row of convecting vortices migrate over a boundary layer at a specific speed range. Such a mechanism is important for noisy flows over bluff bodies, flows inside turbomachinery and flows over helicopter rotor blades. By solving the Navier–Stokes equation, it is shown that this by-pass transition is a consequence of vortices migrating at convection speeds that are significantly lower than the free-stream speed. This situation is commonly found in flows that are affected by the presence of periodic wakes. Whenever the speed of migrating vortices is in a certain critical range, there is a local instability of the underlying shear layer with a very high-growth rate as compared to the growth of pure Tollmien–Schlichting waves created by wall excitation. The above interpretation is supported by solving the linearized and full Navier–Stokes equation for disturbance quantities under the parallel flow approximation in two dimensions. Various ramifications of such a by-pass route of transition are discussed in this paper.     

Empirical slip and viscosity model performance for microscale gas flow

For the simple geometries of Couette and Poiseuille flows, the velocity profile maintains a similar shape from continuum to free molecular flow. Therefore, modifications to the fluid viscosity and slip boundary conditions can improve the continuum based Navier–Stokes solution in the non‐continuum non‐equilibrium regime. In this investigation, the optimal modifications are found by a linear least‐squares fit of the Navier–Stokes solution to the non‐equilibrium solution obtained using the direct simulation Monte Carlo (DSMC) method. Models are then constructed for the Knudsen number dependence of the viscosity correction and the slip model from a database of DSMC solutions for Couette and Poiseuille flows of argon and nitrogen gas, with Knudsen numbers ranging from 0.01 to 10. Finally, the accuracy of the models is measured for non‐equilibrium cases both in and outside the DSMC database. Flows outside the database include: combined Couette and Poiseuille flow, partial wall accommodation, helium gas, and non‐zero convective acceleration. The models reproduce the velocity profiles in the DSMC database within an L₂ error norm of 3% for Couette flows and 7% for Poiseuille flows. However, the errors in the model predictions outside the database are up to five times larger. Copyright © 2005 John Wiley & Sons, Ltd.     

Particle dynamics calculations of wall stresses and slip velocities for couette flow of smooth inelastic spheres

The objective of this work is to characterize the effects of boundary geometry on the flow of dry granular materials composed of smooth, inelastic spheres between parallel, bumpy walls in the absence of gravity. Particle dynamic simulations are done in which wall stresses and slip velocities are computed over a wide range of parameters, including shear gap height, geometry of the wall particles, wall to free particle diameter ratio, and normal restitution coefficient. Calculated wall stresses and slip velocities are found to be highly sensitive to boundary roughness, which is characterized in terms of the mean spacing value for the wall particles. Most noticeable is a pronounced stress drop for dense flows with an associated large slip velocity in a system having a narrow shear gap height.     

Slip yield stress effects in start-up Newtonian Poiseuille flows

Analytical solutions are derived for various start-up Newtonian Poiseuille flows assuming that slip at the wall occurs when the wall shear stress exceeds a critical value, known as the slip yield stress. Two distinct regimes characterise the steady axisymmetric and planar flows, which are defined by a critical value of the pressure gradient. If the imposed pressure gradient is below this critical value, the classical no-slip, start-up solution holds. Otherwise, no-slip flow occurs only initially, for a finite time interval determined by a critical time, after which slip does occur. For the annular case, there is an additional intermediate (steady) flow regime where slip occurs only at the inner wall, and hence, there exist two critical values of the pressure gradient. If the applied pressure gradient exceeds both critical values, the velocity evolves initially with no-slip at both walls up to the first critical time, then with slip only along the inner wall up to the second critical time and finally with slip at both walls.     

Dispersion Equation for Water Waves with Vorticity and Stokes Waves on Flows with Counter-Currents

The two-dimensional free-boundary problem of steady periodic waves with vorticity is considered for water of finite depth. We investigate how flows with small-amplitude Stokes waves on the free surface bifurcate from a horizontal parallel shear flow in which counter-currents may be present. Two bifurcation mechanisms are described: one for waves with fixed Bernoulli’s constant, and the other for waves with fixed wavelength. In both cases the corresponding dispersion equations serve for defining wavelengths from which Stokes waves bifurcate. Necessary and sufficient conditions for the existence of roots of these equations are obtained. Two particular vorticity distributions are considered in order to illustrate the general results.     

Various aspects of fluid vortices

A general review of the current research in vortex dynamics is presented, based on contributions given during a workshop held in May 2003 at Porquerolles, France. This article aims at providing a picture of the work performed on this subject in the French community. Various cases are covered, from 2D vortex patches to 3D vortex tubes; from isolated vortices to shear flows. Different contexts are considered: pure Euler and Navier–Stokes flows as well as stratified, rotating and magnetic flows. To cite this article: I. Delbende et al., C. R. Mecanique 332 (2004).     

A Dynamic-induced Direct-shear Model for Dynamic Triggering of Frictional Slip on Simulated Granular Gouges

This study presents a dynamic-induced direct-shear model to investigate the dynamic triggering of frictional slip on simulated granular gouges. An incident P-wave is generated as a shear load and a normal stress is constantly applied on the gouge layer. The shear stress accumulates in the incident stage and the frictional slip occurs in the slip stage without the effect of the reflected wave. The experimental results show a non-uniform shear stress distribution along the gouge layer, which may be induced by a shear load induced torque and by normal stress vibration along the layer. The shear stress at the trailing edge strongly affects the frictional slip along the P-wave loading direction, while the rebound stress at the leading edge propagates along the opposite direction. The frictional slip is triggered when the maximum shear stress at the trailing edge reaches a critical value. The normal stress influences the maximum shear stress at the trailing edge, the maximum slip displacement and the slip velocity. The advantages and the limitations of this model are discussed at the end.     

General results for complete contacts subject to oscillatory shear

In this paper, we consider the general interfacial characteristics of a square elastic block, pressed onto an elastically similar half-plane by a constant normal force, and subjected to oscillatory shear. It is found that there is a critical coefficient of friction, 0.543, above which the contact is permanently stuck along its entire length for a shearing force below about 55% of that needed to cause sliding. For shearing forces above this, the contact interface will either shakedown to a fully adhered state (depending on the degree of reversal of the shear loading) or will exhibit cyclic slip at an interior point. If the coefficient of friction is below 0.543, the application of normal load alone will produce equal and opposite slip zones attached to the contact edges. The subsequent imposition of a shear force causes the leading edge slip zone to increase in length while the presence of residual slipping tractions at the trailing edge causes the trailing edge to lift off. Under oscillatory loading, the contact edges cycle between slip and separation over a minute region while an interior point may exhibit cyclic slip if the loading history is sufficiently demanding. The results found are of practical relevance to the study of fretting fatigue of complete contacts, such as some types of spline joint.     

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.     

Self-similar regimes of liquid-layer spreading along a superhydrophobic surface

Within the Stokes film approximation, unsteady spreading of a thin layer of a heavy viscous fluid along a horizontal superhydrophobic surface is studied in the presence of a given localized mass supply in the film. The forced (induced by the mass supply) spreading regimes are considered, for which the surface tension effects are insignificant. Plane and axisymmetric flows along the principal direction of the slip tensor of the superhydrophobic surface are studied, when the corresponding slip tensor component is either a constant or a power function of the spatial coordinate, measured in the direction of spreading. An evolution equation for the film thickness is derived. It is shown that this equation has self-similar solutions of a source type. The examples of self-similar solutions are constructed for power and exponential time dependences of mass supply. In the final part of the paper, some of the solutions constructed are generalized to the case of a weak dependence of the flow on the second spatial coordinate, caused by a slight variability of the slip coefficient in the direction normal to that of spreading. The constructed self-similar solutions can be used for experimental determination of the parameters important for hydrodynamics, e.g. the slip tensor components of commercial superhydrophobic surfaces.     

Comparison of low Mach number models for natural convection problems

 We investigate in this paper two numerical methods for solving low Mach number compressible flows and their application to single-phase natural convection flow problems. The first method is based on an asymptotic model of the Navier–Stokes equations valid for small Mach numbers, whereas the second is based on the full compressible Navier–Stokes equations with particular care given to the discretization at low Mach numbers. These models are more general than the Boussinesq incompressible flow model, in the sense that they are valid even for cases in which the fluid is subjected to large temperature differences, that is when the compressibility of the fluid manifests itself through low Mach number effects. Numerical solutions are computed for a series of test problems with fixed Rayleigh number and increasing temperature differences, as well as for varying Rayleigh number for a given temperature difference. Numerical difficulties associated with low Mach number effects are discussed, as well as the accuracy of the approximations. Received on 17 January 2000     

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.     

Simultaneous solution of the Navier-Stokes and elastic membrane equations by a finite element method

Fluid flow through a significantly compressed elastic tube occurs in a variety of physiological situations. Laboratory experiments investigating such flows through finite lengths of tube mounted between rigid supports have demonstrated that the system is one of great dynamical complexity, displaying a rich variety of self-excited oscillations. The physical mechanisms responsible for the onset of such oscillations are not yet fully understood, but simplified models indicate that energy loss by flow separation, variation in longitudinal wall tension and propagation of fluid elastic pressure waves may all be important. Direct numerical solution of the highly non-linear equations governing even the most simplified two-dimensional models aimed at capturing these basic features requires that both the flow field and the domain shape be determined as part of the solution, since neither is known a priori. To accomplish this, previous algorithms have decoupled the solid and fluid mechanics, solving for each separately and converging iteratively on a solution which satisfies both. This paper describes a finite element technique which solves the incompressible Navier-Stokes equatikons simultaneously with the elastic membrane equations on the flexible boundary. The elastic boundary position is parametized in terms of distances along spines in a manner similar to that which has been used successfully in studies of viscous free surface flows, but here the membrane curvature equation rather than the kinematic boundary condition of vanishing normal velocity is used to determine these diatances and the membrane tension varies with the shear stresses exerted on it by the fluid motions. Bothy the grid and the spine positions adjust in response to membrane deformation, and the coupled fluid and elastic equations are solved by a Newton-Raphson scheme which displays quadratic convergence down to low membrane tensions and extreme states of collapse. Solutions to the steady problem are discussed, along with an indication of how the time-dependent problem might be approached.