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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

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

5.4. Boundary Conditions<br />

The boundary conditions for the temperature are as follows:<br />

z = 0<br />

z = 1<br />

→<br />

→<br />

θ = 0<br />

. (5.27)<br />

θ = 0<br />

It is so because <strong>of</strong> the fact that the temperature at the top and the bottom is fixed.<br />

Boundary conditions for the streamfunction - let the shear forces at the top and at<br />

the bottom be neglected:<br />

z = 0<br />

z = 1<br />

→<br />

→<br />

∂υ<br />

x<br />

∂z<br />

∂υ<br />

x<br />

∂z<br />

= 0<br />

= 0<br />

. (5.28)<br />

The follow<strong>in</strong>g expressions satisfy assumed conditions:<br />

Ψ(<br />

x,<br />

z,<br />

t)<br />

= ψ ( t)s<strong>in</strong>(<br />

πz)s<strong>in</strong>(<br />

ax)<br />

θ(<br />

x,<br />

z,<br />

t)<br />

= T1(<br />

t)s<strong>in</strong>(<br />

πz)cos(<br />

ax)<br />

−T2(<br />

t)s<strong>in</strong>(2πz)<br />

where the parameter a is to be determ<strong>in</strong>ed.<br />

, (5.29)<br />

The function Ψ is this part <strong>of</strong> model which is responsible for aris<strong>in</strong>g convective<br />

rolls which can be observed <strong>in</strong> real experiment. The second equation is the<br />

temperature deviation function which consists <strong>of</strong> two parts. The former part<br />

T1<br />

describes the temperature difference between the upward and downward mov<strong>in</strong>g<br />

parts <strong>of</strong> a convective <strong>cell</strong>, while the latter is the description <strong>of</strong> the deviation from the<br />

l<strong>in</strong>ear temperature variation <strong>in</strong> the centre <strong>of</strong> a convective <strong>cell</strong> as a as a function <strong>of</strong><br />

vertical position z .<br />

22

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

Saved successfully!

Ooh no, something went wrong!