Curvature of the interface using interFoam
I am using interFoam to simulate multiphase flow in different geometries and I tried to write a little piece of code (writeCurvature.C that i call as a post processing) to calculate the curvature of the interface following the same method that openFoam normally use in the VOF method. (the calculations that i have done are identical to those in interfaceProperties.C). The problem is that I get strange result. I expected to find high values near the interface and values close to zero in the bulk but it is not!
i attached 2 jpg. one is the colour plot of alpha1 and the other one is the colour plot of the curvature that i got.
someone who has any idea why I find this curvature?
These are the calculations I've done:
|All times are GMT -4. The time now is 22:35.|