一阶线性常系数双曲性方程的有限差分方法的研究2012.5.3.doc

Forpersonaluseonlyinstudyandresearch;mercialuse引言主要讨论双曲性方程及双曲性方程组的差分方法。从简单的一届线性双曲型方程开始,构造差分格式,分析其稳定性及其他性质,然后推广到一届线性双曲性方程组。双曲方程与椭圆方程,抛物方程的重要区别,是双曲方程具有特征和特征关系,其解对初值有局部依赖性质。初值的函数性质(如间断,弱间断等)也沿特征传播,因而解一般无光滑性,迄今已发展许多逼近双曲方程的差分格式,这里只介绍常见的九种方法,讨论了各种求解方法,分析了其性质,最后对初边值问题及二维问题进行了讨论。1一阶线性常系数双曲型方程先考虑线性常系数方程,,t>0其中a为给定常数,这是最简单的双曲型方程,一般称其为对流方程。虽然()式非常简单,但是其差分格式的构造以及差分格式性质的讨论是讨论复杂的双曲型方程和方程组的基础。它的差分格式可以推广到变系数方程,方程组以及拟线性方程和方程组。对于方程()附以初始条件u(x,0)=u(x),在第一章中讨论了初值问题(),()式的解,其解沿方程()的特征线是常数,并可表示为下面讨论双曲性方程的应风格式,Lax-Friedrichs格式,Lax-wendroff格式,Courant-Friedrichs-Lewy条件利用偏微分方程的特征线来构造有限差分格式,蛙跳格式,数值例子。,从而产生了很多好的方法和技巧。迎风各式的基本思想是简单的,就是在双曲型方程中关于空间偏导数用在特征线方向一侧的单边差商来代替,()式的迎风各式是,a>0的截断误差和稳定性:,a>0/(两边乘于),得‚+‚所以截断误差为迎风格式对一阶精度,,()的稳定性:先把差分格式变化为便于计算的形式0其中网格式令则=4当时原差分格式是稳定的。所以迎风格式()是条件稳定。根据Lax等价定理,:,……两边除于,得……(两边乘于)得…‚+‚=…所以=…截断误差为迎风格式的稳定性:将方程改变便于计算的形式:令[]所以()21222122144当差分格式时()是稳定的例讨论差分格式的截断误差和稳定性解截断误因为所以[]‚+‚T(x,t)[]所以差分格式的截断误差为即差分格式是一阶精度的。讨论它的稳定性:先把差分格式改写为令并将其代入时有由于a<0所以取=0差分格式是绝对不稳定的。-Friedrichs格式讨论逼近对流方程()的一个中心差分格式的截断误差和稳定性Lax-Friedrichs格式的截断误差:…(两边乘于)得‚+‚