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函数。