NS方程求解研讨.pptVIP

* * 流过平板的超音速流动 物理模型 超音速平板绕流的类型 10.3.3 有限差分方程 向量形式的控制方程可以表示成： 利用泰勒级数展开，每个网点（i，j）处的流动变量都可以 推进到下一个时间步 其中，U还是表示流场变量，而且在时刻t是已知的，我们定义： 为了求 ，使求解能够推进，采用下列步骤： 由t时刻的已知的流动，用向前空间差分计算控制方程的右 端项，从而得到 的值 用步骤1）的结果，可以得到（t+△t）时刻流场变量的预 估值 合并步骤1），2），预估值的表达式为： 3）将上预估值代入控制方程右端并利用向后空间差分计算， 可以得到预估值的时间导数 4）将上得到的时间导数代入方程（10-6），得到（t+△t）时 刻U的校正量（具有二阶精度），合并步骤3），4），有： 重复1）至4 ）步骤，直到流场变量达到一个定常值，就是所求 的定常解。 为了保持二阶精度，E中关于x的导数项用于 方向相 反的差分离散，而关于y的导数项用中心差分近似。同样的， F中关于y的导数项用于 方向相反的差分离散，而F 中关于x的导数项用中心差分近似。 每个预估步和校正步结束后，由向量U可以得到原始变量 当 确定以后，可以通过10.3.1小节的方程求得其 余的流场变量 μ 和k都是温度T的函数，应用萨瑟兰（Sutherland）公式可 以得到μ，然后由普朗特数为常数的假设可以直接得到k 计算网格为70×70，下面的记号代表流向的网格 （入流边界） （出流边界） 平板的长度已知（LHORI）,所以x方向的步长为 同理 描述了垂直于平板表面的网格。 为了获得准确的解，激波必须位于计算区域之内，根据 布拉休斯的计算，假设计算区域的高度至少为五倍的边 界层厚度式合理的，可以满足计算的要求（图10-4）， 因此计算区域的垂直高度（LVERT）为 10.3.4 空间步长和时间步长的计算 其中边界层厚度由下式给出 因此y方向的步长为 平板长0.00001m，采用上述网格尺寸，x方向与y方向 网格的大小分别为 怎么判断网格的尺寸是否合适呢？ 在每一个时间步和每一网格点，x方向和y方向的网格雷诺数为 网格雷诺数的量级可以衡量计算网格的尺寸是否合适，对于 现在的问题网格雷诺数满足 注：流场变量在垂直平板方向有更大的梯度，因此在y方向的网 格雷诺数要求更高。为了准确地刻画流场，尤其是壁面附近的 流场，通常在垂直壁面方向需要更多的网格点。 其中 为当地音速，单位为m/s 注：K 为柯朗数，它起着调节因子的作用，使解保持稳定，通常取0.5《K《0.8 。 麦考马克方法是显式格式，时间步长受稳定性的限制，可以利用下面形式的CFL条件来确定时间步长 10.3.5 初始条件和边界条件 我们求解的是时间一阶，空间二阶的偏微分方程组，因 此需要给定（速度和温度）初始条件和边界条件。 求解是从初始条件开始推进的，因此需要给定t=0时每个 网格点（i，j）上的流场参数。每个网格点上流场参数的初始 条件，可以取相应的自由来流值。但是在物面上应满足无滑移 的边界条件，并且给定壁面温度TW的值，即 给定初始条件，方程就可以沿时间推进到定常解，推进过 程中，在计算区域的边界上必须指定边界条件。边界条件有四 种类型。 第一类边界条件，在前缘（imin，jmin）也就是（1，1），速 度满足无滑移条件（u（1，1）=v ( 1,1) =0.0）,而温度（T （1，1） ）和 压力(p （1，1）)区自由来流值。 第二类边界条件：在区域的左边界（不包括前缘点）和上边界， 速度的X方向分量u，温度，压力分别取相应的自由来流值，速 度的y方向分量v为零。 第三类边界条件：在平板表面，速度满足无滑移条件，并且假 设温度与壁面温度相同。利用物面上方网格点（j=2和j=3）处 的值，采用外推法计算壁面压力（不包括前缘点），表达式为 第四类边界条件：在区域的右边界，通过j相同的两个内部节点 外推得到边界处的所有参数（不包含jmin=1和jmax=70）。例如， u的表达式如下 上面详细说明了等温壁面的边界条件。CFD最大的优点 就是很容易的改变自由来流或者边界条件，并考察流场有哪 些变化。通过数值试验，可以从物理上更好的解释流动参数 的含义。因此，根据进一步