CFD中文网

应该注意到两个问题：1、volScalarField包含边界场 2、边界上没有体积
因此你应考虑Qc需不需要边界场，如果需要，应考虑边界值怎么计算（该公式无法直接用来计算边界值）；如果不需要，就不用定义为volSaclarField

所以你侧面给的固壁边界条件？这也不物理

@李东岳 多谢大佬点拨

@李东岳 看得懵懵懂懂，还是理不出来相应的数学表达式。感觉是一种线化处理策略，如果没有correction，就变成了纯显式处理

查看OpenFOAM10的代码可以发现，在src/ThermophysicalTransportModels/turbulence/eddyDiffusivity/eddyDiffusivity.C中的divq的实现为：

-correction(fvm::laplacian(this->alpha()*this->alphaEff(), he))
-fvc::laplacian(this->alpha()*this->kappaEff(), this->thermo().T());

而按我的理解，这里应该只是求解 $-\nabla(k_{eff}\nabla T)$，即：

-fvc::laplacian(this->alpha()*this->kappaEff(), this->thermo().T())

那这里多出来的correction应该怎么理解？

在v2012版本中，TEqn.H中有rhoCpPhi的定义为：

     const surfaceScalarField rhoCpPhi
    (
         "rhoCpPhi",
         rhoPhi*(Cp1 - Cp2) + phi*rho2*Cp2
     );

这个定义是如何得到的? 如果Cp1 = Cp2，那rhoCpPhi = phi*rho2*Cp2，这个和 phi*rho*Cp有差距。

如果直接把中心点的rhoCp插值到面上和phi相乘应该表示为 fvc::interpolate(rho1*alpha1*Cp1 + rho2*alpha2*Cp2)*phi。

这个问题跟setFields什么关系？
另外，这很可能是你续算操作本身的问题，应该将续算每一个步骤仔细说明一下

@cfdngu 问题没有描述清楚
1、有冲突是什么冲突
2、现在用的什么方法处理stl边界
3、尝试用toposet+createMesh是怎么尝试的，报错提示又是什么

@yongbo 这里的objectRegistry::template用来修饰 lookupObject<GeometricField>(name)，即在objectRegistry这个类里面有一个函数lookupObject，它是一个template

@wying 都看

仔细看yPlus的定义，它跟当地的流动状态相关

这是两个问题：
1、OpenFOAM中定义的张量，9个分量的顺序是 从左至右、从上到下
2、gradU本身的表达式不要想当然

个人原因，需要在苏州找份工作。但是网上找了一通，也没看到跟CFD相关性比较高的工作，不知有没有人推荐一个，或者告知相关公司、研究单位的名称？

本人在国内一自研CFD软件开发公司有6年的经验，熟悉C++、Fortran，主要侧重不可压流动，对两相流、Lagrangian粒子法有一定的研究，OpenFOAM也相对较为熟悉。

用学校的电脑重头开始算，然后续算看看。

（7）这个公式确实是错的，但是particle里面都用的是(11)这个公式，用来计算“位移”的barycentric坐标，而不是点的barycentric坐标，OpenFOAM/meshes/primitiveShapes/tetrahedron/tetrahedronI.H有计算点的barycentric坐标，你可以比较下。

把patches以外的内容，单独写成一个文件，然后include进来，应该就可以

不要抠字眼。。。 想想具体怎么离散的

试试这个，来源：https://events.hpc.grnet.gr/event/50/sessions/171/attachments/150/223/Paraview.pdf

f 表示 front，r表示rear
Cl(f)就是旋转轴前半部分对Cl的贡献
Cl(r)就是旋转轴后半部分对Cl的贡献

当然，类似的也可以定义Cd(f),Cd(l)

@沐风 没有什么特别的目的，就是写法问题，当然你也可以把它放到整个循环外面，只创建一次。

因为在右端源项里面含有未知数alpha_l，为保证系数矩阵对角占优，当源项大于0时，一般直接做显式处理，而源项小于0时，做隐式处理。

具体到这里，就是为了拆出小于0的隐式项和大于0的显式项。