CFD中文网

盲猜： 构建面是按照单元号从小到大遍历的，面方向为指向单元外侧。故面肯定是由小编号单元先建立的，指向小编号单元外侧（即从owner向neighbour）

@cfdngu 收费也没办法，我之前找了一圈也没有找到比这个好用的

ArmDDT软件，现在好像叫armforge

在OpenFOAM/primitives/Scalar/scalar/scalarFwd.H中：

//- A typedef for float
typedef float floatScalar;

//- A typedef for double
typedef double doubleScalar;


#if defined(WM_SP)

typedef floatScalar scalar;
typedef floatScalar solveScalar;

#elif defined(WM_SPDP)

typedef floatScalar scalar;
typedef doubleScalar solveScalar;

#elif defined(WM_DP)

typedef doubleScalar scalar;
typedef doubleScalar solveScalar;

#else
//    #error "PRECISION must be set to WM_SP, WM_SPDP or WM_DP"
#endif

好像是分内部点（所有面都不是边界）和边界点，内部点根据cell插值（距离倒数分之一），边界点根据face插值。

文件路径OpenFOAM-7/src/fvMotionSolver/fvMotionSolvers/velocity/laplacian/velocityLaplacianFvMotionSolver.C

Foam::tmp<Foam::pointField>
Foam::velocityLaplacianFvMotionSolver::curPoints() const
{
    volPointInterpolation::New(fvMesh_).interpolate     //就是这里
    (
        cellMotionU_,
        pointMotionU_
    );

    tmp<pointField> tcurPoints
    (
        fvMesh_.points()
      + fvMesh_.time().deltaTValue()*pointMotionU_.primitiveField()
    );

    twoDCorrectPoints(tcurPoints.ref());

    return tcurPoints;
}

有没有具体的数学公式，我看了一些书也都没说这些细节的东西。

我记得微信有个CFD招聘的公众号做的挺不错的

自己翻书找到答案了 ，《The FVM in CFD(OpenFOAM, Matlab)》那本书Page 303有讲到：

代码在src/finiteVolume/finiteVolume/gradSchemes/gaussGrad/gaussGrad.C144行：

template<class Type>
void Foam::fv::gaussGrad<Type>::correctBoundaryConditions
(
    const GeometricField<Type, fvPatchField, volMesh>& vsf,   
    GeometricField
    <
        typename outerProduct<vector, Type>::type, fvPatchField, volMesh
    >& gGrad
)
{
    typename GeometricField
    <
        typename outerProduct<vector, Type>::type, fvPatchField, volMesh
    >::Boundary& gGradbf = gGrad.boundaryFieldRef();

    forAll(vsf.boundaryField(), patchi)
    {
        if (!vsf.boundaryField()[patchi].coupled())
        {
            const vectorField n
            (
                vsf.mesh().Sf().boundaryField()[patchi]
              / vsf.mesh().magSf().boundaryField()[patchi]
            );

            gGradbf[patchi] += n *
            (
                vsf.boundaryField()[patchi].snGrad()
              - (n & gGradbf[patchi])
            );
        }
     }
}

里面vsf应该是待求梯度的场吧、然后gGradbf是根据体心梯度插值得到的面心梯度、n是面单位法向、.snGrad()是$\frac{\phi_{f}-\phi_{C}}{d_{Cf}}$ ?
这个面心梯度的修正公式怎么来的我还是没看懂：

gGradbf[patchi] += n *
(
    vsf.boundaryField()[patchi].snGrad()
  - (n & gGradbf[patchi])
);

这几天在看有限体积法的书，大佬的论文、《数值传热学》和 《The FVM inCFD OpenFOAM..Matlab》都看了，逛论坛发现都在推荐《Computational Methods for Fluid Dynamics》这本；去SpringerLink一搜正好今年出了新版，网上资源好像还不多，在这里分享一下： 点击下载

正好最近也有这个疑问，感谢大佬的回答～

@hy1112006 不会吧；我不太记得了，好像是用的https://openfoam.org/的OpenFOAM7，当初因为是租的服务器算，直接用sudo apt-get install openfoam安装的。

补个图片，tif格式好像不支持

可以，但非常复杂（难度不高，知道各点坐标就行），需要编程。推荐这个python库 ；

下面的例子是我用纯blockmesh画的网格，这个网格生成的程序在这个repo可以找到；
这种方法的优点是完全参数化，很方便调整网格。

其实很多人低估了blockmesh的作用，结合这个库使用之后我是觉得已经非常强大了。

把文件里某一个选项的值填为banana（任何一个错误的值都行），让错误信息输出该选项所有的正确可选值

@东岳 现在出了windows termianal，稍微配置下就很好看了~

@Cp_Zhao 找到原因了，创建实例公网带宽记得选按按固定带宽，然后适度调节。
虽然贵一点但至少用户体验好一些

我试了下，注册了账号充了100，然后到最后一步下载openfoam。。。网速只有几k/s，然后就懵逼了。

结论：

1. nut根据计算得来，所以设置成type calculated;然后随便给个值都行。
2. 不同模型间的公式尽量不混着瞎用，老老实实参考NASA里给的建议和https://turbmodels.larc.nasa.gov/ChangesToOpenFOAM.pdf的说明
3. Try not to follow the tutorial values as some (maybe most) of them are not validated per se but more for case file reference only（来自reddit yourstru1y的建议)

不知是否还有需要纠正的地方，感谢~

@东岳 我弄错了，确实稳态的话不会影响最后结果，差别很大是因为还没收敛(我只改了nut和nutTilda的值但没调endTime)。

但我发现合适的初值设置能够很明显的加快收敛速度：
如airFoil2D这个例子，默认是$\nu_{t}=\tilde{\nu}=0.14$，313步就收敛了；
而如果我改成$\nu_{t}=\tilde{\nu}=0$ , 需要371步，这都差别不大；
但如果改得大一点：如$\nu_{t}=\tilde{\nu}=14$ ,需要1018步收敛。

所以感觉有时候一个合理的初值还是挺重要的。