@东岳 谢谢，我是得系统学学。

不好判断，目前只能看出来算错了:zoule:

取决于你对CFD算法的理解，如果是新手的话，比如没怎么做过算法的，应该需要一个博士的流程。对算法中等理解，1-2年，对算法比较熟，1年。这是非常槽的一个估算

东岳老师，那界面不是应该附着在空洞的表面上么，还是说界面长度是指连续的油相前缘的长度，因为我看文献里界面长度都要远远大于这个几何模型的宽度。

dimensionedScalar nu ( "nu", dimViscosity, transportProperties ); wordList patchTypes(mesh.boundaryMesh().size(),"fixedValue"); volScalarField nu0 ( IOobject ( "nu0", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), mesh, nu, patchTypes ); const scalar nuWanted(1000.0); labelUList neiCells(nu0.boundaryFieldRef()[0].patch().faceCells()); forAll(neiCells, cellI) { nu0.primitiveFieldRef()[neiCells[cellI]]=nuWanted; }

@yhdthu 在 如何在domain内加边界条件？ 中说：

域内面积为0的面，然后加边界不？

很抽象，很难想象如何处理的呢

可压缩和不可压缩方程的数学特征的不同，导致需要调用不同的边界条件。比如标量传输方程，是一个双曲方程，如果进出口都给定Dirichlet边界条件就不能借。但对于扩散椭圆方程，如果进出口都给定Dirichlet边界条件很正常。

连续性方程 -> 压力方程，因此压力方程就是你说的这个方程

@dzsmoglai 在 tPEF湍流模型 中说：

oPhaseEulerFoam求解器的时候，总是得

感谢@东岳 老师的答复。目前此问题已经解决。原因也很简单。因为CGKE是基于kEpsilon的。而在PhaseCompressible湍流模型类中，在编译CGKE之前，事先通过PhaseCompressible这一类中的宏makRASModel(kEpsilon)编译了PhaseCompressible版本的kEpsilon模型。因此CGKE中任何函数在调用kEpsilon时就可以调用的到。本问题的解决办法就是同样的在makePhaseCompressibles.C中make一下RNGkE就行了，即makRASModel(RNGkEpsilon)。
如果有其它Foamer对上述解释不理解，我就再补充两句。kEpsilon, RNGkEpsilon这类湍流模型实际上是通用的模型，在incompressible，compressible, phaseCompressible等不同类的湍流模型类中有不同的体现。当然，这个涉及到OpenFoam的湍流模型代码结构，集成关系最多有8层左右，尚不了解的童鞋可以去参考此博客https://marinecfd.xyz/post/openfoam-templated-turbulence-modeling-framework-part-1/。链接是第一部分，一共三个部分。

建议参考这个
https://www.openfoam.com/documentation/guides/latest/doc/guide-turbulence-les-delta-vandriest.html
OpenFOAM里面的LESdelta是基于geometricDelta_的，geometricDelta_（delta关键字后面的那个，比如cubeRootVol）是纯粹与几何有关的一个量，而真正的LESdelta是建立在其之上的。这样的表述实际上更精确，所以才被设计成了两个东西，我猜的。

另外，我不清楚你这个画的是什么，只有一个类似的阶梯状的图。类似这个图我在今年8月份复现过，不过不是你这个算力，很多高阶格式都会把smooth的型线压缩成square的，例如Superbee，即使Minmod也有这种情况存在。

Figure 6.2(b) shows the same test problem as before but with the superbee method. The discontinuity stays considerably sharper.On the other hand, we see a tendency of the smoothhump to become steeper and squared off. This is sometimes a problem with superbee – by choosing the larger of the neighboring slopes it tends to steepen smooth transitions near inflection points.

Leveque 2002 p. 112

不过，我觉得你的问题可能不是这个...只不过图跟我遇到的有点像，给你点启发

看到个帖子不错，记录一下。关键词：householder transform

@Samuel-Tu DynamicField虽然叫Dynamic但声明的时候必须先给定大小，一般先选一个大一点的size初始化，用append添加元素，再用shrink删除多余的空元素。

@东岳 江湖百晓生，CFD万事通。

目前想在这个求解器中添加煤燃烧的反应，将利用反应速率计算需要添加的各个方程的源项。
1 目前的问题的我在src中的thermophysicalModels找到关于固相化学反应的东西，但是找不到有关化学反应速率的求解公式，例如固相阿雷尼乌斯求解反应速率，所以感觉没有地方能下手修改成自己需要的反应速率？
2 如果要直接在求解器中编写一个函数，例如是计算反应速率，我我需要调用到的其他数据（施密特数，直径，雷诺数）那些的需要用到什么代码（原谅我是一个C++沫沫）

在dynamicMeshDict中rigidBody一共有五种type：cuboid；masslessBody；rigidBody；sphere；jointBody，前四种都会用，唯独最后的jointBody没用过，查了一下资料没找着相关说明，有人知道这个jointBody怎么用吗？

只知道是按照从小到大进行排序，这个sort函数有什么作用？

void Foam::SortableList<T>::sort() { sortIndices(indices_); List<T> lst(this->size()); forAll(indices_, i) { lst[i] = this->operator[](indices_[i]); } List<T>::transfer(lst); }

@东岳 老师，这个公式（1）是他们参考另一篇文章
Wu, K., Chen, Z., Li, J., Li, X., Xu, J., & Dong, X. (2017). Wettability effect on nanoconfined water flow. Proceedings of the National Academy of Sciences, 114(13), 3358-3363. https://www.pnas.org/content/114/13/3358.short
我的理解是因为在纳米管内靠近壁面的位置，粘度变化没有统一的解析解，所以看到一些文章就这么分成“多层”粘度，算是粗略简化的方法吧 :xiezuoye:

这个有什么用啊？ cellset和cellset.toc()没有区别吧？

谢谢，知道了。