理论基础2.docVIP

七、差分数值解的耗散(Dissipation)和频散(Dispersion)性质 微分方程解的耗散或频散特性； 考查方程： 的解的性质，———（2－7－1） 设初值条件为： 设解的一般形式为： （2－7－2） 其中 为复数，并且 若写成 其中为实数 则： 将（2－7-2）式的解代入（2－7－1）方程，得： 即 将 代入: （2－7－3） 讨论几种特例； 所有 即 ————行波方程 解为: 这表示分解的振幅将终是 ，所以解是无耗散的；不论k是多少， 波速均为所以波形不发生任何变化，所以也是无频（色）散的解。 ，其余的偏导数项系数均为0。方程为： 解为： 此解表示 仅当时解的振幅才是衰减的，即为（正）耗散解 而 当 时 振幅随时间呈指数增长，解将是无界的，有时也称为负耗散解。 但不论 ，波速与波数k无关，波速恒为，所以是一种无频散的解。 ，其余的偏导数项系数均为0。方程为： 方程（弧立波方程） 解为： 解是无耗散的，但不同的分立波（波数不同），传播速度不一样，其传播速度为 即 时k越大，则波速越大，换言之，波数较大的子波会逐步赶上波数较小的子波，在满足一定的条件下，将形成孤立波。 类似地讨论 iv).时 : 结论: （1）偶阶导数项影响解的耗散，并且对于能被4的整除的偶阶项， 当其系数为负时，是正耗散，为正时是负耗散（解趋于无穷）；而其余的偶阶项（即不能被4整除的偶阶项，例如2阶、6阶、10阶等···）当系数为正时是正耗散，系数为负时是负耗散（解）； （2）奇阶导数项只影响解的频散（色散）特性，不影响解的耗散特性； （3）方程（2－7－1）的解既有耗散，也有频散，其耗散及频散特性与 这两个无穷级数的和有密切的关系。 解的耗散，频散特性的定量讨论方法。 例1； 解为 当时，解可以写成： 对于的每一个分量（即分立波）；引入放大因子G；G是复数 对于复数G，可以考察，以及G的复角G ： 的含义；是相邻时间间隔内解的振幅的改变； G的含义是相邻时间间隔内解的相位差的改变。 例2； 解为； 放大因子 当时 由于只要求G，所以并不一定需要将的仔细形式写出，而重点放在这两个瞬时的解表达式 例3；对于差分方程 例； 的格式， 该差分格式的修正方程为（通式）： 回阅（2－7－3）式，其解为： 初值 若，解可写成： 注意此时 都是复数，其含义不仅仅是据幅了（与的含义不同！！） 从 的改变包含振幅和相位的改变： 因此，若固定空间位置，考虑时间间隔前后的解之比： 另一方面（2－7－3）式的解，当空间位置变化时， 即：当时 ，或时 相应的解的形式改变是： 综合以上： 对于线性差分格式（2－7－4）式： 将此假设代入差分格式（2－7－4），并考虑对于线性差分格式，可以分别讨论每一个Fourier 分量的关系，有： 所有项均有并同除以 有时习惯将写成，所以放大因子： 显而易见，Lax差分格式的解与源方程的解的特性存在差异。 与相比 讨论：对于任意，要求则充分和必要的条件是：或 就是Lax格式的稳定条件。 从另一角度看， 即使，保证了解差分数值解的有界，（稳定了！）但数值解与真解的差别，仍存在着耗散和频散这两个方面的误差。 可通过下列图示表示： 八、差分格式的守恒性质； 如果对一个差分方程在定义域的任意有限空间内作求和运算，（即相当于在连续问题中对微分方程在空间域中作积分运算）所得的表达式仍能满足该区域上物理量的守恒关系，则称该差分格式具有守恒性，或守恒格式。 例；对于连续方程； 有限体积域内的质量守恒律 为简单起见，讨论一维问题；——守恒型 ——非守恒型 对守恒方程用FTCS格式； 若从N至M累加 可见，该格式在离散的概念下，所描述的守恒关系与微分的源方程描述的守恒律是一致的，所以是守恒格式。 例2：Burger方程，其守恒形式是： 格式1（由守恒形式出发）： 格式2（由非守恒形式出发）：