旋翼桨叶跨音速绕流的全势方程计算方法

在随桨叶旋转的坐标系中，用有限差分法求解非守恒格式的全势方程，建立了一个适用于矩形桨尖和后掠桨尖桨叶跨音速流场计算的方法。该方法运用了Ｊａｍｅｓｏｎ的旋转差分技术，并采用等效迎角方法计入旋翼尾迹的影响。作为算例，首先对无升力状态下具有矩形桨尖和后掠桨尖的ＯＮＥＲＡ桨叶进行了计算，取得了与试验数据相吻合的计算结果。然后，对有升力桨叶，也分别就悬停及前飞状态下的桨叶剖面压力分布进行了计算，通过与试验值的对比，进一步验证了该方法的可靠性。     

基于混合网格的直升机旋翼/机身非定常干扰流场数值模拟方法

基于混合网格建立了一个适合于直升机旋翼/机身非定常干扰流场的数值模拟方法。其中,采用混合网格系统生成旋翼/机身干扰流场网格以降低网格生成难度;将二阶MUSCL格式和Roe格式相结合进行空间离散以提高流场的计算精度;同时引入隐式LU-SGS时间格式和OpenMP并行模式以提高干扰流场的计算效率。应用建立的计算方法,对Georgia和Robin旋翼/机身模型进行了数值模拟和分析,与试验结果良好的相关性表明了该方法的有效性。     

基于CFD的直升机旋翼流场及气动力计算

建立了一个基于结构运动嵌套网格的流场求解器,用来精确模拟复杂的旋翼流场,为更好地预测旋翼气动载荷提供一套计算方法.在该求解器中,控制方程为惯性坐标系下的三维非定常Navier-Stokes方程,空间方向上采用低数值耗散的Roe格式结合三阶逆风格式(Monotonic upwind scheme for conservation law,MUSCL),湍流模式采用了一方程的Spalart-Allmaras模型.应用所建立的方法,分别针对Caradonna模型旋翼、UH-60A直升机旋翼的悬停流场和7A旋翼、SA349/2直升机旋翼的前飞流场以及旋翼气动载荷进行了数值模拟.计算结果表明,本文的方法在一定程度上提高了对旋翼流场的模拟能力,进而提高了旋翼气动载荷的计算精度.     

轴对称进气道外罩跨音速非守恒型全位势流AF3算法

本文针对轴对称进气道外罩跨音速绕流问题的特点,采用了贴体坐标网格;并对这类网格拓扑构造了三种Baker的隐式近似因式分解AF3有限差分迭代算法。通过理论分析与数值实验,找到了既能计算亚临界情形,又能计算超临界情形的AF3格式;并就此格式讨论了适合于进气道绕流特性的AF3迭代格式中间变量的提法,具体提出了格式实施的相应方法,研制了可供分析任意轴对称进气道外罩跨音速特性的计算程序。由于采用了任意曲线坐标系上的轴对称非守恒型全位势方程适合于用不同方法生成的计算网格,因此,本文的工作还有可能进一步推广到计算更复杂的轴对称进气道跨音速绕流问题。计算表明,本方法收敛快,结果好。     

基于特征线理论的旋转爆震流场结构特征研究

特征线理论及其计算方法是气体动力学的经典理论与方法，应用于旋转爆震流场分析具有简单高效的特点。将坐标系建立在爆震波上，对旋转爆震流场进行简化，采用特征线理论并结合流场计算单元过程，建立旋转爆震流场计算模型。研究了当量比和喷注参数对氢气/空气、甲烷/空气以及辛烷/空气3种不同预混气的旋转爆震流场结构特征的影响。结果表明:爆震波高度和倾斜角度受混气当量比和喷注总温影响明显；燃料由小分子氢燃料变为大分子碳氢燃料时，爆震波高度和倾斜角度逐渐减小；混气当量比和喷注总温主要通过影响爆震波传播速度、高度和倾斜角度而影响爆震波后宏观流场特征趋势。     

二维非结构网格生成及其数值模拟

提出一种二维非结构网格生成算法以及一种新的绕翼型的数值计算格式。首先，在Delaunay生成算法的框架下结合了阵面推进算法的优点，从覆盖整个计算区域的三角形着手，给出判断准则，由这些三角形的边直接生成新点，让新点与符合要求的三角形的边的两点相连构成三角形，这样生成的三角形网格具有较好的几何性质；其次，针对二维N0012翼型跨声速绕流问题提出了一种新的计算格式。通过建立插值模板的概念，在每一个模板上运用最小二乘和MUSCL思想相结合的办法针对初始变量构造线性插值多项式，改进了原有的权系数的计算准则，将各插值多项式通过加权平均得到在整个模板集上通用的插值多项式，从而构造出高阶的有限体积格式；最后，计算结果表明本文工作具有一定实际应用价值。     

基于三维高分辨率CFD解法的气动优化研究

采用以高分辨率MUSCL格式为核心的三维CFD解法，与自动化网格生成技术等结合，构建了适用于各类亚、超音速气动优化问题的气动造型数值优化系统。分析表明，遗传算法对采用三维高分辨率CFD解法的优化来说工作量太大；通过多激波系超音速进气道优化设计和跨音速翼型阻力最小化的比较显示，基于目标函数梯度的F1etcher-Reeves法的优化速度比不使用梯度信息的Powell法快了近一倍。研究表明．梯度法是相对更适合气动优化系统的非线性优化方法。本文使用梯度法成功地对三维跨音速机翼进行了优化设计，获得了升阻比明显高于原始超临界机翼的设计方案。     

圆转方非轴对称喷管无粘流动的多层网格法求解

本文利用多层网格法求解三维无粘亚音和跨音速圆转方非轴对称喷管的流动,基本格式是把Ni[1]的二维有限体积积分格式推广到三维流场。边界条件采用文[2]中的“预测一修正”处理方法。为了进一步提高格式稳定性,在计算流动变量的一阶时间变化量时,采用单元体积平均法而不是文[1]中的算术平均法。引入当地时间步及多层网格方法,以加快格式的收敛速度。文中还对多层网格构造形式及其收敛效果作了讨论。通过对两种喷管模型流动的计算和与试验数据的对比,验证了该算法的精度。     

飞机亚音速非对称运动流场数值计算

本文以飞机亚音速对称运动（仅小攻角）流场数值计算为基础，采用涡格法理论开发飞机亚音速非对称运动（小攻角、小侧滑角、小角速度）流场特性数值计算。本文特点：从毕－萨公式出发推导的涡格法流场诱速公式比面涡法简便；采用多坐标系，详细推导了各坐标系转换公式；将复杂的干扰流场简化为一阶干扰和二阶干扰效应，并仔细推导了相应的流场诱速计算公式。已在ＩＢＭ－４３４１计算机上建立了流场计算实用性程序，既可单独运算，也可联入我国外挂物投放／发射过程数值仿真系统（ＣＳＳＰ）串联运行，提供所需流场特性。以Ｊ－ＸＸ型号飞机和两个翼身组合作为算例，定性分析了非对称运动流场计算规律；将两个翼身组合体对称动动流场特性计算值与实验值进行了比较。飞机亚音速非对称运动流场特性计算比对称运动情况更为复杂，难度更大。该工作对进一步开发ＣＳＳＰ系统具有重要价值。     

旋转流场测量技术的探讨和实践

本文专门介绍一种测量旋转流场的新技术。作者研制了一合旋转机械编码仪，结合应用在二维激光流速计系统，对几种水泵叶轮的内流场进行测量，获得了成功，得到了平均速度、湍流度、相关系数及雷诺应力等各力学量与相位角度之间的对应分布结果。     

NEW APPROACHES OF MUSCL SCHEMES FOR COMPRESSIBLE FLUIDS

The MUSCL scheme for compressible gas dynamicsis studied in this paper. We propose a new type of Euleri-an MUSCL scheme, which evaluates the intercell flux athalf time step. The intercell flux is comptuted from char-acteristic equations with the data which traced backthrough the approximatecharacteristic from the edge athalf time step. The data is classified in several differenttypes, depending on the characteristic directions of data.We also present a general procedure of oscillation-fre…     

地效翼风洞试验支架干扰数值分析

对地效翼移动地面风洞试验研究中的支架干扰进行了数值模拟和分析。分别对独立地效翼,带支架的地效翼及独立支架进行了数值模拟,计算采用可实现的κ-ε模型,通过求解定常不可压N-S方程,得出地效翼及支架周围流场分布情况。对几组计算结果比较分析了支架和地效翼的空气动力及由于干扰引起的空气动力,发现支架与地效翼之间的相互干扰随着地效翼迎角的增大而增强,如果忽略流动干扰造成的空气动力变化,地效翼升力误差很小,但阻力误差相对较大。同时对有干扰下和没有干扰下的流场进行了对比,分析了支架对翼尖涡流动及绕机翼流动的干扰。翼尖涡在地效翼翼尖附近的发展在0.5犮范围内基本不受支架的干扰;除支架对流场产生干扰外,移动带区域以外的固定地面附近粘性流动也对绕地效翼流动有一定的影响。本研究分析了风洞试验结果的可靠性,为地效翼风洞试验优化设计和地面效应风洞试验研究提供了参考。     

真实火箭燃气流场的测量与数值仿真

对用激光－莫尔偏折法测出的真实火箭燃气射流流场进行了数值仿真,在计算中采用了二阶精度的ＭＵＳＣＬ格式,将三维问题转化为二维问题,结果表明,理论计算与实验测量结果符合较好。     

A New D-GCL for Unidirectional Motion with Large Displacement

Numerical simulations of unsteady flow problems with moving boundaries commonly require the use of geometric conservation law(GCL).However,in cases of unidirectional large mesh deformation,the cumulative error caused by the discrete procedure in GCL can significantly increase,and a direct consequence is that the calculated cell volume may become negative.To control the cumulative error,a new discrete GCL(D-GCL)is proposed.Unlike the original D-GCL,the proposed method uses the control volume analytically evaluated according to the grid motion at the time level n,instead of using the calculated value from the D-GCL itself.Error analysis indicates that the truncation error of the numerical scheme is guaranteed to be the same order as that obtained from the original D-GCL,while the accumulated error is greatly reduced.For validation,two challenging large deformation cases including a rotating circular cylinder case and a descending GAW-(1)two-element airfoil case are selected to be investigated.Good agreements are found between the calculated results and some other literature data,demonstrating the feasibility of the proposed D-GCL for unidirectional motions with large displacements.     

二维双曲型守恒律的一类MmB差分格式

基于通量分裂、单元平均分片线性重构及逆风特性进行空间离散,构造了二维标量非线性双曲型守恒律的一类新的二阶精度的半离散差分格式。进一步地利用二阶TVDRunge-Kutta离散方法对时间进行离散,得到了一类新的时空二阶精度的全离散差分格式,并证明了格式的MmB特性。之后,将格式按分量形式推广到二维非线性双曲型守恒方程组。该方法的一个主要优点是使用分量形式格式计算二维非线性双曲型守恒方程组,无须解黎曼问题且不用进行局部特征分解,因而形式简单、计算量小。通过计算二维可压缩流欧拉方程组的几个算例,数值结果表明,该格式具有高精度、高分辨率及计算简单的特点。     

不对称喷管欠膨胀超声速射流的数值模拟

采用改进的MUSCL格式解三维可压缩平均雷诺纳维尔 斯托克斯方程组,湍流模型为Spalart Allmaras代数模型,数值模拟了四种调整片不对称喷管在喷口压比为2.0、2.8、3.4和4.0条件下的欠膨胀超声速射流场,获得了射流场流谱和参数分布,并给出了单调整片不对称喷管出口压强比与射流偏转角和扩张角的分布曲线,计算结果和实验数据符合良好。     

双曲型守恒律的无振荡中心差分格式

讨论了双曲型守恒律的一类无振荡中心差分格式。Ｈ．Ｎｅｓｙａｈｕ和Ｅ．Ｔａｄｍｏｒ研究了以交错型ＬａｘＦｒｉｅｄｒｉｃｈ格式（ＬｘＦ）为模块的无振荡中心格式的构造与熵不等式。此类格式利用高阶的ＭＵＳＣＬ型插值替代一阶分片常数逼近，减少了ＬａｘＦｒｉｅｄｒｉｃｈ格式的过多数值粘性，建立了一维标量非线性双曲型守恒律的一类高分辨格式。讨论以非交错ＬａｘＦｒｉｅｄｒｉｃｈ格式为模块建立起的差分格式。证明了此格式具有二阶精度、ＴＶＤ性质并在一定条件下满足熵条件。     

Euler方程的动力学通矢量分裂格式

文中基于气体分子动力学理论，研究和发展了一类新型逆风差分格式——动力学通矢量分裂格式（ＫＦＶＳ），并论证了格式的逆风性质和Ｌ１稳定性。最后给出了高精度格式的构造，并就一维激波管问题作了数值试验，还与精确解作了比较，结果令人满意。     

二维非定常化学非平衡流动高精度数值模拟及其与实验的比较

本文成功地将二阶迎风 TVD 数值格式运用到弱电离、化学非平衡斜激波反射流场的数值计算中。气体力学方程与组元质量守恒方程组全耦合求解是目前非平衡流场数值计算中的新方向,文中详细给出了气体力学方程和组元质量守恒方程组全耦合求解过程,最后给出了三种不同情况斜激波反射流场数值结果并进行了详细讨论,本文所得结果与文献〔4〕提供的实验结果和文献〔1〕提供的数值结果相符合。