You are currently browsing the tag archive for the ‘Boussinesq equations’ tag.
The Boussinesq equations for inviscid, incompressible two-dimensional fluid flow in the presence of gravity are given by
where is the velocity field,
is the pressure field, and
is the density field (or, in some physical interpretations, the temperature field). In this post we shall restrict ourselves to formal manipulations, assuming implicitly that all fields are regular enough (or sufficiently decaying at spatial infinity) that the manipulations are justified. Using the material derivative
, one can abbreviate these equations as
One can eliminate the role of the pressure by working with the vorticity
. A standard calculation then leads us to the equivalent “vorticity-stream” formulation
of the Boussinesq equations. The latter two equations can be used to recover the velocity field from the vorticity
by the Biot-Savart law
It has long been observed (see e.g. Section 5.4.1 of Bertozzi-Majda) that the Boussinesq equations are very similar, though not quite identical, to the three-dimensional inviscid incompressible Euler equations under the hypothesis of axial symmetry (with swirl). The Euler equations are
where now the velocity field and pressure field
are over the three-dimensional domain
. If one expresses
in polar coordinates
then one can write the velocity vector field
in these coordinates as
If we make the axial symmetry assumption that these components, as well as , do not depend on the
variable, thus
then after some calculation (which we give below the fold) one can eventually reduce the Euler equations to the system
where is the modified material derivative, and
is the field
. This is almost identical with the Boussinesq equations except for some additional powers of
; thus, the intuition is that the Boussinesq equations are a simplified model for axially symmetric Euler flows when one stays away from the axis
and also does not wander off to
.
However, this heuristic is not rigorous; the above calculations do not actually give an embedding of the Boussinesq equations into Euler. (The equations do match on the cylinder , but this is a measure zero subset of the domain, and so is not enough to give an embedding on any non-trivial region of space.) Recently, while playing around with trying to embed other equations into the Euler equations, I discovered that it is possible to make such an embedding into a four-dimensional Euler equation, albeit on a slightly curved manifold rather than in Euclidean space. More precisely, we use the Ebin-Marsden generalisation
of the Euler equations to an arbitrary Riemannian manifold (ignoring any issues of boundary conditions for this discussion), where
is a time-dependent vector field,
is a time-dependent scalar field, and
is the covariant derivative along
using the Levi-Civita connection
. In Penrose abstract index notation (using the Levi-Civita connection
, and raising and lowering indices using the metric
), the equations of motion become
in coordinates, this becomes
where the Christoffel symbols are given by the formula
where is the inverse to the metric tensor
. If the coordinates are chosen so that the volume form
is the Euclidean volume form
, thus
, then on differentiating we have
, and hence
, and so the divergence-free equation (10) simplifies in this case to
. The Ebin-Marsden Euler equations are the natural generalisation of the Euler equations to arbitrary manifolds; for instance, they (formally) conserve the kinetic energy
and can be viewed as the formal geodesic flow equation on the infinite-dimensional manifold of volume-preserving diffeomorphisms on (see this previous post for a discussion of this in the flat space case).
The specific four-dimensional manifold in question is the space with metric
and solutions to the Boussinesq equation on can be transformed into solutions to the Euler equations on this manifold. This is part of a more general family of embeddings into the Euler equations in which passive scalar fields (such as the field
appearing in the Boussinesq equations) can be incorporated into the dynamics via fluctuations in the Riemannian metric
). I am writing the details below the fold (partly for my own benefit).
Recent Comments