Computer simulation of thermal convection in Rayleigh-Bénard cell ...
Computer simulation of thermal convection in Rayleigh-Bénard cell ...
Computer simulation of thermal convection in Rayleigh-Bénard cell ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Hubert Jopek<br />
<strong>Computer</strong> <strong>simulation</strong> <strong>of</strong> <strong>thermal</strong> <strong>convection</strong> <strong>in</strong> <strong>Rayleigh</strong>-<strong>Bénard</strong> <strong>cell</strong><br />
The next problem is the effect <strong>of</strong> the buoyant force on the pocket <strong>of</strong> fluid<br />
which is proportional to the density different between the pocket and its<br />
surround<strong>in</strong>gs. On the other hand the density difference is proportional to the <strong>thermal</strong><br />
expansion coefficient α and the temperature difference<br />
buoyant force may be calculated as:<br />
∆ T<br />
. Consequently the<br />
δT<br />
F = αρ<br />
0<br />
g∆T<br />
= αρ<br />
0g<br />
∆T<br />
, (4.5)<br />
h<br />
where: ρ0<br />
- the orig<strong>in</strong>al fluid density; g - strength <strong>of</strong> the local gravitational field.<br />
Assum<strong>in</strong>g that buoyant force balances the fluid viscous force the pocket<br />
moves with the constant velocity υ<br />
z<br />
.Hence the displacement through a distance<br />
takes for the pocket a time:<br />
∆ z<br />
τ<br />
d<br />
∆z<br />
= . (4.6)<br />
υ<br />
z<br />
As the viscous force is equal to the viscosity <strong>of</strong> the fluid multiplied by the<br />
Laplacian <strong>of</strong> the velocity, the viscous force may be approximated as follows:<br />
F<br />
v<br />
2<br />
υz<br />
= µ ∇ υz<br />
≈ µ . (4.7)<br />
2<br />
h<br />
Now, by equat<strong>in</strong>g buoyant and viscous force one can obta<strong>in</strong> the υz<br />
expression:<br />
αρ0 ghδT<br />
υz<br />
= ∆z<br />
, (4.8)<br />
µ<br />
10