浙大化工数学.ppt

谢 谢 第六章 近似解析方法——试验函数法 3、边界层动量积分方法 问题：Prandtl边界层方程，非线性PDE方程组 y=0: u=v=0; y→∞: u=U, v=0 x0: u=U, v=0 第六章 近似解析方法——试验函数法 方法要点： 在边界层内用积分形式的动量方程代替微分方程 选择满足边界条件的多项式或其它函数为试验函数 第六章 近似解析方法——试验函数法 1）边界层积分动量方程的推导 边界层厚度?(x)是一个待定的函数 第六章 近似解析方法——试验函数法 2）试验函数的选取 满足以下边界条件 取 第六章 近似解析方法——摄动法 无量纲化 1、先求内解，内解可完全确定 第六章 近似解析方法——摄动法 令 第六章 近似解析方法——摄动法 得到两项近似内解 第六章 近似解析方法——摄动法 2、直接展开求外解，外解不满足初值，含任意常数 3、内、外解匹配确定外解任意常数 得到外解 第六章 近似解析方法——摄动法 第六章 近似解析方法——摄动法 4、合成含有快、慢尺度的统一解 第六章 近似解析方法——摄动法 第六章 近似解析方法——摄动法 4、移动的空间边界层问题 非线性色谱过程的浓度前沿 非线性吸附效应与扩散效应之间的竞争作用 移动的空间边界层的形成 求解思路 外解——非线性色谱问题的激波解 内解——采用跟随激波的移动坐标系，放大边界层 匹配与合成 第六章 近似解析方法——摄动法 问题 第六章 近似解析方法——摄动法 1、外解 由特征线法 浓度激波位置xs由匹配条件确定 第六章 近似解析方法——摄动法 2、边界层内解 积分得 3、匹配 第六章 近似解析方法——摄动法 由以上Prandtl匹配条件得激波间断关系 解得激波轨迹 边界层内解 第六章 近似解析方法——摄动法 第六章 近似解析方法——摄动法 4、0阶近似合成解 第六章 近似解析方法——试验函数法 §2 试验函数方法 Trial Function Method 思想：用已知的、含待定参数的简单函数近似代替准确解，用积分形式的方程或点近似方程代替微分方程，确定不定参数。 以牺牲一些局部的精确性为代价，换取对问题整体规律性的把握，在一定的近似范围内解决问题。 要点：试验函数的选择;残差处理方法 Keypoint: selection of the trial functions treatment of the residual 第六章 近似解析方法——试验函数法 1、试验函数与方程残差 The residual of equations 例1 落石问题 分析：下落速度从零增加到末速度 第六章 近似解析方法——试验函数法 设试验函数为 ? 是待定参数，代入方程得到残差 若要求在t=τ时刻方程成立，R(τ)=0，得 第六章 近似解析方法——试验函数法 由准确解 特点：方程只在一个点满足，近似解“八九不离十” 例2 催化剂颗粒有效系数计算 第六章 近似解析方法——试验函数法 设试验函数 要求方程积分满足，得 第六章 近似解析方法——试验函数法 取s=0，r(y)=y ，得 准确解 ? 1 时，相差甚微（1%左右） ，? 越大相差越大。 原因：快速反应浓度分布空心化，偏离抛物分布。 第六章 近似解析方法——试验函数法 改进，对于快速反应，采用以下蛋白型试验函数 仍要求方程积分满足，确定参数xp 第六章 近似解析方法——试验函数法 准确解 ，说明试验函数越接近真实，结果越准确。 例3 试井问题 拭井：反求地层参数的工业试验方法，压力变化方程 第六章 近似解析方法——试验函数法 第六章 近似解析方法——试验函数法 分析：影响半径R＝R(?) ，漏斗型分布，拟稳态假设 无穷远边值的有限化 积分平均近似 拟稳态试验函数 第六章 近似解析方法——试验函数法 由边界条件 影响半径为待定函数，代入积分的压力方程，得 准确解 第六章 近似解析方法——试验函数法 小结：试验函数法 试验函数的选择 尽可能接近真实 事先满足初始与边界条件 方程残差的处理 点近似 积分平均近似