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做的求解器嘛，那个就可以，还能算凝固过程中的宏观偏析（化工里面好像叫分凝）

@羽之下 thank you，谢谢你的回复。我记得好像在哪看过说这里引用了一个假设，我决定把这部分系统的看一看，有什么新发现我会继续在这里交流~

更新：使用fvm的时候可以编译但计算不行，改成fvc 可以计算接下来就是看能不能收敛了

我想了解下OpenFOAM是否可以批量截取截面？比如轴向有1000个网格点，就截取1000个面，然后计算每个面上的某物理量平均值。我看cfdpost可以联合Python进行批处理。我尝试了foamDataToFluent，但导入fluent流场都错了，就做不了批处理。

谢谢呢，，已经解决了呀，，哈哈

1.JPG
src文件夹下没有fvoptions，但想利用其中的CellSetOption功能，有人知道吗？

@东岳 对，我那个5118是学校的计算集群，我记得我在哪儿看到过说这堆集群都是把超线程关了的，所以不可能使用的线程数量超过物理核心数量；然后睿频的话，应该是说如果你不关他，他默认是开的，但学校集群可能是关了的，你的如果是自己用应该是开着的，就是说如果睿频开着的，CPU温度和功率不超过他的限值的时候他可以运行在更高的频率上。比如你的5120 是14cores 28 threads 2.2G, 但是如果散热比较好功耗不超过某个限制的时候，他是可以整体跑在更高的频率上，这就能更快。
话说我觉得最惊奇的是我的9750H跑6线程居然比服务器CPU慢这么多，理论上说不过去啊，因为我看到跑的时候主频都跑到4G了，服务器才2.3G。得研究一下构架了

因为分散力模型无法找到，https://zhuanlan.zhihu.com/p/90057423

@东岳 好的，我试试