@东岳 谢谢，我也觉得这个方程很简单，这段代码放到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;
}
}
