twoLiquidMixingFoam中添加能量方程温度出现发散

[chengan.wang]

请教各位CFDers
我打算在twoLiquidMixingFoam中添加能量方程
\begin{equation}
\frac{\partial T}{\partial t }+\nabla \cdot (\mathbf{U}T)= \nabla \cdot (a \nabla T)
\end{equation}
我在原始求解器中加入相应的代码如下：

{
    fvScalarMatrix TEqn
    (
        fvm::ddt(T)
        + fvm::div(phi, T)
        - fvm::laplacian(DT, T)
    );

    TEqn.solve();
}

本质上就是利用twoLiquidMixingFoam中的流场数据计算温度场，温度场不影响流场。计算结果发现温度场出现发散，不知道原因是什么？
希望大家能提供一点思路建议。谢谢

[李东岳]

这是一个非常简单的方程，我觉得可能是你的边界条件不对，不应该出现发散

[chengan.wang]

@东岳 谢谢，我也觉得这个方程很简单，这段代码放到icoFoam里面没有出现过问题。但是放在两种液体混合求解器中就不行了。我感觉边界条件没有问题，流体从一侧流入，另一侧流出。

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      T;
}
dimensions      [0 0 0 1 0 0 0];

internalField   uniform 293;

    inlet
    {
        type            groovyBC;
        valueExpression "(pos().y>1) ? (295):(293)"; 			
    }

    outlet
    {
        type            zeroGradient;	
    }

    freeSurface
    {
        type            zeroGradient;
    }

    bottom
    {
        type            fixedValue;
        value           293;			
    }

    front
    {
        type            symmetryPlane;	
    }

    back
    {
        type            symmetryPlane;		
    }

FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (1 0 0); 

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform (1 0 0);
    }

    outlet
    {
        type            inletOutlet;
        inletValue      uniform (0 0 0);
        value           $internalField;		
    }

    freeSurface
    {
        type            pressureInletOutletParSlipVelocity;
        value           uniform (1 0 0);
    }

    bottom
    {
        type            fixedValue;
        value           uniform (0 0 0);			
    }

    front
    {
        type            symmetryPlane;		
    }

    back
    {
        type            symmetryPlane;	
    }

}