CFD-FASTRAN气动热计算模型及网格效应分析

介绍了CFD-FASTRAN软件五种湍流模型。通过B-L湍流模型不同网格气动热计算结果的对比分析,说明近壁面第一层网格间距须保证y＋值处于所选湍流模型要求值范围的1/3左右,才能保证计算结果的准确性。进一步通过直二次圆锥形不同湍流模型计算与试验的对比结果,说明CFD-FASTRAN软件的k-ε及B-L两种模型求解高超声速气动热问题能达到较好的准确性,且B-L模型计算速度快、使用方便,更适宜于整机大型数值计算。     

低雷诺数κ-ε封闭模式在三维湍流边界层计算中的应用

本文采用低雷诺数形式的k-ε湍流模型,计算具有大横向流的不可压缩三维湍流边界层。在湍流模型中,应用了涡粘性的各向异性系数。本文所得结果与用混合长模型所得结果以及量测结果的比较表明:对于边界层的不同特征量,湍流模型和涡粘性各向异性对计算结果的影响程度是不同的。总的来说,应用考虑涡粘性各向异性的k-ε模型,较之应用混合长模型,所得计算结果有所改善,但还不是根本性的改善。     

航空发动机燃烧室流动数值计算中湍流模型的比较

分别采用标准k-ε湍流模型、RNGk-ε湍流模型、Realizablek-ε湍流模型以及雷诺应力模型,对某型航空发动机燃烧室流动进行了数值计算.近壁处理采用标准壁面函数法,计算得到速度矢量分布以及质量流量、湍流粘度比和湍流强度等参数.四种湍流模型计算的总体流动差别较小,但流场的细节有较明显的不同.标准k-ε模型、Realizable k-ε模型和雷诺应力模型的湍流粘性计算结果较为接近,而RNGk-ε模型计算的湍流粘性较小.     

绕翼型低雷诺数流动的数值分析研究

分离泡的产生是绕翼型低雷诺数流动的一个重要特征,分离泡通常都是非定常的,会对整个流场产生极大的影响,由于分离附面层的不稳定性,很快诱发转捩,并产生湍流再附。文中通过求解雷诺平均N-S方程,数值模拟了绕Eppler387翼型的低雷诺数非定常流动,并对两方程SSTk-ω湍流模型,代数B-L模型和层流的计算结果作了比较。N-S方程和两方程湍流模型的控制方程非耦合求解,时间推进均采用近似因子方法(AF);空间离散无粘项采用ROE格式,粘性项用中心差分方法计算;计算非定常流场时,采用伪时间子迭代(τt-s)方法保证二阶时间精度,湍流模型计算时都是在固定点转捩。最后分析了转捩对低雷诺数流动的影响、分离泡的存在引起的流动不稳定性和周期性的脱出涡,数值计算给出的时间平均升力系数、阻力系数和压力分布与实验结果比较吻合。     

运用低雷诺数非线性k-ε模型研究湍流分离流动

给出了一个低雷诺数非线性ｋ－ε湍流模型，其中采用了一个非线性雷诺应力表达式并引入了两个考虑湍流近壁效应的修正函数。应用该模型和标准ｋ－ε模型计算了两种不同膨胀比的后台阶分离流。与实测结果比较表明，低雷诺数非线性模型结果较标准ｋ－ε模型结果有明显改进，它对流场的各特征量包括分离—再附长度、时均速度及各雷诺应力分量均给出了较好的预测。     

利用分区ALE算法数值模拟圆柱湍流涡致振动

用分区算法结合任意拉格朗日-欧拉法(ALE)数值模拟了圆柱湍流涡致振动.求解基于非交错网格系统,压力求解采用压力泊松方程提法,湍流模型采用标准k-ε模型和重整化群(RNG)k-ε模型.计算取中等雷诺数Re=5000、10000、15000、25000、50000等,质量系数M=10.阻尼系数ζ=0.00331,自振频率fn=0.18315、0.1628.计算结果表明:湍流涡致振动下圆柱时均阻力系数大于孤立圆柱绕流,而升力系数(振幅)值都小于孤立圆柱绕流.随着雷诺数增大,湍流粘性系数随之增大,但湍动能和湍流耗散率变化趋势不明显.对孤立圆柱绕流,研究结果与前人实验结果基本相符.     

两种差分格式和两种湍流模型在轴对称冲击射流数值计算中的比较

用延迟修正QUICK和乘方两种差分格式，结合标准k-ε模型和一种基于重正化群（RNG）思想的k-ε模型，对半封闭轴对称涡流冲击射流场进行了数值模拟。通过与实验结果的比较，考察了两种差分格式对计算结果的影响，并讨论了两种湍流模型对冲击射流场的数值预测能力。研究表明：当都采用标准k-ε模型时，QUICK格式得到的结果比乘方格式更接近于实验值，尤其对湍流能分布；当都采用QUICK格式时，RNG k-ε模型的预测值比标准k-ε模型更准确。     

五种湍流涡粘模型在二维方柱绕流数值模拟中的对比研究

为研究不同雷诺时均Navier-Stokes(RANS)模型在求解钝体绕流场的差异,采用五种两方程湍流涡粘模型计算了二维方柱的绕流场(雷诺数Re=22 000),得到了不同湍流模型在不同网格离散方案下的方柱气动力与流场特性。结合以往文献数据,通过对比不同湍流模型计算结果的差异,揭示了不同湍流模型的特点。研究结果表明:网格离散方案对方柱绕流数值模拟结果有重要影响;RNGk-ε湍流模型的计算效率最高,SSTk-ω湍流模型的计算效率最低;Standard k-ε湍流模型的计算结果准确程度整体弱于其余湍流模型;RNGk-ε湍流模型、Realizable k-ε湍流模型与Standard k-ω湍流模型的计算结果大致相当,较接近大涡模拟结果;SSTk-ω湍流模型的模拟结果优于其余湍流模型,其尾流区速度场与试验结果吻合较好;两方程湍流模型计算的二维方柱绕流结果与试验结果以及大涡模拟结果相比,存在不可忽略的差异。     

应用改进的低雷诺数湍流模型预测转捩流动

基于Launder-Sharma(LS)低雷诺数k-ε两方程湍流模型发展了一种改进的具有转捩敏感性的低雷诺数湍流模型.针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,模型进行了改进.改进模型只使用当地物理量,不需要求解壁面距离、y+和边界层积分参数.改进模型能够适用于广泛的流动,且容易应用到通用的计算流体动力学(CFD)程序中.对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,改进模型能够预测转捩流动,并能对自由湍流变化给出合理的响应.      

临近空间螺旋桨低雷诺数高效翼型数值分析

针对临近空间翼型边界层变厚更容易发生分离的特点,在修正了湍流模型后,通过分析比较13种低雷诺数翼型的升阻系数得到了比较适合用于临近空间螺旋桨叶素的高效翼型。结果显示,攻角在-6°～8°、8°～14°、14°～22°三种工况下,分别采用SST k-ω、RNG k-ε、Realizable k-ε湍流模型,可以得到与国外实验比较接近的合理结果;工况条件对翼型气动性能好坏的影响很大,翼型S-1223和FX63-137在所研究的工况内都具有较好的气动性能;因此可以选择这两种翼型作为临近空间螺旋桨用高效翼型。     

RNG k-ε模型数值模拟油雾燃烧流场

采用RNG(R enorm alization G roup)k-ε紊流模型对某环形燃烧室火焰筒内三维两相燃烧流场进行数值模拟。运用偏微分方程法和区域法生成贴体网格,EBU-A rrhen ius紊流燃烧模型估算化学反应速率,六通量热辐射模型估算辐射通量。在非交错网格下求解气相采用S IM PLE算法,液相采用颗粒轨道模型和PS IC算法。近壁区处理分别采用两层和三层壁面函数。计算数据与试验值比较表明,采用三层壁面函数的RNG k-ε模型更适用于模拟三维两相燃烧流场。      

环形燃烧室冷态流场数值模拟中的数学方法

针对带双级径向旋流器环形燃烧室火焰筒,研究不同紊流模型(standardk-εmodel,RNG(renormalization group)k-εmodel和代数应力模型(algebraic stress model))及算法(SI MPLE(semi-i mplicitmethod for pressure-linked equations),PISO(pressure i mplicit split-operator))对复杂流场求解准确性及效率的影响.火焰筒三维贴体网格采用型线定点与偏微分方程相结合的方式生成,用壁面函数法处理近壁区域流动.计算结果表明,PISO算法提高求解效率;不同紊流模型预测结果均与试验数据符合较好,其中代数应力模型最为准确,更适用于复杂流场的数值预测.      

封闭转静盘腔中紊流的数值模拟

对一个旋转盘腔的流场进行了数值计算研究。该盘腔由一个旋转盘,一个静止盘及静止的外围屏组成。在盘腔的核心区采用标准k-ε模型,在近壁区分别采用了两种低雷诺数紊流模型,即(1)单方程模型和(2)Launder-Sharma低雷诺数k-ε模型,并且在第二个模型中考虑了固体旋转引起的耗散率降低的修正。      

计算风工程中 k-ε模型的边界条件研究

从 k-ε湍流模型控制方程入手，对两类满足平衡大气边界层理论要求的入口边界条件进行了详细分析，比较了二者的理论差异。然后以我国公路桥梁抗风规范中建议的 A 类风场为来流条件，定义了相应的 k-ε模型常数与壁面条件，使用两类边界条件进行数值模拟。计算结果表明：两类边界条件结合相应的模型常数与壁面条件，均能在数值模拟中构建基本满足规范要求的平衡大气边界层，但各有其适用范围。在此基础上，对模型常数 C 1和湍动能耗散率与数值模拟结果之间的关系进行了研究，采用修改模型常数与边界条件的方法，改善了数值模拟结果。     

湍流状态下超疏水表面流场减阻特性数值仿真研究

通过对湍流状态下具有特定微观尺寸的超疏水表面流场进行数值仿真计算,对超疏水表面流场的减阻特性进行了分析。针对超疏水表面矩形微观形貌特点,计算域采用结构化网格进行划分,采用VOF多相流模型,Realizable湍流模型,对超疏水表面流场进行仿真。结果表明:受微观形貌的影响,超疏水表面在宏观上的壁面滑移、微观凹坑处的低剪应力和近壁面的低湍流度是其具有减阻特性的重要原因;超疏水表面减阻特性受凹坑内空气体积比影响很大,但是在凹坑内全充满液体条件下,依然具有减阻效果。     

过膨胀状态下轴对称收-扩喷管内外流场计算及分析

轴对称收-扩喷管在过膨胀状态下工作时, 喷管内部将产生强烈的激波以及附面层分离, 流态非常复杂.计算实践表明, 应用广泛采用的重整化群(RNG)k-ε湍流模型, 并结合标准壁面函数对其进行数值模拟时, 误差较大.针对标准壁面函数的不足, 采用增强型壁面函数处理固体壁面附近的流动, 对轴对称收-扩喷管的内外流场进行了数值模拟, 并与试验数据进行了分析比较.结果表明:增强型壁面函数更加适合处理此类存在强逆压梯度导致附面层分离的问题.      

狭小空间内气膜孔流量系数的数值模拟

针对多层层板冷却、双层壁冷却等冷却技术,利用数值模拟分析了这些冷却结构形成的狭小冷却通道中,气膜孔附近的流动特性,重点研究了气膜孔流量系数C_{d随吹风比M(0.5～2.0)、气膜孔雷诺数Re(5000～10000)以及冷却通道高度H和气膜孔直径D之比H/D(0.33～1.0)等参数的变化规律.数值计算中湍流模型为Realizable k-ε模型,近壁处采用非平衡壁面函数,利用SIMPLE算法和二阶迎风格式进行离散求解.计算结果表明:气膜孔流量系数随吹风比的增大而增大,在吹风比M小于1时,影响尤为明显.研究中同时发现在相同吹风比的条件下,Cd随着气膜孔雷诺数的增大而减小,但变化的幅度不大;在相同的气膜孔雷诺数下,Cd随着H/D的减小而降低,特别是当H/D小于0.75时,随着流体在进入气膜孔前逐步受到限制,Cd随着H/D的减小而快速降低.}      

扰流柱直径对倾斜出气孔层板内流动和换热影响的数值研究

用数值模拟的方法, 研究了倾斜出气孔层板结构最小单元体内扰流柱直径对于层板内流阻和换热的影响.采用了非结构化网格, 并对气流转角较大区域做了局部加密;选择k-ε湍流模型方程, 配以壁面函数率, 求解稳态无旋转的层板;采用了流体域与固体域耦合求解法计算耦合换热.给出了计算结果并分析了流动和换热情况.研究结果表明适当增加扰流柱直径有利于降低流阻、提高冷却效率.      

二维非对称喷管数值模拟与验证

采用Fluent软件针对二维非对称喷管进行数值模拟,考察了其中几种常用湍流模型及网格密度对模拟精度的影响,并将计算结果与NASA公布的实验数据进行比较,找出了精度较高、求解快捷的模拟方案,结果表明FLUENT软件可以应用于带高超声速外流的非对称超声速喷管内外流场的数值计算。