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

@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一搜正好今年出了新版，网上资源好像还不多，在这里分享一下： 点击下载

忽然想发个帖结果老用户没币，难受一波。