直接间断galerkin有限元办法求解反应扩散方程组.pdf

2015年 3月 高 等 学 校 计 算 数 学 学 报 第 37卷第 1期 直接间断Galerkin有限元方法 求解反应扩散方程组 张荣培 f辽宁石油化工大学理学院，抚顺 113oo1) THE DIRECT DISCONTINUOUSGALERKIN FINITE ELEM ENT M ETHOD FOR REACTION DIFFUSION EQUATIONS ZhangRongpei (SchO0lofScienceLiaoningShihuaUniversity，Fushun11300) Abstract ThispaperdiscussesthedirectdiscontinuousGalerkin(DDG)finite elementmethod fornonlinearlreactiondiffusion equations．Thebasefuncton n 矗niteelementspaceisconstructedbasedontheGaussintegrationpoint，theIl1lmer— icalfluxisobtainedbysolvingthediffusivegeneralizedRiemannproblem (dGRP)， sowepr0poseadirectdiscontinuousfiniteelementmethod．Theimplicitintegra_ tionfactorfIIF)timediscretizationisappliedforsolvingthenonlnearordnarY di髓rentialequations(ODEs)producedbyDDGspatialdiscretization·umenca exDerimentsarisingfrom theecologymodelarepresentedtoverifythee~ciency ofourmethod． KevWOrds DirectGalerkinmethod，reactiondiffusion equatiOils’mplctnte— grationfactormethod． AMS(2000)subjectclassifications 65M60 中图法分类号 O241．82 国家自然科学基金资助项目 收稿 日期：2012—01_03． 2 · 张荣培：直接间断 Galerkin有限元方法求解反应扩散方程组 第 1期 1 引 言 化学反应流、燃烧以及生态学上的食物链种群演化等模型，在数学上都可以用如下的 反应扩散方程 (组)来描述 OU = DAU+f(u)，(，￡)∈2【×(0， ， (1．1) 其中U—fu1(X，t)，u2(x，t)，… ，札M(x，t)]为各物种密度，M 为参加反应的物种数，D为 关于扩散系数的M ×M 对角矩阵，△是Laplacian算子，f是关于u 的非线性函数．t为 时间， 为一维或二维有界区域．近年来，许多学者对方程 (1．1)提出了多种数值方法，主 要有有限差分方法，有限元方法等 f1-3]． 间断Galerkin(DG)有限元方法是近年来比较流行的数值方法，该方法最初用于求解 双 曲方程 f4]．DG方法通过在原方程两边同时乘以试探函数并在每个单元积分，得到原问 题的弱形式，然后选取适当的试探函数空间，并在单元边界引入数值通量形成连接条件． 选择合适的数值流量是 DG方法的关键．在 DG求解双 曲方程时，可以通过求解黎曼问 题构造各种数值流量 对=F扩散方程的弱形式，定义适当的数值流