larsgeb / rayleigh-bernard

A simple finite difference scheme for the Rayleigh-Bernard equations (buoyancy driven convection).

A finite difference code to simulate 2D advection-convection, or Rayleigh-Bernard, systems. Has the dependencies:

1. Python 3.7
2. Numpy
3. SciPy
4. Matplotlib

Run the code by changing the configuration values in main.py and entering the following on your terminal:

python main.py

Make a video out of the results by going into the simulations directory and using ffmpeg like this:

ffmpeg -framerate 5 -pattern_type glob -i '*.png' \
  -c:v libx264 -r 30 -pix_fmt yuv420p out.mp4

An example simulation with an unstable initial condition and 3:1 aspect ratio (which will favour 3 convection cells) is given below. This is the parameters currently set in main.py.

rayleigh-bernard's Issues

Question on the governing PDE

Hi Lars, great job! Here I have one question about the governing equations (of T, Psi) you referred to. It looks that the equation for T is

What about the equation for Psi? Any reference or link would be appreciated. Thank you.

help to understand

Hi

I find your code very good. I would like to use it for a very specific application of a low-temperature nitrogen slot between -160°C and -100°C, 1m long, and 0.4m high (possibly with a thickness of 4mm). I would like to adapt your code for a 2D case and possibly try to add a third dimension. What do you think?

Steven