基于非结构混合网格的间断Galerkin有限元方法研究

来源 :西北工业大学 | 被引量 : 0次 | 上传用户:jacky1228
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
非结构网格消除了网格节点的结构性限制,适合求解复杂流场,高阶精度方法具有较小的数值色散和耗散,能够精细模拟复杂流场,它是提高数值模拟效率的主要手段之一,而大涡模拟(Large Eddy Simulation,LES)、直接数值模拟(Direct Numerical Simulation,DNS)等数值方法都需要高阶精度格式,因此开展基于非结构混合网格的高阶精度数值方法研究具有重要意义。间断Galerkin有限元方法(Discontinuous Galerkin Finite Element Method,DG)是适用于非结构网格的高精度紧致数值方法的典型代表,已被成功应用于多个领域,然而计算效率低和激波捕捉的不稳定问题严重制约了该方法的发展和应用。为此,本文发展了基于非结构混合网格的高精度间断Galerkin有限元方法,建立了一套实用的高阶精度DG计算程序,通过发展高效率的隐式迭代方法、高效的并行计算技术和激波捕捉技术,重点解决该方法中的计算效率及激波捕捉问题。为下一步建立分离涡模拟(Detached Eddy Simulation,DES)、隐式大涡模拟(Implicit Large Eddy Simulation,ILES)提供平台,并最终将方法用于高超气动力/热、发动机内流、旋翼流场、气动噪声等问题的数值模拟,以期改善数值模拟精准度。论文的研究工作主要包含以下几个方面:1、针对非结构混合网格,发展了基于Taylor基函数的高阶精度间断Galerkin有限元方法。为提高RANS(Reynolds Average Navier Stokes)方程数值模拟的鲁棒性,采用N-S方程与湍流模型方程松耦合方式,发展了DG方法求解主控方程,FV(Finite Volume Method)方法求解湍流模型方程的混合求解技术。通过等熵涡、平板层流边界层、库埃特流动以及平板湍流边界层等典型算例验证了发展的方法的计算精度。2、为提高DG方法计算效率,发展了基于残差Jacobi矩阵精确求解技术的隐式离散方法和并行计算技术。为提高求解定常问题收敛速度,结合PETSc库(Portable,Extensible Toolkit for Scientific Computation),发展了DG方法的GMRES(Generalized Minimum RESidual)隐式时间迭代方法,同时建立了LU-SGS(Lower-Upper Symmetric-Gauss-Seidel)方法;为进一步提高GMRES的收敛速度,针对无粘通量Roe格式和粘性通量BR2格式(Bassi-Rebay 2),通过链式法则,发展了残差Jacobi矩阵精确求解技术,该技术极大提高了隐式GMRES的计算效率;基于MPI(Message Passing Interface)发展了并行计算方法,解决了DG方法计算大网格量时的内存需求以及计算效率问题。通过多个典型算例分析了发展方法的计算效率,研究结果表明:基于Jacobi矩阵精确求解的GMRES能够增大高阶DG计算的CFL数,极大提高计算效率,同时能够改善特定条件下高阶计算的稳定性。3、为处理间断,在kxrcf间断探测器基础上,结合流场分析,根据单元变量梯度关系提出了一类新的kxrcf-grad间断探测器,针对弱间断问题设计了一种具有较好残差收敛特性的限制方法。发展了DG方法激波捕捉技术中的限制器方法,包括斜率类以及重构类限制器、间断探测器、负压负密度修正技术,更加精准地探测“问题单元”,提高了计算精度和效率;发展了Venkatakrishnan与HWENO(Hermite Weighter Essentialli Non-oscillatory)两类限制器,给出Venkatakrishnan限制器中自由参数K的影响规律;分析了限制器影响残差不降的原因,针对弱间断建立了具有较好残差收敛特性的限制方法,改善了限制器带来的残差滞止问题。针对强间断问题中可能出现非物理的负密度和负压力问题,发展了与DG方法相适应的修正方法,提高了计算的稳定性。通过双马赫反射、二维翼型、三维机翼典型算例模拟,验证了发展的激波捕捉技术的有效性。4、基于发展的高阶精度DG方法计算平台,开展了高阶精度DG方法在各类典型流动问题中的初步应用,进一步验证方法的计算精度和数值鲁棒性。算例涵盖了亚跨超全速域的典型算例,具体包括:二维层流流动、激波边界层干扰、30P30N多段翼湍流,三维DPW-W1机翼湍流流动、65°大后掠尖前缘三角翼湍流流动和小展弦比标模等多个复杂计算,研究结果表明:发展的DG方法能够用于层流与湍流、激波边界层干扰和大分离流动等复杂流动问题的精细化数值模拟。
其他文献
一直以来小波分析都是应用数学和工程技术等领域研究的热点问题,近年来,随着信息和计算机技术的快速发展,多小波和小波框架由于其自身的优点已广泛应用于信号处理、数据压缩和大数据分析等领域.对多小波和小波框架的研究已成为小波分析领域的重要研究课题,受到了学者们的高度关注,成为当下的研究热点问题之一.本文基于多分辨率分析(multiresolution analysis,MRA)的思想,主要研究多小波和小波
与传统的电磁兼容测试环境相比,混波室内形成的电磁场是统计均匀、各向同性和随机极化的。由于腔体的高品质因数,较小的输入功率可形成很大的场强,且能模拟受试设备所处的实际电磁环境,非常适用于进行辐射敏感度试验。搅拌模式混波室辐射敏感度试验方法具有更高的测试效率,在航空和军用装备领域应用前景广阔,但其场的变异性较大,可能导致试验结果具有更大的测量不确定度。IEC 61000-4-21:2011等标准尚未对
量子干涉测量可以实现光子量级的信号检测,是微弱信号精密测量的技术手段。作为评估测量性能的关键指标,量子干涉测量的相位灵敏度受到广泛关注,相位灵敏度的理论极限、影响因素和提升方法的研究是国内外重点关注的问题。本论文围绕量子干涉测量的相位灵敏度开展了四个方面的研究:相位灵敏度极限的研究;相位灵敏度内在影响因素的研究;相位灵敏度外在影响因素的研究;相位灵敏度提高的新方法研究。本文首先完成了相位灵敏度极限
随着基于深度学习(Deep learning,DL)的语音识别技术在实际场景中的大规模应用,提升声学模型的鲁棒性(Robustness)至关重要。一个鲁棒的声学模型应当对噪声、信道、说话人以及口音等扰动因素不敏感,在真实复杂的声学环境中依然能够表现出良好的性能。近年来,以神经网络(Deep neural network,DNN)为载体,深度对抗学习(Deep adversarial learnin
随着飞行器对动力要求越来越高,传统涡轮设计技术已经难以满足未来航空燃气涡轮发动机对涡轮性能进一步提升的需求。因此,综合考虑多种对涡轮性能影响因素的精细化设计成为主流,多种影响因素中涡轮轮缘封严气流对主流流动特性及性能的影响成为必须考虑的重要问题,深入研究涡轮轮缘封严气流与主流非定常干涉机理和损失控制方法具有重要的学术价值和工程应用意义。本文以1.5级涡轮为研究对象,对涡轮轮缘封严气流与主流非定常干
数字全息术(Digital holography,DH)是一种能够动态、定量、非接触、全场、高精度地获取原始物光波复振幅分布的干涉测量技术,被广泛用于复杂流场测量、生物细胞检测以及光学元件形貌检测等领域。除了可见光波段,在红外波段、太赫兹波段等,数字全息术也获得了广泛的应用。由于数字全息测量系统再现像的像质与最终的测量结果息息相关,因而如何改善再现像的像质得到了极大的关注。本文从影响数字全息测量系
异喹啉类生物碱是一类具有显著药理活性的含氮杂环化合物。阿朴菲类生物碱作为异喹啉生物碱中的一个重要类群,植物分布广泛、化学结构独特、生物活性多样。DactylicapnosineA(2-27)和 DactyllactoneA(3-1)是从民间药用植物紫金龙中分离出来的两个结构新颖,骨架独特的D环重排阿朴菲类生物碱。这两种高氧化度的骨架重排的生物碱此前从未报道。值得一提的是Dactylicapnosi
当前我国教育最重要的环节之一就是核心素养学习,每一科学科最终的目的都不是为了单科成绩,而是要通过积累知识,提高学生的学习能力和思维能力。小学阶段培养学生小学数学的核心素养也是非常重要的一项任务,其中包含很多内容,比如学生的创新能力、几何能力、运算能力以及推理能力等,所以怎样才能在教学中培养学生的核心素养,成为了小学数学教学中一项非常重要的任务。
图上的游戏和匹配问题是图论研究的两个重要内容,它们不仅对认识图的性质和结构有重要作用,而且在计算机科学、组合最优化和信息论等方面有着广泛的应用.本文首先借助传统的组合博弈理论研究了两类图上的游戏,分别给出了这两类游戏的博弈结果和最优策略.然后,根据匹配覆盖图的性质和耳分解定理,对匹配覆盖图上的非可行边集进行分析,得到了关于匹配覆盖图上非可行边集的相关结论.主要研究内容如下:首先,研究了一类毛毛虫树
周向重叠三分螺旋折流板电加热器壳侧复杂流道内的螺旋流和耦合涡旋扭曲椭圆管换热器壳侧复杂流道内的涡旋流动有助于强化传热。针对工业工程中常用的多圈环形分层布置U形管电加热器提出的“斜日字”单元布管三分螺旋折流板电加热器,具有换热系数高、流动阻力低、无流动死区、防结垢和易抗振等优点。螺旋折流板被认为是最有潜力代替弓形折流板的强化传热结构。轴向分隔则是一种采用小倾角螺旋折流板将相邻折流板适当分开以实现较大