The flow field induced by the torsional oscillations of a spherical cell containing a fluid drop

The flow field induced by the torsional oscillations of a spherical cell containing a fluid drop is examined. It has been found that in addition to the oscillating motion of each fluid particle there exist three standing vortices in every quadrant between the drop and the container. The steady streaming into the drop can be directed either clockwise or counterclockwise, depending on the values of the parameter of the fluids inside and outside the drop. Typical flow fields are shown graphically.     

Investigation of trajectories of inviscid fluid particles in two-dimensional rotating boxes

The particle trajectories of inviscid fluid flow within two-dimensional rotating (elliptic, triangular, and square) boxes are numerically investigated. The source panel method is employed to represent the instantaneous potential interior flow field, and the Runge–Kutta method is used to track the fluid particles. The analytic solutions for the fluid trajectories for the elliptic box are used to verify the numerical accuracy of the method. The numerical error can be reduced to the level of the round-off error if the panels are properly configured and an appropriate number of panels is used. The stagnation of the particles at the corners of the triangular box is successfully predicted with this method. The corner of the square box is found to be a singularity. A logarithmic complex potential is proposed to account for the singularity, using which the stagnation of the particles at the corner in the square box is also captured. The natural frequency of the particles in the rotating elliptic box is constant throughout the flow domain, and the fluid trajectories are epitrochoidal curves. In the triangular box and the square box, the natural frequency strongly depends on the particle position, and the particle trajectories are similar to epitrochoidal curves. In general, the trajectory patterns depend only on the box rotating frequency and the natural frequency of the fluid particle motion.      

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.     

Motion of a sphere down an inclined plane in a viscous flow

The motion of a rigid particle near a wall in a fluid flow is an important element of particle transport by fluids. The aim of this study was to carry out an experimental and theoretical investigation of the gravity-induced motion of a rigid sphere in a viscous fluid in the presence of a transverse flow. The experimental study of this configuration is a way of understanding the specific features of the hydrodynamically constrained particle motion. It is established that the transverse motion of the fluid substantially increases the particle settling velocity, which grows with increase in the transverse flow velocity. This effect is most pronounced for small angles of inclination of the plane. The difference in the particle settling velocities in the presence and absence of the transverse flow could reach a factor of two.     

Calculation of the particle concentration field in the problem of aerosol aspiration from moving air

A mathematical model of a slot sampler is proposed and the aspiration coefficient and the particle concentration field are studied numerically for the problem of aerosol aspiration from a moving gas. In the absence of a particle effect on the gas, the carrier flow is calculated in the viscous incompressible fluid approximation using the FLUENT program package for solving the Navier-Stokes equations. The particle motion equations are supplemented by the equations for determining the particle concentration along the particle trajectories. The particle concentration distribution in the neighborhood of and inside the slot sampler is investigated. Parametric studies of the aspiration coefficient as a function of the ratio of the wind velocity and the aspiration rate are carried out for different Stokes numbers.     

A comparison between a collocation and weak implementation of the rigid‐body motion constraint on a particle surface

In this work, we implemented and compared two different methods to impose the rigid‐body motion constraint on a solid particle moving inside a fluid. We consider a fictitious domain method to easily manage the particle motion. As the solid as well as the fluid inertia are neglected, the particle can be discretized through its boundary only. The rigid‐body motion is imposed via Lagrange multipliers on the boundary. In the first method, such constraints are imposed in discrete points on the boundary (collocation), whereas in the second the constraint is imposed in a weak way on elements dividing the particle surface. Two test problems, that is, a spherical and an ellipsoidal particle in a sheared Newtonian fluid, are chosen to compare the methods. In both cases, the analysis is carried out in 2D as well as in 3D. The results show that for the collocation method an optimal number of collocation points exist leading to the smallest error. However, small variations in the optimal value can generate large deviations. In the weak implementation, the error is only mildly affected by the number of elements used to discretize the particle boundary and by the Lagrange multiplier's interpolation space. A further analysis is carried out to study the effect of an approximated integration of weak constraints. A comparison between the two methods showed that the same accuracy can be achieved by using less constraints if the weak discretization is used. Finally, the rigid‐body motion imposed via weak constraints leads to better conditioned linear systems. Copyright © 2009 John Wiley & Sons, Ltd.     

Suspension of solid particles in a newtonian fluid

We consider the motion of a Newtonian fluid containing solid particles in the case of high concentration (the partial volume of solid is comparable to that of fluid) by using a homogeniztion technique associated with the small parameter ? (the ratio of the particle length to the characteristic macroscopical length). The limit behaviour as ? → 0 is that of fluid with anisotropy properties associated with a microstructure. The evolution equations for the limit flow and for the microstructure are given.     

Slow Motion of an Assemblage of Porous Spherical Shells Relative to a Fluid

The body-force-driven motion of a homogeneous distribution of spherically symmetric porous shells in an incompressible Newtonian fluid and the fluid flow through a bed of these shell particles are investigated analytically. The effect of the hydrodynamic interaction among the porous shell particles is taken into account by employing a cell-model representation. In the limit of small Reynolds number, the Stokes and Brinkman equations are solved for the flow field around a single particle in a unit cell, and the drag force acting on the particle by the fluid is obtained in closed forms. For a suspension of porous spherical shells, the mobility of the particles decreases or the hydrodynamic interaction among the particles increases monotonically with a decrease in the permeability of the porous shells. The effect of particle interactions on the creeping motion of porous spherical shells relative to a fluid can be quite significant in some situations. In the limiting cases, the analytical solution describing the drag force or mobility for a suspension of porous spherical shells reduces to those for suspensions of impermeable solid spheres and of porous spheres. The particle-interaction behavior for a suspension of porous spherical shells with a relatively low permeability may be approximated by that of permeable spheres when the porous shells are sufficiently thick.     

Motion of a spherical solid particle in a nonuniform flow of a viscous incompressible liquid

The effect of a particle on the basic flow is studied, and the equations of motion of the particle are formulated. The problem is solved in the Stokes approximation with an accuracy up to the cube of the ratio of the radius of the sphere to the distance from the center of the sphere to peculiarities in the basic flow. An analogous problem concerning the motion of a sphere in a nonuniform flow of an ideal liquid has been discussed in [1]. We note that the solution is known in the case of flow around two spheres by a uniform flow of a viscous incompressible liquid [2], and we also note the papers [3, 4] on the motion of a small particle in a cylindrical tube.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 71–74, July–August, 1976.     

血液入口段两相流动的层流NS方程数值模拟

本文采用球形刚性颗粒悬浮于牛顿流体的血液模型，基于有相间滑移的两相流动层流Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，采用相间滑移算法（Ｉｎｔｅｒ－ＰｈａｓｅＳｌｉｐＡｌｇｏｒｉｔｈｍ）对圆管内定常轴对称血液入口两相流动问题进行了计算，获得了与实验结果吻合的计算结果，较好地模拟血液流动中红细胞的径向迁移现象，结果表明，采用二相流动模型研究血液流动是一种有前途的方法。     

Full-field measurements of flow through a scaled metal foam replica

Open-celled foam geometries show great promise in heat/mass transfer, chemical treatment, and enhanced mixing applications. Flow measurements on these geometries have consisted primarily of observations of the upstream and downstream effects the foam has on the velocity field. Unfortunately, these observations give little insight into the flow inside the foam. We have performed quantitative flow measurements inside a scaled replica of a metal foam, ϕ = 0.921, D _Cell = 2.5 mm, by Magnetic Resonance Velocimetry to better understand the fluid motion inside the foam and give an alternative method to determine the foam cell and pore sizes. Through these 3-D, spatially resolved measurements of the flow field for a cell Reynolds number of 840, we have shown that the transverse motion of the fluid has velocities 20–30% of the superficial velocity passing through the foam. This strong transverse motion creates and dissipates streamwise jets with peak velocities 2–3 times the superficial velocity and whose coherence length is strongly correlated to the cell size of the foam. This complex fluid motion is described as “mechanical mixing” and is attributed to the geometry of the foam. A mechanical dispersion coefficient, D _M, was formed which demonstrates the transverse dispersion of this geometry to be 14 times the kinematic viscosity and 10 times the thermal diffusivity of air at 20°C and 1 atm.     

Dynamics of particles floating on liquid flowing in a semicircular open channel

The dynamic characteristics of surface-floating particles in liquids flowing in a two-dimensional, semicircular open channel is studied experimentally. For high visibility in the experiments, relatively large particles are employed whose particle-liquid density ratio is either equal to or less than unity. Particles of different size and geometry are tested in a water-glycerin mixture. A video camera traces the pathline of each particle from which the velocity and direction of particle motion are evaluated. Liquid velocity distribution is determined by hot-film anemometry. A modified dynamics (Basset-Boussinesq-Oseen) equation is derived and numerically solved by means of a finite-difference technique to determine fluid velocity. A new dimensionless parameter is disclosed which is pertinent to both particle geometry and fluid flow conditions. It correlates particle trajectory and velocity, trajectory dispersion and fluid-particle velocity ratio.Visiting Scholar on leave from Department of Mechanical Engineering, Fukuyama University, Fukujama, Japan     

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.     

Motion and heat and mass transfer in a disperse admixture in turbulent nonisothermal jets of a gas and a low-temperature plasma

The motion and heat and mass transfer of particles of a disperse admixture in nonisothermal jets of a gas and a low-temperature plasma are simulated with allowance for the migration mechanism of particle motion actuated by the turbophoresis force and the influence of turbulent fluctuations of the jet flow velocity on heat and mass transfer of the particle. The temperature distribution inside the particle at each time step is found by solving the equation of unsteady heat conduction. The laws of scattering of the admixture and the laws of melting and evaporation of an individual particle are studied, depending on the injection velocity and on the method of particle insertion into the jet flow. The calculated results are compared with data obtained with ignored influence of turbulent fluctuations on the motion and heat and mass transfer of the disperse phase. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 95–108, May–June, 2008.     

Slow Motion of a Porous Sphere Translating Along the Axis of a Circular Cylindrical Pore Subject to a Stress Jump Condition

The coupled flow problem of an incompressible axisymmetrical quasisteady motion of a porous sphere translating in a viscous fluid along the axis of a circular cylindrical pore is discussed using a combined analytical–numerical technique. At the fluid–porous interface, the stress jump boundary condition for the tangential stress along with continuity of normal stress and velocity components are employed. The flow through the porous particle is governed by the Brinkman model and the flow in the outside porous region is governed by Stokes equations. A general solution for the field equations in the clear region is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are satisfied first at the cylindrical pore wall by the Fourier transforms and then on the surface of the porous particle by a collocation method. The collocation solutions for the normalized hydrodynamic drag force exerted by the clear fluid on the porous particle is calculated with good convergence for various values of the ratio of radii of the porous sphere and pore, the stress jump coefficient, and a coefficient that is proportional to the permeability. The shape effect of the cylindrical pore on the axial translation of the porous sphere is compared with that of the particle in a spherical cavity; it found that the porous particle in a circular cylindrical pore in general attains a lower hydrodynamic drag than in a spherical envelope.     

Efficient approaches of determining the motion of a spherical particle in a swirling fluid flow using weighted residual methods

The motion of a spherical particle released in a swirling fluid flow is studied employing the least-squares method and method of moments. The governing equations are obtained and solved employing the two methods. The accuracy of the results is examined against the results of a fourth-order Runge–Kutta numerical method. The effects of various parameters, namely the initial radius, initial radial velocity, initial angular velocity, and drag-to-inertia ratio, on the non-dimensional velocity profiles and particle position distribution are considered. The results show that the radial velocity increases over time while the angular velocity decreases, and that an increase in the initial radial velocity increases the particle radial distance and angular velocity but decreases the radial velocity profile.     

Effect of Shape on Inertial Particle Dynamics in a Channel Flow

Particle dynamics in a channel flow are investigated using large eddy simulation and a Lagrangian particle tracking technique. Following validation of single-phase flow predictions against DNS results, fluid velocities are subsequently used to study the behaviour of particles of differing shape assuming one-way coupling between the fluid and the particles. The influence of shape- and orientation-dependent drag and lift forces on both the translational and rotational motion of the particles is accounted for to ensure accurate representation of the flow dynamics of non-spherical particles. The size of the particles studied was obtained based on an equivalent-volume sphere, and differing shapes were modelled using super-quadratic ellipsoid forms by varying their aspect ratio, with their orientation predicted using the incidence angle between the particle relative velocity and the particle principal axis. Results are presented for spherical, needle- and platelet-like particles at a number of different boundary layer locations along the wall-normal direction within the channel. The time evolution and probability density function of selected particle translational and rotational properties show a clear distinction between the behaviour of the various particles types, and indicate the significance of particle shape when modelling many practically relevant flows.     

Effects of vibrations on particle motion near a wall: Existence of attraction force

The effects of small vibrations on a particle oscillating near a solid wall in a fluid cell, relevant to material processing such as crystal growth in space, have been investigated experimentally and theoretically. Assuming the boundary layer around the particle to be thin compared to the particle radius at high vibration frequencies, an inviscid fluid model was developed to predict the motion of a spherical particle placed near a wall of a rectangular liquid-filled cell subjected to a sinusoidal vibration. Under these conditions, a non-uniform pressure distribution around the particle results in an average pressure that gives rise to an attraction force. Theoretical expressions for the attraction force are derived for the particle vibrating normal to and parallel with the nearest cell wall. The magnitude of this attractive force has been verified experimentally by measuring the motion of a steel particle suspended in the fluid cell by a thin wire. Experiments performed at high frequencies showed that the mean particle position, when the particle is brought near a cell wall, shifts towards the same wall, and is dependent on the cell amplitude and frequency, particle and fluid densities.     

Multiphase non‐Newtonian effects on pulsatile hemodynamics in a coronary artery

Hemodynamic stresses are involved in the development and progression of vascular diseases. This study investigates the influence of mechanical factors on the hemodynamics of the curved coronary artery in an attempt to identify critical factors of non‐Newtonian models. Multiphase non‐Newtonian fluid simulations of pulsatile flow were performed and compared with the standard Newtonian fluid models. Different inlet hematocrit levels were used with the simulations to analyze the relationship that hematocrit levels have with red blood cell (RBC) viscosity, shear stress, velocity, and secondary flow. Our results demonstrated that high hematocrit levels induce secondary flow on the inside curvature of the vessel. In addition, RBC viscosity and wall shear stress (WSS) vary as a function of hematocrit level. Low WSS was found to be associated with areas of high hematocrit. These results describe how RBCs interact with the curvature of artery walls. It is concluded that although all models have a good approximation in blood behavior, the multiphase non‐Newtonian viscosity model is optimal to demonstrate effects of changes in hematocrit. They provide a better stimulation of realistic blood flow analysis. Copyright © 2008 John Wiley & Sons, Ltd.     

The motion of rodlike particles in the pressure-driven flow between two flat plates

The motion of freely suspended rodlike particles has been observed in the pressure-driven flow between the two flat plates of a Hele Shaw flow cell at low Reynolds numbers. Data are reported for rodlike particles with aspect ratios of 12.0 suspended in a Newtonian fluid for gap thickness to particle length ratios of 3, 6, and 20; and for rodlike particles with aspect ratios between 5 and 8 in a non-Newtonian fluid (79.25 wt.% water, 20.2 wt.% glycerine, and 0.55 wt.% polyacrylamide). For the Newtonian fluid, the time-dependent orientation of the particles near and far from walls was shown to be in quantitative agreement with Jeffery's theory for ellipsoids suspended in a simple shear flow if an effective aspect ratio is calculated from the experimental period of rotation. Particles aligned with the flow direction and less than a particle half-length from a wall interacted irreversibly with the wall. For the non-Newtonian fluid, the timedependent orientation far from a wall was shown to be in qualitative agreement with Leal's theory for a second-order fluid; however, particles that were aligned with the flow direction and were near walls did not rotate.