圆柱跨声速绕流中的激波/湍流相互作用大涡模拟研究

采用大涡模拟方法研究了圆柱跨声速绕流中的激波/湍流相互作用问题,来流马赫数M∞取为0.75,基于圆柱直径D的雷诺数为2×105。计算结果表明,圆柱分离点处出现一道斜激波,并且以与涡脱泻Strouhal数一致的特征频率向上游传播。激波运动导致流场中出现反对称的流动模态,剪切层中压力信号的功率谱曲线中存在0.4、-1和-5次方的斜率关系,剪切层中的剪切应力角约为0°,脉动速度以流向脉动速度为主,并且沿剪切层的大尺度结构组织性减小。     

趋于临界马赫数的圆柱跨声速绕流特性分析

通过深化认识趋于临界马赫数Ma_cr的圆柱跨声速绕流特性,明确新型飞行器增升减阻设计的空气动力学理论依据。采用大涡模拟方法数值研究了来流马赫数Ma_∞为0.75和0.85、雷诺数Re为2×10⁵的圆柱跨声速绕流。结果表明:当Ma_∞趋于临界马赫数(Ma_cr≈0.9)时,圆柱的阻力下降且升力系数振荡被抑制;通过力的分解,得知圆柱的阻力减小来自旋涡力的影响,而非可压缩性;圆柱的阻力减小与其背压上升有关;剪切层初始阶段的对流马赫数Ma_c随Ma_∞的增加而增大,而增长率相反,这使得剪切层更为稳定、柱体背压更高。此外,由于Ma_∞=0.85时边界层分离点处的激波和尾迹处的激波向下游推移,使得近尾迹处的湍流脉动减弱,进而导致柱体的表面压力振荡和升力系数振荡被抑制。     

战斗机垂尾脉动压力数值模拟

在亚跨超计算流体力学（CFD）软件平台（TRIP）上开发了基于RANS/LES混合思路的IDDES流动模拟技术,并通过NACA0021翼型60°大迎角分离流动与串列圆柱绕流模拟对RANS/LES混合方法的精确度进行了验证,针对某战斗机外形的垂尾脉动压力开展了数值模拟研究。战斗机来流马赫数为0.1,基于全机长度的雷诺数为2×10⁶,模型迎角为20°、30°和40°。分别通过脉动压力系数、脉动压力功率谱密度、空间流动结构以及侧向力响应曲线等对战斗机的垂尾脉动压力进行了分析。脉动压力模拟结果表明：当垂尾完全沉浸在边条翼脱体涡破碎后的宽频湍流脉动气流中时,垂尾翼梢位置的脉动压力会发生明显的增大。     

一种二元进气道起动特性的数值与实验考察

采用Spalart-Allmaras湍流模型对来流马赫数变化引起的二元高超声速进气道模型的起动特性进行了三维数值模拟,并在来流马赫数5.5和5.9两种工况进行了激波风洞实验验证。结果表明,进气道不起动时,模型侧板上的流动分离与外压缩面上的流动分离相融合,堵塞了进气道入口,三维流动特征显著。在预报进气道自起动和不起动时,应当考虑三维流动效应。逐步增大来流马赫数和逐步减小来流马赫数两种路径下,进气道存在起动迟滞现象。对于本文的进气道模型,当内收缩比增大0.2时,自起动马赫数约升高1,不起动马赫数约升高0.5,进气道的起动迟滞环变大,起动性能下降。在激波风洞中获得的进气道自起动以及起动/不起动双解区的实验结果与数值模拟相符合。     

跨声速半球柱绕流流场的激波系结构

本文在1相似文献   

伴流速度对平行喷口射流影响的数值研究

采用大涡模拟方法计算研究平行喷口出口马赫数为0.9,伴流速度比分别为0.1、0.3和0.5时的喷射流场特性。计算时采用高精度的数值模拟方法,并结合Smagorinsky亚网格尺度模型。考察了流场统计平均特性、脉动特性以及射流流场中涡结构的发展演变过程,结果表明:伴流速度的增大使得势流核长度变长,减缓了空间剪切层的发展,转捩延迟。喷射流场速度分布具有自相似性,而湍流强度的分布则不具有相似性。通过分析剪切层中轴向速度脉动、径向速度脉动、压力脉动在空间任意两点上的时空相关性,发现随着伴流速度的增大,脉动量在空间上的相关性减弱,而脉动量向下游的传递速率增加。该研究结果为进一步揭示伴流速度对喷流声场的影响提供基础。     

高马赫数低雷诺数条件下圆球绕流曳力系数

圆球绕流曳力系数通常表示为雷诺数的函数.而随着马赫数增大,曳力系数与马赫数的相关性逐渐增强.本文应用硬球-拟颗粒模型方法模拟研究1≤Re≤20和0.1≤Ma≤3条件下的圆球绕流过程.分析了模拟区域设置和拟颗粒-圆球作用的热适应系数对曳力系数结果的影响,得到合适的区域设置条件和热适应系数,获得高马赫数低雷诺数条件下圆球绕流曳力系数,并依此初步建立1≤Re≤20和0.1≤Ma≤3范围内马赫数对曳力系数的修正模型.     

细长旋成体背风面非对称流场的DES模拟

利用脱体涡模拟DES方法对细长旋成体非对称绕流进行数值模拟,比较了基于S-A湍流模型的DES与RANS方法对背风面非对称分离涡的数值模拟能力,发现在大迎角非对称涡未破裂情况下,DES与RANS均能模拟出与实验相符的非对称绕流;当迎角增大,背风面非对称涡发生破裂时,RANS无法准确捕捉到背风面流场的非定常性,而DES能准确预测非对称涡的飘起与破裂,并与实验值接近。计算结果表明:与RANS方法相比,DES方法具有更好的模拟大攻角,大分离流动的能力,尤其在非对称流场的分离涡破裂模拟方面具有明显优势,能够更真实模拟出细长旋成体背风面分离涡破裂之后的非定常流动特征。     

对涡旋流影响压气机转子性能和稳定性的研究

为研究对涡旋流影响跨声速轴流压气机性能和稳定性的机理,设计了一种叶片式旋流发生器,并对旋流发生器和跨声速单转子进行了联合数值模拟研究。计算结果表明旋流发生器叶片数越少,对涡强度越低;对涡旋流导致压气机总压比、峰值效率、稳定工作范围和堵塞边界流量减小,失速边界流量增大,对涡强度等于60°时峰值效率和堵塞边界流量分别降低1.11%和2.12%,失速边界流量增加4.17%;对涡中的同向涡导致叶尖进口攻角增加,进口相对马赫数减小;反向涡使叶尖进口攻角降低,进口相对马赫数增大;对涡前缘的轴向速度偏低,造成叶尖进口攻角大幅增加,叶尖泄漏流堵塞严重;由于叶尖泄漏流在叶片前缘溢流导致失稳。     

屏蔽式总温探针流动与换热分析

采用定常/非定常数值模拟方法研究了某屏蔽式总温探针的流动与换热特点,讨论了来流马赫数、进出口面积比对其稳态误差的影响,给出了不同条件下屏蔽式总温探针稳态误差的变化规律。结果表明,在来流马赫数0.5条件下,与仅考虑对流换热相比,导热与辐射使屏蔽罩表面最高温度降低约19K;气流在屏蔽室的头部和中部区域形成两个主要的涡区,其中中部的肾型涡尺度最大,是屏蔽室里的主控涡;在高来流马赫数条件下,速度误差和导热误差对进出口面积比的变化非常敏感,总误差随进出口面积比的增大而迅速增大;对于本文研究的探针,在Ma=0.6时,出口半径从1.5mm减小至0.5mm,探针总稳态误差将增大超过12K。     

叶轮机械新流型探索雏议

从“流型”的角度, 对比飞行器空气动力学和叶轮机械气体动力学的发展历程, 认为叶轮机械内的流动仍在追求定常附体流型, 由于其流动及结构特性, 使得“有害”分离问题难以避免。本文提出应探索将外流领域内的定常脱体涡流型或下一代非定常脱体涡流型引入到叶轮机械的可能性与技术途径本文并综述了上述探索性研究的有关背景和现状, 以及所设想的基本途径和影响因素。      

平板扰流片振动诱导涡运动研究

 本文研究了平板上扰流片静止及振动时所诱导的分离涡的形成过程。当拢流片静止时。在拢流片后仅形成不规则的大尺度涡,一旦执流片振动。就诱导出有固定频率的分离涡。拢流片振动频率,对分离涡的尺度、强弱、成对、合并过程有重要影响。文中还对的流片振幅大小的影响进行了讨论。     

超音速绕凹角湍流分离流动的计算

本文用统一的Levy-Lees变换以及正算法与逆算法相结合,求解了超音速绕凹角湍流分离流动。 对附着流区用边界层正算法,压强分布用流过尖劈统一的高超音速与超音速公式,湍流模型取代数涡粘性模型;对凹角分离区用边界层逆算法,给定位移厚度δ~*分布,湍流模型取代数松弛模型;边界层计算采用Cebeci-Keller Box方法;计算成功地算得分离流场,较好地预估了分离点与重附点位置以及壁面压强分布与表面摩擦应力分布。     

跨音速叶栅粘流计算的多级广义有限元法

将 Taylor-Galerkin广义有限元法和多级有限元的思想结合起来 ,并在人工粘性的处理上作了改进 ,构成了收敛速度和稳定性均较好的多级广义有限元法。利用该方法 ,并基于 N-S方程研究了二维跨音速叶栅流动。计算结果与实验结果符合较好。通过计算表明 ,该方法计算稳定、收敛较快 ,是叶轮机械内部跨音速流动强有力的计算手段     

凹半球降落伞模型在加速和定常流中的流场观测

 在北京航空学院水槽中,用氢泡法对绕凹半球模型流动进行了研究。模型前后流速分布是用氢泡时间线测量的,它与在起动开始时由无旋流理论估算值接近,但当起动旋涡变大以后,在尾流中流速差别变大。在模型后部诱导出较大的反流速席。所以由无旋估计的虚质量只在开始时是可用的,并且在起动涡变大时,虚质量应更大。 在定常流中,在模型前有时形成一个大旋涡,这个旋涡是非定常的并且可能引起侧向力和伞的非定常运动。模型上适当开孔可减小甚至消除此旋涡。     

利用机上传感器测量燃油流量的方法研究

主要介绍了在小型商用飞机无法加装流量传感器的情况下,利用机上传感器实现流量测量的方法,通过该方法完成了该飞机在各种状态下燃油流量的测量。     

非定常N-S流中的瞬时临界点

本文用微分方程的相似理论把非定常轨线方程变换成自洽系统,考察变换后临界点邻域在非定常流场中的映象,进而对非定常流中“瞬时临界点”溉念的可用性及其限度作了讨论。     

LTRAN2的一种改进方案

本文提出了LTRAN2的一种改进方案BTRAN2用来分析翼型的各种频率的非定常运动。用E-O调转换格式和ADI方法求解了完全的二维非定常跨音速小扰动位流方程,并用单调转换的AF2格式计算了定常跨音速小扰动方程,以此作为非定常计算的初场。本文给出了带后缘正弦振荡襟翼的NACA 64 A006翼型的绕流和做正弦俯仰振荡的NACA 64 A010翼型的绕流计算结果,它们与Euler方程解或实验数据很吻合。     

文氏效应在叶片根部马蹄涡控制中的应用

为研究文氏管在叶片根部马蹄涡控制作用的应用,针对简化了的叶片-平板结构在叶片上游平板上的某位置开狭缝,在平板的另一侧构造一文氏管(文氏管的喉道位于狭缝处),引射叶片根部的流体,进而控制马蹄涡.数值模拟了不同狭缝位置的湍流流动;结果表明了改变狭缝相对于叶片的位置能很好地控制马蹄涡,甚至能完全消除叶片前马蹄涡.提出的方法可以利用外势流的能量而无需额外消耗能量,为一种被动的控制手段,有很大的实际应用前景.利用文氏效应的控制方法简单实用,因此工程中采用该法是不错的一种选择.      

飞机亚跨超绕流的数值模拟研究

本文以Euler方程为数学模型,采用一种高精度的TVD(Total Variation Dimishing)离散格式及一种含近似因式分解的推进迭代方法,求解亚跨超绕流’流场。通过若干算例的试算,证明方法是可行的,可以用来模拟飞机的复杂流场。