流动传热与燃烧的数值计算5流场求解..ppt

第四章：流场的求解 传热与燃烧过程的求解是在已知流场的条件下进行的，而一般条件下速度场是不能作为已知条件给出的，必须通过如下的控制方程计算。 数值方法求解速度场的真正困难在于不知道压力场，而且确定压力场的途径似乎也不清楚。 用正确的压力场代入动量方程时，得到的速度场满足连续性方程，但是压力场不是由连续性方程直接给出的，直接求解由连续性方程和动量方程产生的离散化方程组得到u,v,w,p也是不可能的。 ： §4.1求解流场的困难（控制容积法、某些差分格式） §4.1.1压力梯度的表示 以一维问题为例： 假定控制容积面e和w位于相应结点的正中间 假定压力是分片线性分布的 ： §4.1.2连续性方程的表示 以一维问题为例： 假定控制容积面e和w位于相应结点的正中间 假定速度u是分片线性分布的 §4.2交错网格 §4.2.1交错网格 交错网格的基本思想是把速度分量安置在不同于所有其它变量所用的网格上。 即：计算速度分量时单独使用一套网格，计算压力和其它所有变量时使用另外一套网格，两组网格相互交错。 在交错网格中，速度分量是对位于控制容积面上的点计算的。 例如：u是在垂直于x方向的面上计算的,u的位置相对主结点仅在x方向错开；位于两相邻主结点在x方向的连线上，且必须在控制容积面上。 交错网格的优点： 1.对代表性控制容积，离散化的连续性方程中包含的是邻近的速度分量之差，可避免波形速度场也满足连续性方程。即，只有合理的速度场才可能被连续性方程接受。 2.两邻近结点之间的压力差现在成了这两个结点之间的速度分量的自然驱动力，使波形和棋盘形速度场不会再被感觉成均匀的压力场，也不会再作为可能的解出现。 3.穿过控制容积面的质量流率不用对有关速度进行插值就可以计算，为建立离散化方程提供了方便。 交错网格的缺点： 基于交错网格的计算机程序必须带有关于速度分量位置的下标和方向的信息，并且仍需进行一定的较为麻烦的插值，但这些麻烦和交错网格的好处相比是微不足道的。 §4.2.2交错网格中的动量方程 交错网格的使用确实使离散化的动量方程稍微不同于其它相对主结点计算速度分量u,v,w的离散化方程，但这并不是本质上的差别，只是细节上的差别，是对动量方程使用交错的控制容积造成的。 X方向动量方程的交错控制容积如图所示，它相对于围绕主结点P的正常控制容积在X方向是错开的。 ： §4.3 SIMPLE算法 §4.3.1压力和速度的校正 1.速度校正方程 2.压力校正方程 P结点的控制容积如图所示： ： §4.3.3 SIMPLE算法(Semi-Implicit Method for Pressure Linked Equations) ： 2.讨论 ： §4.4 SIMPLER算法(Revised) 将(1),(2),(3)代入连续性方程的积分式，可得到压力方程： 2.SIMPLER算法的求解步骤 3.讨论： （1）由压力校正方程能产生合理的速度场，由合理的速度场和压力方程能算出合理的压力场，因此，SIMPLER的收敛速度应该更快。 （2）如果给出的速度场恰好正确，SIMPLER中的压力方程将直接给出正确的压力场，不需要进一步迭代，而SIMPLE 仍要进行多次迭代。 （3）由于压力方程和压力校正方程之间的类似性，关于p’方程边界条件的讨论对压力方程也适用。 （4）虽然SIMPLER比SIMPLE收敛的快，但是SIMPLER的每次迭代包括了更多的计算工作量。 §4.5 SIMPLEC算法(Consistent) ： ? ： ： * * P E W w e u x y P u v E e P E W S N *