Lid driven cavity Matlab example with Re=1000

Lattice Boltzmann Method LBM: Algorithms
#1

Hi everyone,

I’m starting with LBM and I’ve been running some tests with the Lid Driven Cavity problem written in Matlab and available at the wiki page (http://wiki.palabos.org/_media/numerics:cavity.m)

The problem is that when I try to change the Reynolds number from 100 to 1000 in this code I get some errors in the calculations. I have no idea why.

Do I need to change another parameter as well? Is there a limitation in this code? Could someone please give me some hints?

Best regards,
BMR

#2

Hi,
I have a quick check on your code. It was really efficient.
But about your problem, I have to say that this is a common error in this problem. I have also the same error when changing Re from 100 to 1000. As I understood, the main source of this unstability is omega. I do not know exactly, but if you change omega to, say 1.0, you will not see the unstability in the results.