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.

5 Attachment(s)
Quote:
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 pressurefield 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 2D 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?

3 Attachment(s)
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? 
2 Attachment(s)
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.8e5 {m} but question is remained, 1why axisymetric modeling shows dependent result? 2is there any way to solve this problem for axisymetric domain? 3is there any method which can avoid us from this very fine mesh and give us reasonable result in coarser mesh? 
Hi
1Because in the axisymetric case you really simulate a bubble near the wall with 90 degree contact angle. 2AFAIK nope. 3Use a more precise scheme. 
Dear Nima,
I got same problem with interFoam, did you solve it ? please advise me. Thanks. 
All times are GMT 4. The time now is 09:09. 