这个叫函数的重载，在面向对象的编程语言中是个普遍应用的特性。:high:

@队长别开枪 大佬，教程有了吗:140:

是欧拉-欧拉，确实如果引入边界条件的话不知道压力场怎么处理。谢谢分析，我再想想。加多孔介质主要就是物理上好解释，我现在觉得可以在wall的区域再附加一个指向wall外面的源项（根据具体问题可以指定方向），类似重力，以此抵消重力以及其它数值问题造成的穿透问题，相当于在这个源项下，wall区域的固体会被很缓慢地排出wall。我估计一开始如果就加上的话，可能就可以阻止固相进入wall，而且有多孔介质巨大的阻力系数在，排出的固相应该速度非常小，对wall外面的计算应该影响有限。

@dyb 你可以参考自带的案例——airfoil2D，里面有升阻力系数监测的方法。

这个估计很难了。你这个涉及流固耦合的问题。CFD这边做不鸟。但是你可以这么做，流场获得你需要的压强时程曲线，然后你在固体场中把压强加上去。这样是可以的。

@Samuel-Tu 非常感谢您的回复。讲的非常明白，我之前没有想到这些。其实，我想实现的就是一个热源点慢慢加热最终提高出口温度，或许我应该试试改成热源强度再试一试看。非常感谢您。

感谢分享！！！收录在此 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实现的，我在这里还是有点困惑，李老师您能详细的指点我一下吗。

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