求解二维反应扩散方程的局部间断有限元分裂算法研究

来源 :南京大学 | 被引量 : 0次 | 上传用户:youyou306
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文主要讨论高维发展型反应扩散方程的局部间断Galerkin有限元(LDG)方法的时间离散技术。针对抛物方程特性,隐式时间推进是常用的方法。这势必造成每个时间层上大的方程组求解,导致LDG方法的计算效率下降。为此,本文将采用所谓的算子分裂技术来降低LDG格式的数据存储量和计算速度。算法主要包含两部分:首先分裂高维扩散项与反应项部分,然后采用局部一维化(LOD)技术求解高维扩散部分。在时间分裂上,我们讨论了Backward-Euler(SBE)分裂,Strang分裂,以及Richardson外推技巧,在时间方向分别达到1阶、2阶和4阶精度。同时,我们还讨论了非周期边界(Dirichlet以及Neumann边界)的中间时间层的边界设置问题,恢复了简单处理所丢失的时间精度。数值实验表明上述分裂方法获得了相当好的数值计算效果。最后我们数值模拟了一些典型应用问题,如急速冷却问题、二维燃烧问题以及Allen-Cahn方程。
其他文献
广义中心三项式系数Tn(b,c)表示关于x的多项式(x2+bx+c)n其展开式中xn的系数,(其中b,c为整数,n为正整数)。下设p>3为素数,m为正整数且m(≡)0(mod p)。2011年,孙智伟教授观察到∑p-1k=
本文基于复杂网络和投资组合理论相结合来研究投资风险,为解决风险管理问题提供了一个新的视角。主要内容包括:⑴选取道琼斯中国88指数和美国金融市场的股票数据,利用边相关系数
The Legendre orthogonal functions are employed to design the family of PID controllers for a variety of plants. In the proposed method, the PID controller and t
本文主要研究了如下一类带衰减位势的非线性Schrdinger方程变号解的存在性:其中,1
设f(x)是整系数上m(≥2)次不可约多项式,n是大于等于2的整数。根据Siegel等人的一些工作,如果(m,n)≠(2,2),我们知道丢番图方程f(x)=yn, x,y都是整数,只有有限个整数解。而现在研究的
本文研究推广的坐标松弛法xk+1=xk-ωkt1hp1h求解非线性方程组F(x)=0,x∈D的问题,其中F:D()Rn→Rn是非线性映射。本文获得以下3个的结果:   1当F是非梯度映射,且迭代次序i
在这篇文章中,我们试图求解形如  LVI(Ω,H,c)∶(x-x*)T(Hx*+c)≥0,(∨)x∈Ω  的一类线性变分不等式问题。首先我们通过一些例子分析了已知自适应算法[5,11]的不足之处,其次我
农村基层干部如何落实十六届四中全会提出的“加强党的执政能力建设”的号召?根据我在宁武县杨庄村担任党支部书记10多年的经历,我认为,必须按照胡锦涛总书记明确要求的“权
Hopfield神经网络是一个典型的递归神经网络,具有很强的联想记忆和优化计算能力,所以该模型的动态特性已受到广泛关注,并且Hopfield神经网络的研究和应用已经渗透于各项工作之中
本文将基于Hirosawa和Reissig的一篇论文(F.Hirosawa, M.Reissig, Levicondition for hyperbolic equations with oscillating coefficients, J.D.E.223(2006)329-350),主要关