一阶双及曲型方程及方程组的数值边界处理.pptx

5.3 一阶双曲型方程及方程组的数值边界处理 我们知道一阶双曲型方程不是每个边界都给定条件的，因此差分方程所需要的边界条件往往比微分方程的边界条件要多，由此必须加上附加的边界条件，亦称数值边界条件。如果这样的条件处理不好，就会影响到内部网格点的计算，从而导致整个计算的不稳定。在本小节中我们仅以例子作说明而不进行一般性的讨论。例: （5.17）采用蛙跳格式: 初边界条件的离散取 (5.18)(5.19a) 由于蛙跳格式（5.17）是一个三层格式，因此初始条件中还需加上 的值。注意到特别取 t=0, 并假定g(x)一次连续可微，则有这样，我们得（5.19b）我们假定利用(5.19a)及（5.19b）可以得出. 这样利用差分格式（5.17）和边界条件（5.18）可以得出.接下来要计算 ，当计算 时就必须用到 ，这个值不可能由差分格式（5.17）给出。于是，对于差分格式(5.17)来说必须增补x=0处边界条件。数值边界条件（人工边界条件） 外推法: 由内点值外推边界值(5.20)如果选取特殊的初值考虑条件（5.20），那么容易得到其中 为一个常数，，由此看出选取（5.20）作边界条件是不行的。如果我们用数值边界条件（5.21）来代替（5.20），其结果也是不行的，也是不可用的边界条件。下面给出两个可行的附加边界条件（5.22）（5.23）条件（5.22）就是利用迎风格式进行边界处理，一般来说这总是可行的。