[mykkujinu2201]

Dear all,

I have a question about implementation of coupled level-set and volume-of-fluid method (CLSVOF) on OpenFOAM.
I found Sankar's work on the internet. (http://www.tfd.chalmers.se/~hani/kurser/OS_CFD/)
I could get the OF code from the above URL.
However, his work was done on OF 2.3.0 and I am using OF 4.0 so I need to modify the code to run.
I share the codes here. (deleted)
There are two types of code; sclsVOFsfFoam which use level-set method only for the modification of surface tension and sclsVOFFoam which additionally modify density and viscosity with Heaviside function.

I tested the code for the given tutorial as known as 'rising bubble'.
The codes worked but the problem is that I can not get the same result which Sankar presented on his slides and report.
As you may see the figure I share, the result of scslVOFsfFoam looks similar to that of Sankar's presentation.
For scslVOFFoam, the result is weird; as the bubble rises, it shows unrealistic shape as attached.

I modified the code based on interFoam.
I mean that I started from interFoam and added the additional parts including level-set function generation.
It would be appreciate and really helpful if someone could help me for this issue.

Thank you in advance.

[mykkujinu2201]

I tested the code line by line and I found that the density correction part in solveLSFunction.H might be the reason.

rho == limitedH*rho1 + (1.0 - limitedH)*rho2;

If I disable this part so the density update is done as in interFoam,
the result is same as that of sclsVOFsfFoam.
Is there any problem with this part?

[mykkujinu2201]

Dear forum,

Though I couldn't find the direct reason of the wrong result,
I could find the other sclsVOFFoam solver from Amir Taqieddin based on Takuya Yamamoto on github. (https://github.com/amirtaqi/sCLSVOF_mass_transfer)
The code is combined one with sclsVOFFoam and additional equation for mass transfer.
The code was made for OF 4.1 but there is no problem to run the code and the result was satisfied.
The code is a little bit different from that of Sankar's.
I will compare the codes to find the reason.

Best regards,
Jun

[Jasper Z]

Originally Posted by mykkujinu2201

Dear forum,

Though I couldn't find the direct reason of the wrong result,
I could find the other sclsVOFFoam solver from Amir Taqieddin based on Takuya Yamamoto on github. (https://github.com/amirtaqi/sCLSVOF_mass_transfer)
The code is combined one with sclsVOFFoam and additional equation for mass transfer.
The code was made for OF 4.1 but there is no problem to run the code and the result was satisfied.
The code is a little bit different from that of Sankar's.
I will compare the codes to find the reason.

Best regards,
Jun

Dear Jun,

Have you figured out the issue in sclsVOFFoam of Sankar's work? I simulated the same result as your first figure.

Kind regards,
Jasper