计算水力学02剖析.ppt

第二章 离散方法 预测自然界和水利工程中水流和输运现象的问题，可归结为求解形式不同的微分方程组。按照现象本身的复杂程度和工程实践所要求的精细程度不同，待解的微分方程组可能具有各种不同的形式和不同的组合，但一般说来，水流输运问题中遇到的微分方程，总可以化作统一的形式，通用微分方程．本章将研究求解通用微分方程的数值计算方法。 电子计算机问世以来，数值计算技术发展异常迅猛，计算方法千姿百态、纷纭复杂。本章试以加权余量法(Method of Weighted Residuals)为总纲，分析一些主要的计算方法的出发点和基本原理，了解目前比较流行的计算方法(如有限差分法，有限单元法，边界元法等)在计算方法系统中所处的位置及其相互联系，理解各类计算方法的实质。本章的后面部分，着重介绍本课程中采用的计算方法－－有限体积法的基本原理和实现原则。 §2—1 加权余量法 一、加权余量法的基本思想 设有区域 中的微分方程 边界法、区域法和混合法 在区域法中， 可将近似解写为 试函数 应为完整函数组中的线性独立函数 ； 还应具有必要的连续阶数，才可保证余量R不为零。 将试函数代入微分方程，必然产生非零余量R 二、子区域法(Subdomain Method) 例1 求解方程 边界条件为 设 易知其满足边界条件 . 作为第一次近似，取整个计算区域oxl作为唯一的子区域 可解出 作为第二次近似，取两个子区域：ox0.5和0x1，可得 解得 二阶近似解为 原方程的精确解为 子区域法计算结果 三、配置法(Collocation Method) 配置法的效果与有限差分法一致，都是使求解的微分方程在计算区域的若干点上得到满足，而不计及因变量在这些点之间的变化．从这个意义说，有限差分法可解释为没有试函数的配置法。 四、伽辽金法(Galerkin Method) 五、矩法(Method of Moments) 六、最小二乘法(Least-squares Method) §2—2 各类计算方法的联系和比较 这样，总体近似被单元内的近似所代替，总体近似函数被单元内的插值函数(又称形函数)所代替，原先不具有明确意义的未知系数ai被单元节点上因变量的未知数值所代替。这种近似称为局部近似(Local Approximation)，这种方法称为离散方法。 二、各类加权余量法的比较 三、有限差分法、有限单元法和谱方法的比较  §2-3 有限体积法 物理的合理性是指数值解答应与精确解具有同样的定性变化趋势。 物理上合理的和不合理的解答 §2-4 有限体积法的四条基本原则 原则1——控制体积交界面的一致性 当一个表面为相邻的两个控制体积所共有时，在这两个控制体积的离散方程中，通过该表面的通量的表达式必须相同。 原则2--正系数 在大多数水流输运问题中，节点上因变量的数值只通过对流过程和扩散过程受到相邻节点的影响，所以，当其它条件不变时，一个节点上数值的增加，必引起相邻节点数值的增加而不是减少。对于通用方程，可知相邻节点系数和中心节点系数必定同号。系数的数值可以全为正值或全为负值，不妨规定离散方程的系数皆为正值。原则2可叙述为，中心节点系数和相邻节点系数必须恒为正值。 原则3－源项的负坡线性化 原则4－相邻节点的系数之和 李光炽 计算水力学 李光炽 计算水力学 李光炽 计算水力学 李光炽 计算水力学 将权函数取为余量R对未知系数ai的偏导数,就是最小二乘法 其效果等价于令余量R的内积取得极小值： 正是从这个意义考虑，此法被称为最小二乘法。 李光炽 计算水力学 一、谱方法和离散方法 在用加权余量法求出未知系数ai并代入近似解式求出计算区域中任意一点的因变量的数值。如此安排试函数和未知函数，给出的是总体近似(Global Approximation)。如果试函数恰取为所求解方程的本征函数，则这种总体近似的方法与数学物理方法中介绍的分离变量法相一致。采用总体近似的各类方法，又称为谱方法(Spectral Method)。 李光炽 计算水力学 所关心的问题往往不是因变量在整个计算区域中的分布，而是因变量在空间若干特定位的数值。将因变量在给定点的数值直接作为未知系数ai，并求解这些数值，作为满足实际需要的解答，正是一切离散方法的出发点。为了用因变量在特定点的数值ai来描述因变量在整个计算区域的分布，我们可将计算区域划分为许多子区域即单元，假设因变量在单元内的分布规律，据此将因变量在单元内的分布描述为单元节点上因变量数值的函数。 李光炽 计算水力学 李光炽 计算水力学 李光炽 计算水力学 伽辽金法精度高，应用范围广，最小二乘法用于求解椭圆问题，精度与伽辽金法