Numerical calculations (1) -parameter Convection in a rotating cylindrical annulus with fixed heat flux boundaries Neutral curve (2) Numerical calculations (2)
In order to confirm the convection behavior obtained by analyses of the truncated linear system presented so far, we have performed numerical time integrations of the full system (1)-(4). The computational domain extends to 8 horizontally and to 1 vertically. The numbers of grid points are 256 horizontally and 32 vertically. The cyclic condition is adopted in the x direction. 4th order Runge-Kutta scheme is used for the time integration and non-linear terms are discretized by Arakawa's Jacobian scheme (Arakawa 1966). As the initial condition, a temperature disturbance with an amplitude of 0.01 is placed at x=4 and z=1/4 onto the thermal conductive solution.

The Prandtl number and h are fixed at P=1 and h =1000. The numerical computations are performed for the cases (a) R=10000 and (b) R=8000. In the large Rayleigh number case (a), there are unstable modes with large and small horizontal wave numbers, while in the small Rayleigh number case (b), there are only unstable modes with small horizontal wave numbers.

img55a.gif (19404 bytes)

Fig.5 The values of the Rayleigh number used for the numerical calculations.


Numerical calculations (1) -parameter Convection in a rotating cylindrical annulus with fixed heat flux boundaries Neutral curve (2) Numerical calculations (2)