NUMERICAL PREDICTION OF FORCE ON RECTANGULAR CYLINDERS IN OSCILLATING VISCOUS FLOW

Numerical results are presented for an oscillating viscous flow past a square cylinder with square and rounded corners and a diamond cylinder with square corners at Keulegan–Carpenter numbers up to 5. This unsteady flow problem is formulated by the two-dimensional Navier–Stokes equations in vorticity and stream-function form on body-fitted coordinates and solved by a finite-difference method. Second-order Adams-Bashforth and central-difference schemes are used to discretize the vorticity transport equation while a third-order upwinding scheme is incorporated to represent the nonlinear convective terms. Since the vorticity distribution has a mathematical singularity at a sharp corner and since the force coefficients are found in experiments to be sensitive to the corner radius of rectangular cylinders, a grid-generation technique is applied to provide an efficient mesh system for this complex flow. Local grid concentration near the sharp corners, instead of any artificial treatment of the sharp corners being introduced, is used in order to obtain high numerical resolution. The elliptic partial differential equation for stream function and vorticity in the transformed plane is solved by a multigrid iteration method. For an oscillating flow past a rectangular cylinder, vortex detachment occurs at irregular high frequency modes at KC numbers larger than 3 for a square cylinder, larger than 1 for a diamond cylinder and larger than 3 for a square cylinder with rounded corners. The calculated drag and inertia coefficients are in very good agreement with the experimental data. The calculated vortex patterns are used to explain some of the force coefficient behavior.     

High-order resolution Eulerian–Lagrangian simulations of particle dispersion in the accelerated flow behind a moving shock

This paper presents a computational study of the two-dimensional particle-laden flow developments of bronze particle clouds in the accelerated flow behind a moving normal shock. Particle clouds with a particle volume concentration of 4% are arranged initially in a rectangular, triangular and circular shape. Simulations are performed with a recently developed high-order resolution Eulerian–Lagrangian method that approximates the Euler equations governing the gas dynamics with the improved high order weighted essentially non-oscillatory (WENO-Z) scheme, while individual particles are traced in the Lagrangian frame using high-order time integration schemes. Reflected shocks form ahead of all the cloud shapes. The detached shock in front of the triangular cloud is weakest. At later times, the wake behind the cloud becomes unstable, and a two-dimensional vortex-dominated wake forms. Separated shear layers at the edges of the clouds pull particles initially out of the clouds that are consequently transported along the shear layers. Since flows separated trivially at sharp corners, particles are mostly transported out of the cloud into the flow at the sharp front corner of the rectangular cloud and at the trailing corner of the triangular cloud. Particles are transported smoothly out of the circular cloud, since it lacks sharp corners. At late times, the accelerated flow behind the running shock disperses the particles in cross-stream direction the most for the circular cloud, followed by the rectangular cloud and the triangular cloud.     

垂直旋转流场中单颗粒运动状态判别及受力分析

为研究单颗粒在旋转流场中的运动状态及受力情况,以毫米级球形颗粒为例,利用旋转流场颗粒运动装置,通过使用摄像机记录颗粒在流场中的运动轨迹以获取其运动参数,分析了不同转速和颗粒直径条件下颗粒的运动轨迹,拟合得到了颗粒运动状态判别公式以及颗粒运动轨迹公式,分析了颗粒在旋转流场中的受力情况。结果表明,颗粒在旋转流场平衡状态下运动状态主要分为两类,一类是未离开壁面保持静止,另一类是离开壁面保持稳定周向运动;颗粒进行周向运动的轨迹为椭圆形,并且圆心随着转速的增大靠近旋转中心,而随着粒径的增大靠近壁面;颗粒在旋转流场的运动过程中主要受到离心力和旋转科式力作用。     

A numerical study of flows driven by a rotating magnetic field in a square container

The magnetically induced fluid flow in a square container is investigated by means of numerical simulations. Low frequency/ low induction conditions are assumed. The effect of the rotating magnetic field gives rise to a time-independent magnetic body force, computed via the electrical potential equation and Ohm's law and a time-dependent part that is neglected due to the low-interaction parameter. The magnetic body force calculation is verified successfully by comparison with the exact solution. The behavior of the fluid flow in the square container reveals similar features to the flow in the cylindrical container, for instance, in the dependence on the intensities of the magnetic field. However, we did find differences in the velocity field distribution. Particularly, in the finite as well as infinite geometry, the velocity field is influenced by the corner of the container and remains non-axisymmetric in a wide range of Taylor numbers.     

Turbulent transport of particles in a straight square duct

Turbulent flow through a duct of square cross-section gives rise to off-axis secondary flows, which are known to transfer momentum between fluid layers thereby flattening the velocity profile. The aim of this study is to investigate the role of the secondary flows in the transport and dispersion of particles suspended in a turbulent square duct flow. We have numerically simulated a flow through a square duct having a Reynolds number of Re_τ = 300 through discretization of the Navier–Stokes equations, and followed the trajectories of a large number of passive tracers and finite-inertia particles under a one-way coupling assumption. Snapshots of particle locations and statistics of single-particle and particle pair dispersion were analyzed. It was found that lateral mixing is enhanced for passive tracers and low-inertia particles due to the lateral advective transport that is absent in straight pipe and channels flows. Higher inertia particles accumulate close to the wall, and thus tend to mix more efficiently in the streamwise direction since a large number of the particles spend more time in a region where the mean fluid velocity is small compared to the bulk. Passive tracers tend to remain within the secondary swirling flows, circulating between the core and boundary of the duct.     

Three-dimensional, three-component wall-PIV

This paper describes a new time-resolved three-dimensional, three-component (3D-3C) measurement technique called wall-PIV. It was developed to assess near wall flow fields and shear rates near non-planar surfaces. The method is based on light absorption according to Beer–Lambert’s law. The fluid containing a molecular dye and seeded with buoyant particles is illuminated by a monochromatic, diffuse light. Due to the dye, the depth of view is limited to the near wall layer. The three-dimensional particle positions can be reconstructed by the intensities of the particle’s projection on an image sensor. The flow estimation is performed by a new algorithm, based on learned particle trajectories. Possible sources of measurement errors related to the wall-PIV technique are analyzed. The accuracy analysis was based on single particle experiments and a three-dimensional artificial data set simulating a rotating sphere.     

Effects of vortex pairing on particle dispersion in turbulent shear flows

Particle dispersion in large-scale dominated turbulent shear flow is investigated numerically with special emphasis on the effects of the vortex-pairing phenomenon. The particle dispersion is visualized numerically by following the particle trajectories in a flow consisting of large vortices which are undergoing pairing interaction. The flow field is generated by a discrete vortex method. Important global and local fiow quantities from the numerical simulation compare reasonably well with experimental measurements.

For both cases of point sources with continuous particle release and an initially distributed line source, the particle dispersion results demonstrate that the extent of particle dispersion depends strongly on the Stokes number, the ratio of the particle aerodynamic response time to the characteristic time of the vortex-pairing flow field. Particles with relatively small Stokes numbers disperse laterally at approximately the saine rate as that of the fluid particles and particles with large Stokes numbers disperse much less than the fluid particles. Particles with intermediate Stokes numbers (0.5-5) may be dispersed laterally farther than the fiuid particles and may actually be flung out of the vortex structures. Due to the strong particie entrainment power, the flow during the vortex-pairing process seems to produce higher particle lateral dispersion than the pre-pairing and post-pairing flows.     

Turbulence-induced preferential concentration of solid particles in microgravity conditions

A box of near-isotropic, particle-laden turbulence was flown aboard NASA's reduced gravity aircraft in order to measure the turbulence-induced preferential concentration of solid particles in microgravity. Three particle sets of Stokes numbers based on the fluid Kolmogorov time scale of approximately 0.5, 5, and 50 were tested for relative amounts of preferential concentration. Eight fans in each corner of a Lexan box generated fluid turbulence. Particle concentrations were measured using an imaging system consisting of a camera viewing perpendicular to a white light sheet. Post-processing of video images found largest concentration gradients for the intermediate-sized particles of Stokes number 5, closely followed by the Stokes number 0.5 particles. The experimental results agreed well with the trends seen in direct numerical simulations. The quantitative effects of turbulence modulation by the presence of particles were not measured in the experiment, but were most likely present. Received: 10 February 2000 / Accepted: 9 November 2001     

Mathematical modelling of two‐phase non‐Newtonian flow in a helical pipe

Governing equations for a two‐phase 3D helical pipe flow of a non‐Newtonian fluid with large particles are derived in an orthogonal helical coordinate system. The Lagrangian approach is utilized to model solid particle trajectories. The interaction between solid particles and the fluid that carries them is accounted for by a source term in the momentum equation for the fluid. The force‐coupling method (FCM), developed by M.R. Maxey and his group, is adopted; in this method the momentum source term is no longer a Dirac delta function but is spread on a numerical mesh by using a finite‐sized envelop with a spherical Gaussian distribution. The influence of inter‐particle and particle–wall collisions is also taken into account. Copyright © 2005 John Wiley & Sons, Ltd.     

On model equations for particle dispersion in inhomogeneous turbulence

Comparisons are made between the Advection–Diffusion Equation (ADE) approach for particle transport and the two-fluid model approach based on the PDF method. In principle, the ADE approach offers a much simpler way of calculating the inertial deposition of particles in a turbulent boundary layer than that based on the PDF approach. However the ADE equations that have recently been used are only strictly valid for a simple Gaussian process when particle inertia is small. Using a prescribed, but in general non-Gaussian random particle velocity field, it is shown that the net particle mass flux contains a drift term in addition to that from the mean velocity of the particle velocity field, associated with the compressibility of the velocity field. Furthermore the diffusive flux in general depends not only upon the gradient of the mean concentration (true only for a Gaussian random flow field) but also upon higher order derivatives whose relative contribution depends on diffusion coefficients D_ijk… etc. These coefficients depend upon the statistical moments associated with random displacements and compressibility of the particle flow field along particle trajectories which in turn depend upon particle inertia. In contrast the PDF approach offers the advantage of using a simple gradient (Gaussian) approximation in particle phase space which can lead to a non-Gaussian spatial dispersion process when particle inertia is important. Conditions based on the particle mean free path are derived for which a simple ADE is appropriate. Some of the features of particle transport in an inhomogeneous turbulent flow are illustrated by examining particle dispersion in a random flow field composed of pairs of counter rotating vortices which has an rms velocity which increase linearly from a stagnation point.     

Numerical simulation of the accumulation of heavy particles in a circular bounded vortex flow

In present work, an Eulerian–Lagrangian CFD model based on the discrete element method (DEM) and immersed boundary method (IBM) has been developed, validated and used to investigate the accumulation of heavy particles in a circular bounded viscous vortex flow. The inter-particle and particle-wall collisions are resolved by a hard-sphere model. Effects of one-way and two-way coupling, Reynolds number, and particle diameter are systematically explored. Results show that, in case of one-way coupling, the majority of particles will spiral into an accumulation point located near the stagnation point of the flow field. The accumulation point represents a stable equilibrium point as the drag created by the flow field balances the destabilizing centrifugal force on the particle. However, in case of two-way coupling, there does not exist a stable accumulation point due to the strong interaction between the particles and fluid dynamics. Instead most particles are expelled from the circular domain and accumulate on the confining wall. The percentage of accumulated particles on the wall increases with increasing Reynolds number and particle diameter. Moreover, influence of three well-known drag models is also studied and they give consistent results on the particle accumulation behavior, although small quantitative differences can still be discerned.     

An extended finite element method for the simulation of particulate viscoelastic flows

We present an extended finite element method (XFEM) for the direct numerical simulation of the flow of viscoelastic fluids with suspended particles. For moving particle problems, we devise a temporary arbitrary Lagrangian–Eulerian (ALE) scheme which defines the mapping of field variables at previous time levels onto the computational mesh at the current time level. In this method, a regular mesh is used for the whole computational domain including both fluid and particles. A temporary ALE mesh is constructed separately and the computational mesh is kept unchanged throughout the whole computations. Particles are moving on a fixed Eulerian mesh without any need of re-meshing. For mesh refinements around the interface, we combine XFEM with the grid deformation method, in which nodal points are redistributed close to the interface while preserving the mesh topology. Our method is verified by comparing with the results of boundary fitted mesh problems combined with the conventional ALE scheme. The proposed method shows similar accuracy compared with boundary fitted mesh problems and superior accuracy compared with the fictitious domain method. If the grid deformation method is combined with XFEM, the required computational time is reduced significantly compared to uniform mesh refinements, while providing mesh convergent solutions. We apply the proposed method to the particle migration in rotating Couette flow of a Giesekus fluid. We investigate the effect of initial particle positions, the Weissenberg number, the mobility parameter of the Giesekus model and the particle size on the particle migration. We also show two-particle interactions in confined shear flow of a viscoelastic fluid. We find three different regimes of particle motions according to initial separations of particles.     

Numerical study of an inviscid incompressible flow through a channel of finite length

A two‐dimensional inviscid incompressible flow in a rectilinear channel of finite length is studied numerically. Both the normal velocity and the vorticity are given at the inlet, and only the normal velocity is specified at the outlet. The flow is described in terms of the stream function and vorticity. To solve the unsteady problem numerically, we propose a version of the vortex particle method. The vorticity field is approximated using its values at a set of fluid particles. A pseudo‐symplectic integrator is employed to solve the system of ordinary differential equations governing the motion of fluid particles. The stream function is computed using the Galerkin method. Unsteady flows developing from an initial perturbation in the form of an elliptical patch of vorticity are calculated for various values of the volume flux of fluid through the channel. It is shown that if the flux of fluid is large, the initial vortex patch is washed out of the channel, and when the flux is reduced, the initial perturbation evolves to a steady flow with stagnation regions. Copyright © 2008 John Wiley & Sons, Ltd.     

Self-consistent particle simulation of model-stabilized colloidal suspensions

Dynamics of model-stabilized colloidal suspensions were investigated by the self-consistent particle simulation method (SC), a new simulation algorithm that takes into account the interaction between the particles and suspending fluid. In this method, the fluid-particle interaction is introduced self-consistently by combining the finite element method (FEM) for fluid motion with Brownian dynamics (BD) for particle dynamics. To validate the reliability of the proposed algorithm, the shear dynamics of the stable particle suspensions were investigated. Relative viscosity and microstructure as a function of dimensionless shear rate at different volume fractions were in good agreement with previous observations. The robustness of the method was also verified through numerical convergence test. The effect of the fluid-particle interaction was well represented in simulations of two model problems, pressure-driven channel flow and rotating Couette flow. Plug-shaped velocity profile was observed in pressure-driven channel flow, which arised from shear thinning behavior of the stable suspension. In rotating Couette flow, shear banded nonlinear flow profile was observed. Although full hydrodynamic interaction (HI) was not rigorously taken into account, it successfully captured the macroscopic structure-induced flow field. It also takes advantage of the geometrical adaptability of FEM and computational efficiency of BD. We expect this newly developed simulation platform to be useful and efficient for probing the complex flow dynamics of particle systems as well as for practical applications in the complex flow of complex fluids.     

A model for tracking inertial particles in a lattice Boltzmann turbulent flow simulation

A computationally inexpensive model for tracking inertial particles through a turbulent flow is presented and applied to the turbulent flow through a square duct having a friction Reynolds number of Re_τ = 300. Prior to introducing particles into the model, the flow is simulated using a lattice Boltzmann computation, which is allowed to evolve until a steady state turbulent flow is achieved. A snapshot of the flow is then stored, and the trajectories of particles are computed through the flow domain under the influence of this static probability field. Although the flow is not computationally evolving during the particle tracking simulation, the local velocity is obtained stochastically from the local probability function, thus allowing the dynamics of the turbulent flow to be resolved from the point of view of the suspended particles. Particle inertia is modeled by using a relaxation parameter based on the particle Stokes number that allows for a particle velocity history to be incorporated during each time step. Wall deposition rates and deposition patterns are obtained and exhibit a high level of agreement with previously obtained DNS computational results and experimental results for a wide range of particle inertia. These results suggest that accurate particle tracking through complex turbulent flows may be feasible given a suitable probability field, such as one obtained from a lattice Boltzmann simulation. This in turn presents a new paradigm for the rapid acquisition of particle transport statistics without the need for concurrent computations of fluid flow evolution.     

Effect of weak rotation on natural convection in a horizontal porous cylinder

Natural convection in a fluid saturated porous medium confined in a horizontal circular cylinder and rotating about its axis, with isothermal boundary conditions and uniform internal heat sink, is studied by both numerical and perturbation methods. No symmetry with respect to the vertical diameter is expected for the flow and temperature fields and the whole region must be involved in the computation. Only the weak rotation regime, for which the centrifugal force is negligible compared to gravity, is considered. Governing equations for the two-dimensional flow field are solved in both rotating and non-rotating coordinate systems. Results indicate that rotation significantly decreases the radial amplitude of fluid particle trajectories in the radial direction and thus reduces the overall heat transfer.     

Flows near stagnation points in non-orthogonally colliding disperse viscous flows

Within the framework of a two-fluid dusty-gas model, a family of steady self-similar flows near a stagnation point formed by two different non-orthogonally colliding viscous incompressible streams, one of which contains solid inertial particles, is studied. The limiting case of non-orthogonal impingement of a viscous disperse medium on a rigid wall is considered separately. The possibility of the formation of multiple particle accumulation zones on the envelopes of particle trajectories and on the contact (in the limiting case, rigid) surface is demonstrated. The local structure of the particle velocity and concentration fields is studied. The threshold values of the governing parameters corresponding to qualitative flow pattern reconstruction are found.     

荷电颗粒在旋转环形通道内分离性能研究

为研究带电旋转环形通道内荷电颗粒的运动和沉积特性,本文使用计算流体力学两相流离散颗粒法对带电旋转环形通道内的荷电颗粒的运动过程进行了模拟。根据模拟结果分析了不同粒径、电压、入口雷诺数和通道长径比等参数对荷电颗粒运动和沉积的影响,研究了荷电颗粒在旋转通道内离心力与电场力之间的竞争关系,探索了离心力和电场力导致的荷电颗粒运动和沉积变化的规律。结果表明,单个不同粒径颗粒具有不同的颗粒逃逸电压区间,区间的大小随着颗粒粒径的增大而增大,且区间的宽度随着通道长径比的增大将会明显变小;多个颗粒的逃逸率曲线,不同粒径的颗粒将会有不同程度的交叉,随着长径比的增大,颗粒逃逸率曲线的高度与交叉会有明显的减小,而随着转速的增大,颗粒逃逸率曲线的交叉会有一定程度的减小,且高度不会有明显变化。     

On the Particle Paths and the Stagnation Points in Small-Amplitude Deep-Water Waves

In order to obtain quite precise information about the shape of the particle paths below small-amplitude gravity waves travelling on irrotational deep water, analytic solutions of the nonlinear differential equation system describing the particle motion are provided. All these solutions are not closed curves. Some particle trajectories are peakon-like, others can be expressed with the aid of the Jacobi elliptic functions or with the aid of the hyperelliptic functions. Remarks on the stagnation points of the small-amplitude irrotational deep-water waves are also made.