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

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