球锥俯仰阻尼导数的数值计算

本文从非定常Ｅｕｌｅｒ方程出发，应用谐振摄动法计算了球锥体超声速、高超声速下的俯仰阻尼导数和。Ｅｕｌｅｒ方程的求解利用ＳＣＭ方法和ＭａｃＣｏｒｍａｃｋ的两步算法。文中用任意攻角下不同俯仰轴的气动导数转换公式来简化计算，用较粗网格可以得到较精确的结果，程序编制相对简单。一次计算可以得到定常绕流解及一系列钝头率、不同俯仰轴下的多种动导数结果，在ＥＬＸＳＩ６４００机上约需２０～４０分钟。所得结果与试验结果符合程度令人满意。     

小迎角轴对称体跨声速绕流的一种Euler解法

应用小迎角展开的方法将小迎角下旋成体绕流Ｅｕｌｅｒ方程进行摄动展开，对零介和一阶方程分别采用隐式ＮＮＤ格式进行了计算，与实验进行比较得到了较为满意的结果。该方法采用二维的计算得到了三维的结果，大大提高了计算的效率。     

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

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

跨音速压气机动静叶排相干的三维非定常流动数值分析

采用时空二阶精度ＮＮＤ显式差分格式求解三维Ｅｕｌｅｒ方程，对跨音速压气机级中动静叶排相干形成的非定常流动进行了数值分析。除动静叶排流动的控制方程体系、数值方法外，着重讨论了动静叶排计算域、网格以及边界条件的处理。对某压气机第１级动静叶排相干的非定常流场的计算，给出了叶片表面的马赫数分布、压力波动幅值，以及叶片槽道中瞬态密度等值线等，并分析了结果的合理性。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

再入飞行器高超声速绕流的一步和二步有限元计算

本文针对再入飞行器的空气动力问题，采用一步和二步有限元格式求解可缩守恒型Ｅｕｌｅｒ方程，并结合不同的人工粘性模型来捕捉激波，给出了高超声速二维钝头绕流的结果。     

可压缩球涡和平面激波相互作用的数值研究

行进了可压缩球涡和平面激波相互作用的数值研究。一方面，用非定常Ｅｕｌｅｒ方程数值模拟了一种可压缩球涡的自身演化，同时，利用Ｒａｎｋｉｎｅ－Ｈｕｇｏｎｉｏｔ关系，在流场中嵌入运动激波，求解了激波－球涡相干的流动过程。     

可压缩涡环和平面激波相互作用的数值研究

本文从非定常的Ｅｕｌｅｒ方程出发,进行了可压缩涡环和平面激波相互作用的数值研究。     

M∞=1附近轴对称Euler方程的跨声速计算

本文综合应用近似因子分裂法，矢通量分裂法及ＮＮＤ格式的离散思想，对轴对称Ｅｕｌｅｒ方程构造了一种修正的隐式格式，对球柱及多截锥旋成体在Ｍ∞＝１附近跨声速绕流情况进行了计算，与其它方法比较，本方法具有收敛速度快、精度高的优点。尤其对Ｍ∞＝１时的跨声速绕流情况，上述优点就显得更加重要。     

带副翼三维机翼绕流的Euler方程解

提出了一种求解带副翼偏转三维机翼绕的Ｅｕｌｅｒ方程计算方法。采用对接分区网格民分区求解算法的结合,有铲地求解绕引外形的复杂流动,提出了一种满足通量守恒的内边界耦吉选用ＶａｎＬｅｅｒ分裂格式离散玩粘通量项,并构造了一种Ｌｉｍｉｔｅｒ函数以保证ＴＶＤ性质。数值算例表明本方法是求解带操纵纵面偏转的机翼绕流的有效方法。     

振动翼型跨音速绕流的数值模拟

从积分形式的非定常Euler方程出发,在固联于振动冀型的贴体坐标系下,采用有限体积法进行离散,并根据对称型TVD格式构造相应的数值通量,数值模拟了振动翼型的非定常跨音速绕流的变化过程。计算结果表明,该方法能较好地捕捉流场中的运动激波,准确地定出激波位置与强度;并且易于推广到求解任意振动方式的跨音速非定常绕流问题。     

考虑地面效应的翼型动态特性数值模拟

为了研究考虑地面干扰的飞行器非定常气动特性,基于滑移网格技术,通过求解二维非定常Euler方程,对不同近地高度下的NACA0012翼型的动态振动特性进行了数值模拟。分析了存在地面干扰时的升力系数以及俯仰力矩系数迟滞环变化,并结合强迫振动法,计算了翼型各个近地高度下的纵向组合动导数。结果表明,地面效应不仅对定常流场产生影响,更显著地影响了非定常气动力及力矩,近地高度越小,这种效应越明显,使得升力系数的迟滞环面积变小,幅值减小,而俯仰力矩系数的迟滞环变化不规律,力矩系数的动导数随着近地高度减小而增大(绝对值减小),纵向的动态稳定性产生损失。地面效应干扰在飞行器非定常气动研究中应该引起重视。     

超声速平面钝体流动的快速解法

本文介绍了一种平面流动的快速欧拉方程解法。该方法将原参数非定常欧拉方程组重新组合成以广义称曼变量表示的欧拉方程组，再用二点二步迎风格式离散解。针对钝体流动，本文先建立了动网格下的方程，构造了动网格下的算法。提出了一种简单的激波处理方法。计算结果表明，该方法速度快，稳定性好，对初场不敏感。     

叶片通道内二次涡的数值模拟

采用三维Ｅｕｌｅｒ方程数值模拟叶片通道内流向二次涡,及其对叶片出口气流角的影响。算例表明：采用这种方法可以计算出流向二次涡;该涡造成的出口气流偏转角与实验值比较,在远离壁面处相差较小。此方法可用于叶栅通道内流向二次涡的初步分析     

钝体尾迹的稳定性研究及流动控制探讨

系统地回顾了研究钝尾尾迹稳定性的主要方法，并利用局部绝对／对流稳定性分析及整体稳定性分析两种方法，研究了圆柱绕流定常和时均尾迹流场的稳定性。采用低维Ｇａｌｅｒｋｉｎ方法研究了尾迹对圆柱流向振动的响应问题，揭示了频率锁定，涡街的激发和抑制等现象。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

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

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

跨声速无粘流绕振动翼型流动的数值模拟

本文从Euler方程出发,利用LU-ADI格式对绕二元振动翼型跨声速流场进行了数值模拟。文章中给出了绕NACA0012和NACA 64A 010翼型不同状态的几组数值结果,在同某些实验和其它文献结果基本吻合的前提下,研究了振动对翼型表面压力、激波位置的非定常效应以及其它振动参数对流场非定常效应的影响。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.