高精度有限体积法格式的构造方法及应用研究

来源 :哈尔滨工程大学 | 被引量 : 0次 | 上传用户:wangold
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
微分方程的数值方法已经广泛的应用于长期的天气预报、海洋工程、地下天然气的开发、造船业等各类工程问题。现有的有限差分法、边界元法、有限体积法、无网格法等是求解微分方程数值解主要采用的方法。从最近几十年得到的数值计算研究成果上看,有限体积法是一种比较活跃的数值计算方法,在微分方程的求解中得到了广泛地应用,并且有良好地计算效果,而用有限体积法得到的离散方程保持了微分方程的积分守恒性、每项也都有明确的物理意义,并且得到的离散方程具有一定的规范性,方便于在计算机上编程实现计算。但是现有的高精度有限体积法格式很少,满足不了工程上对高精度的要求,因此,构造高精度有限体积法的离散格式是非常有必要的。  本文以对流扩散方程为研究对象,在利用有限体积法得到离散格式的过程中,如果采用Lagrange多项式插值方法得到的函数代替方程中的未知函数,那么再在控制体上做积分就会得到对流扩散方程的离散格式,但是在估计离散格式的空间精度时,会遇到需要估计插值多项式截断误差导数大小的问题,而截断误差导数的大小又难以估计,为了解决这个问题,给出了构造一种高精度有限体积法格式的新方法——多重积分有限体积法,本文的创新之处在于,首先在方程两边对空间变量做变限积分,然后再在控制体内做积分,这样做积分可消除方程中未知函数的导数项,而在控制体内我们又加入了一个控制变量,可以自由的控制积分区域的大小,这样得到的离散方程的方法具有一定的普遍性。其次利用Lagrange多项式插值方法得到插值函数代替方程中的未知函数,结合对时间项的处理(采用中心差分格式,这样得到的时间精度要比一般的差分格式要高,也能与空间精度相吻合)给出了全离散格式,且对离散方程进行了误差估计,最后利用Fourier分析法讨论了全离散格式的稳定性。由此方法得到,对一维对流扩散方程来说,由于作插值函数的时候只用了三个节点,故时间精度是一阶,空间精度是一阶,且全离散格式是无条件稳定的;对二维对流扩散方程来说,作插值函数的时候用了九个节点,故时间精度是二阶,空间精度是四阶,且离散格式是无条件稳定的,而给出相应的算例,进一步验证了全离散格式的稳定性及其精度符合理论推导。应该指出的是利用本文中给出的构造离散格式可以得到更高精度的离散格式。
其他文献
人工神经网络是一种非线性信息处理系统,被广泛应用于组合最优化、模式识别、图像处理、金融预测、信号处理、自动控制、人工智能等领域,具有巨大的发展潜力和研究价值.在目
同宿、异宿轨道作为动力系统理论中一类非常有趣的不变集,曾引起了许多专家学者的关注.人们知道Smale马蹄为我们描述了混沌的动力学行为,那么什么会触发混沌由Birkhoff-Smale
医院感染是当今突出的医疗安全问题,医院感染一旦发生,控制相对困难,影响患者预后与转归,增加医疗费用,给患者和社会带来巨大的经济负担,严重的可以导致患者残疾或死亡。因此对医院
摘 要:电气自动化是研究与电气工程有关的系统运行、自动控制、电力电子技术、信息处理、试验分析、研制开发以及电子与计算机应用等领域的一门学科。电气自动化控制是增强生产、流通、交换、分配等关键一环。  关键词:人工智能 电气 自动化控制  人类智能主要要包括三个力面,即感知能力,思维能力,行为能力,而人工智能是指由人类制造出来的“机器”所表现出来的智能。人工智能主要包括感知能力、思维能力和行为能力。人
Uncertainty on the geological contacts and the block volumes of the models along boundaries is often a major part of the global uncertainty of reserve estimatio
Keller-Segel方程组是生物数学中的一个重要模型,刻画种群发展中的趋化现象,数学上属于具有交叉扩散的抛物-椭圆型或抛物-抛物型非线性偏微分方程组.本文研究一个抛物-抛物型Ke
调和分析于18世纪初形成,调和分析又称Fourier分析.它最原始研究的内容是函数中的Fourier变换和求和法,随着时代的发展,进一步的研究了Hardy-Littlewood极大函数,并且此函数已经成
布尔函数在密码学中占据极为重要的地位,其密码学性质的好坏决定着密码算法的安全性,具有重要的研究价值。  在流密码中,代数攻击的提出和发展对布尔函数的设计提出了更高的要