跨音速非定常二维翼型绕流计算

引言 为了计算跨音速非定常气动载荷和了解其流场的性质,需要对翼型绕流进行数值模拟。在构造一个跨音速非定常流动的计算方法时,重要的是要考虑能够处理流场中不断运动的激波;其次要能够处理各种类型的翼型运动,并具有高的计算效率。     

包括粘性影响的跨音速纵向大扰动非定常翼型绕流计算

本文提出一种计算效率高、并改进小扰动理论的二维跨音速定常和非定常流的计算方法——非定常纵向大扰动流速势方程和边界条件的数值解。本方法还考虑了包括边界层位移厚度以及激波-边界层干扰的粘性影响。文中给出了NACA 0012翼型和NLR 7301超临界翼型绕流的算例,计算结果与实验作了比较。     

跨声速无粘流绕振动翼型流动的数值模拟

本文从Euler方程出发,利用LU-ADI格式对绕二元振动翼型跨声速流场进行了数值模拟。文章中给出了绕NACA0012和NACA 64A 010翼型不同状态的几组数值结果,在同某些实验和其它文献结果基本吻合的前提下,研究了振动对翼型表面压力、激波位置的非定常效应以及其它振动参数对流场非定常效应的影响。     

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

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

二维亚、跨音速数字解的流线座标系方程组

本文导出了适用于亚、跨音速二维绕流问题数字解的流线座标系气体动力学方程组。由于在跨音速流场内可能存在激波,所以也导出了流线座标系二维激波关系式。最后附有算例。     

用单调相关格式计算跨音速定常和非定常翼型绕流

在跨音速非定常势流计算方面,人们一直在努力寻求处理包括激波运动的非线性问题的能力和计算效率的提高。Ballhaus和Goorjian（1977年）用交替方向隐式差分法（ADI）成功地直接求解了时间域的非定常跨音速小扰动方程（此法称为“LTRAN2”）。文献[5]引用坐标变换也做了类似LTRAN2的工作。     

欧拉方程组的隐式近似因子分解差分方法的研究

本文对求解二维欧拉方程组的隐式近似因子分解差分方法进行了简要的分析,并利用它计算了简单直扩张通道内的亚音速、超音速和跨音速管流、绕无限长拐角和对称棱形翼型的超音速绕流以及绕NACA 0012翼型的亚音速、跨音速和超音速绕流等典型定常流动。计算结果表明,它具有良好的适应性和精度,可以正确地捕获强激波和自动满足Kutta条件,Courant数可以取大于1的数值,从而加速收敛。分析和计算都是在贴体坐标下进行的。方法本身可以推广到三维流的计算。     

激波绕流与激波碰撞的数值模拟

本文应用基于Poe平均构造的最小耗散的扩散与反扩散守恒差分格式，数值模拟单激波绕流三我形障碍物，及等强度激波碰撞和不等强度激波碰撞的二维非定常复杂流场。计算结果表明，该格式具有很高的单波精度及较好的激波与接触间断分辨度，可直接用于复杂流场的数值模拟。     

计算绕薄翼型跨音速非定常流的积分方程法

 本文从小扰动方程出发推导出绕薄翼型的跨音速非定常流动的积分方程,为使此方程适用于具有激波的流场引入人工粘性,并对其数值求解。通过算例讨论了人工粘性和计算域大小对计算结果的影响,与实验及其它数值计算结果比较,表明本方法的准确度令人满意。计算量小,收敛性好是此法的特点。     

LTRAN2的一种改进方案

本文提出了LTRAN2的一种改进方案BTRAN2用来分析翼型的各种频率的非定常运动。用E-O调转换格式和ADI方法求解了完全的二维非定常跨音速小扰动位流方程,并用单调转换的AF2格式计算了定常跨音速小扰动方程,以此作为非定常计算的初场。本文给出了带后缘正弦振荡襟翼的NACA 64 A006翼型的绕流和做正弦俯仰振荡的NACA 64 A010翼型的绕流计算结果,它们与Euler方程解或实验数据很吻合。     

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

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

多体分离非定常气动特性数值模拟

发展建立了用于研究多体分离的非定常流场数值模拟方法,该方法采用网格弹性变形的结构网格动网格技术,采用有限体积法求解ALE(Arbitrary Lagrangian-Eulerian)描述形式的三维非定常Euler方程,与刚体运动方程耦合,数值模拟了子弹穿越母弹头部激波的子母弹分离过程,并与采用攻角修正的准定常计算结果进行了对比分析,表明多数情况下两种数值方法得到的结果比较接近。     

绕振动旋成体的跨声速流数值分析

本文采用Ｅｕｌｅｒ方程和ＴＶＤ格式数值求解了小幅振动旋成体的非定常跨声速绕流问题，，假设在绕小幅振旋成体的非定常气流中，物理量的变化为平均定常扰动物理量的高阶小量，然后对Ｅｕｌｅｒ方程进行摄动展开，导出前两阶（零阶及一阶）摄动量的控制方程和相庆的边界条件，并进行数值求解。文中仅给出了振动旋成体的非定常压力分布。     

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

飞机亚跨超绕流的数值模拟研究

本文以Euler方程为数学模型,采用一种高精度的TVD(Total Variation Dimishing)离散格式及一种含近似因式分解的推进迭代方法,求解亚跨超绕流’流场。通过若干算例的试算,证明方法是可行的,可以用来模拟飞机的复杂流场。     

压气机叶片排中三维流场计算的多网格法

以三维有限差分形式的MacCormack格式为基础, 对相对圆柱坐标系下弱守恒型控制方程, 导出了残差在粗、细网格上的传播与回插公式, 形成了一种多网格算法。计算表明, 这种算法大幅度地提高了计算收敛速度, 而且方法本身简便易行。      

跨声速翼型绕流的Euler／边界层方程干扰数值解

本文利用Euler方程和可压缩湍流边界层积分方程研究绕跨声速翼型的有粘与无粘强干扰流动。应用有限差分法在贴体的网格上求解时间相关的Euler方程,以剪功积分方法求解翼面贴附和分离湍流边界层流动,并引入一个松弛方程描述剪应力对上游湍流历程的延迟响应。有粘/无粘干扰采用表面源模型。计算结果表明,对翼面存在强干扰流动情况,获得了与实验值基本吻合的结果。     

基于激波拟合与柔性壁的反设计数值方法研究

将激波和待求激波压缩面作为运动间断边界拟合,研究了超声速来流下的型面反设计数值计算方法.算法秉承拟合法思想,采用二阶格式时间推进欧拉方程组,在运动边界和动网格条件下,按预估的波系结构划分网格块,反求可产生指定压力分布的三维型面,并获得精确的激波型面.通过包括凸包进气道在内的三个反设计算例验证了本文算法的合理性及其对网格质量的适应性,证明了在型面反设计中,采用激波拟合法可大大提高型面精度.本算法可用于超声速进气道压缩面和唇口完全三维反设计.     

求解可压缩N－S方程的无粘－粘性分数步法

用按物理过程作时间分裂的无粘－粘性分数步法求解可压缩Ｎ－Ｓ方程，无粘步用高效欧拉算法解欧拉方程，粘性步用全隐格式解抛物方程；并采取了有效措施消除由分步法产生的时间和空间不相容性，文中给出了轴对称喷管内外流场的数值模拟结果。