icoReactingMultiphaseInterFoam - wrong link to surfaceTensionModel
icoReactingMultiphaseInterFoam links with the surfaceTensionModel from
src/transportModels/interfaceProperties/surfaceTensionModels instead the one from the phaseSystem library located in
As a result, the
type keyword in the
phaseProperties file for specifying the surfaceTensionModel is not read by the solver, the setting has no influence. Only the sigma value is needed. Furthermore, one cannot simply create a new surfaceTensionModel in the
phasesSystem library as one has to use the
interfaceProperties library from transportModels.
The developers are also aware of this issue. Looking at
Make/file in the phasesSystem library reveals the following lines:
/* Ununsed? */ /* surfaceTension = interfaceModels/surfaceTensionModels $(surfaceTension)/surfaceTensionModel/surfaceTensionModel.C $(surfaceTension)/constantSurfaceTensionCoefficient/constantSurfaceTensionCoefficient.C */
So the surfaceTensionModels are not even compiled. Unfortunately, simply uncommenting those lines and removing the links to the
src/transportModels/interfaceProperties does result in other errors during compilation.
How should it behave?
The correct surfaceTensionModels from
src/phaseSystemModels/multiphaseInter/phasesSystem/interfaceModels should be used for
icoReactingMultiphaseInterFoam for easier maintenance and extension of existing surfaceTensionModels.
- OpenFOAM version : v2106
- Operating system : Debian 10
- Compiler : Gcc 9.1.0