飞行器流动仿真讲稿第6章-计算流体力学的基本方法详细分析.ppt

第6.8节 有限体积方法简介 无论哪一种网格系统，Uj都只代表U在控制体内部的某种平均值； 于是，有限体积法一般形式的物理意义是：控制体内部U的平均值随时间的变化是该控制体与相邻控制体通过界面交换通量和控制体内部源项作用的结果。如果通量在某一侧边对该控制体的贡献为正，则它对同一界面的另一个控制体的贡献必为负，且在数值上是相等的。 四、界面通量计算 为了计算界面上的通量，需要规定物理量U及其偏导数关于时间和空间的局部分布形式（型线或插值公式）； 第6.8节 有限体积方法简介 常用的分段线性分布和阶梯式分布型线示意图 U随空间坐标变化的两种型线 U随时间坐标变化的两种型线 U U U U x x t t W w P e E W w P e E t t+?t t t+?t 分段线性分布 阶梯式分布 第6.8节 有限体积方法简介 型线选取仅是一种手段，目的是推导离散方程，一旦离散方程建立起来，型线便不再具有任何意义； 因此，选取型线时主要考虑方便以及最终离散方程的数值特性，而不必追求一致性：同一个控制方程中不同的物理量、同一物理量对不同的坐标可以用不同的型线，甚至同一物理量在方程不同项中对同一坐标也可以用不同的型线； 但是，选择不同的型线将会得到不同的差分格式，且对离散方程的求解方法和计算结果也有很大影响。 第6.8节 有限体积方法简介 五、有限体积法应用举例 以均匀网格为例，对控制体[w，e]，推导守恒型一维模型方程 的有限体积法离散方程。 i-1 i i+1 (?x)e (?x)w x Δx W w P E e 单位质量的某物理量 为简洁起见，假设为不可压缩 控制体 i 第6.8节 有限体积方法简介 对空间和时间积分，有 通过型线对各项 分别进行处理 第1项—非定常项：变量φ沿x的变化取阶梯式分布，在控制体中均等于节点P的值φP（不同控制体有不同的值） i-1 i i+1 (?x)e (?x)w x Δx W w P E e 第6.8节 有限体积方法简介 第2项—对流项：变量φ随时间t的变化取显式阶梯分布，即在[t，t+Δt]时间内， 取t时刻的值φt φ随x的变化用分段线性表示 i-1 i i+1 (?x)e (?x)w x Δx W w P E e 第6.8节 有限体积方法简介 第3项—扩散项：一阶偏导数随t的变化取阶梯式分布 φ随x的变化也用分段线性表示 第4项—源项：通常，源项随x和t的变化均取阶梯分布 t时刻控制体中Q的平均值 第6.8节 有限体积方法简介 将各式代入有限体积法离散方程，整理得最终离散方程 t时间层 t+Δt时间层 观察上式可以发现，非定常项的离散实际就是一阶向前差分，而对流项则为二阶中心差分； 与有限差分方法相比，除了方程形式（守恒型与非守恒型）不同，两者的差分格式是相同的，而这正是选择特定型线的结果。 第6.9节 CFD编程与计算数据的处理 用CFD研究流动问题的一般步骤可以表述为三步曲： 第一步——提出物理模型，选择适当的数值格式和算法，网格生成； 第二步——编程调试与计算，边界条件处理； 第三步——计算数据处理与分析。 第6章 计算流体力学的基本方法 第6.9节 CFD编程与计算数据的处理 一、物理模型 在物理模型选择方面主要有如下一些考虑： 合理确定维数。解决实际问题需要一维、二维还是三维，与求解精度、关心的主要问题和计算条件等有关。维数增加将使问题复杂化，并大大提高对计算机能力的要求； 时间因素。对非定常过程，前一时刻发生的变化可以影响后来的变化，反之则不然，即时间影响是单向的，影响程度除与过程本身的非定常特性有关外，还与计算中能选取的时间步长有关，一般随时间步长增大而减弱； 对流与扩散因素。对流对上、下游的影响与马赫数有关，超声速时影响是单向的，只影响下游，亚声速时却是双向的；扩散（包括热传导、黏性等）在所有方向上都有影响。上、下游之间究竟是相互影响还是单向影响，将直接关系到差分格式的选择； 第6.9节 CFD编程与计算数据的处理 压强因素。对流动的影响与当地马赫数有关，因为压强变化以声速向空间各方向传播。对亚声速流影响是相互的；而对超声速流，影响则从上游向下游传播，是单向的； 边界条件。分析边界条件影响有时很重要。例如考虑无回流管流时，若给定进口质量流率求进出口压差，则影响基本上是单向的，可以推进求解；若给定进出口压差求质量流率，则出口压强直接影响质量流率大小，影响是相互的，必须迭代求解。这与控制方程分为三种类型直接相关：若影响是相互的则为椭圆型；若影响仅向下游传播则为抛物型；若影响仅在下游某个劈锥或楔范围内传播，则属于双曲型，锥或劈顶角取决于当地马赫数。微分方程在不同坐标方向上也可以分为上述三种类型。 第6.9节 CFD编程与计算数据的处理 二、编程方面的考虑 编制计算机程序