《表2 梯度求解方法的数值精度验证(Grid G)》
体心梯度的求解精度决定了求解器的整体精度。梯度求解为一阶精度时,求解器整体有二阶精度。梯度求解为零阶精度时,求解器整体只有一阶精度。图9给出了不同梯度求解方法预测的阻力系数。坐标都是在对数坐标下的表示。横坐标为网格的特征尺度,取为,其中N为总的单元数目。格林高斯法随网格加密有明显的非线性变化,说明此时还未完全进入网格收敛解区域,根据较密的两套网格估算的数值精度也在一阶左右。其他三种方法阻力系数随着网格加密线性变化,且数值精度在二阶左右。在一般的网格下,使用距离倒数权插值到节点的节点格林高斯法求解梯度时只有0阶梯度,但是如果节点周围的近似满足:,则该种方法计算得到的梯度也有一阶精度,此时阻力系数的预测是二阶精度。表2给出了在网格G上通过阻力系数计算得到的空间离散精度,进一步证实了上述的观察。
图表编号 | XD00200623100 严禁用于非法目的 |
---|---|
绘制时间 | 2019.10.01 |
作者 | 赵辉、张耀冰、陈江涛、邓有奇 |
绘制单位 | 中国空气动力研究与发展中心、中国空气动力研究与发展中心、中国空气动力研究与发展中心、中国空气动力研究与发展中心 |
更多格式 | 高清、无水印(增值服务) |