为什么计算域尺度不同会导致计算发散呀？？

[东方白杨]

如题，最近我用SRF求解器算DNS，同样的进出口条件、网格大小，以channel395为例，原始的网格尺度是4m×2m×2m，主流速度0.1335m/s，设置旋转速度rpm=120，这样计算出来使用probe监测任意一点速度，发现相对速度很离谱，随着时间一直升高，最后都10的6、7次方相对速度了。
反而，我在blockMesh里把网格scale设置为0.01，也就是计算域变成了0.04m×0.02m×0.02m，其他的设置都一样，然后计算就能收敛稳定了，如下图是通道中心相对速度随时间变化：

这到底是为什么呀？？？其它条件都一样，就改了网格尺度就算不了了，求大神解答~

[东方白杨]

进出口还都设置的是普通的patch条件，不是周期性条件，已经排除了离心力做功这个因素
@李东岳 老师

[学流体的小明]

DNS计算要考虑到网格尺寸的，你改变计算域大小，网格尺寸也改变了，$ \Delta x^+$这些都改变了。
可能是因为大计算域时候网格尺寸不满足DNS的条件？就是计算域尺度方面的一些问题。

[Afr1yne]

@学流体的小明 老师 buoyantpimplefoam计算大尺度也会发散 您遇到过类似情况吗

[东方白杨]

调了半天还是不行。。。。网格不管加密多大，都算不了。
channel395 网格尺度为0.04m×0.02m×0.02m，旋转中心为x=-0.2m，转速120rpm，设置进口来流0.1335m/s。probe监测x=0.02m，y=z=0.01m处，计算可以稳定，稳定后速度为0.160271m/s
当保持上述设置，只把旋转中心为x=-2m，计算没法稳定，速度在1~2m/s波动，显然不正确。
按理说离心力是一种有势力的话，他只和起点终点有关系，我改了旋转中心应该不影响离心力加功的多少啊？？？上述速度都是相对速度Urel哈

[学流体的小明]

@Afr1yne 我还读在读博😂
我是拿pimpleFoam算的，网格都是精心设计过，和前人文献建议的一致，即使这样，OpenFOAM算出来的结果都不是很好。还是使用了更加精细的网格，得到的才和文献的DNS结果一致。

[学流体的小明]

旋转中心设置为-2 m，计算域还是0到0.04m，计算域距离旋转中心就很远了啊，旋转产生的速度$v=\omega \times r$也应该很大的吧。

[东方白杨]

@学流体的小明 你说的这个是绝对速度，看相对速度的话，和转轴没多大关系
我突然有个问题，目前算槽道流都是按6.4h-2h-3.2h比例算的，那我是不是可以计算域改成0.64m-0.2m-0.32m，和计算域为6.4m-2m-3.2m是一样的吧，就是只要比例一致就行了嘛？？

[学流体的小明]

h大小会影响到槽道的无量纲参数，比如摩擦雷诺数$Re_\tau=u_\tau h / \nu$，你要算湍流的话，网格要跟着$Re_\tau$去画。当然不能简单地缩放。

我不理解这句话

你说的这个是绝对速度，看相对速度的话，和转轴没多大关系

[东方白杨]

@学流体的小明 说错了，我的意思是v=w×r，这个算的是圆周速度，而我们看的是相对坐标系下的相对速度Urel，从NS方程来说它多受到了离心力和哥氏力，哥氏力和主流不在一个方向，而离心力是有势力，如果我进口主流速度都一样的话，不管旋转半径多大，与进口差距一定距离的点，其受到的离心力做功都是一样的，所以理论上改变旋转半径大小，那个点的相对速度应该基本是一样的吧。。。。