注意到simpleFoam源码的createField.H中含有一行代码：

singlePhaseTransportModel laminarTransport(U, phi);

将类singlePhaseTransportModel实例化为laminarTransport，并且传给构造函数两个参数U和phi。
在simpleFoam.C中找到了#include "singlePhaseTransportModel.H"。类singlePhaseTransportModel很有可能在里面定义，因此进入singlePhaseTransportModel.H。
其路径为$FOAM_SRC/transportModels/incompressible/singlePhaseTransportModel/singlePhaseTransportModel.H，发现类singlePhaseTransportModel，及其构造函数的声明。因此进入singlePhaseTransportModel.C中查看其定义。

Foam::singlePhaseTransportModel::singlePhaseTransportModel
(
    const volVectorField& U,
    const surfaceScalarField& phi
)
:
    IOdictionary
    (
        IOobject
        (
            "transportProperties",
            U.time().constant(),
            U.db(),
            IOobject::MUST_READ_IF_MODIFIED,
            IOobject::NO_WRITE
        )
    ),
    viscosityModelPtr_(viscosityModel::New("nu", *this, U, phi))
{}

在类singlePhaseTransportModel构造前，先建立了读取constant/transportProperties的字典，并用
viscosityModel::New返回了一个叫viscosityModelPtr_的指针。
因此需要查看viscosityModel::New的定义。路径为$FOAM_SRC/transportModels/incompressible/viscosityModels/viscosityModel/viscosityModelNew.C。相关代码如下：

Foam::autoPtr<Foam::viscosityModel> Foam::viscosityModel::New
(
    const word& name,
    const dictionary& viscosityProperties,
    const volVectorField& U,
    const surfaceScalarField& phi
)
{
    const word modelType(viscosityProperties.lookup("transportModel"));

    Info<< "Selecting incompressible transport model " << modelType << endl;

    auto cstrIter = dictionaryConstructorTablePtr_->cfind(modelType);

    if (!cstrIter.found())
    {
        FatalErrorInFunction
            << "Unknown viscosityModel type "
            << modelType << nl << nl
            << "Valid viscosityModels :" << endl
            << dictionaryConstructorTablePtr_->sortedToc()
            << exit(FatalError);
    }

    return autoPtr<viscosityModel>
        (cstrIter()(name, viscosityProperties, U, phi));
}

此函数本质是在contant/transportProperties文件中查询关键字transportModel对应的值（以Newtonian为例），并在哈希表中查看是否存在这个Newtonian。返回Newtonian在哈希表中对应的值。以我目前的学习经验来看，这个对应值应该是类Newtonian的一个临时对象，并用cstrIter()(name, viscosityProperties, U, phi)括号里的参数进行了初始化。因此需要找到类Newtonian的构造函数。路径为：$FOAM_SRC/transportModels/incompressible/viscosityModels/Newtonian/Newtonian.C

Foam::viscosityModels::Newtonian::Newtonian
(
    const word& name,
    const dictionary& viscosityProperties,
    const volVectorField& U,
    const surfaceScalarField& phi
)
:
    viscosityModel(name, viscosityProperties, U, phi),
    nu0_("nu", dimViscosity, viscosityProperties_),
    nu_
    (
        IOobject
        (
            name,
            U_.time().timeName(),
            U_.db(),
            IOobject::NO_READ,
            IOobject::NO_WRITE
        ),
        U_.mesh(),
        nu0_
    )
{}

在构造类Newtonian时先构造了类viscosityModel。注意到类viscosityModel是类Newtonian的一个父类。因此需要找到类viscosityModel的构造函数。路径为：$FOAM_SRC/transportModels/incompressible/viscosityModels/viscosityModel/viscosityModel.C

Foam::viscosityModel::viscosityModel
(
    const word& name,
    const dictionary& viscosityProperties,
    const volVectorField& U,
    const surfaceScalarField& phi
)
:
    name_(name),
    viscosityProperties_(viscosityProperties),
    U_(U),
    phi_(phi)
{}

类viscosityModel的构造函数实际是将形参用于name_，viscosityProperties_，U_和phi_初始化了。这里的viscosityProperties_由于被类Newtonian继承，在类Newtonian的构造函数中可以看见nu0_使用了viscosityProperties_进行初始化。而viscosityProperties_的本质实际就是

 IOdictionary
    (
        IOobject
        (
            "transportProperties",
            U.time().constant(),
            U.db(),
            IOobject::MUST_READ_IF_MODIFIED,
            IOobject::NO_WRITE
        )
    ),

返回看nu0_的初始化：

nu0_("nu", dimViscosity, viscosityProperties_)

实际是读取constant/transportProperties文件里的nu关键字对应的dimensionedScalar。
而nu_

    nu_
    (
        IOobject
        (
            name,
            U_.time().timeName(),
            U_.db(),
            IOobject::NO_READ,
            IOobject::NO_WRITE
        ),
        U_.mesh(),
        nu0_
    )

nu_则是以nu0_这个dimensionedScalar构建了一个volScalarField。

此贴基本完结，按着源代码琢磨了一下午，基本弄清楚了是如何计算的。nearestPoint()主要使用了实时碰撞检测算法，可以计算出网格中心点到一个面的最短距离。这个最短距离是通过寻找面上离网格中心点的最近点来计算的。等我过段时间总结好了，到时候再补充到这个贴子下面。

在头文件找到了描述：The low-Re correction is activated by setting the entry lowReCorrection to 'on'; in this mode the model switches between laminar and turbulent functions based on the laminar-to-turbulent y+ value derived from kappa and E. When the lowReCorrection is inactive, the wall function operates in high-Re mode.可能是在constant文件里使用这epsilonWallFunction的时候能够选择打不打开LowRe

遇到一个以前没遇到过的工程问题，就是流场里面有一个没有门的房子，风可以从门洞吹进房子 里面。以前遇到的问题都是全封闭的房子，因此把房子内部的网格挖除，把房子边界面上的面单元设为wall就行。但是现在风可以吹进房子里面，也就是说房子内部也有网格，不能挖除，这个怎么画网格呢。

目前我的想法是，房子内外整体划分网格，把处于房子交界面上的面单元设置为wall，即红色的面单元是房子的屋面和墙壁，设置为wall，蓝色的面单元是门口，还是内部单元。不知道这样是否可行。也就是说房子内部和外部的单元是共节点的，交界面上的单元属性从内部单元，改成了wall

@Vortex 论文里一般设为symmetry

codeFixedValue

确实我看波浪用层流比较多。

以前用OF自带的interpolation，发现用并行也有点问题，当插值点在交接面的时候就会出错。不知道你这个是不是质量源在交接面时出了问题

奇怪了，虽然用起来没什么问题，但是打开的时候还是　会有warning

Invalid $PV_PLUGIN_PATH and -plugin-path= not defined
No supplementary ParaView/OpenFOAM reader modules
Using builtin reader: paraFoam -vtk

Created temporary 'platePimple.foam'

@xpqiu 还真是，我用paraview5.6.3可以完美显示出来了

链接文本，研究了一下用法，大家可以看看

@CFDngu 是不是写了什么代码没有在多核时注意reduce..

@sibo 并行限制2核，网格限制是52万