关于前向-后向热方程的数值方式.pdfVIP

摘 要 本文主要研究一类特殊的热方程，前向．后向热方程的数值方法，包括差分 方程的构造，误差估计，数值求解等问题，研究分别对一维问题和二维问题进行 展开。． 差分方法是解偏微分方程定解问题的常用近似方法之一。文献[22】给出了 一维前向．后向热方程的一种差分格式，即在两个子区域分别应用前向和后向差 分格式，而在交界线应用二阶中心差分格式，本文对此差分格式进行了改进，在 交界线上应用粗网格的中心差分格式。 令让为r．-J题(2．1)的准确解，Z为差分方程组(2．23)-(2．25)的解，定义误差目= u(ih．，歹下)一霉，则有： 定理2．2令Q是Q的闭包。如果云1J酉02uf，壶J象J在Q上有界，且有界常数 为G，则 lIEIIⅣ≤妄G(7-4-h24．-H3)， (1) 其中，lIZll^，=ma,xC{h．jr)“l层|． 因为本文构造的差分格式对前向．后向热方程而言是隐格式，因此我们用基 于区域分解的迭代方法进行数值求解。 解，4(1≤J≤N一1)是差分方程组(2．23)一(2．25)的解，则有 (2) 1臼msaⅣx一1(1名一妒七+1·’1)≤(1一日)19m≤aⅣx一1(1晶一Ck．J1)， 即当七一CO时，矿·七收敛于晶，且收敛率为1一日． 这个结果比文f22】中的结果要好。同时，我们对此方法在二维情况下进行了 推广，得到： 定理3．1假职互1l丽02ul’壶I嘉|’壶I券}在Q的闭包上有界，贝IJ lIEllⅣ≤2Co(7-4-h24-H3)， (3) 其中，岛为有界常数， ㈣Ⅳ2(蚺ma≈rx—E针k 摘 要 ii 定理3．2对1≤J≤M一1，1≤k≤N一1，P=0，l，·．．．，有 Ⅳ ． (4) 显格式容易实现，但稳定性条件限制了时间步长的取值，而隐格式虽然无 条件稳定，但需要在每个时恻步求解代数方程组，Saulyev格式则很好地弥补了 这些缺陷，它是无条件稳定的，而且对一般的热方程它是显式的。 本文详细讨论了将Saulyev格式应用于前向．后向热方程的数值求解，并在 此基础上，构造了前向．后向热方程的分组格式。但是对于前向．后向热方程，这 些格式都是隐式的，所以本文仍考虑用迭代法求解，并且证明： 定理2．4当口∈(0，1)时，给定的区域分解算法收敛。 关键词：前向．后向热方程，差分方法，区域分解，Saulyev格式，分组格式，迭 代方法 Abstract Inthis thenumericalmethodsofthe paper，westudy forward．backwardheat inonedimensionand equation two thedifference dimension，includingmethods， theerrorestimateand numerical solving． Thedifference methodisoneofthe methodsforthe dif- a