DLR—F6翼身组合体跨音速流场CFD应用计算研究

针对跨音速运输机经典算例DLR—F6翼身组合体模型,采用CFD方法对其气动特性进行了黏性流动数值模拟,流动模型为雷诺平均N-S(RANS)方程。首先采用"超立方体"概念生成绕DLR—F6翼身组合体的高质量多块结构拼接网格,研究网格拓扑结构对气动特性的影响。在此基础上通过网格细分和粗分考查了网格密度对计算结果的影响,最后进行了湍流模型的影响研究。通过与实验数据对比分析,得出了适宜DLR—F6翼身组合体跨音速黏性流动的计算网格,并总结出了能较好模拟其跨音速流场特性的湍流模型。结果表明:网格拓扑结构的合理设计会对计算结果产生一定的影响。网格密度对机翼表面压力分布没有明显影响,但对阻力系数影响显著。湍流模型对机翼表面压力系数分布的影响主要体现在激波位置上,对翼根处的分离也有一定的影响。SST模型计算的气动力系数比SA模型接近实验值。     

弧形翼-身组合体绕流流场数值模拟及实验研究

为了研究标准技术协作计划(TTCP)弧形翼-身组合体的空气动力学特性,对其在超声速下的三维绕流流场进行了数值模拟.应用雷诺应力湍流模型和平流迎风劈裂法离散格式,对弧形翼-身组合体在Ma∈[2.0,5.5]范围内的超声速绕流流场结构进行了分析.计算分析了弧形翼-身组合体的各气动力系数随马赫数、攻角的变化规律.阐述了弧形翼在零攻角状态下产生自诱导滚转力矩的机理.利用超声速风洞进行了标准TTCP弧形翼-身组合体在Ma=3.0条件下的吹风实验.实验结果与数值模拟结果的计算误差小于10％,验证了所采用数值模拟方法的可信度.     

DLR-F6翼身组合体跨音速流场

针对跨音速运输机经典算例DLR-F6翼身组合体模型,采用CFD方法对其气动特性进行了粘性流动数值模拟,流动模型为雷诺平均NS(RANS)方程。首先采用“超立方体”概念生成绕DLR-F6翼身组合体的高质量多块结构拼接网格,研究网格拓扑结构对气动特性的影响,在此基础上通过网格细分和粗分考查了网格密度对计算结果的影响,最后进行了湍流模型的影响研究。通过与实验数据对比分析,得出了适宜DLR-F6翼身组合体跨音速粘性流动的计算网格,并总结出了能较好模拟其跨音速流场特性的湍流模型。结果表明：网格拓扑结构的合理设计会对计算结果产生一定的影响;网格密度对机翼表面压力分布没有明显影响,但对阻力系数影响显著;湍流模型对机翼表面压力系数分布的影响主要体现在激波位置上,对翼根处的分离也有一定的影响;SST 模型计算的气动力系数比SA模型接近实验值。     

DLR-F6翼身组合体跨声速阻力计算

采用美国航空航天学会阻力测试小组提供的多块对接网格,结合Spalart-Allmaras、Wilcox’sk-ω和Menter’s k-ω SST三种湍流模型,通过求解雷诺平均Navier-Stokes方程、数值模拟DLR-F6翼身组合体的流场来研究阻力计算精度,考查网格和湍流模型对翼身组合体构型气动特性的影响.结果表明:三种湍流模型得到的机翼表面压力系数分布与实验数据吻合良好,气动力随攻角的变化趋势与实验结果一致;Spalart-Allmaras模型得到了网格收敛结果,所得阻力优于其他软件的结果;网格密度对阻力有影响,对机翼表面压力系数分布无明显影响;湍流模型对机翼表面压力系数分布的影响主要体现在激波位置上,对升力影响较小,对阻力(尤其是摩擦阻力)影响显著,对翼根处的流动分离有一定影响;在跨声速流动中,Menter’s k-ω SST模型的结果最接近实验数据.     

二维非定常有涡流动的数值模拟方法

为有效地模拟二维有分离现象的粘性流动,发展了求解涡量流函数方程的数值方法。对涡量输运方程的时间项运用四阶Runge-Kutta法离散,从而将方程分解为4个计算步分别求解。在每一计算步对方程进行Steger-Warming近似因式分解,从而使涡量的计算在空间的两个方向上分别进行。对流项采用Chakravaythy-Oscher总变差减少(TVD)格式离散。流函数的Poisson方程运用Tschebyscheff SLOR方法交替方向迭代求解。运用该方法对突然起动圆柱和高雷诺数时弯曲薄翼的非定常有涡流动进行了数值模拟,并将计算结果与其它方法的计算结果和试验结果进行了对比,结果表明本方法具有精度高、数值稳定性好和计算效率高的优点。     

船舶Z形操纵试验的数值模拟

建立三维黏性流数值水池,采用CFD方法对船舶Z形操纵试验进行数值模拟研究.利用重叠网格技术实现船舶航行中的姿态变化,通过求解RANS方程和船舶运动方程,并结合单相Level Set方法,对船舶在开敞水域黏性流场Z形操纵试验进行数值模拟.以"育鲲"轮为对象,数值模拟10°/10°Z形试验.结果表明,计算结果和实船试验数据吻合度较高.     

自由剪切流动数值模拟的计算网格设计

从流体力学理论的角度出发,采用理论分析与证明的方法,通过干扰剪切流(ISF)的理论分析表明,在自由剪切流动的计算时存在最佳计算网格. 从ISF理论可以证明最佳计算网格是一个网格线与ISF黏性剪切薄层的流向平行,且在ISF黏性薄层法向在薄层内进行局部加密的正交网格. 对于非最佳网格,即使把网格加得很密,也难以捕捉ISF黏性薄层的物理黏性效应,算出的黏性效应是数值黏性,并通过一个不可压缩自由剪切层流流动的数值模拟证实了这个结论.      

具有运动边界不可压缩黏性流动的CBS有限元解法

基于特征分裂(CBS)算法,提出了一种适用于具有运动边界的不可压缩黏性流动问题的有限元解法,给出了不可压缩黏性流动问题的控制方程,并推导了其无量纲形式.将传统CBS算法与动网格技术相结合,提出了一种动边界流动问题的求解流程,并对所采用的基于弹簧近似的动网格方法进行了说明.采用该方法对静水和均匀来流中的简谐振动圆柱绕流问题进行了求解,计算所得圆柱受力、圆柱表面压力和涡量分布以及流场速度分布等结果均与已有数值或实验结果相吻合,初步验证了所提方法的可行性.     

基于Fluent的多圆柱体绕流场数值模拟

 圆柱绕流一直是海洋工程研究热点,对海洋工程有重要意义。在海洋工程应用中,多圆柱间相互干扰现象普遍存在,多圆柱绕流特性一直是研究的热点和难点。为了分析多圆柱绕流场的形态结构,本文应用Fluent软件,采用数值求解RANS方程的黏性流方法,对典型的多圆柱黏性绕流场进行模拟。结果表明,串列双圆柱绕流场与圆柱间距有关;并列双圆柱尾流存在反向对称涡街;斜置双圆柱与3圆柱绕流场存在后柱抑制前柱涡的分离现象;方形排列4柱绕流场尾涡存在非对称性。同时表明Fluent软件对多圆柱体绕流场的数值模拟是有效的,本文采用的方法能定性正确地模拟多圆柱绕流的流动细节。     

导叶可调式液力变矩器导轮特性研究

液力变矩器内部流动复杂,是非定常、不可压缩、三维、黏性流体的流动.可调式液力变矩器在三元件变矩器基础上增加可调导轮.在可调式液力变矩器可调导轮导叶存在多个开度的基础上,建立不同开度下的变矩器全流道模型;采用分离式求解器、SIMPLE算法及大涡模拟湍流模型对变矩器内部流动进行数值模拟.基于变矩器内部流场数值模拟结果,分析其流动特性,研究导轮的工作特性,得到固定导轮与可调导轮的作用,可为可调式液力变矩器性能改善和导叶的设计提供参考.     

组合的LU分解迎风方法的多重网格收敛

为提高计算收敛速度,以新的高收敛率的LU型隐式格式和高精度、高分辨率的MUSCL TVD迎风格式为基础,对可压缩无粘和粘性流动问题计算的多重网格方法进行了探讨。采用二维凸包通道和VKI LS 59跨音速透平叶栅为算例,通过与单网格上的计算迭代性能做对比,证明这种算法在原有基础上更大幅度地提高了计算收敛速度,节省了CPU时间,而且方法本身也简单易行。     

Numerical Simulation of Rotor-stator Interactions in a Transonic Compressor Stage

ＩｎｔｒｏｄｕｃｔｉｏｎＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｒｏｔｏｒｓｔａｔｏｒｉｎｔｅｒａｃｔｉｏｎｓｈａｓｒｅｃｅｉｖｅｄｗｉｄｅｓｐｒｅａｄａｔｔｅｎｔｉｏｎｉｎｔｈｅｌａｓｔｆｅｗｙｅａｒｓｂｅｃａｕｓｅｏｆｔｈｅｉｎｃｒｅａｓｉｎｇｉｎｔｅｒｅｓｔｉｎｕｎｄｅｒｓｔａｎｄｉｎｇｔｈｅｕｎｓｔｅａｄｙｎａｔｕｒｅｏｆｆｌｏｗｓｉｎｔｕｒｂｏｍａｃｈｉｎｅｒｙ.Ｔｈｅｒｅａｌｆｌｏｗｓｉｎｓｉｄｅａｔｕｒｂｉｎｅｏｒａｃｏｍｐｒｅｓｓｏｒａｒｅｕｎｓｔｅａｄｙａｎｄａｒｅｃｈａｒａ…     

二维翼型绕流计算中预处理和多重网格方法的应用

结合Turkel矩阵预处理方法和多重网格方法,发展一种适合低速粘性流动计算的高效数值方法,通过对Navier-Stokes方程的时间导数项实施Turkel矩阵预处理,使得可压缩Navier-Stokes方程在低速情况下的系统刚性得到改善,为进一步加速收敛,提高计算效率,采用多重网格的3层V循环方式,对RAE2822超临...     

叶栅激波流场的数值模拟

给出了ＡＵＳＭ 格式在任意曲线坐标下的形式,将ＡＵＳＭ 格式和三阶精度ＭＵＳＣＬ格式融合,并结合ＬＵ－ＳＧＳ隐式时间推进法,对典型双圆弧叶栅无粘跨音速和超音速流动进行了数值模拟．本文计算结果与文献计算结果相符     

水驱粘塑性非牛顿原油渗流特性的数值模拟

针对水驱粘塑性非牛顿原油在地层下呈强非线性渗流特征,建立了描述其流动过程的数学方程; 采用全隐式离散、牛顿－ 拉弗松方法迭代线化,稳定双共轭梯度方法求解,对特定初边条件下二维油藏油水两相流动问题作了数值模拟．结果表明,水驱非牛顿流体驱油效率明显低于水驱牛顿流体     

投弃式海流剖面仪探头流场数值模拟

 以三维Navier-Stokes方程为控制方程,结合k-ε两方程湍流模型,采用隐式有限体积法对投弃式海洋剖面仪（XCP）探头进行数值模拟,分析XCP探头周围流场的分布情况。根据XCP探头周围流场分布,对探头的运动特性,重点是探头的头部和尾翼流场进行分析。研究表明,在探头迎流前端小于10cm处,流体层流边界层发展为湍流边界层,因此探头头部的弧线型设计对于保证探头获得稳定的下落速度作用至关重要,探头尾翼虽然会减小下落速度,但有利于保证探头稳定的垂直下落姿态。研究结果验证了模拟计算方法的准确性,同时为同类探头的优化设计提供了参考依据。     

低速湍流中气动噪声预测方法

 为了研究气动噪声的产生机理和传播过程,在考虑介质黏性的影响情况下,采用分解法并结合湍流模型,对低速湍流流动中气动噪声问题在时域上进行数值计算。基于非结构化同位网格和有限体积法,把可压缩N-S方程分解成不可压N-S方程、含黏性项的声扰动方程;为了考虑质点振速和声压耦合,采用SimpleC算法来同步求解不可压N-S方程和声扰动方程。进出口远场边界采用以渐近解为基础的无反射边界条件,并采用与内部区域相对应的有限体积法、时间隐式格式对其进行求解。利用所编制的程序进行层流状态下圆柱绕流气动噪声仿真验证,并与文献结果进行对比,检验本方法的正确性;并结合湍流模型将数值解法推广到湍流状态下气动噪声数值模拟中。结果表明该方法能够很好地反映流场和声场的形态,无反射边界能很好地抑制声波在边界处的反射,适合低速流气动噪声问题模拟,为实际工程中的降噪工作提供预测信息。     

Higher-order Accurate and High-resolution Implicit Upwind Finite Volume Scheme for Solving Euler/Reynolds-averaged Navier-Stokes Equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｅｍｅｒｇｅｎｃｅｏｆＣＦＤａｂｏｕｔ15ｙｅａｒｓａｇｏｐｒｏｖｉｄｅｄａｍａｊｏｒｉｍｐｅｔｕｓｔｏｓｏｌｖｅｔｈｅＥｕｌｅｒａｎｄＮａｖｉｅｒＳｔｏｋｅｓｅｑｕａｔｉｏｎｓｇｏｖｅｒｎｉｎｇｔｈｅｆｌｏｗｆｉｅｌｄｉｎｅｘｔｅｒｎａｌａｎｄｉｎｔｅｒｎａｌｆｌｏｗｓ.ＭａｊｏｒｐｒｏｇｒｅｓｓｈａｓｂｅｅｎｍａｄｅｉｎｔｈｅｆｉｅｌｄｏｆＣＦＤａｌｇｏｒｉｔｈｍｄｅｖｅｌｏｐｍｅｎｔ,ｓｕｃｈａｓｍｕｌｔｉｇｒｉｄａｃｃｅｌｅｒａｔｉｏｎｔｅｃｈｎｉｑｕｅｓ,ｄｉ…     

全三维粘性流场计算的快速收敛方法

以新的高收敛率的 L U型隐式格式和高精度、高分辨率的 MU SCL TVD迎风格式为基础 ,提出了一种快速求解三维粘性流场的计算方法。为有效减少计算工作量 ,近固壁区采用了壁面函数方法。湍流模型采用简单的混合长度模型。为提高迭代收敛速度 ,从多重网格理论出发 ,结合改良型L U隐式和迎风格式 ,对 Reynolds平均三维可压缩 Navier-Stokes方程进行求解。用此方法数值模拟了 NASA L ewis37跨音速压气机转子流场。计算得到的流场与实验结果进行了对比。结果表明 ,此方法可快速得到三维粘性流场的流动特性且计算结果精度较高 ,可望在叶轮机械多叶片排全三维粘性流场校核计算中发挥作用。     

三维非定常粘性流动的数值方法及其应用

提出了适用于模拟非定常三维粘性流动问题的数值计算体系，包括高速收敛稳定性好的新的ＬＵ型隐式格式，捕捉弱间断面和滑移面的高精度、高分辨率的改良型ＭＵＳＣＬＴＶＤ格式。将基于此数值模型而发展起来的非定常粘性流动通用软件应用于叶轮机械气动热力学计算，为研究内流粘性损失、激波和非定常尾涡干涉、三维分离流动，尤其为探明分离泡内部的涡结构打下了基础。此软件亦适用于叶栅设计等方面的工程应用问题