感谢分享！！！收录在此 https://www.cfd-china.com/topic/3499/openfoam小代码

@xpqiu 感谢，我去查一下!不过目前我用break也能顺利终止程序，并行和串行都可以

@东岳 非常感谢东岳老师！按照您的指导，我设置了新的边界条件immersedboundaryEnthalpyFvPatchField，其他部分的设置目前都正常，但是想要使用浸没温度边界上的函数来更新焓值时会报错，我的程序部分如下所示：

const label patchi = patch().index(); wordList type = thermo.T().boundaryField().types(); if(isA<immersedBoundaryFvPatchScalarField>(type[patchi])) { immersedBoundaryFvPatchScalarField& TIb = thermo.T().boundaryField()[patchi]; scalarField & Tref = TIb.refValue(); }

之后程序就会报错，

error: invalid initialization of reference of type ‘Foam::immersedBoundaryFvPatchScalarField& {aka Foam::immersedBoundaryFvPatchField<double>&}’ from expression of type ‘const Foam::fvPatchField<double>’ thermo.T().boundaryField()[patchi];

应该是无法将基类转换为派生类，想和老师请教一下遇到这个问题要怎么解决呢，麻烦老师了！！！:xiexie:

问题背景是在simpleFoam计算中添加被动标量运输，因此在controlDict里添加了function的scalarTransport。但是在设置参数的时候发现resetOnStartUp。
在源代码中的discription是Reset scalar to zero on start-up。
现在比较疑惑的是，如果打开resetOnStartUp:
１．是在每个时间增量计算开始时把标量reset为０
２．还是一整个计算结束后，再次重新计算时把标量设为０
大家有经验吗

有现成的，可以看看6DOF模块的代码，里面又控制方程，也有写好的阻力

@东岳 这位是我后来一步一步检查，用checkMesh得到得具体的问题
30bf2ade-fc6e-412f-9486-f0083a60212b-image.png
如果感觉没有那么复杂的不规则边界，那么基本就是几何拓扑关系出了问题吧

pEqn里面用到了rAUc

@东岳 感谢!！!！!！

@东岳 踩坑踩多了而已:135:

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

@史浩 嗯嗯，好的，很详细了，谢谢您的回复！我看看mpi里的参数设置，琢磨下看看

@东岳 6f4bce5b-d530-4de8-8d67-dc76fc841d83-image.png 这个是出现问题前的最后一步，网格变形太厉害了，我在考虑是不是可以用圆尾的好一些

@东岳 感谢李老师回复！圆柱周围的网格细化我明白，但是在整个水槽的自由液面处网格细化也是用snappyHexMesh做的吗，而且这里的自由液面区域和圆柱周围区域的衔接也很好，看上去像是一个整体加密区域。但是我总感觉自由液面处的加密像是通过blockMesh实现的，我在这里还是有点困惑，李老师您能详细的指点我一下吗。

@东岳 而且我也不要把其他场映射过去，只需要颗粒速度场

@东岳 @Hunahum ，我在cfd-online上看到这样一段代码，想尝试着写出来，不过总是报错，请问你知道问题在哪儿吗1.png 2.png 3.png

@东岳 @xiezhuoyu 谢谢！

@东岳 谢谢！

@Samuel-Tu !就像这样，我显示的是位移，但是总有一层初始的膜（初始时刻的位移图）一样把我的结果给包住了，看不见我想要的结果。。(6ecae2f1-fd37-4b5e-9dc2-4041e9e8ceae-image.png 图片地址)

此处的压力无反射边界条件是哪一个？

各路大神好，想请教个问题。使用dynamicMotionSolverFvMesh网格运动求解器和displacementLaplacian动网格求解器时，设定什么样的时均运动边界条件可以指定为旋转运动，让旋转的角速度随时间均匀变化并且可以中途改变旋转方向。比如：从0度开始顺时针匀速旋转到10度，随后又逆时针匀速旋转到0度，这样往复旋转。