(6.3.1)--5.3-几何变换系数与雅可比行列式讲义.pdfVIP

《计算流体力学基础》讲义 第五章 网格与相应变换 5.3 几何变换系数与雅克比行列式 首先回顾一下上节课学过的变换关系式(5.1a)至(5.1c)，其反映的是计算空间自变量 ,, 与物理空间自变量x，y，z 之间的变换关系。  (x, y ,t ) (5.1a)  (x, y ,t ) (5.1b)  (t) (5.1c) 观察上节课所学的式(5.2)至(5.15)，可以发现变换关系里涉及网格几何形状的项，如  / x ， / y ， / x ， / y 称为几何变换系数metrics。如果式(5.1)给出的变换 是解析的，则可以得到几何变换系数的解析值；如果式(5.1)给出的是数值的对应关系，那么 几何变换系数须采用有限差分法计算。但是，在很多实际应用中，这种变换是以方程(5.1a) 至(5.1c)的反函数形式给出的，那么这种情况下如何求解几何变换系数呢？ x x (,, ) (5.18a) y y (,, ) (5.18b) t t() (5.18c) 我们观察式(5.18)给出的用反函数所表示的变换关系： 式中自变量是,, ，但是式(5.2)至(5.15)中的几何变换关系如 / x，/y 等，是 以x、y、t 为自变量的偏导数。因此，为了由反变换式(5.18)计算这些方程中的几何变换系 数，需要将 / x，/y 等与其逆形式x / , y /  等联系起来,下面将详细讨论。 考虑流动控制方程中的一个因变量，例如 x 方向的速度分量 u,令 u=u(x,y)，其中 , ,  x=x( )，y=y( )。u 的全微分形式可以写成式(5.19)。进一步可以得到u 偏 的表达式 (5.20)，u 偏 的表达式(5.21)。这两个方程可以看作是关于u / x, u / y 两个未知数的方  程组。 u u du dx =+ dy (5.19) x y u u x u y =+ (5.20)  x  y  u u x u y =+ (5.21)  x  y  用克莱姆法则可以很容易求解该方程组，可得u / x 的解：