对流扩散方程的局部间断有限元方法研究

来源 :哈尔滨工业大学 | 被引量 : 0次 | 上传用户:wangming2106
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
作为有限体积法和有限元方法的进一步推广,间断有限元方法是一类用于求解双曲方程的高效的数值解法。对于高阶方程,首先通过建立辅助变量将其转换为一阶方程组,然后对方程组应用间断有限元方法进行计算,我们称这样的方法为局部间断有限元方法。二者均具有高阶收敛精度,是高效稳定的近现代数值算法。我们可以通过选择间断有限元空间中用于逼近真解的多项式的阶数,来使得数值解对精确解的收敛误差具有任意高的阶数。此外,针对解存在间断的情况,可以通过选择恰当的斜率限制器让数值解在间断点附近能够得到更好的控制,即使其呈现出陡峭的状态,同时消除震荡。本文主要研究周期边界条件下的一维变系数线性对流扩散方程的间断有限元方法的稳定性和误差估计。稳定性可以保证由格式得到的数值解是可由其初值控制住的,不会随时间的发展而发生严重偏移;收敛性是衡量一个数值算法是否优秀的重要标准之一,好的算法应该具有较高的收敛精度和数值分辨率,是数值解与精确解之间相似程度的体现。本文给出了一维变系数线性对流扩散方程的局部间断有限元方法的稳定性结论和收敛性结论,对局部间断有限元方法用于求解变系数偏微分方程的理论依据进行了补充。首先,本文针对对流项带变系数函数,扩散项为常系数的对流扩散方程的稳定性以及对流项为常系数,扩散项带变系数函数的对流扩散方程的稳定性进行了证明,综合上述两部分的证明结果,给出了周期边界条件下的一维变系数线性对流扩散方程的稳定性结论。在稳定性证明过程中的关键是数值流通量的选择,本文选择一般数值流通量进行LDG格式的建立。其中令对流通量参数与各剖分单元边界点处变系数函数的取值相匹配是格式保持稳定的关键。然后通过建立分片全局投影,并以定理的形式给出其存在唯一性及最佳逼近性质的证明。结合Gauss-Radau投影与标准局部~2L投影,分别对对流项带变系数函数,扩散项为常系数和对流项为常系数,扩散项带变系数函数的对流扩散方程进行了误差估计,分析结果表明:对流项带变系数函数扩散项为常系数的对流扩散方程的LDG格式可达k阶收敛,而对流项为常系数扩散项带变系数函数的对流扩散方程的LDG格式可达k(10)1阶收敛,综合上述两个部分的证明过程,本文给出周期边界条件下的一维变系数线性对流扩散方程的局部间断有限元方法所得数值解为k阶收敛的结论。最后,在稳定性分析和收敛性分析的理论基础上,本文选择了三个有代表性的数值算例进行数值试验,根据计算结果可以看到,局部间断有限元方法在选择一般数值流通量的情况下用于求解变系数线性对流扩散方程可达k(10)1阶收敛,这是该方法作为数值算法求解高阶方程的最优收敛阶。总而言之,局部间断有限元方法是一种用于求解高阶偏微分方程的,具有高阶数值精度,高分辨率的数值算法。本文主要研究了选择了一般数值流通量的局部间断有限元方法在求解周期边界条件下的一维变系数线性对流扩散方程的稳定性和收敛性,给出了这一方法用于求解变系数方程的理论基础,并进一步说明了该方法在求解高阶方程时的高精度特点。同时利用数值算例为理论结果提供了事实依据,有力地说明了对于变系数方程,局部间断有限元方法具有计算稳定性和结果优越性。
其他文献
通过提出格库铁路南山口站至格尔木站间ITCS系统和CTCS-0系统覆盖范围不同的3种方案,并从ITCS系统和CTCS-0系统工程实施内容、车载系统配置、行车列控模式切换、司乘人员操作
普鲁士蓝(PB)及普鲁士蓝类似物(PBA)作为一种典型的金属有机框架(MOF)材料,具有催化、磁性、光学和传感等性质,已应用于电催化、药物传输、能量储备和生物传感等领域。本论文通过煅烧和与其他材料复合的方法,设计了几种基于双金属PBA的生物传感器用于检测目标分析物。首先我们设计并构建了一种NiFe PBA衍生的,中空双金属NiFe氧化物嵌入多孔碳的复合材料基适体传感器用来检测三磷酸腺苷(ATP)。
英国学徒制在经历了中世纪的辉煌和工业革命后的衰落,如今又以“现代学徒制”的形式,重新出现在英国社会中,并受到英国政府的重视。本文运用文献研究法、历史研究法和比较研究法,梳理与总结英国传统学徒制历史与现代代学徒制治理理念,并从政府、市场和学校三个主体,分析英国现代学徒制的治理结构及运行机制。最后,通过总结分析,对英国现代学徒制未来的治理改革进行展望。英国学徒制有其历史传承。始于12世纪的传统学徒制,