对流扩散方程的迎风广义差分格式

一、引言 考虑对流扩散方程的稳态问题 ?这里Ω为R~2上的有界凸域,Ω为其光滑围道;a为常数,c(x)为上的光滑函数,b(x)=(b_1(x),b_2(x))为上的光滑向量函数,并且满足 |b(x)|>>a>0,c(x)≥0,c(x)-divb(x)≥0,x∈.一般情况下,方程(1)的解u(x)在一窄层内迅速变化,用通常的差分法或有限元法计算,将产生严重的振荡失真现象.本文基于广义差分法构造了一类新的迎风格式,它具     

A PRACTICAL PARALLEL DIFFERENCE SCHEME FOR PARABOLIC EQUATIONS

A practical parallel difference scheme for parabolic equations is constructed as follows: to decompose the domain Ω into some overlapping subdomains, take flux of the last time layer as Neumann boundary conditions for the time layer on inner boundary points of subdomains, solve it with the fully implicit scheme on each subdomain, then take correspondent values of its neighbor subdomains as its values for inner boundary points of each subdomain and mean of its neighbor subdomain and itself at overlapping points. The scheme is unconditionally convergent. Though its truncation error is O(τ h), the convergent order for the solution can be improved to O(τ h2).     

拟线性粘弹性方程新H~1-Galerkin最低阶混合元格式的高精度分析

利用双线性元和零阶Raviart-Thomas元,针对拟线性粘弹性方程建立新的H~1-Galerkin混合元逼近格式.在半离散格式下,给出原始变量u的H~1模和应力=?ut的H(div;?)模的超逼近性和超收敛结果.同时,导出向后欧拉格式和Crank-Nicolson-Galerkin格式的最优误差估计.最后,通过数值算例表明逼近格式是有效的.     

伪压缩族公共不动点的复合隐格式迭代逼近问题

设H是一实Hillber空间,K是H之一非空间凸子集,设{Ti}Ni=1是N个Lipschitz伪压缩映象使得F=∩Ni=1F(Ti)≠0,其中F(Ti)={x∈K:Tix=x}并且{αn}n∞=1,{βn}∞n=1[0,1]是满足如下条件的实序列(i)∑∞n=1(1-αn)2= ∞;(ii)limn→∞(1-αn)=0;(iii)∑∞n=1(1-βn)< ∞;(iv)(1-αn)L2<1,n1;(v)αn(1-βn)2 αn[βn L(1-βn)]2<1,其中L1是{Ti}iN=1的公共Lipschitz常数,对于x0∈K,设{xn}n∞=1是由下列定义的复合隐格式迭代xn=αnxn-1 (1-αn)Tnyn,yn=βnxn (1-βn)Tnxn,其中Tn=TnmodN,则(i)limn→∞‖xn-p‖存在,对于所有的p∈F;(ii)limn→∞d(xn,f)存在,其中d(xn,F)=infp∈F‖xn-p‖;(iii)liminfn→∞‖xn-Tnxn‖=0.本文的结果推广并且改进H-K.Xu和R.G.Ori在2001年的结果和Osilike在2004年的结果,并且在这篇文章中,主要的证明方法也不同与H-K.Xu和Osilike的方法.     

有限个非扩张非自映像隐迭代格式的逼近

假设E为一致凸Banach空间,K为E的非空闭凸子集且为E的非扩张收缩,P为非扩张收缩映像.{Ti:i=1,2,…,N}:K→E为非扩张映像且F(T)=∩ from i=1 to N F(Ti)≠■.定义{xn}如下:x0∈K,xn=P(αnxn-1+(1-αn)TnP[βnxn-1+(1-βn)Tnxn]),n≥1,这里{αn},{βn}为[δ,1-δ]中的实序列,其中δ∈(0,1).若{Ti:i=1,2,…,N}满足条件(B),则{xn}强收敛于x*∈F(T).     

二维抛物型方程的一族高精度显式差分格式

构造了一族解二维抛物型方程的高精度显格式 ,其稳定性条件为r=Δt/Δx2 =Δt/Δy2 <1 /2 ,截断误差为O(Δt3 +Δx4)     

不可压缩核废料污染问题的变网格特征有限元法

1 引言 多孔介质中的核废料污染问题是环境保护领域的重要课题。对于不可压缩二维模型,它是地层中迁移型耦合抛物型方程组的初边值问题:     

关于U_t=aU_(xx)差分格式的进一步研究(英文)

对于热传导方程构造了两个高阶精度的差分格式，一个是三层七点显格式，另一个是三层九点隐格式．证明了差分格式的收敛性和稳定性，最后给出数值计算结果．     

自忆模式中差分格式的稳定性研究

基于大气运动是一种不可逆过程的观点，引进了忆及过去时次资料的记忆函数，导出热传导的自忆性方程，研究了方程分别取Richardson和DuFort-Frankel格式，回溯阶p取1时的稳定性.探讨了多时刻模式中数值积分有时发散的问题，揭示了由过去时次资料动态求取记忆函数，改变了原定设计的差分格式，且它是一个时间平滑因子的本质. 关键词： 回溯阶差分格式 多时次 稳定性 记忆性     

前向爆轰驱动变截面激波管特性的数值模拟

讨论的激波管由氢氧前向爆轰产生高压驱动气体、并配有收缩截面段,它能产生高焓超高速气流.对此进行了模拟.既对变截面角度的作用,两端截面积比的作用又对比正向反向爆轰有何不同效果做了探讨.正向爆轰时主膜处的收缩段产生的汇聚作用既加强主激波又产生反向激波,缓解了爆轰波阵面后紧跟着的稀疏波导致主激波衰减偏快的不利影响.当收缩角度在30°和45°之间时,主激波的强度较高,衰减较小;当β接近90°时在主激波传播过一小段距离之后,主激波后高温高速气流较均匀,且主激波的衰减最小,具有实用价值.驱动段与被驱动段的截面积比越大,主激波的强度越高,但是最初阶段衰减也越快.反向爆轰时产生的主激波衰减最缓,但是同样的主激波强度需要的驱动段与被驱动段初始压力比前向爆轰高1个量级.