@星星星星晴 感觉nSector不能设为1，最后输出的sum(total mass)，sum(average mass flow rate)是nSector>1的4倍左右。

if (nSector_ > 1) { refDir = this->coeffDict().lookup("refDir"); refDir -= normal_[0]*(normal_[0] & refDir); refDir /= mag(refDir); } else { // set 4 quadrants for single sector cases nS = 4; vector tangent = vector::zero; scalar magTangent = 0.0; Random rnd(1234); while (magTangent < SMALL) { vector v = rnd.vector01(); tangent = v - (v & normal_[0])*normal_[0]; magTangent = mag(tangent); } refDir = tangent/magTangent; }

可能跟源代码中nS = 4;有关系

00d4f1b7-3cbe-4541-8fc1-e899bf8ec13a-image.png
代码路径：
channelIndex.H
channelIndexTemplates.C

想请教各位，他在channelIndex.H这里传的是const bool asymmetric = false，但是在模板实现channelIndexTemplates.C中又是const bool asymmetric。那么在初始化的时候，到底是用的哪个呢？

@李东岳 老师，DPMFoam的解析目前您是否已经整理完成了呢，在哪里可以学习到呢

想提取一下运动物体后面的速度，也就是探测点的位置也是随时间变化的，有没有大佬知道怎么做

请教一下研究过粒子的前辈，在openFoam里如何调用硬球模型？跟parcelTypeId有关么？自带的算例没找到相关的

有些网格面还不一定有网格点的定义 :chitang:

@香柏树 好的，多谢楼主分享

https://www.cfd-china.com/topic/2808 你可以去这里问问，我对这个turbulentDFSEMInlet还不太熟悉，那个帖子大佬很多

目前计算一个喷雾燃烧场时，想改变下燃油喷射位置，因此把几何模型的喷嘴结构(5mm圆孔)向上游移动了10mm, map前在system文件夹下创建了mapFieldsDict文件，具体信息如图1所示；然后在对应时刻文件夹下建立了和原始计算场相同时刻的文件（共125个）。在map时使用了1）mapFields /OF6/P25；2）mapFields /OF6/P25 -sourceTime 0.116555 -mapMethod mapNearest -targetRegion region0两种方式均出现phi未能成功map,同时PaSR_kappa和PaSR_Qdot，sprayCloud类文件除了sprayCloudTheta文件，其余均未能map过来，如图2所示。另外在现有map结果基础上试算了下，log和err文件如图3所示，在读取化学反应机理时浮点溢出了，尝试了foamChemistryReader和new format均未能解决。请问这俩问题大家有什么方法解决么？万分感谢！
771f0b91-f3ee-41db-83b0-01decf0d55ac-图片.png
f30a09d4-20a4-45cc-9750-171e698d5778-图片.png

@fangyuanaza 我也遇到这个问题了，貌似是工具里面判定网格中心点距离的问题，太小的网格就提取不到值

你把斜坡做成STL，用snappyHexMesh抠，应该就行了

@myheart 懂的懂的。我这面也经常是，一个bug卡好几天，有个都是几个月，还有按年算的。。。就比如你这个回帖，3行搞定一个问题，其实背后肯定debug好久。辛苦辛苦。

@李东岳 我现在没有好的算法，知识发现零梯度这种边界条件是可以使得场内的浓度出去，不会在边界处积攒，如果编写算法，应该是和零梯度类似的可以。

@李东岳 好的，谢谢李老师

微信图片_20220823200056.jpg 微信图片_20220823200145.jpg 微信图片_20220823200149.jpg
自适应网格达到设置最大的网格数后面的计算就不加密了，
但是我想实现的是网格在压力波传播的附件加密，压力波传播过去后加密过的网格恢复原始大小
请大佬指教是可以怎么解决
（我用的blastFoam求解器）

@李东岳 谢谢李老师！

@星星星星晴 非常感谢！我找找看

@星星星星晴 感谢大佬指点！

看起来像个bug。你可以试试双核的，然后手动把traslationRotation写进processor里面，然后看看报错不。不报错的话可能是bug

@李东岳 谢谢李老师，我再看看