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

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

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

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

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

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

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

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

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

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

Euler及N-S方程的二维四边形网格自适应算法

描述了一种基于贴体网格的Euler和N-S方程自适应算法。采用中心有限体积法,对Euler及N-S方程进行数值求解,N-S方程的计算中加入了基于SA湍流模型的DES模型。在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离、涡等特性的捕捉和分辨能力。采用上述方法,数值分析了单段翼型的绕流问题。文末给出的2个二维问题的对比算例表明本文方法有效。     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

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

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

二维Euler方程SCB格式的并行实现

本文讨论二维双曲型守恒律方程组的高分辨ＳＣＢ格式，并应用这些格式求解二维Ｅｕｌｅｒ方程组。计算过程中在江南－Ⅱ型并行计算机上进行的，计算过程运用了许多并行计算技术，计算结果令人满意。     

固体火箭发动机燃烧室喷管一体化网格生成

给出固体火箭发动机燃烧室喷管一体化网格策略,用二维方法生成平面网格经旋转或堆砌形成空间网格。对网格正交性受几何域形状限制的情况进行了讨论。给出了一种燃烧室－喷管外形的网格结果。     

计算跨音速N─S方程组的二阶隐式矢通量分裂格式

提出了一种二阶精度隐式矢通量分裂单步差分格式并用于求解时间平均Ｎ－Ｓ方程组。与传统方法相比,该方法是Ｊａｃｏｂｉ谱半径格式和无近似因子分解。数值试验表明：该方法具有精度高、计算量小和收敛快等优点。在平面叶栅和扩压器进气道的跨音速流场的求解中,较好地捕获了激波、分离等现象,与实验数据比较结果令人满意     

三维组合药型单室双推力发动机燃面与内弹道计算

基于有限容积的思想，构造任意形状的三维药柱燃面退移通用计算方法，并推广应用于嵌金属丝推进剂的燃面计算。计算过程只关心燃面的始态和终态，可完全避开错综复杂的几何燃面变化规律，实现任意药型下多推进剂、变燃速燃面退移的计算机数值模拟。在此基础上，对三维组合药型单室双推力发动机进行了喉部烧蚀条件下的内弹道计算，同时与试验结果进行比较，取得了较好的效果。     

亚跨超音速喷管内流场的整体矢通量分裂算法

从N－S方程出发，采用先进的矢通量分裂算法，对亚跨超音速喷管内流场进行了整体模拟计算。文中对几种不同计算格式的计算效率进行了比较，给出了超音速占优混合流动的直接分割求解方法。计算采用代数湍流模型，跨音速段的计算结果同其它实验和计算结果进行了比较。计算中发现对出口边界条件进行亚音速、超音速和有回流分别处理很有必要；边界层内的M数沿壁面法向变化非常快；不同壁温条件和不同的流动Re数均会改变边界层的厚度。     

强制偏流喷管的三维网络生成

运用泊松方程实强了制偏流喷管内的三维网络生成，方程中的源项采用了与边界形状和边界网格点分布无关的数学形式，源项内的自由参数由横向网格线与边界的交角及横向网格线在边界处的局部曲率两个约束条件来确定，实现了由边界上的网格分布直接控制内部的网格点，使生成的网格具有良好的正交性和合理性。     

固体火箭发动机有约束流场推力特性计算

为计算带约束流场的固体火箭发动机推力特性，建立了一维准定常变截面的内弹道计算和约束流场流动的数学模型。计算结果表明，约束流场使推力增加14%，冲波进入喷管扩张段。与试验结果很接近，表明计算模型正确，可用于其它型号带约束流场的推力特性计算。     

网格分区构筑中区域边界的隐式处理方法

介绍一种网格分区构筑中区域边界的隐式处理方法──ADMZI法。该方法具有简单实用、数据结构灵活等特点,特别适合于区域边界难以确定的情况。将该方法应用于多级火箭级间分离时的流场物理域网格构筑,结果证明是成功的。     

燃气轮机燃烧室中贴体网格的生成与应用

运用贴体坐标系统分别对燃气轮机燃烧室的冷态和热态流场进行了数值模拟。燃烧室的曲线壁面边界用Ｔｈｏｍｐｓｏｎ的非正交贴体坐标网格处理，并把ＳＩＭＰＬＥ算法应用到曲线坐标系下求解各守恒方程。采用双方程ｋ－ε模型来描述紊流特性；采用按简单的单步化学反应假设的Ａｒｒｈｅｎｉｕｓ－ＥＢＵ模型表述紊流燃烧；采用热通量法辐射模型来估算辐射通量和壁面温度。另外，还研究了冷态时不同旋流数和紊流强度对气流结构的影响。     

三维装药通用计算方法——体积元素法

本文提出一种新的三维装药通用计算方法——体积元素法.这种计算方法使计算过程的数学处理十分简单、方便.根据该方法编制的计算机程序可以对各种类型的三维装药进行快速、准确的计算.用该程序对一台具有复杂燃面变化规律的实际发动机进行了计算,获得了令人满意的结果.     

三维装药几何通用型面模拟计算法

本文对目前广泛采用的“坐标药柱通用模拟计算法”作了有益的改进.模拟装药表面的各基本几何型面更加灵活通用,数学过程也变得简单方便.可对各种类型的三维装药进行快速、准确的计算.通过对某翼柱型装药的几何计算,获得了令人满意的结果.     

喷管潜入段二维两相冷流实验

建立了喷管潜入段二维两相冷流实验设备,进行了相应的实验。实验清晰地显示了在潜入喷管进口背壁凹区旋涡的产生、稳定、消失,以及粒子的冲刷、粘贴、脱落、沉积等反复作用的复杂现象,它们与流动的雷诺数、粒子的含量和尺寸、潜入喷管进口处的形状和潜入深度、背壁型面等因素有关。