Boussinesq approximation (buoyancy)

This article is about the Boussinesq approximation in buoyancy-driven flows. For other uses, see Boussinesq approximation (disambiguation).

In fluid dynamics, the Boussinesq approximation (pronounced: [businɛsk], named for Joseph Valentin Boussinesq) is used in the field of buoyancy-driven flow (also known as natural convection). It ignores density differences except where they appear in terms multiplied by g, the acceleration due to gravity. The essence of the Boussinesq approximation is that the difference in inertia is negligible but gravity is sufficiently strong to make the specific weight appreciably different between the two fluids. Sound waves are impossible/neglected when the Boussinesq approximation is used since sound waves move via density variations.

Boussinesq flows are common in nature (such as atmospheric fronts, oceanic circulation, katabatic winds), industry (dense gas dispersion, fume cupboard ventilation), and the built environment (natural ventilation, central heating). The approximation is extremely accurate for many such flows, and makes the mathematics and physics simpler.

The approximation

The Boussinesq approximation is applied to problems where the fluid varies in temperature from one place to another, driving a flow of fluid and heat transfer. The fluid satisfies conservation of mass, conservation of momentum and conservation of energy. In the Boussinesq approximation, variations in fluid properties other than density $\scriptstyle \rho$ are ignored, and density only appears when it is multiplied by $\scriptstyle g$ , the gravitational acceleration.^[1]^:127–128 If $\scriptstyle \mathbf{u}$ is the local velocity of a parcel of fluid, the continuity equation for conservation of mass is

\frac{\partial\rho}{\partial t} + \nabla\cdot\left(\rho\mathbf{u}\right) = 0.

^[1]^:52

If density variations are ignored, this reduces to

$\nabla\cdot\mathbf{u} = 0.$ ^[1]^:128

(1)

The general expression for conservation of momentum (the Navier–Stokes equations) is

\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u}\cdot\nabla\mathbf{u} = -\frac{1}{\rho}\nabla p + \nu\nabla^2 \mathbf{u} + \frac{1}{\rho}\mathbf{F},

where $\scriptstyle \nu$ is the kinematic viscosity and $\scriptstyle \mathbf{F}$ is the sum of any body forces such as gravity.^[1]^:59 In this equation, density variations are assumed to have a fixed part and another part that has a linear dependence on temperature:

\rho = \rho_0 - \alpha\rho_0\Delta T,

where $\scriptstyle \alpha$ is the coefficient of thermal expansion.^[1]^:128–129 If $\scriptstyle\mathbf{F}=\rho\mathbf{g}$ is the gravitational body force, the resulting conservation equation is

$\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u}\cdot\nabla\mathbf{u} = -\frac{1}{\rho}\nabla p + \nu\nabla^2 \mathbf{u} - \mathbf{g}\alpha\Delta T.$ ^[1]^:129

(2)

In the equation for heat flow in a temperature gradient , the heat capacity per unit volume, $\scriptstyle \rho C_\rho$ , is assumed constant. The resulting equation is

$\frac{\partial T}{\partial t} + \mathbf{u}\cdot\nabla T = \kappa\nabla^2T +\frac{J}{\rho C_\rho},$

(3)

where $\scriptstyle J$ is the rate per unit volume of internal heat production and $\scriptstyle \kappa$ is the thermal diffusivity.^[1]^:129

The three numbered equations are the basic convection equations in the Boussinesq approximation.

Advantages

The approximation's advantage arises because when considering a flow of, say, warm and cold water of density $\rho_1$ and $\rho_2$ one needs only consider a single density $\rho$ : the difference $\Delta\rho= \rho_1-\rho_2$ is negligible. Dimensional analysis shows that, under these circumstances, the only sensible way that acceleration due to gravity g should enter into the equations of motion is in the reduced gravity $g'$ where

g' = g{\rho_1-\rho_2\over \rho}.

(Note that the denominator may be either density without affecting the result because the change would be of order $g(\Delta\rho/\rho)^2$ ). The most generally used dimensionless number would be the Richardson number and Rayleigh number.

The mathematics of the flow is therefore simpler because the density ratio ( $\rho_1/\rho_2$ , a dimensionless number) does not affect the flow; the Boussinesq approximation states that it may be assumed to be exactly one.

Inversions

One feature of Boussinesq flows is that they look the same when viewed upside-down, provided that the identities of the fluids are reversed. The Boussinesq approximation is inaccurate when the nondimensionalised density difference $\Delta\rho/\rho$ is of order unity.

For example, consider an open window in a warm room. The warm air inside is less dense than the cold air outside, which flows into the room and down towards the floor. Now imagine the opposite: a cold room exposed to warm outside air. Here the air flowing in moves up toward the ceiling. If the flow is Boussinesq (and the room is otherwise symmetrical), then viewing the cold room upside down is exactly the same as viewing the warm room right-way-round. This is because the only way density enters the problem is via the reduced gravity $g'$ which undergoes only a sign change when changing from the warm room flow to the cold room flow.

An example of a non-Boussinesq flow is bubbles rising in water. The behaviour of air bubbles rising in water is very different from the behaviour of water falling in air: in the former case rising bubbles tend to form hemispherical shells, while water falling in air splits into raindrops (at small length scales surface tension enters the problem and confuses the issue).

References

1 2 3 4 5 6 7 Tritton, D. J. (1977). Physical fluid dynamics. New York: Van Nostrand Reinhold Co. ISBN 9789400999923.

Boussinesq approximation (buoyancy)

The approximation

Advantages

Inversions

References

Further reading