@多相流FOAM 热计算直接给定热流密度或者用对流换热算比较合适

@wangxiaowei hello，我可能理解错了你的问题，我先用imageJ软件从ct扫描切片中提取出土壤介质，生成了stl文件，然后用snappyHexMesh画出孔隙的网格

哦好像是定义没有弄明白，miu和nu是不一样的····

@yhdthu 好的，谢谢老师

@东岳

我没有找到李老师说的“tutorial: decayIsoTurb”。不过在OF里找到了这个前处理程序7dbae584-0745-4a52-bc5f-ddfbc18a3618-image.png 正在看这篇文章，还看到了一篇中文的“一种改进的均匀各向同性湍流初始化方法”

我看的比较慢，有些基本概念还不太清楚。:136:

@D-Benjamin 具体一共都有哪些函数我也不是很清楚，具体函数和怎么用通常都是看看openfoam自带的边界或者求解器或者其它的应用有类似你想获得功能，就去翻源代码。gsum应该也可以对边界上的网格体积或者面进行求和，应该就能实现你说的功能。

edit -> find data
如果数据只在网格或者点，可以先用filter point data to cell 或者cell data to point筛一下。

感谢分享 :146: :146:

插值只对patch进行，和是否并行无关。并行对local patch进行插值。

@Izumi 这个没添加，我回去试试

@东岳 结果.tar
这个是发散的相分数，速度，压力图

个人理解，如果不改transport property的话，引用的nu应该要加上phase1的关键字，比如可能要输出的应该是nu.air。关键是看nu在你用的solver里注册成了什么名字，这个名字可能和输入的名字完全不一样（注：是注册的名字不是变量名）

@东岳 多谢东岳老师指点:tishizi: 我再改一改边界条件算一下，搞定了再发帖，用这套软件的道友太难找了:zoule:

@东岳 好的，谢谢东岳老师的回答，另外想问一下老师，DPM中颗粒碰撞力的PairMOdel有没有什么参考资料？？谢谢

感谢分享！alias确实更方便！

@xinxinzhou1995 可以交流一下吗，关于动态接触角模型

残差图IMG_20200601_201953.jpg

这个叫函数的重载，在面向对象的编程语言中是个普遍应用的特性。: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实现的，我在这里还是有点困惑，李老师您能详细的指点我一下吗。

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

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

@东岳 @xiezhuoyu 谢谢！

@东岳 谢谢！

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

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

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

我记得构造哪个class的时候就是默认了1

@羽之下 动网格用的自带算例wingMotion为模板的。目前问题找打了，是因为网格变形过大。

研究流体携带的砂粒在一定的工艺参数条件（比如注入速率、携砂液粘度、砂子粒径、砂子密度、裂缝宽度、裂缝角度等等）下的运移和分布规律，除了使用DEM外，有没有什么较为便宜的方法，可以得到相类似的结果呢，DEM跑起来太。。。。。

@hoversoar 好的，我看看，谢谢您

@xpc 这个大佬毕业去商汤了，年薪80万 :146:

:high: :high: :high: :high:

@东岳 感谢东岳老师指点:chitang: