三维扩散方程非正交六面体网格的有限体积差分方法

来源 :清华大学 | 被引量 : 0次 | 上传用户:fine_yhy
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
扩散方程的非正交网格差分方法,要求有二阶精度,正交网格下等价于普通的二阶差分格式,并且计算要守恒,数值解收敛于方程的解析解.论文在二维扩散方程的有限体积差分方法的基础上,研究了在非正交六面体网格下三维扩散方程的有限体积差分方法,提出了一个具有很高计算精度,通量守恒,适应大变形网格的有限体积差分格式.取单元中心做为计算节点,减少向量和单元体积的计算,从而减小了计算量;把单元分成24个小的积分区域,以方便更好的计算单元界面上流量密度;利用通量守恒条件确定界面中心的函数值,保证方法的守恒性;对网格点采用了拉格朗日因子插值法,具有二阶精度,因为考虑了各插值点的相对位置,所以更适应非正交网格的计算;证明了格式在正交网格下等价于普通的七点二阶差分格式.在边界上,采用了不同于内部的处理方法,这样可以更精确的逼近边界,减少边界的离散误差.采用Bi-CGSTAB方法求解线性代数方程组,具有稳定的收敛速度,同时采用了不完全LU分解预处理技巧,提高线性代数方程组迭代收敛的速度.
其他文献
在现代医学领域中,2×2配对样本的非劣性检验是一个非常重要的研究课题。目前,基于2×2配对样本的非劣性检验方法主要有三种类型。一种是基于两率之差(RD)的非劣性检验方法。另一
该文给出了所有有非平凡主幂等元的p阶不可分环的分类及其矩阵表示,其中p是一个素数,特别地,这样的环共有35(p≠2)或33(p=2)个.
1976年,Gross提出了下述问题:能否找到一个集合S,使得对于任意一对非常数整函数f与g,只要满足E(S,f)=E(S,g),就有f≡g?该文研究了与之相关的唯一性问题.在第一章中,我们介绍
在第一阶段临床试验中,为了估计最大容忍剂量(MTD),O’Quigley Pepe和Fisher(1990)突破传统的up-and-down方法,提出了连续重估方法(CRM)。 最近,Gasparini和Eisele(2000)在CRM
该文针对现有汉语语音识别中常用的特征参数对于非平稳部分识别中的不足,对非平稳部分的特性进行探索,根据语音非平稳部分的统计特性,基于平稳性检验提出了新的有效度量语音
本文的目的是研究具有P-正则断面的正则半群.首先建立了正则半群上的P-正则断面的定义及与其相关的一些概念.为了刻画具有P-正则断面的正则半群,讨论了P-正则断面的基本性质.
该文分为两个相对独立的篇章:第一部分主要讨论了Kac-Moody代数中的一类基本问题,即给定一个实根或虚根,其对应的实根向量和虚根向量该如何表示?我们要求给定的广义Cartan矩阵A
笔者试图通过多年的乒乓球选修课教学效果及学生的学习兴趣,探析高职院校乒乓球教学对学生终身体育有意识的培养方法,为日后的乒乓球教学和终身体育意识的养成提供理论参考.
如何贯彻落实“三个代表”重要思想,武宣县在实施“五村、两规范”建设活动中,以“党群带富帮扶联合体”为载体带领群众共同致富,通过组建各种联合体,党员充分发挥先锋模范作
顺应时代发展潮流,在传播新观念、探索新方法上下功夫在改革开放的新形势下,要使宣传思想工作与时俱进,反映社会现实,满足新时期广大干部群众尤其是青少年的要求,必须认真研