计算流体CFD2004-4.pptVIP

第四章 有限差分方程的解法 模型方程: 4.1 椭圆方程 4.1.2 迭代法 4.2 抛物线方程 Von Neumann稳定性分析 4.2.2 隐式格式 4.2.3 交错方向隐式 (ADI) 格式 4.2.4 近似因子分解 4.3 双曲方程 4.3.2 隐格式 4.3.4 非线性问题 4.4 Burgers方程 Von Neumann稳定性分析 4.4.2 Runge-Kutta 法 4.5 代数方程的解法和误差来源 扩散输运方程 对流输运方程 4.6 任意几何形状的坐标变换 * 华中科技大学能源与动力工程学院 计算方法简介： 1、有限差分法(Finite difference method，FDM) 偏微分方程差分，最古老、最普遍的方法； 十分简便而有效，很容易引入对流项的高阶格式 离散方程的守恒特性难以保证，不规则区域的适应性差 2、有限容积法(Finite volume method，FVM) 积分守恒型控制方程的离散，是目前应用最普遍的方法； 具有守恒性，对区域形状的适应性较好， 3、有限元法(Finite element method，FEM) 对积分控制方程离散，流体计算不如有限容积法发展成熟 不规则几何区域的适应性好， 4、有限分析法(Finite analytic method，FAM) 80年代初发展起来的一种数值方法，可克服高Reynolds 数下数值解容易发散或振荡的缺点， 计算工作量较大，对计算区域几何形状的适应性也较差 5、边界元法(Boudary clcmcnt mcthod。BEM) 求解问题的空间维数降低一阶，从而使计算的工作量及所需计算机容量大大减小。 边界元法最大限制是：需要已知所求解偏微分方程的格林函数基本解。 6、谱分析方法(SPcctrZllHlcthod，SM) 被求解的函数用有限项的级数表示，目前用得不普遍 优点是可以获得高精度的解，但不宜编制通用程序 7、数值积分变换法(Integraltransformation method，ITM) 对不具备分析解的非线性偏微分方程，把它的解表示成一个特征值问题的解及一个降维的定解问题解的组合 优点是计算时精度可以较高，但不容易形成通用程序，特征值问题的选取有任意性，对强非线性问题，计算量较大。 8、格子—Boltzmann方法(lattice—Boltzman method，LBM) 基于分子运动论的一种模拟流体流动的数值方法，不再基于连续介质的假设，而是把流体看成是许多只有质量没有体积的微粒所组成，这些微粒可以向空间若干个方向任意运动。通过其质量、动量守恒的原理，建立起表征质点在给定的时刻位于空间某一个位置附近的概率密度函数。再通过统计的方法来获得质点微粒的概率密度分布函数与宏观运动参数间的关系 工程问题的控制方程是Euler方程及N.S.方程。为了反映实际流场，方程是多维的，边界也较复杂。为了弄清数值计算的一些机理性问题，提高计算效率，验证新格式的效果，可先对简化的模型方程进行研究，然后推广到一般问题中去。 解法：针对离散后代数方程组 显式与隐式 迭代解与直接解 Euler方程组(非定常)的模型方程 对N．S．方程组，最简单模型方程是Burgers方程 稳定热传导等的模型方程是泊桑方程或拉普拉斯方程 其中：β=Δx/Δy 令 γ= - 2(1+ β2 ) 5点格式： 5 4 3 2 j=1 u=u4 i=1 2 3 4 5 u=u3 u=u2 u=u1 0 0 0 0 两种解法：迭代解法4.1.2，直接解法4.1.3 直接解法比迭代的解法要耗费更多的时间，容易受到舍入误差的影响；在迭代解法中，每一步产生的误差都会在下一步得到修正，因此，通常不用担心舍入误差。 显式点迭代 隐式线迭代 迭代因子？ 边界值的传递？ 差分方程： (FTCS) 利用Von Neumann稳定性分析，由放大因子|g|?1得到稳定条件： 0?d=αΔt/Δx2 ?1/2； 4.2.1 显式 用分析的办法提出对稳定 性的验证是一个分析差分格式可用性的有效途径。有许多分析稳定性的办法，如能量法、 离散扰动法、矩阵法和Fourier法、Hirt 法等。但应用最广的是Von Neumann提出的Fourier分析 法。 1、 Von Neumann的稳定性分析法 微分方程 离散方程（FTCS） 富氏级数展开： An是波数为kx