×

反应扩散方程组的全隐式时间步长格式和非线性求解器。 (英文) Zbl 1336.65167号

摘要:在本文中,我们提出了反应扩散方程组的鲁棒、高效和精确的全隐式时间步长格式和非线性求解器。反应扩散系统的应用在文献中非常丰富,从发育生物学中的模式形成建模到癌症研究、伤口愈合、组织和骨再生以及细胞运动。因此,建模师、分析师和生物学家能够准确有效地求解复杂平稳区域和曲面上的高度非线性抛物型偏微分方程组,有时甚至是连续演化区域和曲面,这一点至关重要。本文的主要贡献是通过将牛顿法和Picard迭代应用于反向欧拉、Crank-Nicolson(及其修正)和分数阶(theta)方法,研究了全隐式格式。我们的结果表明,分数步长(θ)方法与每个时间步长的单个牛顿迭代相结合,与完全自适应牛顿方法一样精确;两者都优于Picard迭代。特别是,这些结果有力地支持了这样一个观察,即一次牛顿迭代就足以产生与使用自适应牛顿方法获得的结果一样准确的结果。这在求解演化区域和曲面上高度复杂的非线性偏微分方程时尤其有利。为了验证我们的理论结果,我们在静止平面域和静止曲面体上进行了各种适当的数值实验。

MSC公司:

65M60毫米 偏微分方程初值和初边值问题的有限元、Rayleigh-Riz和Galerkin方法
35K51型 二阶抛物型方程组的初边值问题
35K57型 反应扩散方程
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 阿舍尔,U.M。;Ruuth,S.J。;Wetton,B.T.R.,含时偏微分方程的隐式显式方法,SIAM J.Numer。分析。,32,3797-823(1995年)·Zbl 0841.65081号
[2] Atkinson,K.,《数值分析导论》(1989),John Wiley&Sons·Zbl 0718.65001号
[3] 班杰斯,W。;哈特曼,R。;Kanschat,G.,交易。II-通用面向对象有限元库ACM Trans。数学。软质。,33, 4, 24/1-24/27 (2007) ·Zbl 1365.65248号
[4] 巴雷拉,R。;Elliott,C.M。;Madzvmou,A.,《进化生物表面图案形成的表面元素法》,J.Math。生物学,63,1095-1119(2011)·Zbl 1234.92007年
[5] Batchelor,G.K.,《流体动力学导论》(2000),剑桥大学出版社·Zbl 0958.76001号
[6] Britz,D。;Ø斯特比,O。;Strutwolf,J.,《曲柄尼科尔森误差振荡的阻尼》,计算。生物化学。,27, 3, 253-263 (2003) ·Zbl 1039.92047号
[7] 卡斯特,V。;杜洛斯,E。;Boissonade,J。;De Kepper,P.,持续图灵型平衡化学模式的实验证据,Phys。修订稿。,64, 3, 2953-2956 (1990)
[8] Dufiet,V。;Boissonade,J.,《传统和非传统图灵模式》,J.Chem。物理。,96, 1, 664-673 (1992)
[9] Elliott,C.M。;斯蒂纳,B。;Venkataraman,C.,用进化的表面nite元素模拟细胞运动和趋化性,J.R.Soc.Interface,9,76,3027-3044(2012)
[10] Gierer,A。;Meinhardt,H.,《生物模式形成理论》,Biol。赛博。,12, 1, 30-39 (1972) ·兹比尔1434.92013
[11] Glowinski,R.,《不可压缩粘性流的有限元方法》。不可压缩粘性流的有限元方法,《数值分析手册》,第9卷(2003),爱思唯尔科学有限公司·Zbl 1040.76001号
[12] 海尔,E。;诺塞特,S.P。;Wanner,G.,解常微分方程I(1987),Springer-Verlag·Zbl 0638.65058号
[13] Hughes,T.J.R.,《有限元法》(2003),多佛出版社
[14] 俄亥俄州拉基斯。;Madzvamuy,A。;Venkataraman,C.,发展域上rdes的有限元近似隐式显式时间步进,SIAM J.Numer。分析。,51, 4, 2309-2330 (2013) ·Zbl 1280.65111号
[16] Madzvamue,A.,移动有限元的时间步长方案,应用于固定域和增长域上的反应扩散系统,J.Compute。物理。,214, 239-263 (2006) ·Zbl 1089.65098号
[17] Madzvamuy,A。;Maini,P.K。;Wathen,A.J.,应用于模型生物模式生成器的移动网格有限元方法,J.Compute。物理。,190, 478-500 (2003) ·Zbl 1029.65113号
[18] Maini,P.K。;贝克,R.E。;Chong,C.M.,《图灵模型来自分子时代》(特邀视角),《科学》,3141397-1398(2006)
[19] Murray,J.D.,《数学生物学》,第2卷(2003年),斯普林格出版社·Zbl 1006.92002号
[20] 尼尔森,M.P。;Veltman,D.M。;van Haaster,P.J.M。;韦伯,S.D。;麦肯齐,J.A。;Insall,R.H.,《趋化性:一个基于反馈的计算模型稳健地预测真实细胞行为的多个方面》,PLoS Biol。,9,e1000618(2011)
[21] sterby,O.,《减少曲柄尼科尔森振荡的五种方法》,BIT Numer。数学。,43, 4, 811-822 (2003) ·Zbl 1038.65080号
[22] 欧阳(Q.Ouyang)。;Swinney,H.L.,《从均匀状态到六边形和条纹图灵图案的转变》,《自然》,352610-612(1991)
[23] 普里戈金,I。;Lefever,R.,耗散系统中的对称破缺不稳定性。二、 化学杂志。物理。,48, 4, 1695-1700 (1968)
[24] Quarteroni,A。;萨列里,F。;Sacco,R.,《数值数学》(2006),斯普林格出版社·Zbl 0913.65002号
[25] Ruuth,S.J.,模式形成中反应扩散问题的隐式-显式方法,J.Math。《生物学》,34,2,148-176(1995)·Zbl 0835.92006号
[26] Saad,Y.,大特征值问题的数值方法(1992),SIAM·Zbl 0991.65039号
[27] Schnakenberg,J.,《具有极限循环行为的简单化学反应系统》,J.Theor。《生物学》,81,3,389-400(1979)
[28] 生病,S。;Reinker,S。;Timmer,J。;Schlake,T.、Wnt和dkk通过反应扩散机制确定毛囊间距,《科学》,3141447-1450(2006)
[29] Smoller,J.,《冲击波和反应扩散方程》(1994),Kluwer学术出版社·Zbl 0807.35002号
[30] Solnica-Krezel,L.,《脊椎动物发育:驯服节波》,《现代生物学》。,13, 7-9 (2003)
[31] (Thomas,D.;Kernevez,J.,《固定化酶系统的分析和控制》(1976),爱思唯尔),113-150,(关节酶膜、运输、记忆和振荡现象一章)
[32] Thomée,V.,抛物问题的Galerkin有限元方法(2006),Springer·Zbl 1105.65102号
[33] 图灵,A.M.,《形态发生的化学基础》,菲洛斯。事务处理。R.Soc.伦敦。序列号。B、 生物。科学。,641, 33-72 (1952) ·Zbl 1403.92034号
[34] Venkataraman,C。;Sekimura,T。;加夫尼,E.A。;Maini,P.K。;Madzvamuy,A.,parr mark pat的模型。形态。在阿玛戈鳟鱼的早期发育时期,Phys。版本E,84,4,041923(2011)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。