Skip to content

icoReactingMultiphaseInterFoam - wrong link to surfaceTensionModel

Summary

icoReactingMultiphaseInterFoam links with the surfaceTensionModel from src/transportModels/interfaceProperties/surfaceTensionModels instead the one from the phaseSystem library located in src/phaseSystemModels/multiphaseInter/phasesSystem/interfaceModels/surfaceTensionModels .

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.

Environment information

  • OpenFOAM version : v2106
  • Operating system : Debian 10
  • Compiler : Gcc 9.1.0