openfoam中压力梯度计算不均
-
最近发现在openfoam计算压力梯度时,用
blockMesh
生成的均匀网格,给定压力场为均匀初始场,当压力初始场为零时,全场压力梯度为零,但是当压力初始场为101325
时,梯度在x,y
(2D)方向上绝对值大小为1e-6
(由于网格质量引起的数值误差?),虽然这个值量级比较小,但是目前计算的系统网格长度量级为1e-5
,会对最终的结果产生误差,我分别在openfoam2.2.x
以及openfoam6.0
上测试都存在该问题,压力梯度计算方法为Gauss linear
,不知道有没有人碰到类似的情况,或者有没有推荐的绘制网格质量较高的软件,我试试在不同网格下计算看看
下图为压力初始值为101325
的压力梯度分布
(另附测试算例链接https://pan.baidu.com/s/1uFbx_xqbO-5QAbcmTSr-fA)