偏微分方程课件

1、浙江大学实用数值计算方法17.1 偏微分方程的一般介绍 Partial Differential Equations(PDEs) 自变量数 至少2个 阶数 方程中导数的最高阶数 性态 以一阶方程为例浙江大学实用数值计算方法27.1浙江大学实用数值计算方法37.1 类型 一阶栓区型方程流动方程Advection Equation(AE) 二阶线性方程浙江大学实用数值计算方法47.1 求解方法有限差分法 Method of Finite Differences (MFD)特征线法 Method of Characteristics (MOC)线上求解法 Method of Lines (MOL)有

2、限元素法 Method of Finite Elements (MFE)加权残差法 Method of Weighled Residuals (MWR) 问题 收敛性 Convergence 当采取的步骤趋于无限时，数值结果是否趋于理论值？ 稳定性 Stability 在某一步引入的误差，经多步数值计算后，会扩大或抑制？浙江大学实用数值计算方法57.2 一阶双曲型方程的差分求解法或称流动方程 Advective Advection Equation (AE) v为流速因子该方程的介折解求具体解时需要提供2个辅助条件浙江大学实用数值计算方法67.2 assuming the forcing fu

3、nction is a Rump The solution of is shown below.图 7.1 Propagation of the Wave Front浙江大学实用数值计算方法7 7.2.1 最简单的差分化格式构想图 7.2浙江大学实用数值计算方法87.2.1 以上方法称为 时间镶嵌空间中心 的差分表达 Forward Time Centered Space FTCS represetation 实际上这个方法不能用：不稳定的方法 Unstable Method 考虑数据误差 r 由于原方程为线性，故误差的传播关系 是与原方程完全相同的差分方程 差分方程独立解的一般形式 Inde

4、pendent Solutions of Difference Equations浙江大学实用数值计算方法97.2.1 应为补充解和特殊解之和 补充解系由下式求出补充解系由两个独立解组成浙江大学实用数值计算方法107.2.1 差分方程的解 可用算符运算方法 Operator Calculus 导出 差分算符 Difference Operator它和微分算符一样，是一种线性算符用于线性二阶差分方程和微分方程类似，它的补充解可由下式得到浙江大学实用数值计算方法117.2.1故补充系由两个独立解组成 （Independent Solutions）两个独立解为差分方程的一个独立解（Eigenmode

5、）浙江大学实用数值计算方法127.2.1 差分方程独立解的一般形式用于本题的情况将独立解代入差分表达式得到浙江大学实用数值计算方法137.2.2 差分格式的改进Courant Condition图 7.3图 7.4浙江大学实用数值计算方法147.2.2Courant 条件的物理意义波形传递系沿x=vt线t节点的选取 当节点取在线上： 当节点取在线外： 当节点取在线内： Lax差分格式也写成以下形式可以看成为以下偏微分方程的FTCS差分式dissipative term 耗散项Numerical Viscosity 数值黏度图 7.5浙江大学实用数值计算方法157.3 一阶双曲型方程的特征线求解

6、法Method of Characteristics (MOC) 这是原方程的转换方程，它们的解相同。 为原方程的特征线方程在特征线上，满足 的为解。浙江大学实用数值计算方法167.3.1 Method of Characteristics (MOC)图 7.6浙江大学实用数值计算方法17 7.3.1 Method of Characteristics (MOC)图 7.7浙江大学实用数值计算方法187.3.1浙江大学实用数值计算方法197.4 一阶双曲型方程的线上求解法Method of Lines (MOL) 有限差分法：偏微分方程完全离散成为 一组差分方程 用线性代数方程组求解 线上求解

7、法：偏微分方程部分离散成为 一组常微分方程 用常微分方程积分方法求解浙江大学实用数值计算方法20线上求解法 Method of Lines (MOL)线间距积分步长7.4图 7.8浙江大学实用数值计算方法217.5 二阶椭圆型方程的差分求解法称为稳态热传导方程，通式为 Dirichlet 问题 Neumann 问题浙江大学实用数值计算方法22u(xm,y)=f2(y)u(x0,y)=f1(y)Laplace 方程的 Dirichlet 边界条件和 Neumann 边界条件和 Poisson 方程边界条件也需4个，有3类给定方法 Dirichlet 边 界 条 件 Neumann 边 界 条 件

8、 混合 边 界 条 件7.5图 7.9浙江大学实用数值计算方法23 7.5.1 Laplace算符的差分表达用于Laplace算符浙江大学实用数值计算方法247.5.1图 7.10浙江大学实用数值计算方法25例：Laplace 方程的Dirichlet 边界问题7.5.1图 7.11浙江大学实用数值计算方法26为了提高精度需要加密网络7.5.1图 7.12浙江大学实用数值计算方法27Laplace 方程 Dirichlet边界问题的差分求解 消去法 直接迭代 Liebmann 方法 相继松弛 S.O.R. 方法 交替方向A.D.I.方法7.5.1浙江大学实用数值计算方法28 7.6 二阶椭圆型

9、方程的有限元素法求 Method of Finite Elements (MFE)以Laplace 方程的Dirichlet 问题为例根据变分原则VariationalPrinciples等价性定理以上方程的解将使以下泛函为最小。图 7.13浙江大学实用数值计算方法297.6将D进行剖分，常用的是三角剖分法对任何一个元素用二原线性函数近似在三个顶点上可得到其中浙江大学实用数值计算方法307.6Ui=Wi Uk=WkUj=Wj 图 7.14浙江大学实用数值计算方法317.6所以其中既然顶点坐标均为规定，所以并有浙江大学实用数值计算方法327.6使泛函最小的问题，即对近似为对求极值，或因此得到：可

10、解得n为内部节点数边界上的W为给定浙江大学实用数值计算方法337.6对于更为一般性的情况需要极小化的泛函将是也可剖分为有限个元素后求解图 7.15浙江大学实用数值计算方法34 7.8 二阶抛物型方程的差分求解法动态扩散方程对于一维空间用差商代替微商，可以有各种选择，例如所以有需要另有更方便的方法浙江大学实用数值计算方法357.8显式方法得到或者：则有：图 7.16浙江大学实用数值计算方法367.8示例：取得到的数值解与以下解析解比较饱和蒸汽C2H5OH空气图 7.17浙江大学实用数值计算方法37Number of time stepsAnalytical SolutionsNumerical SolutionsAnalytical versus Numerical SolutionsDiffusion Dynamics r0.257.8图 7.18浙江大学实用数值计算方法38Number of time