数值模拟有外挂物的空腔流动

数值模拟了有外挂物存在的开式空腔的超音速流动。数值计算采用全隐式无分裂方法,湍流模型采用修正的B-L代数模型。研究了有外挂物存在的空腔流动的剪切层自持性振荡,并与无外挂物存在时的现象进行类比。压强与时间的变化曲线、空腔底部流向的静压系数分布以及空腔后壁面法向的静压系数分布表明,空腔内的压力脉动与外挂物所处位置是有关的。对于Rossiter振荡模态,数值计算结果能够与实验结果相吻合。     

高升力系统外形的数值模拟计算

分析和讨论了高升力系统外形无黏和黏性流动的物理特征,表明了流动现象的复杂性和求解的艰难性。深入研究了二维和三维高升力外形空气动力数值模拟计算的雷诺平均N-S(RANS)方程方法、计算结果、以及它们和实验数据的比较。二维结果表明压强分布、速度型等计算值和实验值吻合很好,但最大升力系数两者仍有差异。三维计算最大升力系数的结果与二维类似。还分析了影响计算结果的各种因素,包括计算网格、湍流模型、外形的几何模拟的逼真度等。特别指出了转捩模型研究对提高计算准确度的重要性。文中解释了研究高升力系统外型的数值模拟方法,以及同时发展能估算C_Lmax的工程半经验方法的重要性和迫切性。     

腔内平板对空腔自激励振荡的影响及预估振荡频率方程的改进

运用全隐式无分裂方法,紊流模型采用修正的B-L代数模型,数值模拟了超音速流通过腔内存在平板的空腔流动.研究了腔内有不同长度平板存在的空腔流动的剪切层自持性振荡,并与无平板存在时的现象进行类比.腔内平板将剪切层与空腔隔开,减弱了剪切层振荡与空腔压强脉动的耦合作用;随着腔内平板长度的增加,剪切层振荡的程度有所减小.基于特征长度可变的思想,本文改进了Heller的频率预估方程,改进后的频率方程能更好地预估空腔的振荡频率.与Heller的方程相比,修改后的方程预估所得到的频率值更接近于闭式空腔声学模型方程的预估结果;各个方程预估空腔振荡频率值的比较进一步验证了改进后的方程的正确性及合理性.     

导弹垂直发射系统三维超声速射流的数值模拟

对导弹垂直发射系统(VLS,Vertical launching system)内的三维湍流流场进行了数值模拟.采用了单元中心有限体积法和三阶迎风TVD差分格式离散Navier-Stokes方程,采用双时间步和预处理过程提高计算收敛速度,引入对角化的双时间ADI算法求解控制方程.仿真结果准确地模拟了发射系统内的激波形成过程,真实地反映出VLS中的复杂流动.研究表明,压力室底部射流中心的波后压强以大约1 200Hz的频率振荡.      

直接供气预旋转静系流动和换热数值模拟

采用CFD商业软件FLUENT分别从维度模型和湍流模型两个角度对接近真实的高压涡轮旋转盘腔进行数值模拟,研究了一种直接供气预旋转静系旋转盘腔模型内的流动和换热特性.通过计算研究发现:低位腔内的流动仅由中心入流控制,二维计算不能反映高位腔内中心入流与预旋入流的掺混;二维和三维计算所得的腔内压力分布相似,实验值介于两者之间;主盘面的换热主要受中心入流控制,二维和三维计算差别不大,Realizable k-ε湍流模型在定性和定量上均能更好地反映换热实验结果.      

模型固体火箭发动机内涡脱落过程的大涡模拟

为研究后向台阶结构引起的涡脱落对固体火箭发动机内流动稳定性的影响，采用大涡模拟技术，对模型发动机内的流动进行了数值模拟，获取了发动机燃烧室内压强随时间的变化曲线及其振荡频率，并和已有的计算结果及实验结果进行了对比分析。结果表明采用LES方法的计算结果明显优于采用k-ε湍流模型的计算结果；不考虑燃烧时模型发动机内由涡脱落引起的压强振荡的振幅较小，但频率较高，且随来流速度的增加而增大；在后向台阶这种结构中，涡一般会在台阶的边沿形成，并在其下游周期性地脱落；在保证y^＋较小的情况下，网格大小对计算结果的影响不明显。     

空腔流动特性数值模拟研究

本文采用数值方法模拟研究了空腔的长深比（L／D）对流动特性的影响。控制方程采用Navier-Stokes方程，计算方法采用有限体积法，湍流模型采用Baldwin-Lomax模型。通过分析三种典型空腔的流线图谱、压力云图、底部压力分布曲线，得到空腔流动特性与其几何参数和来流参数间的相应关系。从数值模拟的结果来看开式空腔流动的空腔内压力分布均匀，流场结构相对稳定，有利于内埋式武器安全平稳地分离。     

考虑空腔的高超声速多流动区域同步数值模拟

先进的高超声速飞行器具有薄壁空腔结构,在飞行过程中受热会产生空腔内气体流动现象,从而影响流场和结构的温度分布。采用数值方法准确模拟高超声速流场、结构温度场和空腔内流动对热结构分析是很有必要的。以研究空腔流动对结构温度分布影响为目的,发展了一种适用于多流动区域流场/结构温度场耦合问题的同步计算方法,并以高超声速带空腔结构物体为例,数值研究了其外部气动热/结构热传导引起的空腔热对流问题。以已发展的高超声速外流场/结构温度场同步计算方法为基础,为了进一步考虑空腔内低速流场,采用了预处理矩阵方法。在流场与结构温度场的交界面两侧分别引入虚拟单元,从而高效地实现相邻场之间物理信息交换。首先通过标准算例验证了方法在求解单独气动热/结构热传导问题以及空腔自然对流问题中的准确性。进而对封闭和带有开孔的两种高超声速运动圆环分别进行多流动区域同步数值模拟。计算结果表明,由于结构温度不均匀引起的空腔内热对流反之也会对结构温度场分布产生轻微的影响。在空腔内气体流动的影响下,封闭圆环的前缘温度在35s内最多下降0.8%左右。对于带开孔空腔的圆环,其孔壁周边温度在0.5s内能够超过外流前缘驻点温度。     

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

低雷诺数翼型层流分离现象大涡模拟方法

为准确预测翼型低雷诺数条件下出现的层流分离流动现象,发展了基于结构化拼接网格的大涡模拟方法。控制方程为Favre滤波Naver-Stokes方程,并选用多种亚格子模型。空间离散采用AUSM格式以及高阶WENO格式,时间推进采用显式方法和隐式方法。以SD7003翼型在雷诺数60 000及4°迎角下的层流分离流动为研究对象,对比分析了数值格式、亚格子模型、网格尺度对流场预测结果的影响。数值结果表明:划分的计算网格能够有效解析小尺度流动结构,基于隐式亚格子模型、采用AUSM格式和双时间步推进的大涡模拟方法能够准确预测流动分离、转捩、再附等复杂流动现象,计算得到的平均压强系数与雷诺应力分布与文献数据吻合较好。与转捩模式计算结果对比进一步表明:发展的大涡模拟方法能够准确预测翼型低雷诺数层流分离流动非定常演化过程,为下一步的研究工作建立了有效的数值模拟手段。     

超声速空腔流动数值模拟研究

采用Roe格式、S-A一方程湍流模型和LU-SGS隐式算法对三个不同类型的空腔进行了模拟,并把数值计算结果和试验结果进行了比较,二者一致.通过对流过不同L/D空腔的可压缩流动进行模拟,得到了不同L/D空腔的底部压力分布、空腔内流线、等压线和等马赫数云图,分析了空腔流动类型随L/D增大转变的机理.     

超燃冲压发动机不同燃烧室结构流场数值模拟

对带不同长深比凹腔结构（燃烧室）的二维超燃冲压发动机流场进行数值模拟,采用迎风二阶精度NND格式求解二维欧拉方程,对流场进行冷态数值计算,能够更加清楚地捕捉、分析激波的特点,反映流场的状态和旋流特点;结果证明凹腔可以使流场中形成低速高压回流区,以达到燃料混合燃烧的目的。凹腔长深比对流场特性、总压恢复、燃烧室阻力影响显著。随着长深比的越大,凹腔旋流强度越强,总压损失越大,燃烧室阻力越大。     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

三维不可压流动有限元数值模拟

用原始变量带涡量的伽辽金有限元法计算三维定常不可压Ｎ－Ｓ问题，克服了过去由于速度和压力导数不同阶而不得不采用不同阶插值函数的困难。引入牛顿法线化Ｎ－Ｓ方程，使得用有限元离散后的大型代数方程组线性化，用高斯－赛德尔迭代法求解。算例为三维方穴驱动和带潜入锥的固体火箭发动机燃烧室轴对称加质流动。三维方穴驱动的计算结果与文献值相近。     

长深比对空腔流动与声学特性的影响分析

长深比(L/D)是影响空腔流动的最重要外形参数。为了对空腔的流动机理有更深入的认识,采用基于Realizable k-ε湍流模型的脱体涡模拟(DES)数值方法,研究空腔L/D的变化对空腔流动与声学特性的影响。选取的来流马赫数为0.85,基于空腔长度的雷诺数为6.84×106。结果表明:随着L/D的增加,空腔前部的压强系数逐渐降低,后部的压强系数逐渐升高,导致空腔内的静压梯度逐渐增大;整个空腔内的压强脉动强度减弱,噪声环境有所改善。研究结果可为空腔流动类型的划分提供一定的依据。     

烧蚀外形三维粘性流场数值模拟

本文通过求解层流Navier-Stokes方程,数值模拟了三维奶头型烧蚀外形粘性流场。差分格式为Yee和Harten的迎风隐式TVD格式,在空间方向具有二阶精度。计算得到了高质量的激波和合理的流场结构。物面压力分布与实验结果进行比较,符合较好。     

不同去旋角度进气共转盘腔内流动与换热研究

以带有去旋喷嘴的共转盘腔为研究对象,运用RNG k-ε湍流模型,分别对3种喷嘴角度（90°、60°、30°）的3维模型进行了数值计算,揭示了去旋角度对盘腔内气体切向速度、盘腔内压降和转盘表面换热的影响规律。研究结果表明：随着去旋角度的减小,盘腔内气体的切向速度减小,盘腔内压降减小,转盘表面的局部努赛尔数Nur增大。     

二维高亚声速空腔流激振荡的数值模拟研究

采用基于P .L .Roe的近似Riemann解的修正Osher&Chakravarthy (MOC)三阶TVD有限差分格式 ,数值求解二维雷诺平均全Navier Stokes方程 ,并用Cebeci Smith代数湍流模型 (对腔内区域作修正 )来模拟湍流效应 ,时间方向的积分采用四阶Runge Kutta方法 ,通过对GAMM超音速前台阶绕流的计算 ,验证了格式及程序的有效性 ,对高亚声速来流下的空腔流动作了数值模拟研究 ,取得了较好的效果     

三维高超声速喷流干扰流场的数值模拟

采用ＬＵ隐式方法求解了有限体积法离散的完全Ｎ－Ｓ方程,数值模拟了三维高超声速喷流干扰流场。计算中采用了Ｊａｍｅｓｏｎ的当地极值递减（ＬＥＤ）格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,与实验结果相比较,数值方法对附体流动计算较准确,对分离区附近流场计算有一定误差。     

含侧喷流弹丸三维绕流干扰流场数值分析

从全Navier-Stokes方程出发，应用三维隐式有限体积TVD格式，在超声速和零攻角条件下，分别对无侧喷流和含侧喷流的低阻远程弹丸三维粘性绕流流场进行了数值模拟，为避免分区计算而用代数方法生成了弹丸绕流场O型网络，计算结果得到了清晰的流场波系结构，并从流场波系结构、弹体和弹底表面压力分布等角度将两者进行了对比与分析，表明侧喷流不仅提供增程推力，而且改变了弹丸阻力特性，其结果可用于研究弹丸的火箭增程技术。