有限差分法中的贴体坐标变换

摘要

cqvip:分析了有限差分法曲线贴体坐标系下守恒型控制方程的推导过程,认为在离散条件下所采用的数学恒等式不成立,推断目前CFD广泛采用的齐次方程是原始Descartes直角坐标系下方程的近似,提出增加源项的非齐次方程作为离散等价方程.采用数值实验研究了源项,结论是大部分情况下源项不等于0,且对数值解的影响大于差分格式的截断误差,在分析了引起源项非0的原因和推导过程后,提出源项离散的相容性准则.利用坐标变换系数和守恒型方程对流通量的特性,建立了源项隐式计算的耦合算法,通过数值实验证明耦合算法有效消除了坐标变换引起的误差.