请教扩散现象不随时间变化的原因？

[merlinizlf]

大家好：
请教一个问题，添加温度源的时候，发现一段时间后就不再扩散，明显不符合常识啊，会是什么原因？在我的例子中t=6之后就不再变化，如图：

温度源是写在system/fvOptions里的，SRC是topoSet中分出来的一个区域，附上代码：

TSource
{
    type            scalarSemiImplicitSource;
    active          true;
    scalarSemiImplicitSourceCoeffs
    {
        selectionMode   cellSet; // all, cellSet, cellZone, points
        cellSet         SRC;
        volumeMode      absolute; // absolute or specific
        injectionRateSuSp
        {
            T     (10 0);
        }
    }
}

求解器是基于pisoFoam写的，其中添加的温度方程参照scalarTransportFoam写的，代码：

fvScalarMatrix TEqn
(
    fvm::ddt(T)
    + fvm::div(phi, T)
    - fvm::laplacian(DT, T)
    == fvOptions(T)
);
TEqn.relax();
fvOptions.constrain(TEqn);
TEqn.solve();
fvOptions.correct(T);

请大家不吝解惑，欢迎任何评论

[Samuel-Tu]

T是瞬态还是稳态，是不是速度稳定了，热源也稳定，求解达到稳态了。另外你的TEqn放的位置在哪里，是在pisoLoop外面还是里面

[merlinizlf]

@Samuel-Tu 非常感谢您的回复。T是瞬态，但是，TEqn跟U只是单向耦合。我的TEqn放在了pisoLoop的外面，这一点我没有考虑过，请问会造成什么影响？因为我看openfoamwiki上icoFoamTemp修改是放在了runTime.loop外面，所以我也放在了同一位置
非常期待您的回复
附上ofwiki的链接链接文本

[Samuel-Tu]

@merlinizlf 看了下wiki的代码，TEqn是在pisoLoop里面。只是在压强迭代外面。pisoFoam速度场也可能是稳定的啊，你想一个水管，里面没有阻挡物，入口边界不变，层流，一直放水，那个速度场基本会很快稳定。你看看你的速度场６s后是不是基本稳定。速度场稳定了，源项也是稳态的。最终温度分布也会稳定吧

[merlinizlf]

@Samuel-Tu 非常感谢您的回复。讲的非常明白，我之前没有想到这些。其实，我想实现的就是一个热源点慢慢加热最终提高出口温度，或许我应该试试改成热源强度再试一试看。非常感谢您。