计算水动力学重点.doc

◎ 非守恒型控制方程的推导： 几点说明： 从微元体的角度看，控制方程的守恒型与非守恒型是等价的，都是物理守恒定律的数学表示。但对有限大小的计算体积，两个形式是有区别的。 ??不论节点布置的疏密程度如何，根据控制方程导出的离散方程也具有对任意大小容积守恒的特性。 ??非守恒型方程便于对离散方程进行理论分析，而守恒型控制方程能保持物理量守恒的性质，便于克服对流项非线性引起的问题。 第七章 双边界法 第八章 SIMPLE算法自1972年问世以来，在计算流体力学及传热学中得到了广泛的应用，同时也得到了不断的改进与发展。近年来，SIMPLE改进算法的研究成果主要有：Patankar于1980年提出的SIMPLER算法，Spalding于1981年提出的SIMPLEST算法和Doormal与Raithby于1984年提出的SIMPLEC算法[29]等。 （一）（一）SIMPLER算法 由SIMPLE算法得出的值对修正速度而言是相当好的，但对修正压力则是过分了。虽然对采用了亚松驰处理，也未必能恰到好处。这样就使速度场的改进与压力场的改进不能较好的同步进行，最终影响了整个流场的迭代收敛速度。于是就产生了这样的想法：只用来修正速度，压力场的改进则另谋更合适的方法。此外，在SIMPLE算法中，为了确定动量离散方程的系数，一开使就假设了一个速度分布，那么与这一速度分布相协调的压力场即可由动量方程计算而得，不必在单独假定一个压力场。把以上两个思想结合起来，就构成了SIMPLER算法。 在SIMPLER算法中，初始的压力场和速度场是相协调的，算出的压力场不必亚松驰，使迭代层次数减少。但另一方面，每一层次计算中所花的时间则较SIMPLE多，因为SIMPLER算法中要多解一个Poisson方程。但就总的计算时间来说，SIMPLER算法一般较SIMPLE少。 （二）（二）SIMPLEST算法 与SIMPLE算法相比，它主要有以下两个特点： （1）对流项采用迎风格式，因为这是一个绝对稳定的格式，且扩散项与对流项的影响系数可以分离开来，不像指数（或乘方）格式那样综合在一起。至于由迎风差分所引起的假扩散问题，则采取逐步加密网格，以获得与网格疏密程度无关的解这种方法加以克服。 （2）把邻点的影响系数表示成对流分量及扩散分量之和，并把对流部分全部归入源项，于是对的动量方程为： ????? （3-34） 由此可见，当扩散项忽略不计时，动量方程实际上采用了Jacobi的点迭代求解。因此，在这种算法中扩散项采用线迭代而对流项采用点迭代。点迭代的收敛速度是比较慢的，但是由于对流项与压力之间的耦合关系等原因，正希望利用这一特性以防止迭代发散。这种混合格式的计算方法有利于促进强烈非线性问题的迭代过程收敛。SIMPLEST算法的计算步骤与SIMPLE算法基本相同。 （三）（三）SIMPLEC算法 在SIMPLE算法中，为了求解的方便，略去了速度修正值方程中的及项，从而把速度的修正完全归结为有压力差项的直接作用所致。当在略去等项不计时，实际上犯了一个“不协调一致”的错误，因为略去相当于使，而根据系数计算公式，，在令等号后面的时，没有同时令等号前的也等于零。为了能略去等项而同时又能使方程基本协调，在方程的等号两端同时减去，即： ?????? ??????（3-35） 可以预期，与其邻点的修正值具有相同的量级，因而略去所产生的影响较在原方程中不计所带来的影响要小得多。于是得： ???????????????????????????? （3-36a） 类似地， ???????????????????????????? （3-36b） ?? ?这就是协调一致的SIMPLE算法，简称SIMPLEC。SIMPLEC算法的收敛特性远优于SIMPLE，有时甚至优于前面介绍的SIMPLER。从实质上来说，它是对SIMPLE的校正公式稍作改进，改善其协调性，而计算步骤是完全一样的，从而工作量增加较少。 第六章 小结 1、有限分析法FAM (1)基本思想：把求解区域划分成许多矩形网格，四个网格组成一个单元；在局部单元内将微分方程线性化，在单元边界上为一近似函数；局部单元内求解微分方程的解析解；建立单元中心的和其周围八个结点之间的迭代关系式。 (2)格式推导： 对象：对流扩散方程----标准的椭圆型方程：2Afx +2Bfy =fxx +fyy +G 步骤：1)、线性化 2)、齐次化 (3)优、缺点 优点： 1)、稳定性好。有限分析系数的大小取决于单元雷诺数Ah、Bk及网格长宽比k/h=b；有限分析系数CSC、CNC、CWC、CEC、CSW、CSE、CNW、CNE 及P1、P2、E2的值都在0和1之间；八个系数之和等于1。这一特性能保证迭代格式具有较好的稳定性。 2)、对称性好。由