[überschwupper]

Peace Foamers,


I'm currently implementing an evaporation model and for that, I need the interface Area between my two phases alpha1 and alpha2.


Does anyone know how to do that? I assumed that one of the interface reconstruction methods would deliver a function that could be called or the class surfaceInterpolation but I havent found anything suitable there. (btw. I'm not greatly familiar with PLIC, isoAdvector, etc.)


My first thought of an apporach looked like that:
1.) Identifying interface cells (e.g. alpha = 0.5)
2.) Then calculate the volume occupied by alpha for each of these interface cells (mesh.V()*alpha1)

3.) calculate normal vector with nHat = grad(alpha)/mag(grad(alpha)) to get information about the orientation of that approximated interface plane
______________ Now the part where I dont know how to proceed
Here I thought about an algorithm that is starting from a suitable point. I try to explain it in a simple case: 2D Square. Starting from a corner then follow the diagonal and calculate a orthogonal line where both ends of this line touch the sides of that square. This will be continued until the areas (equivalent to the volume in 2.)) will match. The resulting length(area) is what can be used for the calculation.


Has someone an idea how to realize that - or even better: Is there something which gives me the interface area between both phases without a complicate implementation?


Kind regards and a thank you in advance!