基于风洞试验和数值模拟的超临界机翼雷诺数修正方法研究

结合风洞试验方法和数值模拟,对采用超临界机翼的大型飞机进行雷诺数影响规律研究。对比分析了不同雷诺数下的试验结果和数值模拟结果,在此基础上研究了基于数值模拟结果的雷诺数修正方法,将低雷诺数试验结果向高雷诺数进行修正。修正结果与相应雷诺数试验结果相比,阻力系数相差不超过0.0004,升阻比最大误差约为0.2。针对于力矩系数修正误差问题进行了修正方法改进,改进后的修正误差从0.01降为0.001,表明了修正结果在飞行雷诺数下的适用性。     

固体燃料冲压发动机燃烧效率建模与数值分析

为了获得固体燃料冲压发动机的燃烧效率计算模型,采用理论分析的方法创建了燃烧效率模型,初步推导出燃烧效率的表达式,并运用数值模拟的方法进行了深入研究。研究结果表明,所建数值计算模型准确度较高,与试验结果的相对误差均在3%以内,可用数值模拟代替试验验证研究结果的可信度。燃烧室燃烧效率与进气流量、燃烧室压强燃烧室长度成指数关系,同燃料内径以D=2为分界点成分段指数关系。改变工况后运用拟合公式计算所得结果与数值模拟结果契合度较好。     

基于耦合隐式方法的复杂外形绕流计算研究

飞机复杂外形的气动特性数值模拟在飞机型号研制中有重要的作用,也是数值模拟中比较复杂的问题。笔者采用基于Octree技术的方法生成非结构四面体网格,用耦合隐式离散方法进行有限体积离散,利用此方法对某型飞机带外挂和舵面偏转情况的气动特性进行了数值模拟,将计算结果与实验结果作了比较,收到较好效果。     

反潜导弹分离舱张开特性数值方法研究

反潜导弹弹箭分离时,分离舱张开投放战斗载荷。采用重叠网格方法对分离舱动态张开过程进行非定常数值模拟,并通过火箭橇地面试验验证了数值方法的正确性。基于典型工况非定常数值模拟结果,考虑分离侧滑角、分离舱张开角和张开角速度等因素,通过多维插值获取分离舱实时气动力求解不同分离条件下分离舱张开特性,插值求解结果与直接非定常数值模拟结果吻合良好。这表明,基于非定常数值模拟结果的插值求解方法可以对分离舱张开特性进行快速预报,适用于弹箭分离方案初步设计和评估。     

轴流压缩系统失速可恢复性的预测

基于描述轴流压气机及其系统的动态气动响应的一维非定常单元容积模型 ,通过数值求解并利用可视化的编程方法 ,模拟了单级、三级和九级压缩系统的失速可恢复性。单级压缩系统的数值模拟结果取得了实验结果的验证 ,九级压缩系统的数值模拟结果与该压缩系统在台架上的实际工作情况符合得很好 ,数值模拟结果揭示出了影响压缩系统失速恢复性的一些因素     

基于混合RANS/LES方法与FW-H方程的气动声学计算研究

在气动噪声的数值计算过程中,非定常流动的求解精度对声学计算结果有着重要的影响.以机体噪声计算的标准算例双圆柱绕流气动噪声问题为研究对象,采用基于非线性k-ε湍流模型的限制数值尺度(LNS)方法对双圆柱绕流进行了数值模拟,将计算得到的气动特性和流动特征与相应的试验结果进行了对比分析.为了求解远场观测点处的气动噪声,在精确求解双圆柱绕流流动的基础上结合基于FW-H(Ffowcs Williams-Hawkings)方程的声类比方法进行数值计算,并通过圆柱体的展向相关性对计算结果进行了修正,将得到的最终结果与相应的声学试验结果进行了对比,两者吻合良好,表明该数值方法是准确、可靠的.     

弹射座椅减速性能的数值仿真计算

减速性能是火箭弹射座椅的一个重要性能指标,而且一直不能得到很好的评估预测,只能根据地面弹射试验和经验粗略估计座椅的减速性能。采用数值仿真的方法,在获得弹射座椅的气动参数的基础上求解空间六自由度方程,得到了各气动参数以及座椅的减速性能。方法基于计算流体力学以及matlab的simulink仿真工具箱。最后数值仿真的结果与地面弹射试验的数据进行比较,其平均误差小于5%,说明数值计算的结果是比较准确的。数值仿真的方法在弹射座椅研究中的应用是完全可行的,有广泛的应用前景。     

气泡和液滴在上升或下降时大变形的数值模拟

本文在考虑界面张力的条件下用涡片方法导出了二维无粘气泡和液滴在另一种静止无粘流体中上升或下降时运动和变形的支配方程，提出了求解支配方程的数值方法，并在数值模拟的例子中研究了空气泡和煤油滴在水中上升时及水滴在煤油中下降时的变形规律，研究了数值模拟中界面张力对界面形状的影响，气泡上升的数值模拟结果与相应的试验结果基本一致。数值模拟结果表明本方法可正确模拟饣泡和液滴在上升或下降过程中的大变形，直到界面发     

基于浸入式边界法的振荡转子叶片数值模拟

对于流固耦合问题,传统的数值模拟方法由于需要不断的网格重构使得问题变得复杂,本文通过利用浸入式边界法建立了在单一坐标系和网格下求解振荡转子叶片的快速计算模型,从而避免了传统方法中由于需要不断地重构贴体网格造成的数值模拟的复杂性。为了验证该方法的正确性,分别对低KC(Keulegan-Carpenter)数下的振荡圆柱和两自由度振荡圆柱进行了数值模拟,计算结果与以往的试验结果和数值结果吻合得很好,证明了浸入式边界方法的可靠性。在此基础上,对振荡转子叶片进行了数值模拟。结果表明折合速度和叶栅稠度是影响振荡转子叶片的重要因素。并且,值得注意的是这种耦合过程并没有生成贴体网格,因此大大减小了计算时间,这样可以更加快速、准确地模拟真实的叶栅流动情况。     

基于能量法和特征值法的颤振预测数值方法研究

为研究适用于工程设计的颤振数值预测方法,采用能量法和特征值法对出现颤振现象的真实叶片进行预测。通过将2种数值方法的预测结果与试验结果进行对比,验证数值预测方法的精度和工程适用性。结果表明:2种方法都较准确的判断了某典型风扇转子叶片的气弹稳定性,与试验结果基本吻合,具有较高的工程实用价值。     

高超声速球锥粘性激波层压力分裂法数值求解

本文将压力分裂法引入细长钝体高超声速粘性激波层数值求解中,较好地解决了这种流场远后身部分的推进求解问题。针对非解析球锥体绕流数值求解的困难,作者提出了实用型初始激波形状计算方法,达到了简化数值求解步骤,提高数值计算效率的目的,数值求解结果与实验结果相符合。     

计算跨音速三维机翼非定常载荷的核函数法

 本文给出了在跨音速流动中,计算三维机翼非定常载荷分布的核函数法和计算结果。此方法是基于线化升力面理论,利用在亚、超音速情况下,联系载荷分布与下洗的积分方程的核函数,采用局部线化法来处理跨音速的非定常绕流问题。在计算中引进了线分布的冲波偶极子,并嵌入适当的正冲波边界条件。 我们从统一的超、亚音速核函数表达式出发,对亚音速区数值积分作了新的简便的处理,对于亚、超音速区相互的诱导下洗给出了自己的数值积分过程。我们算了三个例题,和文献上的数值符合较好。     

跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

用弱非线性理论研究二次稳定性问题

本文把研究三波共振的弱非线性理论应用于平板边界层和平面Poiseuille流三维扰动的二次稳定性问题。其数值结果与文献[1]、[2]用Floquet形式的二次稳定性理论得出的结论是一致的。     

超音速气流绕局部透气凹角流动的数值模拟

以具有压力分裂形式的简化Ｎ Ｓ方程为控制方程,数值模拟了超音速来流条件下的激波 边界层干扰被动控制（ｐａｓｓｉｖｅｃｏｎｔｒｏｌｏｆｓｈｏｃｋ ｂｏｕｎｄａｒｙｌａｙｅｒｉｎｔｅｒａｃｔｉｏｎ）。模拟是以预先给定激波前吹气和激波后吸气的流量来实现的。为了定性地确定吹气或吸气对激波 边界层干扰的影响,首先计算了单独吹气和单独吸气两种情况。数值计算时采用了多重扫描法对控制方程差分离散,以反映亚音速区压力对流场的椭圆性影响。     

超声速主流中横向喷流场的激波—旋涡结构的数值模拟

本文利用NND格式,通过求解NS方程,对二维超声速主流中横向喷流干扰流场进行了数值模拟,计算清楚地给出了激波结构、回流区和混合层。本文计算得到的激波结构和实验相当一致。最有兴趣的是由于喷流的干扰,主流在喷口前发生主涡分叉,观察到三个流向旋转涡和两个反流向旋转涡;在喷口后的背风区,存在具有低压和回流区的尾迹。     

不同壁温及差分格式下超燃冲压发动机的仿真

针对圆截面超燃冲压发动机相关实验进行数值仿真,采用了不同的壁面温度条件和差分格式,并探究其对仿真结果的影响规律。结果表明采用常用的绝热壁面边界条件会导致激波边界层交互位置处的回流区变大,进而导致流道收缩,壁面压力升高。采用高精度格式不仅不能更为精确地计算回流区的大小,反而会使回流区面积变得更大。通过改变壁面边界条件,采用300 K的等温壁面,可以使激波边界层交互位置处的回流区变得极小,使仿真结果与实验结果十分吻合。计算结果表明,针对当前构型与实验条件,盲目提高数值格式的精度并不一定会得到更好的结果,相反可能会使仿真与实验更加偏离,而适当修改壁面温度条件,即使数值格式精度较低,仍然可以得到很好的仿真结果。最后,针对与实验结果吻合的仿真结果,分析了圆截面超燃冲压发动机的流场特性,重点研究了凹腔剪切层及其质量交换特性。     

NND格式在航天飞机头部段N—S方程求解中的应用

本文采用实质上为二阶精度的无波动、无自由参数的耗散差分格式(NND格式)求解了非定常的二维和三维完全N-S方程,对类似于航天飞机头部外形的超声速粘性气体绕流流场进行了模拟。在计算中,使用了代数网格生成技术,对脱体激波采用了装配法。作为验证算例的二维圆柱绕流计算结果同文献[4]的数据进行了比较。二维凹陷外形绕流的计算结果同无粘流结果进行了对比,由于粘性流动出现了分离,二者的壁面压力是有差异的。最后,对三维外形的绕流也进行了模拟,计算结果提供了流场的细节。     

时间推进积分方程法求解非定常跨声速流动

本文提出了一种求解非定常跨声速流动的新方法——时间推进积分方程法,此法克服了时间线化积分方程法的限制,能较好地模拟激波的运动。本文首先用一维波(?)问题——模型问题阐明此法的基本思想,然后将它应用于二维非定常跨声速流动中。本文还首次引入拟速度位的概念,使时间推进积分方程式得到简化,尾涡条件和Kutta条件更易处理。数值计算表明时间推进积分方程法是合理可靠的。     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...