特型钝体绕流动态特性的低速实验研究

在低速自由射流风洞上 ,利用热线及相关设备对由半球头部、颈部及旋成体构成的特型钝体周围非定常流动的基本特性、联合特性及湍流度分布进行了测试与分析 ,实验工况为迎角α=- 1 5°～ 1 5°、侧滑角 β =0°～ 1 0°。结果表明 ,在上述工况范围内 ,该钝体周围流场中速度脉动的能量分布平坦 ,属宽频带随机信号 ,流场中没有发生明显的流动分离 ,流场动态品质良好 ;凹陷区内气流的脉动以不同的速度向下游空间传播 ;在上述α、β变化范围内 ,钝体颈部凹陷区的最大湍流度高达 1 0 .9～ 2 4.8%。     

有侧滑角椭圆弹身大攻角状态绕流的数值模拟

用数值模拟方法研究了椭圆细长体在大侧滑角、大攻角状态下的流动特性。数值模拟的出发方程为三维可压缩流的全Ｎ－Ｓ方程，数值格式为ＴＶＤ格式。通过数值模拟，研究了椭圆截面细长体在Ｍ∞＝２．５、侧滑角β＝４５°状态下，攻角从１０°￣３５°范围内流场中的分离特征。     

应用GAO-YONG湍流模式数值模拟三维激波/湍流边界层干扰

应用GAO-YONG可压缩湍流模式数值模拟了三维激波/湍流边界层干扰算例之一——单鳍流动.攻角20°,来流马赫数2.93,雷诺数9.8×105.对流项和扩散项分别采用Roe格式和二阶中心差分格式计算.Runge-Kutta显示时间推进方法求解了半离散的控制方程.包括壁面压力分布,边界层内流动偏移角等在内的计算值与试验数据进行了比较.准确地预测出了三维激波/湍流边界层干扰流场的主要流动特性——λ波结构,主分离涡核,膨胀区,滑移线等.计算与Alvi等提出的单鳍流动的理论模型符合很好,得到了平板表面压力以及分离线、再附线等在单鳍流动中所独有的半圆锥特性.      

熵层分析与钝体绕流

本文对钝头旋转体的流动特点进行了分析, 在较严密的量级分析基础上讨论了熵层流动特点, 提出了简单的流动模型和计及熵层效应的计算方法。计算了球头、椭球头等两种钝锥绕流的物面流动参数、流场及小攻角气动力系数, 与特征线法和实验结果的比较都具有很好的定性趋势与定量的结果。根据理论分析和给出的计算结果对熵层特点进行了讨论, 另外对钝锥体气动力特性也作了一些分析与讨论。     

主流攻角对双射流孔气膜冷却特性的影响

应用压力敏感漆技术,在平板上测量了不同主流攻角(i=-30°,-20°,-10°,0°,10°,20°,30°)下双射流孔的气膜冷却效率,并利用计算流体动力学(CFD)计算得到的流场对气膜冷却效率的规律进行了解析。所研究的双射流孔结构的孔间无量纲横向距离为0.5,孔间无量纲流向距离为3;射流与主流密度比为1.0,吹风比分别为0.5、1.0、1.5、2.0。结果表明小的主流攻角(i=-10°,0°,10°)下,流场中存在反肾型涡对或挤压作用,气膜层与壁面贴附良好,气膜冷却效率最高;大正值攻角(i=20°,30°)下,虽然气膜覆盖面积大,但反肾型涡对退化,气膜冷却效率下降;大负值攻角(i=-20°,-30°)下,流场中有肾型涡对,且气膜横向覆盖受限,气膜冷却效率最低。     

细长体大攻角绕流的数值模拟

用数值模拟方法研究了细长体大攻角的绕流特性（求解在周向３６０°范围内进行）。数值模拟的出发方程是三维可压流的全Ｎ－Ｓ方程,格式为二阶迎风ＴＶＤ格式。通过数值模拟研究了附近,攻角１９°、３２．５°和５０°状态下,细长体绕流场中的激波、分离流动及表面压力分布等重要特征。     

高精度差分格式在机翼分离流计算中的应用

用一个模型方程分析了边界层方程的数值稳定性 ,指出稳定性问题会随计算雷诺数减小而变得严重 ,因而对较低的Re数 ,计算分离边界层流动采用高阶精度差分格式十分必要。本文同时给出用四阶精度差分格式求解在流动分离情况下边界层方程的技术方法 ;对一个大展弦比后掠翼在攻角等于 1 4°时计算了机翼上表面的分离线 ,并与文献上公布的用别的测算方法测出结果做了比较 ,两者总体上符合良好     

高超声速再入飞行器气动特性的快速预示—局部方法的推广应用

本文将计算高超声速稀薄气流过渡领域中气动特性的局部方法,推广应用到连续介质中弹头型高超声速再入飞行器气动力特性的快速估算。由激波风洞中M_∞=9.9时,一个8°钝锥的气动力测量结果,导出这一实验条件下的领域系数,并以此来估算不同锥角、不同钝度比及不同外形弹头型再入飞行器的气动力和力矩系数,其结果与无粘数值解及实验结果作了比较,在攻角2°～14°范围内吻合得很好。局部方法可用于弹头型高超声速再入飞行器气动特性的快速预示。     

高超声速无粘流场的推进计算

本文将文[1]中用于时间相关法计算的NND格式推广到定常超声速流动的空间推进计算,采用二步的预测、校正方法保证了推进方向的二阶精度,可以证明,这种二阶精度的NND格式具有TVD性质,是MacCormack二步显式格式的推广。本文首先将格式应用于二维平板上斜激波反射流场的推进计算,以检验格式捕捉激波的能力,同时研究了不同的通量分裂方法对格式捕捉激波能力的影响,得到了相当满意的结果。在此基础上,计算了航天飞机简化外形的身部超声速流场,给出了M_∞=10,α=0°,和M_∞=5,α=5°两种状态的部分结果,计算结果清楚地描绘了由于气流在机翼附近受到强烈压缩而产生的内嵌激波与外激波相交的复杂流场结构,与文[7]相比,流场结构更为清晰。     

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。     

用空间推进法来解钝锥三维粘性绕流

本文将文献［１］中提出的方法推广应用于三维可压缩Ｎ－Ｓ方程，文中给出了钝锥超声速三维粘性绕流的结果。得到了背风面分离的流动图画，计算结果与实验相吻合。     

使用RANS/LES混合方法对钝前缘三角翼进行数值模拟(英文)

阐述了建立在一方程湍流模型基础上的分离涡模型(DES)和延迟脱体涡模型(DDES),它们很好地模拟了壁面附近的小尺度流动。该方法在大攻角分离流区域,使用了亚格子应力的Smagorinsky大涡数值模拟,模拟亚音速旋涡流动。使用大攻角下的钝前缘65°大后掠三角翼模型,研究了CFD对涡的发展、破裂以及复杂涡的发展和演变的模拟能力。研究结果说明了方法对模拟分离流动是可行的。     

三角翼大攻角分离流开缝吸气效应研究

分别采用数值方法和实验方法研究了大后掠角三角翼大攻角条件下,背风面开缝吸气对流场结构和气动力特性所造成的影响。数值模拟采用 Ｈａｒｔｅｎ Ｙｅｅ 的二阶精度隐式 Ｔ Ｖ Ｄ 格式和 Ｎ-Ｓ方程;实验采用激光蒸汽屏流场显示和应变天平测力技术。通过在三角翼背风面开缝吸气,抽掉低能气流,实现增升的效果不理想,同二维吸气涡控制增升效果相差很远,这主要是由于二维同三维分离产生机理上的差异所造成的。     

再人体复杂外形无粘流场数值计算研究

本文在文[1,2,3]的基础上研究了非定常Euler方程的推进迭代方法,并将无波动、无自由参数的耗散差分格式(下称NND格式)发展为隐式、迎风格式,用以计算有差曲控制翼的再入复杂飞行器的育攻角、有侧滑角的超声速无粘绕流流场。其结果是准确的。本文的方法具有以下特点: 1.可以在同一计算程序中实现亚跨超流场的计算:在超声速区域,采用推进技术;在亚跨声速区,则自动选择推进迭代求解技术。 2.NND格式不仅具有公式简单,无自由参数的特点,而且过激波时满足熵条件,不会产生非物理解,也不发生振荡。 3.隐式差分格式的无条件稳定和Gauss-Seidel迭代的快速收敛,计算省时。 本文的方法曾用于烧蚀头部凹陷外形和航天飞机简化外形的计算,能自动选择迭代和推进区域。可以预言,如计算机条件许可,本方法可以用于大攻角的无粘流场计算。     

超声速高超声速球锥绕流的边界层稳定性特点初探

本文对一个球锥组合体,在超声速来流及零迎角条件下,对其边界层的稳定性特点进行了分析并讨论了转捩预测问题。结果表明转捩总是最先出现在来流马赫数不太大,或来流雷诺数较高并且具有较大的周向速度或周向波数的情况下。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

钝前缘三角翼无人机气动特性研究

三角翼布局因其优良的气动特性在军用飞机和无人机上获得了广泛应用.为了研究钝前缘三角翼无人机的气动特性,首先采用求解雷诺平均N-S方程的方法对NASA钝前缘三角翼标模进行对比计算,以验证计算方法的可靠度;然后对无人机四个升降舵偏角的气动力和流场特性进行分析研究.结果表明:三角翼无人机在升力系数较小时具有较高的升阻比,当迎角小于1 5°时,钝前缘三角翼前缘气流附体、吸力较高,翼面的横向流动不明显,使飞机的升阻比提高;当迎角大于15°后,涡流特征起主导作用,使得飞机在直到40°迎角范围内没有出现大面积气流分离,具有良好的俯仰稳定性,升降舵效率较高.钝前缘三角翼气动布局在翼展受限、翼载较小的条件下具有一定的气动特性优势.     

抛物化Navier-Stokes方程在求解高超声速流动稳定性分析 的基本流中的应用

通过在边界层内部使用边界层上缘的流向压力梯度代替该处的流向压力梯度,对传统抛物化Navier-Stokes(PNS)方程求解中的流向压力梯度的处理方法进行改进,从而使它可以求解用于稳定性分析的高超声速流动的基本流,并分别计算了平板、零攻角钝锥和小攻角钝锥三种典型算例.对于流向没有分离,且横向流动不强的流动,使用PNS方程计算高超声速流动稳定性分析用的基本流是可靠的.特别是改进后的方法求得的基本流的稳定性分析结果与直接求解Navier-Stokes(N-S)方程的结果非常吻合,但计算的时间消耗和空间占用都减少了一个数量级.     

方形截面飞行器上仰机动对滚转特性影响的数值模拟

飞行器从中小迎角至大迎角范围内,由于背风区流动分离形态的演化,静态气动特性特别是横侧向气动特性也随迎角显著变化,可能诱发复杂的滚转运动。但飞行器一般是上仰机动时,才从平飞状态快速拉起至大迎角,此机动过程对横侧向气动特性和滚转运动可能产生较大影响。本文发展了刚体动力学方程和Navier-Stokes方程的松耦合求解技术,并通过数值模拟航天飞机脱落碎片的六自由度运动轨迹进行了验证。针对背风区涡流形态及横侧向气动特性复杂的方形截面飞行器,数值模拟研究了其不同迎角下的静态滚转气动特性、自由滚转运动特性,以及上仰机动时不同拉起速率对滚转运动特性的影响。结果表明,对于此飞行器,静态时存在临界迎角约为13°,当迎角小于临界迎角时,滚转方向是静不稳定的,诱发快速滚转运动;当迎角大于临界迎角时,滚转方向是静稳定的,其滚转运动是收敛的。但上仰机动时,滚转运动的形态还与拉起速率相关,即使拉起的终止迎角大于临界迎角,如果拉起速率较慢,也可能出现快速滚转运动。