implement level set to compressibleInterFoam
Dear Foamers,
Does anybody implement level set method into compressibleInterFoam? Could you please tell me some theory guide? Best Regards CYW |
I have been working on it in my free time so I haven't had a significant progress yet.
Here is a good starting point. Look for "Coupled Level-Set with VOF interFoam, Sankar Menon (2014)" in the link and you can find the report as well as the files. |
level set in compressibleInterFoam
Quote:
|
I will resume the work next month, will let you know if I have any achievement in this regard and if interested we can work together.
|
Quote:
|
Sorry, I've got my hands full this month. If you have any success please let me know. If you have any specific errors regarding implementation also maybe I can help.
|
Quote:
|
Parasitic currents are a result of the way surface tension is implemented in OF so you have pay extra attention to the way compressibleInterFoam deals with surface tension modeling as it is different than interFoam. So implementing LS is not as straight forward as interFoam.
|
surface tension in compressibleInterFoam
Quote:
which is interface.surfaceTensionForce(). |
In case of a quiescent flow without gravity, lack of balance in surface tension and pressure gradient forces, and errors in curvature calculations is the source of the currents.
So I would suggest to look into this case and see how compressibility makes a difference in this simple case. At least that's what I was trying to do and understand. Just put a square bubble in a domain where the boundaries are 4 times the radius away from the bubble, remove the gravity and try it with the models that you have and figure out what are the differences. |
Quote:
|
I see.
In general I don't think coupling with LS provide a significant advantages in this case, otherwise they would have implemented it officially. The CIF solver is complex, adding more complexity wouldn't help I suppose. There are some alternatives though. If you need sharp interface you can try using isoAdvector solver. It's already a part of OF v1706, here. There are some nice comparison between interFoam and CLSVOF in this paper. interFoam performs very well. Also, you can try implementing a force balancing algorithm which are specifically for this purpose, for example this paper. |
Have you simulated the dambreak case in the tutorial with sclsVOFFoam? What is the result?
In the report, the sclsVOFFoam can only handle zerogradient boundary condionts? What does that mean? How to handle other boundary conditions? |
Quote:
Did you manage to find out how to handle other boundary conditions in the sclsVOFFoam solver? |
All times are GMT -4. The time now is 14:32. |