LU—ADI方法在求解跨声速翼型绕流问题中的应用

本文利用LU-ADI方法求解了跨声速范围内的Euler、NS方程。研究了边界的处理和某些求解方法的改进。     

提高飞行器升力的若干方法研究

本研究用数值模拟方法研究了四种增升方法。数据模拟的出发方程是可压流的Ｎ－Ｓ方程，格式为Ｂｅａｍ－Ｗａｒｍｉｎｇ格式的改进型。通过上述方法分别研究了：（１）动界效应对提高升力的作用；（２）涡的有利干扰对主升力面的影响；（３）吸气的驻涡增效应；（４）升力面下表面向下喷气的增升作用。研究结果表明，采用上述四种方法可有效地提高飞行器的升力，对改善飞行器的性能有着十分重要的意义。     

翼型大攻角绕流的数值模拟

 以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

旋涡对翼型的气动干扰

 用数值模拟方法研究了旋涡对翼型的气动干扰。出发方程是N-S方程。数值格式为Beam-Warming格式。结果表明,旋涡位置改变,其对翼型的气动(激波和升力)干扰随之有显著变化;涡的有利干扰可使翼面上流动分离受到一定程度的控制,导致翼型升力增大;涡的位置和强度的合理组合能有效地提高升力。计算还表明,单一旋涡扫过翼型时平均升力将提高50％左右。     

三角翼大攻角分离流开缝吸气效应研究

 分别采用数值方法和实验方法研究了大后掠角三角翼大攻角条件下,背风面开缝吸气对流场结构和气动力特性所造成的影响。数值模拟采用 Ｈａｒｔｅｎ Ｙｅｅ 的二阶精度隐式 Ｔ Ｖ Ｄ 格式和 Ｎ-Ｓ方程;实验采用激光蒸汽屏流场显示和应变天平测力技术。通过在三角翼背风面开缝吸气,抽掉低能气流,实现增升的效果不理想,同二维吸气涡控制增升效果相差很远,这主要是由于二维同三维分离产生机理上的差异所造成的。     

粘性流绕俯仰振动翼型流动的数值模拟

本文直接从N-S方程出发,利用LU-ADI格式和Baldwin-Lomax代数湍流模型研究了粘性流绕俯仰振动翼型的流动。数值实验表明,本文的数值结果同实验吻合较好。通过数值模拟手段研究了振动对流场中激波和分离这两大主要特征的影响,结果发现;(1)激波滞后于攻角的变化,如当NACA0012翼型在负攻角状态下上翼面存在激波而下翼面无激波。(2)由于翼型的振动,分离被减轻升力相对提高。(3)随翼型攻角和振幅的变化,翼型振动的升力回线走向可不同。     

分布存储结构并行处理机上隐式有限差分算法并行尝试

 通过对传统的Beam-Warming隐式格式的并行化尝试，实现了无粘和粘性空气动力学方程的并行求解。计算中采用了SPMD并行编程模式，并将静态非重叠网格区域分裂技术用于建立计算区域和并行处理系统处理器节点之间的映射。结果表明原等价串行算法的稳定性和收敛性被基本保留，获得的加速比和并行效率较为满意。     

细长体大攻角绕流的数值模拟

 用数值模拟方法研究了细长体大攻角的绕流特性（求解在周向３６０°范围内进行）。数值模拟的出发方程是三维可压流的全Ｎ－Ｓ方程,格式为二阶迎风ＴＶＤ格式。通过数值模拟研究了附近,攻角１９°、３２．５°和５０°状态下,细长体绕流场中的激波、分离流动及表面压力分布等重要特征。     

利用转动前缘的动边界效应控制翼面分离的数值模拟

本文使用数值模拟方法研究了活动翼面对翼型背风面分离的控制作用。数值实验的模型为一头部装有一可旋转圆柱的对称Joukowsky翼型。使用的数值方法为在贴体网格上用有限差分法求解可压流的NS方程。数值实验表明,当头部圆柱以一定角速度旋转时,其升力比不旋转时有明显的增加,分离也受到一定程度的控制。     

大钝体三维粘性流场数值模拟

本文用有限差分法数值求解薄层Ｎ－Ｓ方程，模拟了高超声速绕大钝头倒锥体的三维粘性流场。来流Ｍ∞＝７．０，Ｒｅ∞４．５×１０５，差分格式为二阶迎风ＴＶＤ格式。给出了包括复杂三维分离流动现象在内的流场结构，及其整体物面压力分布。研究了攻角对流场的影响，发现在一定攻角以后，流场结构发生很大的变化。