扩散方程扩散.ppt

计算物理 有限差分方法 物理问题和数学方程(1/5) 物理问题和数学方程(2/5) 物理问题和数学方程(3/5) 物理问题和数学方程(4/5) 物理问题和数学方程(5/5) 有限差分原理(1/3) 有限差分原理(2/3) 有限差分原理(3/3) 矩形区域中的泊松方程(1/1) 迭代解法(1/6) 迭代解法(2/6) 迭代解法(3/6) 迭代解法(4/6) 迭代解法(5/6) 迭代解法(6/6) 非矩形区域中的泊松方程(1/3) 非矩形区域中的泊松方程(2/3) 非矩形区域中的泊松方程(3/3) 一维扩散方程(1/3) 一维扩散方程(2/3) 一维扩散方程(3/3) 二维扩散方程(1/3) 二维扩散方程(2/3) 二维扩散方程(3/3) 一维波动方程(1/3) 一维波动方程(2/3) 一维波动方程(3/3) 作业 * * 3/lesson/ComputationalPhysics 有限差分方法 物理问题和数学方程 有限差分原理 矩形区域中的泊松方程 迭代解法 非矩形区域中的泊松方程 一维扩散方程 二维扩散方程 一维波动方程 √ 物理问题 连续体(气体、液体、固体)和场(电磁场、温度场) 状态用关于时间和空间的函数描述：u = u(x, y, z, t) 系统遵循特定的物理规律，即 u 的变化满足特定的(偏微分)方程 数学方程的建立 稳定过程：泊松方程 静电场：电场的散度正比电流密度，等于势场u的梯度 静磁场：类似于静电场 √ 输运过程：扩散方程 扩散：流体由于不均匀而发生扩散，扩散密度正比于密度的梯度，同时满足质量守恒 热传导：类似于扩散 振动传播：波动方程 交变电磁场 机械振动 √ 方程的分类 二阶线性偏微分方程的一般形式 椭圆型(B2-4AC0)：二维泊松方程 抛物型(B2-4AC=0)：一维扩散方程 双曲型(B2-4AC0)：一维波动方程 √ 边界条件和初始条件 物理：过程的具体特征与初始状态和边界(受外界影响)有关 数学：偏微分方程有无限个解，定解需要初始和边界条件 边界条件 第一类 第二类 当 u 是电磁场的势，则 代表场强 当 u 是密度/温度/位移，则 是流量/热流/应力 当 q0=0 时称为第二类齐次边界条件 √ 第三类 例：热传导，系统通过表面与外界交换热量：表面热流 正比于表面温度 u 与外界温度 u0 之差，即 初始条件 初始瞬间待求函数 u 在各处的值： 初始瞬间 u 在各处对时间的变化率： 求解方法 大多数没有解析解，只能数值求解 求解方法：变分解法、有限元法、有限差分法，边界元法 √ 差商格式 基础：用差商代替微商(离散化) √ 误差为O(h)的差商公式 √ 误差为O(h2)的差商公式 差商格式的收敛性和稳定性 收敛性：当步长 h?0 时，差分方程的解是否收敛于微分方程的解 稳定性：误差 Du 在运算过程中是否失控，即累计误差是否无限增加 √ 五点差商格式 二维泊松方程 收敛性和稳定性：在给定边界条件下具有唯一解，当 h?0 时趋于解析解 拉普拉斯方程 方程( f(x,y)=0 ) 五点差商格式 √ 差分方程组的特点 方程个数等于内点数，每条方程最多含5个未知项 系数矩阵是稀疏和带状的 跌代法求解：同步法、异步法和逐次超松弛法 同步法 用第 k 步的 u 值，代入方程的右边，计算得到新的 u 值，用在第 k+1 步，…；直到新旧的 u 差值小于设定的误差 特点：需要两套内存(存放 u )，收敛慢 异步法 计算第 k+1 步的 ui,j 时， ui-1,j 和 ui,j-1 已经知道 特点：需要一套内存(存放 u )，收敛较快 √ 逐次超松弛法 将第 k 步的 u 值，与异步法中计算得到的 u 值，加权平均，用在第 k+1 步 特点：加权因子 w 对跌代次数的影响很大。经验公式 √ 矩形区域的第二和三类边界条件 当 a 和 b 是 x, y 的函数时，应 a=a(xi, yj) 和 b=b(xi, yj) 对第二边界条件，令 a=0 √ 不规则区域 第一类边界(不对称网格方法) 第二类边界 结点在边界上 结点不在边界上：过结点 P 向边界作垂线，交于 P 点，以 P 代替 P 第三类边界 前两类边界条件的组合 √ 例：求解拉普拉斯方程 √ 例：求解拉普拉斯方程，边界条件如右图所示 取 h=5，有三个内点：uA, uB, uC 转化为线性方程组 方程的解 收敛性 26.095 5.491 1.095 5/32 26.097 26.107 26.144 26.289 26.786 uC 5.498 5.525 5.632 6.019 7.143 uB 1.097 1.107 1.144 1.289 1.786 uA 5/16 5/8 5/4 5/2 5 h √ 二维圆