基于各向异性非结构网格的超声速流动自适应计算

基于各向异性非结构网格,实现了自适应求解技术,将其应用于超声速楔形体绕流及超声速横向喷流问题的计算研究,对各向异性网格自适应计算过程的收敛性和求解精度进行了对比分析,展示了各向异性网格自适应算法在降低问题求解规模。提高各向异性流场分辨率等方面的优势及存在的问题。超声速楔形体无粘绕流场的自适应计算,清晰地捕捉到了激波附近的流场信息,激波前后的马赫数、压力、密度。温度等物理量与理论分析解吻合很好。超声速横向喷流流场存在激波、分离涡、边界层等流场结构的相互干扰。计算研究表明,单纯基于Hes-sian度量张量的各向异性网格生成及自适应算法不能有效模拟边界层内的流动情况,是将来需要进一步开展的研究挑战。     

空气节流对乙烯燃料超燃冲压发动机流场结构影响研究

空气节流可以有效地实现超燃冲压发动机燃烧室燃料的稳定燃烧,采用非定常数值模拟和地面试验研究了空气节流作用下乙烯燃料超燃冲压发动机的流场特性。结果表明:在燃烧室入口马赫数2、静温530K,静压0.1MPa,845mm处节流条件下,冷流流场达到稳定所需时间较短,约为2ms;在空气节流开始后的6ms注入燃料实现发动机的起动点火是合适的;28ms后燃烧流场稳定,流场呈现短周期的小幅振荡,周期约为0.86ms,振荡对流场的影响很小;壁面压力数据的数值模拟结果与试验结果匹配良好,数值模拟更加全面地给出了流场结构信息。     

三维贴体坐标系下燃烧室热态流场的大涡模拟

大涡模拟三维贴体坐标系下环形燃烧室火焰简热态紊流瞬态流场。利用椭圆方程方法生成三维贴体网格，计算中采用k方程亚网格尺度模型估算亚网格紊流粘性；亚网格EBU燃烧模型估算化学反应速率；热通量辐射模型估算辐射通量。并在非交错网格系下采用sIMPLE算法和混合差分格式求解离散方程，利用壁面函数处理固壁边界条件。计算结果与实验结果的比较表明，采用大涡模拟方法能更真实反映环形燃烧室火焰简内紊流化学反应流气流结构和燃烧过程。     

螺旋桨飞机滑流非定常数值模拟研究

采用非定常数值模拟方法,基于滑移网格技术,分析研究螺旋桨滑流对全机气动特性的干扰影响。用多块结构网格数值模拟软件,从搭接边界建立、搭接网格生成、流场插值传递3个方面对螺旋桨滑移搭接网格的数值模拟方法进行全面详细的介绍,并且阐述基于该技术的非定常数值模拟方法。在此基础上,开展螺旋桨飞机模型验证计算,分析研究螺旋桨滑流对全机流场和气动性能的影响,通过与风洞实验数据的对比分析,验证滑移搭接网格技术的可行性与准确性。研究结果表明,基于滑移搭接网格技术的螺旋桨滑流数值模拟,可以真实反映螺旋桨滑流流场特性,为螺旋桨类飞机的设计评估提供一种有效的数值模拟手段,在工程方面具有重要的实用价值。     

壁温对氢燃料超燃燃烧室流动特性的影响研究

为了研究超燃燃烧室壁面换热与超声速燃烧之间的关系,运用FLUENT软件,采用RNG k-ε湍流模型、有限速率/涡耗散燃烧模型、密度基AUSM+隐式算法对中国空气动力研究与发展中心的双模态燃烧室模型开展三维冷态和热态流场计算,模拟条件来流马赫数为2.05,总温T_t为1870.9K,分别模拟了壁面温度为500K,700K,900K,1300K以及绝热条件下的超燃燃烧室的燃烧场。燃烧室壁面压力数值模拟结果与实验结果吻合较好,壁面温度为500K,700K,900K,1300K和绝热时,平均误差分别为8.89%,5.78%,14.41%,13.97%,16.53%。通过对比分析发现:随着壁面温度的降低,壁面压力趋势大致不变,但壁面压力值降低,同时壁面压力的压升起始点大幅后移;燃烧所产生的激波串逐渐向燃烧室下游移动,激波串结构发生改变,但激波串前端均为X形激波;燃烧室内马赫数有所升高;燃烧场高温区域面积减小;燃烧室燃烧模态由亚燃模态逐渐向超燃模态转换。     

双旋流燃烧室主燃区流动特性PIV测量和分析

针对采用斜切径向双级旋流器的环形燃烧室单头部矩形模型,利用非接触式测量方法粒子成像测速仪(PIV)测量了主燃区的300K冷态速度场。采用Realizableκ-?湍流模型和稳态层流火焰面燃烧模型对燃烧室的冷流和燃烧流场进行数值模拟,得到燃烧室流场的速度分布和流场结构,并与试验测量数据进行对比验证。结果表明:瞬态流场结构变化剧烈,旋流和主燃射流的边界形成大量小尺度漩涡结构,回流区具有强烈的搅拌作用,回流区下游滞止点位置是随机变化的;反向旋流器比同向旋流器产生的回流区尺寸更小,燃烧状态的回流区尺寸比冷流的小,但主要受火焰筒壁面和主燃射流的约束;外旋流在距离头部5mm距离内控制内旋流,保持旋向相同;燃烧增大主燃射流穿透深度,改变流场的对称性。     

实际进口条件下超声速燃烧室全流场研究

针对自行设计的超燃冲压发动机进气道/隔离段+燃烧室一体化模型,进行了燃烧室实际超声速来流条件下的冷态全流场试验研究与数值模拟,获得了模型在不同来流总压和燃烧室出口反压下的流场结构及沿程静压力分布。试验真实模拟了冷流条件下燃烧室实际的进口环境,结果表明:当燃烧室反压系数较低时,激波串在燃烧室扩张段附近形成;当燃烧室反压系数较高时,激波串在燃烧室进口处形成,并且由于凹槽的影响,激波串向燃烧室上壁面偏斜。     

三维直角叉树切割网格Euler方程自适应算法

基于叉树数据结构，实现了一种用于三维直角叉树切割网格的自适应算法，包括对几何外,形以及对流场计算的自适应网格加密技术。在初始网格的生成过程中，根据相邻网格的物面法向向量间的差值，进行针对外形的自适应网格加密；在流场计算中，根据相邻网格间选定物理量梯度的变化，进行针对流场的自适应网格加密。详细地描述了三维直角叉树切割网格的生成过程，以及以任意网格的切割细分算法。在自适应过程中，分别采用了八叉树和全叉树的数据结构，八叉树是基本的数据结构，而全叉树的采用，使网格具有了各向异性的特征，从而大大的减少了自适应网格的数量。采用中心有限体积法，求解Euler方程，并运用上术方法，完成了对外形和流场的逢适应网格加密算法，获得了较好的数值计算结果，证明了自适应算法的正确性，体现了直角叉树切割网格自适应技术的有效性和实用性。     

非结构混合网格鲁棒自适应技术

数值格式、湍流模型和计算网格是影响CFD数值模拟精度的3个主要因素。结合流场信息的网格自适应技术具备动态优化计算网格的能力，被NASA列为未来CFD发展的一项关键技术。本文针对非结构混合网格，发展了网格单元分布优化、表面网格几何投影和空间网格协调匹配3项关键技术，建立了高鲁棒性几何保真的网格自适应系统。首先，为了提高自适应方法的鲁棒性和通用性，发展了基于标准面网格的多面体网格单元分布优化方法。其次，发展了仅依赖表面网格信息的局部曲面重构技术，采用参数点映射方法实现了新增表面网格点的几何投影，消除了自适应系统对几何CAD系统的依赖。再次，采用改进的距离函数方法实现了空间网格与投影后表面网格的快速匹配。最后，结合基于流场特征的自适应探测器，采用二阶格式的有限体积方法，开展了30P30N三段翼绕流和三角翼大迎角绕流的网格自适应数值模拟。结果表明，通过网格自适应对网格单元的分布进行优化后，流场求解的收敛性和模拟精度都得到了显著提高。     

超燃冲压发动机隔离段流动特性研究

隔离段是双模态超燃冲压发动机实现双模态和模态转换的一个重要部件,同时它把进气道和燃烧室隔离开,以防止燃烧室工作对进气道干扰,引起进气道不启动。由于隔离段内流场的复杂性及其广泛的工程应用前景,隔离段内的流动特性引起了人们广泛的关注。本文利用计算流体力学软件Fluent,采用k-ε湍流模型,近壁面采用非平衡壁面函数方法处理,控制方程的离散采用二阶迎风格式,对隔离段的流场进行了数值模拟,并分析了出口反压、来流马赫数等参数对隔离段内流场的影响。     

A mesh optimization method using machine learning technique and variational mesh adaptation

Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation. In light of the introduced large amount of computational costs for many adaptive mesh methods, moving mesh methods keep the number of nodes and topology of a mesh unchanged and do not increase CFD computational expense. As the state-of-the-art moving mesh method, the variational mesh adaptation approach has been introduced to CFD calculation. However, quickly estimating the flow fi...     

气体动理论BGK格式的网格自适应方法

为了提高气体动理论BGK（Bhatnagr-Gross-Krook）格式在超声速流动问题计算时激波捕捉的准确性与计算效率,提出了一种适用于气体动理论BGK格式的网格自适应加密方法。该方法采用基于四边形的链表技术来描述网格的拓扑结构,在物理量重构过程中,使用了在四边形网格中表现优异的van Leer限制器,以保证粗细网格过渡处物理量重构的精度。用跨声速翼型绕流（马赫数Ma=0.85）、超声速前台阶流（Ma=3）和高超声速圆柱绕流（Ma=8.03）等多个典型算例验证了BGK自适应网格方法。计算结果表明,自适应网格BGK方法在保证数值精度的前提下,可大幅度提高计算效率。这为该方法用于高效地解决复杂问题提供了一种选择。     

基于自适应网格的超声速液体射流破碎过程研究

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

带主动冷却的超声速燃烧室传热分析

介绍了流体、固体传热耦合的一维分析方法,对带主动冷却系统的马赫数2.5超声速燃烧室进行了传热分析。该分析以实验测量的燃烧室壁面静压以及超临界煤油换热特性数据为基础,考虑了燃气的高温离解效应,燃烧特性以及碳氢燃料的高温热物理特性,对不同燃烧状态、冷却条件下的主动冷却过程进行了分析。结果表明有燃烧时壁面热流可高达1MW/m2以上,是无燃烧时的2～3倍。当煤油流量较小时(当量比为0.45),冷却后的壁面温度仍偏高,而且冷却壁内温度分布不均匀。随着煤油流量的增加,冷却效果明显提高,冷却壁内温度分布趋于均匀;并且煤油的出口温度也显著减小。     

基于自适应网格的钝头体激波诱导燃烧现象数值模拟研究

钝头体激波诱导燃烧是爆震研究的一个基本问题。针对化学恰当量比的H2/Air预混气在Ma=4.79和Ma=6.46时的激波诱导燃烧现象开展数值模拟研究，采用基于有限体积法的块结构自适应网格加密程序AMROC对带化学反应源项的轴对称Euler方程解耦求解，考察了数值模拟中不同形式的MUSCL重构格式、限制器类型以及化学反应机理等重要因素对模拟结果的影响。结果表明，程序能够根据设定的加密判据较好地实现网格自适应加密，减小总网格量，实现高效数值模拟。通过与实验数据的对比，表明非定常激波诱导燃烧算例的准确程度不仅取决于化学反应机理，也取决于限制器类型，而采用两种不同形式的MUSCL重构格式获得的振荡频率则几乎一致，与试验结果的误差分别为1.17%和0.97%。模拟对比经典的Jachimowski机理和近年来新发展的几种包含压力相关反应步的氢/氧反应机理，模拟结果表明：对于Ma=4.79时的非定常激波诱导燃烧模拟，经典的Jachimowski机理仍然是能够给出与实验结果最接近的反应机理；而对于Ma=6.46时的定常激波诱导燃烧模拟，几种反应机理均能给出与实验吻合较好的结果。     

NLR-7301高升力翼型绕流评估计算

探讨了在二维高升力多段翼型绕流计算中附面层网格采用各向异性三角形网格的技术。通过对经典的NLR-7301两段翼型绕流的计算,从数值计算角度评估了自研混合网格雷诺平均N-S方程求解软件对该类型附面层网格的解算能力。通过与确认试验数据对比,在附面层采用各向异性三角形网格进行二维高升力多段翼型绕流计算,其结果具有较高的可信度。在此基础上,采用同样策略,评估了多段翼型绕流计算中翼型后缘的几何处理、湍流模型选取等关键CFD要素对计算结果的影响。     

支板喷射超声速湍流燃烧的大涡模拟

为了发展可行的超声速湍流燃烧大涡模拟方法,将设定型PDF(Probability Density Function)模型与LES(Large Eddy Simulation)相结合以封闭亚格子湍流-燃烧相互作用,并将模型用于支板喷射超声速湍流燃烧流场的数值模拟。分别对冷流及燃烧流场进行了模拟,计算结果与实验测量符合较好,表明了所采用方法及模型的可行性。冷流条件下,大尺度湍流涡通过卷吸、拉伸运动主导支板尾迹区的近场混合,并通过破碎过程影响远场混合。燃烧条件下,回流区尺度扩大,剪切层中形成的高温燃烧产物通过大涡卷吸以及回流区末端对流作用进入回流区并与其中的燃料喷流相互作用,使部分燃料预热升温并进入回流区两侧剪切层与主流空气混合、燃烧,从而实现火焰稳定。在薄反应剪切层及大尺度反应涡的边界区域,LES网格不足以直接求解湍流与燃烧的相互作用,PDF模型给出了较强的亚格子脉动。     

三维非结构自适应多重网格技术

发展了一种基于几何外形的非结构自适应多重网格技术。根据流场参数的变化梯度确定加密边,由加密准则对粗网格进行自适应剖分得到分布合理的较密网格。通过预先生成的初始极密表面网格来将边界网格的加密点投影到边界上,使得边界保持初始外形。为了提高流场求解速度,将自适应加密后的一系列网格作为多重网格,采用了多重网格算法进行计算。M6机翼的流场计算表明该方法可以准确捕捉激波位置,并且表面的压力系数分布与实验值相吻合。     

超声速燃烧凹腔质量交换特性的混合RANS/LES模拟

对超声速冷流及燃烧流条件下超燃冲压发动机火焰稳定凹腔的质量交换特性进行了研究.采用混合RANS/LES(Reynolds-averaged Navier-Stokes/large-eddy simulation)方法对非定常流场进行数值模拟,系统研究了凹腔结构参数和横向喷流燃烧放热对凹腔质量交换特性的影响.首次采用"原子追踪法"对化学反应流条件下的凹腔质量交换过程进行了刻画.结果表明:在冷流条件和燃烧流条件下,凹腔驻留时间均随长深比的增大而增加;凹腔后缘倾角对驻留时间的影响在冷流条件下比较显著,而在燃烧流条件下明显减弱;相同结构的凹腔,在燃烧流条件下的驻留时间相对于冷流条件明显减小.      

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。