Cyclic周期性边界条件

[Calf.Z-DNS]

@cxzhmg 我之前跟同学讨论的是泄压可能是因为设置cyclic的效果就像是一道虚设的面，面上可以流进流出，而真正有阻挡作用的是wall壁面，因此channel就像是上下两个平板之前存在一个自由流动的空间一样，于是留着留着就泄压。但如果像圆管那样四周都是wall，就不会泄压。另外我采用了低马赫数假设，压力的变动不影响物性计算。因此无论压力是否泄压，我的物性计算还是正常的。

[cxzhmg]

@Calf-Z-DNS 圆管不会泄压吗？我计算二维的box，左右壁面，也是会出现泄压的问题。如果cyclic确实是周期边界的话，我觉得即使是3维无限大平板间流动，或者三维全周期流动，也都不应该出现泄压或者质量不守恒的问题，因为从任何边界流出去的都会从对应的neighbour流进来，所以我怀疑openFOAM里面的cyclic并不是准确的periodic，当然具体还是要看代码，我最近正在看。

[XieXiaoyang]

@气固两相流 您好，我在simpleFoam里加了能量方程想测试一下fixedJump边界，是三维槽道算例，边界上会出现如下图的情况，在槽道出入口处有些网格上的结果感觉比较异常。

该算例中采用fvoptions添加动量源项实现动量方程的周期边界，因此速度的边界就是cyclic，在出入口处云图就很平滑：

除此以外，我在对一个横掠管束的算例中对压力采用fixedJump边界的时候，也出现了这种现象(压力场)：

但是在采用fvoptions添加动量源项时边界处结果就比较正常(虽然对管束这么实现周期边界是不正确的)：

还有就是采用这种边界，温度场好像会因为数值误差的累计(存疑)使得整体计算区域温度越算越高，下图分别是采用simpleFoam刚开始计算和很久之后的温度场分布，明显发现计算很久之后的温度场温度整体偏高，按理说稳态计算应该是不变的才对。(另：这个网格原来是用来算DNS的，所以我觉得算很久之后出现的这种温度场波动应该是算出湍流脉动了)：

这是我最近调试fixedJump边界遇到的一些问题，烦请大佬有时间帮忙瞅瞅~
如有描述不清的地方，还请多多指教！

[李东岳]

@xiexiaoyang simpleFoam怎么加入的温度场

[XieXiaoyang]

@李东岳 就是在simpleFoam求解器中添加了一个温度的标量输运方程

[李东岳]

@xiexiaoyang https://openfoam.org/release/2-2-0/boundary-conditions/ 我看了下这个边界条件的介绍，看起来更倾向于在计算域内部添加baffle然后实现某种jump。类似计算域内部添加的fan边界条件。所以你目前这种前后面肯定不是baffle，经过你测试，如果用起来有问题，那应该就是不能用这个fixedJump边界条件。

simpleFoam加温度，为何不用rhoSimpleFoam。

[XieXiaoyang]

@李东岳 但是就计算结果来说，确实实现了温度边界的阶跃，而且楼上有老师对二维算例的压力边界计算结果，效果也都还好。
就是不太明白为什么会出现这样部分网格计算结果不好

[XieXiaoyang]

@李东岳 没用rhosimpleFoam是因为它解出来的好像是内能e，还涉及到热物性相关内容，不敢乱来hhh

[XieXiaoyang]

@李东岳 东岳老师，jump边界上的问题解决了，可能是因为在并行计算的时候，没有将关联的周期边界放在一起引起的，加了：

这几句以后的计算结果种jump关联边界上的问题基本就没有了。

[李东岳]

@xiexiaoyang 感谢反馈！

[李东岳]

我测试了一下这个fixedJump，会一直一直的导致速度上升，你多算几个周期看看，

[XieXiaoyang]

@李东岳 是对压力边界采用fixedJump边界吗？
我先用层流稳态算例尝试一下看看~

[李东岳]

@xiexiaoyang 是的

[XieXiaoyang]

@李东岳 李老师，对于周期边界下的压力场，piso循环中有pEqn.setReference(pRefCell,pRefValue);用来确定压力场。那在计算周期边界下的被动标量时，应该用什么来对温度场进行修正呢？