喷管分离流动与侧向载荷定常数值模拟

利用商业CFD软件,进行了三维、有粘、定常计算,模拟了某液体火箭大面积比喷管地面条件下的分离流动.计算了多个不同入口总压下的工况.其结果预示了该喷管的分离流动在不同入口总压下的激波模态变化,并获得了非对称分离及其造成的侧向载荷的分布情况.计算为进一步相关研究打下了基础.      

面向三维激波问题的装配方法

给出了一种基于非结构动网格技术的三维激波装配方法。在该方法中,三维激波面由被标记为激波属性的网格点连接构成,标记为激波属性的网格点称为激波点。激波点具有两组参数分别代表激波的上下游,利用激波点上下游参数求解R-H关系式获得激波点运动速度。非结构动网格技术的使用允许激波大幅度运动,降低了对初始激波位置的要求。通过引入网格属性定义避免了对计算网格进行分区,增加了装配激波的灵活性。通过球柱体绕流问题验证了该三维装配方法的合理性,针对三维激波装配中比较困难的交点装配问题,通过对三维激波反射以及三维激波相交等算例进行研究找到了可用的三维激波交点运动速度的确定方法,保证了激波运动过程中交点运动与流场求解之间的相容性,获得了相应的装配结果。     

AUSM+格式的改进

为了提高AUSM 格式对流场计算中激波和剪切流的分辨率,减小数值振荡和激波后数值过冲问题,本文在原有AUSM 格式的基础上经过数值研究分析,对AUSM 格式进行了声速、特征马赫数的处理方式,对网格界面马赫数等进行了改进,并用改进后的AUSM 格式求解了含有脱体激波、斜激波、剪切流、膨胀扇、及湍流边界层等流动特征的复杂流动.经典型算例检验证实,改进后的AUSM 格式对激波和剪切流的分辨率得到了加强,求解的准确性和收敛性有一定的提高.改进后的AUSM 格式继承了原AUSM 格式优点,计算处理问题能力增强,可以有效地应用于亚声、跨声和高超声流动的工程实际问题中.     

非结构网格下可压缩黏性流体的隐式算法

采用牛顿线化和Roe近似Jacobians矩阵隐式算法,实现基于非结构网格的三维定常黏性流动问题的求解.分别对平板、跨声涡轮叶片Mark-Ⅱ和CW-22这三种具有代表性的流动进行数值模拟,并详细分析了流动细节,模拟结果与理论以及实验值吻合良好,准确预测了流场内激波、二次流等重要的流动现象.验证了基于非结构网格的可压缩流体隐式计算方法对于黏性流动问题的计算是正确的.      

吸气式飞行器进气道唇口三维激波/激波干扰

对吸气式飞行器进气道唇口处三维曲面激波/弓形激波干扰流场进行数值模拟,利用典型三维气动干扰试验对采用的数值计算方法进行验证。利用拼接网格技术及逆距离加权插值方法获得入口处流场的非守恒变量,作为激波干扰研究的入口边界条件。数值模拟表明,唇口处激波干扰流动的三维效应十分显著,曲面激波与弓形激波产生斜交,尽管唇口前缘半径很小,但Edney提出的6类激波干扰类型可能沿唇口展向方向同时存在；第Ⅲ和Ⅳ类激波/激波干扰的诱导使得唇口热流分布异常严酷；激波相交处形成斜向“伤疤”状局部高热流条带,峰值热流达到参考热流的4~6倍,可能引起唇口结构的局部烧蚀或破坏。      

进气道内激波/边界层干扰及控制研究进展

进气道作为高速航空航天飞行器的重要气动部件,对飞行器的气动力特性、结构重量、隐身性能等有显著影响。激波/边界层干扰现象是高速进气道内普遍存在的一类流动现象,对进气道的性能有突出的影响。发生于进气道内的激波/边界层干扰现象主要可分为正激波/边界层干扰、斜激波/边界层干扰以及三维激波/边界层干扰几类,由于受到侧壁壁面和进气道内背景波系的影响,这些干扰现象偏离了传统基于简化模型的研究结果,具有显著的耦合干扰特征,干扰区间内三维特征明显。概述了发生于进气道内的激波/边界层干扰特性及相关研究进展,并对目前进气道内激波/边界层干扰现象的控制方法进行了总结。     

模拟存在运动激波流场的无网格自适应算法

将无网格自适应算法发展用于求解存在运动激波的非定常流动问题。为了降低初始布点要求,同时提高对运动激波的分辨率,提出了一种结合点云加密和点云粗化的动态点云技术。采用该技术,在非定常计算的每一个时间步,首先通过基于压强梯度的流场探针识别出激波所在位置,然后对激波附近点云进行加密,从而提高对激波的分辨率,同时对远离激波的点云进行粗化,尽可能减少冗余节点,提高计算效率。采用所提出的无网格自适应算法求解了N-S方程,首先对一维运动激波问题进行了数值模拟,通过将数值模拟结果与精确解比较,验证了所提算法求解存在运动激波的非定常流动的可行性,在此基础上,将方法推广用于二维非定常流动问题的求解,分别对瞬时激波碰撞圆柱以及双弧翼型激波振荡流动进行了数值模拟,计算结果表明,方法可以显著提高对激波的分辨率,同时相对于直接采用较密布点的非自适应算法,自适应算法在计算效率上具有优势。     

高超声速喷管三维热化学非平衡流场的数值模拟

利用数值求解三维热化学非平衡Navier-Stokes方程的方法和多块网格技术,数值模拟了典型高焓激波风洞锥型喷管的三维非平衡流场,分析了热化学非平衡效应对喷管流动的影响,给出了平动温度、组分N2和O2振动温度、马赫数和组分质量分数在轴对称线上的分布,也给出了喷管出口附近的速度和温度沿径向分布等结果,获得了喷管流场的详细信息,并与已经存在的采用轴对称Navier-Stokes方程计算的结果,进行了比较,二者吻合得很好.计算结果表明:目前的代码是可以模拟多块结构网格下的热化学非平衡流动的.     

热化学非平衡来流条件下热化学模型影响研究

为研究热化学非平衡来流条件下热化学模型等计算设定对斜激波压缩流动计算结果的影响,针对尖劈构型和相应的前缘钝化构型的高焓激波风洞实验,采用多种计算设定开展详细的数值模拟研究。计算结果表明,计算采用不同热化学模型,以及来流设定为振动冻结/平衡/非平衡状态,会导致斜激波激波角等参数存在一定差别,其中激波角差别可达约2%。当来流速度一定时,过斜激波后分子内能增量在平动转动能和振动能上的分配方式的差别决定了激波角的差别。前缘钝化情形下,采用不同计算设定所得激波角之间的关系和尖前缘构型的规律一致;但是,采用不同计算设定所得斜激波到壁面距离之间的关系和尖前缘构型的规律有差别,这源于钝化前缘的激波脱体距离的影响。对于自由来流下的斜激波压缩流动问题,若考虑了分子振动能激发但未考虑热力学非平衡(例如热完全气体模型、考虑空气反应的单温度模型等),就斜激波激波角等参数而言,计算误差比量热完全气体模型计算误差更大。     

斜激波串非对称特性研究

针对马赫数2.0～4.0均匀来流条件下斜激波串流场的非对称特征开展了数值仿真研究和建模分析.通过改变出口背压和入口马赫数条件分析了不同状态的激波串结构形态,利用前缘激波起始点位置差定义了激波串非对称度系数,分析得到了激波串的非对称度系数变化规律.结果表明,斜激波串出现非对称现象时,前缘激波起始点位置和前缘激波角度均不相...     

一种新型多维通量函数在求解平面叶栅跨音速流动的应用

采用一种改进的Ｒｏｅ格式对平面叶栅的无粘流动进行了数值模拟,并与Ｒｏｅ格式的计算结果进行了比较。采用四阶Ｒｕｎｇｅ－Ｋｕｔａ时间差分格式,空间三阶有限体积法。网格生成采用贴体网格近壁面修正的方法。计算结果表明,该格式在减少格式对网格的依赖性,抑制数值粘性方面具有明显的作用。     

一种高阶精度大涡模拟技术用于激波/火焰相互作用

基于Smagorinsky亚格子模型和增厚火焰技术,开发了一种高阶精度的反应流大涡模拟求解器,运用于数值研究边界层对激波与火焰相互作用的影响。该求解器的核心在于加入了基于超黏性的激波捕捉技术和时空三阶精度的两步Taylor-Galerkin紧致（TTGC）有限元格式,并通过对一维Shu-Osher问题和二维激波/气泡相互作用问题的计算,验证了求解器对激波、接触间断和湍流脉动等流动细节的捕捉精度,计算结果与实验数据吻合良好。通过对激波管内激波、火焰与边界层相互作用问题的数值模拟,发现由于激波与边界层的相互作用会产生不稳定的激波分叉现象,激波三分叉点传播速度的发展经历了水平匀速运动、小斜率线性增长和大斜率迅速增长3个阶段,由此揭示了激波分叉促进火焰加速的机理。当火焰面传入激波分叉区后,流场不均匀的回流区起到了稳定火焰的作用,一方面分叉结构内火焰面能够为激波的运动持续地供应热量,另一方面局部超声速区域为火焰的快速传播提供动力,使其能紧跟激波。通过对比相同条件下甲烷和乙烯燃烧的数值结果,发现两者爆震点触发的位置都出现在马赫杆后面,热量释放率的变化趋势也大致相同,但乙烯出现爆震的时间比甲烷早。     

跨音速压气机级三维紊流流场数值模拟

在对跨音速压气机级的动／静双叶排中的三维定常紊流流场数值的模拟中，利用三阶高分辨率ＮＮＤ格式和ＬＵ－ＳＧＳ隐式推进迭代法，既保证了流场中激波的模拟质量和粘性流动特征的正确预估，也实现了求解过程的高效率。对某单级跨音速压气机的计算及其与实验数据的对比证实了本文方法的有效性。将计算推广至多级轴流叶轮机械的情况将会是十分简便的     

复杂旋翼流场的耦合欧拉-拉格朗日数值方法

针对影响旋翼流场求解精度的关键因素“桨叶复杂近体流动”和“尾迹涡畸变”,结合计算流体力学(CFD)方法和黏性涡方法,发展了一套适合于复杂旋翼涡流场分析的耦合欧拉-拉格朗日数值方法:为捕捉桨尖三维效应、激波等细节流场特征,在桨叶近体区域采用CFD方法对其进行求解;针对高雷诺数旋翼流场中桨尖涡的紧凑结构特点,引入黏性涡方法建立了高分辨率的尾迹求解模型;两计算域间的信息交换采用了集中涡源法和边界修正法.应用所建立的计算方法,以旋翼CFD标准验证试验(Caradonna-Tung旋翼)为算例,对尾迹影响明显的悬停状态进行了数值模拟,通过对比耦合边界处流场特征及桨叶表面压力系数分布,验证了方法的有效性.此外,还从旋翼尾迹捕捉精度、涡量耗散特征及计算时间等方面对不同计算方法进行了对比分析,结果表明耦合方法可充分发挥CFD和黏性涡方法各自的优点,在旋翼流场数值模拟方面具有独特的优势.     

NASA共同研究模型翼身组合体数值模拟研究

采用结构网格求解器WiseManPlus软件,对DPW5提供的NASA共同研究模型CRM翼身组合体构型进行数值模拟研究.采用雷诺平均N-S方程,选择SA一方程湍流模型和SST两方程湍流模型,开展了网格收敛性分析及抖振分析计算研究,评估了该软件对此类民机构型的阻力预测能力.按照DPW5的要求,采用会议统一提供的基准结构网格进行计算分析,研究结果表明,计算的巡航设计点的阻力及大攻角下的气动力系数与试验值吻合较好,达到会议统计分析平均水平.计算结果表明,所采用的结构网格求解器WiseManPlus软件计算精度较高,适用于大型民机高速巡航及抖振现象的数值模拟研究.     

Reflection of a moving shock wave over an oblique shock wave

The reflection of a moving shock wave over a wedge immersed in a still gas and the reflection of a wedge induced steady shock wave over symmetrical and asymmetrical reflecting surfaces have received intensive considerations since more than 70 years ago. Here we consider a different shock reflection problem—reflection of a moving shock wave over an initially steady oblique shock wave induced by a wedge immersed in supersonic flow. For the flow condition we considered, five moving triple points, with each connecting an incident shock wave, a reflected shock wave and a Mach stem, are identified. By using the reference frame co-moving with each triple point, the type of each shock wave of this triple point is clarified. The present study is significant in that it treats a new shock reflection problem leading to a new shock reflection configuration and showing potential applications in supersonic flow with unsteady shock interaction.     

高升力标模确认计算研究

针对高升力标模NASA Trapwing全展襟翼构型展开确认计算研究。采用N-S方法,预测高升力模型的流动状态,分析高升力的流动现象及机理,评估自主研发的WiseManPlus软件及WoF90软件预测高升力流动的能力。通过对网格生成、计算策略及数据分析等各个环节进行分析,制定出一套适用于高升力模型工程应用的CFD计算方法,得到与实验值吻合较好的计算结果。结果表明,目前所采用的结构网格求解器WiseManPlus软件及混合网格求解器WoF90软件计算精度相当,适用于模拟高升力流动。     

低速增升装置失速特性预测研究

民机低速增升装置失速特性是具有挑战性的空气动力学难题之一,决定着飞机起、降阶段机翼处于流动分离或尾迹干扰状态时的操纵稳定性,是气动设计必须关注的重要方面。利用基于N-S方程的自研多块结构网格计算软件和混合网格计算软件,对NASA高升力梯形机翼标准模型开展了失速特性计算研究,系统评估了计算软件、计算网格、湍流模型等主要因素对失速特性的影响,获得了有效预测高升力构型失速特性的计算方案,并将其成功应用于工程实际,为增升装置失速特性计算评估提供了有意义的参考。     

跨声速压气机机匣处理的数值研究

介绍了所开发的压气机内流场三维N-S方程计算程序.在利用典型算例验证了计算程序的描述能力和精度的基础之上, 对带有顶隙和机匣处理的跨声速压气机内的流动进行了三维数值模拟.研究了环向槽机匣处理方式, 比较了不同开槽方式下, 机匣处理对顶隙泄漏流和压气机性能的影响.      

Review of shock wave detection method in CFD post-processing

In the present computational fluid dynamics (CFD) community, post-processing is regarded as a procedure to view parameter distribution, detect characteristic structure and reveal physical mechanism of fluid flow based on computational or experimental results. Field plots by contours, iso-surfaces, streamlines, vectors and others are traditional post-processing techniques. While the shock wave, as one important and critical flow structure in many aerodynamic problems, can hardly be detected or distinguished in a direct way using these traditional methods, due to possible confusions with other similar discontinuous flow structures like slip line, contact discontinuity, etc. Therefore, method for automatic detection of shock wave in post-processing is of great importance for both academic research and engineering applications. In this paper, the current status of methodologies developed for shock wave detection and implementations in post-processing platform are reviewed, as well as discussions on advantages and limitations of the existing methods and proposals for further studies of shock wave detection method. We also develop an advanced post-processing software, with improved shock detection.