A VoF-based multiphase solver including heat transfer, marangoni flow and laser interaction in OpenFOAM intended for modelling complex surface-tension driven flows
doesn't seem to really compute the correct divergence of the capillary stress tensor, the surface tension itself stays homogeneous. In a temperature gradient field, this can't be right.
This section of
divCapillaryStress()
https://github.com/pzimbrod/thermocapillaryInterFoam/blob/20f3a8574dabb849e723d43e04b45c2d88351df3/phasesSystem/phaseSystem/phaseSystem.C#L1083-L1118
doesn't seem to really compute the correct divergence of the capillary stress tensor, the surface tension itself stays homogeneous. In a temperature gradient field, this can't be right.