@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()没有区别吧？

谢谢，知道了。

https://www.cfd-online.com/Forums/openfoam-meshing/182460-creating-baffles-blockmesh.html

本人使用的是rigidBody，和sixDoFRigidBody不同的是，它有一项是transform，个人理解是用于从全局坐标系转到局部坐标系（动坐标系）的指令，所以如果动坐标系是固定在物体的重心处，transform中的向量项应该和centreOfMass是一致的才对，但是一算就死，后面试了下，把transform中的向量项写成（0 0 0）就能算，不过计算结果中会缺少部分自由度，有人遇到过这个问题吗？

感谢鸭蛋分享

其实我只是想实现tecplot里面的ddx函数，但是我找不到tecplot讨论区，所以先发这里了。
我得到的数据变量在储存在节点，而不是cell中心。网格是二维四边形网格。因为是非结构的，所以不能像笛卡尔网格那样直接求差商。我看了tecplot帮助文档里面ddx的实现方法是Moving Least-Squares，说是假定变量关于节点是二次函数，之后最小二乘法拟合出二次函数，然后将二次函数在节点的梯度近似为该变量在节点的梯度。我不明白对于二维的二次函数，tecplot是否考虑了交叉项xy?
如果不考虑的话，那我是否可以偷懒:假设需要求f的梯度，求对x的偏导的时候，完全不管y坐标。对单元的四个节点的x坐标x1，x2，x3，x4，和他们的函数值f1，f2，f3，f4进行一维的二次函数拟合。如果得到的拟合结果是f(x，y)=ax^2+bx+c，那第一个节点的x方向梯度就是把x1带入f'(x)=2ax+b中？

手机打字不太方便，请见谅。

Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: sixDoFRigidBodyMotion Selecting sixDoFSolver symplectic Translational constraint tensor (1 -0 -0 -0 1 -0 -0 -0 -624) Rotational constraint tensor (0 0 0 0 0 0 0 0 0) PIMPLE: No convergence criteria found PIMPLE: No corrector convergence criteria found Calculations will do 2 corrections PIMPLE: Operating solver in transient mode with 2 outer correctors Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Selecting laminar stress model Stokes No MRF models present No finite volume options present Constructing face velocity Uf Courant Number mean: 0.000625613 max: 0.10552 fieldAverage fieldAverage1: Starting averaging at time 0 forces forces1: Not including porosity effects forces forceCoeffs1: Not including porosity effects forceCoeffs forceCoeffs1: Not including porosity effects Starting time loop Courant Number mean: 0.000625613 max: 0.10552 deltaT = 0.005 Time = 0.005 PIMPLE: Iteration 1 forces forces: Not including porosity effects #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::sixDoFSolvers::symplectic::solve(bool, Foam::Vector<double> const&, Foam::Vector<double> const&, double, double) at ??:? #4 Foam::sixDoFRigidBodyMotion::update(bool, Foam::Vector<double> const&, Foam::Vector<double> const&, double, double) at ??:? #5 Foam::sixDoFRigidBodyMotionSolver::solve() at ??:? #6 Foam::motionSolver::newPoints() at ??:? #7 Foam::dynamicMotionSolverFvMesh::update() at ??:? #8 ? in "/home/z/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/pimpleFoam" #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #10 ? in "/home/z/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/pimpleFoam" 浮点数例外 (核心已转储)

@悟空空 我也是这样的问题，直接复制并行命令就出错，但输入并行命令就没问题，这是什么原因啊？

第一个。如果二维平面，墙面为x轴。无滑移边界条件，速度为0，说的是延x轴方向速度为0，没规定沿y轴方向。如果墙面沿y轴方向有个加速度。就不是第一类边界条件了，但还是无滑移的。一般墙不动，所以fixedvalue也没啥毛病。

第二个，可以找找fluent中symmetry和axis两种边界条件的区别。求解选项中axisymmetric和axisymmetric swirl的区别。比较绕，我还没弄清楚完，但大概和这个相关。

问题解决了，是因为我将Omega设置为零了，所以在计算的时，湍流模型在算湍流黏性系数是nut=k/Omega 出现除零错误。

@东岳 是的，这个operator[]不会每个时间步都去读，感觉这样改起来是个无底洞呢。不知是否能设置个贴体坐标，所有的点都相对其运动起来

没有指导意见，我现在用的sctoch划分，说的是划分后交界面最小，好像也没比simple快

@东岳 李老师下午好，我去看了simplifiedSiwek输出文件下拉格朗日文件里的d文件，我发现里面的颗粒粒径从反应开始到结束粒径都没有变化，但是颗粒反应过程中质量确实损失了，粒径应该要变小，所以想问一下李老师是不是因为coalChemistryFoam暂时还无法得出粒径的变化只有质量损失的变化？:xinxin:

@东岳 比较通用的就是多孔介质模型，好用好收敛，算的还快。之前不是在公众号上发过一个普度大学基于OpenFOAM做的求解器嘛，那个就可以，还能算凝固过程中的宏观偏析（化工里面好像叫分凝）