计算流体力学汇编.pptVIP

物理模型是指把实际的问题，通过相关的物理定律概括和抽象出来并满足实际情况的物理表征。比如，我们研究管道内的流体流动，抽象出来一个直管，和粘性流体模型，或者我们认为管道内的液体是没有粘性的，使用一个直管和无粘流体模型.还有，我们根据热传导定律，认为固体的热流率是温度梯度的线形函数，相应的傅立叶定律就是导热问题的物理模型。因此，不难理解物理模型是对实际问题的抽象概念，对实际问题的一种描述方式，这种抽象包括了实际问题的几何模型，时间尺度，以及相应的物理规律。数学模型就好理解了，就是对物理模型的数学描写。比如N-S方程就是对粘性流体动力学的一种数学描写，值得注意的是，数学模型对物理模型的描写也要通过抽象，简化的过程。;建立控制方程;确定边界条件与初始条件 初始条件与边界条件是控制方程有确定解的前提，控制方程与相应的初始条件、边界条件的组合构成对一个物理过程完整的数学描述。 初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况。对于瞬态问题，必须给定初始条件。对于稳态问题，不需要初始条件。 边界条件是在求解区域的边界上所求解的变量或其导数随地点和时间的变化规律。对于任何问题，都需要给定边界条件。例如，在锥管内的流动，在锥管进口断面上，我们可给定速度、压力沿半径方向的分布，而在管壁上，对速度取无滑移边界条件。 对于初始条件和边界条件的处理，直接影响计算结果的精度。 ;划分计算网格 采用数值方法求解控制方程时，都是想办法将控制方程在空间区域上进行离散，然后求解得到的离散方程组。要想在空间域上离散控制方程，必须使用网格 不同的问题采用不同数值解法时，所需要的网格形式是有一定区别的，但生成网格的方法基本是一致的 ; 在同一种离散化方法中，如在有限体积法中，对式中的对流项所采用的离散格式不同，也将导致最终有不向形式的离散方程。 对于瞬态问题，除了在空间域上的离散外，还要涉及在时间域上的离散。 ;建立离散方程 对于在求解域内所建立的偏微分方程，理论上是有真解（或称精确解或解析解）的。但由于所处理的问题自身的复杂性，一般很难获得方程的真解。因此，就需要通过数值方法把计算域内有限数量位置(网格节点或网格中心点)上的因变量值当作基本未知量来处理，从而建立一组关于这些未知量的代数方程组，然后通过求解代数方程组来得到这些节点值，而计算域内其他位置上的值则根据节点位置上的值来确定。由于所引入的应变量在节点之间的分布假设及推导离散化方程的方法不同，就形成了有限差分法、有限元法、有限元体积法等不同类型的离散化方法。 ;离散初始条件和边界条件 前面所给定的初始条件和边界条件是连续性的，如在静止壁面上速度为0，现在需要针对所生成的网格，将连续型的初始条件和边界条件转化为特定节点上的值，如静止壁面上共有90个节点，则这些节点上的速度值应均设为0。这样，连同在各节点处所建立的离散的控制方程，才能对方程组进行求解。 。 ;给定求解控制参数 在离散空间上建立了离散化的代数方程组，并施加离散化的初始条件和边界条件后，还需要给定流体的物理参数和紊流模型的经验系数等。此外，还要给定迭代计算的控制精度、瞬态问题的时间步长和输出频率等。 ;求解离散方程 在进行了上述设置后，生成了具有定解条件的代数方程组。对于这些方程组，数学上已有相应的解法，如线性方程组可采用Guass消去法或Guass-Seidel迭代法求解，而对非线性方程组，可采用Newton-Raphson方法。;;显示和输出计算结果 线值图：在二维或三维空间上，将横坐标取为空间长度或时间历程，将纵坐标取为某一物理量，然后用光滑曲线或曲面在坐标系内绘制出某一物理量沿空间或时间的变化情况。 矢量图：直接给出二维或三维空间里矢量（如速度）的方向及大小，一般用不同颜色和长度的箭头表示速度矢量。矢量图可以比较容易地让用户发现其中存在的旋涡区。 等值线图：用不同颜色的线条表示相等物理量(如温度)的一条线。; 对流-扩散方程的混合格式及乘方格式;中心差分与一阶迎风格式的比较;二、混合格式(Spalding,1971);三、指数格式;三、指数格式;四、乘方格式;五、5种3点格式系数汇总;1、二阶迎风格式（SUD）;3、QUICK格式;提出SIMPLE算法的缘由;SIMPLE算法的假设条件;SIMPLE算法的计算步骤;湍流/紊流;湍流是一种高度复杂的三维非稳态、带旋转的不规则流动。湍流流体的各种物理参数