跨音速小扰动方程的多重解

本文利用Engquist-Osher格式,离散跨音速小扰动方程(TSD)和边界条件,选择不同初场计算NACA0012翼型二维定常位势流场,得到了与试验值较符合的结果,同时用Murman-Cole非守恒和守恒两格式计算并进行了比较。着重用数值实验进一步探讨多重解问题,发现并总结出了几点很有参考价值的规律和现象。     

近壁湍流相干结构的数值模拟

提出了近壁湍流对称单个相干结构的理论模型,采用紧致有限差分和Fourier谱展开相结合的数值方法研究它的演化问题,得到雷诺应力分布及高剪切层的形成过程与实验结果一致。     

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

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

大扰动跨音速流积分方程的导出及亚临界流的数值解法

本文对一方向上大扰动的平面跨音速势流微分方程提出了新的近似简化方案,从数学上完成了相应积分方程的推导;其中关于边界及激波条件作了小扰动近似处理;最后提供了关于亚临界流的算例。     

典型航空分离流动的雷诺应力模型数值模拟

先进飞行器设计对CFD方法的边界层分离模拟能力提出了更高要求,传统雷诺平均纳维-斯托克斯(RANS)涡黏性模型因建模和构造层面的理论缺陷导致其分离流动预测可信度较低。雷诺应力模型由于未对雷诺应力及其分量关系进行建模,对湍流非平衡、旋转以及雷诺应力各向异性等流动现象具有天然的理论优势。为验证与确认雷诺应力模型对典型航空分离流动的预测能力,基于SSG/LRR-g模型,以NACA4412翼型大攻角分离、M6机翼跨声速分离以及F6翼身结合区分离流动为例,探讨了雷诺应力模型对逆压梯度、激波诱导分离、二次流动分离等典型航空分离流动预测的适应性。通过与k-ω剪应力输运(SST)模型模拟结果对比发现,雷诺应力模型对分离泡大小、速度型分布、雷诺应力分布和激波位置等关键特征的模拟精度较涡黏性模型显著提升,基本验证了雷诺应力模型可在翼身接合区角区流动和三维强激波诱导分离等问题中得到正确的流动特征,而SST模型在此类流动中基本失效,显示了雷诺应力模型在典型航空分离流动中较涡黏模型的优势。同时,发现k-ω SST模型所包含的Bradshaw假设在三维激波诱导分离较强时严重影响了模型预测的准确性,是预测结果偏离试...     

超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

壁面湍流模型对湍流分离流动数值模拟的影响

本文用三种近壁湍流模型计算了二个二维的湍流分离流动：单包流动和多包流动。结果表明：二层模型和低Reynolds数模型具有类似的特性。它们基本上给出合理的分离流动结构。壁函数的方法由于对数律的假定基本上不适合于计算湍流的分离流动。     

壁温对压缩拐角流动影响的数值模拟研究

压缩拐角是高超声速飞行器上的典型非连续区域,其分离/再附结构对局部热环境有较大影响。本文采用基于Navier-Stokes方程的自研程序,开展了典型压缩拐角外形的气动热环境数值模拟研究,分析了不同壁温条件下压缩拐角的热环境分布规律。计算结果表明:随着壁温升高,流场整体结构变化不大,但由于近壁面流体物性参数变化较大,导致拐角分离涡分离点向前移动、再附点向后移动,拐角处干扰区扩大;压缩拐角大部分区域的热流随着壁温的升高而减小,且干扰区内热流减小的幅度比无干扰区更大,但热流并不完全遵循随壁温升高而减小的规律。另外,通过对比变壁温计算热流和热壁修正公式修正热流,发现热壁修正公式在干扰区存在精度降低、适用性不足的问题。本文的研究进一步加深了壁温对压缩拐角流动影响的认识,缩小了热壁修正公式的适用范围。     

测量三维湍流场的热线响应方程及热线校正方程

从Ｊｏｒｇｅｎｓｅｎ等效冷却速度公式出发推出了高阶精度的热线响应方程,它适用于密度不变低速条件下求解湍流场中时均速度和雷诺应力张量。应用该热线响应方程时不严格要求热线探头轴线与流场速度矢量完全一致,因此可简化热线测量的支撑和传动机构。导出的显式热线校正方程可保证校正精度。     

钝锥有攻角分离流动的数值模拟及其分析

本文定性分析了开式分离的性状,并对钝锥有攻角超声速绕流的开式分离作了数值模拟。分析指出,开式分离可能存在两种形态,第一种分离线的起点为正常点,第二种分离线的起始为鞍、结点(包括螺旋点)的组合。对于第一种形态,分离线的起点是横向分离的起始点,除分离线外,分离面上的流线不是从分离线的起点发出的。对文中计算的情况,流动属第一种开式分离。计算证实了定性分析的结论。计算和分析均指出,对第一种开式分离,在分离的起始区域,分离流面尚未卷曲,但在下游,则变成卷曲面。文中还研究了围绕物体的流管在分离诱导下的变形情况。     

非定常分离流动的Lagrange方法数值模拟

本文对Ｌａｇｒａｎｇｅ坐标系下的边界层方程进行了数值模拟。研究了静止流场中的线涡诱导的非定常边界层流动。数值结果表明，随着边界层接近分离，确实出现了有限时间奇性。流动在驻定点附近很窄的区域内堆积，并在垂直壁面方向迅速增长，形成一脉冲似的突跃。基于此现象讨论了非定常分离的性质和三种分离模式。采用Ｌａｇｒａｎｇｅ方法能很好地说明有限时间奇性，便干精确地给出分离的定义和判据，而且物理常清楚。Ｌａｇｒａｎｇｅ方法有着明显的解析和计算上的优点。     

运用SIMPLE方法计算大曲率弯道中的湍流流动

在三维任意贴体非正交曲线坐标系下应用同位网格技术,发展了适于求解复杂流动的压力修正算法,采用速度的协变物理分量为动量方程中速度的求解变量,湍流模型选用标准的κ-ε模型.计算表明,方法具有良好的收敛性、准确性和稳定性,计算值与试验数值吻合良好,而且能预示通道内复杂的流动现象.     

二维旋转湍流边界层内气固两相流动的数值模拟

本文用修正的ＤＯＲＯＤ－ＦＥＭ方法对不可压流体二维旋转湍流边界层进行了数值模拟，得到了与实验数据较符合的边界层形状因子和壁面摩擦系数的计算结果。采用单向耦合（Ｏｎｅ－ｗａｙ Ｃｏｕｐｌｉｎｇ）方法，计算了固粒在湍流边界层和无粘流中的运动轨迹，给出了壁面被固粒撞击后的冲蚀磨损分布，分析了近壁流场性态与壁面冲蚀磨损的关系。     

裙后向台阶超声速绕流特性数值研究

本文研究了在中部带有裙体后向台阶的轴对称体的超声速绕流特性。采用空间推进法数值模拟超声速无粘流场，对于轴对称后向台阶的分离区采用Ｃｈａｐｍａｎ－Ｋｏｒｓｔ理论模型进行处理，计算结果给出了与实验值吻合的物面压力分布以及超声速流场的流动图画。     

脱体涡流型翼面非定常气动力计算方法

本文研究了细长三角翼在脱体涡流型时变攻角所产生的非定常气动力的机理,并采用了一种时间历程的,考虑翼面前、后缘涡面随局部气流拖出并由时间发展逐步形成的数值方法,称之为涡面生成法,计算结果表明该方法是稳定和有效的。     

大曲率弯道和叶栅内的湍流数值模拟

本文在贴体曲线坐标系下应用非交错网格技术，发展了适于求解复杂流动的压力修正算法，采用速度的协变物理分量为动量方程的求解变量，推导得出一个适于任何曲线坐标系下控制单元体界面流通量的插值模式，可以有效地抑制不合理的波动压力场，计算表明，算法稳定收敛，计算值与试验值吻合良好，同时能预示通道内复杂的流动现象。     

二维分离流中的NS方程数值计算方法

本文采用显式多步龙格－库塔方法求解了二维可压非常完全ＮＳ方程，为了使计算适用于一般的风格，采用有限体积法进行几何离散，研究定常问题计算中的一些加速收敛措施，如：当地时间步长、残值光顺技术等。讨论了模拟分离流问题中的技术问题，计算结果表明本文的计算结果与国外文献是吻合的。     

反向喷流与主流干扰数值模拟

反向喷流是改变钝体迎风面压力和温度分布的有效手段之一，所以了解反向喷流场对再入体防热及减阻设计具有重要意义，本文用有差分法求解轴称Ｎ－Ｓ方程模拟了球柱体部反向喷流与超声速自由来流的干扰流场。研究了喷流出口压力对流场的影响，物面压力分布及喷流Ｍａｃｈ盘、弓形激波的脱体距离与实验值进行了比较，取得了一致的结果，喷流附近分离区的再附位置及再附产生的激波位置也与实验相符，但高压比下再附激波略强于实验。     

钝锥高超声速绕流空间分离形态的数值模拟及拓扑分析

本文通过数值模拟和理论分析,系统研究了钝锥高超声速有迎角绕流时垂直于物面的截面流态沿轴向的演变过程。得出了高超声速情况下,背风分离区内的流态在大多数情况下貌似对称,其实是非对称的这一重要结论。考察了周向网格疏密程度对流动结构分辩率的影响,发现对分离流动的数值模拟,必须在分离区内布置足够密的网格才能给出正确的结果。     

吸气条件对圆柱非定常分离流影响的数值研究

本文根据二维可压非定常Ｎ－Ｓ方程，基于Ｂｅａｍ ａｎｄ Ｗａｒｍｉｇｎ及Ｓｔｅｇｅｒ的方法，在圆柱绕流出现周期分离后，施加吸气条件，研究气对于圆柱非定常绕流的影响，计算结果表明，在恰当地区域施加恰当的吸气条件，可以希望通过减小压差阻力来减小总阻力，吸气还有抑制尾迹中非对称流动的作用，从而使圆柱所受的垂直于来流方向的升力的变化幅度减小，变化频率减小，但是吸气会增加磨擦阻力，过份的吸气会导致摩阻增大太