Modeling incompressible flows using a finite particle method

This paper describes the applications of a finite particle method (FPM) to modeling incompressible flow problems. FPM is a meshfree particle method in which the approximation of a field variable and its derivatives can be simultaneously obtained through solving a pointwise matrix equation. A set of basis functions is employed to obtain the coefficient matrix through a sequence of transformations. The finite particle method can be used to discretize the Navier–Stokes equation that governs fluid flows. The incompressible flows are modeled as slightly compressible via specially selected equations of state. Four numerical examples including the classic Poiseuille flow, Couette flow, shear driven cavity and a dam collapsing problem are presented with comparisons to other sources. The numerical examples demonstrate that FPM is a very attractive alternative for simulating incompressible flows, especially those with free surfaces, moving interfaces or deformable boundaries.     

Large shear rate behavior for the Hébraud-Lequeux model

The Hébraud-Lequeux model is a model describing the flow of soft glassy material in a simple shear flow configuration.It is given by a kinetic/Fokker-Planck-type equation whose coefficients depend on the shear rate of the experiment.In this paper we want to study what happens to the stationary solutions of this model when the shear rate is asymptotically large.In order to do that,we expand the solution of the equation using singular perturbation tools.In the end,we rigorously prove the estimate of Hébraud and Lequeux that the material asymptotically behaves as a Newtonian fluid.     

非牛顿流体内流传热研究

在本文中,研究了注入轴对称模腔非牛顿流体非定常流动.本文的第二部份研究了上随体Maxwell流体管内热流动.对于注入模腔流动.其本构方程采用幂律流体模型方程.为了避免在表现粘度中温度关系引起的非线性.引进了一特征粘度的概念.描述本力学过程的基本方程是,本构方程、定常状态的运动方程、非定常能量方程及连续方程.该方程组在空间是二维问题,在数学上是三维问题.采用分裂差分格式求得本方程组的数值解答.分裂法曾成功应用于求解牛顿流体问题.在本文中,首次将分裂法成功地应用解决非牛顿流体流动问题.对于圆管内热流,给出了差分格式,使基本方程组化为一个三对角方程组.其结果,给出了不同时刻的模腔内二维温度分布.     

Flow of polymer melts in a rotary extruder head

The flow of a non-Newtonian fluid in an annular extruder channel formed by two coaxial cylinders is investigated theoretically. Complex shear conditions are created by the pressure difference between the ends of the channel and the rotation of the cylinders in opposite directions at constant speeds. An expression is obtained for the rate of flow of a polymer melt regarded as a non-Newtonian fluid. The flow of a polyethylene melt at 140° C through an extruder head with a rotating mandrel has been experimentally investigated. The experimental data serve to confirm the theoretical conclusions.Mekhanika Polimerov, Vol. 4, No. 3, pp. 531–539, 1968     

Finite amplitude long-wave instability of a thin viscoelastic fluid during spin coating

This paper investigates the stability of a thin incompressible viscoelastic fluid designated as Walters’ liquid B″ during spin coating. The long-wave perturbation method is proposed to derive a generalized kinematic model of the film flow. The method of normal mode is applied to study the linear stability. The amplitude growth rates and the threshold conditions are characterized subsequently and summarized as the by-products of the linear solutions. Using the multiple scales method, the weakly nonlinear stability analysis is studied for the evolution equation of a film flow. The Ginzburg–Landau equation is determined to discuss the threshold conditions of the various critical flow states. The study reveals that the rotation number and the radius of the rotating circular disk generate the destabilizing effects. Moreover, the viscoelastic parameter k indeed plays a more significant role in destabilizing the film flow than a thin Newtonian fluid during spin coating [27].     

Influence of Mass Loading on Particle-Laden Turbulent Pipe Flow

A CFD code in the framework of OpenFOAM was validated for simulations of particle-laden pipe and channel flows at low to intermediate mass loadings. The code is based on an Eulerian two-fluid approach with Reynolds-averaged conservation equations, including turbulence modeling and four-way coupling. Pipe flow simulations of particles in air against gravity were conducted at Reynolds numbers up to 50000. The particle mass loading was varied and its effect on the mean velocities and turbulent fluctuations of the two phases was studied. Special attention was paid to the influence of mass loading on the centerline velocity and the wall shear velocity of the fluid phase for various flow parameters and particle properties. Empirical correlations were established between these two quantities and the flow Reynolds number, particle Reynolds number, Stokes number and particle to fluid density ratio for a range of particle mass loadings. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Weakly dissipative solutions and weak–strong uniqueness for the Navier–Stokes–Smoluchowski system

This article deals with a fluid–particle interaction model for the evolution of particles dispersed in a fluid. The fluid flow is governed by the Navier–Stokes equations for a compressible fluid while the evolution of the particle densities is given by the Smoluchowski equation. The coupling between the dispersed and dense phases is obtained through the drag forces that the fluid and the particles exert mutually. The existence of weakly dissipative solutions is established under reasonable physical assumptions on the initial data, the physical domain, and the external potential. Furthermore, a weak–strong uniqueness result is established via the relative entropy method yielding that a weakly dissipative solution agrees with a classical solution with the same initial data when such a classical solution exists.     

A note on the unsteady flow of a generalized second-grade fluid through a circular cylinder subject to a time dependent shear stress

The unsteady flow of a generalized second-grade fluid through an infinite straight circular cylinder is considered. The flow of the fluid is due to the longitudinal time dependent shear stress that is prescribed on the boundary of the cylinder. The fractional calculus approach in the governing equation corresponding to a second-grade fluid is introduced. The velocity field and the resulting shear stress are obtained by means of the finite Hankel and Laplace transforms. In order to avoid lengthy calculations of residues and contour integrals, the discrete inverse Laplace transform method is used. The corresponding solutions for ordinary second-grade and Newtonian fluids, performing the same motion, are obtained as limiting cases of our general solutions. Finally, the influence of the material constants and of the fractional parameter on the velocity and shear stress variations is underlined by graphical illustrations.     

Coloured noise for dispersion of contaminants in shallow waters

In this article, we explore the application of a set of stochastic differential equations called particle model in simulating the advection and diffusion of pollutants in shallow waters. The Fokker–Planck equation associated with this set of stochastic differential equations is interpreted as an advection–diffusion equation. This enables us to derive an underlying particle model that is exactly consistent with the advection–diffusion equation. Still, neither the advection–diffusion equation nor the related traditional particle model accurately takes into account the short-term spreading behaviour of particles. To improve the behaviour of the model shortly after the deployment of contaminants, a particle model forced by a coloured noise process is developed in this article. The use of coloured noise as a driving force unlike Brownian motion, enables to us to take into account the short-term correlated turbulent fluid flow velocity of the particles. Furthermore, it is shown that for long-term simulations of the dispersion of particles, both the particle due to Brownian motion and the particle model due to coloured noise are consistent with the advection–diffusion equation.     

Motion of a distant solid particle in a shear flow along a porous slab

The motion of a solid and no-slipping particle immersed in a shear flow along a sufficiently porous slab is investigated. The fluid flow outside and inside of the slab is governed by the Stokes and Darcy equations, respectively, and the so-called Beavers and Joseph slip boundary conditions are enforced on the slab surface. The problem is solved for a distant particle with length scale a in terms of the small parameter a/d where d designates the large particle–slab separation. This is achieved by asymptotically inverting a relevant boundary-integral equation on the particle surface, which has been recently proposed for any particle location (distant or close particle) in Khabthani et al. (J Fluid Mech 713:271–306, 2012). It is found that at order O(a/d) the slab behaves for any particle shape as a solid plane no-slip wall while the slab properties (thickness, permeability, associated slip length) solely enter at O((a/d)²). Moreover, for a spherical particle, the numerical results published in Khabthani et al. (J Fluid Mech 713:271–306, 2012) perfectly agree with the present asymptotic analysis.     

Effects of slip on steady Bödewadt flow of a non-Newtonian fluid

The steady flow arising due to the rotation of a non-Newtonian fluid at a larger distance from a stationary disk is extended to the case where the disk surface admits partial slip. The constitutive equation of the non-Newtonian fluid is modeled by that for a Reiner–Rivlin fluid. The momentum equation gives rise to a highly nonlinear boundary value problem. Numerical solution of the governing nonlinear equations are obtained over the entire range of the physical parameters. The effects of slip and non-Newtonian fluid characteristics on the momentum boundary layer are discussed in details. It is observed that slip has prominent effect on the velocity field, whereas a predominant influence of the non-Newtonian parameter is observed on the moment coefficient.     

The Maxey–Riley equation: Existence,uniqueness and regularity of solutions

The Maxey–Riley equation describes the motion of an inertial (i.e., finite-size) spherical particle in an ambient fluid flow. The equation is a second-order, implicit integro-differential equation with a singular kernel, and with a forcing term that blows up at the initial time. Despite the widespread use of the equation in applications, the basic properties of its solutions have remained unexplored. Here we fill this gap by proving local existence and uniqueness of mild solutions. For certain initial velocities between the particle and the fluid, the results extend to strong solutions. We also prove continuous differentiability of the mild and strong solutions with respect to their initial conditions. This justifies the search for coherent structures in inertial flows using the Cauchy–Green strain tensor.     

Effects of slip on steady Bödewadt flow and heat transfer of an electrically conducting non-Newtonian fluid

The steady flow and heat transfer arising due to the rotation of a non-Newtonian fluid at a larger distance from a stationary disk is extended to the case where the disk surface admits partial slip. The constitutive equation of the non-Newtonian fluid is modeled by that for a Reiner–Rivlin fluid. The fluid is subjected to an external uniform magnetic field perpendicular to the plane of the disk. The momentum equation gives rise to a highly nonlinear boundary value problem. Numerical solution of the governing nonlinear equations are obtained over the entire range of the physical parameters. The effects of slip, non-Newtonian fluid characteristics and the magnetic interaction parameter on the momentum boundary layer and thermal boundary layer are discussed in detail and shown graphically. It is observed that slip has prominent effects on the velocity and temperature fields.     

Is Landau Damping Possible in a Shear Fluid Flow?

Asymptotic analysis for small long‐wave perturbations of a given stationary shear flow of an ideal fluid with free boundary as is performed. It is shown that small disturbances of the flow are attracted to periodic solution in the case where the governing equations are hyperbolic on the main shear flow solution. A class of shear flows for which Landau damping is realizable, is described. Analytical results obtained are validated by numerical calculations.     

Solutions to a class of transport problems with radially dominant convection

For fluid systems dealing with drops and bubbles, there are many situations in which the flow is dominated by a radial field. An analysis is carried out for a general class of problems, in which the primary flow is a purely radial type in a spherical geometry and the secondary flow is a perturbation on it. In particular, the flow solutions are obtained for a particle in extensional flow, rotating particle, and a particle in a linear shear flow. In addition, the steady state heat/mass flow equations with radial convection are solved in a fairly general form for spherical boundaries. The solutions lead to a new class of polynomials for the radial functions of the separated solutions. Some of the fundamental properties of these polynomials have also been derived.     

Resonant Flow of a Rotating Fluid Past an Obstacle: The General Case

We consider the flow of a rotating fluid past an antisymmetric obstacle placed on the axis of a cylindrical tube, for the case when the upstream flow is nearly resonant, or critical, so that the speed of a free linear long wave is nearly zero in the frame of reference of the obstacle. The perturbed flow is dominated by the resonant mode, whose amplitude satisfies a forced Korteweg—de Vries equation in this general case when the upstream flow contains radial shear andjor radially dependent angular velocity.     

A sphere in a second degree polynomial creeping flow parallel to a wall

Comprehensive results are provided for the creeping flow arounda spherical particle in a viscous fluid close to a plane wall,when the external velocity is parallel to the wall and variesas a second degree polynomial in the coordinates. By linearityof Stokes equations, the solution is a sum of flows for typicalunperturbed flows: a pure shear flow, a ‘modulated shearflow’, for which the rate of shear varies linearly inthe direction normal to the wall, and a quadratic flow. Solutionsconsidered here use the bipolar coordinates technique. Theycomplement the accurate results of Chaoui and Feuillebois (2003)for the pure shear flow. The solution of Goren and O'Neill (1971)for the quadratic flow is reconsidered and a new analyticalsolution is derived for the ambient modulated shear flow. Theperturbed flow fields for these two cases are presented in detailand discussed. Results for the force and torque friction factorsare provided with a 5 x 10^–17 accuracy as a reference.For the quadratic flow, there is a force and a torque on a fixedsphere. A minimum value of the torque is found for a gap ofabout 0·18a, where a is the sphere radius. This minimumis interpreted in term of the corresponding flow structure.For the modulated shear flow, there is only a torque. The freemotion of a sphere in an ambient quadratic flow is also determined.     

Optical Stokes Flow Estimation: An Imaging-Based Control Approach

We present an approach to particle image velocimetry based on optical flow estimation subject to physical constraints. Admissible flow fields are restricted to vector fields satifying the Stokes equation. The latter equation includes control variables that allow to control the optical flow so as to fit to the apparent velocities of particles in a given image pair. We show that when the real unknown flow observed through image measurements conforms to the physical assumption underlying the Stokes equation, the control variables allow for a physical interpretation in terms of pressure distribution and forces acting on the fluid. Although this physical interpretation is lost if the assumptions do not hold, our approach still allows for reliably estimating more general and highly non-rigid flows from image pairs and is able to outperform cross-correlation based techniques. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Gliding motion of bacteria on power‐law slime

The present investigation deals with an undulating surface model for the motility of bacteria gliding on a layer of non‐Newtonian slime. The slime being the viscoelastic material is considered as a power‐law fluid. A hydrodynamical model of motility involving an undulating cell surface which transmits stresses through a layer of exuded slime to the substratum is examined. The non‐linear differential equation resulting from the balance of momentum and mass is solved numerically by a finite difference method with an iteration technique. The manner in which the various exponent values of the power‐law flow affect the structure of the boundary layer is delineated. A comparison is made of the power‐law fluid with the Newtonian fluid. For the power‐law fluid with respect to different power‐law exponent values, shear‐thinning and shear‐thickening effects can be observed, respectively. Copyright © 2004 John Wiley & Sons, Ltd.     

Evolution Equations for Long,Nonlinear Internal Waves in Stratified Shear Flows

Evolution equations for long, nonlinear internal waves are derived when the basic stratified shear flow has a slow temporal and spatial variation as well as the usual dependence on the vertical coordinate. When the horizontal waveguide has a limited vertical extent the evolution equation is a variable coefficient Korteweg-deVries equation, while in the deep fluid case the evolution equation is a variable coefficient Benjamin-Davis-Ono equation. Explicit expressions are obtained for the coefficients of these equations.