流固问题和麦克斯韦方程组的数值方法研究及超收剑分析

来源 :山东大学 | 被引量 : 0次 | 上传用户:icekingfly
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
基于多孔弹性问题和Darcy-Forchheimer问题的研究背景,以及麦克斯韦方程超收敛分析的研究背景,本文分别研究了多孔弹性问题的耦合弱有限元和混合元方法,全解耦的弱有限元方法和非一致网格上的交错有限差分方法,Darcy-Forchheimer问题的二重网格稳定化混合元方法和时谐麦克斯韦方程高阶矩形棱元的超收敛分析。具体内容为:  第一章我们给出了耦合弱有限元和混合元方法求解多孔弹性问题,用弱有限元逼近位移,用混合元逼近流体的压力和速度。首先,我们介绍了多孔弹性问题的数学模型及其混合元变分形式;之后,我们定义了弱有限元空间和混合元空间,引入了一些投影算子,给出了离散变分格式,并证明了格式的稳定性和存在唯一性。随后,我们得到了半离散和全离散格式下位移、速率和压力关于时间最大模的最优误差估计。误差估计不需要假设约束的单位储存系数c0大于0,并且误差分析中没有用到Gronwall不等式。最后,通过三个数值算例来验证理论分析的准确性和有效性。前两个算例表明在c0>0和c0=0两种情况下,我们都得到了位移、速率和压力关于时间最大模的最优误差估计。第三个算例是cantilever bracket问题,通过计算,我们发现耦合连续有限元和混合元方法得到的压力成振荡形态,而耦合弱有限元和混合元方法得到的压力光滑稳定,这验证了我们的方法能有效的消除由c0=0引起的压力振荡现象。  第二章首先,我们回顾了多孔弹性模型,并给出了变分形式。接着,我们引入了最低阶弱有限元空间,并介绍了向量函数和标量函数的弱梯度和弱散度的定义。然后,我们给出了多孔弹性问题的全离散的分离的弱有限元格式。格式分两步分别求解压力和位移,第一步,已知位移的前两层的值,利用流体方程求解压力,第二步,利用求得的压力和力学方程求解位移。这样,我们就把多孔弹性模型完全解耦了。通过格式可知,压力和位移的第一层的值未知,因此,我们又给出了初始层分离的求解格式,在格式中,我们去掉了散度项,因此需要假设初始时间步非常小,保证位移关于时间的变化非常小以致于可以忽略。随后,我们证明了格式的存在唯一性,得到了压力和位移的能量模关于时间和空间的最优误差估计。格式中的稳定化参数是准确估计的。最后,数值实验说明了当c0大于0和c0接近0时,压力和位移的能量模关于时间和空间都是1阶的,这表明数值实验结果与我们的理论分析相吻合。  第三章首先,我们引入了固体压力和流体速度,将多孔弹性模型等价的转化为多孔弹性问题的四场模型,接着,给出了块中心有限差分和MAC有限差分方法的一些基本符号,以及两个逼近空间和两个等价的逼近格式。格式中,位移和速度的x-分量在单元竖边中点处逼近,位移和速度的y-分量在单元的横边中点处逼近,流体压力和固体压力在单元中心处,这与文献[62]中的交错格式不同,我们的格式未知量数量更少。然后,我们介绍了离散LBB条件,得到了格式的稳定性。随后,通过引入位移的插值,证明了当c0>0时,位移和压力的离散H1模在非一致网格下具有一阶最优收敛性;在一致网格下,压力的离散H1模、位移的x-分量沿x-方向的差分算子和位移的y-分量沿y-方向的差分算子的离散L2模具有二阶超收敛性,位移的x-(y-)分量沿y-(x-)方向的差分算子在不包含边界项时具有二阶超收敛性,在包含边界项时为1.5阶收敛性。当c0≥0时,我们得到了压力的离散L2模在一致网格下具有二阶超收敛,在非一致网格下具有一阶最优收敛性。最重要的是,稳定性分析和误差估计关于Lamé常数λ∈(0,+∞)都是一致成立的。因此,我们的格式能有效的消除压力振荡和泊松闭锁现象。最后,通过数值实验验证了理论结果。  第四章我们考虑了Darcy-Forchheimer模型,利用格林公式,得到了混合变分形式,给出了速度u和压力p满足的正则性假设。然后,我们定义了P12-P1元的有限元空间。显然,同阶混合元空间不满足LBB稳定性条件,因此,我们引入了压力映射稳定项,得到了离散变分形式。我们证明了速度和压力的有限元空间满足一个弱的离散LBB条件,从而证明了离散变分形式解的存在唯一性。由于在速度u的值为0处,|u|的导数不存在,因此,我们用光滑函数√ε2+u2的一阶导数近似的代替|u|的一阶导数。利用牛顿修正法,我们给出了二重网格算法。第一步,在粗网格下,求解非线性问题,第二步,利用求得的粗网格上的解,在细网格上,求解线性化的问题,从而得到逼近解。为得到误差估计,我们引入了椭圆投影,得到了它的逼近结果,并证明了在粗网格上,速度和压力的L2误差。然后,我们证明了二重网格方法速度和压力的L2误差估计,得到了粗网格、细网格和参数∈只需满足关系H=O(h1/2)和ε=O(h)。最后,通过三个数值算例来验证二重网格方法的准确性和有效性,前两个算例说明了二重网格方法和稳定的混合元方法得到的误差几乎相同,但是二重网格方法的计算时间更少。第三个算例是注入-产出问题,通过它,也说明了二重网格方法能得到和稳定的混合元方法几乎相同的数值逼近解。  第五章我们首先回顾了时谐麦克斯韦方程,给出了各个记号的意义。介绍了Nédélec元插值算子和标准的L2投影算子,以及他们的性质,证明了任意阶棱元的超逼近结果。然后,我们分析了二阶Nédélec元的插值超收敛性质,得到了电场E、磁场H和电场的一阶导数和二阶导数的插值超收敛结果,利用四个高斯点,我们定义了标量和向量函数的离散l2范数,利用插值超收敛结果,我们证明了二阶Nédélec元电场E、磁场H和curl(E)具有三阶超收敛性,最后,非一致网格和各向异性网格上的数值结果验证了理论分析。数值结果还表明了电场的一阶导数和二阶导数具有二阶超收敛性。然后,通过相似的分析,我们得到了三阶Nédélec元的插值超收敛性质,利用九个高斯点,我们定义了标量和向量函数的离散l2范数,证明了三阶Nédélec元电场E、磁场H和curl(E)具有四阶超收敛性,数值结果与理论分析相吻合。
其他文献
哈密顿函数在平衡点处的二次部分定号时,哈密顿函数在平衡位置有局部极值,可得这个平衡点是李雅普诺夫稳定。哈密顿系统在平衡点处非谱稳定时,必然非李雅普诺夫稳定,但谱稳定不能
行为金融学的发展,为我们研究金融市场中一些无法用传统金融学解释的情况,以及金融市场中存在的异象,提出了新的思路。投资者情绪属于行为金融学范畴,它会导致投资者对股票价格的认知偏差,从而导致股票价格偏离其基本价值,即投资者情绪对股票价格决策及其在股票市场上的变化产生重大影响。因此,对投资者情绪的研究有助于投资者更好地理解自己的情绪,使投资者能够更理性地面对金融市场的波动,从而做出合理的投资决策。本文将
在此论文中,主要介绍图中一定条件的不交团及独立的圈在一些二分图中的相关结果.  本文令G代表一个图,它的顶点集和边集分别用V(G)和E(G)来表示.设v∈V(G),则点v在G中的度数表
随着分数阶微分的发展以及它所具备的优点,它逐渐被广泛应用到各个领域,且在图像处理领域也得到了发展并取得了一些成果。图像处理在信息表达、交流等方面占有很重要的地位。
随着科技的快速发展,图像越来越成为人类生活和工作中用于传递信息和进行交流的重要数据载体,图像的有效表示是进一步处理图像的重要基础。用各种观察系统取得的图像很多都是纹
在语言型多属性决策中,处理语言信息的方法主要有扩展原理法,符号法,语言变量直接运算法和二元语义方法.虽然二元语义的应用较为广泛,但是二元语义的研究理论还不是很成熟.二
图论是数学的一个分支.它以图为研究对象.通常通过研究一个图存在什么样的子图以及这个图具有什么样的性质来了解这个图的结构.本文主要研究图中点不交的星图,小阶完全图,圈以
RNA聚合酶的主要功能是利用DNA来制造RNA。在转录过程中,RNA聚合酶使用DNA作为模板并使用腺嘌呤脱氧核苷酸(A)和胸腺嘧啶脱氧核苷酸(T),胞嘧啶脱氧核苷酸(C)、尿嘧啶脱氧核苷酸
本文的主要内容是基于代数曲线理论来研究孤子方程的Riemann theta函数表示形式的拟周期解,分别研究了Hu族,耦合修正Korteweg-de Vries族,Vakhnenko方程和一族新的耦合非线性演
由于地域差异和时间变化的存在,产生了空间面板数据,这类数据能同时反映时间和空间两个方向上的变化规律,比仅考虑地域差异的截面数据和只关注时间变化的时间序列数据适用性