2024年-CFD程序设计PPT（精编）.ppt

6.1制订一个特殊程序的必要性我们推导了已知流场通用微分方程的公式，但是除了某些特殊的情况之外，流场是不可能已知的，我们必须由相应的控制方程计算局部的速度场与密度场。速度分量受动量方程控制，而动量方程则是对的通用微分方程的特殊情况〔取，等）。于是从一个估计的速度场开始，我们可以迭代求解动量方程，从而得到速度分量的收敛解。6.1-1：主要的困难*Companyname6.1制订一个特殊程序的必要性计算速度场的真正困难在于：未知压力场，压力梯度构成动量方程中源项的一部分，然而还没有一个可以用来求得压力场的明显方程；压力场间接地通过连续性方程规定，当正确的压力场代入动量方程时，所得到的速度场满足连续性方程。但是对于我们所要达到的目的而言，这种间接的规定是不太有用的，除非我们打算去直接求解由动量方程与连续性方程所推得的整个离散化方程组。6.1-1：主要的困难*Companyname6.1制订一个特殊程序的必要性6.1-2：以涡量为基础的方法为了解决确定压力所带来的困难人们提出若干从控制方程中消去压力的方法。涡量法：在二维的问题中，通过交叉微分从两个动量方程中消去压力，就可以导出一个涡量输运方程，把涡量的定义与二维稳态情况下的流函数的定义相结合就是众所周知的“流函数、涡量法”。*Companyname6.1制订一个特殊程序的必要性6.1-2：以涡量为基础的方法流函数／涡量法具有某些吸引人的优点：压力消失了，代替处理连续性方程和两个动量方程的问题，解两个方程可求得流函数和涡量；某些边界条件可以很容易地规定：当一个无旋的外流位于计算域的邻近时，可以方便地取边界的涡量为0。比较突出的缺点：在壁面上的涡量值很难给定，而这往往成为取得收敛解的麻烦的原因所在。*Companyname6.1制订一个特殊程序的必要性6.1-2：以涡量为基础的方法压力，这个量已经被消去了。对于三维的问题，一个以涡量为基础的方法采用六个因变量，即三个涡量向量和三个速度势向量分量，其复杂性实际上超过直接处理三个速度分量和压力本身。此外，涡量向量与速度势向量所包含的概念要比速度分量与压力的意义难以理解和接受。寻求一种采用所谓的基本变量（即速度分量和压力）的方法。*Companyname6.2若干有关的困难6.2-1：压力梯度项的表述对图6.1所示的一维情况进行构造x方向的动量方程的离散化公式，那么唯一的新问题是项的表示式，对整个控制容积的积分表达式。*Companyname6.2若干有关的困难项对离散化方程的贡献是压力降Pw–Pe，它是作用在单位面积的控制容积上的净压力。假设压力是分片线性分布的，如果把控制容积面e和w选在相应网格点之间的中点，就有：6.2-1：压力梯度项的表述*Companyname6.2若干有关的困难这意味着动量方程将包含有两个相间而不是相邻网格点之间的压力差。即压力取的是较实际所采用的网格为粗的网格上的值。这一点导致解的精度下降。由图6.2看出，假设压力场用网格点处的压力值表示，不能认为这样一个锯齿形的压力场是合乎实际的。6.2-1：压力梯度项的表述*Companyname6.2若干有关的困难动量方程对这样一个波形压力场的“感受”竟然与对均匀的应力场的“感受”一样。在二维的情况下，可以更为明显地看清这个困难。x方向的动量受Pw–Pe的影响，y方向的动量则受Ps–Pn的影响；可以得出结论：在图6.3中所示的压力场由布置成棋盘形的四个任意的压力值所组成，它不会在x或y方向产生压力合力。于是一个高度不均匀的压力场就被动量方程的这种特殊离散形式处理成为一个均匀的压力场了。6.2-1：压力梯度项的表述*Companyname6.2若干有关的困难6.2-1：压力梯度项的表述如果得到一个相当光滑的压力场的解，那么在该解上附加上一个棋盘型的压力场就可以构成任意数量的附加解。由于棋盘场包含着没有压力合力的作用的意思，动量方程将不会受该附加压力场的影响，得出这种荒谬的解的数值方法是不理想的。*Companyname6.2若干有关的困难对于一维稳态的常密度流场，连续性方程简化为：如果我们在图6.1所示的控制容积内积分这一方程，就有：6.2-2连续方程的表述*Companyname6.2若干有关的困难与前面对Ｐ一样，对u采用分段线性分布并取控制容积面于中点位置就导致：离散化的连续性方程要求相间