2-D Rayleigh-Benard convection

There are input files

    • ray1.rea which has 10 element mesh
    • ray2.rea which has 9 element mesh

The explanations in the examples say that it is better to use odd number of elements. The Rayleigh number is set in the variable param(76). Actually inside the code the Rayleigh number is calculated as

Ra = param(76) * ( 1 + param(75))

The critical value is about 1707 so we can try the simulation at 1700 and 1750 for example. The usr file is ray0.usr so to run with 9 elements

cp ray2.rea ray0.rea
cp ray2.map ray0.map
makenek ray0
nekb ray0
(Wait to finish)
visnek ray0
visit