Periodic solution of turbulent oscillating channel flows

The time-dependent turbulent Navier–Stokes equations are solved numerically by a finite element method with an algebraic eddy viscosity model (Baldwin–Lomax formulation) for oscillating turbulent channel flows. The method of averaging is used to analyse the resulting periodic motion of the fluid. Numerical results are obtained for various Strouhal numbers and relative amplitudes. A comparison is made between the numerical and published experimental results. It appears that for low relative amplitudes in a certain range of frequencies the agreement is satisfactory.     

Direct numerical simulation of particle behaviour in the wall region of turbulent flows in horizontal channels

The motion of small particles in the wall region of turbulent channel flows has been investigated using direct numerical simulation. It is assumed that the particle concentration is low enough to allow the use of one-way coupling in the calculations, i.e. the fluid moves the particles but there is no feedback from the particles on the fluid motion. The velocity of the fluid is calculated by using a pseudospectral, direct solution of the Navier-Stokes equations. The calculations indicate that particles tend to segregate into the low-speed regions of the fluid motion near the wall. The segregation tendency depends on the time constant of the particle made non-dimensional with the wall shear velocity and kinematic viscosity. For very small and very large time constants, the particles are distributed more uniformly. For intermediate time constants (of the order 3), the segregation into the low-speed fluid regions is the highest. The finding that segregation occurs for a range of particle time constants is supported by experimental results. The findings regarding the more uniform distributions, however, still remain to be verified against experimental data which is not yet available. For horizontal channel flows, it is also found that particles are resuspended by ejections (of portions of the low-speed streaks) from the wall and are, therefore, primarily associated with low-speed fluid. The smaller particles are flung further upwards and, as they fall back towards the wall, they tend to be accelerated close to the fluid velocity. The larger particles have greater inertia and, consequently, accelerate to lower velocities giving higher relative velocities. This velocity difference, as a function of wall-normal distance, follows the same trend as in experiments but is always somewhat smaller in the calculations. This appears to be due to the Reynolds number for the numerical simulation being smaller than that in the experiment. It is concluded that the average particle velocity depends not only on the wall variables for scaling, but also on outer variables associated with the mean fluid velocity and fluid depth in the channel. This is because fluid depth in combination with the wall shear velocity determines how much time a particle, of a given size and density, spends in the outer flow and, hence, how close it gets to the local fluid velocity.     

Überlegungen zur Scherviskosität von Suspensionen aufgrund einer Dimensionsanalyse

Dimensional analysis of the motion of solid particles suspended in a fluid phase shows that the macroscopic relative shear viscosity of suspensions generally depends not only on the volume concentration and particle shape but also on two Reynolds numbers and a dimensionless sedimentation number. These dimensionless numbers are formed using parameters characterizing the structure and motion of the suspension at the microscopic level. The analysis was based on the assumptions that the dispersed particles are rigid and sufficiently large that Brownian motion may be neglected, that the continuous fluid phase is Newtonian and that the interactions between particles and between particles and fluid phase are only hydrodynamic. The Reynolds numbers describe the influence of the inertial forces at the microscopic level, and the sedimentation number the influence of gravity. The dimensionless numbers can be neglected if their values are much smaller than one. For each of the dimensionless numbers both the shear rate and the particle size influence the shear viscosity. Thus sedimentation number is large for low shear rates, whereas the Reynolds numbers are large for high shear rates. The viscosity function for one suspension can be transformed into the viscosity function for another suspension with geometrically similar particles but of a different size. The scale-up rules are derived from the requirement that the relevant dimensionless numbers must be constant. The influence of non-hydrodynamic effects at the microscopic level on the shear viscosity can be detected by deviations from the derived scale-up rules.     

活性流体流变行为的布朗动力学模拟研究

活性流体在用于开发新材料方面具有巨大潜力, 满足这一需求就要定量掌握活性流体所表现的特殊力学行为, 特别是流变行为. 扩展布朗运动方程, 建立自驱动活性粒子的运动模型, 基于反向非平衡法确定活性流体的黏度, 考察活性粒子体积分数、直行速度和转向扩散系数对活性流体流变行为的影响规律, 确定活性流体特殊流变行为的形成机理. 结果表明, 活性流体的流变曲线可被划分为黏度下降区、过渡区和牛顿区; 活性粒子体积分数越高, 活性流体的非牛顿特性越显著, 活性粒子的直行运动引起活性流体在低剪切速率区域黏度下降, 直行运动和转向运动的耦合作用导致中剪切速率区域流变曲线非单调变化, 活性粒子频繁发生转向运动会导致活性流体非牛顿特性受到抑制; 活性流体的宏观流变学特性和粒子的涨落直接相关, 活性粒子体积分数越高、直行速度越快和转向扩散系数越小, 活性流体中活性粒子越容易产生显著的涨落; 低剪切速率区域内活性粒子涨落明显, 随着剪切速率增大, 活性粒子的涨落逐渐被削弱, 粒子的聚集结构不断被破坏, 最终体系的流变行为类似一般被动流体.      

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.     

Turbulence modulation for gas - particle flow in vertical tube and horizontal channel using four-way Eulerian-Lagrangian approach

Turbulence modulation of gas-solid flow in vertical tube and horizontal channel in dilute and moderately dense suspensions is investigated numerically using a four way Eulerian-Lagrangian approach. Low Reynolds number k-l model is used for analyzing the fluid phase motion. A new model is presented based on a source-term formulation, which can predict fluid phase turbulence augmentation due to the presence of large particles and damping of turbulence due to small particles in the core of the channel and tube. Particle-particle and particle-wall collisions are simulated based on a deterministic approach, and coupling terms representing the fluid-particle interactions are also taken into account. The predicted fluid mean velocity and turbulence intensity profiles are in good agreement with the available experimental data. Additional numerical simulation results for variation of the eddy viscosity, turbulence production and dissipation are also presented for different values of loading ratios.     

Particle collisions in a turbulent flow

An equation for the two-point probability density function of the two-particle the coordinate and velocity distribution is obtained. A closed system of equations for the first and second two-point moments of the velocity fluctuations of a pair of particles with allowance for the turbulent flow inhomogeneity is given. Boundary conditions for the equations of the particle concentration and the intensity of the relative random velocity during particle collision are obtained. A unified formula describing the interparticle collision process as a result of turbulent motion and the average relative particle velocity slip is obtained for the kernel of the coagulation equation. The effect of the average velocity slip of the particles and the carrier phase on the parameters of motion of the dispersed admixture and its coagulation is investigated on the basis of a two-point two-time velocity fluctuation autocorrelation function with two time and space scales representing the energy-bearing and small-scale motion of the fluid phase.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 104–116, March–April, 1996.     

Instabilities during batch sedimentation in geometries containing obstacles: A numerical and experimental study

Batch sedimentation of non‐colloidal particle suspensions is studied with nuclear magnetic resonance flow visualization and continuum‐level numerical modelling of particle migration. The experimental method gives particle volume fraction as a function of time and position, which then provides validation data for the numerical model. A finite element method is used to discretize the equations of motion, including an evolution equation for the particle volume fraction and a generalized Newtonian viscosity dependent on local particle concentration. The diffusive‐flux equation is based on the Phillips model (Phys. Fluids A 1992; 4 :30–40) and includes sedimentation terms described by Zhang and Acrivos (Int. J. Multiphase Flow 1994; 20 :579–591). The model and experiments are utilized in three distinct geometries with particles that are heavier and lighter than the suspending fluid, depending on the experiment: (1) sedimentation in a cylinder with a contraction; (2) particle flotation in a horizontal cylinder with a horizontal rod; and (3) flotation around a rectangular inclusion. Secondary flows appear in both the experiments and the simulations when a region of higher density fluid is above a lower density fluid. The secondary flows result in particle inhomogeneities, Rayleigh–Taylor‐like instabilities, and remixing, though the effect in the simulations is more pronounced than in the experiments. Published in 2007 by John Wiley & Sons, Ltd.     

A model of the steady-state motion of suspensions

The authors examine the steady-state one-dimensional motions of suspensions whose particles have a density equal to that of the corresponding dispersion medium. As a whole, the mechanical behavior of such suspensions is described by equations of motion that coincide in form with the Navier-Stokes equations for a certain incompressible fluid whose viscosity is a known function of the particle concentration in the suspensions. To close these equations, the authors postulate a principle of minimum energy dissipation for steady-state motion, which plays the paxt of an equation of state for the suspension. This new equation permits the determination of the spatial distribution in the concentration of solids. Exact solutions are presented for certain variational problems associated with the Poiseuille flow of a fluid of this kind in circular tubes and Couette flows between concentric cylinders and parallel planes. It is shown that in most cases separation of the suspension takes place.     

The dynamic equations of motion for a two-component system

The motion of solid particles in a fluid flow is represented as a random process with independent increments. The resulting kinetic equation for the particle distribution has the form previously proposed [1]. The solution to this equation provides a system of equations for the hydrodynamics of the assembly of solid particles. These equations differ from ones previously proposed [2, 3] in having additional terms related to relative motion of the components, whose presence is due to anisotropy in the distribution of the normal stresses in the pseudogas.I am indebted to V. G. Levich for valuable discussions and for constant interest in the work.     

Motion of deformable drops in pipes and channels using Navier–Stokes equations

The motion of deformable drops in pipes and channels is studied using a level set approach in order to capture the interface of two fluids. The interface is described as the zero level set of a smooth function, which is defined to be the signed normal distance from the interface. In order to solve the Navier–Stokes equations, a second‐order projection method is used. The dimensionless parameters of the problem are the relative size of the drop to the size of the pipe or channel cross‐section, the ratio of the drop viscosity to the viscosity of the suspending fluid and the relative magnitude of viscous forces to the surface tension forces. The shape of the drop, the velocity field and the additional pressure loss due to the presence of the drop, varying systematically with the above‐mentioned dimensionless parameters, are computed. Copyright © 2000 John Wiley & Sons, Ltd.     

Slow motion of a slip spherical particle perpendicular to two plane walls

The problem of the quasisteady motion of a spherical fluid or solid particle with a slip-flow surface in a viscous fluid perpendicular to two parallel plane walls at an arbitrary position between them is investigated theoretically in the limit of small Reynolds number. To solve the axisymmetric Stokes equation for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both circular cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the plane walls by the Hankel transform and then on the particle surface 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 relative viscosity or slip coefficient of the particle and of the relative separation distances between the particle and the confining walls. For the motions of a spherical particle normal to a single plane wall and of a no-slip sphere perpendicular to two plane walls, our drag results are in good agreement with the available solutions in the literature for all relative particle-to-wall spacings. The boundary-corrected 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 geometry, but there are exceptions. For a specified wall-to-wall spacing, the drag force is minimal when the particle is situated midway between the two plane walls and increases monotonically when it approaches either of the walls. The boundary effect on the particle motion normal to two plane walls is found to be significant and much stronger than that parallel to them.     

Stability of plane poiseuille flow of a viscous anisotropic fluid

Experimental investigations show that the presence in a fluid of fibers and rigid asymmetric particles leads to a greater stability of flow in tubes and lowers the turbulent frictional resistance in a certain range of Reynolds numbers [1]. In the present paper, the anisotropic structure of a fluid with additives is described by Ericksen's rheological model [2]. The parameters of the model are particularized in accordance with the paper [3] of Pilipenko, Kalinichenko, and Lemak, and in the limiting case of weak Brownian motion allowance is made for the effect of the predominant orientation of the particles and the influence of additives on the longitudinal and shear viscosity. The stability of the Poiseuille flow is considered in the linear formulation. In an anisotropic viscous fluid, an equation of Orr-Sommerfeld type has a singular point. A rule for choosing the path of integration avoiding the singular point is obtained on the basis of a generalization of the method of Dikii [4] proposed in an investigation of the stability of the flow of an ideal fluid. The results of numerical calculations of the neutral stability curve for two-dimensional perturbations are given.     

Single-Fluid Model of a Mixture for Laminar Flows of Highly Concentrated Suspensions

A model of laminar flow of a highly concentrated suspension is proposed. The model includes the equation of motion for the mixture as a whole and the transport equation for the particle concentration, taking into account a phase slip velocity. The suspension is treated as a Newtonian fluid with an effective viscosity depending on the local particle concentration. The pressure of the solid phase induced by particle-particle interactions and the hydrodynamic drag force with account of the hindering effect are described using empirical formulas. The partial-slip boundary condition for the mixture velocity on the wall models the formation of a slip layer near the wall. The model is validated against experimental data for rotational Couette flow, a plane-channel flow with neutrally buoyant particles, and a fully developed flow with heavy particles in a horizontal pipe. Based on the comparison with the experimental data, it is shown that the model predicts well the dependence of the pressure difference on the mixture velocity and satisfactorily describes the dependence of the delivered particle concentration on the flow velocity.     

The slow motion of two touching fluid spheres along their line of centers

The creeping motion along their line of centers of two fluid spheres in contact is analyzed. An exact solution is presented. Corrections to the Hadamard—Rybezynski equation are tabulated for various particle radii ratios and particle fluid to external fluid viscosity ratios. In the limit of infinite particle viscosity, these corrections are shown to agree with previous calculations for rigid spheres.     

A boundary element method and spectral analysis model for small‐amplitude viscous fluid sloshing in couple with structural vibrations

A boundary element method (BEM) is presented for the coupled motion analysis of structural vibrations with small‐amplitude fluid sloshing in two‐dimensional space. The linearized Navier–Stokes equations are considered in the frequency domain and transformed into a Laplace equation and a Helmholtz equation with pure imaginary constant. An appropriate fundamental solution for the Helmholtz equation is provided. The conditions of zero stress are imposed on the free surface, and non‐slip conditions of fluid particles are imposed on the walls of the container. For rigid motion models, the expressions for added mass and added damping to the structural motion equations are obtained. Numerical examples are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Generalized Plain Couette Flow and Heat Transfer in a Composite Channel

An analytical study of fluid flow and heat transfer in a composite channel is presented. The channel walls are maintained at different constant temperatures in such a way that the temperatures do not allow for free convection. The upper plate is considered to be moving and the lower plate is fixed. The flow is modeled using Darcy–Lapwood–Brinkman equation. The viscous and Darcy dissipation terms are included in the energy equation. By applying suitable matching and boundary conditions, an exact solution has been obtained for the velocity and temperature distributions in the two regions of the composite channel. The effects of various parameters such as the porous medium parameter, viscosity ratio, height ratio, conductivity ratio, Eckert number, and Prandtl number on the velocity and temperature fields are presented graphically and discussed.     

相对运动中一类质点运动方程的同伦摄动解

相对运动中质点的运动方程大多数是非线性的,其精确显式求解是一个难点。本文基于同伦摄动分析方法,研究了相对运动中一类质点运动非线性微分方程的近似显式解。先构造了系统的同伦方程,再结合Lindstedt–Poincare方法和系统的初始条件,推导了系统自由振动的固有频率,求解了系统的位移近似响应,并通过数值仿真验证了理论分析解的正确性,为相对运动的质点运动方程求解提供了一种新的求解方法。