:chitang:

好奇老师最后问题是怎么解决的？今天我也遇到了类似问题？

@韬智tz 请问您找到ccm向OF9网格转换的方法了嘛

求推荐OpenFOAM编程方面的资料或课程。
最好是整个流程讲解得比较全面的。

目前想实现自定义格式的流场数据输出。后续打算开发自定义dsmcFoam求解器。

已解决，个人发现在controldict中的 “ writeformat ” 为binary，更改为ascii即可

请问下东岳老师还有各位大佬，icoFOAM的qDNS案例是可以直接运行的吗，想要作\(u^+\)-\(y^+\)图，参考Calf.Z-DNS的回答eb0a816e-e80f-4390-902e-fb6da872db51-image.png 以及
https://www.cfd-online.com/Wiki/Y_plus_wall_distance_estimation的公式，
算了1000s后结果如图
dfacee60-97d4-45d2-8bf8-55dd702b48c8-image.png
发现几个问题：
1.case给的是\(Re_\tau\)=180，所给文献的数据也是\(Re_\tau\)=180，但是transportProperties文件给的是运动粘度nu的值为2e-5，注释却是\(Re_\tau\)=395，按照2e-5计算后\(Re_\tau=\frac{0.1335*0.001332}{0.00002}\)=8.8911，不知道雷诺数以哪一个为准；
8436cb60-fdbc-408b-80d9-f5919d41eba0-image.png
2.做出的\(u^+\)-\(y^+\)图取了\(y^+\)小于180的部分，大概距离边界层0.2，\(u^+=\frac{u_{mean}}{u_*}\)，其中\(u_{mean}\)的值为第1000s的平均值，提取中线上x方向的速度，位置如图
1a511cf7-4222-4bf0-a0b2-a8993b054d12-image.png
计算后得到\(u^+\)的值偏小，结果如图
a52d1d15-c00c-4bca-ac0b-9a702045aa02-image.png
希望大家给予解答，谢谢。

@shrine 你多试几个，我也是新手，之前遇到过类似的，换几个总有能用的

@number44 好的，已经按照这个思路在做了

@bestucan老师我这个算发散了，后续算一个正确的再上传！

@李东岳 东岳老师，我了解了以下结构化网格的排序，然后刚才通过postProcess -func writeCellCentres 写出了channel395的网格，这里网格是按照先x在y后z的顺序排列的，时间步中的结果应该是和这个对应的请问下东岳老师，有办法把结果输出成三维列表来么？

@bestucan 谢谢前辈，版本是一样的，后来和学校超算人员联系发现是少了一个环境配置加载：module load gcc/7.4.4，加上以后就可以了

@冠竹 按道理来说不应该啊 你这个都属于是神bug级别了 如果有时间你就多debug一下 确定是bug可以去官方提bug 没时间就算了这样吧

三台电脑已经通过网线建立了小集群

最开始看到这个话给我整蒙了，我说3个电脑连起来不通过交换机咋搞呢？

就是普通的网线，水晶头的那种，都连到普通交换机上，网卡什么级别我不太清楚，他们的型号都一样

我仔细看了一下，原来楼主通过交换机连起来的。

各位大神请问一下，我最近在模拟气泡与自由液面之间的影响，用的Depthcharge2d改编的，我想问我设置了初始压力和初始半径，还定了气泡位置，怎么得到最大气泡半径呀，还有有没有什么什么办法知道计算结果中的每个时间步的射流速度和气泡半径，像下面图所示。
微信图片_20220523110946.jpg 微信图片_20220523111112.jpg

@fubianhanshu
在particleIO.C下添加一行代码，注释原先的输出可以输出颗粒全局坐标。

os << position() << token::SPACE << celli_;

先采用最笨的方法，把保存时间步设置小，遍历读取坐标，其他的再慢慢实现，看着lagrangian粒子程序太复杂了:136:

OpenFOAM紧靠壁面网格结点上的物理量（比如湍动能）计算值怎么用代码提取出来？

请问老师前辈们，我在网格文件中设置的matchTolerance数值为1
efa9a684-0340-4883-86d8-734c6180ab43-b812433ad0b8c91dae5ea5ced31c96b.png
但是最后报错始终显示为0.001导致的网格问题
a9c5ebd8-72c9-4e43-9225-c7ff5fd797bc-1653209642(1).png
是什么原因呢，请问要如何处理。

@浪迹天大 请问下这个怎么使用啊

采用PisoFoam计算得到的水翼升阻力曲线不稳定，怎么回事
0f5142e4-4ca8-4d51-9100-f44de4d62bf8-屏幕截图 2022-05-22 110540.jpg

@李东岳 @bestucan
感谢两位老师，我大概明白了，马上就去详细看一看:xiexie:

求助各位，目前of的动网格技术可不可以实现下图功能：在计算域顶端和计算域内的物块之间连接弹簧，同时给整个计算域施加简谐运动？
9dbaecc8-35ff-49cb-a874-3db65bdce331-1653113557(1).png
目前学习of五个月有余，从我的了解来看，of自带动网格库中的弹簧只支持一端连接固定点，一端连接运动的物体。想问问大家有没有做过类似的东西，能提供些思路，不动底层代码能够解决这个问题，不胜感激。

paraview 的操作可以写 python 脚本。直接写比较麻烦，可以宏录制之后复制粘贴改参数

tools -> start trace

@coolhhh 感谢回复！

@李东岳 李老师，我用的是extend中的fsifoam求解器，在open foam中画出网格用extend计算，会这样报错。
433cf721-c2ed-4e1f-8490-c3cbeba84ec8-ac9b3169084b49719ca8976178b2dc4.png
但是我的几何边界中没有命名global的边界。请问这是什么问题呢。

@coolhhh 厉害厉害，我就觉着是变量定义格式问题，但是自己不太懂。感谢～

@yfleung 请问下用paraview该怎么做流向和展向的平均呢？

请教一下各位，我做的圆环柱的传热，流向为周期性边界条件，内壁面恒定热流，外壁面绝热，内外壁面无滑移，fvOptions加了一个meanVelocityForce，（0.7 0 0），0.7就是流向，可是最终的速度云图，产生了一块速度较大的区域，并且该区域还会沿环向流动

如图，可以明显看到红色高速区在缓慢的逆时针移动，请问下为什么会产生这个高速区，以及为什么高速区会移动？

压力与速度是耦合的 重力是个源项啊 没法取梯度吧

@zhe :146: :146:

@hhh 解决了，需要现在ICEM中设置周期性边界，之后就可以了
我参照的是https://blog.csdn.net/muyangzixue/article/details/108631943

@tens 感谢您的提醒，我再去看看

@星星星星晴 感觉还是不太对，网格需要再细一些。

把 return 挨个屏蔽，看看屏蔽了哪个报错。。。:yes:
就知道调用哪个了

我想知道它是什么意思呢？

好像是用 u 和 dim 初始化了一个新的 fvVectorMatrix 类型的变量

@李东岳 东岳老师您好，感谢您的回复。我重新安装gcc 9.4.0和gsl 2.7之后，使用root账户尝试编译发现错误，如下：

make[1]: *** [/home/user1/OpenFOAM/user1-v2012/platforms/linux64GccDPInt32Opt/bin/faceSetToSTL] Error 1 make[1]: Leaving directory `/home/user1/OpenFOAM/user1-v2012/applications/utilities/waves2Foam/applications/utilities/preProcessing/faceSetToSTL' make: *** [faceSetToSTL] Error 2

想询问一下这种错误是来自 gsl 版本问题吗？或者是否需要重新编译OpenFOAM-v2012（之前使用的gcc 4.8.5 编译的v2012）？

在高版本Paraview中我通常是采用contour读取界面，然后选中该表面save data，进行输出界面上的各种数据，但是由于实验室的服务器过于老旧，只能使用paraview4.4版本，同样的操作在输出数据时，得到的数据文件及其繁琐乱糟糟的，一个时刻下竟然生出了5个文件，请教一下大家有没有更好的方法简化输出数据的方法

你有热通量输入到计算域里，但是没有边界让输入的热量离开，系统的温度无法达到统计稳态那就自然随着时间不断提高了

@vien锁 您好，dsmcFoam输出哪些流场文件是由什么控制的。注意到像icoFoam这种求解器是通过fvSolution来设置的，但是dsmcFoam好像并不是。

而且我在使用dsmcFoamPlus时发现，提供的几个案例输出的流场参数也是不一样的。按理说相同求解器输出的哪些流场文件应该是相同的，除非设置了字典文件，然而我在字典文件中并未发现这种设置。

@hellofoam 请问一下,自动画图的温度范围区间不合适，是300-1200，我想把他调整到800-1200该怎么整？

@zhanghan 您好，我最近在用buoyantPimpleFoam做管道的传热，可是温度分布云图却几乎和层流的一样，请问您方便把您之前的算例文件借我借鉴借鉴么？或者说一下您的温度是否出现过我这种问题以及如何解决的？

都不可压了，计算中密度也不会变。可以看看东岳流体对基本求解器的解析。

密度是设定的，在 constant/transportProperties 里

@vbcwl 米

@bestucan 好的，非常感谢大佬的指点！

大家好，我想请教一下KinematicParcelTrackingDataI.H中的TrackCloudType& cloud里的cloud指的是什么？在如下代码中：

rhoInterp_ ( interpolation<scalar>::New ( cloud.solution().interpolationSchemes(), cloud.rho() ) ),

cloud直接调用了函数rho()，为了弄清楚cloud是哪个类，我改变了rho()为任意的字母，比如a()，但是重新编译依然可以编译成功，这是为什么呢?

@cccrrryyy 老师，我这个结果跑完。后续贴上，麻烦！

@hotairballon1997 再coneinjection中 你可以通过控制flowrate来控制你喷射的parcel的速度

case ftFlowRateAndDischarge: { const scalar A = 0.25*pi*(sqr(dOuter_) - sqr(dInner_)); const scalar massFlowRate = this->massTotal()*flowRateProfile_.value(t)/this->volumeTotal(); const scalar Umag = massFlowRate/(parcel.rho()*Cd_.value(t)*A); parcel.U() = Umag*dirVec; break; }

@zhy2022 在 controlDict文件中关于自动调节库朗数 中说：

@hitsc30 稳定状态感觉没什么影响，但是瞬态的话影响很大。我一般都还是保持在1以下。个人愚见库朗数30+还不如把网格画粗一点。不过结果对得上最重要:146:

目前来看确实是这样。不过网格很粗的情况下，我的y+估计就有点大了。一般要是壁面的网格小于1的情况下，保证库郎数小于1的话，感觉时间步小到无法接受。计算成本太高了。网格画的厚一点的代价就是近壁面的y+估计要上升了。不知道怎么权衡了
蚌埠住了，我做的瞬态模拟

@hitsc30 自抛自扣一下，我太蠢了。中间缺少一个文件夹，细心细心再细心

@number44 好的