[Saideep]

Thanks both!!


I have tried cfMesh earlier but instead of prism cells I have diamond cells and the behavior is similar to Fig 5 as attached in the earlier post. As you see in Fig5, adding layers pushes the prism cells into the domain and causes that sort of an interface (in simple words, the prism cells are trying to divide or layer the flow).

Ofcourse, this looks like a terribly simple geometry to model (and do doubt that sHM is doing its bit to capture the shape of stl) and solve but as you see it is not straight forward. I am not sure if you can mesh with blockmesh as you cannot capture the curves accurately any longer.

There could be several possibilities where it is going wrong here:
1. implementation of the contact angle bc in interfaceProperties.C

2. computation of interface curvature,
3. Prism cells.


I think the first two are not the case as they can be mathematically derived so only these few cells are problematic.



No clue how to deal with these.