问题已解决。

发现原因可能是points文件需要的是一条入口inlet面的边上的线。
比如下图，我本来用的是线CD上所有cell的坐标，然后报错如主楼所说，我稍微看了一下代码，发现 if (!fullMatched) 就会输出错误信息如上。

然后我试了试把points中的数据换成AB线上的cell坐标就可以了。

具体原因本工程专业的代码小白现在先不深究了，等有空了再看看。

微信图片_20210924201829.jpg

额 好像发现个问题，质量守恒方程应该是这个
7973d26e-efef-4af5-9e55-686a447c4996-image.png

其中𝑆(𝜌)与 𝑆(𝛾)形式一样，只不过自变量由相体积分数变成了密度，不知道这样对不对？

感觉有点弄乱了，求大佬指点，不胜感激:135:

@wwzhao 教授，您好！能给予鄙人一些建设性意见吗！感激不尽

@hongjiewang 在 OpenFOAM小代码 中说：

@李东岳 目前我需要将温度边界和压力边界定义为随相含量变化的一个场。我只可以做到变化之后的也是同一类边界条件，但是像下图中的气相区压力是第一类边界条件，液相区温度是第二类边界条件。请问这种变边界类型的边界应该怎么定义~1632279639(1).png

已解决压力的~但是温度的第二类和第三类的转换还没有解决~

9818ab58-ba59-4644-af8a-2b91924553cc-image.png
模拟环形截面流道内往复振荡流动，边界条件设置为正弦速度入口,余弦压力出口。理论及文献中曲线在180度时分界，而我的模拟结果发现曲线随周期呈整体右移趋势，在15周期分界点移动到192度。请教各位老师这是什么原因？

@gemini 嗯，上面那个图就是设置的oil相0.5，water相0.5，水的体积分数云图显示出来就是上面那个图的底部蓝色变为了红色，其余为蓝色。现在就是知道计算的时候怎么区分各相，因为有两相的颜色是一样的

@小考拉 自问自答来解决这个问题。

首先定义一个全为0的流场命名为c_flow

然后利用下面的命令

c_flow.boundaryField()[patchID] = c_flow_origin.boundaryField()[patchID]

最后输出c_flow流场即可:chitang:

大家好！

最近我在OpenFOAM中使用fvOptions设置了一个动量方程的源项，定义为
\begin{equation}
S=-C_z \times |U| \times U
\end{equation}
$C_z$是一个常数，$|U|$是速度的大小，$U$是速度，我在fvOptions中是直接写了

eqn += -Cz*mag(U)

我看到源项的表达式应该是写成，是乘以网格体积了的
\begin{equation}
S=S_U \times V_P + S_P \times\ V_P \times T_P
\end{equation}
请问我在fvOptions中写的方式是需要乘以网格的体积吗

ANSYA进行双向流固耦合后，在流场当中给定固体加载udf后，在最后的系统耦合界面中只有一条数据转换，也就是力与力的关系，并不是力与力，位移与位移的关系，这样操作无法跑出来，如何才能正确的设置，请教各位师兄，老师！麻烦大家了！（最简单的一个说法就是某块平板在某流场中按照预定的规律运动，想要的结果是平板弦向的压力云图，不可忽略平板对流场的作用力）

求助模拟大神给看下，这个case中结晶设置是否有问题？想模拟液固流化床氯化铵冷却结晶过程，不知道传质那部分怎么设置？没找到帮助文件中mass mechanism 下拉列表中的population balance机制

@jqbu 好兄弟，of 9不就推出这个功能了吗？你肯定是漏掉东岳老师的公众号推送了，去李老师公众号看:chigua:

@鲸落 请问你解决了吗~动态读取cellZone~

@wying 不好意思明天我就要去杭州了... 国庆后回来 只能看看有没有其他大佬能给你看看了，:135:

@李东岳 谢谢谢谢李老师！！这个是我20年的暑假整的，但是当时没整出来，今年才发了小论文，下午注册的账号，我可太激动了！马上就发了一贴。
补充一下就是气泡那个模拟，两个求解器是用的同一套网格。

众所周知，当我们在进行并行计算的时候，需要把网格进行剖分，OF中可以使用第三方库进行区域分解，可是我们知道，OF中的区域分解是不带重叠层的（或者说没有halo单元），也就是，同一个网格仅可能属于一个处理器。
那么问题来了，我现在想实现带有halo层的重叠型的区域分解，我该怎么实现呢？或者说，有没有那位大佬做过类似的工作或者相应的资料吗？跪谢了！

@hongjiewang 在 温度边界条件 中说：

1632279639(1).png
气液界面是不断变化的，请问红色方框内的边界条件如何添加~

可以识别出气液界面，但是想知道如何在同一个边界里使用第二类和第三类边界条件

@prometheus10 然后，文中提到的另一个链接:papa:
https://www.cnblogs.com/liusuanyatong/p/11259653.html

@马乔 前辈你好，谢谢您的回复，也祝您中秋快乐。

看了您的回复，我的理解是，这个DFSEM算法，是根据我输入的R、L文件的值来计算出脉动速度（也就是您说的波动速度），然后加上U文件的值来生成相应的风场。

我看了tutorials中的case后，现在还有一些细节的地方不知道怎么理解，麻烦您指导一下：
1、对于U、R、L和point四个文件中数据的个数，可以随便定义若干个点就行？比如我inlet面(y×z)有5000个点，我就挑y方向一列800个点的数据，来输入对应的四个文件就行？

2、这个U文件应该是输入的平均风速吧？比如我想要生成的流场是平均速度为（1 0 0），波动速度为0.5，那么U文件中所有的数据都应该是(1 0 0)对么？

3、相应的R和L文件，雷诺应力R应该就是[ui'uj']的一个对称张量，湍流尺度L应该怎么定义呢，看case中是一个列向量，那假如一个三维各向异性的湍流场是否需要输入一个n行3列的数据，如

{ L11 L22 L33 L11' L22' L33' L11'' L22'' L33'' L11''' L22''' L33''' ...... }

再次感谢您的回复，希望您空闲的时候能再指点一二，谢谢！

请教各位同学和老师，我在尝试用MPPIC在一个存在斜面的三维复杂网格中模拟固液两相流时，计算一直没办法收敛。尝试不同的边界条件和不同的k和epsilon值，都没有任何改善。以下是stack trace：

Lookup solver for epsilon.waterFinal [16] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [16] #1 Foam::sigFpe::sigHandler(int) at ??:? [16] #2 ? in "/lib64/libc.so.6" [16] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:? [16] #4 Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > Foam::operator*<Foam::volMesh>(Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > const&, Foam::DimensionedField<double, Foam::volMesh> const&) at ??:? [16] #5 Foam::RASModels::kEpsilon<Foam::PhaseIncompressibleTurbulenceModel<Foam::singlePhaseTransportModel> >::correct() at ??:? [16] #6 ? at ??:? [16] #7 __libc_start_main in "/lib64/libc.so.6" [16] #8 ? at ??:? [cas019:410749:0:410749] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2)) ==== backtrace ==== 0 /lustre/opt/cascadelake/linux-centos7-skylake_avx512/gcc-8.3.0/ucx-1.6.0-dacxhp7gepme2j7jmvpcmbm2qrptc6k7/lib/libucs.so.0(+0x1bf44) [0x2b968ef64f44] 1 /lustre/opt/cascadelake/linux-centos7-skylake_avx512/gcc-8.3.0/ucx-1.6.0-dacxhp7gepme2j7jmvpcmbm2qrptc6k7/lib/libucs.so.0(+0x1c2ca) [0x2b968ef652ca] 2 /lib64/libc.so.6(gsignal+0x37) [0x2b9689945337] 3 /lustre/home/acct-naozzn/naozzn-user1/openfoam/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6_+0xd5) [0x2b9688a0dca5] 4 /lustre/home/acct-naozzn/naozzn-user1/openfoam/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4FoammlINS_7volMeshEEENS_3tmpINS_16DimensionedFieldIdT_EEEERKS6_RKS5_+0x164) [0x2b9683ba5d74] 5 /lustre/home/acct-naozzn/naozzn-user1/openfoam/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libDPMTurbulenceModels.so(_ZN4Foam9RASModels8kEpsilonINS_34PhaseIncompressibleTurbulenceModelINS_25singlePhaseTransportModelEEEE7correctEv+0x547) [0x2b9687558727] 6 MPPICFoam() [0x433d76] 7 /lib64/libc.so.6(__libc_start_main+0xf5) [0x2b9689931505] 8 MPPICFoam() [0x4352e0]

边界条件如下表

f3ef028e-2057-43c0-8e87-fe42f54246e7-image.png

感觉应该是颗粒太大或者是网格划的不太好，有没有大佬能帮我看看怎么调整更准确一些

我用tec读openfoam数据二进制格式无法读取，你要输出十进制。而且就算十进制，一些复杂网格的openfoam数据我用tecplot读进来就闪退，所以还是用paraview了

@tens 好的，谢谢你，我试一下

各位老师，大家好，想请教一下，在拥有管理员账号的前提下，怎样在联想集群的LSF作业调度系统中添加新的普通账户？
之前是按照这样的方式，也忘记了是在哪里找到的了：:135:

//在管理节点上添加新用户(root)： useradd -d /gpfsdata/home/cfdem -m cfdem //在管理节点上给新用户添加密码(root) passwd cfdem //添加到用户组 vim /gpfsdata/apps/lsf/conf/lsbatch/cluster1/configdir/lsb.users //添加完后运行 badmin ckconfig //在计算节点上添加新用户 ssh c01n01 useradd -m cfdem passwd cfdem并输入密码 //同理在c01n02~c01n03同样上述操作 //ssh到管理节点mgt01并切换到cfdem ssh mgt01 su cfdem ssh-keygen//然后一直回车

但是这种方式感觉有点繁琐。:136:
后来了解到PBS管理系统里面可以直接通过命令

clusconf -au 账号名

就可以完成普通账户的添加。
现在想问一下，LSF管理系统添加账户的时候，有没有类似于PBS管理系统的这种一气呵成的简单命令？:146:

各位大佬，我想在一个长宽高为2的domain的内部设置一个面作为inlet，其他的6个面均为wall。。

看了好多topoSet，createPatch，CreateBaffle的教程，还是没整明白。。
求指点到底该怎么做。。

blockmesh如下，
在y=0.5的地方建立一个，x方向长0.5 y方向长0.5 的矩形平面作为inlet

convertToMeters 0.001; vertices ( (-1 0 -1) (1 0 -1) (1 1 -1) (-1 1 -1) (-1 0 1) (1 0 1) (1 1 1) (-1 1 1) );

@李东岳 感谢李老师的交流平台，希望以后有机会向学习OPENfoam，加深对模拟的学习:high: :high: :high:

2021-09-17 15-46-22屏幕截图.png

@冠竹 这个是$\overline{u'u'}$, $\overline{u'v'}$等一共6个，就是解析的雷诺应力除以密度

@五好青年 你可以参考一下Holzmann CFD官网的算例, 关于Magnus Effect的算例中有加入圆柱转动函数条件的相关设置.

专业：流体力学、传热学相关专业
主要职责：用C++语言，进行CUDA和GPU计算开发CFD计算软件。负责流动、燃烧及传热等仿真软件的开发及封装；
工作地点：长沙
公司性质：央企中国电子信息集团所属子公司，协同NUDT研发
薪酬待遇：本科年薪15W左右，硕士年薪20W左右，博士25W左右
其中一个招聘链接
https://www.zhipin.com/job_detail/fe4f05de98c982591nN62Ni7GVZT.html?ka=hot-job-1

@liujuncfd
本科待遇年薪15w左右，硕士年薪20w左右，博士年薪25w。
https://www.zhipin.com/job_detail/fe4f05de98c982591nN62Ni7GVZT.html?ka=hot-job-1

目前我只会做到在边界某段添加固定温度，整个边界都使用热流密度条件也是可以的，请问如何在边界的某段添加热流密度边界~比如下面这段0.000333<=x<=0.000666使用热流边界 应该如何修改~

OUTLET { type codedFixedValue; name heated; value uniform 300; codeInclude #{ #include "fvCFD.H" #}; codeOptions #{ -I$(LIB_SRC)/finiteVolume/lnInclude \ -I$(LIB_SRC)/meshTools/lnInclude #}; //libs needed to visualize BC in paraview codeLibs #{ -lmeshTools \ -lfiniteVolume #}; code #{ /*const IOdictionary& d = static_cast<const IOdictionary&> ( dict.parent().parent() ); const fvMesh& mesh = refCast<const fvMesh>(d.db()); const label id = mesh.boundary().findPatchID("OUTLET"); const fvPatch& patch = mesh.boundary()[id];*/ scalarField T(patch().size(), 300); forAll(T, i) { const scalar x = patch().Cf()[i][0]; //const scalar y = patch.Cf()[i][1]; //const scalar z = patch.Cf()[i][2]; if(x>=0.000333||x<=0.000666) { T[i] = 400; } else { T[i] = 300; } } (*this)==T;//T.writeEntry("", os); #}; }

@队长别开枪 谢谢老师指点，我好好研究一下

@七辆战车
你好！ 请问一下，利用rigidbody求解器，有实现双浮体吗，是直接在dynamicdict文件下给两个浮体的参数吗？ 对于外部导入的水平圆柱浮体，刚体类型可以给rigidbody吗？ 期待你的回答，不胜感激。

@gemini 哦哦，好的，谢谢您

@队长别开枪 谢谢老师您的指点，已经安装上了。但我是用hypremesh画的斜坡网格，是否可以不使用stl文件来进行相初始化。

我看了下这个文章。这个文章非常适合练手，挺有意思。通过植入3个传输方程，处理传热流动。openfoam植入非常straightforward。做个标记。或许下次可以当做课程资料。Fluent那面不清楚了。

你为什么不直接用传热流而是一定要自己植入一个呢？

@袁宝强 啊我知道了 谢谢 :xiexie:

:134:

标准的限制方法是保证插值后网格顶点（vertice）的值不会出现极大值或极小值（比周围网格的最大值大或比周围网格的最小值小）因为梯度插值最容易出现越界值的就是在顶点处，从而用一个标量来scale梯度gradient；另外一种方法是保证网格的面心midponts of cell faces处不出现极大值或极小值即可，就认为满足了插值后不存在越界的问题。我想的是关于第二种方法，面心处插值后不出现最大值或最小值不能代表顶点处不出现最值，对吗

感谢大家 已更新新数据

@李东岳 老师 ，侧面我是没有进行面关联的只是进行了线关联，vertex设置为周期性节点，请问老师取消关联是将edge取消关联吗？

我旋转网格的时候，勾选了merge nodes，保证节点能相互传输数据。因为我把侧面定义成了一个part，我想生成全体网格后再删除这个part，也就删除了面以及面网格，不知道这样行不行？

各位大佬，
我最近在做一个温度场模拟，模型简单来说：是一个密闭空间，空间外部是自然环境，空间内部是真空，内部有一固体热源，热量通过真空辐射传热至封闭空腔的壁面，再经壁面对流换热传至外界。
我的问题是，内部空腔的真空区域怎么设定材料属性？网格导入后，真空区域默认为流体区域，设置空气材料显然是不正确的。
另外，若该真空区域未达到理想的完全真空，真空度为90%(可认为是存在稀薄的空气)，这种情况又该如何模拟？
望答复，祝好！

@tens 谢谢啦

@李东岳 谢谢李老师

@jinlinna 感谢！

@lts您好阿，最近我在使用 simpleFoam 时也遇到了同样的问题，请问您是如何解决的。谢谢！

管内气体流动，管壁定温加热，我在计算的时候设置了report：热流密度，对流换热系数，Nu等参数。发现用热流密度反推的对流换热系数和Nu和报告的对不上，相差很大，请教各位老师这是什么原因？谢谢

之前接触过CFD++和STAR CCM+的噪声计算，里面都有专门的模块。openfoam这方面的功能没试过，不过搜到了一个第三方库，里面有源码和一篇文章，你可以参考下：
https://github.com/unicfdlab/libAcoustics
实在不行就把声源面的压力脉动都导出来，按定义写个程序去积分，虽然麻烦了点，但应该也是能做的

@jacob-z 请问，您找到正确处理办法了吗？

@李东岳 我觉得李博整个一套下来思路蛮清晰的，先在学校搞几年基础，再青椒几年体验一下教学相长，然后创业做项目咨询吸收行业应用经验，现在响应号召开始自主研发，祝成功！