h-型自适应网络方法在计算流体力学中的应用

来源 :北京大学 | 被引量 : 0次 | 上传用户:lolotang530
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
计算流体力学(Computational Fluid Dynamics,简称CFD)是当代迅速发展的一门学科,是利用高速计算机求解流体流动的偏微分方程组,目的是为了更好地从定性上和定量上了解流体流动的物理现象,改进设计的一门学科.目前在航空航天、造船、气象、海洋、水利、液压和石油化工等工程领域都有广泛的应用。   本文主要研究了h-型自适应方法在计算流体力学中的应用.对于不可压两相流模型和机翼绕流模型,作者在非结构的三角形网格上离散求解,并把网格自适应方法应用到这两个模型的求解中。   对于两相流模型,作者选用不可压Navier-Stokes方程作为控制方程,并采用水平集方法描述界面位置.根据两相流模型的控制方程和水平集方程的不同特点,作者把两套网格自适应方法[44]应用到该模型的求解中.在本文的算法中,Navier-Stokes方程和水平集方程分别是在两套不同的自适应网格上离散求解,以便得到更高的计算效率和更好的逼近精度.水平集方法中的一个重要操作是重新初始化,就是在求解水平集方程捕捉界面的过程中为保持水平集函数是一个符号距离函数而引入的操作.重新初始化操作往往计算量比较大,而且会在一定程度上移动界面的位置.本文针对重新初始化的这些缺点,在前人算法的基础上提出了一个新的高效重新初始化算法,这个算法在重新初始化过程中能更好地保持界面不动.对于水平集方程的求解,作者采用特征线法结合高阶Runge-Kutta方法,在求解的过程中使用Y.-J.LIU提出的向前向后的误差修正方法[46]修正.为了更好地节省计算量,作者把Narrow-Band方法应用到水平集方程的求解中.对于Navier-Stokes方程的求解,本文使用有限元中的混合元方法。   对于机翼绕流模型,作者选用定常Euler方程作为控制方程.该问题的计算过程中常常会出现较强的激波.为了能够精确的捕捉激波的位置和形态,需要在非常稠密的网格上对控制方程离散求解,而在远离激波的地方,由于解的光滑性比较好,可以在较稀疏的网格上对控制方程离散求解.因此,作者把网格自适应技术应用到机翼绕流模型的求解中.本文还提出了一个基于有限体积方法的在非结构网格上的定常Euler方程的高效稳定算法.这种新的算法是在牛顿迭代的基础上建立起来的.在每步牛顿迭代中使用以块LU-SGS迭代为光滑子的线性多重网格迭代.在算法中,作者还采用了依赖于局部残量的松弛因子来正则化牛顿迭代中的Jaccobi矩阵,代替传统的基于当地CFL数的时间步长的正则化方式。
其他文献
指数跟踪在一些被动投资基金管理中是非常重要的技术。但是,很多时候,在建立跟踪复制指数表现的策略时,无法买入该指数所有的成分股,这样导致跟踪的效果变差。因此,如何低成
软件复杂程度与软件模块中的缺陷分布有着显著的关系,一些软件复杂性度量元的度量结果与软件测试和验证过程中发现缺陷的数目密切相关。   本文研究了利用软件复杂性度量
学位
本文利用微分方程定性理论和分叉理论对一类非线性传染性疾病模型的动力学性质进行了研究,包括稳定性和分叉情况,对传染病模型的研究有重要的意义。   全文由四章构成  
中心-焦点的判别问题是常微分方程定性理论研究中最古老和最重要的问题之一,它与Hilbert第16问题有密切的联系,并且在上个世纪得到了极大的发展。Poincare和Liapunov提出了最基
本文主要考虑了随机矩阵理论中的几个问题:随机Toeplitz及相关矩阵的极限谱分布和谱涨落的中心极限定理;限制迹的Gauss正交、酉、辛系综;量子混沌中的随机纯态的双粒子纠缠模型,
学位
今年以来,舒城县在全县推行了发展党员会审制。该制度规定,各级党组织在发展党员时,对拟发展对象在培养考察期间,每半年由组织、纪检、信访、计生、综治等部门对入党积极分子
学位
课堂导入是教师引导学生参与学习的过程和手段,是课堂教学的必需环节.因此,恰当的导入有利于营造良好的教学情境,便于集中学生的注意力,更能激发学生的学习兴趣,启迪积极的思
作为小学数学教师,要俯下身子看孩子,从学生的成长进步、全面发展出发,把握新课程改革的精髓,与学生沟通交流,掌握学生的身心发展特点、数学感知能力,倾听学生对数学学科教学
近二十年来,多智能体系统已迅速发展为分布式人工智能的一个重要分支,也是目前国际上人工智能研究的前沿学科,其研究的目的是将大而复杂的现实系统建设成小的、彼此互相通信