基于嵌套网格和计入尾迹影响的旋翼流场数值模拟

给出了三维Euler方程数值模拟悬停状态旋翼流场的方法和模型。为充分考虑旋翼尾迹对流场的影响,采用了由动量理论导出的远场边界条件和由尾迹诱导速度代入滑移条件而构造的新翼面边界条件。为减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。针对嵌套网格中关键的贡献单元搜寻问题,本文采用一种伪贡献单元搜寻法。应用上述方法进行了算例计算,给出了有无尾迹修正时的旋翼桨叶表面压力分布和沿展向升力分布的计算结果,并与可得到的实验数据进行了对比。此外,从计算的等涡量线图上分析了旋翼尾迹流场畸变的发展趋势,表明了Euler方程法在旋翼涡尾迹捕捉方面的能力。     

自适应网格局部加密方法在可压缩流场模拟中的应用

本文基于M.J.Berger和J.Oliger发展的自适应网格局部加密算法,它采用多个组合网格的思想,用于求解双曲型方程组。它结合显式有限差分格式,采用Richardson外推技术自动进行局部截断误差估计,对精度低的区域产生新细网格局部加密,或去除无需再加密的旧细网格,以最小的运算量达到给定的精度要求。网格可一层层加密下去,按层次覆盖,每个网格是有任意方向的矩形均匀网格。这套算法独立于求解所用的差分格式,很容易和各种格式结合。 我们从M.J.Berger和J.Oliger程序出发,参考M.J.Berger和P.Colella,对守恒型方程组,实现了粗细网格交接面上数值通量守恒,从而能计算间断解。对网格边界处理、误差估计和网格重新生成等做了一些改进,并且可以灵活地在任何感兴趣区域指定加密。 我们用不定常Euler方程组计算了带前台阶的二维管道中Mach-3流动,采用MacCormack显式附加人工粘性差分格式。共用了四层网格结构:Δ=1/10、1/20、1/80和1/320。其中基本网格是Δ=1/10;到台阶上第一次激波反射区域从Δ=1/20加密至Δ=1/80,计算到定常所需运算量仅是同等均匀网格△=1/80的28％;我们还用Δ=1/320网格对台阶拐角和切向间断处指定加密,结果表明可明显改善台阶上的激波Mach反射,切向间断清晰可辨,所需的计算量是上述均匀网格Δ=1/     

飞行器复杂外形流场NS方程数值模拟

本文采用分块对接和“O”型网格技术，首先在飞行器表面长出“一层”分块结构网格来，再将计算区域用若干分块结构网格填充，这样既保证了网格的质量和模拟复杂外形的能力，又大大提高了计算效率。在计算中，采用有限体积法，MuScL型迎风差分格式求解Ns方程，根据多块网格的排序，研究了计算方法。最后，对某飞行器外形进行了亚跨超流场数值计算，计算结果与试验符合较好。     

一种跟踪电网频率的新方法

基于正交变换,提出了一种新的电网频率跟踪方法.针对两种模型进行仿真,表明该方法能够较为精确与快速地跟踪电网频率,尤其针对突变频率的信号,具有较好的跟踪效果;同时算法克服了同步采样条件的限制,通过迭代,在不同初始频率下均具有较好的收敛性.     

网格生成和数值模拟的讨论

讨论了多种网格生成方法及其数值模拟。给出一种生成多部件复杂外形飞行器分区网格和Ｅｕｌｅｒ方程分区解。讨论了前缘分离涡为主体的旋涡流场。作了跨声速流分区的并行计算的初步尝试。给出一种二维非结构网格的生成方法。     

螺旋桨飞机俯仰力矩特性改进方法

螺旋桨滑流对飞机各部件的气动干扰造成飞机的俯仰力矩特性恶化，危害飞行安全。开展螺旋桨滑流对螺旋桨飞机俯仰力矩特性的影响机理研究和螺旋桨滑流作用下螺旋桨飞机俯仰力矩特性改进方法研究十分重要。采用动态重叠多块结构网格，通过求解非定常雷诺平均可压缩Navier-Stokes方程，数值模拟了某螺旋桨飞机不同拉力系数下降落构型的绕流流场。结果显示螺旋桨滑流对机翼、平尾的气动干扰是导致飞机俯仰力矩特性恶化的主要原因，改进螺旋桨飞机俯仰力矩特性的关键是改进其平尾的升力特性。以拉力系数等于0.4时的降落构型为优化对象，开展了俯仰力矩特性的改进方法研究，包括降低平尾高度、改变平尾上反角、抬高螺旋桨轴线等方法。研究发现在不增大机翼对平尾整体下洗强度的前提下，减小平尾和螺旋桨轴线的垂向距离，可以明显地改善螺旋桨飞机的俯仰力矩特性。     

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.     

一种基于N—S方程的CFD／CSD耦合计算方法

对于气动弹性问题,在跨音速和大变形情况下,模态叠加法失效,为此需要探讨全场求解耦合方程的新方法。利用耦合计算流体动力学和计算结构动力学（CFD／CSD）方法,数值模拟了机翼的气弹响应问题。气动力的求解采用三维Navier—Stokes方程的有限体积法,实现了一种多块结构网格的动网格方法生成气动网格,采用有限元方法计算应...     

基于有限体积格式的自适应笛卡尔网格虚拟单元方法及其应用

基于自适应笛卡尔网格,结合虚拟单元方法处理浸入边界,发展了一种有限体积格式实现模式,数值模拟了二维无粘可压缩流动。本文发展的方法相对于传统的有限差分方法具有执行效率快,易于获得守恒性等优点。最后,利用一些标准算例和多物体流场的计算问题验证了算法的有效性。     

基于动态混合网格的多体分离数值模拟方法

在以往发展的动态混合网格技术的基础上,耦合六自由度运动方程计算和非定常流场计算,建立了一种多体分离问题的非定常数值模拟方法。首先采用刚体运动的六自由度运动方程,建立刚体运动轨迹计算模块,然后以该模块为纽带将以往建立的非定常计算方法与改进的基于Delaunay背景网格插值方法和局部网格重构方法相结合的动态混合网格生成方法耦合起来,形成了精度较高、效率较快的多体分离数值模拟方法。利用该方法对典型多体分离问题进行了数值模拟,计算结果与风洞实验结果吻合较好。