contactAngleForce parallel bug
When using the finiteArea to study thin films, it is possible to use the 'perturbedTemperatureDependantContactAngle' force, but when using parallelization, a problem occurs at the boundary between processors; the contactAngle boundary condition between processors is set to 'uniform (0 0 0)' instead of 'nonuniform List...'.
Steps to reproduce
Using the 'perturbedTemperatureDependantContactAngle' force with the finiteArea method and using parallelization.
Open the tutorial 'incompressible/pimpleFoam/laminar/inclinedPlaneFilm', set the Ccf in 0/U to 1 then use Allrun-parallel.
What is the current bug behaviour?
Some fluid seems to be imprisoned at the boundary between processors, but if you observe the 'contactAngle:contactForce' in Paraview, you will notice it is set to 0 at the boundary.
What is the expected correct behavior?
The film behavior should not be changed by the communication between processors.
Relevant logs and/or images
This picture is the simulation of a falling thin film, the domain is decomposed in 6 regions you can clearly see. Red/blue is the film height scale and dark is the contactAngleForce scale.
- OpenFOAM version : v2112
- Operating system : centos
- Hardware info :
- Compiler :
I don't know if the problem comes from the angle itself or from the force, but changing the boundary condition between processor should be the solution.