化工问题的建模与数学分析方法化工数学.pptxVIP

化工问题的建模 与数学分析方法——ModellingandAnalyticalMethodsforProblemsinChemicalEngineering

第六章近似解析方法奇异摄动法试验函数法正交配置法

01解析解——由简单函数关系式直接给出的对应关系02结构简单，计算代价小03结果可靠，直观，便于应用04对一般问题难以得到05数值解——以大量数字对应方式给出的函数关系06适用性广，可处理复杂问题和大规模问题07依赖于计算工具和特定算法，代价较大解析解与数值解的比较第六章近似解析方法——概论

近似解析解——准确解的近似解析表达式01局部精确性较差，但整体规律性好02形式简单而满足工程应用03容易得到04数学问题的求解原则05首先求准确解析解06其次求近似解析解07最后采用数值解08第六章近似解析方法——概论

摄动法——将问题对小参数进行级数展开的求解方法1正则摄动：小参数直接展开的方法2奇异摄动：直接展开失效后采用的专门方法或改进方 法31摄动法第六章近似解析方法——摄动法

第六章近似解析方法——摄动法正则摄动与奇异摄动例1最高次项含小参数的非线性代数方程的求解设得

第六章近似解析方法——摄动法正则摄动只能得到一个根，因为直接展开失去了问题的非线性性质。

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

第六章近似解析方法——摄动法准确解为当?→0时，其两个根分别趋于y→a和y→?－1，对应的两个摄动解分别称为正则摄动解与奇异摄动解。

第六章近似解析方法——摄动法例2小参数位于非导数项中的情况设得

第六章近似解析方法——摄动法近似解与准确解极为接近，这种情况下正则摄动法是奏效的。

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

第六章近似解析方法——摄动法直接展开得到取x=1处的边界条件y0(1)=?，y1(1)=0，得到

在x=0处因此，近似解不满足x=0处的边值。第六章近似解析方法——摄动法

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

概念：渐近级数与收敛级数 收敛级数：按变量展开的级数，如泰勒级数，三角级数，幂级数等，级数的精度随项数的增加而提高；渐近级数：按参数展开的级数系数yn(x)是由展开后的问题顺序解出的，因此级数不一定收敛，一般只取级数的2～3项。第六章近似解析方法——摄动法

第六章近似解析方法——摄动法边界层方法基本思想：放大镜——将空间边界层放大，使分布变平缓，突出边界层内的作用；慢镜头——将时间尺度放大，使变化减缓，突出快速变化的过程。历史来源与发展：Prandtl边界层方程，Blasuis匹配方法，PLK方法

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

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

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

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

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

第六章近似解析方法——摄动法匹配Prandtl匹配原理——0阶近似的匹配方法得0阶内解

加法合成法0101020304合成解＝外解＋内解－公共部分高阶近似的匹配——VanDyke匹配原理n项外解的m项内部展开＝m项内解的n项外部展开0203044、合成第六章近似解析方法——摄动法

第六章近似解析方法——摄动法匹配后的两项近似内解合成后的两项近似解

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

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

第六章近似解析方法——摄动法快时间尺度解——定常态近似

第六章近似解析方法——摄动法合成与匹配——VonDyke匹配原理例：催化剂的平行失活问题反应快、失活慢，二者均需要考虑

第六章近似解析方法——摄动法无量纲化先求内解，内解可完全确定

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

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

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

第六章近似