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



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

    温度源是写在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);
    

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



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



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



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



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2