Boussinesq approximation (buoyancy)

From Wikipedia, the free encyclopedia
Jump to: navigation, search
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[edit]

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






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






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},






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.


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.


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).


  1. ^ a b c d e f g Tritton, D. J. (1977). Physical fluid dynamics. New York: Van Nostrand Reinhold Co. ISBN 9789400999923. 

Further reading[edit]