也行也行

@东岳 谢谢老师，确实是这样，setFields从新设置就好了。太感谢老师了。

@r07522121 你最好看一下模型的.H，看看都有什么东西可以定义

有做爆轰的大佬吗，请教如何得到处理出如下的胞格图QQ截图20200519151140.png

请问这个滤波为什么会报错有什么结果吗? 我现在也碰到报错了，但是其他mesh又可以用，在怀疑是不是mesh的问题？

@Jato 我的是在OpenFOAM-2.4.0版本下编译的，没有出现你这个问题

想请教下湍流模型中0文件夹下参数的取值是如何选择的呢，比如k, epsilon, omega， nut等，谢谢！

大家好，请教一个问题，我想用buoyantBoussinesqSimpleFoam求解器求解层流问题。将湍流参数字典设置为层流，我发现仍然需要在transportProperties字典中提供湍流Prt数？这里的Prt数参与计算吗？
同时，0字典中，也仍需要提供alphat文件。

感谢分享

@东岳 老师你好，非常抱歉回复晚了。我以前只尝试过cavity、dambreak、sloshingTank2D、flange、sloshingTank3D6DoF这几个算例，其余的算例还没试过，除了cavity以外都会出现计算很慢然后莫名其妙停止的情况。现在更改为双精度浮点数之后计算效率与Linux基本差不多。

研究了下manualInjection的构造函数，直径和速度分别有构造函数，构造了List，说明可以分别设置，有了更系统的发现再跟大家汇报

更正：LES model没有kOmegaSSTIDDES，应该是kOmegaSSTDES，然后LESdelta type应该使用cubeRootVol

思路：如果是内部场的值，调取a最大值的cell label，然后拿出这个label的b值就可以了。对于边界场同理。

@CFDngu 而且forAll操作可以对空容器进行吧。。实在不行你可以给每个空容器加入第一个数-1，而网格编号不可能小于0，你循环的时候判断下网格编号是不是小于0就行了

楼主做波浪的啊。我最近也读了一点文章，准备开始interFoam，想加一点浮冰啊，空气摩擦的影响进去，感觉可以交流一下。不过我看有些文章里面老是提SWAN，对这个你有了解吗

很厉害很厉害！！

就是全场吧， 因为1-alpha1 是对于水为0， 所以

gsum（（1-alpha1） × mesh.V()）

就是bubble的volume；下面的那个速度也是一样的， 但是那个速度只是一个平均速度。

@hurricane007 脑残了，之前写错地方了。现在找到了

感谢分享！

好的，谢谢东岳老师指点！

@东岳 感谢老师回复，耦合求解器后面放弃了😂

目前想做一个俯仰翼型的算例，使用dynamicMotionSolverFvMesh网格运动求解器和displacementLaplacian动网格求解器。看到一个圆柱振荡的算例，想跟着这个算例进行修改，振荡算例的pointDisplacement文件的边界条件如图：69ecd26a-81fe-4165-9daa-fcad778d8ba4-image.png 如何把type进行修改，变为转动的类型，对应需要的参数又有哪些，请大神指点一下

@东岳 :140:

@东岳 好的

我在计算圆柱涡激振动，想要设置随时间变化的弹簧刚度（在dynamicMeshDict中设置stiffness）。

我采用timeActivatedFileUpdate实现这一目的，即在不同时间段使用不同dynamicMeshDict文件，比如
time = 0 - 200，使用 dynamicMeshDict_1 (stiffness_1)
time = 200- 400，使用 dynamicMeshDict_2 (stiffness_2)

输出结果表明，程序确实按照预先的设置，在不同时间段读取了不同的dynamicMeshDict文件，但是结果似乎不太对。
以圆柱振幅为例，正确的结果应该是：

使用stiffness_1时，圆柱在0-200期间已经产生周期振动，且最大振幅稳定在 a1； 当在 time=200 改变刚度为stiffness_2时，圆柱振幅有一个突变，并逐渐趋于周期振动，最大振幅稳定在a2.

但目前使用timeActivatedFileUpdate的结果是，圆柱振幅逐渐增大并最终稳定在a2，0-200的过程没有表现出来。

请高手指教！！

@东岳 上次问完问题，就没用并行计算，这两天要用了，哈哈哈

我要刷火箭

矩阵求解器发散了，原因很多，不好判断

@程迪 scalar t 没有用到？

可以尝试自己写代码试试啊

这种由发布者发布的视频资料，类似文章的通讯作者，应该具备解答问题的责任，所以你最好联系视频发布者，在这里问，我不好评价别人的工作。
小提示：你的看起来比较粘。

@东岳 瞬态问题的话，是不是就是内残差作为收敛标准，初始残差不看的？

问题背景：将sixDoFRigidBodyMotion中linearSpring的参数stiffness从常数改成随时间变化的函数 stiffness = f(t)
f(t) 可以是分段函数也可以是任意函数，t为目前的计算时间。

因此我对$FOAM_SRC/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion/restraints/linearSpring目录中的.C和.H文件进行改动。但是在获得当前计算时间这一步遇到了问题。。尝试了多种方法都无法编译，以下是我尝试的一些语句和报错：

(1) scalar current_time = db().time().value()
Error: ‘db’ was not declared in this scope
(2) scalar current_time = runTime().value()
scalar current_time = runTime_.time().timeOutputValue();
Error: ‘runTime’ was not declared in this scope
(3) scalar current_time = timeDirs.last().value()
Error: ‘timeDirs’ was not declared in this scope
(4) scalar current_time = U.runTime().value();
scalar current_time = U_.thisDb().time().value();
Error: ‘U’ was not declared in this scope
(5) scalar current_time = mesh_.time().value();
Error: ‘mesh_’ was not declared in this scope
(6) scalar current_time = this->db().time().timeOutputValue();
Error: ‘const class Foam::sixDoFRigidBodyMotionRestraints::linearSpring’ has no member named ‘db’

挺简单一个问题，在网上找了很多方法，都失败了。。。不知道是不是因为没有包含需要的头文件？或者除了修改库文件，是否还有其他方法修改常数参数为time-dependent参数？请大神指教！！

库朗数如果放在80年代可以引一引，不过现在库朗数基本成为了一个必备内容，不引用应该也很安全

:144: Yes 我同意 @huoyanxiadebing 说的。每个湍流模型有自己适用的流型。这方面需要特定的研究领域的专家才能有很深的理解。

@kimy 在 Code production of k ? 中说：

@东岳 Many thanks. In fact, I understood this procedure. One confusion is that the formula implemented in openfoam is not consistent with that mentioned in the Spalding' paper.WeChat Screenshot_20200504110942.png

Nobody is interested in this issue?

@wwzhao 好东西！谢谢

小白刚接触openfoam，因为做振荡浮子有关的课题，所以在学习floatingObject算例。
有几个问题想向诸位前辈请教：
1、constant文件夹下存在两个动网格文件，是两个文件同时起作用吗？还是只有dynamicMeshDict文件控制动网格？
77152af7-c361-4b22-a3bd-a9778793fe8b-image.png
2、dynamicMeshDict中joint部分的参数含义是什么？我感觉与六自由度的constraint有关，不知道对不对？
f4c7d673-3d16-46df-b663-e0acf8965d52-image.png
3、如何通过pointDisplacement获取浮体的运动位移？是否与forces一样可以在controldict文件中设置functionObject来获取？

还有一个myFoam，github可以搜到，里面有个lusgsFoam求解器

我理解的欧拉方法就是网格不变形，材料在网格间流动。但是stressFoam的基本变量是位移。那么stressFoam还是欧拉的范畴吗，也就是说计算过程中网格会有变形(位移)吗

后续有什么进展么？

@东岳 谢谢老师，我调了设置之后kepsilon模型的计算结果基本一致了。因为在对比不同湍流模型对喷流计算结果的影响，发现还是有很多问题呀。。。愁

非常感谢老师们的指点，我去试试看

987.png

横向加密5倍试试

替代文字

@东岳 非常感谢！那我试着推推看这里:jingya:

@麦迪文 你好，其实openfoam-5.x就可以在状态方程中直接使用boussinesq，亲测可以，因为我最近也在用。我有一个疑惑，请问chtMultiRegionFoam共轭传热边界的换热原理是什么？（关于这个问题，我写了个求租贴，在这里），该原理在求解器中好像没有说明，也不知道这个原理是放到哪个头文件中去了，若你知道的话，请不吝赐教，谢谢！
另外，祝你五一快乐！

@BlookCFD 谢谢！已经解决:ok2:

不写的话默认是yes

网格分辨率与时间步长调小试试