自己写的求解器和case,运行时候出错，啥原因啊？

[mark]

[CFD中文网]

只能看出来phaseChangeTwoPhaseMixure::vDotAlphal()函数有问题 :slight_smile:

[mark]

错误提示发散了，该从哪里找原因呢

[wwzhao]

时间步长9e-9？另外这才刚开始算第一步，不是发散:astonished:

[mark]

错误显示分母为0；
Foam::phaseChangeTwoPhaseMixture::vDotAlphal() const
{
volScalarField alphalCoeff(1.0/rho1() - alpha1_*(1.0/rho1() - 1.0/rho2()));
Pair<tmp<volScalarField> > mDotAlphal = this->mDotAlphal();
rho1() rho2() 这个地方，公式没错，该怎样修改呢

[mark]

我做的是可压缩空化，对于密度是变化的，case里面需要定义初始密度？

[wwzhao]

@mark 能否提供测试代码和算例？