Update Meeting notes authored by Maric, Tomislav's avatar Maric, Tomislav
......@@ -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.
`//Jun: reference to reconstructionSchemes
```
//Jun: reference to reconstructionSchemes
const reconstructionSchemes& surf = advector.surf();
const point x0 = surf.centre()[faceOwner[faceI]];
25 // cell interface unit normal
26 vector n0 = surf.normal()[faceOwner[faceI]];
`
```
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!
`
```
// cell interface centre
const point x0 = surf.centre()[faceOwner[faceI]];
// cell interface unit normal
......@@ -41,7 +42,7 @@ Loop over all boundaries of alpha, ask if a boundary patch is a wall, loop over
n0 /= (mag(n0)+SMALL);
cutFace.calcSubFace(faceI, n0, x0);
alphaf[faceI] = mag(cutFace.subFaceArea()) / sfArea;
`
```
A global object that you initialize outside the boundary field loop
......
......