一类特殊对流扩散方程有限差分法的研究

来源 :新疆大学 | 被引量 : 0次 | 上传用户:aa283488665
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
对流扩散方程是描述粘性流体运动的非线形方程的线性化模型方程,而且它本身也描述了许多自然现象,例如在水中或大气中污染物质浓度的扩散,沿海盐度,温度扩散等等。因此求解对流扩散方程的计算方法引起了充分的重视。 考虑对流扩散方程在区域D={x|0≤x≤1,0≤t≤T}上的初边值问题初始条件 u(x,0)=u0(x) 0<x<1边界条件 u(0,t)=u(1,t)=0 0<t≤T其中ε>0为任意常数。 本文对问题(1)的有限差分方法进行了研究。从一般差分格式的构造,并行解法的构造,特征-差分法的构造三个方面出发讨论了对于对流项系数无界,扩散项系数趋近于零情况的处理方法。 目前,求解常系数线性对流扩散方程已有了许多方法和技巧,理论上也很成熟。一般地,这些方法和技巧都可以推广到求解非线性对流扩散方程,在对流项的处理上以及对计算小ε值的情形处理成为一个关键。对于对流项系数有界且对流占优的扩散方程的讨论较多,采用的方法很多,主要以有限元与特征法的结合为主。而对于对流项系数无界的情形,讨论的较少,所采用的方法也只是一般差分方法。无论是哪种情形,若对于小的ε值处理不当,计算结果则会出现非物理振荡。所以对于此类问题的研究工作一直方兴未艾。 本文从三个方面将一般差分方法,并行算法和特征-差分法求解常系数线性对流扩散方程已有的一些方法和技巧推广到了求解变系数线性对流扩散发展方程。在求解问题(1)的一般差分方法中,我们得到了两种绝对稳定的加权隐格式,分别是4点加权两层隐格式和6点加权两层隐格式。得到了三种条件稳定的显格式,分别是4点加权两层显格式,4点加权迎风显格式,和5点加权三层显格式。文中分别给出了这五种差分格式的稳定性条件,并利用能量分析法给出了4点加权两层显格式的稳定性证明。 虽然建立了这些差分格式,但在并行算法的推广中遇到了问题,显式方法很适合于并行计算,但由于稳定性条件限制,必须采用非常小的时间步长 进行计算。隐式格式一般无稳定性条件,但在每一个时间层要求解线性代数 方程组,实现并行计算有一定的困难。本文将Evans和Abd毗lah首先建立的 GER方法,GEL方法,。4GE方法(参见文献p6I)推广到了求解问题(1)。 这三种格式同样能进行显式计算,在实际中具有好的稳定性和计算精度,但 在理论上不能保证其绝对稳定性质,原因在于问题n是关于时间t的发展方 程,若问题(1)的系数与时间t无关,则AGE方法是绝对稳定的(参见文献 [25])。 另外,文中也给出了两种特征-差分方法,它们是差分法与特征法的结 合,即扩散项用差分法,对流项用特征法,大体上都是迎风格式在不同意义 下的推广。其中第一种特征-差分方法为通常所采用的一般格式,形式比较复 杂,数值求解比较困难。第二种特征-差分方法是在第一种方法基础上的改 进,截断误差较小,形式简单,数值求解比较容易。文中给出了第二种方法的 数值结果。 最后,我们针对以上讨论,进行了以下初边值问题的数值实验 r 纱一卫犁=三染.0<、。<1.才>0 D & Tof“a22·——、、-。,。-- lull.0)=1一1“.0<上<1 】犁=0。=0.亡>0 Ill=0.工=1。亡>0 这里。是充分小的常数。 为了便于说明问题,我们以文中的三种构造方法,分别进行数值实验。在 每组数值实验中,我们分别取了4个不同的。值,并将数值结果绘制成图形。 在第一种构造方法中我们将分别利用以【“三种差分格式进行求解 川含有加权因子的两层4点迎风显格式 间两层4点隐格式 (3)含有加权因子的两层6点隐格式 在第二种构造方法中我们将分别利用以下三种差分格式进行求解 (l)肌ul’yev格式,采用从右向左求解. (二)GER格式求解对流扩散方程. (3)AGE格式求解对流扩散方程. 在第三种构造方法中我们采用文中式(4印的差分方程组进行求解。 综上,数值结果表明了文中三种构造方法的可行性。
其他文献
早在上世纪初,伴随着生物学与数学的发展产生了一门新的边缘学科-生物数学.所谓生物数学,顾名思义就是利用数学学科中的一些理论去解决生物学中的某些实际问题.其中,微分方程的
该文试用鞅论与随机微积分中的相关理论,研究了有交易费用情况下的欧式期权定价问题,得多维期权定价公式,可以用数值方法得到期权定价问题,得到了多多定价公维数为1时,文中的
本研究如下几乎临界增长的半线性椭圆型方程组解的渐近行为:{-△u=| x|βvqε, x∈Ω,-△v=| x|α upε, x∈Ω,(0.1)u=v=0, x∈(6)Ω,其中Ω(C) RN为单位球,ε>0,α,β>0,pε=p-(p
该文研究的是动力学系统伪谱与扰动系统谱之间的关系,侧重于研究非正规系统的谱问题,该文对常微分方程系统和微分--代数方程系统进行了较为松弛算子的谱进行了分析.数值实验
具有许多优良特性的k元n方体是应用非常广泛的互连网络之一.k元n方体Qkn(k≥2,n≥1)的顶点集V(Qkn)={u0u1…un-1:0≤ui≤k-1,0≤i≤n-1},两个不同的顶点u=u0u1…un-1和v=v0v1…vn-
该文第一部分通过对期权的标的资产(通常中股票)的价格行为分析,引入在布郎 运动和泊松过程共同驱动下的支付股标红利的新股票价格过程,即就是改变了经典black-Scholes模型其
本文主要讨论非线性时滞微分方程的概周期解的稳定性.全文分为三章,所得结果推广和改进了文献中的相关结论.  第一章,主要介绍了非线性时滞微分系统的概周期解的性质的研究发
该论文研究两上问题,一个是科技进步的数量度量方法,另一个是收入分配不均等性数量度量方法.在前一个问题中,该文讨论了许多经济学家在应用中均使用的一个基本假设(规模收益
该文较系统的研究了一个带有关闭延迟期和启动期的GI/M/1排队系统.这是一个全新的休假排队模型,田乃硕(1992,1997)及Dukhovny(1997)等所研究的模型都是该文中模型的特例.使用
该文针对当今电力营销系统存在的主要问题,构造了一种基于Java的三层Browser/Server模型,在此基础上,依据新型电力营销系统设计规范,将电力营销系统划分为三个层次,即客户服