Stability Analysis of Bioconvection of Gyrotactic Motile Microorganisms in a Fluid Saturated Porous Medium

Despite a large number of publications on bioconvection in suspensions of motile microorganisms, bioconvection in a fluid saturated porous medium is a relatively new area of research. This paper is motivated by experimental research by Kessler (1986) who established that a porous medium prevents the development of convection instability in algal suspensions. This suggests that there may exist a critical value of the permeability of a porous medium. If the permeability is smaller than critical, the system is stable and bioconvection does not develop. If the permeability is larger than critical, bioconvection may develop. This paper presents a model of bioconvection of gyrotactic motile microorganisms in a fluid saturated porous medium. The focus of this research is the determination of the critical value of permeability of a porous medium by a linear stability analysis. A simple but elegant analytical solution for the critical Darcy number is obtained.     

Instability Analysis of Gyrotactic Microorganisms: A Combined Effect of High-Frequency Vertical Vibration and Porous Media

A linear stability analysis is carried out to predict the instability analysis in a dilute suspension of gyrotactic microorganisms in horizontal fluid-saturated porous layer influenced by high-frequency vertical vibration. The governing equations, describing the mean flow, are the time-averaged Boussinesq equations and the analytical solution of the problem has been obtained using Galerkin method. A secular relation involving bioconvection Rayleigh number and its vibrational analogs and other parameters have been established. The graphical interpretations for dependence of bioconvection Rayleigh–Darcy number and corresponding wave number, on gyrotactic number and bioconvection Péclet number in the presence of vibration are utilized to understand the problem.     

Fibre suspension rheology: effect of concentration, aspect ratio and fibre size.

Viscosity data for fibre suspensions are produced using cone-and-plate geometry of enhanced dimensions for the reduced influence of fibre-wall interactions. Semi-concentrated suspensions of monodisperse polyamide fibres in silicone oil, with a variety of fibre concentrations (2, 5 and 8%), lengths and diameters, were studied. The suspension viscosity was measured in a range of shear stress in order to study the stress dependence. The study here focuses on the nature of the forces and interactions that contribute to the suspension viscosity. The results show that at sufficiently high stress levels, the suspension viscosity tends to reach a steady-state. At very low stress levels the suspension viscosity increases over time, most likely due to structures formed by adhesive forces. At higher concentrations, the viscosity depends on the absolute size of the fibres, again indicating the presence of non-hydrodynamic interactions.     

Direct numerical simulation of settling of a large solid particle during bioconvection

Settling of a large solid particle in bioconvection flow caused by gyrotactic microorganisms is investigated. The particle is released from the top of the bioconvection chamber; its settling pattern depends on whether it is released in the centre of the bioconvection plume or at its periphery. The Chimera method is utilized; a subgrid is generated around a moving particle. The method suggested by Liu and Wang (Comput. Fluid 2004; 33 :223–255) is further developed to account for the presence of a moving boundary in the streamfunction‐vorticity formulation using the finite‐difference method. A number of cases for different release positions of the particle are computed. It is demonstrated that bioconvection can either accelerate or decelerate settling of the particle depending on the initial position of the particle relative to the plume centre. It is also shown that the particle impacts bioconvection plume by changing its shape and location in the chamber. Copyright © 2005 John Wiley & Sons, Ltd.     

Overstability Analysis of Thermo-bioconvection Saturating a Porous Medium in a Suspension of Gyrotactic Microorganisms

A linear stability analysis is performed to analyze bioconvection in a dilute suspension of gyrotactic microorganisms in horizontal shallow fluid layer cooling from below and saturated by a porous medium, in the rigid boundary case. It is established that due to cooling from below thermally stratified layer is stabilized, which opposes the development of bioconvection and the situations for oscillatory convection may take place. The stability criterion is obtained in terms of thermal Rayleigh number, bioconvection Rayleigh number, gyrotactic number, bioconvection Peclet number, measure of cell eccentricity, Prandtl number, and Lewis number. It is observed that the presence of porous medium results in decrease of the magnitude of critical bioconvection Rayleigh number in comparison with its non-existence; hence due to porous effect, the system becomes less stable.     

The Dilute Rheology of Swimming Suspensions: A Simple Kinetic Model

A simple kinetic model is presented for the shear rheology of a dilute suspension of particles swimming at low Reynolds number. If interparticle hydrodynamic interactions are neglected, the configuration of the suspension is characterized by the particle orientation distribution, which satisfies a Fokker-Planck equation including the effects of the external shear flow, rotary diffusion, and particle tumbling. The orientation distribution then determines the leading-order term in the particle extra stress in the suspension, which can be evaluated based on the classic theory of Hinch and Leal (J Fluid Mech 52(4):683–712, 1972), and involves an additional contribution arising from the permanent force dipole exerted by the particles as they propel themselves through the fluid. Numerical solutions of the steady-state Fokker-Planck equation were obtained using a spectral method, and results are reported for the shear viscosity and normal stress difference coefficients in terms of flow strength, rotary diffusivity, and correlation time for tumbling. It is found that the rheology is characterized by much stronger normal stress differences than for passive suspensions, and that tail-actuated swimmers result in a strong decrease in the effective shear viscosity of the fluid.     

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall considering the interaction between squirmer and wall are numerically studied with an immersed boundary-lattice Boltzmann method. The power-law index, Reynolds number, initial orientation angle of squirmer, and initial distance of squirmer from the wall are all taken into account to investigate the swimming characteristics for pusher (β?<?0), neutral squirmer (β?=?0), and puller (β?>?0) (three kinds of swimmer types) near the no-slip boundary. Four new kinds of swimming modes are found. Results show that, for the pushers and pullers, the wall displays an increasing attraction with increasing power-law index n, which differs from the neutral squirmer who always departs from the wall after the first collision with the wall. Both the initial orientation angle and initial distance from the wall only affect the moving situations rather than the moving modes of the squirmers. However, the squirmers depart from the wall as the Reynolds number increases and chaotic orbits appear for some squirmers at Re?=?5. Several typical flow fields are analyzed and the power consumption and torque for different kinds of flows are also studied. It is found that, as the absolute value of β increases, the power consumption generally increases in shear-thinning (n?=?0.4), Newtonian (n?=?1), and shear-thickening (n?=?1.6) fluids. Moreover, the pushers (β?<?0) and the pullers (β?>?0) expend almost the same power if the absolute value of β remains the same. In addition, the power consumption of the squirmers is highly dependent on the power-law index n.     

Numerical simulation of flow kinematics and fiber orientation for multi-disperse suspension

The development of flow kinematics and fiber orientation distribution from the parabolic velocity profile and isotropic orientation at the channel inlet was computed in multi-disperse suspension flow through a parallel plate channel and their predictions were compared with those of mono- and bi-disperse suspensions. A statistical scheme (orientations of a large number of fibers are evaluated from the solution of the Jeffery equation along the streamlines) was confirmed to be very useful and feasible method to analyze accurately the orientation distribution of fibers in multi-disperse fiber suspension flow as well as mono- and bi-dispersions, instead of direct solutions of the orientation distribution function of fibers or the evolution equation of the orientation tensor which involves a closure equation. It was found that the flow kinematics and the fiber orientation depend completely on both the fiber aspect-ratio and the fiber parameter for multi-disperse suspension when the fiber–fiber and fiber-wall interactions are neglected. Furthermore, the addition of large aspect-ratio fibers as well as an increase in the fiber parameter related to the large aspect-ratio fibers could suppress the complex velocity field and stress distributions which are observed in suspensions containing small aspect-ratio fibers. From a practical point of view, therefore, the mechanical and physical properties of fiber composites should be improved with an increase in the volume fraction of large aspect-ratio fibers.     

Simulation of the rheological properties of suspensions of oblate spheroidal particles in a Newtonian fluid

A simulation algorithm was developed to predict the rheological properties of oblate spheroidal suspensions. The motion of each particle is described by Jeffery’s solution, which is then modified by the interactions between the particles. The interactions are considered to be short range and are described by results from lubrication theory and by approximating locally the spheroid surface by an equivalent spherical surface. The simulation is first tested on a sphere suspension, results are compared with known experimental and numerical data, and good agreement is found. Results are then presented for suspensions of oblate spheroids of two mean aspect ratios of 0.3 and 0.2. Results for the relative viscosity η _r, normal stress differences N ₁ and N ₂ are reported and compared with the few available results on oblate particle suspensions in a hydrodynamic regime. Evolution of the orientation of the particles is also observed, and a clear alignment with the flow is found to occur after a transient period. A change of sign of N ₁ from negative to positive as the particle concentration is increased is observed. This phenomenon is more significant as the particle aspect ratio increases. It is believed to arise from a change in the suspension microstructure as the particle alignment increases.     

Hydrodynamic interactions in squirmer motion: Swimming with a neighbour and close to a wall

We analyze the hydrodynamic coupling of pairs of squirmers and the impact it has on their short and long-time behavior. The study combines an analytic analysis of the hydrodynamic interactions between pairs of squirmers with computer simulations to elucidate the quantitative capabilities of the theoretical approach. The numerical study allows us to address the motion of simple geometries of squirmers on long times and perform a complete discussion of the effective repulsive interactions in squirmer ensembles. The contrast between analytic and numerical results identifies the features of active motion responsible for such effective interactions. The framework developed also allows for an analysis of the hydrodynamic coupling between a squirmer and a solid wall and shows the possibility of bounded motion next to a solid wall.     

Influence of suspension viscosity on Brownian relaxation of filler particles

Brownian relaxation caused by Brownian movement of particles in suspensions can macroscopically be probed by small-amplitude oscillatory shear experiments. Phenomenological considerations suggest a direct proportionality between suspension viscosity and Brownian relaxation times. To verify this relation experimentally, a set of nanocomposite suspensions with viscosities varying over five decades is presented. The suspensions are chosen in a way to ensure that particle-particle interactions and average particle-particle distances are identical so that they can be used as a model system to study the mere influence of suspension viscosity on Brownian relaxation. The suggested linear relationship between suspension viscosity and Brownian relaxation time can be confirmed. Moreover, a verification of a recently introduced characteristic timescale for Brownian relaxation is presented.     

Investigation of the onset of thermo-bioconvection in a suspension of oxytactic microorganisms in a shallow fluid layer heated from below

This paper investigates the combined effect of density stratification due to oxytactic upswimming and heating from below on the stability of a suspension of motile oxytactic microorganisms in a shallow fluid layer. Different from traditional bioconvection, thermo-bioconvection has two destabilizing mechanisms that contribute to creating the unstable density stratification. This problem may be relevant to a number of geophysical applications, such as the investigation of the dynamics of some species of thermophiles (heat loving microorganisms) living in hot springs. By performing a linear stability analysis, we obtained a correlation between the critical value of the bioconvection Rayleigh number and the traditional, “thermal” Rayleigh number. It is established that heating from below makes the system more unstable and helps the development of bioconvection.     

Shear-induced changes of electrical conductivity in suspensions

The effect of shear on electrical conductivity (rheo-conduction) is studied to give information about particle behaviour in suspensions. Past work is reviewed, and expressions are derived for the rheo-conduction of a suspension of nonconducting spheroids in a conducting matrix for current flow, parallel and normal to the suspension flow direction. A simple apparatus to study rheo-conduction in pipe flow is described, and measurements of steady and time-dependent effects are reported for various suspensions of colloidal particles. Suspensions of anisometric rod- and platelike particles at low concentrations showed rheo-conductive changes of sign, magnitude and relaxation that were consistent with the particle shape, concentration and interactions. The rheo-conductive response decreased with increasing volume fraction for platelike kaolinite particles, attributed to orientational jamming. Spherical latex particles gave unexpected rheo-conductive changes consistent with shear disruption of a conductive network of particles. It is concluded that rheo-conduction measurements are a useful adjunct to conventional rheometry.     

Strain-thickening transition in ferric-oxide suspensions under oscillatory shear

Shear-strain-thickening transition under oscillatory flow was observed in flocculated ferric-oxide suspensions in mineral oil. The value of the dynamic modulus of the suspensions that was measured at small strain amplitude after cessation of shear also became higher when the strain amplitude of the applied shear had been within or above the transition region.The ferric-oxide powders used were an acicular submicron maghemite (magnetic) and the hematite (non-magnetic) that was converted from the maghemite by heat treatment. The powders were treated with a dispersing agent and the suspensions were prepared in 33% by particle weight. The strain-thickening transition was observed in both the magnetic and the non-magnetic suspensions. However, the onset of the strain-thickening in the magnetic suspension was found at about one decade larger strain amplitude than that in the non-magnetic analog suspension, indicating particle interactions affect to the appearance of the phenomenon.A qualitative interpretation was made in view of site percolation for the enhancement of modulus at rest after the application of the large-amplitude oscillatory shear, where the process of the strain-thickening transition under shear and the development of the modulus after stopping the shear was described with a floc model in which the flocculation phase dilates as a result of the reduction of the particle linkages under higher shear.     

Dynamics of large solid particles in bioconvective sedimentation

Settling of one or two large solid particles in a bioconvection flow induced by gyrotactic motile microorganisms is investigated using a 2D numerical model. The results of varying the initial positions of large particles on the bioconvection flow pattern are investigated. The Chimera method is utilized to generate subgrids around the moving particles. It is demonstrated that the introduction of a single large particle displaces bioconvection plume and changes its shape. The introduction of two particles on the same side of the bioconvection plume further displaces the plume while the introduction of two particles on opposite sides reduces this displacement. The influence of the bioconvection plume on the particles' settling paths and particles' settling velocities is investigated. Copyright © 2006 John Wiley & Sons, Ltd.     

The Onset of Convection in a Suspension of Gyrotactic Microorganisms in Superimposed Fluid and Porous Layers: Effect of Vertical Throughflow

The purpose of this paper is to investigate the effect of vertical throughflow on the onset of bioconvection in a suspension of gyrotactic microorganisms. A dilute suspension of gyrotactic microorganisms in a shallow system that consists of superimposed fluid and porous layers is considered. A linear instability analysis of this problem is performed and the Galerkin method is utilized to solve the eigenvalue problem. The analysis leads to an equation for the critical Rayleigh number. It is shown that the vertical throughflow stabilizes the system.     

Stable and contact-free time stepping for dense rigid particle suspensions

We consider suspensions of rigid bodies in a two-dimensional viscous fluid. Even with high-fidelity numerical methods, unphysical contact between particles occurs because of spatial and temporal discretization errors. We extend a time stepping method that avoids overlap by imposing a minimum separation distance between all pairs of bodies. In its original form, the method discretizes interactions between different particles explicitly. Therefore, to avoid stiffness, a large minimum separation distance is used. In this paper, we introduce a new implicit time stepping method that is able to simulate dense suspensions with large time step sizes and a small minimum separation distance. The method is tested on various unbounded and bounded flows, and rheological properties of the resulting suspensions are computed.     

Direct simulation of flexible fibers

A numerical simulation of multiple flexible fibers in suspension in Newtonian simple shear flow is presented. The method used is similar to those of previous recent simulation works by Fan et al. [J. Non-Newtonian Fluid Mech. 74 (1998) 113] and Yamane et al. [J. Non-Newtonian Fluid Mech. 54 (1994) 405], however, the method has been modified to allow a small amount of bending and torsion in the fibers. A restoring moment acts to straighten the fibers as they interact in the flow.It is demonstrated that this simulation can be used to extract basic rheological information about the suspension including fiber orientations and suspension viscosity. The viscosity of semi-concentrated to concentrated flexible fiber suspensions are shown to increase by a magnitude of the order 7–10% greater than the equivalent rigid fiber suspension tested. This is in qualitative agreement with previous experimental work by Goto et al. [Rheologica Acta 25 (1986) 119] and Blakeney [J. Colloid Interface Sci. 22 (1966) 324]. The implication is that any constitutive relation involving particulate suspensions described by orientation vectors may quantitatively underestimate suspension viscosity, particularly for fibers of large aspect ratio, or low Young’s modulus, whereby the tendency to flex is greater [Rheologica Acta 25 (1986) 119]. If particulate deformation were accounted for (by whatever means) in the existing constitutive relationship, predictions of bulk suspension parameters such as viscosity should be noticeably improved. A method is developed to modify an existing rigid-fiber viscosity to an equivalent flexible fiber viscosity, hence improving viscosity prediction ability.     

Cover Image

A hybrid computational method coupling the lattice-Boltzmann (LB) method and a Langevin-dynamics (LD) method is developed to simulate nanoscale particle and polymer (NPP) suspensions in the presence of both thermal fluctuation and long-range many-body hydrodynamic interactions (HIs). Brownian motion of the NPP is explicitly captured by a stochastic forcing term in the LD method. The LD method is two-way coupled to the nonfluctuating LB fluid through a discrete LB forcing source distribution to capture the long-range HI. To ensure intrinsically linear scalability with respect to the number of particles, a Eulerian-host algorithm for short-distance particle neighbor search and interaction is developed and embedded to LB-LD framework. The validity and accuracy of the LB-LD approach are demonstrated through several sample problems. The simulation results show good agreements with theory and experiment. The LB-LD approach can be favorably incorporated into complex multiscale computational frameworks for efficiently simulating multiscale multicomponent particulate suspension systems such as complex blood suspensions.     

Yield stress and maximum packing fraction of concentrated suspensions

The yield stress and features of the structure of concentrated suspensions based on silica flour, with particles of average diameter around 4 m, were investigated in terms of a phenomenological model. The yield stress of a concentrated suspension of known solid volume concentration is estimated by employing a shear-dependent maximum packing fraction _m which is obtained by model fitting equilibrium viscosity data, and by incorporating a first-order kinetic equation. The model proposed was examined by using several mineral suspensions in which silica flour was mixed with metal oxide particles so that microstructural features of the suspensions could be adjusted. A cocoa fat suspension was also used as a test sample having radically different chemistry. The agreement between the model prediction and independently obtained experimental evidence is acceptable. Furthermore, a qualitative explanation is obtained by a scaling analysis in an effort to relate the model parameters with the suspension structure that stems from interactions among the suspension constituents.