资源描述
单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,化工系统工程,第三章 数学模型求解方法,授课教师：杜可杰,化学化工学院,化工系统工程第三章 数学模型求解方法,3.4,非线性代数方程组迭代解法,*,（,不要求，成熟软件,）,3.4.1,直接迭代法（,direct substitution,direct iteration,）,3.4.2,部分迭代法,(partial substitution),3.4.3,牛顿（,Newton,）算法,F(x)=0,道理相似，理论相通，计算复杂！,3.4 非线性代数方程组迭代解法*（不要求，成熟软件）3,牛顿法也可写成,改进：,Broyden,法；修正的牛顿法等等,p77-78,Quasi-Newton Methods,避免了牛顿,每轮迭代需要求导数,和求逆的复杂运算,牛顿法也可写成改进：Broyden法；修正的牛顿法等等,3.5 迭代加速技术,*,(p78,不要求,),3.5.1 DEM,（,Dominant Eigenvalue Method,）,加速法,3.5.2,王德人加速法,(,略,p79),思想与策略,:,步长和方向的选择,-,考虑相邻两次迭代的误差向量的任意一种范数比,;,只利用变量及相邻两次误差计算范数信息,.,OR,适用范围,:,单调收敛,的过程,其他：多维,wegstein,法本质上,F(x)=0,代表了迭代收敛与改进的定量信息！,3.5 迭代加速技术*(p78,不要求)3.5.1 D,3.6 方程组求解的切割技术,简介,，,背景,化工流程模拟是一项非常复杂的任务，往往涉及到,成千上万个非线性方程,，而且这些方程涉及的,函数定义域狭窄、间断点多、连续性和可微性很差，各变量之间交互作用很强,，因此如简单地套用现成的算法往往导致求解效率低下，甚至根本无法求解。,在化工流程模拟实践中，总是需要先对数学模型进行一定的处理，将复杂问题分解为一系列较简单的问题，然后再根据具体情况，调用适当的算法进行求解。切割技术就是其中最重要的一种。,(,基于化工系统工程那部分内容,?),3.6 方程组求解的切割技术简介，背景,案例,1,，,f,1,f,3,f,2,f,4,x,2,x,4,x,1,x,3,X,3,*,信息流图,案例1，f1f3f2f4x2x4x1x3X3*信息流图,案例,2,，,f,1,，,f,2,，,f,3,为关于,x,2,，,x,3,，,x,4,的线性方程组，迭代策略？,f,1,f,2,f,3,f,4,X,1,*,信息流图,案例2，f1，f2，f3为关于x2，x3，x4的线性,Remarks,：,关于方程组切割技术,(p80,，,p82),，,1,.,选取部分变量作为,“,切割变量,”,则整个迭代过程从大局看变为仅对切割变量的迭代过程，使待求解问题的维数降低。,2,.,运用切割法可将难解的多个变量联立求解的问题分解转化为一系列易解的较少变量联立求解的问题。,3,.,切割法的作用是将原迭代问题,转化,为另一类迭代问题，其,本身并非迭代求根算法,。因此，使用了切割法后，仍需使用方程(组)迭代求根算法。,(,简化,),4,.,若将切割变量计算值直接用于下一轮对切割变量的迭代，则实际上对切割变量是在使用直接迭代法。,Remarks：关于方程组切割技术(p80，p82)，1,，,5.,使用切割技术，将,n,个未知数的联立求解问题转化为对少数,m,个切割变量的求解问题，将复杂的问题分解为一系列简单的问题。,关键,在于：适当或巧妙地选取切割变量与中间变量，并将原方程组适当地分解为中间变量方程组和切割变量方程组；所选切割变量（组）最宜是那些物理意义明确、易给初值并密切相关的一组变量。,（关于切割技术将在第四章结合化工过程模拟系统（序贯模块法）进行介绍。,，5.使用切割技术，将n个未知数的联立求解问题转化为对少,方程组切割技术,，,对切割技术运用得当，可从以下许多方面改善计算效率。,容易选取初始值，扩大收敛域。往往选取了切割变量的初值后，中间变量的初值就变得非常容易选取，甚或可直接计算出来。,提高数值稳定性。如中间变量方程组的形式适当，比如是线性方程组，就极为利于求解，不易发生误差传递和放大的现象。,减少内存占用。由于切割后回避了大规模方程组的联立求解，故通常许多大型矩阵和向量的工作单元都省去了。,最易被忽视而最重要的是，提高求解过程的适定性。化工过程数学模型的定义域往往非常狭窄、苛刻，各变量之间取值关系受到物理规律的严格约束。所有变量同时迭代极其容易使得迭代中间结果远超出定义域，违背物理规律，因而导致迭代混乱。正确选取切割变量后就能保证后续的计算在数学模型的定义域内进行，迭代中间结果也能保持物理意义，避免迭代失败。,方程组切割技术，对切割技术运用得当，可从以下许多方面改善计,方程组切割技术,，,提高了算法的收敛性。经过切割后，复杂问题被分解为若干基本简单问题，有利于选取收敛性好的迭代算法，加快计算速度。并且，使用收敛性高的算法，比较容易设置安全而合理的收敛判据和容差。,使用切割技术后，整体上或从外层看是仅对切割变量进行迭代，减少了问题的维数，但在求解中间变量方程组时也常常需要迭代。故此时形成外层迭代中嵌套着内层迭代的情况，似乎比较复杂。其实不然。切割法一般都是针对具有明确物理意义的大型复杂问题，巧妙分解原方程组，使得内层迭代适定性好、求解极其顺利，因而效率很高。如此时采用经典的联立求解，是难于解决上面论述的稳定性、收敛域、适定性、收敛速度等几方面问题的。,切割技术是适用于大规模复杂问题的求解技术，绝不是用来求解一般问题甚或习题的。许多方法都有类似的特点，适用于复杂问题的解法与适用于简单问题的解法是截然不同的。假如对极其简单的问题也采用切割技术，反而是生搬硬套，事倍功半。,方程组切割技术，提高了算法的收敛性。经过切割后，复杂问题,方程组切割技术,，,巧妙运用切割技术的关键在于适当选取切割变量（组）和对应的切割变量方程（组）。当方程组不太复杂时，变量数和方程数都较少，往往容易直接观察出较佳的切割方案。但对于维数稍多的问题，由于,“,组合爆炸,”,的缘故，基本上无法直接观察出合理的切割变量及其对应的方程。即使不考虑问题的物理意义，仅从纯数学的角度，实际上是能够给出选择切割变量及其对应方程的算法的，只不过对此类问题，数学上的合理解往往不是唯一的，也难于从数学角度判定何种方案是最佳的。,方程组分解就是将大型复杂方程组转化为若干简单方程组求解的方法。方程组分解的内容大致包括,“,分块,”,和,“,切割,”,两部分，目的就是解决切割变量选择及对应的方程选择问题。通俗地讲就是要解决由哪些子方程组去求解哪些未知数以及求解的可行次序。,方程组切割技术，巧妙运用切割技术的关键在于适当选取,3.7 双层法,(,不要求，闪蒸，精馏,）,严格模型,简化模型,3.7 双层法(不要求，闪蒸，精馏）严格模型 简化模型,3.8,常微分方程组初值问题与动态模拟（简介）,引入,(,概念,),欧拉法,龙格库塔法,3.8 常微分方程组初值问题与动态模拟（简介）引入(概念,案例,1:,连续流动加热水槽动态模拟,符号说明：,F,T,流量,(kg/hr),与温度,(K),，,T,0,环境温度,(K),，,C,物性，比热,(kJ/kg.K),与密度,(kg/m,3,),，,K,A,传热系数,(kJ/h K m,3,),与传热面积,(,底面积,)(m,2,),，,L,V,液位,(m),与持液体积,(m,3,),，,t,时间,(h),建模简化,假设：,K,，,C,，,为常数，理想搅拌（集中参数体系），流体不可压缩，忽略设备热容与传递滞后，焓衡算用显热衡算代替。,-,液位与温度随时间的动态变化关系,3.8.1,引入,(,定义,),（p92）,案例1:连续流动加热水槽动态模拟 符号说明：F,T 流量,连续流动加热水槽动态模型,首先考虑,质量衡算,：，即,(1),热量衡算,(,显热衡算,),：,即,:,(2,),利用,(2-12),，代入上式，消去,整理得：,(3),式,(1),与式,(3,即构成了所需的数学模型（机理模型）。此,模型描述了液位与温度随时间的动态变化关系,。如给定适当初始条件，即可算出液位与温度随时间的变化曲线,(,初值问题,),。,列出数学表达式,数学推演,连续流动加热水槽动态模型 首先考虑质量衡算：,案例,2:,物料冷却的数学模型,它是含有自变量,t,时间,),、位置函数,T,（随时间变化的物料温度）、,To,（环境温度）、,k(,降温速率）以及温度的一阶导数 ，是一个常微分方程。,案例2:物料冷却的数学模型 它是含有自变量t时间)、位置,3.8.1,定义（p92）,一阶常微分方程：,（,aX b,）,(X,为纯量,;y,为纯量或向量,),数值解法就是寻求方程,(1),、（,2,）的近似解,即求区间（,a,，,b,）上各个分点序列,X,n,，,n=1,，,2,，,m,的数值解,y,n,化工计算中，只有一些特殊形式的,F(x,y),才能计算它们的解析解。大多数情况下只能计算数值解。如果,x,为时间，则表示动态过程。其求解过程即动态模拟,。,数学上成为,”,初值,”,问题,3.8.1定义（p92）一阶常微分方程：化工计算中，只有一些,3.8.2,欧拉法,（,Euler,）,欧拉,:18,世纪科学界的代表人物，是那个时代的巨人。他是历来最有才华、最博学的人物之一，也是历史上最多产的一位数学家,x,0,x,1,x,2,y,(,x,),y,2,y,1,y,0,斜率,向前欧拉法,改进欧拉法（预测校正法精度较高）,y,x,基本思想,：在求解区间上作等距分割，步长为,h,，用差商代替导数计算常微分方程。,作,f(,x,),在,x,i,处的,1,阶向前差商：,d,y,/d,x,y,(,x,(i+1),）,-,y,(,x,i,)/h,d,y/,d,x=f(x,y),f(x,i,y,i,）,向后欧拉法,y,(i+1)=,y,(i),+h,f,（,x,(i+1,),y,(i+,1,),),（,p96,）,预报值,y=f(x,y),3.8.2欧拉法（Euler）欧拉:18世纪科学界的代表人物,向前欧拉公式,应用实例,例,1,：假定某物体的温度,w,因自热而产生的热量可以使物体在每秒钟内以,4%,的速度增长，同时该物体由于散热可使其温度在每秒种内下降,100,k,，则物体温度随时间变化的微分方程：,（,t,以秒为单位,),分别以初始温,x,(0)=1500k,，,y,(0)=2500k,，,z,(0)=3500k,用欧拉公式预测,24,秒后的物体温度趋势。,向前欧拉公式 应用实例例1：假定某物体的温度w因自热而产生的,向前欧拉公式,应用实例,解：,w,0,分别以,x,0,=1500,，,y,0,=2500,，,z,0,=3500,代入。计算结果见如表所示。,图,4-3,三种初始值的温度变化曲线,表,4-1,Note,：表与图的使用,向前欧拉公式 应用实例解：图4-3 三种初始值的温度变化曲线,4.2.1,向前欧拉公式,实例,当自热引起物体温度升高的速度小于散热引起温度下降的速度，物体的温度随时间而逐渐降低；当自热引起物体温度升高的速度与散热引起温度下降的速度平衡时，物体的温度保持不变；当自热引起物体温度升高的速度大于散热引起温度下降的速度，物体的温度随时间而升高。在图中,L,1,，,L,2,，,L,3,分别表示初始值,3500k,，,2500k,和,1500k,的三条温度变化趋势曲线。,4.2.1 向前欧拉公式 实例当自热,Remarks,方 法,显式欧拉,隐式欧拉,改进欧拉法,简单,精度低,稳定性最好,精度低,计算量大,精度提高,计算量大,Remar
展开阅读全文