您好,欢迎访问三七文档
当前位置:首页 > 建筑/环境 > 工程监理 > 高等计算流体力学讲义(1)
1高等计算流体力学讲义(1)第一章计算流体力学基本原理第1节流体力学基本方程一、非定常可压缩Navier-Stokes方程不计质量力的情况下,在直角坐标系中,守恒型N-S方程可以写为下列向量形式:()()()0vvvtxyzUFFGGHH,(1)其中uvwEU2()uupuvuwEpuF2()vvuvpvwEpvG2()wuwvwwpEpwH,0xxxyvxzxxxyxzTuvwkxF0xyyyvyzxyyyyzTuvwkyG,0xzzyvzzxzzyzzTuvwkzH。如果忽略N-S方程中的粘性和热传导,得到的简化方程为Euler方程:20txyzUFGH。(2)方程(1)、(2)称为向量守恒型方程。其重要特点是:连续、动量和能量方程被写为统一形式。其中,,,,,,,vvvUFGHFGH均为列向量,U是方程的解向量,称为守恒变量;,,,,,vvvFGHFGH称为通量(flux),具体说,,FGH为无粘通量,,,vvvFGH为粘性通量。所谓守恒型方程,是空间导数项为散度的形式的方程。(1),(2)式所示的向量型守恒方程,实际上仍然是散度形式。显然,(1),(2)式的另一种等价形式为:0tUE,(3)其中()()()vvvEFFiGGjHHk,或EFiGjHk,通量张量,,,ijk为直角坐标系三个坐标轴方向的单位基矢量。把(3)式在任意固定的控制体上积分,并利用Gauss公式,有0SddStUEn。(4)这就是守恒积分型方程。可见,守恒的微分、积分型方程之间有直接的联系。(4)式是我们以后将要讲到的有限体积方法的出发方程,而(1)、(2)或(3)是则是有限差分方法的出发方程。二、流体力学方程的简化形式根据具体流动状态,N-S方程可以进行各种简化。简化的形式及其适用条件是理论流体力学的重要研究内容之一。这里我们对于各种简化方程作一归纳,见下图:3图1.N-S方程的简化形式三、曲线坐标系中的基本方程当求解域的形状比较复杂时,计算流体力学方法通常在曲线坐标系中实施。因此,有必要得到曲线坐标系中流体力学基本方程的形式。在曲线坐标中,矢量可以采用在直角坐标中的分量形式,也可以采用协变或逆变分量,基本方程也将因此呈现出不同的形式。最简单,应用也最普遍的形式是矢量分量为直角坐标系中的分量。下面,我们讨论这种情况下的流体力学基本方程。直角坐标到曲线坐标的变换及其逆变换关系为:(,)(,)xyxy(5)(,)(,)xxyy(6)1、导数的变换对于一阶偏导数,根据链式求导法则,有xxx。(7)同理可得4yyy。(8)对于二阶偏导数,有222222222222222()()()[][]()2()xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx。同理可得222222222()2()yyyyyyyyy,222222()xyxyxyxyyxxyxy。把导数的变换关系代入微分方程,就可以得到微分方程在计算平面中的形式。以直角坐标系中的Laplace方程22220xy为例,把上述二阶导数的变换关系代入上述Laplace方程,得222222222[()()]2[][()()]()()0xyxxyyxyxxyyxxyy。(9)2、度量系数及其计算方法在导数的坐标变换公式中涉及到下列坐标变换系数:,,,xyxy。这些系数称为坐标变换公式(5)对应的度量系数(metrics)。我们看到,为了求解计算平面中的偏微分方程,如(9)式,必须确定度量系数(有时还包括,,,,,xxxyyyxxxyyy等)的离散值。那么,这些度量系数如何计算呢?由于一般情况下,我们只知道坐标变换关系(5)、(6)的离散表达式,度量系数一般也要通过有限差分方法近似计算。但是,直接构造,,,xyxy的差分近似是不容易的。以x为例,根据偏导数的意义,x为y保持不变时随x的变化,如图2所示,网格点P处的x的计算公式应为:5)QPxPQPxx。由于Q一般不是网格点,因此,QQx是未知的,只能通过插值方法确定。另一方面,我们可以定义逆变换(6)式的度量系数,,,xxyy。在贴体坐标系中,这些度量系数的有限差分离散非常简单。如果采用中心差分离散,有1,1,,,1,1,1,1,,,1,1,)2)2)2)2ijijijijijijijijijijijijxxxxxxyyyyyy。(10)这就提示我们,如果能够找到,,,xyxy和,,,xxyy之间的关系,我们就可以得到,,,xyxy等的计算方法。下面,我们推导二者之间的关系。由变换关系(5)式,有xyxyddxdyddxdy,写成矩阵形式xyxyddxddy。(11)由逆变换(6)式,有xxdxdyydyd。(12)(11)、(12)式中的矩阵称为正变换和逆变换的Jacobi矩阵。由(11)、(12)易知图2x的计算PQ61xyxyxxyy,(13)即1xyxyyxyxJ,(14)其中1xyxyxxJxyxyyy(15)为坐标变换的Jacobi行列式(jacobian)。因此,1111xxyyyJyJxJxJ。(16)根据(10)、(15)、(16)式,可以得到,,,xyxy的差分离散形式。如何计算,,,,,xxxyyyxxxyyy呢?考虑10xyxyxyxy。(17)根据(17)式,我们同样可以得到11xyyJxJ。现在,把(17)式分别对,求导:722()()()()()()2()0()()()()()xyxyxyxxxyxyyyxyxxxyyyxyxyxyxyxxxyxyyyxxyxyxyxyxxyyxyxxyyxyxyxyxyxyxxyy()0()()()()()()0()()yxxxyyyxyxyxyxyxxxyxyyyxyxxxyyyxyxyxxyxxxyxyyyxyxyxyxyxyxxyyxyxxxyxyyyxyxy22()()()()2()0yxyxxxyxyyyxyxxxyyyxyxyxyxyxxyyxyxxyyxy,上面四个关系中,只有三个是独立的,写成矩阵形式,有:2222()2()()()2()xyxxxyxyyyxyxyxxyyxxxyxyyyxyxxyyxy。(18)所以12222()2()()()2()xyxxxyxyyyxyxyxxyyxxxyxyyyxyxxyyxy,222222()2()1()[()()]()2()xyxxxyxyyyxyxyJyJyyJyJxyJxyxyJxyxyJxyxyJxJxxJxxy。(19)同理222222()2()1()[()()]()2()xyxxxyxyyyxyxyJyJyyJyJxyJxyxyJxyxyJxyxyJxJxxJxxy。(20)对(19)、(20)式的右端进行有限差分离散,就可以算出,,,,,xxxyyyxxxyyy。8三、任意曲线坐标系中流体力学方程组的守恒形式考虑直角坐标系中的二维守恒型Navier-Stokes方程:()()0vvtxyUFFGG。(21)利用一阶导数的变换公式,有()()()()0vvvvxyxytUFFGGFFGG。(22)(22)式称为(,)平面上Navier-Stokes方程的弱守恒形式。在(22)式两侧乘以JacobianJ,有()()()()0vvvvxyxyJJJJJtUFFGGFFGG。上面的方程中的各项可以进一步整理:()()(){[()()]}()()()()()(){[()()]}()()()(),vvxyvxvyvxvyvvxyvxvyvxvyJJttJJJJJJJJJJUUFFGGFFGGFFGGFFGGFFGGFFGG得(){[()()]}{[()()]}()[()()]()[()()]0vxvyvxvyvxxvyyJJJtJJJJUFFGGFFGGFFGG。(23)注意到[()()]()()0[()()]()()0xxyyJJyyJJxx,所以,(23)式可化简为(){[()()]}{[()()]}0vxvyvxvyJJJtUFFGGFFGG。(24)(24)式称为(,)平面上Navier-Stokes方程的强守恒形式,一般记为:()()()0vvtUFFGG,(25)9其中[][][][]xyxyvvxvyvvxvyJJJJJ
本文标题:高等计算流体力学讲义(1)
链接地址:https://www.777doc.com/doc-4684319 .html