- Run a single case with and without dynamicAlphaContactAngle BC (standard one in OpenFOAM-v2012) and compare the results (spreading area / time).
- Check [`reconstructionSchemes`](https://www.openfoam.com/documentation/guides/latest/api/classFoam_1_1reconstructionSchemes.html)(specifically[`plicRDF`](https://www.openfoam.com/documentation/guides/latest/api/classFoam_1_1reconstruction_1_1plicRDF.html) - not sure about class names, start with `isoAdvection`) class to see if it does anything special for the boundary patch where we use dynamicAlphaContactAngle -
```
//- Interface area normals
volVectorField normal_;
```
in `reconstructionSchemes.C` : cell-centered normals are used, that are specific to the scheme! What does `dynamicAlphaContactAngle` do with the interface normals? Does it fetch the normal field by name using a dictionary, or modify BC for alpha using fixedGradient?