一个简单的方程组实验模型及应用.pdfVIP

一个简单的方程组实验模型及其应用 孙红兵 方燕 刘儒勋 (云南工业大学数学系，昆明650051)(中国科技大学数学系，合肥230026) 摘 要本文提供了～个二元的，并且给出了它的精确解的方程组模型。通过对它的模拟 实验和精确解的比较，可以更有效地改造、改进和选择合理的差分格式，同时，还可以利用模型 方程直接进行有针对性的差分格式的设计。最后，利用它作了两个方波的碰撞模拟实验，效果很好。 关键词 有限差分格式，实验模型方程，精确解 l 引 言 针对一个实际问题，设计、构造和选择了一个有限差分格式，一般说来需要对它的 可行性和效能进行模型实验，以求作到心中有数。也就是说，在进行数值模拟之前，常 常先要做数值模拟实验，这就需要～个实验的模型方程。目前，～般说来所采用的模型 方程，存在着简单化的问题。 众所周知，在现实中的数值模拟问题多半是方程组形式的，诸如拽水波方程，水锤 模型Navier-Stokes方程，等等．在国内外，不例外地都是利用了一元的各种简化模型 Ⅳl+Q／,／，=0 坼+辑)，=o 甜f+till￡=矿Ⅳ矗 它们分别是匀流方程(又称单向波方程)，无粘Burgers方程(也称激波模型方程)和 Burgers方程。然而，它们终究不同于实际问题的多元形式。这种差异很可能保证不了模 型实验的目的。本文利用特征线理论，提出了一个简单的二元的方程组模型，并给出了 它的解析解。从而解决了多元的实验模型问题。 在论文的第2节，讨论了这种模型方程的基本特征，特征化形式，以及精确解；第3 节里，我们利用这个模型，进行了方波的相互碰撞的数值模拟实验，数值结果表明了数 值解与精确解的吻合程度，以及应当对格式进行修改的要求。譬如，差分格式的数值耗 散性和数值色敬性，等等。 2 特征化理论和模型方程 首先，我们来考察简单的单行波方程 ·48· Il Sf ≥ Of 0 (21) 甜 +0x I| f 、J “ C0 JI (22) 0 甜 Ⅳ ，，L r”) ^“扛功 利用特征线方法，我们不难证明它有准确解(参见【l】)。 1)斗2+2)适定，则有解析解 至埋l 设Cauchy问题(2 甜(x，f)=甜。(x—df)+J：：白(x一口(f—f)，f)dr (23) 1)吖2．2)+(24) 地2设初边值问题(IBVP)(2 玎(O，f)=甜f(，) (2．4) x≥af f“。(X--at)+￡矗(r一口(，一r)，r)dr (2．5) r谢 ”‘x，7’。{“，(≠一言)+正：向(z—d(，一r)，r)dr 【证明】如图1所示， 定义特征线 (26) c：idx=口甜X--aI=C07ISt 州方程(2、1)可以表示为