本文分析了interFoam求解器根目录下的createFields.H文件，未注释的代码可参考本人的博文

//显示输出动压
Info<< "Reading field p_rgh\n" << endl;
volScalarField p_rgh
(IOobject("p_rgh",runTime.timeName(),mesh,IOobject::MUST_READ,IOobject::AUTO_WRITE),mesh
);Info<< "Reading field U\n" << endl;
volVectorField U
(IOobject("U",runTime.timeName(),mesh,IOobject::MUST_READ,IOobject::AUTO_WRITE),mesh
);#include "createPhi.H"//显示读取transportProperties字典文件
Info<< "Reading transportProperties\n" << endl;
//不互溶不可压缩两相流混合
immiscibleIncompressibleTwoPhaseMixture mixture(U, phi);volScalarField& alpha1(mixture.alpha1());// 相1的相分数
volScalarField& alpha2(mixture.alpha2());//相2的相分数const dimensionedScalar& rho1 = mixture.rho1();//相1密度
const dimensionedScalar& rho2 = mixture.rho2();//相2密度// Need to store rho for ddt(rho, U)
volScalarField rho
(IOobject("rho",runTime.timeName(),mesh,IOobject::READ_IF_PRESENT //存在，则读取),alpha1*rho1 + alpha2*rho2 //根据相分数及密度进行计算
);
rho.oldTime();// Mass flux
surfaceScalarField rhoPhi
(IOobject("rhoPhi",runTime.timeName(),mesh,IOobject::NO_READ,IOobject::NO_WRITE //不进行输出),fvc::interpolate(rho)*phi
);// Construct incompressible turbulence model
//autoPtr--智能指针，可实现类className的指针操作。指针会自动释放，因此无需手动释放
autoPtr<incompressible::turbulenceModel> turbulence
(incompressible::turbulenceModel::New(U, phi, mixture)
);#include "readGravitationalAcceleration.H"
#include "readhRef.H"
#include "gh.H"//总压
volScalarField p
(IOobject("p",runTime.timeName(),mesh,IOobject::NO_READ,IOobject::AUTO_WRITE),p_rgh + rho*gh
);
//设置参考网格和压力参考值
label pRefCell = 0;
scalar pRefValue = 0.0;
setRefCell
(p,p_rgh,pimple.dict(),pRefCell,pRefValue
);if (p_rgh.needReference())
{p += dimensionedScalar("p",p.dimensions(),pRefValue - getRefCellValue(p, pRefCell));p_rgh = p - rho*gh;
}mesh.setFluxRequired(p_rgh.name());
mesh.setFluxRequired(alpha1.name());// MULES flux from previous time-step
surfaceScalarField alphaPhi
(IOobject("alphaPhi",runTime.timeName(),mesh,IOobject::READ_IF_PRESENT,IOobject::AUTO_WRITE),phi*fvc::interpolate(alpha1)
);// MULES Correction
tmp<surfaceScalarField> talphaPhiCorr0;#include "createMRF.H"

几点说明：
关于智能化指针模板类autoPtr和瞬态对象操作模板类tmp可查看

1. https://blog.csdn.net/hanbingchegu/article/details/107122685
2. http://blog.sina.com.cn/s/blog_5fdfa7e60100e867.html