Steady Radial Diverging Flow of a Particle-Laden Fluid with Particle Migration

1. Introduction

Plane radial flow of an incompressible viscous fluid in a diverging channel with two non-parallel straight walls, called “Jeffery-Hamel (JH) flow” [1,2,3,4,5], remains an active research topic with significant practical application [6,7,8,9,10,11,12]. Remarkably, the exact solution of Navier-Stokes equations in this case admits the so-called “similarity solution” and the problem is reduced to a simpler 2nd-order nonlinear ordinary differential equation with constant coefficients. Recently, the research on JH channel flow has been extended to nanofluids with dispersed nanometer particles. However, to the best of our knowledge, almost all related works on JH channel flow of nanofluids (see e.g. [13,14,15]) have adopted the single-phase model [16,17,18] which assumes that the dispersed particles and the carrier fluid share the same velocity field and therefore cannot explain many important multiphase flow phenomena such as particle migration. In spite of active researches on various JH-like radial flows in a diverging/converging channel (see e.g. [19,20,21]), the steady spatial distribution of non-neutrally buoyant particles in radial flow of a viscous or inviscid particle-laden suspension is rarely studied in the literature.

In an attempt to study multiphase particulate radial flow in a diverging channel, it turns out that the governing equations of the multiphase model for particle-laden viscous fluids (with a larger number of equations than the Navier-Stokes equations for a clear fluid without dispersed particles) do not admit a similarity solution of JH type due to the required no-slip wall conditions. As a matter of fact, the similarity solution of JH type is not admitted even in the diverging pipe flow of a clear fluid without particles [22].

Therefore, the present work will focus on the JH-like plane radial flow of a particle-laden suspension when the no-slip wall conditions are not applied (see the two types of problems shown in Figure 1 below), with specific interest in the steady particle distribution of a particle-laden suspension attained eventually as a result of particle migration.

The general equations of the two-fluid model for the present problem are given in section 2. In section 3, the initial particle velocity field of a particle-laden suspension with uniform particle distribution is studied with an emphasis on the velocity shift between the particles and the suspension. The steady particle distribution of plane radial flow due to particle migration is studied in section 4 with demonstrated results for various values of the Stokes number of particles heavier or lighter than the carrier fluid, and an asymptotic expression is given for the steady particle distribution for the vanishingly small particle Stokes number. Finally, the main results are summarized in section 5.

2. Equations of the Present Model

Let us consider an incompressible (viscous or inviscid) fluid with initially uniform distribution of identical rigid spherical particles of radius r_S, as an incompressible particle-fluid suspension.

2.1. General Equations of the Model with Particle Migration

With the present model, hydrodynamics of an incompressible suspension with dispersed solid spheres is governed by the modified form of Navier-Stokes equations (the gravity not involved)

$$\rho {\displaystyle \frac{d{\mathit{v}}_m}{dt}}=\rho \left[{\displaystyle \frac{\partial {\mathit{v}}_m}{\partial t}}+\left({\mathit{v}}_m\bullet \nabla \right){\mathit{v}}_m\right]=-\nabla p+\nabla \bullet \left(\mu \left[\nabla \mathit{v}+{\left(\nabla \mathit{v}\right)}^T\right]\right),$$

(1)

$$\mathrm{d}\mathrm{i}\mathrm{v}\mathit{v}=0,$$

(2)

$${\displaystyle \frac{\partial }{\partial t}}\rho +\mathrm{d}\mathrm{i}\mathrm{v}\left[\rho {\mathit{v}}_m\right]=0.$$

(3)

where x and t are the spatial coordinates and time, p(x, t) is the pressure field of the suspension, v(x, t) is the velocity field of the suspension (defined as the velocity field of the geometrical center of the representative unit cell of suspension), v_m(x, t) is the velocity field of the mass center of the representative unit cell defined by (A3) in Appendix, the effective density ρ (per unit volume) of the suspension is given by (A2) in Appendix, ρ_S and ρ_f are the densities of the particles and the carrier fluid, respectively. Here, δ is the volume fraction of particles, μ is the effective viscosity of the suspension which can be estimated by Einstein formula $\mu ={\mu }_f\left(1+2.5\delta \right)$ with the viscosity μ_f of the carrier fluid in the dilute limit, and $\nabla$ and ${\nabla }^2$ are the gradient and Laplacian operators. In general, if the particle volume fraction δ changes with the spatial position and time due to particle migration, the mass density ρ and the effective viscosity μ can depend on the spatial position and time.

As explained in Appendix, the Newton’s second and third laws imply that the resultant external force acting on the representative unit cell, given by the terms on right-hand side of (1), equates to the mass of the unit cell multiplied by the acceleration dv_m/dt of its mass center (rather than the acceleration field dv/dt of its geometrical center), which leads to the above modified form of Navier-Stokes equations (1). Clearly, for a homogenous clear fluid (δ=0), v_m(x, t)=v(x, t) and the equation (1) reduces to the classical form of Navier-Stokes equations.

For a suspension with non-neutrally buoyant particles (ρ_S≠ρ_f), we have v_m(x, t)≠v(x, t). As explained in Appendix, with the gravity not involved, an additional relationship between v_m(x, t) and v(x, t) is given as

$${\mathit{v}}_m+a{\displaystyle \frac{d{\mathit{v}}_m}{dt}}+C_L{\displaystyle \frac{2{\rho r}_S^2}{9\mu }}\left[{\mathit{v}}_m\times \left(\nabla \times \mathit{v}\right)\right]=\mathit{v}+b{\displaystyle \frac{d\mathit{v}}{dt}}+C_L{\displaystyle \frac{2{\rho r}_S^2}{9\mu }}\left[\mathit{v}\times \left(\nabla \times \mathit{v}\right)\right],$$

(4)

$$a=\left(1+C_a{\displaystyle \frac{\rho }{{\rho }_S}}\right){\displaystyle \frac{2{\rho }_S{r}_S^2}{9\mu }},b=a\left({\displaystyle \frac{{\rho }_f}{\rho }}+{\displaystyle \frac{\left(1+C_a\right)\left({\rho }_S-{\rho }_f\right)\delta }{\left(C_a\rho +{\rho }_S\right)}}\right).$$

(5)

Here, d/dt denotes the material derivative of the associated velocity field along its own streamlines, C_a and C_L are the added mass and the lift force coefficients (C_a=C_L=0.5 is often adopted in literature), respectively, the coefficients a and b are derived by considering the Stokes drag, the forces acting on particles due to added mass and fluid acceleration [23,24], and the lift force [25,26,27,28], although the lift force vanishes for the present problem of r-dependent radial flow with $\left(\nabla \times \mathit{v}=0\right)$. Here, it should be stated that the Stokes drag-based models of particle-laden inviscid fluids have been widely adopted for inviscid particulate flows [29,30,31,32,33]. For instance, the inviscid version of Stokes drag-based Saffman model [29] was used by Michael [30] to study Kelvin-Helmholtz instability of particle-laden inviscid flows.

It is stated that the second terms inside the brackets in the expressions of a and b in (5) will be absent (then $a={\displaystyle \frac{2{\rho }_S{r}_S^2}{9\mu }},b=\left({\displaystyle \frac{{\rho }_f}{\rho }}\right)a$) if only the Stokes drag is considered, and a=b and v_m(x,t)=v(x,t) when either δ=0 or ρ_S=ρ_f (“neutrally buoyant particles”) and then the present model reduces to the single-phase models [16,17,18] (see (A1) in Appendix).

2.2. Equations for Steady Plane Radial Flow

It can be verified that, unlike a clear viscous fluid (without dispersed particles) which admits the similarity solution of JH-type for radial flow in a diverging channel, the particle-laden viscous suspensions do not admit such a similarity solution for radial flow in a diverging channel. Therefore, in the present paper we shall focus on the following two problems of steady plane radial flow of an incompressible viscous or inviscid particle-laden fluid shown in Figure 1 whose flow fields depend solely on the radial coordinate r:

(a) axisymmetric plane viscous radial flow from a point source;

(b) inviscid radial flow in a diverging channel.

To our knowledge, radial flow of a particle-laden (viscous or inviscid) suspension shown in Figure 1 has been rarely studied in literature, although radial flow of a clear fluid (without dispersed particles) from a point source has been the topic of several known works [34,35,36,37]. Clearly, the case (a) is an r-dependent axisymmetric flow, and the inviscid radial flow in a diverging channel shown in the case (b) depends solely on the radial coordinate r because the inviscid flow is free to slip on the straight walls.

The steady r-dependent radial flow field (u(r), u_m(r), p(r), δ(r)) in the cylindrical coordinate (r, θ, z) system are given by

$$\mathit{v}=\left(u\left(r\right),0,0\right),{\mathit{v}}_{\mathit{m}}=\left(u_m\left(r\right),0,0\right),p\left(r\right),\delta \left(r\right).$$

(6)

Note that

$$\nabla \mathit{v}{+\left(\nabla \mathit{v}\right)}^{\mathrm{T}}=2\left[\begin{array}{ccc}{u}_{,r}& 0& 0\\ 0& {\displaystyle \frac{u}{r}}& 0\\ 0& 0& 0\end{array}\right],\nabla \bullet \left(\mu \left(r\right)\left[\nabla \mathit{v}+{\left(\nabla \mathit{v}\right)}^T\right]\right)=2\left[\begin{array}{c}\mu \left(r\right){\displaystyle \frac{\partial }{\partial r}}\left({\displaystyle \frac{1}{r}}{\displaystyle \frac{\partial }{\partial r}}\left(ru\right)\right)+u_{,r}{\displaystyle \frac{\partial \mu }{\partial r}}\\ 0\\ 0\end{array}\right].$$

(7)

where the subscript “, “ denotes the partial derivative with respect to r. Thus, when the particle migration is involved, it is verified that equations (1-4) give the following 4 nonlinear equations for (u, u_m, p, δ) as the functions of the single variable r

$$u_m{\displaystyle \frac{\partial u_m}{\partial r}}=-{\displaystyle \frac{1}{\rho \left(r\right)}}{\displaystyle \frac{\partial p}{\partial r}}+{\displaystyle \frac{2}{\rho \left(r\right)}}{u}_{,r}{\displaystyle \frac{\partial \mu }{\partial r}},$$

(8)

$${\displaystyle \frac{\partial }{\partial r}}\left(ru\right)=0,$$

(9)

$$u_m{\displaystyle \frac{\partial \rho }{\partial r}}+\rho \left(r\right){{\displaystyle \frac{1}{r}}\left(r{u}_m\right)}_{,r}=0.$$

(10)

$$u_m+a\left(u_m{\displaystyle \frac{\partial u_m}{\partial r}}\right)=u+b\left(u{\displaystyle \frac{\partial u}{\partial r}}\right),$$

(11)

Let us assume that the particles and fluid have two independent inlet velocities (u_0S, u_0f) with the constant inlet particle fraction δ₀ at the entrance r=r₀ shown in Figure 1, thus the inlet values of (u_m, u) are given by

$${\left.u\right|}_{r=r_0}=u_0={\delta }_0{u}_{0S}+\left({1-\delta }_0\right)u_{0f},{{\left.u_m\right|}_{r=r_0}=\left(u_m\right)}_0={\displaystyle \frac{{{\rho }_S\delta }_0{u}_{0S}+{\rho }_f\left({1-\delta }_0\right)u_{0f}}{{\left({\rho }_S-{\rho }_f\right)\delta }_0+{\rho }_f}}.$$

(12)

For the conciseness of mathematical analysis, let us confine ourselves to case when the inlet velocity of fluid is not zero (u_0f>0). The results derived could offer a qualitative understanding of the limiting case (u_0f=0) by considering sufficiently small inlet velocity of the fluid.

It follows from (9) that $u\left(r\right)={\displaystyle \frac{f}{r}}$, where f is a constant. In the present work, with $u\left(r\right)={\displaystyle \frac{f}{r}}$, we shall focus on equations (10, 11) for the steady velocity field u_m(r) and the particle distribution δ(r), and the pressure p(r) can be determined from (8) once u_m(r) and δ(r) are known.

3. Initial Velocity Field with the Uniform Particle Distribution

In this section, to illustrate why the particulate radial flow with initially uniform particle distribution leads to particle migration, let us first study the initial particle velocity field of a particle-laden viscous suspension with uniformly distributed particles under the assumption that the particle migration is slow enough so that the initial flow field is nearly steady with the constant particle volume fraction and equation (10) associated with particle migration can be ignored.

Thus, the parameters (δ=δ₀, ρ, μ) are all constants in this section, and the unknown radial velocity u_m(r) is determined by (11) which gives

$$u_m+a\left(u_m{\displaystyle \frac{\partial u_m}{\partial r}}\right)={\displaystyle \frac{f}{r}}-b{\displaystyle \frac{1}{r^3}}{f}^2,u\left(r\right)={\displaystyle \frac{f}{r}}.$$

(13)

For a clear fluid (δ₀=0) without dispersed particles, it follows from the definition (A3) that v_m(x, t)=v(x, t), and we have a=b and

$${\delta }_0=0:u_m\left(r\right)=u\left(r\right)=u_0\left(r\right)={\displaystyle \frac{f_0}{r}},f_0=r_0{u}_{0f}>0.$$

(14)

For the particle-laden fluid (δ₀>0), let us write the particle-disturbed flow field u_m(r) as

$${\delta }_0>0:u_m\left(r\right)=u\left(r\right)+∆\left(r\right),f=f_0+r_0{\delta }_0\left(u_{0S}-u_{0f}\right).$$

(15)

Here, let us focus on the dilute limit of particle-laden fluids when the volume fraction of particles is much smaller than the unity. With the Einstein formula $\mu ={\mu }_f\left(1+\alpha \delta \right)$, up to the first powers of δ, the δ–dependent coefficients (μ, ρ, a, b) are expanded as

$$\rho ={\rho }_f\left[1+{\displaystyle \frac{\left({\rho }_S-{\rho }_f\right)}{{\rho }_f}}\delta \right],{\displaystyle \frac{\mu }{{\mu }_f}}=1+\alpha \delta ,$$

$$a=a_0\left(1-\left[\alpha -{\displaystyle \frac{C_a\left({\rho }_S-{\rho }_f\right)}{{C_a\rho }_f+{\rho }_S}}\right]\delta \right),a_0=\left(1+C_a{\displaystyle \frac{{\rho }_f}{{\rho }_S}}\right){\displaystyle \frac{2{\rho }_S{r}_S^2}{9{\mu }_f}},$$

$$b=a_0\left[1-\left(\left[\alpha -{\displaystyle \frac{C_a\left({\rho }_S-{\rho }_f\right)}{{C_a\rho }_f+{\rho }_S}}\right]+{\displaystyle \frac{{\left({\rho }_S-{\rho }_f\right)}^2}{{\rho }_f\left(C_a{\rho }_f+{\rho }_S\right)}}\right)\delta \right].$$

(16)

where $\left({\displaystyle \frac{2{\rho }_S{r}_S^2}{9{\mu }_f}}\right)$ is the dimensional relaxation time of suspended particles, and $a_0=C_a{\displaystyle \frac{2{\rho }_f{r}_S^2}{9{\mu }_f}}$ is the modified relaxation time for massless gas bubbles due to the added mass. The coefficient α depends on the nature of dispersed particles. Typically, α=2.5 is commonly adopted for rigid spheres, and α=1 is suggested for spherical gas bubbles in a liquid [38,39].

For the dilute particle-fluid suspensions with the small number δ<<1, the disturbed velocity field Δ(r)=u_m(r)-u(r) scales with the number δ, and therefore Δ(r) is of the order δ. Substituting (14, 15) into (12) and ignoring all nonlinear terms of Δ(r) and δ, the linear equation for Δ(r) gives

$${\displaystyle \frac{\partial ∆}{\partial r}}+\left({\displaystyle \frac{r^2-a_0{f}_0}{a_0{f}_{0}r}}\right)∆={\displaystyle \frac{{\delta }_0}{r^2}}\left[{\displaystyle \frac{f_0{\left({\rho }_S-{\rho }_f\right)}^2}{{\rho }_f\left(C_a{\rho }_f+{\rho }_S\right)}}\right],$$

(17)

with the boundary condition at the entrance

$${\left.∆\right|}_{r=r_0}={\left(u_m\right)}_0-u_0={\displaystyle \frac{{\left({\rho }_S-{\rho }_f\right)\delta }_0\left({1-\delta }_0\right)\left(u_{0S}-u_{0f}\right)}{{\left({\rho }_S-{\rho }_f\right)\delta }_0+{\rho }_f}}.$$

(18)

The homogeneous solution of (17) is of the form $\left(r{e}^{\left({\displaystyle \frac{-r^2}{{2a}_0{f}_0}}\right)}\right)$. Here, in the case f₀≠0, using the method of variation of constant for (17), explicit solution of Δ(r) for (17) is given by

$$∆\left(r\right)=\left[\begin{array}{c}{\displaystyle \frac{{\left({\rho }_S-{\rho }_f\right)\delta }_0\left({1-\delta }_0\right)\left(u_{0S}-u_{0f}\right)}{\left[{\left({\rho }_S-{\rho }_f\right)\delta }_0+{\rho }_f\right]r_0}}{e}^{\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)}\\ +{\delta }_0\left[{\displaystyle \frac{f_0{\left({\rho }_S-{\rho }_f\right)}^2}{{\rho }_f\left(C_a{\rho }_f+{\rho }_S\right)}}\right]{\int }_{r_0}^r{\displaystyle \frac{e^{\left(\frac{t^2}{{2a}_0{f}_0}\right)}}{t^3}}dt\end{array}\right]r{e}^{\left({\displaystyle \frac{-r^2}{{2a}_0{f}_0}}\right)}.$$

(19)

We are particularly interested in the relative velocity of the particles with respect to the fluid. Based on the general relation (A4), we have

$$\left(u_S-u\right)={\displaystyle \frac{\rho \left(u_m-u\right)}{\delta \left({\rho }_S-{\rho }_f\right)}},\left(u_S-u_f\right)={\displaystyle \frac{\rho \left(u_m-u\right)}{\delta \left(1-\delta \right)\left({\rho }_S-{\rho }_f\right)}}.$$

(20)

Thus, up to the lowest order of δ, the velocity difference between the particles and the fluid normalized by the volume-averaged velocity of the suspension is given by

$${\displaystyle \frac{\left(u_S-u_f\right)}{u}}=\left[\begin{array}{c}{\displaystyle \frac{\rho \left(u_{0S}-u_{0f}\right)r_0}{\left[{\left({\rho }_S-{\rho }_f\right)\delta }_0+{\rho }_f\right]f}}{e}^{\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)}\\ +\left[{\displaystyle \frac{\rho \left({\rho }_S-{\rho }_f\right)f_0}{{\rho }_f\left(C_a{\rho }_f+{\rho }_S\right)f}}\right]{\int }_1^{r^{*}}{\displaystyle \frac{e^{\left(\frac{r_0^2}{{2a}_0{f}_0}{t}^2\right)}}{t^3}}dt\end{array}\right]r^{*2}{e}^{\left({\displaystyle \frac{-r_0^2}{{2a}_0{f}_0}}\right)r^{*2}},r^{*}={\displaystyle \frac{r}{r_0}}\ge 1.$$

(21)

It is seen from (21) that the velocity shift between particles and the fluid vanishes for the neutrally buoyant particles (ρ_S=ρ_f), this is consistent with the fact that the lift force [25,26,27,28] responsible for migration of the neutrally buoyant particles vanishes for the r-dependent radial flow with $\left(\nabla \times \mathit{v}=0\right)$.

3.1. Lighter Particles with Higher Inlet Velocity (u_0S>u_0f)

Let us first discuss the case when the radial flow is driven by the high-speed injection of lighter particle with (u_0S>u_0f). This problem is of major interest in the literature on bubble-driven gas-liquid two-phase flow [40,41,42,43,44,45,46]. It should be stated that the physical concepts and mathematical equations formulated in two-fluid models for dispersed solid spherical particles can be largely applied to fluids with dispersed small spherical gas bubbles when the effects of deformation and non-uniform size distribution of gas bubbles can be ignored [47,48,49,50].

For massless bubbles with (ρ_S<<ρ_f) and Ca=0.5 [47,48,49,50], it follows from (21) that

$${\displaystyle \frac{\left(u_S-u_f\right)}{u}}=\left[{\displaystyle \frac{\left(u_{0S}-u_{0f}\right)}{u_{0f}}}{e}^{\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)}-2{\int }_1^{r^{*}}{\displaystyle \frac{e^{\left(\frac{r_0^2}{{2a}_0{f}_0}{t}^2\right)}}{t^3}}dt\right]r^{*2}{e}^{\left({\displaystyle \frac{-r_0^2}{{2a}_0{f}_0}}\right)r^{*2}},r^{*}={\displaystyle \frac{r}{r_0}}\ge 1.$$

(22)

Here, let us consider two cases when the inlet velocity of bubbles is moderately or much higher than the inlet velocity of the fluid at the entrance of the flow, with u_0S=2u_0f and u_0S=10u_0f, respectively. The normalized velocity difference ((u_S-u_f)/u) between the massless bubbles and the fluid given by (22) are shown in Figure 2 and Figure 3 for several typical values of the ratio $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$ which is considered to be inversely proportional to the (modified) Stokes number of bubbles with $a_0=C_a{\displaystyle \frac{2{\rho }_f{r}_S^2}{9{\mu }_f}}$ , where u_S , u_f and u are the velocities of the bubbles, the fluid and the suspension, respectively.

It is seen from Figure 2 and Figure 3 that although the bubbles move faster than the fluid within a finite distance from the entrance (r*=r/r₀=1), the velocity of bubbles becomes lower than the velocity of fluid beyond that distance. This distance is determined by the bubble-to-fluid inlet velocity ratio and the Stokes number of bubbles. For example, it is seen from Figure 2 for u_0S=2u_0f that the velocity of bubbles of a moderate Stokes number presented by $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$=1 becomes lower than the velocity of fluid beyond a distance slightly above (r*=2), while the velocity of bubbles of a larger Stokes number presented by $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$=0.01 remains faster than the fluid within the distance above (r*=6).

For the higher bubble-to-fluid inlet velocity ratio u_0S=10u_0f, it is seen from Figure 3 that the velocity of bubbles of a moderate Stokes number presented by $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$=1 remains faster than the fluid within the distance above (r*=3), while the velocity of bubbles of a larger Stokes number presented by $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$=0.01 remains faster than the fluid within the large distance abound (r*=33). In particular, it is seen from Figure 3 that the normalized velocity difference given by (22) in the case $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$=0.01 can be very high (beyond 100, too high to be shown entirely in Figure 3) because the suspension velocity diminishes quickly while the bubbles of large Stokes number respond to the decelerating flow field much slower than the suspension.

Since most experiments and numerical simulations on high speed bubble-driven decelerating gas-liquid flow are limited to a diverging channel or pipe of limited finite length [41,42,43] and the liquid is nearly stationary, our results shown in Figure 2 and Figure 3 could suggest that the bubbles of even a moderate Stokes number will move faster than the fluid within a sufficiently long distance when the bubble-to-fluid inlet velocity ratio is very high, consistent with some known experimental observations and numerical simulations.

In addition, it is worthy mention that when the inlet velocities (and therefore the coefficient f) change their signs simultaneously, the normalized velocity difference given by the right-hand side of (21) remains unchanged. This could suggest that the results derived here for a diverging channel may be qualitatively valid for a converging channel. Consistent with this, it is noted that the bubbles accelerate faster than the fluid in the accelerating bubbly flow in a converging channel of finite length, as reported by Auton et al. [40] and reviewed by Magnaudet & Eames [47].

3.2. Particles and Fluid Have the Same Inlet Velocity

When the particles and the fluid get into the decelerating flow field (du/dr<0) in a diverging radial flow with the same inlet velocity (u_0S=u_0f) at the entrance, because the Stokes drag vanishes there, the heavier (lighter) particles of larger (less) inertia will respond to the decelerating flow field slower (faster) than the fluid. Therefore, the particles heavier (lighter) than the fluid will move faster (slower) than the fluid at least within a certain distance from the entrance. In particular, for the massless bubbles with C_a=0.5, it is readily seen from (A7) of Appendix that the deceleration of the bubbles is 3 times the deceleration of the fluid at the entrance when the bubbles and the fluid have the same inlet velocity and the Stokes drag vanishes at the entrance, which implies that the velocity of bubbles is slower than the velocity of fluid at least within a certain distance from the entrance.

Beyond the entrance, the motion of bubbles is determined by the three terms on the right-hand side of (A7) of Appendix, although the motion of heavy particles (ρ_S>>ρ_f) is dominated by the Stokes drag and the other two terms on the right-hand side of (A7) can be ignored. Actually, with the condition (u_0S=u_0f), it follows from (19) that

$$u_{inf}=u_{inS}:∆\left(r\right)={\delta }_0\left[{\displaystyle \frac{f_0{\left({\rho }_S-{\rho }_f\right)}^2}{{\rho }_f\left(C_a{\rho }_f+{\rho }_S\right)}}\right]\left[{\int }_{r_0}^r{\displaystyle \frac{e^{\left(\frac{t^2}{{2a}_0{f}_0}\right)}}{t^3}}dt\right]r{e}^{\left({\displaystyle \frac{-r^2}{{2a}_0{f}_0}}\right)}.$$

(23)

It follows from (21) that

$${\displaystyle \frac{\left(u_S-u_f\right)}{u}}={\displaystyle \frac{\rho \left({\rho }_S-{\rho }_f\right)}{{\rho }_f\left(1-{\delta }_0\right)\left(C_a{\rho }_f+{\rho }_S\right)}}\Omega \left(r^{*}\right),\Omega \left(r^{*}\right)\equiv r^{*2}{e}^{\left({\displaystyle \frac{-r_0^2}{{2a}_0{f}_0}}\right)r^{*2}}\left[{\int }_1^{r^{*}}{\displaystyle \frac{e^{\left(\frac{r_0^2}{{2a}_0{f}_0}\right)t^2}}{t^3}}dt\right],r^{*}={\displaystyle \frac{r}{r_0}}\ge 1.$$

(24)

With the dimension r₀ of the present problem and the inlet velocity ${\displaystyle \frac{f_0}{r_0}}$ at the entrance (r=r₀), the ratio $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$ is inversely proportional to the Stokes number of particles.

The dimensionless function Ω(r*) defined in (24), which determines the velocity difference between the particles and the fluid, is plotted in Figure 4 for three different values of $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$=1, 3 and 10, respectively. It is seen from Figure 4 and (24) that the particles heavier (or lighter) than the fluid will move faster (or slower) than the fluid velocity. In addition, because the ratio $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$ is inversely proportional to the Stokes number of particles, it is seen from (24) that the velocity difference between the particles and the fluid decreases with increasing value of $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$ (or decreasing Stokes number of particles), which suggests that the particle migration can be slow for nanofluids of nanometer particles although the long-term particle migration of nanofluids can be relevant and cannot be ignored for a sufficiently long period of time.

It should be stated that all results derived in section 3 are based on the assumed steady initial flow field with constant parameters (δ, ρ, μ), which only serve to explain the particle migration at the initial stage of the radial flow but cannot offer any detailed data on the steady flow attained eventually as a result of particle migration beyond the initial stage of flow.

4. Steady Particle Distribution of Plane Radial Diverging Flow

Since the radial flow of a particulate fluid with initially uniform distribution of particles cannot remain the uniform particle distribution due to particle migration, it is of major interest to study the steady particle distribution attained eventually as a result of long-term particle migration.

Here, to determine the steady volume fraction δ(r) of particles, substituting (15, 16) into (10, 11) and ignoring all nonlinear terms of Δ(r) and δ(r), the linear equations for Δ(r) and $\delta \left(r\right)$ give

$$\left({\rho }_S-{\rho }_f\right){{\displaystyle \frac{f_0}{r}}{\displaystyle \frac{d\delta }{dr}}+{\displaystyle \frac{{\rho }_f}{r}}\left(r\Delta \right)}_{,r}=0,$$

(25)

$${\displaystyle \frac{d\mathsf{\Delta }}{dr}}+\left({\displaystyle \frac{r^2-a_0{f}_0}{a_0{f}_{0}r}}\right)\mathsf{\Delta }=\left[{\displaystyle \frac{f_0{\left({\rho }_S-{\rho }_f\right)}^2}{{\rho }_f\left(C_a{\rho }_f+{\rho }_S\right)}}\delta \left(r\right)\right]{\displaystyle \frac{1}{r^2}},$$

(26)

with the boundary conditions

$${\left.{\left.∆\right|}_{r=r_0}={\displaystyle \frac{{\left({\rho }_S-{\rho }_f\right)\delta }_0\left({1-\delta }_0\right)\left(u_{0S}-u_{0f}\right)}{{\left({\rho }_S-{\rho }_f\right)\delta }_0+{\rho }_f}},\delta \right|}_{r=r_0}={\delta }_0.$$

(27)

It follows from (25) and the conditions (27) that

$$\left[{\displaystyle \frac{{\left({\rho }_S-{\rho }_f\right)\delta }_0\left({1-\delta }_0\right)\left(u_{0S}-u_{0f}\right)}{{\left({\rho }_S-{\rho }_f\right)\delta }_0+{\rho }_f}}{r}_0-r\Delta \left(r\right)\right]=\left({\rho }_S-{\rho }_f\right){\displaystyle \frac{f_0}{{\rho }_f}}\left(\delta \left(r\right)-{\delta }_{in}\right).$$

(28)

Using (28) to eliminate Δ(r) in (26), the following first-order linear equation for δ(r) can be verified

$${\displaystyle \frac{d\delta }{dr}}+\left[{\displaystyle \frac{r}{a_0{f}_0}}-{\displaystyle \frac{\left({\rho }_S+\left(2C_a+1\right){\rho }_f\right)}{\left(C_a{\rho }_f+{\rho }_S\right)r}}\right]\delta \left(r\right)=B{\displaystyle \frac{{\delta }_{0}r}{a_0{f}_0}}\left(1-2{\displaystyle \frac{a_0{f}_0}{r^2}}\right),$$

(29)

With the constant B is given by

$$B=\left[1+{\displaystyle \frac{\left({1-\delta }_0\right)\left(u_{0S}-u_{0f}\right){\rho }_f}{\left[{\left({\rho }_S-{\rho }_f\right)\delta }_0+{\rho }_f\right]u_{0f}}}\right].$$

(30)

The homogeneous solution of (29) is of the form

$$\delta \left(r\right)\propto r^{\left[{\displaystyle \frac{{\rho }_S+\left(2C_a+1\right){\rho }_f}{\left(C_a{\rho }_f+{\rho }_S\right)}}\right]}{e}^{{\displaystyle \frac{{-r}^2}{{2a}_0{f}_0}}}.$$

(31)

On using the variation of constant, explicit solution of the non-homogeneous equation (29) with the boundary condition (27) gives

$$\begin{gathered} {\displaystyle \frac{\delta \left(r^{*}\right)}{{\delta }_0}}=\left[e^{{\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}}+B{\displaystyle \frac{r_0^2}{a_0{f}_0}}{\int }_1^{r^{*}}\left(1-2{\displaystyle \frac{a_0{f}_0}{t^2{r}_0^2}}\right)t^{\left[1-\left({\displaystyle \frac{{\rho }_S+\left(2C_a+1\right){\rho }_f}{\left(C_a{\rho }_f+{\rho }_S\right)}}\right)\right]}{e}^{\left({\displaystyle \frac{r_0^2{t}^2}{{2a}_0{f}_0}}\right)}dt\right] \\ \times {r^{*}}^{\left[{\displaystyle \frac{{\rho }_S+\left(2C_a+1\right){\rho }_f}{\left(C_a{\rho }_f+{\rho }_S\right)}}\right]}{e}^{\left({\displaystyle \frac{-r_0^2}{{2a}_0{f}_0}}\right)r^{*2}},r^{*}\equiv {\displaystyle \frac{r}{r_0}}\ge 1. \end{gathered}$$

(32)

As expected, it can be verified from (32) that δ(r^*)/δ₀ ≡1 for the neutrally buoyant particles (ρ_S=ρ_f).

4.1. Light Particles with Higher Inlet Velocity

Let us first discuss the two cases discussed in section 3.1 when the inlet velocity of bubbles are higher than the inlet velocity of the fluid. In these case with Ca=0.5 for bubbles, up to the lowest order of δ, we have

$$\begin{gathered} {\displaystyle \frac{\delta \left(r^{*}\right)}{{\delta }_0}}=\left[e^{{\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}}+B{\displaystyle \frac{r_0^2}{a_0{f}_0}}{\int }_1^{r^{*}}\left(1-2{\displaystyle \frac{a_0{f}_0}{t^2{r}_0^2}}\right)t^{-3}{e}^{\left({\displaystyle \frac{r_0^2{t}^2}{{2a}_0{f}_0}}\right)}dt\right]{r^{*}}^4{e}^{\left({\displaystyle \frac{-r_0^2}{{2a}_0{f}_0}}\right)r^{*2}},r^{*}\equiv {\displaystyle \frac{r}{r_0}}\ge 1; \\ B={\displaystyle \frac{\left(u_{0S}-u_{0f}\right)}{u_{0f}}}. \end{gathered}$$

(33)

The dimensionless steady volume fraction of bubbles given by (33), with u_0S=2u_0f and u_0S=10u_0f, are shown in Figure 5 and Figure 6, respectively, for three values of $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$. It is seen from Figure 5 and Figure 6 that the volume fraction of bubbles attains its maximum at a location near the entrance of the flow and after then monotonically decreases with increasing radial coordinate and converges to a finite value determined by the inlet velocity ratio of the bubbles and the fluid, consistent with the conservation of bubbles without considering the breakup of bubbles. In addition, the maximum volume fraction and its location approach the inlet volume fraction multiplied by the inlet velocity ratio and the entrance location of the flow, respectively, as the Stokes number of bubbles approaches zero (or equivalently, as $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$ tends to infinity).

4.2. Particles and Fluid Have the Same Inlet Velocity

On the other hand, when the particles and the fluid have the same inlet velocity (u_0S=u_0f) with C_a=0.5, we have B=1 and

$$\begin{gathered} {\displaystyle \frac{\delta \left(r^{*}\right)}{{\delta }_0}}=\left[e^{{\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}}+{\displaystyle \frac{r_0^2}{a_0{f}_0}}{\int }_1^{r^{*}}\left(1-{\displaystyle \frac{{2a}_0{f}_0}{t^2{r}_0^2}}\right)t^{\left[1-{\displaystyle \frac{2\left({\rho }_S+2{\rho }_f\right)}{\left({\rho }_f+2{\rho }_S\right)}}\right]}{e}^{\left({\displaystyle \frac{r_0^2}{{2a}_{0}f}}\right)t^2}dt\right]{r^{*}}^{\left[{\displaystyle \frac{2\left({\rho }_S+2{\rho }_f\right)}{\left({\rho }_f+2{\rho }_S\right)}}\right]}{e}^{\left({\displaystyle \frac{-r_0^2}{{2a}_{0}f}}\right)r^{*2}}, \\ {r}^{*}\equiv {\displaystyle \frac{r}{r_0}}\ge 1. \end{gathered}$$

(34)

The dimensionless steady volume fraction of particles given by (34) is plotted in Figure 7 for heavy particles (${\displaystyle \frac{2\left({\rho }_S+2{\rho }_f\right)}{\left({\rho }_f+2{\rho }_S\right)}}\approx 1$) with three values of $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$=1, 3 and 10, respectively. It is seen from Figure 7 that the volume fraction of heavy particles attains its minimum at a location nearby the entrance of the flow, and after then the particle volume fraction monotonically increases with increasing radial coordinate and converges to a finite value. Particularly, the minimum value of the particle volume fraction and its location approach the inlet value δ₀ and the entrance location r=r₀ as the value of $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$ approaches infinity, or equivalently as the Stokes number of particles approaches zero.

On the other hand, the dimensionless steady volume fraction of particles given by (34) is plotted in Figure 8 for massless bubbles (${\displaystyle \frac{2\left({\rho }_S+2{\rho }_f\right)}{\left({\rho }_f+2{\rho }_S\right)}}\approx 4$) with three values of $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$=1, 3 and 10, respectively. It is seen from Figure 8 that the volume fraction of massless bubbles attains its maximum at a location nearby the entrance of the flow, and after then the bubble volume fraction monotonically decreases with increasing radial coordinate and converges to a finite value. Particularly, the maximum value of the particle volume fraction and its location approach the inlet particle volume fraction δ₀ and the entrance location r=r₀ as the value of $\left({\displaystyle \frac{r_0^2}{{2a}_0{f}_0}}\right)$ approaches infinity, or equivalently as the modified Stokes number of bubbles approaches zero.

To our current knowledge, the existence of the location with the minimum (or maximum) volume fraction of heavier (or lighter) particles in the steady particle distribution of a diverging radial flow has not been addressed in literature, and therefore a comparison of this interesting prediction with known data cannot be made here due to the lack of available related results in literature.

5. Conclusions

Steady spatial distribution of particles in various flow problems of particle-laden viscous or inviscid fluids is not extensively addressed in literature. The present work focuses on the diverging plane radial flow of a particle-fluid viscous or inviscid suspension when the velocity field depends on the radial coordinate solely, with particular interest in the steady spatial particle distribution attained eventually as a result of particle migration. Our main results include

(1)

In the initial flow field of a particle-fluid suspension with uniformly distributed particles, the relative velocity of particles with respect to the fluid depends on their inlet velocity ratio, the mass density ratio and the Stokes number of particles. For example, when their inlet velocities are equal (then Stokes drag vanishes at the entrance), the particles heavier (or lighter) than the fluid will move faster (or slower) than the fluid. On the other hand, the particles lighter than the fluid can remain faster than the fluid within a sufficiently long distance provided that the inlet velocity of lighter particles is much higher than the inlet velocity of the fluid. This result is qualitatively consistent with some known simulations and experiments on gas-liquid bubbly flow in a diverging channel of finite length driven by high-speed injection of gas bubbles into a nearly stationary liquid.

(2)

An explicit expression is obtained for the steady spatial distribution of particles attained eventually as a result of particle migration. In particular, for massless gas bubbles with the inlet velocity higher than the inlet velocity of the fluid, our results show that the volume fraction of bubbles attains its maximum at a location close to the entrance of the flow and after then monotonically decreases with increasing radial coordinate and converges to a finite value determined by the inlet velocity ratio of the bubbles and the fluid. In addition, the maximum volume fraction and its location approach the inlet volume fraction of the bubbles multiplied by the inlet velocity ratio and the entrance location of the flow, respectively, as the Stokes number of bubbles approaches zero.

(3)

When the particles and the fluid have the same inlet velocity, our results show that the particles heavier than the fluid attains its minimum at a location close to the entrance of the flow and after then monotonically increases with increasing radial coordinate and converges to a finite value, and the minimum volume fraction and its location approach the inlet particle volume fraction and the entrance location of the flow as the Stokes number of heavy particles approaches zero. On the other hand, the volume fraction of light particles attains its maximum at a location close to the entrance of the flow and after then monotonically decreases with increasing radial coordinate and converges to a finite value, and the maximum volume fraction of light particles and its location approach the inlet particle volume fraction and the entrance location of the flow as the Stokes number of light particles approaches zero