MPPICInterFoam中的bug?



  • 在ESI的版本中,MPPICInterFoam中的pEqn:

    fvScalarMatrix p_rghEqn
    (
        fvm::laplacian(Dp, p_rgh)
       ==
       fvc::ddt(alphac) + fvc::div(alphacf*phiHbyA)
    );
    

    其中的fvc::ddt(alphac)为0。 Jasak也表示:

    Thus, if you try to evaluate the ddt the way it is written, it will recognise that the field does not have the oldTime() level stored and will use the current value to initialise it. Basically, fvc::ddt(magSqr(U)) will always give you zero, which is wrong.

    是个bug?将fvc::ddt(alphac)去掉之后求解结果没变化。



  • @李东岳 有另外一个疑问,在MPPICInterFoam的UEqn中的第7行,为什么没有alphac项呢,0_1512089257527_UEqn.JPG ,即是否应为:
    fvm::Sp(fvc::ddt(alphacRho))+fvc::div(alphaRhoPhic)?



  • @linhan-ge 抱歉上面的公式有误,应为:
    fvm::Sp(fvc::ddt(alphacRho)+fvc::div(alphaRhoPhic),U)?



  • 不确定,我还没有细看这个求解器。你可以推一下看看,哈哈。


登录后回复
 

与 CFD中文网 的连接断开,我们正在尝试重连,请耐心等待