InterFoam mesh dependency
 

Hello Friends
i made a simulation by InterFoam for bubble rising in stagnant liquid,
it seems solvers result is mesh dependent and make mesh refined does not make it free from mesh?
has any body similar experience with interFoam?

Hi
Do you have contact angle?

This is probably due to spurious currents. It is known that they get worse on mesh refinement.

no i dont have

however to make it much more clear i put here the result of an spherical bubble rising in stagnant liquid for different mesh,shape is the same but final position in the same time is deferent!

Also, to calculate the terminal velocity of bubble, i calculate the mass gravity of bubble center, then i calculate the slope of mass gravity position in different time,

terminal velocity for different mesh is:
case: terminal velocity (m/s)
40: 0.036
80: 0.030
160: 0.024
320: 0.016

Hi Nima,

I have several (!) cases where I experience the same problems. Droplets, bubbles with or without contact to the wall, hex tet or poly mesh, you name it.

I did not yet find out where it exactly comes from, but I hope to figure it out some time :)

Just want to let you know, that you are not the only person having these problems! I'll let you know if I have anything new.

Edit: Check also the pressure-field in some cases, they differed extremely from each other due to mesh refinement!

regards
Nicklas

Hi Nima
Is this a symmetry BC? If yes this is numerically similar to normal contact angle.Would you please examine with the complete 2-D domain and let me know the results.

nope it is axisymetric modeling, so frontAndBack are wedge and axis is empty!
you can find the case in above attachment, but i will try it with 2D simulation

As suggested by Ata, could you try your case considering a planar 2D simulation with the whole bubble, and see if the problem persists?

Hello
i made a simulation of 2D bubble rising, it seems problem is there:
case: terminal velocity {m/s}
40: 0.0486
80: 0.0411
160: 0.0377

Hi
Would you please attach picture of bubbles on the mesh simultaneously.

whats your mean ata? do you want surface plot with edge to see mesh?
the mesh is uniform hexahedral

Hi
Yes. I want to see your mesh resolution.
An other issue. How much is your density and viscosity ratio? How much is the capillary number?

Hi
it has been solved with finer mesh for 2D case, this time i examined (320x1600) and terminal velocity was 0.0362 {m/s}, so based on previous result in above post, the geometry with 160x800 is suitable for simulation, so the cell size is about 2.8e-5 {m} but question is remained,
1-why axisymetric modeling shows dependent result?
2-is there any way to solve this problem for axisymetric domain?
3-is there any method which can avoid us from this very fine mesh and give us reasonable result in coarser mesh?

Hi
1-Because in the axisymetric case you really simulate a bubble near the wall with 90 degree contact angle.
2-AFAIK nope.
3-Use a more precise scheme.

Dear Nima,

I got same problem with interFoam,

did you solve it ?
please advise me.

Thanks.