Skip to content

contactAngleForce parallel bug

Summary

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.

Example case

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. Angle_contact

Environment information

  • OpenFOAM version : v2112
  • Operating system : centos
  • Hardware info :
  • Compiler :

Possible fixes

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.