A Boussinesq fluid is governed by the equations of heat, motion and continuity (e.g. Kimura, 1983):
where, T is temperature,
is velocity,
is density, P is pressure,
is the thermal expansion coefficient,
is thermal diffusivity,
is kinematic viscosity,
is acceleration of gravity.
We introduce the stream function
![]()
as
Subsituting this into (14) and (15), we obtain(17)
where J(A,B) is the Jacobian
(20) The governing equations are non-dimensionalized using the thickness of the fluid layer b as the length scale, the thermal diffusion time
as the time scale, and the temperature difference between the boundaries
as the temperature scale, yielding,
Here
![]()
denotes a non-dimensional variable. Equations (18) and (19) then become
where,is vorticity,
is the Prandtl number, and
is the Rayleigh number. The subscirpts
denoting a non-dimensional variable are omitted for simplicity.
A fixed temperature condition is specified at the top and the bottom boundaries:
The kinematic and dynamic boundary conditions are impermeable, and free-slip, such that
(23)
(24) Horizontal boundaries are periodic.
(25) These kinematic and dynamic conditions are described by the voriticity and the stream function. The impermeable condition becomes
Since the initial condtion is one of no motion and
(26) , stream function always vanishes at the boundaries,
Since w=0 at the horizontal plane of the boundaries, the free-slip condition becomes
(27) that is,
Further, since
(28) at the horizontal plane of the boundaries, we have,
By following the similar procedure for vertical plane of the boundaries, we obtain (29).
(29)