nonManifoldPoints与多域网格问题

[李东岳]

这个算例我打算做成多计算域求解器的案例放到OKS编程课来讲了。很有意思。外圈算流体，内圈算固体导热 :-)

[李东岳]

嘿！还用不了。这个网格是一套的。就是传统的单计算域网格。我自己做一个吧

[是嘿嘿牛啊]

@李东岳
李老师，我这个网格有大大的问题啊。结构物上方的初始液面呈现锯齿状, 然后我单独计算Openfoam可以计算，但是流固耦合的时候会发散。这应该是网格导致的问题把, 李老师可以给我一些建议来解决这个问题吗？
我的算例放在了github(https://github.com/NUGUY2022/OpenFOAM_Calculix_membrane.git)

[李东岳]

我没有你那个求解器撒

你固体算什么？变形？导热？还是啥。

[是嘿嘿牛啊]

@李东岳
算变形
我想知道我的流体求解器设置是不是有问题。

[李东岳]

那就用你的单独的流体计算域算流体试试就可以，用interfoam

[是嘿嘿牛啊]

@李东岳
流体域只有一个的时候是可以很好的计算的，我之前一直计算的单一流体域。

[李东岳]

把interfoam的算例发我看下

[是嘿嘿牛啊]

@李东岳
之前计算的单一流体域的是计算的弹性板

您的意思是让我只计算流体域的上半部分看看吗？把圆弧结构内部填实吗？

[李东岳]

就是你在一楼提供的网格，设置好interfoam求解器的算例，跑一下看看。你有尝试过么？

流体域只有一个的时候是可以很好的计算的

你这么回我的我以为你尝试过

[是嘿嘿牛啊]

@李东岳
我用interfoam计算，算例会很快发散。查看结果发现结构物附近的p_rgh值有异常的变化, 而且结构物上的液面有流速，位置正是液面锯齿处。但是我用wave2foam来计算的话，计算不会发散，可以计算完20s.我把两个算例放在了github(https://github.com/NUGUY2022/OpenFOAM_revise.git). 请李老师给点意见

[是嘿嘿牛啊]

@李东岳
李老师，如果我采用分区域耦合的话，整个模拟分成Fluid-outer, Solid, Fluid-inner三部分。那我的Fluid-inner的边界条件该如何给定？Fluid-inner有两个边界弧(interface)和底部边界。弧(interface)是需要用来耦合交换数据的。

[李东岳]

边界地方就是处理流体固体耦合的，涉及到信息的传递。这个我不知道你怎么处理的。如果用extend版本的求解器，他应该已经处理好了。我没用过extend

[是嘿嘿牛啊]

@李东岳
内部我就想模拟为一个装满水的容器，弧面我视作固体壁面，然后通过映射把力映射到结构求解器中，在结构求解器中给予材料属性来计算结构变形，然后把变形映射回OpenFOAM。底部边界我也看作壁面，就如同水槽中底部边界那样设置。我现在把边界条件设置成这样alpha.water p_rgh pointDisplacement U, 这样的边界条件无法计算，报错是这样的log.interFoam 。

[李东岳]

[14] Unable to set reference cell for field p
Please supply either pRefCell or pRefPoint

这个是非常普适性的CFD问题。跟你的算例没啥关系。你要在算例fvSolution里面设置pRefCell 0;pRefValue 0;

http://www.dyfluid.com/theory.pdf : 5.10.7 setReference() 函数

由此来看，我觉得你在普适性CFD上应该先走一遍，确保各类求解器算例都了解什么意思，然后再深入到你的课题。这是正常的学习思路。急于求成反而会导致困惑。

[是嘿嘿牛啊]

@李东岳
是的，我确实还需要多加学习！谢谢李老师