too dense mesh deteriorates the simulation
 

hi,
I'm using pimpleFoam to simulate the flow around a circular cylinder. I found that many times that if a dense mesh is used for a low Reynolds number simulation, say Re=100, the simulation is not physically correct(in terms of the drag or lift force coefficients). However, if I maintain the same mesh and raise the reynolds number to like 1000, the simulation seems to go well.

does someone else also have the same experience with me?

Hi kkpal, this is very different to the normal results improvement due to mesh elemnt size decrease and it is not aligned to numerical theory...have you checked everything in your model? are you using wall resolved boundary layers in particular?

y+=1
 

Right, when using lowRe Models without wall functions, the refinement of the wall depends on the Reynolds number. You have to achieve a y+=1 for the first layer. When changing the Reynolds number from 100 to 1000, the first layer thickness to achieve y+=1 varies with round about factor 5.

Cheers

Fabian

Thank you, Andrea and Fabian,
I am not using any wallfunctions. Just pure DNS for low Re flows.
I think I've found the root to the problem.
It turned out a larger value of nCorrectors(in the fvSolution under PIMPLE directory) needs to be used.
While in the case of coarse mesh and small Re(100), a value of nCorrectors=1 suffices, when the mesh is refined, a larger value, say nCorrectors=3 needs to be used.
Another detail: I'm using fixed courant about 2 to dynamic adjust the timestep. So I think my previous problem is not a result of the timestep.
To be honest I do not know the reason behind this:confused::confused:. I would appreciate it if someone could give me some hints.

It is very out of theory that a fine mesh gives worst result than a coarse mesh...it sounds to me like what you get with the coarse mesh is a lucky result, while the fine grid tells that something is not resolved properly... I don't have experience with DNS, so I cannot help further in detail. But, are you comparing cfd results with experiment at same Reynolds number? are you comparing lift and drag coefficients only or do you have any other data available?

Which kind of mesh do you use, structured, unstructured ... ?

See the post Weird results on Cavity and unstructured meshes.