论文部分内容阅读
本文用谱方法讨论如下形式的具常迁移率的四阶抛物型CahnHilliard方程的初边值问题：
{(e)u/(e)t+γ(e)4u/(e)x4=(e)2ψ(u)/(e)x20＜t＜T,-π≤x≤π(1){u(x-π，t)=u(x+π,t)0＜t＜T(2){u｜t=0=u0(x)-π≤x≤π(3)
其中u=u(x，t)为相互扩散的两种物质之一的密度，γ＞0为迁移率，假设其为常数。ψ(u)≡H(u)，H(u)为双井位势的典型形式：
H(u)=-1/2u2+1/3γ1u3+1/4γ2u4
本文只讨论γ1=0的情形，此时ψ(u)≡H(u)=-u+γ2u3本文的主要目的是要在γ2＞0时用谱方法建立问题(1)，(2)，(3)半离散近似并讨论其稳定性和收敛性。
Cahn-Hilliard方程是一类重要的四阶非线性扩散方程。最初是由Cahn和Hilliard于1958年在研究热力学中两种物质(如合金、聚合物等等)之间相互扩散现象时提出。后来在描述生物种群的竞争与排斥现象，河床迁移过程，固体表面上微滴的扩散等许多扩散现象的研究中也提出了同样的数学模型。从本世纪八十年代以后，人们开始系统的研究Cahn-Hilliard方程，近年来由于Cahn-Hilliard方程在化学、化工和材料科学等方面的重要背景，关于Cahn-Hilliard方程的理论研究成果甚多，本文主要从数值方面来研究Cahn-Hilliard方程。所谓谱方法就是以三角函数为基底的有限元法。三角函数具有很好的逼近性质，这使得谱方法成为一种高精度的数值方法，选取三角多项式作为逼近空间的基函数，用谱方法构找出相应的拟谱近似方程。
本文首先对Cahn-Hilliard方程之周期边值问题(1)，(2)，(3)建立了与之等价的变分形式。证明了如果u：[0，T]→C4[0,π]是方程(1)，(2)，(3)的解，则u一定是方程((e)u/(e)t,v)+γ((e)2u/(e)x2,(e)v/(e)x2)=(ψ(u),(e)2v/(e)x2)(A)v∈HE(4)的解.如果u：[0，T]→HE是方程(4)的解，对任意的t有u(t，·)∈C4(0，π)且u满足边界条件(2)及初始条件(3)则它一定是方程(1)，(2)，(3)的解.然后用谱方法对空间变量x作离散近似得到了以下半离散问题：求解方程(3.1)的谱方法为求uN：[0，T]→SN使之满足((e)uN/(e)t,v)N+γ((e)2uN/(e)x2,(e)2v/(e)x2)N=(ψ(uN),(e)2v/(e)x2)N(A)v∈SN(5)uN(x-π，t)=uN(x+π，t)(6)uN|t=0=uN0(7)其中uN=uN(x，t)=a0(t)+N∑j=1aj(t)cosjx+N∑j=1bj(t)sinjx,uN0=PNu0为u0(x)在SN上的投影。接下来证明了半离散问题(5)，(6)，(7)的解的存在唯一性。即证明了当γ＞0，γ2＞0时半离散问题(5)，(6)，(7)在[0，T]上必存在唯一的解。在此基础上证明了半离散问题(5)，(6)，(7)解的稳定性和收敛性。即证明了如果u∈HkE是方程(5)，(6)，(7)的解，则(A)t∈[0，T]有‖e‖+∫t0‖exx(τ)‖2Ndτ≤CN-k+‖e(0)‖其中C是仅依赖于γ，γ2和u0(x)的常数.‖e(0)‖是初始逼近误差。设u是方程(1)，(2)，(3)的解，uN是方程(5)，(6)，(7)的解，我们得到‖u-uN‖≤C‖u0-uN0‖+CN-k.