31.03.2014 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!