数值数学(海德堡)。2015; 131(1): 1–31.
广义二级流体的Rayleigh–Stokes问题分析
,1 ,2 ,1,三和三 埃米利亚·巴兹列科娃
1保加利亚科学院数学与信息研究所。保加利亚索非亚1113号Bl.8 G.Bonchev街
金邦迪(Bangti Jin)
2英国WC1E 6BT伦敦Gower街伦敦大学学院计算机科学系
雷奇·拉扎罗夫
1保加利亚科学院数学与信息研究所。保加利亚索菲亚1113号G·邦切夫街8号楼
三得克萨斯农工大学数学系,美国得克萨斯州大学城,邮编77843
Zhi Zhou公司
三德克萨斯农工大学数学系,美国德克萨斯州大学城,邮编77843
1保加利亚科学院数学与信息研究所。保加利亚索非亚1113号Bl.8 G.Bonchev街
2英国WC1E 6BT伦敦Gower街伦敦大学学院计算机科学系
三德克萨斯农工大学数学系,美国德克萨斯州大学城,邮编77843
通讯作者。 2014年4月15日收到;2014年10月29日修订
开放式访问本文是根据知识共享署名许可证的条款分发的,该许可证允许在任何媒体上使用、分发和复制,前提是原创作者和来源得到了认可。
摘要
我们研究了广义二级流体的Rayleigh–Stokes问题,其中包含Riemann–Liouville分数阶导数,并对连续、空间半离散和全离散形式的问题进行了分析。我们建立了光滑和非光滑初始数据齐次问题的Sobolev正则性,包括发展了一种使用连续分段线性有限元的空间半离散Galerkin格式,并导出了有限元近似的相对于初始数据正则性的最优误差估计。进一步,基于后向欧拉方法和二阶后向差分方法以及相关的卷积求积,发展了两种全离散格式,并导出了光滑和非光滑初始数据的全离散近似的最优误差估计。给出了具有光滑和非光滑初始数据的一维和二维算例的数值结果,以说明该方法的有效性,并验证了收敛理论。
数学学科分类:65M60、65M15
介绍
本文研究具有分数阶导数模型的广义二级流体的齐次Rayleigh–Stokes问题。让是边界为的凸多面体域、和成为一个固定的时间。然后通过以下公式给出数学模型
哪里是一个固定常数,是初始数据,、和是Riemann–Liouville分数阶导数由定义[11,24]:
瑞利-斯托克斯问题(1.1)近年来受到了相当大的关注。分数导数在模型中,用于捕捉流体的粘弹性行为;参见示例[5,28]获取派生详细信息。模型(1.1)在描述某些非牛顿流体的行为方面起着重要作用。
为了深入了解该模型解的行为,人们对推导特殊情况下的闭式解有着浓厚的兴趣;参见,例如[5,28,32]. 例如,Shen等人[28]利用傅里叶正弦变换和分数拉普拉斯变换得到了问题的精确解。赵和杨[32]在齐次初始和边界条件的情况下,使用矩形域上的本征函数展开导出精确解。在这些研究中获得的解本质上是形式的,尤其是解的正则性尚未得到研究。在Sect。 2下面,我们填补了这个空白,并建立了光滑和非光滑初始数据解的Sobolev正则性。我们想提一下Girault和Saadouni[7]分析了一类密切相关的含时二级流体模型弱解的存在唯一性。
这些研究中获得的精确解涉及无穷级数和特殊函数,例如广义Mittag–Leffler函数,因此不便于数值计算。此外,封闭式解决方案仅适用于有限类别的问题设置。因此,必须为问题开发高效且最精确的数值算法(1.1). 这是在早些时候考虑的[1,2,12,21,31]. Chen等人[1]基于空间有限差分法和时间分数阶导数的Grünwald–Letnikov离散化,发展了隐式和显式格式,并使用傅里叶方法分析了它们的稳定性和收敛速度。同样的味道是工作[2],其中考虑了基于傅里叶级数展开的方案。吴[31]通过变换问题发展了一种隐式数值逼近格式(1.1)转化为一个积分方程,并通过能量论证证明了其稳定性和收敛性。林和江[12]描述了一种基于再生核希尔伯特空间的方法。最近,Mohebbi等人[21]比较了紧致有限差分法和径向基函数法。然而,在所有这些研究中,误差估计是在以下假设下获得的:(1.1)足够平滑且域是一个矩形。因此,不包括非光滑数据(初始数据或右侧数据)和一般域的有趣情况。
在过去的十年中,关于分数阶导数微分方程数值方法的理论研究受到了相当大的关注。麦克林和穆斯塔法[18,22]对分段常数法和分段线性间断Galerkin法进行了时间分析,推导了光滑初始数据的误差估计;另请参见[23]获得相关的超收敛结果。在[8,10]问题的空间半离散Galerkin有限元法和集中质量法具有(与是一个椭圆算子,并且是Caputo衍生物)进行了分析。为初始数据建立了几乎最优的误差估计,(参见第2(定义见下文)。注意,这包括弱(非光滑),以及非常微弱的数据,.英寸[19,第4节],McLean和Thomée研究了以下方程(与为的Riemann-Liouville积分和导数算子和(分别是)和派生的-空间半离散格式的误差估计和(并适当平滑)并讨论了一些基于拉普拉斯变换的全离散格式。相应的数据估计和衍生于[20]. Lubich等人[15]为该问题开发了两个完全离散的方案具有和基于分数阶导数项的卷积求积,导出了非光滑初始数据和右侧的最优误差估计。Cuesta等人[4]考虑了具有卷积求积的模型的半线性对应项,它还涵盖了分数扩散情况,即。,,并提供了一个统一的框架,用于在抽象的Banach空间设置中进行具有最佳误差估计的误差分析。
本文针对这一问题开发了Galerkin有限元法(1.1)对于光滑和非光滑的初始数据,导出了关于数据正则性误差估计的最优解。近似值基于有限元空间形状正则拟均匀划分族上的连续分段线性函数域的进入之内-单纯形,其中是最大直径。问题的半离散Galerkin有限元法(1.1)是:查找这样的话
哪里、和是初始数据的近似值。我们的默认选项是投影,假设和Ritz投影,假设此外,我们基于后向Euler方法和二阶后向差分方法以及分数阶导数项的相关卷积求积,开发了两个全离散格式,分别在时间上达到了一阶和二阶精度。对于半离散和全离散格式,都提供了关于数据正则性的最佳误差估计。
我们的主要贡献如下。首先,在定理中2.1,使用运算符方法[25]通过建立问题解的平滑性和衰减性,我们为我们的研究奠定了理论基础(1.1). 其次,对于平滑的初始数据和非光滑初始数据,我们导出了空间半离散格式的误差估计,cf.定理3.1和3.2:
对恶化为方法。误差估计是根据Fujita和Suzuki的方法得出的[6]. 接下来,在定理中4.1和4.2我们建立最佳两种完全离散格式的误差估计。该证明受到了Cuesta等人的基本工作的启发[4]它依赖于卷积求积的已知误差估计和卷积核的界。例如,我们证明了离散解通过反向欧拉方法(在时间步长为1的均匀网格上)满足以下先验误差界
二阶后向差分法也有类似的估计。
论文的其余部分组织如下。在Sect。 2我们建立了解的Sobolev正则性。在Sect。 三,我们分析了空间半离散格式,并导出了光滑和非光滑初始数据的最优误差估计。然后在帮派。 4基于分数阶导数的卷积求积逼近,我们发展了两种全离散格式。两种方案都提供了最佳误差估计。最终进入门派。 5给出了一维和二维算例的数值结果,以说明收敛理论。从头到尾,符号表示在不同情况下可能不同的常数,但它始终独立于解,网格大小和时间步长.
溶液的规律性
在本节中,我们建立了(1.1)在同质情况下.我们首先回顾关于椭圆算子和函数空间的预备知识。然后我们导出了适当的解表示,证明了弱解的存在性,并建立了齐次问题解的Sobolev正则性。主要工具是在[25]. 此外,我们通过特征函数展开给出了另一种解表示,并导出了时间相关分量的定性性质。
前期工作
首先我们介绍一些符号。对于,我们表示为范数诱导的希尔伯特空间
具有表示中的内积和和分别是的Dirichlet特征值和特征函数在域上。像往常一样,我们识别一个函数在里面具有功能在里面由定义,对于所有人。然后设置和形式正交基和分别是。因此是标准的是标准的和相当于中的规范什么时候在[29]. 请注意形成插值空间的希尔伯特尺度。因此,我们表示作为插值尺度上的范数和对于在间隔中和作为插值尺度上的范数和什么时候在中然后和规范与任何通过插值,同样地和规范与任何.
对于和我们引入轮廓由定义
其中圆弧的方向是逆时针的,而两条射线的方向是增加虚部。此外,我们表示为行业
我们重铸问题(1.1)带有通过将控制方程的两边积分为Volterra积分方程(1.1)
其中内核由提供
和操作员由定义使用域. The椭圆问题的正则性对于我们的讨论至关重要,它来自于域上的凸性假设众所周知,操作员生成角的有界解析半群,即,对于任何
同时,将拉普拉斯变换应用于(2.1)收益率
即。,,带内核由提供
哪里是函数的拉普拉斯变换因此,通过拉普拉斯逆变换,我们推断出解算子由提供
哪里.
首先我们陈述一个关于核的基本估计.
引理2.1
修复,并让在中定义(2.3). 然后
证明
让,即。然后通过注释,
为了证明(2.5)我们注意到
让.由于功能达到最低值,具有最小值,它来自(2.7)那个
自,这导致了第一个断言
发件人(2.7)由此可见
因此,我们得到
这就完成了引理的证明。
解的先验估计
现在我们可以说明问题的规律性(1.1)带有.
定理2.1
对于任何和有一个独特的解决方案到问题(1.1)和
此外,以下稳定性估计适用于和:
哪里是常数取决于和,和常数也取决于.
证明
按引理2.1和(2.2)我们获得
我们根据(2.3)和(2.11)那个
然后由[25,定理2.1和推论2.4],对于任何有一个独特的解决方案第页,共页(2.1)它是由
仍需显示估计值。
让。我们选择并简称为
首先我们推导(2.9)的和.来自(2.4)和(2.12)我们推断
接下来我们证明估计(2.9)的和.通过应用运算符到的两侧(2.4)以及我们达到的差异化
使用身份
它是从(2.12)和引理2.1那个
通过采取,我们从(2.14)
这显示了估计值(2.9). 证明估算(2.10)带有我们观察到
现在通过记录身份
事实上对于,我们有
由(2.11)我们获得
因此,使用这个估计,我们得到
自对于,我们推断
具有最后,请注意(2.10)带有等于(2.9)带有和替换为.
备注2.1
我们注意到这个论点适用于任何扇形算子,包括空间中的Riemann–Liouville分数导数算子[9].
此外,定理中的估计2.1通过插值暗示以下结果。
备注2.2
解决方案到问题(1.1)带有满足
在哪里和或和.
关于解决方案行为的进一步讨论
估计值(2.9)保留任何然而,在这种情况下和我们可以大大改进这个估计也就是说,如果我们应用边界来自(2.15)在估计(2.14),我们得到以下更清晰的大边界:
备注2.3
对于我们有以下界限
它比(2.9)对于大型.此与绑定在一起(2.9)带有,暗示以下问题解决的先验估计(1.1):
此外,通过应用特征函数展开,求解了瑞利-斯托克斯问题(1.1)可以写在表格中
哪里和满足以下方程式:
要解决(2.17)我们应用拉普拉斯变换并使用恒等式
适用于功能,连续,并且这样是有限的[16,方程式(1.15)]。这样,对于拉普拉斯变换,一个到达
基于这种表示,在下一个定理中,我们总结了时间相关分量的一些性质,这对于研究解的行为,包括非齐次问题是有用的。
回想一下,一个函数被称为完全单调的当且仅当
定理2.2
功能具有以下属性:
其中常量不依赖于和.
证明
我们介绍辅助功能由它们的拉普拉斯变换定义
通过拉普拉斯变换的性质我们获得和进一步,采用拉普拉斯逆变换(2.17),我们得到
哪里是Bromwich路[30]. 积分下的函数有一个分支点,所以我们切断了实轴的负部分。注意,函数在Riemann曲面的主片中(包括其切割边界)没有零。的确,如果,使用,然后
自从和有相同的标志和因此,可以通过将Bromwich路径弯曲到Hankel路径来找到,从开始沿着负实轴的下侧,围绕圆盘逆时针方向,结束于沿着负实轴的上侧。通过采取我们获得
哪里
自、和,有个保持为所有人因此,根据伯恩斯坦定理,是完全单调的函数。特别是,它们是正的且单调递减的。这显示了前两个断言。
以同样的方式,我们证明了是完全单调的,因此.签署人(2.18),以及(2.20),
在进行拉普拉斯逆变换时,意味着现在第三个断言是
最后,使用表示法
具有
其中函数定义为(2.3),最后一个断言是通过应用定理证明中的论点得出的2.1具有替换为并使用以下类似于(2.15):
这就完成了命题的证明。
根据定理2.1,对于任何,解决方案运算符在二阶空间中具有平滑特性。在极限情况下然而,它没有任何平滑特性。为此,我们考虑本征函数展开:
在这种情况下我们从中推断(2.17)和(2.18)
这表明该问题不具有平滑特性。
备注2.4
我们观察到,如果,然后行为类似作为此行为与细分扩散方程解的行为相同;参见[17,定理4.1]和[26,定理2.1]。然而,作为腐烂得像,如标准扩散方程的情况。解决方案第页,共页(1.1)腐烂得像对于。这比,细分扩散方程解的衰减[26,推论2.6],但比扩散方程的指数衰减慢得多。
我们可以推广定理2.1对于非常弱的初始数据,即。,具有显然,对于任何人函数满足等式(1.1)在某种意义上.然后我们呼吁扩大(2.21). 重复定理的论证2.1收益率根据勒贝格的支配收敛定理,我们推导出
因此,函数满足(1.1)和用于收敛到在里面即。,表示解决方案。此外,定理的论证2.1收益率对于任何.
半离散Galerkin有限元法
在本节中,我们考虑空间半离散有限元近似,并导出齐次问题的最佳误差估计。
半离散Galerkin格式
首先我们回顾一下-正交投影和里兹投影分别由定义
对于对于,的-投影定义不明确。尽管如此,人们可能会认为对于作为空间之间的对偶和并定义以同样的方式。
里兹投影和-投影具有以下属性。
引理3.1
让网格是准均匀的。然后是操作员和满足:
此外,在上稳定对于.
引入离散拉普拉斯算子由定义
和,我们可以写下空间离散问题(1.2)至于发现这样的话
哪里是初始条件的适当近似值相应地,解决方案操作员对于半离散问题(1.2)由提供
哪里是中定义的轮廓(2.13)和此外,利用本征对离散拉普拉斯算子,我们定义了离散范数在空间上对于任何
操作员的稳定性如下所示。这个证明与定理的证明相似2.1,因此省略。
引理3.2
让由定义(3.2)和.然后
在哪里和或和.
现在我们导出了半离散Galerkin格式的误差估计(3.2)使用操作员技巧,遵循Fujita和Suzuki的有趣工作[6]. 我们注意到,类似的估计也来自于[10],但以牺牲额外的对数因子为代价在非光滑初始数据的情况下。
以下引理在推导误差估计中起着关键作用。
引理3.3
对于任何和对于,有个保持
证明
由[6,引理7.1],我们对任何
或者,它遵循不等式
带着选择和.通过引理2.1,为所有人,这就完成了证明。
下一个引理显示了之间的误差估计及其离散模拟.
引理3.4
让、和.然后是保持
证明
根据定义,和分别满足
减去这两个恒等式得到以下误差的正交性关系:
这个和引理3.3暗示对任何人
通过采取,的拉格朗日插值利用柯西-施瓦兹不等式,我们得出
再次吸引引理3.3选择,我们获得
因此
鉴于(3.8),一个绑定可以导出
它由此而来(3.7)那个
这就产生了
这给出了所需的界限下一步,我们推导出通过二元论证。对于,通过设置
我们有二元性
然后,期望的估计值如下所示(3.6)和(3.9)由
这就完成了引理的证明。
半离散格式的误差估计
现在我们可以说明非光滑初始数据的误差估计.
定理3.1
让和是问题的解决方案(1.1)和(3.2)带有和分别是。然后针对,持有:
证明
错误可以表示为
具有和.通过引理3.4定理证明中的论证2.1我们有
类似的论点也产生了-估计。
接下来我们讨论平滑初始数据的情况,即。,和.我们再次测量轮廓。然后是错误可以表示为
通过平等
我们可以获得
哪里和然后我们得出以下误差估计。
定理3.2
让和成为问题的解决方案(1.1)和(3.2)带有和分别是。然后针对,持有:
证明
让和.然后是引理3.1和3.4和身份给
现在它是从这个和表示(3.10)那个
因此,我们获得-误差估计。这个-误差估计如下。
备注3.1
对于平滑的初始数据,我们也可以取近似值。然后可以将错误拆分为
定理3.2给出了第一项的估计。第二项的界来自引理3.1和3.2
因此,误差估计(3.11)保持初始近似值它由此而来,定理3.1和插值和,有个保持
备注3.2
如果初始数据很弱,即。,具有,然后是[8,定理2]给出了半离散有限元近似的以下最佳误差估计(1.2)
完全离散方案
现在我们为这个问题开发了两个完全离散的方案(1.1)基于卷积求积(参见[4,13–15]以供详细讨论),并导出平滑和非平滑初始数据的最佳误差估计。
卷积求积
首先,我们简要描述了中的抽象框架[4第2节和第3节],该节有助于开发和分析全离散方案。让是在扇区内解析的复值或算子值函数并以
对于一些实数和.然后是分布的拉普拉斯变换在实线上消失,其单数支持为空或集中于,是的分析函数。对于,解析函数由反演公式给出
哪里是位于分析性扇区的轮廓,平行于其边界并以增加的虚部为方向。使用作为时间差,我们定义作为与核的(分布式)卷积的算子对于函数具有适当的平滑度。
卷积求积近似通过离散卷积具体来说,我们划分时间间隔进入之内具有时间步长的等子间隔,并定义近似值:
其中正交权重由生成函数决定
在这里是稳定一致线性多步方法生成多项式的商。在这项工作中,我们考虑了反向欧拉(BE)方法和二阶反向差分(SBD)方法,其中
现在我们专门研究半离散问题的构造(3.2). 通过集成(3.2)来自到,我们得到了半离散解的表示
哪里,表示Riemann–Liouville积分左边是一个卷积,我们近似于具有通过
其中符号和参考由相应的线性多步方法生成的相关卷积求积。为了便于数值实现,我们以时间步进的形式重写它们。
反向欧拉(BE)方法
BE方法如下:Find对于这样的话
用卷积求积和由BE方法生成。通过应用加入计划(4.2)以及卷积的结合性,我们推断它可以改写为:和,查找对于这样的话
备注4.1
在方案中(4.3),术语位于在里面可以省略而不影响其收敛速度[15,27].
二阶后向差分(SBD)方法
现在我们转向SBD方案。众所周知,如果例如,对于[13,定理5.1][4,第3节]。在数值上也观察到问题的一阶收敛性(1.1). 因此,需要纠正方案,我们遵循[4,15]. 使用身份
我们可以重写半离散解进入之内
哪里和这导致了卷积求积
保留操作员的目的在中完好无损(4.4)是为了达到二阶精度,参见引理4.4如下所示。出租并注意到身份在栅格点和卷积的结合性(4.4)可以重写为
接下来通过应用运算符,我们获得
因此,我们得出了一个时间步进方案:,查找这样的话
和用于
其中卷积求积由提供
用砝码由SBD方法生成。
全离散格式的误差分析(4.3)和(4.5)为了这个案子将按照[4,第4节]。
反向欧拉方法的误差分析
调用函数时来自(2.3)并表示
我们可以写出和作为
为了进行误差分析,我们需要以下估计[13,定理5.2]。
引理4.1
让分析和(4.1)保持。然后针对,基于BE的卷积求积满足
现在我们可以说明非光滑初始数据的误差估计.
引理4.2
让和成为问题的解决方案(3.2)和(4.3)带有和分别是。然后就有了
证明
由(2.2)和身份对于,有个保持
然后(4.7)和引理4.1(与和)给予
所需结果直接来自的稳定性.
接下来我们转向平滑初始数据,即。,.
引理4.3
让和是问题的解决方案(3.2)和(4.3)带有和分别是。然后就有了
证明
以身份
并表示,错误可以表示为
发件人(2.2)和引理2.1我们推断
现在引理4.1(与和)给予
所需的估计值直接来自等式.
备注4.2
按引理4.3,误差估计表现出奇异的顺序行为作为,即使对于平滑的初始数据尽管如此,作为,问题(1.1)简化为标准抛物线方程,因此光滑数据的奇异性消失,这与抛物线方程相一致[29].
现在我们可以说明全离散格式的误差估计(4.3)在光滑和非光滑初始数据下,利用三角形不等式,定理3.1和3.2,引理4.2和4.3分别针对非光滑和光滑的初始数据。
定理4.1
让和是问题的解决方案(1.1)和(4.3)带有和分别是。然后,以下估计成立。
如果和,然后
如果和,然后
备注4.3
对于,我们也可以选择.让是全离散格式的对应解通过该方案的稳定性,引理的一个直接结果4.3,我们有
因此定理中的估计4.1(a) 仍然保持然后通过插值估计,我们推断
备注4.4
如果初始数据很弱,即。,具有,通过引理4.2,逆不等式[三第140页]和引理3.1我们有
这个和备注3.2得出以下误差估计
二阶后向差分法的误差分析
使用,我们有
类似引理4.1,以下估计成立(参见[13,定理5.2][14,定理2.2])。
引理4.4
让在中进行分析和(4.1)保持。然后针对,基于SBD的卷积求积满足
现在我们可以说明非光滑初始数据的误差估计.
引理4.5
让和是问题的解决方案(3.2)和(4.5)带有和分别是。然后就有了
证明
由(2.2)和身份
有个等待
然后(4.8)和引理4.4(与和)给予
所需结果直接来自的稳定性.
接下来,我们转向平滑初始数据.
引理4.6
让和成为问题的解决方案(3.2)和(4.5)带有和分别是。然后就有了
证明
通过设置可以表示为
发件人(2.2)和引理2.1我们推断
现在引理4.4(与和)给予
所需的估计值来自等式.
然后我们对全离散格式进行以下误差估计(4.5).
定理4.2
让和成为问题的解决方案(1.1)和(4.5)带有和分别是。然后,以下误差估计成立。
如果、和,有个保持
如果、和,有个保持
备注4.5
通过格式的稳定性,引理的一个直接推论4.6和备注中的论点4.3,定理中的估计4.2(a) 仍然保持然后通过插值我们得到
备注4.6
如果初始数据很弱,备注中的论点4.4收益率
数值结果
在这一部分中,我们给出了数值结果来验证Sects中的收敛理论。 三和4我们将考虑具有光滑、非光滑和非常弱初始数据的一维和二维示例。在一维情况下,我们取,在二维情况下。这里我们使用符号对于集合的特征函数考虑以下四种情况。
平滑:哪个在.
非光滑:; 跳跃和导致; 但对于任何.
非常微弱的数据:哪一个是狄拉克-功能集中于.根据索波列夫嵌入定理,对于.
二维示例:哪个在对于任何.
在我们的实验中,我们固定了参数英寸(1.1)适用于所有情况。我们分别考察了对于非光滑初始数据的情况,我们特别关注以下方面的误差接近于零。这些例子的精确解可以用广义Mittag–Leffler函数表示,但这很难计算,因此我们在非常精细的网格上计算参考解。我们报告标准化误差和对于平滑和非平滑数据。
在我们的计算中,我们划分单位间隔进入之内等距子间隔,具有网格大小.有限元空间由连续的分段线性函数组成。类似地,我们采用具有时间步长的统一时间网格,使用是感兴趣的时间。
数值结果,例如(a)
首先,我们确定网格大小在因此,空间离散化产生的误差可以忽略不计,这使我们能够检查时间收敛速度。在表中,我们显示-误差范数对于不同的值。在表中,BE和SBD分别表示反向欧拉方法和二阶反向差分方法,速率是指当时间步长为(或网格大小)两半,括号中的数字表示理论收敛速度。在图中我们为以下对象绘制结果对数刻度。阶收敛速度和对于BE方法和SBD方法,分别观察到了,这与我们的收敛理论非常吻合。此外,我们观察到误差随着分数阶的增加而减小增加。
表1
这个-误差标准,例如(a):和
|
|
|
|
|
|
| 费率 |
---|
比利时 |
| 6.75电子三 | 2.42秒三 | 1.00秒三 | 4.55秒4 | 2.15秒4 |
1.15 (1.00) |
| 3.68秒三 | 1.73当量三 | 8.42当量4 | 4.13秒4 | 2.03秒4 |
1.04 (1.00) |
| 4.12秒4 | 2.03秒4 | 1.00秒4 | 4.96秒5 | 2.43秒5 |
1.03 (1.00) |
SBD公司 |
| 5.59秒三 | 4.82秒4 | 1.18秒4 | 2.77秒5 | 6.66秒6 |
2.06 (2.00) |
| 1.05秒三 | 2.39秒4 | 5.33秒5 | 1.28秒5 | 3.14秒6 |
2.08 (2.00) |
| 7.62秒5 | 1.64秒5 | 3.86秒6 | 9.48秒7 | 2.46秒7 |
2.06 (2.00) |
误差图,例如(a),使用和
在表中和图,我们显示-和-误差标准BE方案。我们设置并检查空间收敛速度。数值结果表明和分别针对-和-误差范数,充分证明了定理3.2此外,经验收敛速度几乎与分数阶无关.
表2
错误,例如(a):和
|
|
|
|
|
|
| 费率 |
---|
|
-规范 | 6月16日4 | 1.59当量4 | 4.00秒5 | 9.90秒6 | 2.38秒6 |
2.01 (2.00) |
-规范 | 1.19秒2 | 5.99秒三 | 2.99秒三 | 1.49秒三 | 7.26秒4 |
1.01 (1.00) |
|
-规范 | 1.58秒三 | 4.00秒4 | 1.00秒4 | 2.48秒5 | 5.95秒6 |
2.01 (2.00) |
-规范 | 3.92秒2 | 1.98秒2 | 9.88秒三 | 4.91秒三 | 2.40秒三 |
1.01 (1.00) |
|
-规范 | 1.38秒三 | 3.47秒4 | 8.67欧元5 | 2.15电子5 | 5.16页6 |
2.01 (2.00) |
-规范 | 3.56秒2 | 1.79秒2 | 8.96秒三 | 4.45电子三 | 2.17电子三 |
1.01 (1.00) |
错误,例如(a):和
数值结果,例如(b)
在表格中和我们给出了示例(b)的结果。时间收敛速度为和对于BE和SBD方法,分别参见表,并且空间收敛速度是有序的在里面-规范和在里面-标准,参见表对于非光滑的初始数据,我们对以下方面的误差特别感兴趣接近于零。因此,我们也在和在表中.数值结果完全证实了预测的速率。
表3
这个-误差范数,例如(b)at,使用
|
|
|
|
|
|
| 费率 |
---|
比利时 |
| 2.82秒2 | 1.42秒2 | 7.13秒三 | 3.56秒三 | 1.76秒三 |
1.00 (1.00) |
| 8.67欧元三 | 4.18秒三 | 2.05秒三 | 1.01秒三 | 4.97秒4 |
1.02 (1.00) |
| 9.06秒4 | 4.47秒4 | 2.21当量4 | 1.09秒4 | 5.42秒5 |
1.02 (1.00) |
SBD公司 |
| 7.14秒三 | 1.61秒三 | 3.92秒4 | 9.63电子5 | 2.38电子5 |
2.05 (2.00) |
| 2.46秒三 | 5.05秒4 | 1.17秒4 | 2.82秒5 | 6.91欧元6 |
2.06 (2.00) |
| 1.67秒4 | 3.58秒5 | 8.40秒6 | 2.04秒6 | 5.11e条7 |
2.08 (2.00) |
表4
错误,例如(b):和
|
|
|
|
|
|
| 费率 |
---|
|
-规范 | 1.63秒三 | 4.09秒4 | 1.02秒4 | 2.55秒5 | 6.30秒6 |
2.00 (2.00) |
-规范 | 4.04秒2 | 2.02秒2 | 1.01秒2 | 5.04秒三 | 2.51秒三 |
1.00 (1.00) |
|
-规范 | 5.87电子三 | 1.47秒三 | 3.66秒4 | 9.13秒5 | 2.26秒5 |
2.00 (2.00) |
-规范 | 1.62当量1 | 8.08电子2 | 4.04秒2 | 2.02秒2 | 1.00秒2 |
1.00 (1.00) |
|
-规范 | 1.47秒2 | 3.66秒三 | 9.15秒4 | 2.28秒4 | 5.65秒5 |
2.00 (2.00) |
-规范 | 4.48秒1 | 2.24秒1 | 1.12秒1 | 5.60秒2 | 2.78秒2 |
1.00 (1.00) |
此外,在表中和图我们展示了-示例(a)和(b)的误差范数和。为了检查空间离散化误差,我们固定了时间步长在并使用SBD方法,使得时间离散化误差可以忽略不计。我们观察到,在光滑情况下,即示例(a),空间误差基本上保持不变,而在非光滑情况下(即示例(b)),其恶化为在示例(b)中,初始数据对于任何和备注4.5,错误增长如下作为表中的经验比率和图与理论预测吻合良好,即。,对于.
表5
这个-示例(a)和(b)的误差范数、和
| 1个三 | 1个4 | 1个5 | 1个6 | 1个7 | 1个8 | 费率 |
---|
(a) | 2.48秒4 | 3.07秒4 | 3.27秒4 | 3.46秒4 | 3.55电子4 | 3.58电子4 |
0.02 (0) |
(b) | 2.28秒4 | 5.07秒4 | 1.22秒三 | 2.89秒三 | 6.78秒三 | 1.56秒2 |
0.37 (0.37) |
示例(a)和(b)的误差图对于
数值结果,例如(c)
根据备注,如果数据很弱4.4和4.6,我们只能期望小时间步长的空间收敛表中的结果用速率表示超收敛现象在中-规范和在中-规范。这是由于在一维中Dirac的解-函数作为初始数据从支撑点和有限元空间的两侧平滑具有良好的逼近性。当奇点未与网格对齐,表显示了一个和价格-和-误差的范数。
表6
错误,例如(c):、和
|
|
|
|
|
|
| 费率 |
---|
|
-规范 | 1.19秒4 | 2.98秒5 | 7.45秒6 | 1.86当量6 | 4.62秒7 |
2.00 (1.50) |
-规范 | 5.35秒三 | 2.69秒三 | 1.35秒三 | 6.72当量4 | 第3.34页4 |
1.00 (0.50) |
|
-规范 | 2.41秒三 | 6.04秒4 | 1.51秒4 | 3.77秒5 | 9.31秒6 |
2.00 (1.50) |
-规范 | 3.98秒2 | 1.99秒2 | 9.92秒三 | 4.95秒三 | 2.46秒三 |
1.00 (0.50) |
|
-规范 | 1.25秒2 | 3.12秒三 | 7.80秒4 | 1.94秒4 | 4.83秒5 |
2.00 (1.50) |
-规范 | 5.00秒1 | 2.50秒1 | 1.25秒1 | 6.23秒2 | 3.09秒2 |
1.00 (0.50) |
表7
错误,例如(c):和
|
|
|
|
|
|
| 费率 |
---|
|
-规范 | 5.84秒三 | 2.22秒三 | 8.15秒4 | 2.93秒4 | 1.04秒4 |
1.50(1.50) |
-规范 | 1.79秒1 | 1.29秒1 | 9.16秒2 | 6.44秒2 | 4.45秒2 |
0.52 (0.50) |
|
-规范 | 2.42秒2 | 9.54秒三 | 3.57秒三 | 1.30秒三 | 4.63秒4 |
1.48 (1.50) |
-规范 | 7.77秒1 | 5.68英1 | 4.07秒1 | 2.87秒1 | 1.98秒1 |
0.51 (0.50) |
|
-规范 | 8.01秒2 | 3.27秒2 | 1.25秒2 | 4.57秒三 | 1.64秒三 |
1.46 (1.50) |
-规范 | 2.65欧元 | 1.97e0 | 1.43e0 | 1.02欧元 | 7.05秒1 |
0.49 (0.50) |
数值结果,例如(d)
这里我们考虑单位正方形上的二维示例对于非光滑的初始数据。为了离散化这个问题,我们划分了单位区间进入之内具有网格大小的等距子间隔从而将域划分为小正方形。我们通过连接每个小正方形的对角线来获得域的对称三角剖分。表显示了BE和SBD方法的一阶和二阶时间收敛速度。空间误差和如表所示,这意味着收敛速度为在中-规范和在中-规范。在图中和我们绘制了表中所示的结果和分别是。所有数值结果都证实了我们的收敛理论。
示例(d)的误差图具有和
示例(d)的误差图:和在和
表8
这个-误差范数,例如(d)at,使用和
|
|
|
|
|
|
| 费率 |
---|
比利时 |
| 4.53秒三 | 2.15秒三 | 1.04秒三 | 5.17页4 | 2.56秒4 |
1.03 (1.00) |
SBD公司 |
| 1.33秒三 | 2.80秒4 | 6.48秒5 | 1.56秒5 | 3.79秒6 |
2.11 (2.00) |
表9
错误,例如(d):和
|
|
|
|
|
|
| 费率 |
---|
|
-规范 | 1.95秒三 | 5.02秒4 | 1.26秒4 | 3.12秒5 | 7.61秒6 |
2.01(2.00) |
-规范 | 3.29秒2 | 1.63秒2 | 8.11秒三 | 4.03秒三 | 1.97秒三 |
1.00 (1.00) |
|
-规范 | 7.79秒三 | 2.00秒三 | 5.03秒4 | 1.25秒4 | 2.98秒5 |
2.02 (2.00) |
-规范 | 1.43秒1 | 7.09秒2 | 3.53秒2 | 1.75秒2 | 8.56秒三 |
1.01 (1.00) |
|
-规范 | 1.97秒2 | 5.09秒三 | 1.28秒三 | 3.19秒4 | 7.05秒5 |
2.00 (2.00) |
-规范 | 4.44秒1 | 2.22当量1 | 1.11秒1 | 5.52秒2 | 2.69秒2 |
1.01 (1.00) |
结束语
在这项工作中,我们研究了二级广义流中Rayleigh–Stokes方程的齐次问题。利用算子理论方法建立了解的Sobolev正则性。提出了一种基于Galerkin有限元方法的空间半离散格式,以及两种基于后向Euler方法和二阶后向差分方法以及相关卷积求积的完全离散格式,并给出了半离散和完全离散的数据正则性误差估计的最优解离散方案。大量的数值实验充分证明了我们收敛分析的敏锐性。
致谢
作者感谢克里斯蒂安·卢比奇教授对该论文早期版本的有益评论,这使Section中的表述得到了显著改进。 4以及一位匿名的裁判,为许多建设性意见提供支持。B.Jin的研究得到了NSF拨款DMS-1319052的支持,R.Lazarov的部分研究得到了国家科学基金拨款DMS-1016525的支持,还获得了阿卜杜拉国王科技大学(KAUST)颁发的编号为KUS-C1-016-04的奖助。
工具书类
1Chen C-M,Liu F,Anh V.带分数导数的加热广义二级流体Rayleigh–Stokes问题的数值分析。申请。数学。计算。2008;204(1):340–351. doi:10.1016/j.amc.2008.06.052。[交叉参考][谷歌学者] 2Chen C-M,Liu F,Anh V.具有分数导数的加热广义二级流体Stokes第一问题的Fourier方法和外推技术。J.计算。申请。数学。2009;223(2):777–789. doi:10.1016/j.cam.2008.03.001。[交叉参考][谷歌学者] 三。Ciarlet PG公司。椭圆问题的有限元方法。费城:SIAM;2002[谷歌学者] 4Cuesta E,Lubich C,Palencia C.分数阶扩散波方程的卷积求积时间离散化。数学。计算。2006;75(254):673–696. doi:10.1090/S0025-5718-06-01788-1。[交叉参考][谷歌学者] 5Fetecau C,Jamil M,FetecauC,Vieru D。广义Oldroyd-B流体中边的Rayleigh–Stokes问题。Z.安圭。数学。物理学。2009;60(5):921–933. doi:10.1007/s00033-008-8055-5。[交叉参考][谷歌学者] 6Fujita,H.,Suzuki,T.:进化问题。参见:《数值分析手册》,第二卷,Handb。数字。分析。,II、 第789–928页。荷兰北部,阿姆斯特丹(1991年)
7Girault V,Saadouni M.关于二维时变二级流体模型。计算。数学。申请。2007;53(3–4):347–360. doi:10.1016/j.camwa.2006.02.048。[交叉参考][谷歌学者] 8Jin,B.,Lazarov,R.,Pasciak,J.,Zhou,Z.:分数阶抛物型方程的Galerkin有限元法.LNCS 8236。在:过程。第五个Conf.Numer。分析。申请。(2012年6月15日至20日)。柏林施普林格出版社,第24-37页(2013年)
9Jin B,Lazarov R,Pasciak J,Zhou Z。空间分式抛物方程有限元方法的误差分析。SIAM J.数字。分析。2014;52(5):2272–2294. doi:10.1137/13093933X。[交叉参考][谷歌学者] 10Jin B,Lazarov R,Zhou Z。分数阶抛物型方程半离散有限元方法的误差估计。SIAM J.数字。分析。2013;51(1):445–466. doi:10.1137/120873984。[交叉参考][谷歌学者] 11Kilbas A、Srivastava H、Trujillo J。分数阶微分方程的理论与应用。阿姆斯特丹:爱思唯尔;2006[谷歌学者] 12Lin Y,Jiang W.带分数导数的加热广义二级流体Stokes第一问题的数值方法。数字。方法部分差异。埃克。2011;27(6):1599–1609. doi:10.1002/num.20598。[交叉参考][谷歌学者] 13Lubich C.卷积求积和离散运算演算。一、。数字。数学。1988;52(2):129–145. doi:10.1007/BF01398686。[交叉参考][谷歌学者] 14卢比奇C.卷积求积重温。比特。2004;44(3):503–514. doi:10.1023/B:BITN.0000046813.23911.2d。[交叉参考][谷歌学者] 15Lubich C,Sloan IH,Thomée V.具有正记忆项的演化方程近似的非光滑数据误差估计。数学。计算。1996;65(213):1–17. doi:10.1090/S0025-5718-96-00677-1。[交叉参考][谷歌学者] 16Mainardi F,Gorenflo R.放松过程中的时间分数导数:一项辅导调查。分形。计算应用程序。分析。2007;10(3):269–308. [谷歌学者] 17McLean W.时间分数阶扩散方程解的正则性。ANZIAM J。2010;52(2):123–138. doi:10.1017/S1446181111000617。[交叉参考][谷歌学者] 18McLean W,Mustapha K.子扩散方程的不连续Galerkin方法的收敛性分析。数字。阿尔戈。2009;52(1):69–88. doi:10.1007/s11075-008-9258-8。[交叉参考][谷歌学者] 19McLean W,Thomée V.通过拉普拉斯变换和分数阶发展方程求积对数值解进行最大范数误差分析。IMA J.数字。分析。2010;30(1):208–230. doi:10.1093/imanum/drp004。[交叉参考][谷歌学者] 20McLean W,Thomée V.分数阶演化方程的拉普拉斯变换数值解。J.整合。埃克。申请。2010;22(1):57–94. doi:10.1216/JIE-2010-22-1-57。[交叉参考][谷歌学者] 21Mohebbi A,Abbaszadeh M,Dehghan M。紧凑有限差分格式和RBF无网格方法,用于求解带有分数导数的加热广义二级流体的二维Rayleigh–Stokes问题。计算。方法应用。机械。工程师。2013;264:163–177. doi:10.1016/j.cma.2013.05.012。[交叉参考][谷歌学者] 22Mustapha K,McLean W.分数扩散方程的分段线性间断Galerkin方法。数字。阿尔戈。2011;56:159–184. doi:10.1007/s11075-010-9379-8。[交叉参考][谷歌学者] 23Mustapha K,McLean W.分数扩散和波动方程的间断Galerkin方法的超收敛性。SIAM J.数字。分析。2013;51(1):491–515. doi:10.1137/120880719。[交叉参考][谷歌学者] 24波德鲁布尼一世。分数微分方程。圣地亚哥:学术出版社;1999[谷歌学者] 25普吕斯·J。进化积分方程及其应用,数学专著,第87卷。巴塞尔:Birkhäuser Verlag;1993[谷歌学者] 26Sakamoto K,Yamamoto M.分数阶扩散波方程的初值/边值问题及其在一些反问题中的应用。数学杂志。分析。申请。2011;382(1):426–447. doi:10.1016/j.jmaa.2011.04.058。[交叉参考][谷歌学者] 27Sanz-Serna JM。偏积分微分方程的数值方法。SIAM J.数字。分析。1988;25(2):319–327. doi:10.1137/0725022。[交叉参考][谷歌学者] 28Shen F,Tan W,Zhao Y,Masuoka T。分数阶导数模型下加热广义二级流体的Rayleigh–Stokes问题。非线性分析。真实世界应用。2006;7(5):1072–1080. doi:10.1016/j.nonrwa.2005.09.007。[交叉参考][谷歌学者] 29托梅·五世。抛物型问题的Galerkin有限元方法,计算数学中的Springer级数,第25卷。2.柏林:施普林格;2006[谷歌学者] 30Widder DV(抖动DV)。拉普拉斯变换。普林斯顿数学系列,第6卷。普林斯顿:普林斯顿大学出版社;1941[谷歌学者] 31Wu C.具有分数导数的加热广义二阶流体的Stokes第一问题的数值解。申请。数字。数学。2009;59(10) :2571–2583。doi:10.1016/j.apnum.2009.05.009。[交叉参考][谷歌学者] 32赵C,杨C.矩形微通道中粘弹性流体电渗流动的精确解。申请。数学。计算。2009;211(2):502–509. doi:10.1016/j.amc.2009.01.068。[交叉参考][谷歌学者]