interFoam水槽试验验证,明渠流速断面现象随网格增加异常变化问题
-
@李东岳 老师这是之前在不同湍流模型下的结果:
值得注意的是,这些网格数是在97万网格数的结果,并且网格与最开始展现的模型存在一些差异,由以上不同湍流模型的流速断面图可知下游两个断面测点流速完全是跟物模相反的现象(物模顶部流速大,底部流速小,而这些都是底部流速大顶部流速小),后来经过调试发现当把水流进入逆坡段部分添加10层边界层网格(之前只加了1层)后,现象会有所好转,在低网格数下会出现与物模比较相似的现象(最开始的那几个模型结果),但是当网格数量增加后现象却与物模相差反而越来越大。
并且增加了边界层后,我对比过RNGke与komegaSST的流速结果,基本没什么差别。我也尝试更改添加边界层的设置(边界层尺寸是绝对尺寸或者相对尺寸),但是结果都没有得到很好的改善。 -
@李东岳 李东岳老师,当我采用流速出口时,当出口的流量大于进口流量时,出口的过流面积会变大,出口的流量也就变大了,所以能够保证当水位稳定时进出口的流量也稳定了,并且在后处理过程中我统计了不同断面的质量流量通量(流速出口和压力出口),发现是守恒的(各断面的流量误差在基本在2%以内)。但是下游的两个监测断面流速分布与物模对不上,物模试验测出来的是沿底面向重力反方向流速逐渐增大(每个断面中轴线共五个测点,第一个断面:从0.97m/s增长到1.61m/s,第二个断面:从1.02m/s增长到1.43m/s),但是数模计算出的结果最开始的时候时底部流速大顶部流速小,下图是之前数模逆坡后第一个监测断面的流速分布
,后来我在逆坡段添加了边界层后情况有所好转,相同位置断面流速分布
,但是这个模型的主要问题目前在于两点,1、虽然整体流速分布现象正确但是数值上误差还较大(物模是0.97m/s-1.61m/s,数模结果1.21m/s-1.43m/s),误差在20%左右。2、当前现象能够匹配物模的模型在当前网格数下误差较小,但是当我仅加密blockMesh背景网格时则断面流速分布则会变化(变为底部流速大,顶部流速小,1.48m/s至1.24m/s),各个测点的误差最大增加到52%)
还有一点令我困惑的现象,我目前所提到的模型背景网格都是严格沿着xyz轴方向的,背景网格示意图如下:这样构建的模型基本能够保证逆坡段的底部流速大,但是明渠水平段的流速分布不对
然后我尝试变更了背景网格的画法,让逆坡段的网格是倾斜的,示意图如下,这样构建的模型结果会有完全不一样的结果,速度云图如下
-
@李东岳 感谢您的回复,正如您所讲,我这个是两相流,我为了判断进出口的流量是否守恒,我在controlDict的function里面添加了
“inletFlux { type surfaceFieldValue; libs (fieldFunctionObjects); writeControl runTime; writeInterval 1; log false; // Output field values as well writeFields false; regionType patch; name inlet; operation sum; fields ( rhoPhi ); } outletFlux { $inletFlux; name outlet; }”
内容以监测进出口断面的质量流量(主要排除空气的干扰,空气部分的质量流量与水的质量流量相比可以忽略不记),结果如下
,所以我判断整个水槽试验在进行到最后,系统流量是守恒,但是就是下游的流速分布现象不对。