SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Nanofluids.

Open Access Nano Express

Nanofluid bioconvection in water-based suspensions containing nanoparticles and oxytactic microorganisms: oscillatory instability

Andrey V Kuznetsov

Author affiliations

Dept. of Mechanical and Aerospace Engineering, North Carolina State University, Campus Box 7910, Raleigh, NC 27695-7910, USA

Citation and License

Nanoscale Research Letters 2011, 6:100  doi:10.1186/1556-276X-6-100

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/6/1/100

Received:20 September 2010
Accepted:25 January 2011
Published:25 January 2011

© 2011 Kuznetsov; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The aim of this article is to propose a novel type of a nanofluid that contains both nanoparticles and motile (oxytactic) microorganisms. The benefits of adding motile microorganisms to the suspension include enhanced mass transfer, microscale mixing, and anticipated improved stability of the nanofluid. In order to understand the behavior of such a suspension at the fundamental level, this article investigates its stability when it occupies a shallow horizontal layer. The oscillatory mode of nanofluid bioconvection may be induced by the interaction of three competing agencies: oxytactic microorganisms, heating or cooling from the bottom, and top or bottom-heavy nanoparticle distribution. The model includes equations expressing conservation of total mass, momentum, thermal energy, nanoparticles, microorganisms, and oxygen. Physical mechanisms responsible for the slip velocity between the nanoparticles and the base fluid, such as Brownian motion and thermophoresis, are accounted for in the model. An approximate analytical solution of the eigenvalue problem is obtained using the Galerkin method. The obtained solution provides important physical insights into the behavior of this system; it also explains when the oscillatory mode of instability is possible in such system.


The term "nanofluid" was coined by Choi in his seminal paper presented in 1995 at the ASME Winter Annual Meeting [1]. It refers to a liquid containing a dispersion of submicronic solid particles (nanoparticles) with typical length on the order of 1-50 nm [2]. The unique properties of nanofluids include the impressive enhancement of thermal conductivity as well as overall heat transfer [3-7]. Various mechanisms leading to heat transfer enhancement in nanofluids are discussed in numerous publications; see, for example [8-12].

Wang [13-15] pioneered in developing the constructal approach, created by Bejan [16-19], for designing nanofluids. Nanofluids enhance the thermal performance of the base fluid; the utilization of the constructal theory makes it possible to design a nanofluid with the best microstructure and performance within a specified type of microstructures.

Recent publications show significant interest in applications of nanofluids in various types of microsystems. These include microchannels [20], microheat pipes [21], microchannel heat sinks [22], and microreactors [23]. There is also significant potential in using nanomaterials in different bio-microsystems, such as enzyme biosensors [24]. In [25], the performance of a bioseparation system for capturing nanoparticles was simulated. There is also strong interest in developing chip-size microdevices for evaluating nanoparticle toxicity; Huh et al. [26] suggested a biomimetic microsystem that reconstitutes the critical functional alveolar-capillary interface of the human lung to evaluate toxic and inflammatory responses of the lung to silica nanoparticles.

The aim of this article is to propose a novel type of a nanofluid that contains both nanoparticles and oxytactic microorganisms, such as a soil bacterium Bacillus subtilis. These particular microorganisms are oxygen consumers that swim up the oxygen concentration gradient. There are important similarities and differences between nanoparticles and motile microorganisms. In their impressive review of nanofluids research, Wang and Fan [27] pointed out that nanofluids involve four scales: the molecular scale, the microscale, the macroscale, and the megascale. There is interaction between these scales. For example, by manipulating the structure and distribution of nanoparticles the researcher can impact macroscopic properties of the nanofluid, such as its thermal conductivity. Similar to nanofluids, in suspensions of motile microorganisms that exhibit spontaneous formation of flow patterns (this phenomenon is called bioconvection) physical laws that govern smaller scales lead to a phenomenon visible on a larger scale. While superfluidity and superconductivity are quantum phenomena visible at the macroscale, bioconvection is a mesoscale phenomenon, in which the motion of motile microorganisms induces a macroscopic motion (convection) in the fluid. This happens because motile microorganisms are heavier than water and they generally swim in the upward direction, causing an unstable top-heavy density stratification which under certain conditions leads to the development of hydrodynamic instability. Unlike motile microorganisms, nanoparticles are not self-propelled; they just move due to such phenomena as Brownian motion and thermophoresis and are carried by the flow of the base fluid. On the contrary, motile microorganisms can actively swim in the fluid in response to such stimuli as gravity, light, or chemical attraction. Combining nanoparticles and motile microorganisms in a suspension makes it possible to use benefits of both of these microsystems.

One possible application of bioconvection in bio-microsystems is for mass transport enhancement and mixing, which are important issues in many microsystems [28,29]. Also, the results presented in [30] suggest using bioconvection in a toxic compound sensor due to the ability of some toxic compounds to inhibit the flagella movement and thus suppress bioconvection. Also, preventing nanoparticles from agglomerating and aggregating remains a significant challenge. One of the reasons why this is challenging is because although inducing mixing at the macroscale is easy and can be achieved by stirring, inducing and controlling mixing at the microscale is difficult. Bioconvection can provide both types of mixing. Macroscale mixing is provided by inducing the unstable density stratification due to microorganisms' upswimming. Mixing at the microscale is provided by flagella (or flagella bundle) motion of individual microorganisms. Due to flagella rotation, microorganisms push fluid along their axis of symmetry, and suck it from the sides [31]. While the estimates given in [32] show that the stresslet stress produced by individual microorganisms have negligible effect on macroscopic motion of the fluid (which is rather driven by the buoyancy force induced by the top-heavy density stratification due to microorganisms' upswimming), the effect produced by flagella rotation is not negligible on the microscopic scale (on the scale of a microorganism and a nanoparticle).

In order to use suspensions containing both nanoparticles and motile microorganisms in microsystems, the behavior of such suspensions must be understood at the fundamental level. Bio-thermal convection caused by the combined effect of upswimming of oxytacic microorganisms and temperature variation was investigated in [33-36]. Bioconvection in nanofluids is expected to occur if the concentration of nanoparticles is small, so that nanoparticles do not cause any significant increase of the viscosity of the base fluid. The problem of bioconvection in suspensions containing small solid particles (nanoparticles) was first studied in [37-41] and then recently in [42]. Non-oscillatory bioconvection in suspensions of oxytactic microorganisms was considered in Kuznetsov AV: Nanofluid bioconvection: Interaction of microorganisms oxytactic upswimming, nanoparticle distribution and heating/cooling from below. Theor Comput Fluid Dyn 2010, submitted. This article extends the theory to the case of oscillatory convection in suspensions containing both nanoparticles and oxytactic microorganisms.

Governing equations

The governing equations are formulated for a water-based nanofluid containing nanoparticles and oxytactic microorganisms. The nanofluid occupies a horizontal layer of depth H. It is assumed that the nanoparticle suspension is stable. According to Choi [2], there are methods (including suspending nanoparticles using either surfactant or surface charge technology) that lead to stable nanofluids. It is further assumed that the presence of nanoparticles has no effect on the direction of microorganisms' swimming and on their swimming velocity. This is a reasonable assumption if the nanoparticle suspension is dilute; the concentration of nanoparticles has to be small anyway for the bioconvection-induced flow to occur (otherwise, a large concentration of nanoparticles would result in a large suspension viscosity which would suppress bioconvection).

In formulating the governing equations, the terms pertaining to nanoparticles are written using the theory developed in Buongiorno [43], while the terms pertaining to oxytactic microorganisms are written using the approach developed by Hillesdon and Pedley [44,45].

The continuity equation for the nanoparticle-microorganism suspension considered in this research is

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M1">View MathML</a>


where U = (u,v,w) is the dimensionless nanofluid velocity, defined as U*H/αf; U* is the dimensional nanofluid velocity; αf is the thermal diffusivity of a nanofluid, k/(ρc)f; k is the thermal conductivity of the nanofluid; and (ρc)f is the volumetric heat capacity of the nanofluid. The dimensionless coordinates are defined as (x,y,z) = (x*, y*, z*)/H, where z is the vertically downward coordinate.

The buoyancy force can be considered to be made up of three separate components that result from: the temperature variation of the fluid, the nanoparticle distribution (nanoparticles are heavier than water), and the microorganism distribution (microorganisms are also heavier than water). Utilizing the Boussinesq approximation (which is valid because the inertial effects of the density stratification are negligible, the dominant term multiplying the inertia terms is the density of the base fluid that exceeds by far the density stratification), the momentum equation can be written as:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M2">View MathML</a>


where <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M76">View MathML</a> is the vertically downward unit vector.

The dimensionless variables in Equation 2 are defined as:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M3">View MathML</a>


where t is the dimensionless time, p is the dimensionless pressure, ϕ is the relative nanoparticle volume fraction, T is the dimensionless temperature, n is the dimensionless concentration of microorganisms, t* is the time, p* is the pressure, μ is the viscosity of the suspension (containing the base fluid, nanoparticles and microorganisms), ϕ* is the nanoparticle volume fraction, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M77">View MathML</a> is the nanoparticle volume fraction at the lower wall, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M78">View MathML</a> is the nanoparticle volume fraction at the upper wall, T* is the nanofluid temperature, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M79">View MathML</a> is the temperature at the upper wall (also used as a reference temperature), <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M80">View MathML</a> is the temperature at the lower wall, n* is the concentration of microorganisms, and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M81">View MathML</a> is the average concentration of microorganisms (concentration of microorganisms in a well-stirred suspension).

The dimensionless parameters in Equation 2, namely, the Prandtl number, Pr; the basic-density Rayleigh number, Rm; the traditional thermal Rayleigh number, Ra; the nanoparticle concentration Rayleigh number, Rn; the bioconvection Rayleigh number, Rb; and the bioconvection Lewis number, Lb, are defined as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M4">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M5">View MathML</a>


where ρf0 is the base-fluid density at the reference temperature; ρp is the nanoparticle mass density; g is the gravity; β is the volumetric thermal expansion coefficient of the base fluid; Δρ is the density difference between microorganisms and a base fluid, ρmo - ρf0; ρmo is the microorganism mass density; θ is the average volume of a microorganism; and Dmo is the diffusivity of microorganisms (in this model, following [44,45], all random motions of microorganisms are simulated by a diffusion process).

The conservation equation for nanoparticles contains two diffusion terms on the right-hand side, which represent the Brownian diffusion of nanoparticles and their transport by thermophoresis (a detailed derivation is available in [43,46]):

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M6">View MathML</a>


In Equation 6, the nanoparticle Lewis number, Ln, and a modified diffusivity ratio, NA (this parameter is somewhat similar to the Soret parameter that arises in cross-diffusion phenomena in solutions), are defined as:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M7">View MathML</a>


where DB is the Brownian diffusion coefficient of nanoparticles and DT is the thermophoretic diffusion coefficient.

The right-hand side of the thermal energy equation for a nanofluid accounts for thermal energy transport by conduction in a nanofluid as well as for the energy transport because of the mass flux of nanoparticles (again, a detailed derivation is available in [43,46]):

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M8">View MathML</a>


In Equation 8, NB is a modified particle-density increment, defined as:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M9">View MathML</a>


where (ρc)p is the volumetric heat capacity of the nanoparticles.

The right-hand side of the equation expressing the conservation of microorganisms describes three modes of microorganisms transport: due to macroscopic motion (convection) of the fluid, due to self-propelled directional swimming of microorganisms relative to the fluid, and due diffusion, which approximates all stochastic motions of microorganisms:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M10">View MathML</a>


where V is the dimensionless swimming velocity of a microorganism, V*H/αf, which is calculated as [44,45]:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M11">View MathML</a>


In Equation 11 <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M82">View MathML</a> is the Heaviside step function and C is the dimensionless oxygen concentration, defined as:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M12">View MathML</a>


where C* is the dimensional oxygen concentration, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M83">View MathML</a> is the upper-surface oxygen concentration (the upper surface is assumed to be open to atmosphere), and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M84">View MathML</a> is the minimum oxygen concentration that microorganisms need to be active. Equation 11 thus assumes that microorganisms swim up the oxygen concentration gradient and that their swimming velocity is proportional to that gradient; however, in order for microorganisms to be active the oxygen concentration need to be above <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M84">View MathML</a>. Since this article deals with a shallow layer situation, it is assumed that <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M85">View MathML</a> throughout the layer thickness, and the Heaviside step function, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M86">View MathML</a>, in Equation 11 is equal to unity.

Also, the bioconvection Péclet number, Pe, in Equation 11 is defined as:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M13">View MathML</a>


where b is the chemotaxis constant (which has the dimension of length) and Wmo is the maximum swimming speed of a microorganism (the product bWmo is assumed to be constant).

Finally, the oxygen conservation equation is:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M14">View MathML</a>


The first term on the right-hand side of Equation 14 represents oxygen diffusion, while the second term represents oxygen consumption by microorganisms.

The new dimensionless parameters in Equation 14 are

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M15">View MathML</a>


where Le is the traditional Lewis number, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M87">View MathML</a> is the dimensionless parameter describing oxygen consumption by the microorganisms, DS is the diffusivity of oxygen, and γ is a dimensional constant describing consumption of oxygen by the microorganisms.

According to Hillesdon and Pedley [45], the layer can be treated as shallow as long as the following condition is satisfied:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M16">View MathML</a>


Equation 16 gives the maximum layer depth for which the oxygen concentration at the bottom does not drop below <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M84">View MathML</a>.

The boundary conditions for Equations 1, 2, 6, 8, 10, and 14 are imposed as follows. It is assumed that the temperature and the volumetric fraction of the nanoparticles are constant on the boundaries and the flux of microorganisms through the boundaries is equal to zero. The lower boundary is always assumed rigid and the upper boundary can be either rigid or stress-free. The boundary conditions for case of a rigid upper wall are

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M17">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M18">View MathML</a>


The fifth equation in (18) is equivalent to the statement that the total flux of microorganisms at the upper surface is equal to zero: the microorganisms swim vertically upward at the top surface but (because their concentration gradient at the top surface is directed vertically upward) they are simultaneously pushed downward by diffusion; the two fluxes are equal but opposite in direction).

If the upper surface is stress-free, the second equation in (18) is replaced with the following equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M19">View MathML</a>


Basic state

The solution for the basic state corresponds to a time-independent quiescent situation. The solution is of the following form:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M20">View MathML</a>


In this case, the solution of Equations 6, 8, 10, and 14 subjects to boundary conditions (17) and (18) is (the particular form of hydrodynamic boundary conditions at the upper surface is not important because the solution in the basic state is quiescent):

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M21">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M22">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M23">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M24">View MathML</a>


where A1 is the smallest positive root of the transcendental equation

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M25">View MathML</a>


The solutions given by Equations 23 and 24 were first reported in [44].

The pressure distribution in the basic state, pb (z), can then be obtained by integrating the following form of the momentum equation (which follows from Equation 2):

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M26">View MathML</a>


Equations 21 and 22 can be simplified if characteristic parameter values for a typical nanofluid are considered. Based on the data presented in Buongiorno [43] for an alumina/water nanofluid, the following dimensional parameter values are utilized: <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M88','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M88">View MathML</a>, αf = 2 × 10-7m2/s, DB = 4 × 10-11m2/s, μ = 10-3 Pas, and ρf0 = 103 kg/m3. The thermophoretic diffusion coefficient, DT, is estimated as <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M89">View MathML</a>, where, according to Buongiorno [43], τ is estimated as 0.006. This results in DT = 6 × 10-11m2/s. The nanoparticle Lewis number is then estimated as Ln = 5.0 × 103. The modified diffusivity ratio, NA, and the modified particle-density increment, NB, depend on the temperature difference between the lower and the upper plates and on the nanoparticle fraction decrement. Assuming that <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M90">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M91','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M91">View MathML</a>, and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M92','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M92">View MathML</a>, gives the following estimates: NA = 5 and NB = 7.5 × 10-4. This suggests that the exponents in Equations 21 and 22 are small and that these equations can be simplified as:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M27">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M28">View MathML</a>


Linear instability analysis

Perturbations are superimposed on the basic solution, as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M29">View MathML</a>


Equation 29 is then substituted into Equations 1, 2, 6, 8, 10, and 14, the resulting equations are linearized and the use is made of Equations 27 and 28. This procedure results in the following equations for the perturbation quantities:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M30">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M31">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M32">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M33">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M34">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M35">View MathML</a>


Equations 30 to 35 are independent of Rm since this parameter is just a measure of the basic static pressure gradient. In order to eliminate the pressure and horizontal components of velocity from Equations 30 and 31, Equation 31 (see [46]) is operated with <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M93','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M93">View MathML</a> and the use is made of the identity curl curl ≡ grad div - ∇2 together with Equation 30. This results in the reduction of Equations 30 and 31 to the following scalar equation which involves only one component of the perturbation velocity, w':

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M36">View MathML</a>


where <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M94','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M94">View MathML</a> is the two-dimensional Laplacian operator in the horizontal plane and ∇4w' is the Laplacian of the Laplacian of w'.

Equations 17 and 18 then lead to the following boundary conditions for the perturbation quantities for the case when both the lower and upper walls are rigid:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M37">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M38">View MathML</a>


If the upper boundary is stress-free, the second equation in Equation 38 is replaced by

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M39">View MathML</a>


The method of normal modes is used to solve a linear boundary-value problem composed of differential Equations 32 to 36 and boundary conditions (37), (38) (or (39)). A normal mode expansion is introduced as:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M40">View MathML</a>


where the function f(x,y) satisfies the following equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M41">View MathML</a>


and m is the dimensionless horizontal wavenumber.

Substituting Equation 40 into Equations 36 and 32 to 35, utilizing Equation 41, and letting <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M95','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M95">View MathML</a> (so that the resulting equation for amplitudes would depend on the product <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M96">View MathML</a> rather than on Pe and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M87">View MathML</a> individually), the following equations for the amplitudes, W, Θ, Φ, N, and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M97','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M97">View MathML</a>, are obtained:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M42">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M43">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M44">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M45">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M46">View MathML</a>


where Equation 25 for A1 is reduced to

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M47">View MathML</a>


In Equations 42 to 46 s is a dimensionless growth factor; for neutral stability the real part of s is zero, so it is written s = , where ω is a dimensionless frequency (it is a real number).

For the case of rigid-rigid walls, the boundary conditions for the amplitudes are

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M48">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M49">View MathML</a>


If the upper surface is stress-free, the second equation in (49) is replaced by

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M50">View MathML</a>


Equations 42 to 46 are solved by a single-term Galerkin method. For the case of the rigid-rigid boundaries, the trial functions, which satisfy the boundary conditions given by Equations 48 and 49, are

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M51">View MathML</a>



<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M52">View MathML</a>


and A1 is given by Equation 47.

If the upper boundary is stress-free, W1 is replaced by

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M53">View MathML</a>


and the rest of the trial functions are still given by Equation 51. W1 given by Equation 53 satisfies the boundary condition given by Equation 50.

Results and discussion

Rigid-rigid boundaries

For the case of the rigid-rigid boundaries the utilization of a standard Galerkin procedure (see, for example [47]), which involves substituting the trial functions given by Equation 51 into Equations 42 to 46, calculating the residuals, and making the residuals orthogonal to the relevant trial functions, results in the following eigenvalue equation relating three Rayleigh numbers, Ra, Rn, and Rb:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M54">View MathML</a>


where functions F1, F2, F3, and F4 are given in the appendix [see Equations A1 to A4], they depend on Lb, Le, Ln, Pr, NA, ϖ, ω, and m. It is interesting that Equation 54 is independent of NB at this order (one-term Galerkin) of approximation.

In order to evaluate the accuracy of the one-term Galerkin approximation used in obtaining Equation 54 the accuracy of this equation is estimated for the case of non-oscillatory instability (which corresponds to ω = 0) for the situation when the suspension contains no microorganisms (this corresponds to <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M98','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M98">View MathML</a>, which leads to Rb = 0) and no nanoparticles (this leads to Rn = 0).

In this limiting case Equation 54 collapses to

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M55">View MathML</a>


The right-hand side of Equation 55 takes the minimum value of 1750 at mc = 3.116; the obtained critical value of Ra is 2.5% greater than the exact value (1707.762) for this problem reported in [48]. The corresponding critical value of the wavenumber is 0.03% smaller than the exact value (3.117) reported in [48].

Based on the data presented in [44,45] for soil bacterium Bacillus subtilis, the following parameter values for these microorganisms are used: Dm = 1.3 × 10-10 m2/s, Ds = 2.12 × 10-9 m2/s, Δρ = 100 kg/m3, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M99','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M99">View MathML</a>, θ = 10-18 m3, and H = 2.5 × 10-3 m (or 2.5 mm, this is a typical depth of a shallow layer; this size is also typical for a microdevice). Also, according to Hillesdon et al. [45], for Bacillus subtilis dimensionless parameters can be estimated as follows: Pe = 15H, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M100','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M100">View MathML</a>, where the layer depth, H, must be given in mm. Based on [43], the following parameter values for a typical alumina/water nanofluid are utilized: <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M101','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M101">View MathML</a>, ρf0 = 103 kg/m3, ρp = 4 × 103 kg/m3, (ρc)p = 3.1 × 106 J/m3, αf = 2 × 10-7 m2/s, DB = 4 × 10-11 m2/s, DT = 6 × 10-11 m2/s, and μ = 10-3 Pas. It is also assumed that <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M102','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M102">View MathML</a>, β = 3.4 × 10-31/K, (ρC)f = 4 × 106J/m3, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M103','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M103">View MathML</a>, and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M104','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M104">View MathML</a>.

The parameter values given above result in the following representative values of dimensionless parameters: Lb = 1.5 × 103, Le = 94, Ln = 5.0 × 103, Pr = 5.0, NA = 5, NB = 7.5 × 10-4, Pe = 37, <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M105','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M105">View MathML</a>, ϖ = 17, Ra = 2.7 × 103, Rb = 1.2 × 105, Rm = 8.0 × 105, and Rn = 2.3 × 103. The values of Ra and Rb can be controlled by changing the temperature difference between the plates and the microorganism concentration, respectively, and Rn depends on nanoparticle concentrations at the boundaries.

For Figure 1a,b,c, the following values of dimensionless parameters are utilized: Lb = 1500, Le = 94, Ln = 5000, Pr = 5, NA = 5, ϖ = 17, and Rb = 0 (which corresponds to the situation with zero concentration of microorganisms). Rn is changing in the range between -1.2 and 1.2. In Figure 1a, the boundary for non-oscillatory instability (shown by a solid line) is obtained by setting ω to zero in Equation 54, solving this equation for Ra and then finding the minimum with respect to m of the right-hand side of the obtained equation. The boundary for oscillatory instability (shown by a dotted line) is obtained by the following procedure. Two coupled equations are produced by taking the real and imaginary parts of Equation 54. One of these equations is used to eliminate ω, and the resulting equation is then solved for Ra; the critical value of Ra is again obtained by calculating the minimum value that the expression for Ra takes with respect to m.

thumbnailFigure 1. The case of rigid upper and lower walls, Rb = 0 (no microorganisms): (a) Oscillatory and non-oscillatory instability boundaries in the (Rac, Rn) plane. (b) Critical wavenumber in the (Rac, Rn) plane. (c) Square of the oscillation frequency, ω2, versus the nanoparticle concentration Rayleigh number (for oscillatory instability to occur, ω2 must be positive so that ω remains real).

Figure 1a shows that for Rb = 0 the curve representing the instability boundary for non-oscillatory convection (solid line) is a straight line in the (Rac, Rn) plane. Rn is defined in Equation 5 in such a way that positive Rn corresponds to a top-heavy nanoparticle distribution. Therefore, the increase of Rn produces the destabilizing effect and reduces the critical value of Ra. A comparison between instability boundaries for non-oscillatory (solid line) and oscillatory (dotted line) cases indicates that in order for the oscillatory instability to occur, Rn generally must be negative, which corresponds to a bottom-heavy (stabilizing) nanoparticle distribution. In this case the destabilizing effect of the temperature gradient (positive Ra corresponds to heating from the bottom) and destabilizing effect from upswimming of oxytactic microorganisms compete with the stabilizing effect of the nanoparticle distribution.

Figure 1b shows that the critical value of the wavenumber, mc, is independent of Rn and for the case displayed in Figure 1a (Rb = 0) is equal to 3.116; also, it is almost independent of the mode of instability (non-oscillatory versus oscillatory).

Figure 1c shows the square of the oscillation frequency, ω2, versus the nanoparticle concentration Rayleigh number, Rn. The value of ω2 for the oscillatory instability boundary is obtained by eliminating Ra from the two coupled equations resulting from taking the real and imaginary parts of Equation 54 and solving the resulting equation for ω2. The solution is presented in terms of ω2 rather than ω because the resulting equation is bi-quadratic in ω. For oscillatory instability to occur, ω2 must be positive so that ω is real. Figure 1c shows that for Rb = 0 ω is real when Rn is negative.

Figure 2a,b,c is computed for the same parameter values as Figure 1a,b,c, but now with Rb = 120000. Figure 2a,b,c thus shows the effect of microorganisms. By comparing Figure 2a with 1a, it is evident that the presence of microorganisms produces the destabilizing effect and reduces the critical value of Ra. For example, at (NA + Ln) Rn = -5000 in Figure 1a the value of Rac corresponding to the non-oscillatory instability boundary is 6750 and in Figure 2a the corresponding value of Rac is 6437. At (NA + Ln) Rn = 5000 in Figure 1a the value of Rac corresponding to the non-oscillatory instability boundary is -3250 and in Figure 2a the corresponding value of Rac is -3563. The destabilizing effect of oxytactic microorganisms is explained as follows. These microorganisms are heavier than water and on average they swim in the upward direction. Therefore, the presence of microorganisms produces a top-heavy density stratification and contributes to destabilizing the suspension.

thumbnailFigure 2. Similar to Figure 1, but now with Rb = 120000.

The comparison of Figure 2b with 1b shows that the presence of microorganisms increases the critical wavenumber (in Figure 1b it was 3.116 and in Figure 2b it is 3.441).

Figure 2c brings an interesting insight. Apparently, if the concentration of microorganisms is above a certain value, the oscillatory mode of instability is not possible. Indeed, ω2 in Figure 2c is negative for the whole range of Rn (-1.2 ≤ Rn ≤ 1.2) used for computing this figure. This means that ω is imaginary and oscillatory instability does not occur for the value of Rb used in computing Figure 2.

Rigid-free boundaries

For the case when the upper boundary is stress-free, the eigenvalue equation is

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M56">View MathML</a>


where functions F5, F6, F7, and F8 are given in the appendix [see Equations A10 to A13].

Again, to evaluate of the accuracy of the one-term Galerkin approximation in this case, the accuracy of Equation 56 is estimated for the case of non-oscillatory instability (which corresponds to ω = 0) for the situation when the suspension contains no microorganisms (Rb = 0) and no nanoparticles (Rn 0). In this limiting case Equation 56 collapses to

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M57">View MathML</a>


The right-hand side of Equation 57 takes the minimum value of 1139 at mc=2.670; the obtained value of Rac is 3.48% greater than the exact value (1100.65) for this problem reported in [48]. The corresponding critical value of the wavenumber is 0.45% smaller than the exact value (2.682) reported in [48].

For Figures 3a,b,c and 4a,b,c, which show the results for the rigid-free boundaries, the same parameter values as for Figures 1 and 2 are utilized. Figure 3a, which is computed for Rb = 0 (no microorganisms), shows boundaries of non-oscillatory and oscillatory instabilities. This figure is similar to Figure 1a, but since now the case of the rigid-free boundaries is considered, the values of the critical Rayleigh number in Figure 3a are smaller than those in Figure 1a. Again, the comparison between the non-oscillatory and oscillatory instability boundaries indicates that in order for oscillatory instability to occur Rn must be negative; in this case at the instability boundary the effect of the nanoparticle distribution is stabilizing and the effect of the temperature gradient is destabilizing; the presence of these two competing agencies makes the oscillatory instability possible.

thumbnailFigure 3. The case of a rigid lower wall and a stress-free upper wall, Rb = 0 (no microorganisms): (a) Oscillatory and non-oscillatory instability boundaries in the (Rac, Rn) plane. (b) Critical wavenumber in the (Rac, Rn) plane. (c) Square of the oscillation frequency, ω2, versus the nanoparticle concentration Rayleigh number (for oscillatory instability to occur, ω2 must be positive so that ω remains real).

The critical wavenumber shown in Figure 3b (mc = 2.670) is smaller than the corresponding critical wavenumber for the rigid-rigid boundaries shown in Figure 1b. Again, it is independent of Rn and almost independent of the mode of instability (non-oscillatory versus oscillatory).

Figure 3c, similar to Figure 1c, shows that ω is real when Rn is negative, which means that for negative values of Rn oscillatory instability is indeed possible.

Figure 4a,b,c shows the results for rigid-free boundaries computed with Rb = 120000, meaning that the difference with Figure 3a,b,c is the presence of microorganisms. As in the case with rigid-rigid boundaries, the presence of microorganisms produces a destabilizing effect and reduces the critical value of the Rayleigh number (compare Figures 4a and 3a).

thumbnailFigure 4. Similar to Figure 3, but now with Rb = 120000.

Also, the presence of microorganisms increases the critical value of the wavenumber (compare Figures 4b and 3b).

Figure 4c again shows that for the range of Rn used for this figure the presence of microorganisms makes the oscillatory mode of instability impossible (corresponding values of ω are imaginary).


The possibility of oscillatory mode of instability in a nanofluid suspension that contains oxytactic microorganisms is investigated. Since these microorganisms swim up the oxygen concentration gradient, toward the free surface (which is open to the air), and they are heavier than water, they always produce the destabilising effect on the suspension. The destabilizing effect of microorganisms is larger if their concentration in the suspension is larger. The concentration of microorganisms is measured by the bioconvection Rayleigh number, Rb, which by definition is always non-negative (the zero value of Rb corresponds to a suspension with no microorganisms). The increase of Rb thus destabilizes the suspension. It is also shown that the presence of microorganisms increases the critical wavenumber.

The effect of the temperature distribution can be either stabilizing (heating from the top, negative thermal Rayleigh number Ra) or destabilizing (heating from the bottom, positive Ra). The effect of nanoparticles can also be stabilizing (bottom-heavy nanoparticle distribution, negative nanoparticle concentration Rayleigh number Rn) or destabilizing (top-heavy nanoparticle distribution, positive Rn).

The results obtained in this article indicate that in order for the oscillatory instability to occur, Rn generally must be negative, which corresponds to a bottom-heavy (stabilizing) nanoparticle distribution. In this case the destabilizing effect of the temperature gradient (positive Ra) and destabilizing effect from upswimming of oxytactic microorganisms compete with the stabilizing effect of the nanoparticle distribution.

In order for the oscillatory mode of instability to occur, the dimensionless oscillation frequency, ω, must be real. Since increasing Rb pushes ω2 to negative values, oscillatory instability is possible only if the concentration of microorganisms is below a certain value.

The results for the rigid-rigid and rigid-free boundaries are similar, but the critical Rayleigh number for the rigid-free boundaries is smaller. The critical wavenumber for the rigid-free boundaries can be either smaller or larger, depending on the concentration of microorganisms. For Rb = 0 the critical wavenumber is smaller for the rigid-free boundaries but for Rb = 120000 it is larger than for the rigid-rigid boundaries.


The functions F1, F2, F3, and F4 defining the eigenvalue equation for the layer with the rigid-rigid boundaries [given by Equation 54] are

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M58">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M59">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M60">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M61">View MathML</a>


The integrals I1 to I5 in Equations A1 to A4 are functions of Le and ϖ. The expressions for these integrals for the rigid upper boundary case are given below:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M62">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M63">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M64">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M65">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M66">View MathML</a>


The functions F5, F6, F7, and F8 defining the eigenvalue equation for the layer with the rigid-free boundaries [given by Equation 56] are

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M67">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M68">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M69">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M70">View MathML</a>


The integrals <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M106','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M106">View MathML</a> to <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M107','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M107">View MathML</a> in Equations A10 to A13 are functions of Le and ϖ. The expressions for these integrals for the stress-free upper boundary case are given below:

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M71">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M72">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M73">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M74">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/100/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/100/mathml/M75">View MathML</a>


Competing interests

The author declares that he has no competing interests.

Authors' contributions

AVK carried out all the work regarding the development of the model, performing simulations, writing and revising the paper and approving the final manuscript.


  1. Choi SUS: Enhancing thermal conductivity of fluids with nanoparticles. In Developments and Applications of Non-Newtonian Flows. Volume 99. Edited by Siginer DA, Wang HP. New York: ASME; 1995. OpenURL

  2. Choi SUS: Nanofluids: From vision to reality through research.

    J Heat Transf Trans ASME 2009, 131:033106. Publisher Full Text OpenURL

  3. Lee S, Choi SUS, Li S, Eastman JA: Measuring thermal conductivity of fluids containing oxide nanoparticles.

    J Heat Transf Trans ASME 1999, 121:280. Publisher Full Text OpenURL

  4. Choi SUS, Zhang ZG, Yu W, Lockwood FE, Grulke EA: Anomalous thermal conductivity enhancement in nanotube suspensions.

    Appl Phys Lett 2001, 79:2252. Publisher Full Text OpenURL

  5. Eastman JA, Choi SUS, Li S, Yu W, Thompson LJ: Anomalously increased effective thermal conductivities of ethylene glycol-based nanofluids containing copper nanoparticles.

    Appl Phys Lett 2001, 78:718. Publisher Full Text OpenURL

  6. Choi SUS, Zhang Z, Keblinski P: Nanofluids. In Encyclopedia of Nanoscience and Nanotechnology. Volume 757. Edited by Nalwa H. New York: American Scientific Publishers; 2004. OpenURL

  7. Das S, Choi SUS, Yu W, Pradeep T: Nanofluids Science and Technology. Hoboken, NJ: Wiley; 2008. OpenURL

  8. Jang SP, Choi SUS: Role of brownian motion in the enhanced thermal conductivity of nanofluids.

    Appl Phys Lett 2004, 84:4316. Publisher Full Text OpenURL

  9. Jang SP, Choi SUS: Effects of various parameters on nanofluid thermal conductivity.

    J Heat Transf Trans ASME 2007, 129:617. Publisher Full Text OpenURL

  10. Vadasz JJ, Govender S, Vadasz P: Heat transfer enhancement in nano-fluids suspensions: Possible mechanisms and explanations.

    Int J Heat Mass Transf 2005, 48:2673. Publisher Full Text OpenURL

  11. Vadasz P: Heat conduction in nanofluid suspensions.

    J Heat Transf Trans ASME 2006, 128:465. Publisher Full Text OpenURL

  12. Wu C, Cho TJ, Xu J, Lee D, Yang B, Zachariah MR: Effect of nanoparticle clustering on the effective thermal conductivity of concentrated silica colloids.

    Phys Rev E 2010, 81:011406. Publisher Full Text OpenURL

  13. Bai C, Wang L: Constructal design of particle volume fraction in nanofluids.

    J Heat Transf Trans ASME 2009, 131:112402. Publisher Full Text OpenURL

  14. Bai C, Wang L: Constructal allocation of nanoparticles in nanofluids.

    J Heat Transf Trans ASME 2010, 132:052404. Publisher Full Text OpenURL

  15. Fan J, Wang L: Constructal design of nanofluids.

    Int J Heat Mass Transf 2010, 53:4238. Publisher Full Text OpenURL

  16. Bejan A, Lorente S: Constructal theory of generation of configuration in nature and engineering.

    J Appl Phys 2006, 100:041301. Publisher Full Text OpenURL

  17. Bejan A, Lorente S: Design with Constructal Theory. Hoboken, NJ: Wiley; 2008. OpenURL

  18. Bejan A, Lorente S:

    Philos Trans Roy Soc B Biol Sci. 2010, 365:1335. Publisher Full Text OpenURL

  19. Bello-Ochende T, Meyer JP, Bejan A: Constructal multi-scale pin-fins.

    Int J Heat Mass Transf 2010, 53:2773. Publisher Full Text OpenURL

  20. Wu X, Wu H, Cheng P: Pressure drop and heat transfer of Al2O3-H2O nanofluids through silicon microchannels.

    J Micromech Microeng 2009, 19:105020. Publisher Full Text OpenURL

  21. Do KH, Jang SP: Effect of nanofluids on the thermal performance of a flat micro heat pipe with a rectangular grooved wick.

    Int J Heat Mass Transf 2010, 53:2183. Publisher Full Text OpenURL

  22. Ebrahimi S, Sabbaghzadeh J, Lajevardi M, Hadi I: Cooling performance of a microchannel heat sink with nanofluids containing cylindrical nanoparticles (carbon nanotubes).

    Heat Mass Transf 2010, 46:549. Publisher Full Text OpenURL

  23. Fan X, Chen H, Ding Y, Plucinski PK, Lapkin AA: Potential of 'nanofluids' to further intensify microreactors.

    Green Chem 10:670.


    Publisher Full Text OpenURL

  24. Li H, Liu S, Dai Z, Bao J, Yang Z: Applications of nanomaterials in electrochemical enzyme biosensors.

    Sensors 2009, 9:8547. Publisher Full Text OpenURL

  25. Munir A, Wang J, Zhou HS: Dynamics of capturing process of multiple magnetic nanoparticles in a flow through microfluidic bioseparation system.

    IET Nanobiotechnol 2009, 3:55. PubMed Abstract | Publisher Full Text OpenURL

  26. Huh D, Matthews BD, Mammoto A, Montoya-Zavala M, Hsin HY, Ingber DE: Reconstituting organ-level lung functions on a chip.

    Science 2010, 328:1662. PubMed Abstract | Publisher Full Text OpenURL

  27. Wang L, Fan J: Nanofluids research: Key issues.

    Nanoscale Res Lett 2010, 5:1241. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Sokolov A, Goldstein RE, Feldchtein FI, Aranson IS: Enhanced mixing and spatial instability in concentrated bacterial suspensions.

    Phys Rev E 2009, 80:031903. Publisher Full Text OpenURL

  29. Tsai T, Liou D, Kuo L, Chen P: Rapid mixing between ferro-nanofluid and water in a semi-active Y-type micromixer.

    Sensors Actuators A Phys 2009, 153:267. Publisher Full Text OpenURL

  30. Shitanda I, Yoshida Y, Tatsuma T: Microimaging of algal bioconvection by scanning electrochemical microscopy.

    Anal Chem 2007, 79:4237. PubMed Abstract | Publisher Full Text OpenURL

  31. Pedley TJ: Instability of uniform micro-organism suspensions revisited.

    J Fluid Mech 2010, 647:335. Publisher Full Text OpenURL

  32. Pedley TJ, Hill NA, Kessler JO: The growth of bioconvection patterns in a uniform suspension of gyrotactic microorganisms.

    J Fluid Mech 1988, 195:223. PubMed Abstract | Publisher Full Text OpenURL

  33. Kuznetsov AV: Thermo-bioconvection in a suspension of oxytactic bacteria.

    Int Commun Heat Mass Transf 2005, 32:991. Publisher Full Text OpenURL

  34. Kuznetsov AV: Investigation of the onset of thermo-bioconvection in a suspension of oxytactic microorganisms in a shallow fluid layer heated from below.

    Theor Comput Dyn 2005, 19:287. Publisher Full Text OpenURL

  35. Kuznetsov AV: The onset of thermo-bioconvection in a shallow fluid saturated porous layer heated from below in a suspension of oxytactic microorganisms.

    Eur J Mech B Fluids 2006, 25:223. Publisher Full Text OpenURL

  36. Avramenko AA, Kuznetsov AV: Bio-thermal convection caused by combined effects of swimming of oxytactic bacteria and inclined temperature gradient in a shallow fluid layer.

    Int J Numer Methods Heat Fluid Flow 2010, 20:157. Publisher Full Text OpenURL

  37. Kuznetsov AV, Avramenko AV: Effect of small particles on the stability of bioconvection in a suspension of gyrotactic microorganisms in a layer of finite depth.

    Int Commun Heat Mass Transf 2004, 31:1. Publisher Full Text OpenURL

  38. Geng P, Kuznetsov AV: Effect of small solid particles on the development of bioconvection plumes.

    Int Commun Heat Mass Transf 2004, 31:629. Publisher Full Text OpenURL

  39. Geng P, Kuznetsov AV: Settling of bidispersed small solid particles in a dilute suspension containing gyrotactic micro-organisms.

    Int J Eng Sci 2005, 43:992. Publisher Full Text OpenURL

  40. Kuznetsov AV, Geng P: The interaction of bioconvection caused by gyrotactic micro-organisms and settling of small solid particles.

    Int J Numer Methods Heat Fluid Flow 2005, 15:328. Publisher Full Text OpenURL

  41. Geng P, Kuznetsov AV: Introducing the concept of effective diffusivity to evaluate the effect of bioconvection on small solid particles.

    Int J Transp Phenom 2005, 7:321. OpenURL

  42. Kuznetsov AV: Non-oscillatory and oscillatory nanofluid bio-thermal convection in a horizontal layer of finite depth.

    Eur J Mech B Fluids 2011, 30(2):156-165. Publisher Full Text OpenURL

  43. Buongiorno J: Convective transport in nanofluids.

    J Heat Transf Trans ASME 2006, 128:240. Publisher Full Text OpenURL

  44. Hillesdon AJ, Pedley TJ, Kessler JO: The development of concentration gradients in a suspension of chemotactic bacteria.

    Bull Math Biol 1995, 57:299. PubMed Abstract OpenURL

  45. Hillesdon AJ, Pedley TJ: Bioconvection in suspensions of oxytactic bacteria: Linear theory.

    J Fluid Mech 1996, 324:223. Publisher Full Text OpenURL

  46. Nield DA, Kuznetsov AV: The onset of convection in a horizontal nanofluid layer of finite depth.

    Eur J Mech B Fluids 2010, 217:052405. OpenURL

  47. Finlayson BA: The Method of Weighted Residuals and Variational Principles. New York: Academic Press; 1972. OpenURL

  48. Chandrasekhar S: Hydrodynamic and Hydromagnetic Stability. Oxford: Clarendon Press; 1961. OpenURL