一类双曲型方程的摄动有限差分格式及其数值分析-perturbation finite difference scheme for a class of hyperbolic equations and its numerical analysis.docxVIP

第一章引言1. 1 问题的引入自然界存在大量复杂的流动现象，如大气运动，潮汐涨落，石油开采，台风，龙卷 风，原子弹爆炸，超新星爆发以及人类的日常生等等。人类的生产生活无不与流体流动 息息相关，随着人类对流体流动现象认识的深入，人类开始利用流动规律改造自然界， 如航空航天，武器研制，工业制造和农业灌溉等。流体运动的规律由一组控制方程描述，在计算机发明前，流体力学家们在对方程经 过大量简化后通过解析方法求得其解，但实际的流动问题大都是复杂的强非线性问题， 无法求得精确的解析解，计算机的出现及其迅速发展使人们直接求解控制方程组的梦想 渐渐得到实现，从而产生了计算流体力学这门交叉学科。计算流体力学（CFD，Computational Fluid Dynamics）是一门用数值计算方法直接 求解流动主控方程（Euler 方程或 Navier-Stokes 方程）以发现各种流动现象规律的学科， 它综合了计算数学，计算机科学和流体力学等多种学科，包括计算水动力学，计算空气 动力学，计算燃烧学，计算传热学，计算化学反应流动以及数值天气预报等[17~22]。 随 着 CFD 技术的飞速发展并广泛应用于各领域，大大减少了实物试验，缩短了研发周期， 节约了研究经费，利用 CFD 进行反复的设计，分析和优化已成为许多应用领域中标准 的必经步骤和手段，同时，各应用领域对 CFD 的需求也极大刺激了 CFD 的发展。在航空、气象、海洋、石油勘探等方面的流体力学问题，在很多情况下都归结成双 曲型偏微分方程（国外文献称为守恒律）的问题，由于这类方程的解中常常出现间断和 振荡（通常称为激波或振荡问题），所以有关的数值解的研究一直引起人们的关注。我们以一维问题为例，简要说明间断和振荡产生的原因，考虑如下双曲型方程?u ??a ?u ??0 ， ??∞???x ???∞，0 ??t ≤ T(1.1)?t?xu(x,0) ? ?(x)(1.2)其中?(x) 给定的充分光滑函数。容易验证当 a 为常数时，方程的解为u(x, t) ???(x ??at) ， ??∞???x ? ?∞也就是说， xt 平面上沿着， 0 ??t ≤ T(1.3)x ??at ? c(c为常数)(1.4)这样的直线， u 的值保持不变。如果把(1.2)看作在初始时刻 t=0 时的“波形”， 则(1.3)表示这个波以速度|a|传播，当 a0 时沿 x 轴正方向传播，当 a0 时沿 x 轴负方向传播，而波形保持不变。dx在 xt 上半平面上，方向为??a 的直线为(1.4)，由于dtdu ???u ???u dx ???u ??a ?u ??0dt?t?x dt?t?x所以 u 在 xt 上半平面上任一斜率为 a 的直线(1.4)上其值为常数。当然 c 不同，即直线不同，u 的值也可以不一样，这簇直线覆盖了整个 xt 上半平面，它可以刻画初值问题（1.1）（1.2）的某些特征，通常称此直线为特征线。在 xt 上半平面上过任一点 (x0 , t0 ) 的特征线 x ??at ??x0 ??at0 和直线 t=0 交点为(x0 ??at0 ,0) , 在这一点上 u 的值是已知的，即0u(x0 ??at,0) ???(x0 ??at0 )由于 u 沿着这条特征线的值为常数，所以u(x0 , at0) ???(x0 ??at0 )因为 (x0 , t0 ) 是任一点，所以u(x, t) ???(x ??at)就是（1.1）（1.2）的精确解。某些方程可求得精确解可以帮助我们验证数值解法的优越性，如果 a 不是常数， 则往往不能得到精确解，所以只能用数值解法。但当 a 不是常数时，而是变量 u 时， 无论初始函数多么光滑，初值问题（1.1）（1.2）的解不再是连续可微的，这个事实在流 体力学中有重要意义，它反映了自然现象从连续到间断的转折。因此，如何计算间断解 是人们关注的课题，文献[26]中指出，只能用一阶的显式格式计算间断解，如果采用高 阶显式格式（包括二阶），那么数值解在间断点处会出现振荡现象，因此，在间断点如 何提高解的精确度，引起人们极大的兴趣。探讨这类问题始于二战时期，Join Von Neumann 试图用所谓激波捕捉法去“捕捉” 激波，虽然并不很成功，但他的思想对解决这类问题起了重大影响。激波捕捉法即无论 双曲问题是否有间断解，都不加区别地按统一的格式进行计算。采用激波捕捉法间断解 是作为解的一部分进行计算的，间断解被表示为具有一定过渡层的连续解，这个连续解 基本上能反映出本来的间断特征，这类方法最大的优点是简单。对一维双曲型方程：??u ??u ?u???t?x（1. 5）?0?u( x, 0) ??u ( x)众所周知，即使 u0 (x) 是 x 的充分光滑函数，式（