前缘对进气道性能影响的数值模拟

应用质量平均的N-S方程和B-L代数湍流模型，在不同压缩面前缘半径条件下，以高超声速进气道的流场进行了数值模拟，对比研究了进气道压前缘半径对流场结构及性能的影响，在计算过程中，对方程中的对流项采用了空间为二阶精度的TVD格式，扩散项则采用了二阶中心差分离散。     

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

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

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

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

翼型扑翼运动流场的数值模拟

基于非定常N-S方程,利用动网格方法和预处理方法对翼型在低速度下扑翼运动进行数值模拟。采用Jameson中心格式的有限体积法、双时间推进法求解非定常N-S方程;紊流模型采用基于SA紊流模型的DES模型。本文计算模拟翼型在低雷诺数下扑翼运动过程,计算所得推力系数结果与实验参考结果吻合较好,对扑翼飞行器设计有一定参考价值。     

机翼弹性变形对气动特性的影响

在非结构混合网格框架下,采用有限体积法求解了Euler方程和N-S方程。使用结构影响系数法计算机翼结构的弹性变形,通过动态网格技术实现了气动力与结构变形之间的耦合运算,计算和分析了在气动载荷作用下的机翼变形和此时的气动特性,算例分析表明：典型后掠机翼的弹性变形会使机翼产生了正的挠度和负的剖面扭转角,在弯曲和扭转的相互作用下,机翼的有效迎角减小,从而导致机翼的升力系数减小。在Euler方程和N-S方程算例中,考虑弹性变形条件下的机翼升力系数比刚性机翼的升力系数分别减小了25.3%和24.1%。     

翼型和机翼阵风响应的数值模拟

采用有限体积法来数值求解非定常Euler和NS方程,并在非定常计算中通过引入"网格速度"来计入阵风的影响,实现了对刚性翼型和机翼阵风响应的有效数值模拟。首先采用该方法来对翼型攻角突然增大的阵风、陡边型阵风这两种典型阵风情况下的翼型气动力响应进行了计算,结果表明:计算结果与现有的一些理论结果都符合较好,与文献中计算结果也吻合,NS方程的计算结果与Euler方程计算结果在变化趋势上是一致的,但具体的数值有差别。最后针对攻角突然增大的阵风,对不同马赫数下不同展弦比平直机翼的阵风响应进行了数值模拟,得出的计算结果与文献中的计算结果一致。     

蛇形进气道影响下发动机进口部件水撞击特性的数值研究

为了了解亚声速蛇形进气道影响下航空发动机进口段迎风表面处的水滴撞击特性,采用欧拉/欧拉法对蛇形进气道和发动机进口段在来流马赫数为0.3及环境温度为266.5K的大气结冰条件下进行了空气/过冷水滴两相流场的计算,再由已获得的水滴流场进行支板表面水撞击特性的计算。整流罩表面水滴的撞击计算时采用插值的方法将原两相流场的计算结果插值到新生成的整流罩的近壁网格处,之后采用新生成的整流罩面网格进行水撞击特性的计算。通过计算获得了发动机进口结冰参数、支板以及整流罩表面水收集系数的分布。计算结果表明整流罩和支板表面水滴撞击的主要区域与进气道的弯曲方向有关且位于其下表面一侧,最大水收集系数达到0.98。与之相对比的整流罩和支板其他区域的水收集系数则显得非常微小,在数值上约等于0。此外,对于与进气道对称面平行的竖直支板,其与机匣相交的区域存在水收集系数为0遮蔽区。     

高速进气道抛罩过程内流场演化特性研究

为探究进气道整流罩打开过程中内流场的演化规律,针对下颌式进气道及其配套整流罩,基于重叠网格法求解N-S方程,详细分析了抛罩过程中内流场演化过程与特征流场结构。研究表明,堵罩时,整流罩前端发展出V形激波,整流罩尾部形成后向台阶且产生一对涡矢量方向相反的流向涡,开罩过程中,整流罩前端的V形激波与尾部的后向台阶均会对内流场流动产生影响,捕获流量在中间某时刻达到峰值且大于稳定状态。对于开罩过程流动准定常且内流道无大分离区的流场,进气道在V形激波打在唇口时完成起动且整流罩在此后不影响内流场流动。假设的匀加速转动与重构出的真实转动吻合较好,对于更精确的运动过程模拟可采用三次函数的角速度变化律。     

跨声速机翼操纵面非定常气动力Euler方程计算

采用代数法与椭圆型方程优化相结合的结构生成方法为带操纵面F5机翼生成面对接块结构网格;弹性变形技术生成动态网格;有限体积法、双时间推进法求解非定常Euler方程.给出一种适于面对接块结构网格计算的守恒量插值方法.数值计算了F5机翼操纵面作简谐运动时的非定常气动力,计算值与试验值吻合较好.     

超音进气道流场的数值计算

用Jameson的有限体积法,对某型号飞机在不同飞行状态下的进气道流场进行了数值模拟。求解的控制方程为三维非定常N-S方程。在对流项的计算上用二阶中心差分有限体积法,粘性项的计算则采用作者在［1］中发表的一种有限体积框架下的新的离散化技术,时间方向用多步Runge-Kutta方法进行推进,以增大CFL条件数。湍流模型则采用Baldwin-Lomax代数模型。求解过程中还采用了当地时间步长,隐式残量平均、人工粘性等加速收敛技巧。最后给出了各种飞行状态下进气逼出口截面的平均总压恢复与总压畸变,得出了比较满意的结果。     

任意回转面亚、跨音流场计算方法

本文所设计的任意回转面亚、跨音流场的计算软件 ,能根据叶型和进出口条件分析流动特性 ,并自动选择精确和节省机时的合适的数值解法。用户在给出叶型坐标和定解条件后 ,可以得到叶面马赫数分布或压力分布 ,全流场的等马赫数或等压图     

风扇进气畸变三维非定常数值模拟技术研究

采用三维欧拉方程加体积力模型的方法, 发展了一套三维非定常进气畸变数值模拟程序, 以实现风扇/压气机在进气畸变条件下的全环面三维非定常数值模拟.利用该程序模拟计算了某航空发动机风扇进口稳态总压畸变条件下特性和流场结构, 并与试验结果进行了对比, 结果表明该程序用较小的计算机资源即可以实现对风扇全环面三维非定场畸变流场的模拟, 能够以较高的精度预估气流畸变对压气机/风扇特性的影响.      

黏性彻体力建模方法及其应用

研究了将Navier-Stokes方程中黏性项简化为黏性彻体力的建模方法,并将该模型加入已有的计算流体力学(CFD)程序中,实现了用三维欧拉方程加黏性彻体力模型快速模拟叶轮机三维非定常流动的方法.使用该方法对进气周向总压畸变条件下的某双级压气机非定常流场进行了模拟,并与雷诺平均的Navier-Stokes(RANS)程...     

复合振动翼型跨声速非定常流计算

在文献[1]、[2]的基础上,本文把绕翼弦中点作旋转振动的翼型跨声速流计算推广到翼型作水平、垂直方向的平移振动以及平移与旋转振动的叠加。本方法主要特点是网格与翼型一起运动,对无穷远及翼面边界条件作了精确的处理以及计算网格数很少。本文结果对进一步研究旋翼、叶轮转子等的非定常气动力将有重要的意义。     

熵变量欧拉方程组跨超声速流SUPG有限元数值解

通过变量代换将拟线性守恒变量欧拉方程组转换为具有对称性的熵变量欧拉方程组，对该方程组构造SUPG的弱解式，其中权函数摄动项中包含流向算子和间断捕获算子。对翼型的跨超声速流动算例表明，本法具有较好的稳定性和能有效地抑制激波两侧振荡。应用带预处理的GMRES迭代法求解非对称方程组能显著地提高收敛速度。     

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

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

低于自起动马赫数时高超进气道的非定常流动特性

采用非定常数值仿真的方法研究了低于自起动马赫数时高超声速进气道的非定常流动特性.研究表明:低于进气道自起动马赫数时,进气道处于不起动状态,流场发生喉道壅塞性振荡现象,其流场振荡频率为250Hz.流场振荡主要发生在喉道之前,对其后流场影响相对较小,扰动信号由喉道以当地气流速度向下游传播.隔离段长度对喉道壅塞性流场振荡几乎没有影响.飞行马赫数较小时流场未出现振荡现象,当飞行马赫数靠近自起动马赫数时流场出现周期性振荡现象,并且随着飞行马赫数的增大,此类流场振荡趋于强烈;进气道压差阻力随着时间推进呈现周期性变化,振荡频率同样为250Hz.      

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

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

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

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

Numerical Simulation of Unsteady Flow Around Forward Flight Helicopter with Coaxial Rotors

Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of the two coaxial rotors both become worse because of the aerodynamic interaction between them, and the influence of the top rotor on the bottom rotor is greater than that of the bottom rotor on the top rotor. The downwash velocity at the bottom rotor plane is much larger than that at the top rotor plane, and the downwash velocity at the top rotor plane is a little larger than that at an individual rotor plane. The downwash velocity and thrust coefficient both become larger when the collective angle of blades is added. When the spacing between the two coaxial rotors increases, the thrust coefficient of the top rotor increases, but the total thrust coefficient reduces a little, because the decrease of the bottom rotor thrust coefficient is larger than the increase of the top rotor thrust coefficient.