翼-身组合体跨声速流场计算

本文把Holst的AF2因式分解格式推广到求解机翼-机身跨声速绕流的全速势方程。利用求解椭圆型偏微分方程组来数值生成贴体的计算网格。算例表明本方法运算可靠,费时较少,可以向工程实用部门推广使用。     

翼—身—尾组合体绕流的Euler方程计算

本文采用重迭网格技术和Ｅｕｌｅｒ方程计算翼－身－尾组合体绕流。对翼－身与尾－身部分采用各自的Ｈ－Ｏ型网格。Ｅｕｌｅｒ方程求解采用Ｊａｍｅｓｏｎ有有限体积法，即中心差分挖和显式Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进。采用前后区交替迭代使前后两区通过重迭区交换信息。本文用ＮＡＣＡＴＮ４０４１翼－身－尾模型为例，计算的空气动力特性与实验符合较好。     

超声速前机身及翼—身组合体的全速势方程数值计算

通过数值求解全速热方程，计算了超声速来流的前机身及翼－身组合体。当流场存在亚声速区时，在此区域内采用中心差分格式，迭代求解，并引入多重网格技术，加快收敛；当流场中某些一区域沿某一方向是超声速时，在此区域内采用沿该方向的推进解法。计算结果表明，本文的方法可靠，结果准确，可以向工程应用方面推广。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

翼－身－尾组合体绕流的Euler方程计算

本文采用重迭网格技术和ＥＪｅｒ方程计算翼－身－尾组合体绕流。对算身与尾。身部分采用各自的Ｈ。Ｏ型网格，Ｅ＊ｅ／方程求解采用１ａｍｅｓｏｎ的有限体积法，即中心差分近似和显式Ｒｕｎｇｅ·Ｋｕｔｔａ时间推进。采用前后区交替迭代使前后两区通过重迭区交换信息。本文用ＮＡＣＡＴＮ４０４１翼－身，尾模型为例，计算的空气动力特性与实验符合较好。     

大展弦比飞机翼身Euler绕流计算

给出了用Ｅｕｌｅｒ方程对大展弦比翼身组合体绕流的数值模拟，将Ｋｒｏｌｌ的求解欧拉方程的有限体积法，推广到求解机翼机身流动问题，基本思想是使用保角变换、代数变换和椭圆方程变换混合方法生成机翼机身贴体网格。     

基于欧拉方程的跨声速翼型设计

给出了一种基于欧拉方程的跨声速翼型设计方法。方法以Ｔａｋａｎａｓｈｉ提出的“正反迭代余量修正”设计原理为基础，在气动力分析模块中，以欧拉方程为基本控制方程，并采用Ｗａｌｚ边界层方法对其进行粘性修正；反设计模块采用经过改进的二维翼型设计方法。方法的程序经过几个设计实例主题证明是有效而实用的。     

跨声速欧拉方程的一种隐式有限元解法

本文用有限体积Ｇａｌｅｒｋｉｎ法求解了跨声速二维欧拉方程。求通量运用了Ｏｓｈｅｒ格式，通过外插把原来只有一阶精度的逆风格式改进成二阶精度。并采用基于一性化流量的隐式格式和非结构化网格，使程序具有较高的收敛效率。     

翼-身-尾组合体的网格生成及 Euler 方程数值模拟

采用代数方法快速生成翼-身-尾组合体三维O-C型贴体网格。该代数网格生成方法使用双边界法,可保证物面附近近似正交,真实地模拟了机身头部及与机翼、尾翼交界处的形状,提高了代数方法生成网格的质量。并用中心格式有限体积法求解跨声速Euler方程,以某翼-身-尾组合体模型为例,计算结果与实验符合良好。     

超声速翼—身组合体Euler方程数值解

本文采用MacCormack二步显式差分格式,用空向推进法给出了超声速翼-身组合体Euler方程的数值解。应用薄翼假设生成简易的多区网格。采用Kentzer方案处理边界。计算中通过对横流速度型的修正改善了背风面压力分布特性,解决了推进中断问题,同时模拟出背风涡。     

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

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

欧拉方程多重网格方法及三维跨声速大迎角流动数值计算

本文给出一种粗、细网格时间步长可调整匹配至最优组合的多重网格方法以获得在最细网格上最大的加速收敛效益，二重网格中节省机时可达３５％。文中的数值方法可用于飞行器跨声速大迎角无粘有旋的涡干扰流场的模拟研究。     

机翼机身跨声速绕流的计算

本文采用变换与数值网格生成方法，对由椭圆抛物而描述的每一个展向曲面进行网格生成，引入一种用于控制网格疏密分布的函数变换方法，其特点是易于对不同族网格线分别控制，根据物面梯度自适应地调整网格的疏密。用有限体积法计算守恒型全速势方程跨声速绕流，对M-100翼身组合体进行计算。理论计算结果与实验结果表明，本文方法有效。     

一种求解欧拉方程改进的通量分裂方法

本文以Steger-Warming提出求解Euler方程的通量分裂方法为基础,利用改进的特征值分裂方法,建立了一种连续的通量分裂方法。该方法克服了Steger-Warming通量分裂法存在的分裂后的通量在个别点导数不连续的问题,同时方法简单,便于推广。对于通量分裂后的Euler方程,采用迎风格式,建立了相应的二阶隐式有限体积方程。为了减少求解所得块矩阵方程组所需要的工作量,本文利用特征变换,在不改变原方程组离散化精度的前提下,将所得块矩阵方程组化为分离的代数方程组进行求解,从而大大减少了计算量。文中分别对一维管流,二维绕NACA-0012和NACA64A-10翼型在亚声速和超声速来流的跨声速流场进行了计算,并将所得结果与有关实验和计算结果进行了比较,均符合良好。 本文所建立的方法具有占机内存小,计算时间短的优点。在普通中小型计算机上亦可使用。     

三维非结构网格的欧拉方程聚合多重网格法

介绍一种基于三维非结构网格的聚合多重网格技术,用于求解Euler方程定常解,加速收敛效果十分明显。求解Euler方程的基本算法采用有限体积中心差分格点格式。聚合多重网格技术按一定的规则将细网格的若干控制体积聚合在一起,组成一个较粗的网格。数值实验表明本文多重网格法十分有效。