Basset–Boussinesq–Oseen equation

Short description: Equation describing particle motion in unsteady, low-Reynolds fluid flow

In fluid dynamics, the Basset–Boussinesq–Oseen equation (BBO equation) describes the motion of – and forces on – a small particle in unsteady flow at low Reynolds numbers. The equation is named after Joseph Valentin Boussinesq, Alfred Barnard Basset and Carl Wilhelm Oseen.

Formulation

The BBO equation, in the formulation as given by (Zhu Fan) and (Soo 1990), pertains to a small spherical particle of diameter [math]\displaystyle{ d_p }[/math] having mean density [math]\displaystyle{ \rho_p }[/math] whose center is located at [math]\displaystyle{ \boldsymbol{X}_p(t) }[/math]. The particle moves with Lagrangian velocity [math]\displaystyle{ \boldsymbol{U}_p(t)=\text{d} \boldsymbol{X}_p / \text{d}t }[/math] in a fluid of density [math]\displaystyle{ \rho_f }[/math], dynamic viscosity [math]\displaystyle{ \mu }[/math] and Eulerian velocity field [math]\displaystyle{ \boldsymbol{u}_f(\boldsymbol{x},t) }[/math]. The fluid velocity field surrounding the particle consists of the undisturbed, local Eulerian velocity field [math]\displaystyle{ \boldsymbol{u}_f }[/math] plus a disturbance field – created by the presence of the particle and its motion with respect to the undisturbed field [math]\displaystyle{ \boldsymbol{u}_f. }[/math] For very small particle diameter the latter is locally a constant whose value is given by the undisturbed Eulerian field evaluated at the location of the particle center, [math]\displaystyle{ \boldsymbol{U}_f(t)=\boldsymbol{u}_f(\boldsymbol{X}_p(t),t) }[/math]. The small particle size also implies that the disturbed flow can be found in the limit of very small Reynolds number, leading to a drag force given by Stokes' drag. Unsteadiness of the flow relative to the particle results in force contributions by added mass and the Basset force. The BBO equation states:

[math]\displaystyle{ \begin{align} \frac{\pi}{6} \rho_p d_p^3 \frac{\text{d} \boldsymbol{U}_p}{\text{d} t} &= \underbrace{3 \pi \mu d_p \left( \boldsymbol{U}_f - \boldsymbol{U}_p \right)}_{\text{term 1}} - \underbrace{\frac{\pi}{6} d_p^3 \boldsymbol{\nabla} p}_{\text{term 2}} + \underbrace{\frac{\pi}{12} \rho_f d_p^3\, \frac{\text{d}}{\text{d} t} \left( \boldsymbol{U}_f - \boldsymbol{U}_p \right)}_{\text{term 3}} \\ & + \underbrace{\frac{3}{2} d_p^2 \sqrt{\pi \rho_f \mu} \int_{t_{_0}}^t \frac{1}{\sqrt{t-\tau}}\, \frac{\text{d}}{\text{d} \tau} \left( \boldsymbol{U}_f - \boldsymbol{U}_p \right)\, \text{d} \tau}_{\text{term 4}} + \underbrace{\sum_k \boldsymbol{F}_k}_{\text{term 5}} . \end{align} }[/math]

This is Newton's second law, in which the left-hand side is the rate of change of the particle's linear momentum, and the right-hand side is the summation of forces acting on the particle. The terms on the right-hand side are, respectively, the:

1. Stokes' drag,
2. Froude–Krylov force due to the pressure gradient in the undisturbed flow, with [math]\displaystyle{ \boldsymbol{\nabla} }[/math] the gradient operator and [math]\displaystyle{ p(\boldsymbol{x},t) }[/math] the undisturbed pressure field,
3. added mass,
4. Basset force and
5. other forces acting on the particle, such as gravity, etc.

The particle Reynolds number [math]\displaystyle{ R_e: }[/math]

[math]\displaystyle{ R_e = \frac{\max\left\{ \left| \boldsymbol{U}_p - \boldsymbol{U}_f \right| \right\}\, d_p}{\mu/\rho_f} }[/math]

has to be less than unity, [math]\displaystyle{ R_e \lt 1 }[/math], for the BBO equation to give an adequate representation of the forces on the particle.

Also (Zhu Fan) suggest to estimate the pressure gradient from the Navier–Stokes equations:

[math]\displaystyle{ -\boldsymbol{\nabla} p = \rho_f \frac{\text{D} \boldsymbol{u}_f}{\text{D} t} - \mu \nabla^2 \boldsymbol{u}_f, }[/math]

with [math]\displaystyle{ \text{D} \boldsymbol{u}_f / \text{D} t }[/math] the material derivative of [math]\displaystyle{ \boldsymbol{u}_f. }[/math] Note that in the Navier–Stokes equations [math]\displaystyle{ \boldsymbol{u}_f(\boldsymbol{x},t) }[/math] is the fluid velocity field, while, as indicated above, in the BBO equation [math]\displaystyle{ \boldsymbol{U}_f }[/math] is the velocity of the undisturbed flow as seen by an observer moving with the particle. Thus, even in steady Eulerian flow [math]\displaystyle{ \boldsymbol{u}_f }[/math] depends on time if the Eulerian field is non-uniform.

References

• Zhu, Chao; Fan, Liang-Shi (1998). "Chapter 18 – Multiphase flow: Gas/Solid". in Johnson, Richard W.. The Handbook of Fluid Dynamics. Springer. ISBN 9783540646129. 
• Soo, Shao L. (1990). Multiphase Fluid Dynamics. Ashgate Publishing. ISBN 9780566090332.