五点差分格式算法在数学物理方程中的应用.pdfVIP

第 35卷第 4期 辽宁师范大学学报 (自然科学版) Vo1．35 No．4 2012年 12月 JournalofLiaoningNormalUniversity (NaturalScienceEdition) Dec． 2O12 文章编号 ：1000—1735(2012)04—0467—04 五点差分格式算法在数学物理方程中的应用 林 爽 ， 张 杰 ， 姜行健。 (1．大连工业大学 信息科学与工程学院，辽宁 大连 116034；2．辽宁师范大学 数学学院，辽宁 大连 116029； 3．大连理工大学 机械工程与材料能源学部，辽宁 大连 116024) 摘 要：许多物理现象中的稳定过程都 归结为椭圆型微分方程，差分法是解椭 圆型微分 方程的重要方法．首先给出了解椭 圆型微分方程五点差分格 式的算法框架，然后对数学 物理 问题 中的热传导方程及附加条件作差分逼近，求 出其数值解，最后将理论应用于实 际中去，在实践中得到可行性的检验． 关键词 ：椭圆方程 ；五点差分格式；热传导方程 中图分类号：O241．82 文献标识码 ：A 许多物理现象中的稳定过程都可以用椭圆型微分方程来描述 ，如辐射热传导方程 ，静电场的电势， 稳定温度场的温度分布，流体的流动等．通常情况下我们只能求得其数值解 ，即求出函数的精确解在 自 变量一系列离散点处的近似值．差分法是解椭圆型微分方程的主要数值方法，如五点差分格式、九点差 分格式等．本文主要介绍五点差分格式算法在数学物理方程中的应用． 1 五点差分格式数值解法 在许多工程问题 ，例如热传导、电磁场等问题都可以归结为 Poisson方程边值问题 f—Au=f(x，)，(z，)∈n， (1) 【“一 (z，3，)，(，)∈r， (2) 其中是二维有界区域 的边界，“是z，Y的函数“(z，3，)，△一寿+ 是二维Laplace算子P·。is— son方程 (1)是最简单最典型的椭 圆形方程． 在平面上取一个 固定点(z。，Y。)，过该点分别以步长 ^和r作垂直于z轴和 轴的两族平行线 z —z。+ih，一0，1，… ，N1，— 。+_r『，一0，1，… ，N2，将区域分成矩形 区域 Q∈((，3，)Izo≤z≤ o+Nlh，Y。≤ ≤ 。+N2r)，记 (z，YJ)为 (i，)． (，)是正则内点，应用 Taylor公式 ，得到其五点差分格式 一 寿[㈩，』--2u+H．]一专[地 1--2u+ 广1]= ， (3) 其中 “d是 “(，_)『的一个近似值 ， 一厂(i，)． 如果 r=h，则差分方程 (3)化为 1 壶(4 一“卜I，一f+1，一“一1一“+1)= ， (4) 收稿Et期 ：2012—07—01 基金项 目：国家 自然科学基金项 目 作者简介 ：林爽 (1981一)，女 ，辽宁大连人 ，大连工业大学讲师 ，硕士 468 辽宁师范大学学报(自然科学版) 第 35卷 若 _厂三0，则为 Laplace方程或调和方程 [2] “===÷(一1，+Ui+1，J+ ，，』一1+ul，+1)． (5) (，_)『不是正则 内点时，可通过转移或线性插值推导其差分格式． 2 五点差分格式数值解法的应用 2．1 热传导方程的应用 我们考察三维空间 一(z，z：，X。)中由曲面 I1所界的某体积G内的定常热分布．热传导过程由傅 里叶定律确定 ：热流密度 向量 W 正比于温度 “的梯度，W一一