twoLiquidMixingFoam中添加能量方程温度出现发散
-
请教各位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
中的流场数据计算温度场,温度场不影响流场。计算结果发现温度场出现发散,不知道原因是什么?
希望大家能提供一点思路建议。谢谢 -
这是一个非常简单的方程,我觉得可能是你的边界条件不对,不应该出现发散
-
@东岳 谢谢,我也觉得这个方程很简单,这段代码放到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; } }