多相流如何提取扭矩

[派大星huangxiaojie]

各位老师大家好，我正在使用openfoam10中的multiphaseEulerFoam求解器，模拟的为气液两相流旋转机械，目前我在扭矩的提取上碰到了问题。以下是我的functions：

forces
{
type forces;
libs ("libforces.so");

writeControl timeStep;
writeInterval 1;

rho rhoInf;

patches (WALLIMP);

rhoInf 1000;

CofR (0 0 0);

pitchAxis (0 0 -1);

U U.water;
p p;
}

当我运行计算时，出现了如下提示：

--> FOAM FATAL ERROR:
No valid model for viscous stress calculation

From function Foam::tmp<Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> > Foam::functionObjects::forces::devTau() const
in file forces/forces.C at line 288.

FOAM exiting

我该如何解决？谢谢各位老师！