提高k-ωSST模型对翼型失速特性的模拟能力

采用CFD软件fluent,k-ω SST模型对S825翼型进行了二维数值模拟.针对k-ω SST模型对翼型分离失速特性模拟不准的情况,对分离区域内湍流的强非平衡输运特性进行了分析.研究表明:原始模型由于不能准确模拟分离区内湍流的强非平衡输运特性,导致其对翼型失速特性模拟失效;提出了通过修正模型系数a₁及β_∞^*,提高k-ω SST模型预测非平衡湍流输运特性的方法,从而提高对翼型失速特性的模拟精度.     

高超声速二方程湍流模型的数值模拟对比

对基于Reynolds和Favré(密度加权)混合平均的二方程湍流模型进行了修正,同时根据数值模拟高超声速流动时必须具有高分辨率捕捉间断面与在边界层内抑制数值粘性的能力的要求,提出了新的总变差减小(TVD)格式熵修正函数.在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对基于Reynolds和Favré混合平均的二方程湍流模型,以及其它不可压缩模型及可压缩性修正模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了对高超声速压缩拐角型流动,湍流模型可压缩修正的必要性,得到了基于Reynolds和Favré混合平均的二方程湍流模型的计算结果最接近实验结果.      

多块网格计算跨音压气机带叶尖间隙流场

开发了三维数值模拟程序研究轴流跨音转子叶尖间隙流动,应用高雷诺数k-ε湍流模型加壁面函数的方法,计算了轴流跨音转子NASA Rotor37在设计转速下的流场.叶尖间隙采用分区的H型网格和主流区连续对接耦合计算,没有用间隙模型,也没有考虑Vena收缩效应而减小间隙量.在用有限体积法对Navier-Stokes方程和湍流方程进行空间离散的过程中采用了交错网格的方法将N-S方程与湍流方程紧密地耦合在一起,从而提高了计算精度.计算结果和实验数据进行了详细的比较和分析.结果表明,中部叶展具有与实验结果非常一致的流场特征,根、尖区流场则因涡粘假设和激波问题的存在而使流动细节与实验结果略有偏差.     

某高空螺旋桨气动特性数值模拟与风洞试验

以某高空螺旋桨为研究对象,利用k-ε湍流模型,采用三维Navier-Stokes方程,并通过数值计算模拟该螺旋桨的复杂绕流,分析其气动特性.按照等雷诺数等前进比准则,进行了该螺旋桨地面缩比风洞试验.经过比较发现,计算结果与试验结果吻合良好.验证了数值模拟计算高空螺旋桨气动特性的正确性,为高空螺旋桨的设计和试验提供参考依据.     

近壁区湍能耗散率输运方程的理论模型

提出在湍流边界层的近壁区采用三维不稳定波来描述湍流扰动速度,然后根据理论模型对湍能耗散率的生成、粘性破坏以及扩散等特性进行了系统定量的计算和分析.理论计算结果与直接数值模拟结果符合较好,特别是在y<10的区域中明显优于传统湍流模型的计算结果.这不仅在理论上有益于对湍流物理机制的了解,而且可能为湍流模型计算开辟一条新的途径.     

使用GAO-YONG湍流方程组对扩压器流动的计算

采用基于同位网格的SIMPLE方法求解GAO-YONG不可压湍流方程组,对二维扩压器流动进行了数值模拟.通过界面速度动量插值法解决了压力锯齿波问题,并且采用了正交贴体网格和二阶QUICK格式离散对流项以提高计算精度.与实验结果以及BL(Baldw-Lomax)模型计算结果的比较表明,不需要任何经验系数及壁面函数的GAO-YONG不可压湍流模型方程组能够对有压力梯度的湍流流动做出很好的预测.计算发现,在机械能方程中引入平均流压力梯度的作用,对GAO-YONG湍流方程组正确模拟逆压力梯度流动起到了关键作用.      

均匀各向同性湍流的脱体涡数值模拟

采用脱体涡模拟方法对均匀各向同性湍流进行了数值模拟,并与Comte-Bellot实验测量结果进行了对比,验证了该文的脱体涡模拟方法对均匀各向同性湍流模拟的可靠性.初始速度场的生成采用Rogallo所提出的构造方法,初始湍能谱满足Von Karman波谱分布;在对流输运项的选择方面,分别采用二阶中心型格式、四阶偏斜对称型中心格式和迎风型低耗散通量分裂格式,考察它们在均匀各向同性湍流模拟中的计算精度和适用性.同时,通过改变计算域大小以及脱体涡模拟方法中的模型常数,达到能谱截断波数的改变,考察它对各向同性湍流计算的能谱以及能谱截断处小尺度涡能量积累问题的影响.      

飞行器表面三维流场与固壁温度场的耦合分析

分析了飞行器表面温度场形成的物理机制,建立了飞行器蒙皮温度场耦合计算的理论模型,较为完整的考虑了各种热源对飞行器蒙皮温度的作用.采用有限体积法和重整化群RNG(Renormalization Group)k-ε湍流模型对流场求解N-S方程,蒙皮固壁中采用热传导模型,采用离散坐标法DOM(Discrete Ordinates Method)计算蒙皮固体壁面对外界大气的辐射.通过数值计算研究了飞行高度、马赫数和内热源对蒙皮温度分布的影响,以及高度及马赫数对飞行器表面气动加热的作用.结果表明,该方法适合考虑飞行器外部流场、固壁等多种热源作用的温度场数值模拟.     

DES方法对返回舱气动热环境数值模拟

脱体涡模拟方法(DES,Detached Eddy Simulation)由于结合了湍流模型(RNAS,Reynolds Averaged Navier-Stokes)和大涡模拟(LES,Large Eddy Simulation)两者各自优势,是模拟非定常大分离流动的有效方法.采用基于SA(Spalart-Allmaras)湍流模型的DES方法,对高超声速返回舱外形进行数值模拟,计算所得的返回舱表面压力及热流分布与实验结果吻合,证实了DES方法的优越性.最后,在DES方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.     

DES方法对返回舱气动热环境数值模拟

脱体涡模拟方法(DES,Detached Eddy Simulation)由于结合了湍流模型(RNAS,Reynolds Averaged Navier-stokes)和大涡模拟(LES,Large Eddy Simulation)两者各自优势,是模拟非定常大分离流动的有效方法.采用基于SA (Spalart-Allmaras)湍流模型的DES方法,对高超声速返回舱外形进行数值模拟,计算所得的返回舱表面压力及热流分布与实验结果吻合,证实了DES方法的优越性.最后,在DES方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.     

曲率对旋转态气膜冷却效率影响的数值模拟

通过对旋转状态下曲率叶片模型上气膜冷却现象的流动和换热进行数值模拟,得到了不同主流雷诺数、吹风比和旋转数情况下吸力面和压力面上的冷却效率分布.计算选用κ-ω和SST(Shear-Stress Transport)湍流模型,主流雷诺数Re=3 198.4~6 716.6,吹风比M=0.2~1.2,旋转数Rt=0~0.015 9.结果表明:旋转数的增大导致气膜孔下游中心区域的冷却效率下降,但使压力面整场的冷却效果略有提高;吹风比的增大使得吸力面和压力面上的冷却效率逐渐降低,主流雷诺数的变化对壁面整体冷却效果则影响不大.此外,相同工况下吸力面上的冷却效率要高于压力面上的对应值.     

高超声速流动湍流模式评估

选取了超声速二维膨胀-压缩拐角和高超声速双椭球绕流作为基准流动,考察了几种当前CFD(Computational Fluid Dynamics)工程应用较流行的湍流模式:BL(Baldwin-Lomax)模式、SA(Spalart-Allmaras)模式、k-ω模式以及SST(Shear-Stress Transport)模式,通过数值计算结果和实验结果的对比及分析,对有关的湍流模式进行了评估,得到一些有意义的结论,并指出受流动可压缩效应的影响,计算高超声速湍流壁面热流时要对现有的湍流模式进行相应的可压缩修正.      

改进SA模型对翼型分离流动的数值模拟

采用(CFD,Computational Fluid Dynamics)软件Fluent,选取(SA,Spalart-Allmaras)模型,对风力机翼型S825进行了二维数值模拟,并与实验数据进行比较.针对SA模型不能准确预测翼型尾部流动分离的情况,分析了分离区域内湍流的强非平衡输运特性.研究表明,SA模型中系数C_b1与湍流的非平衡输运特性密切相关,进而提出了修正模型系数C_b1改进SA模型的方法,结果表明该方法能够更准确地模拟风力机翼型分离流动.     

基于空间分层组合设计的火箭落点实时计算模型

为保证航天发射任务航区安全,实现火箭落点的高精度实时快速计算,提出一种基于空间分层的数学模型,对落点计算几何模型与动力学模型进行组合设计.该模型将弹道空间分层为自由段和再入段,自由段采用地球扁率修正的落点计算几何椭圆模型,而再入段采用考虑了空气阻力的动力学模型,并加入地球摄动及旋转等因素的修正.经与实际落点数据对比验证...     

前飞旋翼三维湍流场的数值模拟

采用计算流体力学CFD(Computational Fluid Dynamics)方法,分析正常前飞和近地前飞时直升机的流场与性能.在直角坐标系中,采用有限体积法和压力耦合方程的半隐式法SIMPLER(Semi-Implicit Method for Pressure-Linked Equations Revised)算法求解三维定常不可压的湍流N-S(Navier-Stokes)方程.在分析中,旋转的螺桨被描绘成沿螺桨桨叶展向分布、与本地流动参数相关以及时间平均的动量源项.计算方法还包括K-ε湍流模型和壁面函数法等措施.流场分析和性能预测同实验测量数据的良好的一致性表明,这种CFD方法可以有效的分析直升机的具体设计问题.      

带肋变截面回转通道内流动与换热的数值模拟

开发了三维流动换热的通用计算程序,数值研究了带肋变截面回转通道内流动与换热的特性.湍流模型采用低雷诺数k-ε模型.通道肋间距为25mm,肋高分别为1mm,1.5mm,2mm,冷气进口雷诺数Re分别为7500,12500,18500,25000.计算结果表明:①通道的平均努赛尔数均随进口雷诺数的增大而增大;②对于Re=7500和12500,肋高越高,换热越强;对于Re=18500和25000,肋高为1.5mm的通道换热最强;③局部雷诺数的不同和离心力的影响导致通道内各区域的局部换热随肋高的变化趋势并不一致;在进口段,肋高越高,换热越强;在出口段,当Re=7500和12500时,肋高越高,换热越强,而当Re=18500和25000时,存在最佳肋高1.5mm.     

旋转状态下曲率对气膜与主流掺混区域的影响

气膜冷却是应用于航空燃气轮机上的冷却技术,旋转及表面曲率是影响气膜与 主流掺混区域的重要因素,通过数值计算方法对旋转状态下曲率对气膜与主流掺混区域的影 响进行了研究,湍流模型选取了k-ω模型.增加旋转速度,会引起吸力面气膜的分离; 固定转速,降低表面曲率半径,压力面气膜发生分离,吸力面气膜冷却效果得到改善.当动 量流量比在小于1的范围内变化时,旋转只改变压力面气膜与主流掺混区域的分布,而对吸 力面没有影响.      

用GAO-YONG湍流模型计算激波/边界层干扰

采用GAO-YONG可压缩湍流方程组,模拟了平板激波/湍流边界层干扰现象.运用SIMPLE算法求解方程组,并分别采用三阶精度的QUICK格式和中心格式离散对流项和扩散项.计算结果较好预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征: 分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波.对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布和速度型与实验值比较吻合很好.结果表明GAO-YONG可压缩湍流方程组能够高精度模拟平板激波/湍流边界层干扰流动.      

适用于UH模型的单元破坏修正方法

在采用UH模型的用户材料子程序(UMAT)进行有限元计算过程中时常会出现局部单元破坏的情况,主要包括拉裂和剪坏,这种破坏应力状态的存在不仅会使得计算结果不合理,也会降低有限元计算的稳定性。针对采用UH模型UMAT进行有限元计算时遇到的局部单元破坏问题,根据一定的假设条件,并结合不同坐标系下的应力变换关系,推导得到适用于UH模型的三维问题单元破坏修正公式;采用FORTRAN语言编写出相应的子程序,将其嵌入UH模型UMAT中,以消除采用UH模型进行有限元计算时出现的不合理的破坏应力状态,提高UH模型UMAT有限元计算的稳定性,并通过有限元模拟基坑开挖的例子验证本文所提的单元破坏修正方法的有效性和合理性。     

飞机装配坐标系公共基准点粗差检测与修正方法

针对因厂房地基下沉等因素引起的飞机数字化装配体系中公共基准点产生偏差的问题,提出一种飞机装配坐标系公共基准点粗差检测与校正方法.该方法通过泰勒展开布尔沙-沃尔夫空间坐标转换模型形成了适用于大旋转角的坐标系转换线性模型,并结合拟准检定法检测并估算偏差较大的基准点,进而给出了公共基准点在装配全局坐标系下的坐标修正模型.采用激光跟踪仪对公共基准点进行测量并用上述方法修正公共基准点坐标.实验证明,该方法能有效检测粗差并修正,提高飞机装配坐标系的精度.