OpenFOAM常用头文件(一)——createPhi.H

Info<< "Reading/calculating face flux field phi\\n" << endl;

surfaceScalarField phi
(
    IOobject
    (
        "phi",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT, 
        IOobject::AUTO_WRITE
    ),
    fvc::flux(U)  
);

OpenFOAM的求解器都是基于有限体积法进行求解的。有限体积法的一大特点就是利用散度定理,将体积分转化成面通量之和来处理梯度、散度项。通量的构造与重构在求解器中经常能看到,其中最常见的就是createPhi.H了。phi是速度通量。这个H文件生成由速度场构造而成的通量场:

fvc::flux(U) 相当于 fvc::interpolate(U) & mesh.Sf()

而当需要重构时,则需要用到reconstruct函数。

© 版权声明
THE END
喜欢就支持一下吧
点赞0 分享
评论 抢沙发
头像
欢迎您留下宝贵的见解!
提交
头像

昵称

取消
昵称表情代码图片