@@ -20,20 +20,21 @@ with a droplet advected on the bottom wall (Mathis et al.)
...
@@ -20,20 +20,21 @@ with a droplet advected on the bottom wall (Mathis et al.)
Look at [reconstructionSchemes](https://www.openfoam.com/documentation/guides/latest/api/classFoam_1_1reconstructionSchemes.html) member functions `center` and `normal` : those are the PLIC plane point and normal. Use those to calculate contact angle.
Look at [reconstructionSchemes](https://www.openfoam.com/documentation/guides/latest/api/classFoam_1_1reconstructionSchemes.html) member functions `center` and `normal` : those are the PLIC plane point and normal. Use those to calculate contact angle.
Loop over all boundaries of alpha, ask if a boundary patch is a wall, loop over all faces of the boundary patch, use the code snippet above to fetch center and normal, do some geometry to get contact angle.
Loop over all boundaries of alpha, ask if a boundary patch is a wall, loop over all faces of the boundary patch, use the code snippet above to fetch center and normal, do some geometry to get contact angle.
**FYI**: that's how we could calculate alphaf without interpolation!
**FYI**: that's how we could calculate alphaf without interpolation!
`
```
// cell interface centre
// cell interface centre
const point x0 = surf.centre()[faceOwner[faceI]];
const point x0 = surf.centre()[faceOwner[faceI]];
// cell interface unit normal
// cell interface unit normal
...
@@ -41,7 +42,7 @@ Loop over all boundaries of alpha, ask if a boundary patch is a wall, loop over
...
@@ -41,7 +42,7 @@ Loop over all boundaries of alpha, ask if a boundary patch is a wall, loop over