第一章有限差分法.pptVIP

定解条件的离散化 h ? L 3 ? 3 第三类边界条件 的差分离散化 ? 0 0 n y D ? ? ? ? ? ? ? ? ? ? f ( r , t ) ? f ( r , t ) 1 2 ? ? ? n ? S ? ? 第一种情况，当结点 刚好着落 于边界线 L 上时，这还取决于边 界结点处的 外法线与网格线 重合， x ? 0 ? f 1 ( r 0 ) ? 0 ? ? 3 h ? f 2 ( r 0 ) 21 定解条件的离散化 ? 外法线与网格线不重合情况，边界结点上的外向法向方向 与水平夹角为 ā ，其法向导数显然是在 x 和 y 方向的导数在 法向的 投影组合 ， L ? ? ? ? ? ? ? ? ? ? cos ? ? sin ? ? ? ? ? ? ? ? y ? ? n ? 0 ? ? x ? 0 ? 0 ? ? 3 ? 0 ? ? 2 ? ? cos ? ? sin ? h h 3 0 ? 0 ? n h y D 2 ? 0 ? ? 2 ? ? 3 ? ? 0 ? ? 0 ? f 1 ( r 0 ) ? cos ? ? sin ? ? ? f 2 ( r 0 ) h ? h ? 22 x 定解条件的离散化 h y L 3 o o ? ? 0 n 2 D x ? 第二种情况，当结点不落于边界线 L 上时，只需要 引入 于 结点 0 相关的边界结点 O‘ ，点的外方向 n 作为结点 0 处的 “外方向 n” ，且 近似地认为 边界条件中给定的函数和均在 O 点上的取值。这样，此种情况下的 第三类边界条件 的离 散格式于式相似， ? 0 ? ? 2 ? ? 3 ? ? 0 ? ? 0 ? f 1 ( r 0 ) ? cos ? ? sin ? ? ? f 2 ( r 0 ) h ? h ? 23 定解条件的离散化 第二类边界条件 的差分离散化 ? ? / ? n ? 0 ? 第二类齐次边界条件为第三类边界条件的特殊情况，即。 我们这里讨论最常见的一种情况 f 1 ( r ) ? f 2 ( r ) ? 0 加一层虚拟边界 ? ? 1 ? ? 3 1 2 ? 0 ? ? 2 ? 1 ? ? 2 ? ? 4 ? h F ? 4 ? L ? 2 3 h D 2 0 ? 3 ? 0 4 1 ? 1 上面也是 对称边界条件 的离散公式 ? 4 24 有限差分法的求解 ? ? 综上所述，对 场域 D 内各结点 （包括所有场域内结点和边 界结点） 逐一列出对应的差分计算格式 ，即构成以这些离 散结点上的 位函数 为待求量的差分方程组 （代数方程 组）。求解这些代数方程组，得到场域中的电位值 计算步骤通常是： ? 离散场域 ，采用一定的网格剖分方式离散化计算区域。 ? 离散化场方程 ，即基于差分原理的应用，对场域内场的 偏微分方程以及定解条件进行差分化处理，得到方程的 差分格式。 ? 计算离散解 ，建立的差分格式（与原定解问题对立的离 散数学模型 — 代数方程组），选用合适的代数方程组解 法，编写相应的计算程序，算出待求的结点上场值。 ? 25 有限差分法的求解 电 磁 场 问 题 前 处 理 离 散 化 场 域 （ 网 格 剖 分 ） 离 散 化 方 程 （ 差 分 原 理 ） 差 分 方 程 组 （ 代 数 方 程 组 ） 数 据 计 算 计 算 方 程 组 （ 迭 代 法 ） 离 散 解 插 值 计 算 其 他 值 或 可 视 化 显 示 结 果 后 处 理 26 有限差分法格式特点 仔细分析离散的差分方程组，例如泊松方程，从离散方程 式不难看出，该方程组的系数一般是 有规律 的，且方程都 很简单，每个方程的项数不多（待求量最多 不超过 5 项 ） ? 各离散结点上的方程组形式 （结点顺序按坐标 先从 y 轴增加 (i-1,j+1) (i,j+1) (i+1,j+1) j+1 、再 x 轴增加 （ 从下到上、从左 到右，即先列后行 ）排列 ? ( i , j ) j ? j increase (i-1,j) (i,j) (i+1,j) j-1 (i-1,j-1) (i,j-1) (i+1,j-1) i-1 i 27 i+1 i increase 有限差分法格式特点 ? ? ? ? 2 ? ? i ? 1, j ? 1 ? ? i ? 1, j ? 1 ? ? i , j ? ? i , j ? 2 ? 4 ? i , j ? 1 ? h F i , j ? 1 ? 2 ? j i ncrease ? ? i ? 1, j ? ? i ? 1, j ? ? i , j ? 1 ? ? i , j ? 1 ? 4 ? i , j ? h F i , j ?