化工问题的建模与数学分析方法化工数学6-课件

化工问题的建模 与数学分析方法——ModellingandAnalyticalMethodsforProblemsinChemicalEngineering化工问题的建模1第六章近似解析方法1、奇异摄动法2、试验函数法3、正交配置法第六章近似解析方法2第六章近似解析方法——概论解析解与数值解的比较

解析解——由简单函数关系式直接给出的对应关系 结构简单，计算代价小 结果可靠，直观，便于应用 对一般问题难以得到 数值解——以大量数字对应方式给出的函数关系 适用性广，可处理复杂问题和大规模问题依赖于计算工具和特定算法，代价较大第六章近似解析方法——概论解析解与数值解的比较3第六章近似解析方法——概论 近似解析解——准确解的近似解析表达式 局部精确性较差，但整体规律性好 形式简单而满足工程应用 容易得到数学问题的求解原则 首先求准确解析解 其次求近似解析解 最后采用数值解第六章近似解析方法——概论 近似解析解——准确解的近似解析表4第六章近似解析方法——摄动法§1摄动法

摄动法——将问题对小参数进行级数展开的求解方法 正则摄动：小参数直接展开的方法 奇异摄动：直接展开失效后采用的专门方法或改进方 法第六章近似解析方法——摄动法§1摄动法5第六章近似解析方法——摄动法1、正则摄动与奇异摄动例1最高次项含小参数的非线性代数方程的求解 设

得

第六章近似解析方法——摄动法1、正则摄动与奇异摄动6第六章近似解析方法——摄动法正则摄动只能得到一个根，因为直接展开失去了问题的非线性性质。 第六章近似解析方法——摄动法7第六章近似解析方法——摄动法如果作变换y=u/，得

然后对u直接展开，得到另一个根第六章近似解析方法——摄动法如果作变换y=u/8第六章近似解析方法——摄动法

准确解为

当→0时，其两个根分别趋于y→a和y→－1，对应的两个摄动解分别称为正则摄动解与奇异摄动解。第六章近似解析方法——摄动法 9第六章近似解析方法——摄动法例2小参数位于非导数项中的情况 设 得第六章近似解析方法——摄动法例2小参数位于非导数项中的情10第六章近似解析方法——摄动法近似解与准确解极为接近，这种情况下正则摄动法是奏效的。

第六章近似解析方法——摄动法11第六章近似解析方法——摄动法

例3方程最高阶导数乘小参数的情况 当＝0时，方程由二阶退化成一阶方程，近似解只能满足一个边值而难以同时满足两个边值。第六章近似解析方法——摄动法 例3方程最高阶导数乘小参数12第六章近似解析方法——摄动法直接展开得到

取x=1处的边界条件y0(1)=，y1(1)=0，得到

第六章近似解析方法——摄动法直接展开得到13第六章近似解析方法——摄动法

在x=0处 因此，近似解不满足x=0处的边值。第六章近似解析方法——摄动法 14第六章近似解析方法——摄动法分析：x=0处存在一个边界层 边界层的存在是小参数乘最高阶导数问题的特征

第六章近似解析方法——摄动法分析：x=0处存在一个边15第六章近似解析方法——摄动法概念：渐近级数与收敛级数 收敛级数：按变量展开的级数，如泰勒级数，三角级数，幂级数等，级数的精度随项数的增加而提高； 渐近级数：按参数展开的级数 系数yn(x)是由展开后的问题顺序解出的，因此级数不一定收敛，一般只取级数的2～3项。第六章近似解析方法——摄动法概念：渐近级数与收敛级数 16第六章近似解析方法——摄动法2、边界层方法 基本思想：放大镜——将空间边界层放大，使分布变平缓，突出边界层内的作用；慢镜头——将时间尺度放大，使变化减缓，突出快速变化的过程。历史来源与发展： Prandtl边界层方程，Blasuis匹配方法，PLK方法第六章近似解析方法——摄动法2、边界层方法17第六章近似解析方法——摄动法边界层方法的求解步骤 1、外解——直接展开

2、内解——边界层放大 3、匹配——内解与外解的衔接 4、合成——内解与外解的组合第六章近似解析方法——摄动法边界层方法的求解步骤18第六章近似解析方法——摄动法例3

1、外解

第六章近似解析方法——摄动法例319第六章近似解析方法——摄动法

2、内解

边界层放大，定义内部坐标第六章近似解析方法——摄动法 20第六章近似解析方法——摄动法

取＝1以保留二阶导数项，得 令 得第六章近似解析方法——摄动法 取＝1以保留二阶导数项，得21第六章近似解析方法——摄动法

解出0阶近似 常数C由匹配条件确定第六章近似解析方法——摄动法 22第六章近似解析方法——摄动法 3、匹配

Prandtl匹配原理——0阶近似的匹配方法 得0阶内解第六章近似解析方法——摄动法 3、匹配23第六章近似解析方法——摄动法 4、合成

加法合成法 合成解＝外解＋内解－公共部分

高阶近似的匹配——VanDyke匹配原理 n项外解的m项内部展开＝m项内解的n项外部展开

第六章近似解析方法——摄动法 4、合成24第六章近似解析方法——摄动法匹配后的两项近似内解合成后的两项近似解第六章近似解析方法——摄动法匹配后的两项近似内解25第六章近似解析方法——摄动法3、时间边界层——刚性问题（stiffequs)

刚性问题：具有不同时间尺度的变化问题； 特点：快步骤与慢步骤共存 拟稳态近似与定常态近似 计算难点：数值振荡，多步Gear方法 奇异摄动：慢镜头分析，给出完整的结果第六章近似解析方法——摄动法3、时间边界层——刚性问题（s26第六章近似解析方法——摄动法 例

慢时间尺度解（＝0）——拟稳态近似第六章近似解析方法——摄动法 例27第六章近似解析方法——摄动法

快时间尺度解——定常态近似第六章近似解析方法——摄动法28第六章近似解析方法——摄动法合成与匹配——VonDyke匹配原理 例：催化剂的平行失活问题

反应快、失活慢，二者均需要考虑第六章近似解析方法——摄动法合成与匹配——VonDyke29第六章近似解析方法——摄动法无量纲化

1、先求内解，内解可完全确定

第六章近似解析方法——摄动法无量纲化 30第六章近似解析方法——摄动法

令

第六章近似解析方法——摄动法 31第六章近似解析方法——摄动法得到两项近似内解

第六章近似解析方法——摄动法得到两项近似内解 32第六章近似解析方法——摄动法2、直接展开求外解，外解不满足初值，含任意常数3、内、外解匹配确定外解任意常数

得到外解第六章近似解析方法——摄动法2、直接展开求外解，外解不满足33第六章近似解析方法——摄动法

第六章近似解析方法——摄动法 34第六章近似解析方法——摄动法4、合成含有快、慢尺度的统一解

第六章近似解析方法——摄动法4、合成含有快、慢尺度的统一解35第六章近似解析方法——摄动法

第六章近似解析方法——摄动法36第六章近似解析方法——摄动法4、移动的空间边界层问题

非线性色谱过程的浓度前沿 非线性吸附效应与扩散效应之间的竞争作用 移动的空间边界层的形成求解思路 外解——非线性色谱问题的激波解 内解——采用跟随激波的移动坐标系，放大边界层 匹配与合成第六章近似解析方法——摄动法4、移动的空间边界层问题37第六章近似解析方法——摄动法问题

第六章近似解析方法——摄动法问题 38第六章近似解析方法——摄动法1、外解

由特征线法浓度激波位置xs由匹配条件确定

第六章近似解析方法——摄动法1、外解39第六章近似解析方法——摄动法2、边界层内解

积分得3、匹配第六章近似解析方法——摄动法2、边界层内解40第六章近似解析方法——摄动法由以上Prandtl匹配条件得激波间断关系解得激波轨迹

边界层内解第六章近似解析方法——摄动法由以上Prandtl匹配条件得41第六章近似解析方法——摄动法

第六章近似解析方法——摄动法42第六章近似解析方法——摄动法4、0阶近似合成解

第六章近似解析方法——摄动法4、0阶近似合成解43第六章近似解析方法——试验函数法§2试验函数方法

思想：用已知的、含待定参数的简单函数近似代替准确解，用积分形式的方程或点近似方程代替微分方程，确定不定参数。 以牺牲一些局部的精确性为代价，换取对问题整体规律性的把握，在一定的近似范围内解决问题。 要点：试验函数的选择 残差处理方法第六章近似解析方法——试验函数法§2试验函数方法44第六章近似解析方法——试验函数法1、试验函数与方程残差

例1落石问题 分析：下落速度从零增加到末速度第六章近似解析方法——试验函数法1、试验函数与方程残差45第六章近似解析方法——试验函数法设试验函数为

是待定参数，代入方程得到残差若要求在t=τ时刻方程成立，R(τ)=0，得第六章近似解析方法——试验函数法设试验函数为46第六章近似解析方法——试验函数法由准确解

特点：方程只在一个点满足，近似解“八九不离十” 例2催化剂颗粒有效系数计算

第六章近似解析方法——试验函数法由准确解47第六章近似解析方法——试验函数法

设试验函数要求方程积分满足，得第六章近似解析方法——试验函数法48第六章近似解析方法——试验函数法取s=0，r(y)=y，得

准确解<1时，相差甚微（1%左右），越大相差越大。 原因：快速反应浓度分布空心化，偏离抛物分布。第六章近似解析方法——试验函数法取s=0，r(y)=y，49第六章近似解析方法——试验函数法改进，对于快速反应，采用以下蛋白型试验函数

仍要求方程积分满足，确定参数xp第六章近似解析方法——试验函数法改进，对于快速反应，采用以50第六章近似解析方法——试验函数法

准确解，说明试验函数越接近真实，结果越准确。

例3试井问题

拭井：反求地层参数的工业试验方法，压力变化方程第六章近似解析方法——试验函数法51第六章近似解析方法——试验函数法

第六章近似解析方法——试验函数法 52第六章近似解析方法——试验函数法分析：影响半径R＝R()，漏斗型分布，拟稳态假设

无穷远边值的有限化

积分平均近似

拟稳态试验函数第六章近似解析方法——试验函数法分析：影响半径R＝R()53第六章近似解析方法——试验函数法由边界条件

影响半径为待定函数，代入积分的压力方程，得准确解第六章近似解析方法——试验函数法由边界条件54第六章近似解析方法——试验函数法小结：试验函数法试验函数的选择 尽可能接近真实 事先满足初始与边界条件方程残差的处理 点近似 积分平均近似 加权积分近似第六章近似解析方法——试验函数法小结：试验函数法55第六章近似解析方法——试验函数法2、空间平均近似

例：球形颗粒上的不定常扩散

采用抛物型试验函数:第六章近似解析方法——试验函数法2、空间平均近似56第六章近似解析方法——试验函数法代入方程，令空间积分为0，得

系数A由初始条件确定，定义空间平均浓度，得由初值为0

第六章近似解析方法——试验函数法代入方程，令空间积分为0，57第六章近似解析方法——试验函数法近似解与准确解的比较： 长时间后准确，短时间内偏离。 原因：渗透区的存在，偏离抛物型试验函数。第六章近似解析方法——试验函数法近似解与准确解的比较：58第六章近似解析方法——试验函数法第六章近似解析方法——试验函数法59第六章近似解析方法——试验函数法改进——取渗透型试验函数由空间平均近似第六章近似解析方法——试验函数法改进——取渗透型试验函数60第六章近似解析方法——试验函数法短时间解准确解第六章近似解析方法——试验函数法61第六章近似解析方法——试验函数法3、边界层动量积分方法 问题：Prandtl边界层方程，非线性PDE方程组

y=0:u=v=0;y→∞:u=U,v=0

x<0:u=U,v=0第六章近似解析方法——试验函数法3、边界层动量积分方法62第六章近似解析方法——试验函数法

方法要点： 在边界层内用积分形式的动量方程代替微分方程 选择满足边界条件的多项式或其它函数为试验函数第六章近似解析方法——试验函数法63第六章近似解析方法——试验函数法1）边界层积分动量方程的推导 边界层厚度(x)是一个待定的函数第六章近似解析方法——试验函数法1）边界层积分动量方程的推64第六章近似解析方法——试验函数法2）试验函数的选取 满足以下边界条件取第六章近似解析方法——试验函数法2）试验函数的选取65第六章近似解析方法——试验函数法代入动量积分方程得到确定边界层厚度(x)的方程准确解第六章近似解析方法——试验函数法66第六章近似解析方法——试验函数法4、加权余量法（Galerkin方法）

残差加权积分为0的近似方法权函数的选择 点近似 积分平均近似 矩量积分近似 最小二乘近似第六章近似解析方法——试验函数法4、加权余量法（Galer67第六章近似解析方法——试验函数法Galerkin方法

Galerkin方法应用最广，其物理基础为变分原理第六章近似解析方法——试验函数法Galerkin方法68第六章近似解析方法——试验函数法例：催化剂有效系数计算取试验函数方程残差第六章近似解析方法——试验函数法例：催化剂有效系数计算69第六章近似解析方法——试验函数法由Galerkin方法n=2时第六章近似解析方法——试验函数法由Galerkin方法70第六章近似解析方法——试验函数法 准确解＝0.446比较：空间平均近似法 =0.384

第六章近似解析方法——试验函数法71第六章近似解析方法——正交配置法§3正交配置法

Galerkin法的特点精度高积分计算量大改进用高斯积分公式进行Galerkin法的积分运算——正交配置法（Finlayson1972,Villadsen1978)第六章近似解析方法——正交配置法§3正交配置法72第六章近似解析方法——正交配置法1、以待定参数为未知量的正交配置法

例：催化剂颗粒问题 设

残差 Galerkin法确定参数a的方程第六章近似解析方法——正交配置法1、以待定参数为未知量的正73第六章近似解析方法——正交配置法采用Gauss求积公式来计算Galerkin积分取N个Jacobi正交多项式的根为节点，得到

第六章近似解析方法——正交配置法74第六章近似解析方法——正交配置法上式是关于N个残差项RNj的线性代数方程组，由于

得到充分必要条件上式为N个节点上的配置方程，可定出N个参数a。因此，当采用高斯积分公式来计算Galerkin积分时，Galerkin方法就成为高斯节点上的配置方法，二者的误差就是高斯积分误差，具有2N－1阶代数精度。相应的方法称为正交配置法。第六章近似解析方法——正交配置法上式是关于N个残差项RNj75第六章近似解析方法——正交配置法Jaccobi正交多项式的概念 (0,1)区间上的特征值问题的非零解

N次正交多项式在(0,1)区间的N个根就是配置节点。例：催化剂有效系数的计算

第六章近似解析方法——正交配置法Jaccobi正交多项式的76第六章近似解析方法——正交配置法 用高斯积分公式来代替Galerkin积分公式计算残差，单点近似 z1=1/3，仍取n=2,=3，解出

第六章近似解析方法——正交配置法 用高斯积分公式来代替Ga77第六章近似解析方法——正交配置法正交配置解Galerkin方法

＝0.463

差别源于高斯积分的微小误差。小结：

Galerkin法：通过残差加权积分为零确定参数 正交配置法：残差在正交节点上直接为零确定参数 二则差别仅为高斯积分近似误差第六章近似解析方法——正交配置法正交配置解78第六章近似解析方法——正交配置法2、以节点函数值为未知量的正交配置法取正交节点上的函数值yi=yi(x)为待定参数，代替ai，试验函数取为N次多项式Lagrangian插值函数第六章近似解析方法——正交配置法2、以节点函数值为未知量的79第六章近似解析方法——正交配置法以节点函数值为待定参数，意义明确，插值函数已知，便于求导导数的离散第六章近似解析方法——正交配置法以节点函数值为待定参数，意80第六章近似解析方法——正交配置法因此，对于任意的微分方程L(y)=0，采用正交配置法以后，就可以将其离散为配置节点xi上的N个代数方程

L(yi)=0,i=1,2,…,N 求出N个yi后，得解。例：圆柱形催化剂颗粒的内扩散问题，取

第六章近似解析方法——正交配置法因此，对于任意的微分方程L81第六章近似解析方法——正交配置法试验函数节点残差第六章近似解析方法——正交配置法试验函数82第六章近似解析方法——正交配置法代入边界条件yN+1=1，得到

确定N个未知量的N个线性代数方程，数值求解。

第六章近似解析方法——正交配置法代入边界条件yN+1=1，83第六章近似解析方法——正交配置法小结 正交配置法是一种解析解与数值解相结合的试验函数方法，它以正交多项式的根为插值节点，以节点函数值为参数，通过节点残差为零的条件将微分方程化为N个代数方程，确定N个未知函数值。 节点上的导数离散公式保证了正交配置法可应用于大规模的数值模拟问题，其计算量小，精度高，得到的结果是解析形式的多项式解。

第六章近似解析方法——正交配置法小结84第六章近似解析方法——正交配置法辅助程序 1、Jaccobi多项式的根和各阶导数计算

JCOBI(ND,N,N0,N1,ALFA,BETA,DIF1,DIF2, DIF3,ROOT） 2、插值多项式的微分矩阵和高斯求积系数第六章近似解析方法——正交配置法辅助程序85第六章近似解析方法——正交配置法DFOPR(ND,N,N0,N1,I,ID,DIF1,DIF2,DIF3,ROOT,VEC)第六章近似解析方法——正交配置法86谢谢谢谢87化工问题的建模 与数学分析方法——ModellingandAnalyticalMethodsforProblemsinChemicalEngineering化工问题的建模88第六章近似解析方法1、奇异摄动法2、试验函数法3、正交配置法第六章近似解析方法89第六章近似解析方法——概论解析解与数值解的比较

解析解——由简单函数关系式直接给出的对应关系 结构简单，计算代价小 结果可靠，直观，便于应用 对一般问题难以得到 数值解——以大量数字对应方式给出的函数关系 适用性广，可处理复杂问题和大规模问题依赖于计算工具和特定算法，代价较大第六章近似解析方法——概论解析解与数值解的比较90第六章近似解析方法——概论 近似解析解——准确解的近似解析表达式 局部精确性较差，但整体规律性好 形式简单而满足工程应用 容易得到数学问题的求解原则 首先求准确解析解 其次求近似解析解 最后采用数值解第六章近似解析方法——概论 近似解析解——准确解的近似解析表91第六章近似解析方法——摄动法§1摄动法

摄动法——将问题对小参数进行级数展开的求解方法 正则摄动：小参数直接展开的方法 奇异摄动：直接展开失效后采用的专门方法或改进方 法第六章近似解析方法——摄动法§1摄动法92第六章近似解析方法——摄动法1、正则摄动与奇异摄动例1最高次项含小参数的非线性代数方程的求解 设

得

第六章近似解析方法——摄动法1、正则摄动与奇异摄动93第六章近似解析方法——摄动法正则摄动只能得到一个根，因为直接展开失去了问题的非线性性质。 第六章近似解析方法——摄动法94第六章近似解析方法——摄动法如果作变换y=u/，得

然后对u直接展开，得到另一个根第六章近似解析方法——摄动法如果作变换y=u/95第六章近似解析方法——摄动法

准确解为

当→0时，其两个根分别趋于y→a和y→－1，对应的两个摄动解分别称为正则摄动解与奇异摄动解。第六章近似解析方法——摄动法 96第六章近似解析方法——摄动法例2小参数位于非导数项中的情况 设 得第六章近似解析方法——摄动法例2小参数位于非导数项中的情97第六章近似解析方法——摄动法近似解与准确解极为接近，这种情况下正则摄动法是奏效的。

第六章近似解析方法——摄动法98第六章近似解析方法——摄动法

例3方程最高阶导数乘小参数的情况 当＝0时，方程由二阶退化成一阶方程，近似解只能满足一个边值而难以同时满足两个边值。第六章近似解析方法——摄动法 例3方程最高阶导数乘小参数99第六章近似解析方法——摄动法直接展开得到

取x=1处的边界条件y0(1)=，y1(1)=0，得到

第六章近似解析方法——摄动法直接展开得到100第六章近似解析方法——摄动法

在x=0处 因此，近似解不满足x=0处的边值。第六章近似解析方法——摄动法 101第六章近似解析方法——摄动法分析：x=0处存在一个边界层 边界层的存在是小参数乘最高阶导数问题的特征

第六章近似解析方法——摄动法分析：x=0处存在一个边102第六章近似解析方法——摄动法概念：渐近级数与收敛级数 收敛级数：按变量展开的级数，如泰勒级数，三角级数，幂级数等，级数的精度随项数的增加而提高； 渐近级数：按参数展开的级数 系数yn(x)是由展开后的问题顺序解出的，因此级数不一定收敛，一般只取级数的2～3项。第六章近似解析方法——摄动法概念：渐近级数与收敛级数 103第六章近似解析方法——摄动法2、边界层方法 基本思想：放大镜——将空间边界层放大，使分布变平缓，突出边界层内的作用；慢镜头——将时间尺度放大，使变化减缓，突出快速变化的过程。历史来源与发展： Prandtl边界层方程，Blasuis匹配方法，PLK方法第六章近似解析方法——摄动法2、边界层方法104第六章近似解析方法——摄动法边界层方法的求解步骤 1、外解——直接展开

2、内解——边界层放大 3、匹配——内解与外解的衔接 4、合成——内解与外解的组合第六章近似解析方法——摄动法边界层方法的求解步骤105第六章近似解析方法——摄动法例3

1、外解

第六章近似解析方法——摄动法例3106第六章近似解析方法——摄动法

2、内解

边界层放大，定义内部坐标第六章近似解析方法——摄动法 107第六章近似解析方法——摄动法

取＝1以保留二阶导数项，得 令 得第六章近似解析方法——摄动法 取＝1以保留二阶导数项，得108第六章近似解析方法——摄动法

解出0阶近似 常数C由匹配条件确定第六章近似解析方法——摄动法 109第六章近似解析方法——摄动法 3、匹配

Prandtl匹配原理——0阶近似的匹配方法 得0阶内解第六章近似解析方法——摄动法 3、匹配110第六章近似解析方法——摄动法 4、合成

加法合成法 合成解＝外解＋内解－公共部分

高阶近似的匹配——VanDyke匹配原理 n项外解的m项内部展开＝m项内解的n项外部展开

第六章近似解析方法——摄动法 4、合成111第六章近似解析方法——摄动法匹配后的两项近似内解合成后的两项近似解第六章近似解析方法——摄动法匹配后的两项近似内解112第六章近似解析方法——摄动法3、时间边界层——刚性问题（stiffequs)

刚性问题：具有不同时间尺度的变化问题； 特点：快步骤与慢步骤共存 拟稳态近似与定常态近似 计算难点：数值振荡，多步Gear方法 奇异摄动：慢镜头分析，给出完整的结果第六章近似解析方法——摄动法3、时间边界层——刚性问题（s113第六章近似解析方法——摄动法 例

慢时间尺度解（＝0）——拟稳态近似第六章近似解析方法——摄动法 例114第六章近似解析方法——摄动法

快时间尺度解——定常态近似第六章近似解析方法——摄动法115第六章近似解析方法——摄动法合成与匹配——VonDyke匹配原理 例：催化剂的平行失活问题

反应快、失活慢，二者均需要考虑第六章近似解析方法——摄动法合成与匹配——VonDyke116第六章近似解析方法——摄动法无量纲化

1、先求内解，内解可完全确定

第六章近似解析方法——摄动法无量纲化 117第六章近似解析方法——摄动法

令

第六章近似解析方法——摄动法 118第六章近似解析方法——摄动法得到两项近似内解

第六章近似解析方法——摄动法得到两项近似内解 119第六章近似解析方法——摄动法2、直接展开求外解，外解不满足初值，含任意常数3、内、外解匹配确定外解任意常数

得到外解第六章近似解析方法——摄动法2、直接展开求外解，外解不满足120第六章近似解析方法——摄动法

第六章近似解析方法——摄动法 121第六章近似解析方法——摄动法4、合成含有快、慢尺度的统一解

第六章近似解析方法——摄动法4、合成含有快、慢尺度的统一解122第六章近似解析方法——摄动法

第六章近似解析方法——摄动法123第六章近似解析方法——摄动法4、移动的空间边界层问题

非线性色谱过程的浓度前沿 非线性吸附效应与扩散效应之间的竞争作用 移动的空间边界层的形成求解思路 外解——非线性色谱问题的激波解 内解——采用跟随激波的移动坐标系，放大边界层 匹配与合成第六章近似解析方法——摄动法4、移动的空间边界层问题124第六章近似解析方法——摄动法问题

第六章近似解析方法——摄动法问题 125第六章近似解析方法——摄动法1、外解

由特征线法浓度激波位置xs由匹配条件确定

第六章近似解析方法——摄动法1、外解126第六章近似解析方法——摄动法2、边界层内解

积分得3、匹配第六章近似解析方法——摄动法2、边界层内解127第六章近似解析方法——摄动法由以上Prandtl匹配条件得激波间断关系解得激波轨迹

边界层内解第六章近似解析方法——摄动法由以上Prandtl匹配条件得128第六章近似解析方法——摄动法

第六章近似解析方法——摄动法129第六章近似解析方法——摄动法4、0阶近似合成解

第六章近似解析方法——摄动法4、0阶近似合成解130第六章近似解析方法——试验函数法§2试验函数方法

思想：用已知的、含待定参数的简单函数近似代替准确解，用积分形式的方程或点近似方程代替微分方程，确定不定参数。 以牺牲一些局部的精确性为代价，换取对问题整体规律性的把握，在一定的近似范围内解决问题。 要点：试验函数的选择 残差处理方法第六章近似解析方法——试验函数法§2试验函数方法131第六章近似解析方法——试验函数法1、试验函数与方程残差

例1落石问题 分析：下落速度从零增加到末速度第六章近似解析方法——试验函数法1、试验函数与方程残差132第六章近似解析方法——试验函数法设试验函数为

是待定参数，代入方程得到残差若要求在t=τ时刻方程成立，R(τ)=0，得第六章近似解析方法——试验函数法设试验函数为133第六章近似解析方法——试验函数法由准确解

特点：方程只在一个点满足，近似解“八九不离十” 例2催化剂颗粒有效系数计算

第六章近似解析方法——试验函数法由准确解134第六章近似解析方法——试验函数法

设试验函数要求方程积分满足，得第六章近似解析方法——试验函数法135第六章近似解析方法——试验函数法取s=0，r(y)=y，得

准确解<1时，相差甚微（1%左右），越大相差越大。 原因：快速反应浓度分布空心化，偏离抛物分布。第六章近似解析方法——试验函数法取s=0，r(y)=y，136第六章近似解析方法——试验函数法改进，对于快速反应，采用以下蛋白型试验函数

仍要求方程积分满足，确定参数xp第六章近似解析方法——试验函数法改进，对于快速反应，采用以137第六章近似解析方法——试验函数法

准确解，说明试验函数越接近真实，结果越准确。

例3试井问题

拭井：反求地层参数的工业试验方法，压力变化方程第六章近似解析方法——试验函数法138第六章近似解析方法——试验函数法

第六章近似解析方法——试验函数法 139第六章近似解析方法——试验函数法分析：影响半径R＝R()，漏斗型分布，拟稳态假设

无穷远边值的有限化

积分平均近似

拟稳态试验函数第六章近似解析方法——试验函数法分析：影响半径R＝R()140第六章近似解析方法——试验函数法由边界条件

影响半径为待定函数，代入积分的压力方程，得准确解第六章近似解析方法——试验函数法由边界条件141第六章近似解析方法——试验函数法小结：试验函数法试验函数的选择 尽可能接近真实 事先满足初始与边界条件方程残差的处理 点近似 积分平均近似 加权积分近似第六章近似解析方法——试验函数法小结：试验函数法142第六章近似解析方法——试验函数法2、空间平均近似

例：球形颗粒上的不定常扩散

采用抛物型试验函数:第六章近似解析方法——试验函数法2、空间平均近似143第六章近似解析方法——试验函数法代入方程，令空间积分为0，得

系数A由初始条件确定，定义空间平均浓度，得由初值为0

第六章近似解析方法——试验函数法代入方程，令空间积分为0，144第六章近似解析方法——试验函数法近似解与准确解的比较： 长时间后准确，短时间内偏离。 原因：渗透区的存在，偏离抛物型试验函数。第六章近似解析方法——试验函数法近似解与准确解的比较：145第六章近似解析方法——试验函数法第六章近似解析方法——试验函数法146第六章近似解析方法——试验函数法改进——取渗透型试验函数由空间平均近似第六章近似解析方法——试验函数法改进——取渗透型试验函数147第六章近似解析方法——试验函数法短时间解准确解第六章近似解析方法——试验函数法148第六章近似解析方法——试验函数法3、边界层动量积分方法 问题：Prandtl边界层方程，非线性PDE方程组

y=0:u=v=0;y→∞:u=U,v=0

x<0:u=U,v=0第六章近似解析方法——试验函数法3、边界层动量积分方法149第六章近似解析方法——试验函数法

方法要点： 在边界层内用积分形式的动量方程代替微分方程 选择满足边界条件的多项式或其它函数为试验函数第六章近似解析方法——试验函数法150第六章近似解析方法——试验函数法1）边界层积分动量方程的推导 边界层厚度(x)是一个待定的函数第六章近似解析方法——试验函数法1）边界层积分动量方程的推151第六章近似解析方法——试验函数法2）试验函数的选取 满足以下边界条件取第六章近似解析方法——试验函数法2）试验函数的选取152第六章近似解析方法——试验函数法代入动量积分方程得到确定边界层厚度(x)的方程准确解第六章近似解析方法——试验函数法153第六章近似解析方法——试验函数法4、加权余量法（Galerkin方法）

残差加权积分为0的近似方法权函数的选择 点近似 积分平均近似 矩量积分近似 最小二乘近似第六章近似解析方法——试验函数法4、加权余量法（Galer154第六章近似解析方法——试验函数法Galerkin方法

Galerkin方法应用最广，其物理基础为变分原理第六章近似解析方法——试验函数法Galerkin方法155第六章近似解析方法——试验函数法例：催化剂有效系数计算取试验函数方程残差第六章近似解析方法——试验函数法例：催化剂有效系数计算156第六章近似解析方法——试验函数法由Galerkin方法n=2时第六章近似解析方法——试验函数法由Galerkin方法157第六章近似解析方法——试验函数法 准确解＝0.446比较：空间平均近似法 =0.384

第六章近似解析方法——试验函数法158第六章近似解析方法——正交配置法§3正交配置法

Galerkin法的特点精度高积分计算量大改进用高斯积分公式进行Galerkin法的积分运算——正交配置法（Finlayson1972,Villadsen1978)第六章近似解析方法——正交配置法§3正交配置法159第六章近似解析方法——正交配置法1、以待定参数为未知量的正交配置法

例：催化剂颗粒问题 设

残差 Galerkin法确定参数a的方程第六章近似解析方法——正交配置法1、以待定参数为未知量的正160第六章近似解析方法——正交配置法