(6.4.1)--5.4-计算平面的强守恒型方程推导讲义.pdfVIP

《计算流体力学基础》讲义 第五章 网格与相应变换 5.4 计算平面的强守恒型方程推导 为方便说明问题起见，我们以二维非定常流动为例，其强守恒型控制方程可以写成： U F G + + 0 (5.37) t x y 那么将以上控制方程变换到( , )计算平面后，方程是否还是强守恒型形式？即，它 在变换之后能否写成： U F G 1 + 1 + 1 0 (5.38) t   式中，F1，G1 是原始通量F 和G 的适当组合。 如果是可行的话，我们就能够在计算平面上利用2.10 节中描述的CFD计算中强守恒形 式方程的所有优点。答案是肯定的，下面我们来具体看看是如何实现的。 首先，根据式(5.2)和式(5.3)给出的偏导数变换公式将式(5.37) 的空间变量变换到计算平 面，得到： U F  F  G  G  + ( ) + ( ) + ( ) + ( ) 0 (5.39) t  x  x  y  y 然后，将上式两端同乘以雅克比行列式J，得到： U F  F  G  G  J +J ( ) +J ( ) +J ( ) +J ( ) 0 (5.40) t  x  x  y  y  / x  先不处理式(5.40)，我们先考虑将JF( )项对 求偏导数，其展开式为： [JF ( x )]  F   J ( ) =+F (J ) (5.41)  x   x 将式(5.41)进行重新整理，可以得到： F  [JF ( x )]   J ( ) =−F (J ) (5.42)  x   x 式左端刚好是式(5.40)的左端第二项，我们稍后再代入。同样道理，将 JF(  / x )对 求  偏导数并进行整理，得到： F  [JF ( x )]   J ( ) =−F (J ) (5.43)  x   x