流体传热及数值模拟第九讲.ppt

流体流动与传热过程数值模拟讲义之九 主要内容 常微分方程概念 数值求解方法 隐式法 显式法 欧拉法 （重点） 龙格-库塔法 （重点2） 应用举例 欧拉及R-K法编程求解 常微分方程与空间微分方程区别 对时间或演变过程因变量进行微分-推进或演变概念。 初值与边界条件 初值问题的Euler方法 * 问题 /* Numerical Methods for Ordinary Differential Equations */ 常微分方程数值解 例：对下图所示的二维非稳态金属块的无源导热问题，其初始温度为室温27度，四周边界温度为定值。求解金属块内的非稳态二维热传导。导热系数k=1000 w/(mK)。 控制方程： 边界条件： x Y 常微分方程提出 初始条件： t=0, T=27 定义 得到 常微分方程 （对时间微分） （对空间微分） 初值条件 t=0, T=27 常微分方程数值解 /* Numerical Methods for Ordinary Differential Equations */ ? 待求解的问题：一阶常微分方程的初值问题 /* Initial-Value Problem */: 解的存在唯一性（“常微分方程”理论）：只要 f (x, y) 在[a, b] ? R1 上连续，且关于 y 满足 Lipschitz 条件，即存在与 x, y 无关的常数 L 使 对任意定义在 [a, b] 上的 y1(x) 和 y2(x) 都成立，则上述IVP存在唯一解。 解析解法：（常微分方程理论） 只能求解极少一类常微分方程；实际中给定的问题不一定是解析表达式，而是函数表，无法用解析解法。 如何求解 计算解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 节点间距 为步长，通常采用等距节点，即取 hi = h (常数)。 数值解法: 求解所有的常微分方程 步进式：根据已知的或已求出的节点上的函数值计算当前节点上的函数值，一步一步向前推进。因此只需建立由已知的或已求出的节点上的函数值求当前节点函数值的递推公式即可。 --------Euler’s Method §1 欧拉方法 /* Euler’s Method */ §1 Euler’s Method Taylor展开法 几何意义 亦称为欧拉折线法 /* Euler’s polygonal arc method*/ 几何直观是帮助我们寻找解决一个问题的思路的好办法哦 定义 　　　在假设 yn = y(xn)，即第 n 步计算是精确的前提下，考虑公式或方法本身带来的误差: Rn = y(xn+1) ? yn+1 , 称为局部截断误差 /* local truncation error */。 说明 显然，这种近似有一定误差， 而且步长越大，误差越大， 如何估计这种误差y(xn+1) ? yn+1 ？ §1 Euler’s Method 截断误差: 实际上，y(xn) ? yn， yn 也有误差，它对yn+1的误差也有影响，见下图。但这里不考虑此误差的影响，仅考虑方法或公式本身带来的误差，因此称为方法误差或截断误差。 局部截断误差的分析：由于假设yn = y(xn) ，即yn准确，因此分析局部截断误差时将y(xn+1) 和 yn+1都用点xn上的信息来表示，工具：Taylor展开。 ? 欧拉法的局部截断误差： Rn+1 的主项 /* leading term */ §1 Euler’s Method 定义 　　　若某算法的局部截断误差为O(hp+1)，则称该算法有p 阶精度。 ? 欧拉法具有 1 阶精度。 在xn点用一阶向前差商近似一阶导数 在讨论差分公式时 介绍了差商的概念和性质， 各阶差商可以近似各阶导数，具有不同的精度， 且可以用函数值来表示。 Euler’s method §1 Euler’s Method §1 Euler’s Method ? 欧拉公式的改进： 隐式欧拉法或后退欧拉法 /* implicit Euler method or backward Euler method*/ xn+1点向后差商近似导数 隐式或后退欧拉公式 由于未知数 yn+1 同时出现在等式的两边，故称为隐式 /* implicit */ 欧拉公式，而前者称为显式 /* explicit */ 欧拉公式。隐式公式不能直接求解，一般需要用Euler显式公式得到初值，然后用Euler隐式公式迭代求解。因此隐式公式较显式公式计算复杂，但稳定性好（后面分析）。 收敛性 §1 Eu