《计算流体力学》结课作业.docVIP

2012～20学年第学期 是研究流体(液体和气体)的力学运动规律及其应用的学科的研究，、飞行器和核电站的设计及其运行，可燃气体或，都广泛地用到流体力学知识。许多现代所关心的问题既受流体力学的指导，同时也促进了流体力学不断发展。1950年后计算机的发展给予流体力学极大的推动。 同物理学、化学等学科一样，流体力学离不开实验，尤其是对新的流体运动现象的研究。实验能显示运动特点及其主要趋势，有助于形成概念，检验理论的正确性。二百年来流体力学发展史中每一项重大进展都离不开实验。: 几何相似、运动相似、动力相似、边界条件、初始条件相似。两个流场动力相似，则两个流场所有的动力相似准则应分别相等。但要做到两个独立的动力相似准则同时分别与原型的同名准则相等是不可能的，所以只能部分相似，即近似模型实验。模拟实验在流体力学中占有重要地位。根据模型实验所得的数据可以用像换算单位制那样的简单算法求出原型的数据。流体力学中最常用的基本模型有：连续介质、牛顿流体、不可压缩流体、理想流体、平面流动等。用数学语言将质量守恒、动量守恒、能量守恒等定律表达出来，从而得到连续性方程、动量方程和能量方程。此外，还要加上某些联系流动参量的关系式(例如状态方程)，或者其他方程流体力学基本方程组。求出方程组的解后，可结合具体流动，解释这些解的物理含义和流动机理。 数值研究的一般过程是：对流体力学数学方程作简化和数值离散化，编制程序作数值计算，将计算结果与实验结果比较。 在流体力学理论研究和工程应用中，描述流体运动的数学方程是非线性偏微分方程组，只对极少数的简化模型可以通过数学方法，获得理论分析解，多数情况下，只能通过数值计算的途径进行求解。这里说的“数值计算”，是指利用高速电子计算机，对描述流体力学具体问题的偏微分方程初边值问题进行离散化计算，从而获得流动区域中离散点上的流体物理量的求解方法。这种通过数值计算获得流动区域中离散点的数值解的方法，通常称为流体力学数值解法，也可称之为计算流体力学。随着高速电子计算机的发展与普及，数值方法越来越受到重视，已成为流体力学理论研究和工程应用的重要手段。计算流体力学是以计算机为工具、以流体力学的基本方程（纳维－斯托克斯方程）为理论依据，采用离散化的数值方法对流体力学问题进行数值模拟和分析的流体力学分支学科。在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了“计算流体力学”。解决流体力学问题时、实验、理论分析和数值是相辅相成的。实验需要理论指导，才能从分散的、表面上无联系的现象和实验数据中得出规律性的结论理论分析和数值也要依靠实验给出物理或数据，以建立流动的力学模型和数学模式；最后，还须依靠实验来检验这些模型的完善程度。此外，实际流动往往异常复杂，理论分析和数值计算会遇到巨大的数学和计算方面的困难，得不到具体结果，只能通过实验进行研究。数值模拟已经成为人类改造世界的第三种手段按时间步长和空间步长将时间和空间区域FTCS格式（时间向前差分、空间中心差分） FTFS格式（时间向前差分、空间向前差分） FTBS格式（时间向前差分、空间向后差分） 根据冯诺依曼稳定性分析，FTBS格式，在和的条件下稳定，而FTFS格式，在和的条件下稳定。这里，当的符号改变时，为了使差分格式稳定，空间差分的方向也做了相应的变化。 3、逼近求解，这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程。 求解差分方程组一般采用Gauss消去法、追赶法、迭代法、交替方向隐式差分法（ADI法）、隐式近似因式分解法（AF法）等，上述消去法和追赶法对求解离散后的代数方程组没有特别的优势，采用迭代法来求解方程组在收敛速度上有一定的优势。迭代法基本思路为：首先对求解的未知量给一个预测值，代入代数方程组，它一定不满足方程组。利用一些特性对预测值进行修正，并把修正后的预测值再代入方程组，它仍不满足方程组。再修正预测值，再代入方程组，通过不断迭代过程，直到收敛于数值解。 迭代法中的高斯—赛德尔迭代法，简称为G-S迭代法，具有形式简单，收敛速度较快的特点。假设求解过程是按x和y增长方向进行，于是在求点的值时，在和点上的值实际上已经求出。G-S迭代法基本思路是把已经求的的值，立即代入迭代式中去。它的迭代差分格式为： 松弛迭代法是对G-S迭代法的一种