计算流体力学结课.docVIP

计算流体力学结课作业 晴哥哥有其它需要，请加微信： 二、概述3种数值方法，介绍其基本思想、技术要点、应用步骤，并就几种方法的特点和区别做比较分析。 （1）有限差分法（Finite Difference Method，FD）是偏微分方程（partial differential equation简称PDE）数值求解的最为古老的方法，是Euler在18世纪提出的，主要用来求解与时间有关的问题，该方法概念清晰，方法简单、直观，是应用最多的一种数值方法。 有限容积法（Finite Volume Method）又称为控制体积法，有限容积法（FVM）是计算流体力学（CFD）和计算传热学（NHT）中应用应用广泛的数值离散方法。 基本思想：为了求解由偏微分方程定解问题所构造的数学模型，有限差分法将定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程，即有限差分方程组，解此方程组就可以得到原问题在离散点商的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 有限容积法式将计算区域划分为网格，并使每个网格点周围有一个互不重复的控制体积；将待解的偏微分方程对每一个控制体积积分，从而得出一组离散方程，其中的未知量是网格点的特征变量。为了求出控制体积的积分，必须假定特征变量值的网格点之间的变化规律。从积分区域的选取方法看来，有限体积法属于加权余量法中的子域法；从未知解的近似方法看来，有限体积法属于采用近似的离散方法。 （2）技术要点： 有限差分方法重要的技术就是将偏微分方程转换成差分方程组，然而构造这样的方程组有很多方法，我的目的就是寻找一个最接近偏微分方程的解的差分方程。构造差分方程就要知道怎么构造每一项的差分形式。 有限差分法的主要内容包括：如何根据问题的特点将定解区域做网格剖分；如何把 原微分离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算n+1层的近似值时要用到第n层的近似值，直到与初始值有关。前面各层若有舍入误差，必然影响到后面隔层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有一下三种方法。最常用的方法是数值微分方法，比如用差商代替微商等等。另一种方法叫差值法，以为在实际问题中得出的微分方程常常反应物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式；从差分的空间形式来考虑，可分为中心差分格式和逆风差分格式；考虑时间因子的影响，差分格式还可以分为显式格式、隐式格式、显隐交替格式等；目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式；差分方法主要适用于结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等等，其中前两种格式为一阶计算精度，后两种为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限体积法的基本思路易于理解，并能得出直接的物理解释。离散方程的物理意义，就是因变量在有限大小的控制体积中的守恒原理，如同微分方程表示因变量在无限小的控制体积中的守恒原理一样。有限体积法得到的离散方程，要求因变量的积分守恒对任意一组控制体积都得到满足，对整个计算区域，自然也得到满足。这是有限体积法吸引人的优点。有一些离散方法，例如有限差分法，仅当网格及其密时，离散方程才能满足积分守恒；而有限体积法即使在粗网格的情况下，也显示出准确的积分守恒。有限体积法和有限差分方法一样，也需要对计算域进行离散，将其分割成有限大小的离散网格。区域离散的实质是用有限个离散点来代替原来的连续空间，即生成计