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 , stream function always vanishes at the boundaries,
(26) 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)