多重网格数值求解不可压流体的局部Fourier分析

来源 :昆明理工大学 | 被引量 : 1次 | 上传用户:helloliuhh
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
多重网格算法是偏微分方程数值求解的一种快速算法。主要针对离散微分方程后所得的代数方程组进行数值求解,在椭圆型偏微分方程的数值解中已被证明是最优的数值算法,其收敛性与网格尺度的大小无关,且计算成本与问题的规模成正比。由于多重网格算法的优越性,使得它成为计算流体力学中一种高效的数值方法而受到广泛关注和研究。本文依托高等学校博士学科点专项科研基金(优先发展领域)(20135314130002)项目、国家自然科学基金面上项目(51279071),研究多重网格法在水力机械内部流数值模拟方面的理论和应用,重点是多重网格光滑理论中的局部Fourier分析方法,对数值求解不可压缩流体控制方程的多重网格方法进行收敛性分析。主要研究内容和创新如下:(1)结合水力机械流道湍流的流动特点,提出了多重网格算法及其误差迭代的格式。基于局部Fourier分析理论,分别定义了离散算子和松弛迭代算子的椭圆率和光滑因子,并利用不同粗、细网格层Fourier组分之间的关系,定义了新的不变子空间,分析了不同粗化方式下网格转化算子的Fourier表述方式,研究了多重网格算法渐进收敛因子的理论计算方法,创新了两色松弛在两种不同的Fourier模态函数不变子空间中的光滑分析方法,得到了基于多色松弛矩阵的Fourier分析的理论表示,并以泊松方程为例给出了相应的分析结果。研究表明,基于多色松弛的多重网格光滑分析过程具有一般的迭代格式,所得结果具有代表性和应用前景。(2)基于交错网格和非交错网格提出了求解Stokes流的离散格式,并对该离散系统实施两种不同多重网格的松弛算法进行了光滑分析:即聚松弛和分布松弛光滑分析。在交错网格的离散系统中实现了多重网格分布松弛,发现该离散系统的光滑性取决于Laplace算子,并得到了相应的光滑因子。其次,在非交错网格离散系统中,分别实施了多重网格分布松弛和聚松弛,在两色松弛的Fourier谐波空间中,讨论了这两种松弛的光滑性质,得出光滑因子关于附加人工压力项参数的表达式。结果表明:松弛方法的收敛性与网格尺度无关,而依赖于附加人工压力项参数。(3)基于最优红黑Jacobi逐点松弛方法,从理论上分析了Possion方程两层网格算法的收敛性。给出了对流扩散方程的一阶上迎风离散格式,分析了对流占优参数和扩散参数对该离散格式的椭圆率影响,探索了对流扩散方程各参数对多重网格光滑性和两层网格收敛性的影响。在提出的理论方法基础上,利用Riemann解的通量差分分裂法-Godunov方法处理Oseen流控制方程的离散,得到了基于一阶上迎风格式的离散方程,并分析了使用多重网格方法求解该离散方程的V-循环算法和W-循环算法的收敛性,并通过局部Fourier分析方法,对获得的离散方程的聚对称交替线Gauss-Seidel松弛的光滑性质进行了系统研究。结果表明:使用多重网格的两层网格及三层网格算法求解具有不同Reynolds数的Oseen流,即便是在较高Reynolds数情况下,聚对称交替线Gauss-Seidel松弛仍然具有很好的光滑性质,且W-循环算法收敛性比V-循环算法好。(4)首次对基于非定常不可压缩流体的NS方程进行基于交错网格离散系统实施多重网格分布松弛。通过局部Fourier分析,发现该离散系统的光滑性质由时间依赖的对流扩散算子决定,并对两种处理时间依赖问题的多重网格松弛,时空松弛和波形松弛进行了系统研究。在交错网格上,提出了非定常不可压缩流体NS方程仅对空间变量进行离散的半离散格式,并对该离散系统实施分布松弛,使得离散系统多重网格松弛的光滑性质仅取决于时间依赖的对流扩散算子。通过局部Fourier分析,对时间依赖的对流扩散问题所使用的时空多重网格方法和波形多重网格方法进行了光滑性分析。另一方面,在时空多重网格方法的光滑分析中,采用了时空离散格式,其中时间离散采用一阶Euler向后格式,而空间离散采用一阶上迎风格式。提出了多重网格的粗化仅对空间粗化的半粗化方法以及与时空多重网格对应的各种松弛的局部Fourier分析方法。而在波形多重网格方法中,首先利用Laplace变换将时间依赖问题转化为带有复参数的定常问题,然后对应用于波形多重网格方法的各种松弛进行局部Fourier光滑分析。通过提出的两种多重网格方法的光滑分析,研究了对流占优参数和雷诺数对各种松弛算子光滑性的影响,给出了相应的最优光滑因子和最佳松弛参数的选取方法。提出的理论和方法部分用于了由导师负责的国家基金面上项目“水轮机旋转湍流全欧拉并行多层网格模拟研究”等项目的算法设计和代码开发应用中,并获得成功。
其他文献
通过对铝基板及不锈钢基板上条纹沟槽表面接触角的测试,研究了条纹沟槽对材料表面润湿性的影响。试验用条纹沟槽表面采用准LIGA技术制备。该技术以波长为300~400μm紫外光作为
当前我国高职教育取得长足的发展,毕业生就业率稳步提升。摆脱过时的毕业生跟踪方式,建立一套完备的跟踪服务系统至关重要。在功能上,系统分为前台用户操作模块和后台系统管
改革开放以来,流入我国的FDI规模迅速增长,而近年来就业日益成为中国理论界和政府部门关注的一个焦点问题。本文以中国东、中、西部地区为横截面单元,利用各横截面单元1985年~200
边坡土体中的膨胀力是导致膨胀土边坡失稳的主要原因之一。文章基于湿度应力场理论,采用ANSYS软件的热传导分析功能模拟分析边坡的降雨入渗以及膨胀变形,并采用有限元强度折
纺织服装业属于典型的劳动密集型产业,由于其要素密集度的特性与中国和东盟各主要国家现阶段要素禀赋结构相吻合,因而近年来得到了飞速的发展:但从另外一个角度来看,中国和东
制备均分散纳米复合陶瓷的关键是使基质中的纳米相不团聚而均匀分布在其它组分中。介绍了纳米颗粒的分散以及与陶瓷基质均匀混合的新工艺
硫酸在有机化学反应中的作用甘肃省纺织工业学校邢桂芳在有机化学教学中涉及到很多用到硫酸的反应。例如,烯烃间接水合制醇的反应,醇分子内脱水生成烯烃和醇分子间脱水生成醚,醇
目的:了解非甾体抗炎药(Non—steroidAnti-inflammatoryDrugs,NSAIDS)在眼科围手术期的使用情况,评价该类药物对眼科手术术前及术后非感染性炎症的预防和抗炎效果。方法:通过随机抽
2018年11月9日,国家综合性消防救援队伍授旗仪式在人民大会堂举行。中共中央总书记、国家主席、中央军委主席习近平向国家综合性消防救援队伍授旗并致训词,这是全国消防救援