关于重叠网格的疑问

[秋山结月]

本人先是跑了一下tutorials的overInterDyMFoam的算例，一切正常，然后将模型换成了自己的一个算例，然后问题来了。
openfoam识别不了两套网格的交叉区域，两边好像都是在独自在跑，没有插值，具体问题如下：

Creating cellMask field to block out hole cells

inverseDistance : detected 2 mesh regions
    zone:0 nCells:1331200  voxels:(35 35 35) bb:(-10.0000079951 -5.00000799512 -6.00000799512) (-4.87499200488 -1.62499200488 -0.874992004884)
    zone:1 nCells:238837  voxels:(35 35 35) bb:(-62.0000295973 -25.0000295973 -26.0000295973) (-35.9999704027 -14.9999704027 -15.9999704027)
Overset analysis : nCells : 1570037
    calculated   : 238837
    interpolated : 0 (interpolated from local:0  mixed local/remote:0  remote:0)
    hole         : 1331200

Creating interpolatedCells field 

从计算开始到结束都没有变化，interpolated的数量始终为0

主要有几个问题:

1. 该怎么改.....
2. 设置重叠网格跟setFields里设置的zoneID有关吗，如果无关，控制重叠网格的字典应该是哪一个？

[秋山结月]

@wwzhao 请教一下老师，非常感谢

[wwzhao]

看起来两个 zone 没有 overlap？

[秋山结月]

@wwzhao 对 是这样 两套网格似乎是在独立算没有交互
我后来又换了个算例，结果又有overlap了，然后今天又做了个，又变成零了。。
这个规律到底是闹哪样。。。

[秋山结月]

@wwzhao 想请教您一个问题

现在我在计算中遇到了remote单元不断上升的情况，然后接着就是某一步速度或者压力突然迭代次数剧增，然后计算发散。

我想请问一下您，一般背景跟overlap的网格差距设定为多大合适?

另外就是PIMPLE各个设定，看您在别的帖子提到需要提高外循环的次数，一般多少比较合适呢?

[wwzhao]

@秋山结月

保证插值区域处不同网格块的单元尺寸相当就可以。

提高PIMPLE外循环可以适当增加时间步长。你的发散可能不是这个问题，网格的可能性大些，不过你可以尝试增加PIMPLE外循环次数，看看能不能解决问题。

[秋山结月]

@wwzhao 在 关于重叠网格的疑问 中说：

处不同网格块的单元尺寸相当就可以。
提高PIMPLE外循环可以适当增加时间步长。你

我目前的办法是为了保证插值区相差不大，在背景网格里用snappyHexMesh选了个refineBox来局部加密
pimple的话我是设定一个高次数的外循环，再用一个relTol控制次数（差不多4次左右就会满足relTol）
但是还是会发散，可能真的是本身网格有些问题
总而言之先谢谢老师了 我再试试调整一下网格

[秋山结月]

@wwzhao 另外想请教您一个问题，一般重叠网格（贴体的那个网格），尺寸应该取多大比较合适呢。
我看tutorials里的构件都比较小，可能实体是 0.2×0.3×0.4的大小，重叠网格大概就是0.4×0.6×0.8

但是我要计算的物体几何尺寸可能有几m到十几米，这样也重叠网格的尺寸也要翻倍？
还是说，只要能够保证重叠区包裹住实体就可以了？

[wwzhao]

@秋山结月 重叠网格的block size没有特殊要求，包住物体就行。

[秋山结月]

@wwzhao 好的好的 谢谢老师~~

[秋山结月]

@wwzhao 老师我想请教一个问题。我看到一些论文把重叠网格分为插值边界单元，贡献单元以及孤点单元，我想问着是否跟OpenFOAM的重叠单元的interpolated from local , mixed local/remote 以及 remoted相对应呢？

我计算的时候遇到了remoted单元不断上升的情况（初始0，计算100s后可能上升到几千），但是我的计算确实也没有发散，能够进行下去，如果remoted对应孤点单元，是不是代表我这个计算设定是不准确的呢？

[wwzhao]

@秋山结月 http://cfd-china.com/post/11197

[秋山结月]

@wwzhao 看起来显示的只是个并行计算的信息？
总而言之谢谢老师的指点了

[Yingdong]

@秋山结月 背景网格使用sanppy加密部分区域好像真的会发散，而且p_rgh的求解矩阵也会变为对称，不加密是非对称的，请问您最终怎么解决的这个问题

[z597288]

@秋山结月 您好，请问您解决这个问题了吗，我目前也遇到了 应该怎么改呢