亚、跨音速时尖头旋成体在大迎角时的法向力特性

本文对一组尖头旋成体在大迎角时的法向力特性进行了实验研究,并对实验结果作了分析,研究表明,目前常用的以横流理论为基础的计算方法未能合理地考虑到粘性和压缩性的影响,因此不能反映出流动的真实情况。 此外,对不同安装方位角时法向力系数的实验研究表明,法向力系数和背风面上的旋涡谱有着密切的关系。由于大迎角时尖头旋成体背风面上的旋涡谱难于预计,由此可以推断,尖头旋成体在大迎角时的气动特性不应该是唯一的,而可在一定范围内变化。     

细长旋成体背风面非对称流场的DES模拟

利用脱体涡模拟DES方法对细长旋成体非对称绕流进行数值模拟,比较了基于S-A湍流模型的DES与RANS方法对背风面非对称分离涡的数值模拟能力,发现在大迎角非对称涡未破裂情况下,DES与RANS均能模拟出与实验相符的非对称绕流;当迎角增大,背风面非对称涡发生破裂时,RANS无法准确捕捉到背风面流场的非定常性,而DES能准确预测非对称涡的飘起与破裂,并与实验值接近。计算结果表明:与RANS方法相比,DES方法具有更好的模拟大攻角,大分离流动的能力,尤其在非对称流场的分离涡破裂模拟方面具有明显优势,能够更真实模拟出细长旋成体背风面分离涡破裂之后的非定常流动特征。     

细长体空间涡系的观测

本文用流态显示技术和空间总压测量,对于钝头和尖拱形头部两个细长旋成体模型在中等和大迎角状态的复杂背涡系进行了实验研究。侧重于“二次分离区”的流动观察和测量。实验表明:在中等迎角时细长体背风面上存在有一对二次涡,它们的旋转方向相反;靠外侧的一个涡尺度较大,其旋转方向与同侧主涡相同,当同侧主涡破裂时,它也发生破裂。随迎角进一步增加此涡呈现出从物面上间断地形成和脱落。     

机身前体横截面形状对三角翼涡破裂位置的影响

通过求解三维定常雷诺时均方程,采用剪切应力输运(SST)湍流模型,在亚声速范围内,分别对融合体前体-三角翼组合体和旋成体前体-三角翼组合体流场中翼涡破裂现象进行了数值模拟.模拟结果显示:与旋成体前体相比,在中、大迎角时,融合体前体的分离涡,涡量集中、强度高,进入机翼上方流场后,能够与翼涡密切耦合,彻底地改变翼涡强度的分...     

低速大迎角三角翼旋涡流动的有粘与无粘数值模拟

通过人工压缩性方法求解了三维定常不可压Euler方程和雷诺平均Nayierstokes方程，应用Beam－Warming近似因子分解格式及其对角化形式，采用BaldwinLomax代数湍流模型。计算模型为70°三角翼，分析了有粘与无粘流动旋涡特性的影响，尤其是对Euler方程模拟大迎角分离涡、涡破裂的能力及其局限性提供了数值验证，计算结果与实验结果相一致。     

旋成体轴对称跨声速全位势流的高效差分算法

本文用守恒型全位势方程,贴体坐标网格,对旋成体轴对称跨声速绕流的差分数值计算方法进行了研究;根据最佳收敛准则,提出了轴对称情形的AF2迭代算法,并将此算法应用于半球头柱体、弹体等各种外形的旋成体。与一般方法仅适用于亚声速自由流不同,本计算可从亚声速、跨声速自由流一直到低超声速自由流。计算结果表明,本文方法收敛快,与实验及其它方法的结果符合较好。     

基于RANS-LES混合方法的翼型大迎角非定常分离流动研究

使用雷诺平均Navier-Stokes方程-大涡模拟（RANS-LES）混合方法中的延迟分离涡模拟（DDES）方法,模拟了NACA 0015翼型在大迎角下的静态绕流和强迫振荡运动并与实验值进行了比较。在大迎角静态翼型大分离流动模拟中,DDES方法捕获了非定常RANS计算未能获得的机翼背风面的涡脱落现象。在所采用的RANS和DDES模型中,基于剪切应力输运（SST）湍流模型的SST-DDES混合方法给出的时均压力系数分布与实验吻合得最好。在大迎角强迫振荡翼型绕流模拟中,DDES方法得到的非定常气动载荷与实验值吻合得很好,正确地反映了最大迎角处阻力和俯仰力矩的阶跃性突变,而非定常RANS计算则给出了完全错误的趋势。     

旋成体大迎角侧向力的控制方法计算研究

提出了运用充气气囊控制旋成体大迎角非对称分离涡的思路,并通过求解雷诺平均Navier-Stokes方程模拟了不同形状充气气囊作用下的流场,得到了充气气囊形状对旋成体侧向力的影响规律,计算结果表明充气气囊完全可以形成期望的流动干扰,为飞行器大迎角侧向力的控制探索了一条新的途径。     

细长旋成体大迎角俯仰振荡的数值研究

为研究俯仰振荡对细长旋成体流场与气动特性的影响,利用雷诺平均方法与大涡模拟的混合方法相结合的脱体涡模拟方法对大迎角下细长旋成体俯仰运动进行了数值模拟研究。通过对背风面分离涡的强度和位置、细长体绕流形态、截面压强系数和侧向力系数的观察和分析,发现施加小振幅、特定频率的俯仰振荡对大迎角细长旋成体背风面流场有明显的控制能力,...     

垂直于物面的横截面上流态的拓扑

对于旋成体有迎角绕流,本文研究了垂直于物面的横截面上的流态结果。利用解析方法,给出了物面轮廓线上、迎风和背风对称线上以及全横截面上截面流线方程奇点的分布规律以及奇点沿流动方向的演变规律,从而可定性地确定截面流态的结构。通过求解ＮＳ方程数值模拟了钝锥有迎角高超声速绕流流场,计算获得的垂直于物面的横截面流线的拓扑结果和理论分析完全一致。     

大迎角体涡,翼涡干扰的Euler方程数值模拟

本文用Ｅｕｌｅｒ方程的简化区算法计算翼－身组合体绕流。当大迎角时，计算能自动捕捉到大后掠前缘的分离涡，机射的分离涡由强加Ｋｕｔｔａ条件生成，形成体涡与翼涡的相互干拢。本文对加长前机身的ＡＧＡＲＤ－Ｂ翼－射组合体模型，给出机射有分离和无分离时，横流截面的速度矢量分布、法向力和压力中心，分析体涡的影响。     

用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。     

用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。     

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

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

绕跨声速三角翼的激波/涡干扰流场数值模拟

 在绕三角翼的跨声速流动中,随着迎角的增加,三角翼上的涡破裂位置会出现突然前移的现象。针对这一与亚声速下不同的流动现象,采用带曲率修正的Spalart-Allmaras(SAR)湍流模型,求解定常雷诺平均Navier-Stokes(RANS)方程,对不同迎角下绕65°后掠尖前缘三角翼的跨声速流动进行数值模拟,并在此基础上,采用基于SAR湍流模型的脱体涡模拟(DES)方法,对由激波干扰导致的前缘涡破裂位置的运动规律进行了初步探讨。模拟结果与试验结果对比表明:SAR湍流模型能准确地模拟出三角翼上的激波系统和旋涡结构,并能准确模拟出由于激波干扰导致的涡破裂位置突然前移的现象。此外,对涡破裂后流场的非定常数值研究发现,支架前端正激波的干扰作用使得涡破裂位置向下游移动比较突然,而向上游移动则相对缓慢。     

再人体复杂外形无粘流场数值计算研究

本文在文[1,2,3]的基础上研究了非定常Euler方程的推进迭代方法,并将无波动、无自由参数的耗散差分格式(下称NND格式)发展为隐式、迎风格式,用以计算有差曲控制翼的再入复杂飞行器的育攻角、有侧滑角的超声速无粘绕流流场。其结果是准确的。本文的方法具有以下特点: 1.可以在同一计算程序中实现亚跨超流场的计算:在超声速区域,采用推进技术;在亚跨声速区,则自动选择推进迭代求解技术。 2.NND格式不仅具有公式简单,无自由参数的特点,而且过激波时满足熵条件,不会产生非物理解,也不发生振荡。 3.隐式差分格式的无条件稳定和Gauss-Seidel迭代的快速收敛,计算省时。 本文的方法曾用于烧蚀头部凹陷外形和航天飞机简化外形的计算,能自动选择迭代和推进区域。可以预言,如计算机条件许可,本方法可以用于大攻角的无粘流场计算。     

旋翼桨-涡干扰气动特性计算及参数影响研究

建立了一个适用于旋翼桨-涡干扰(BVI)气动特性分析的数值方法。在该方法中,控制方程采用惯性坐标系下的非定常Euler方程,以适合于分析BVI流场的特点;为便于前飞流场分区求解和信息传递,使用了嵌套网格方法;对于流场中涡线的模拟,建立了一种适用于有限体积格式的涡引入方法——广义网格速度法,以简化计算。应用上述方法对旋翼BVI流场进行了计算,并与可得到的试验数据进行对比,验证了方法的有效性。通过对比桨叶弦向位置压强的变化得出,在旋翼BVI过程中,气动载荷主要来自桨叶10%弦长内的前缘部分的压强突变。文中还进一步分析了涡强、干扰距离和干扰角度对BVI气动特性的影响。本文的BVI数值计算表明:当涡接近桨叶前缘时,升力达到最大;而涡位于后缘位置时,诱导速度改变桨叶环量分布,破坏了后缘的库塔条件,但随着涡的远离,桨叶表面环量开始重构。     

超声速条件下内埋式武器分离特性的数值分析

采用非结构动网格技术,耦合求解欧拉方程和六自由度(6-DOF)运动方程,对内埋式导弹在亚声速和超声速条件下的分离轨迹和姿态变化进行了比较,分析了不同马赫数下气动力对导弹分离特性的影响,重点讨论了超声速条件下导弹分离过程中相对于载机的位置和导弹的姿态变化情况,分析了超声速飞行环境下导弹安全分离的条件.     

超音速多圆弧(MCA)叶型叶珊的研究

本文提出了含脱体激波的叶栅有旋超音速流场的特征线起始线计算方法。采用激波嵌入法解完全的欧拉方程,计算了含脱体激波的多圆弧叶栅的超音速进口流场。本文发展了一种将超音速流场计算与亚音速扩散损失计算的尾迹法相关联的计算超音速叶栅损失的方法,可分别确定激波损失和亚音速扩散损失。      

变来流速度下旋翼翼型非定常气动特性分析

为分析变来流速度状态下的旋翼翼型气动特性,提出了利用翼型平移来模拟来流速度变化的数值方法.在此方法基础上,采用基于隐式LU-SGS(lower upper symmetric Gauss-Seidal)方法的非定常雷诺平均N-S(Navier-Stokes)(RANS)方程,模拟了SC1095旋翼翼型在定迎角 变来流速度及变迎角 变来流速度状态下的非定常气动特性.通过对比分析发现:翼型在变速度-定迎角状态下会表现出明显的非定常现象,产生了前缘分离涡,气动特性会出现明显的迟滞效应及波动现象,脉动速度越大,非定常效果越明显.并且基准速度越大,翼型气动特性的峰值越大;翼型迎角越大,非定常涡出现的也越早.考虑直升机旋翼翼型实际工作环境,在变速度-动态失速状态下,翼型最大迎角处的气动力会得到一定程度的削弱,在小迎角下的气动力得到一定程度的增强,且脉动速度越大,翼型的非定常特性也越强.