三维Euler方程的自适应非结构网络计算

介绍一种三维自适应非结构网格的Euler方程求解技术。采用有限体积中心差分的格点格式对Euler方程进行空间离散,多步龙格-库塔时间推进,并采用当地时间步长、焓阻尼修正等加速收敛。采用Delaunay三角剖分技术生成三维流场的四面体网格。网格自适应技术中的网格局部加密也采用Delaunay三角剖分技术。最后用ONERAM6机翼的跨音速绕流计算说明本文方法的有效性。     

三维Euler方程非结构自适应网格投影和光顺技术研究

研究了三维Euler方程非结构自适应网格的投影和光顺技术.以ONERA M6机翼为计算模型,非结构网格为计算网格,守恒变量的范数为自适应加密准则,对采用一般自适应计算策略得到的网格和数值结果进行了分析,并发展了一种考虑网格投影和光顺的自适应计算策略.数值试验结果表明:通过对自适应所生成网格的投影与光顺,实现了自适应计算过程中的几何保真性,获得了相对光滑的压力分布曲线.     

Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

二维Euler方程的自适应非结构网格生成方法

发展了一种用于数值模拟二维Euler方程的自适应非结构网格的生成方法。这一方法根据背景网格提供的信息直接确定内部结点的位置，然后依照Delaunay准则，将新结点加入到已有网格划分中。用于平观流及管道流的实际计算中，结果表明，该方法不受流场形状及流动特征的限制，具有较好的通用性。     

二维Euler方程的非结构网格再生成自适应方法

给出了提高二维Ｅｕｌｅｒ方程定常解质量的非结构网格自适应方法和模拟结果。流场算法为结点型有限体积法和Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式，网格用阵面推进法生成，生成初始网格时使用均匀的直角坐标背景网格。     

三维非结构混合网格生成与N-S方程求解

 提出了一套较为通用的, 完全自动化的非结构混合网格生成方法, 在物面粘性作用区, 运用一种改进的推进层方法生成三棱柱形和金字塔形网格; 在其他流动区域采用阵面推进方法生成四面体网格。采用一种改进精度的格心有限体积法和Spalart-Allmaras 一方程湍流模型对三维N-S 方程进行了求解, 并以M6 机翼、DLR-F4 翼身组合体和某运输机外形的粘性流场作为数值算例, 验证了上述网格生成和流场求解方法的正确性和实用性。     

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

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

二维欧拉方程的非结构网格自适应方法

给出提高二维欧拉定常解质量的非结构网格自适应方法。网格的生成采用Ｅｅｌａｕｎｅｙ方法,流场的算法为节点型有限体积法,能量的计算采用Ｖａｎｌｅｅｒ格式,自适应采用网格的增加和删减。对流场特性的探测采用密度物质导数的绝对值。最后计算了Ｎａｃａ００１２翼型的跨音速绕流,结果证明自适应效果很好。     

三维直角叉树切割网格Euler方程自适应算法

基于叉树数据结构，实现了一种用于三维直角叉树切割网格的自适应算法，包括对几何外,形以及对流场计算的自适应网格加密技术。在初始网格的生成过程中，根据相邻网格的物面法向向量间的差值，进行针对外形的自适应网格加密；在流场计算中，根据相邻网格间选定物理量梯度的变化，进行针对流场的自适应网格加密。详细地描述了三维直角叉树切割网格的生成过程，以及以任意网格的切割细分算法。在自适应过程中，分别采用了八叉树和全叉树的数据结构，八叉树是基本的数据结构，而全叉树的采用，使网格具有了各向异性的特征，从而大大的减少了自适应网格的数量。采用中心有限体积法，求解Euler方程，并运用上术方法，完成了对外形和流场的逢适应网格加密算法，获得了较好的数值计算结果，证明了自适应算法的正确性，体现了直角叉树切割网格自适应技术的有效性和实用性。     

二维Euler方程的自适应非结构直角网格算法

本文用四叉树型非结构直角网格和有限体积法求解了二维Ｅｕｌｅｒ方程,在基本网格生成过程中,采用了以几何外形为基础的网格加细方法,包括物面边界网格加细和曲率加细;并提出了一个用自适应的方法在激波处及其附近局部加细基本网格过程,提高了激波分辨率,在计算中将结构网格经常采用的以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法做了推广,减少了计算工作量,节省了机时,为今后推广到三维创造条件,本文对单个翼型,多段翼     

非结构网格下Euler方程的高分辨率高精度解

提出了一种非结构网格下求解Ｅｕｌｅｒ方程的高分辨率高精度迎风格式。以Ｒｏｅ的矢通量差分分裂为基础,吸收了ＮＮＤ格式的优点,使其具有捕捉激波和滑移线的良好性能;在时间方向上采用Ｊａｍｅｓｏｎ的Ｒｕｎｇｅ－Ｋｕｔａ积分,并结合局部最大时间步长和残差光滑技术加速收敛。最后成功地完成了二维平板激波反射、跨音速Ｌａｖａｌ喷管内的流动和ＧＡＭＭ超音速前台阶绕流等算例,显示了该方法的有效性     

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

用非结构直角网格和欧拉方程计算飞机绕流

 采用八叉树结构 ,生成复杂外形绕流计算的非结构直角网格。物面附近用投影方法 ,使网格贴体。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。对歼击机模型的绕流计算表明 ,网格生成的机时花费很少 ,总体质量好 ,因而欧拉方程解算的收敛质量也好。     

二维欧拉方程在非结构网格上的自适应有限元算法

发展了一种易于和自适应网格加密技术相结合的非结构三角形网格生成技术。这种方法的特点是域内节点可任意规定;根据欧拉方程解的误差估计,在流场某些区域增加节点后,能重新快速生成新的具有高质量三角形单元的非结构网格。欧拉方程采用Ｇａｌｅｒｋｉｎ空间离散和Ｒｉｃｈｔｍｙｅｒ两步显式时间推进相结合的有限元法求解,并以沿流线熵的变化量作为物理判据进行网格自适应加密,提高解的精度。     

用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。     

三维Euler方程自适应八叉树结构各向异性直角网格算法

采用各向异性直角网格对相应的三维复杂流场进行了模拟 ,提出了物面三角化方法及网格与物面的相交判断过程。用以中心差分为基础的 Jameson的有限体积法对 Euler方程进行求解。计算结果表明 ,各向异性直角网格算法对相应流场进行模拟 ,计算时间较短而效果较好     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

二维可压无粘流的自适应流量修正有限元解

从非定常Ｅｕｌｅｒ方程组出发,采用非结构三角形网格,研究二维无粘流动问题的自适应流量修正有限元解。将ＦＥＭ－ＦＣＴ与网格自适应相结合计算了ＮＡＣＡ００１２翼型的三个算例。计算结果与实验及其他计算结果比较,亚音速和超音速绕流符合很好;跨音速绕流情况存在一些差别。     

二维非定常欧拉方程的摄动解在超声速矩形翼上的应用

本文利用超声速线化理论给出的矩形翼载荷系数解析解作为权函数,对沉浮和俯仰谐振的二维平板翼压强系数进行加权修正,获得矩形翼谐振气动力的近似解,计算结果与线化理论数值结果符合得很好,本方法是超声速和高超声速矩形翼非定常气动力的一种统一和快速的计算方法。