Allen—Cahn方程弱形式的勒让德配置法及基于张量积的矩阵形式.pdfVIP

2013年 9月 四川师范大学学报(自然科学版) Sept．，2013 第36卷 第5期 JournalofSicbuanNormalUniversity(NaturalScience) VD1．36．No．5 Allen—Cahn方程弱形式的勒让德配置法 及基于张量积的矩阵形式 吴夭云z顺，l贝 (贵州师范大学数学与汁算机科学学院，贵州 贵阳550001) 摘要：给出了Allen—Cahn方程弱形式的勒让德配置法的全离散格式，并对稳定性进行了分析．在适当 的稳定性条件下，证明了给出的格式是能量稳定的．然后利用勒让德微分矩阵和Gauss—Lobatto积分的性 质，全离散的勒让德配置格式被转化为矩阵方程 ，再利用张量积转化为线性代数方程组． 关键词：Allen—Cahn方程；弱形式；勒让德配置法；稳定性分析；张量积 中图分类号：0241．8 文献标志码：A 文章编号：1001—8395(2013)05一‘0691—05 doi：10．3969／j．issn．1001—8395．2013．05．009 在本文中考虑的是下述Allen—Cahn方程弱形 述实际工程问题时，会包含一些复杂的初值函数和 式的勒让德配置法及其稳定性分析 内能项，要找到它们的精确解是不容易的．因此，设 一 计一些有效和精确的数值方案，使之满足相应的离 △ + Ⅱ)-0’ 散能量律，即能量稳定是非常重要的． (，)，，)∈ ×(0，]， (1) 主要考虑能量函数F(2／)的导数 u)=F(M) Im =0， t∈(0， 满足以下条件：存在一个常数 ￡使得 maxl厂(“)I≤L· (4) J：0= 0， (，Y)∈。， 注意到，通过限制 F(2／)的增长，即当ll≥ 时，用 其中，=(一1，1)。，2／被定义为混合物中2种成分 2／的二次式来代替F()，这个条件被许多相关的 之间的浓度差 (2／)=F (2／)，F(“)是一个给定的 电位满足．众所周知，Allen—Cahn方程满足最大值 势能．这里小的正常数 是与交面能量相关的梯度 原理，可以截断F(3／,)，使之在区间[一 ， ]的外 能量系数u，n是外法向， 是总时间．而Allen— 部按二次增长，如果初始条件 ‰ 的最大范数不超 Cahn方程可以被视为李雅普诺夫能量函数E(2／) 过 ，则不会对解造成影响_2J． 在 ￡中的梯度流 ． 关于Allen—Cahn方程的数值分析已经有大量 E(M)：=I(÷lVMI+ F(“))dxdy，(2) 的工作 卜”J．这些工作都是对于有限 有限差分 Jtl ‘- ￡ 1 法或带着周期性边界条件的傅里叶谱方法进行分 更确切地说，利用 (1)式与 一△ + 厂(u)作 内积，