高阶精度DG/FV混合方法在二维粘性流动模拟中的推广

DG/FV混合方法因其具有紧致性、易于推广至高阶及相比同阶DGM计算量、存储量小等优点,已成功应用于一维/二维标量方程和Euler方程的求解。在此基础上,将该方法推广于二维三角形/矩形混合网格上的NavierStokes方程数值模拟,将格式形式精度提高至4~5阶。物理量的空间重构及离散使用DG/FV混合重构方法;无粘通量计算采用Roe格式;粘性通量计算采用BR2格式;时间方向离散采用高阶显式R-K方法或隐式方法。利用该方法计算了有解析解的Couette流动问题以验证几种格式的数值精度阶,并计算了层流平板流动和定常、非定常圆柱绕流问题等经典算例。计算结果表明DG/FV混合方法达到了设计的精度阶,在较粗的网格上亦能得到高精度的计算结果;定性分析和数值结果表明相比同阶DG方法单步计算量减少约40%。     

点火冲击波与药柱裂纹流固耦合作用机理数值模拟研究

流体计算采用基于非结构动网格的有限体积方法,结构计算采用分段精细时程积分方法,求解结构动力学有限元方程,流体计算和固体计算之间采用松耦合在不断运动的交界面上进行作用力和运动特性的传递,对固体火箭发动机点火形成的冲击波与药柱缺陷形成的弹性裂纹之间相互耦合作用过程进行了模拟。研究表明,作用在弹性裂纹上的冲击波载荷引起的壁面变形随时间变化明显,开始表现为撑开裂纹出口,随后出现了收缩出口撑开底部的瓶状;固体变形对冲击波在弹性狭缝内的压力波振荡影响较小,以前采用刚体模型得到的研究结论仍有指导意义。     

空间曲面的非结构网格生成

介绍一种空间曲面的 Delaunay三角形网格生成技术。采用参数表示法将空间曲面变换为平面参数域,引入相关的度量矩阵后将平面参数域进行三角剖分,最后得到空间曲面的三角剖分。空间曲面的边界点作为已知点被给定,内点采用一类自动生成的方法计算,例如按边插点法,Voronoi边插点法等。该方法可以很好地控制曲面网格的质量,使用方便,可以直接用于物理问题数值模拟中的曲面网格生成,也可以用于三维非结构网格生成中提供边界面三角剖分。文中给出的算例说明了本方法的应用情况。     

基于点簇合并的多边形表面模型简化算法研究

 在计算机图形的动态实时显示技术中, 将高度密集的多边形表面模型用尽可能少的多边形逼近成为近年研究的焦点, 也出现了一批简化算法。针对存在的问题, 文中以累进网格为基础, 并将点簇合并引起的几何逼近作为研究重点, 将多种度量方法通过权因子统一到基于二次误差的模型逼近度量方法中, 并提出了改进的基于二次误差的多边形表面简化算法。实例表明该算法可快速、高质的实现逼近。该算法也适用于实时环境下的多边形网格的动态多分辨率表示。     

复杂外形的动态混合网格生成方法

本文将层推进方法、四分树方法和阵面推进法相结合,建立了二维复杂外形的动态混合网格生成方法.首先用层推进方法在运动物体附近生成贴体的四边形网格,然后用四分树方法在外场生成自适应的矩形网格,而中间由三角形网格过渡.三角形网格由阵面推进法生成,并由节点松弛和"对角线交换"方法进行优化.当物体运动时,贴体的四边形网格随物体运动而运动,而外场的矩形网格保持静止.中间过渡层的三角形网格随物体运动而变形.当物体运动位移较大导致出现质量较差的单元甚至导致网格相交时,则在局部重新生成网格.利用上述动态混合网格生成方法,生成了多种含相对运动的复杂外形的动态混合网格,网格质量令人满意.     

用非结构动网格技术计算确定再入飞行器配平攻角

提出了用数值计算直接确定飞行器配平攻角的自由振动法.该方法采用非结构动网格技术,使飞行器在气动力作用下绕通过质心的轴线自由振动,最终衰减到配平攻角.构造了一种MUSCL类型有限体积格式,求解考虑了动网格效应的三维非定常Euler方程,采用弹簧近似技术实现非结构网格的变形运动.作为例证,计算确定了载人飞船返回舱的配平攻角,计算结果与实验结果接近.该方法外形适应能力强,效率较高,具有较好的应用前景.     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法．在背景网格建立和三维初始推进面生成的前提下，利用推进阵面法生成了翼身组合体的三维空间网格，并应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

二维非结构网格上的高阶间断有限元方法研究

对二维非结构网格上的求解欧拉方程的高阶间断有限元方法进行了研究.运用间断有限元基本理论,采用多项式基函数对解进行近似描述,欧拉方程的数值通量项运用高斯积分公式和Roe格式求取,时间推进采用传统四步Runge-Kutta方法.运用发展的间断有限元方法对NACA0012翼型在亚跨音速情况下的无粘流动和圆柱的低速无粘流动进行了数值计算,结果表明:间断有限元方法具有良好的收敛特性、很小的数值耗散和很好的激波捕捉能力.     

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。