关于重叠网格的几点疑问,想请教一下相关大神
-
我先是用flunet meshing分别绘制了旋转区网格与背景网格,然后在fluent里将两套网格一起输出为mesh文件,再用fluent3DMeshToFoam转到openfoam中,再用topoSet文件选取所有旋转域网格,然后通过反选选取流场域网格,然后在Setfields里将旋转域定为zoneID 1(差值区),并将流场域定为zoneID 0 (计算区),并在zoneID字典里定义旋转域的overset边界为1。想请教一下处理重叠网格这样的做法是否正确
同时我在paraview里面看zoneID,发现整个旋转区的zoneID为1,流场为0,
然而我在计算的时候发现,其实旋转区里面的也是作为calculated,而差值区interpolated只有4万网格,这也挺符合常理的,但是却与paraview里的显示不大一致所以不知道自己这样处理重叠网格是否正确,以前也没有相关经验,就单单学了overPimpleDyMFoam的几个例子就拿来实践了,但里面的例子也是过于简单,所以想请问一下相关重叠网格专家,我的方法正确吗?
-
没用过openfoam的重叠网格 不知道 @wwzhao 懂不懂这块
-
zoneID
只是用来标记background网格和component网格的,所以图示中paraview的显示结果是正确的,之后的网格挖洞的操作由overset库自己完成,位于createCellMask.H
,setCellMask.H
,createInterpolatedCells.H
和setInterpolatedCells.H
这四个文件里,所以要在paraview中查看单元属性应该显示cellTypes
场,cellTypes
场不存在于初始文件夹0
,位于其余时间文件夹内。
-
@队长别开枪 感谢大神解答,还有一点疑问,关于计算过程中的inerpolated项后面的local与remote,这两项是何意。运行了算例之后,发现差值区都是local那一项的。可是当具体实践的时候,自己导入的网格计算时,显示的都是remote。是否是因为组件网格与背景网格之间的差异导致的。
-
@队长别开枪
这是我的fvsolution的设置
这是fvschemes的设置
对于这两个文件,想请教一下有没有问题很大的地方,因为最终得到的结果跟实际差异有点大,不知道是不是网格的原因还是这块求解设置的原因
-
@gentle OpenFOAM有嵌套网格求解器吗?我也想学习嵌套网格。
-
@gentle local和remote是
inverseDistance
插值方法输出的信息,具体的我也不是很清楚,是它内部的一个插值过程,源代码在/src/overset/cellCellStencil/inverseDistance/inverseDistanceCellCellStencil.C
,如果只是使用重叠网格的应用,这个我认为可以暂时不管它。如果不对,应该是网格的问题,我是用STAR CCM+
或者Pointwise + cfmesh
准备重叠网格的,没有问题,也没关心过这个local和remote,你可以参考一下。
-
@gentle 你把
constant/polyMesh
这个文件夹的文件树,constant/dynaminMeshDicttr
文件,constant/polyMesh/cellZones
文件贴出来看看。
-
@gentle
constant/dynaminMeshDict
文件
-
local: stencil 都在本进程(并行)
remote: stencil 都在其他进程(并行)
mixed: 以上二者的混合
-
-
@wwzhao 明白了,感谢大神回复,想请问下我前面的fvsolution与fvschemes设置的是否合理,因为所计算的风机阻力系数与实际情况差的有点多,不知道问题是否出在这里
-
@gentle 大体上应该没问题,建议增加
fvSolution
里面PIMPLE
内外修正次数,尤其是对于3D算例。dynaminMeshDict
里面转轴axis
设置有点奇怪啊,一般都是平行于某一个主坐标轴,还有omega
注意一下单位(rad/s)。
-
@gentle 阻力系数不对的原因太多了,不一定是重叠网格的问题。你试过滑移网格吗?
-
@队长别开枪 因为我的转子有一个5度的俯仰角,所以我这个是转子的法线方向,不知道我这么理解是不是正确
-
@wwzhao 还没试过滑移网格,因为刚上手就是重叠网格,感觉可以先学一下滑移网格哦
-
@队长别开枪 如何对重叠网格计算的结果进行后处理?导入tecplot有问题,洞网格还在。如何在parafoam中进行重叠网格的后处理?可以不显示洞网格吗?
-
@我是河滩 两个软件都需要额外操作才能正确显示重叠网格。tecplot可以参考链接
http://2016.oversetgridsymposium.org/assets/presentations/tutorials/ogs2016-fowler-tecplot.pptx
和https://fun3d.larc.nasa.gov/session9_july.pdf
,paraview可以参考https://twitter.com/openfoam/status/883256187390750720
,里面说了The cellTypes field can be used to subset only the active set of cells - in ParaFoam/ParaView: Filter->Threshold
,最小值0,最大值1。
-
@队长别开枪 太感谢了,我下来试试。
-
@队长别开枪 你好。overPotentialFoam中有一个cylinder算例,但计算结果并没有cellTypes变量,因而无法区别网格的类型,在哪可以定义cellTypes并让其输出啊?
我试了overSimpleFoam和overPimpleFoam求解器均有cellTypes,但不知道这个是在那定义以及如何定义的。希望能告知一下,谢谢!
-
@我是河滩 There is no Chinese input tool in my Linux workstation, so I replied this in Eng. Take
plus 1812
version for instance. It is simple if you prefer to writecellTypes
field in the solver ofoverPotentialFoam
, just addmesh.write()
to the end ofoverPotentialFoam.C
, i.e.,p.write(); } mesh.write(); // Add this runTime.functionObjects().end();
mesh.write()
will call the member functionFoam::dynamicOversetFvMesh::writeObject()
wherecellTypes
andzoneID
fields have been defined. They will be written to0
time folder inoverPotentialFoam
.
-
@队长别开枪 谢谢你耐心的回复,明白了。
-
@队长别开枪 感谢分享,楼主好人
-
@队长别开枪 你好,在利用overLaplacianDyMFoam 求解器计算heatTransfer算例时,我找了一下/system /constant /0下的文件中,并没有指定物面hole的平移速度,但计算过程中,屋面hole是运动的,我感觉很奇怪,屋面hole的运动定义到底是在什么地方啊?
-
@我是河滩 在
pointDisplacement
里面,"(free|hole)" { patchType overset; type uniformFixedValue; uniformValue table ( (0.0 (0 0 0)) (1.0 (0.2 0 0)) (2.0 (0 0 0)) ); }
x
方向速度分量在0~1.0s
时间段内均匀由0
加速到0.2
,1.0s~2.0s
时间段内均匀减为0
。
-
哦,明白了,谢谢。