连续双扫掠激波/湍流边界层干扰流动特性研究

为了探究连续双扫掠激波/湍流边界层干扰的流动特性,采用仿真方法对一双尖鳍/平板物理模型进行研究。结果表明:双扫掠激波/湍流边界层干扰形成的两个干扰区存在明显的相干现象,虽然第一道扫掠激波/边界层干扰流动仍具有典型的准锥形相似特性,但受其干扰所形成的非均匀流的影响,第二道扫掠激波/边界层干扰却不再具有准锥形相似特性,同时第二个干扰区将影响其上游临近气流的运动甚至影响第一个干扰区的再附线和分离线等结构。两个干扰区形成各自的λ波结构,并且沿着流向两个干扰区内的激波结构相互汇聚,最终合并为单个更强的λ波结构;不仅如此,两个干扰区内还形成了复杂的旋涡结构,包括一级主旋涡和二级主旋涡,这些旋涡向下游运动,最终融合成一个尺度更大的锥形主旋涡。     

一种基于动模态分解的翼型流动转捩预测新方法

考虑自由转捩的定常/非定常流动Navier-Stokes方程数值求解,对于翼型流动细节的精确模拟和气动力的精确预测均具有十分重要的意义。采用动模态分解(DMD)方法进行流动稳定性分析,再结合e N方法,提出了一套适用于翼型绕流的转捩预测新方法,称为DMD/e N方法。相比于传统的线性稳定性分析方法,DMD方法不需要求解附面层方程和线性稳定性方程,也没有引入平行流假设,具有更好的理论适用性和算法鲁棒性。开展了NLF0416、S809和SD7003等翼型的转捩预测数值验证研究,通过与实验结果以及与传统的基于线性稳定性分析的e N方法的比较,验证了本文所发展的转捩预测新方法在预测翼型的定常流动和非定常流动转捩方面的正确性,也表明了该方法具有解决含层流分离泡的翼型绕流转捩预测的能力。     

NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

自然层流超临界翼型的设计研究

阐述了研究发展我国自然层流超临界翼型ＮＰＵ－Ｌ７２５１３的设计思想,设计要求和转捩位置的判别技术和设计方法,首次提出了利用弱激波形成足够顺压梯度的自然层流超临界翼型设计思想,解决了维持层流所需的有一定顺压梯度压力分布形态和无激波超临界翼型的屋顶状压力分布要求的矛盾。风洞实验结果表明,所设计自然层流超临界翼型达到了设计要求。     

近壁湍流相干结构传播速度的研究

近壁湍流相干结构传播速度的问题,目前尚无理论结果。本文利用流动稳定性理论中的共振三波概念建立单个相干结构模型,采用紧致有限差分和Ｆｏｕｒｉｅｒ谱展开相结合的方法对近壁湍流相干结构进行了数值模拟,计算得到的相干结构流向传播速度与直接数值模拟的结果一致。     

两种湍流模式对绕椭球体三维分离流计算的比较

将新近的Shih Lumley(S L)模式应用于雷诺数为7 2×106的6∶1椭球绕流的全N S方程计算,作为参考也计算了Jonse Launder(J L)模式。在没有发生分离时,两种湍流模式计算的平均流场与实验结果符合较好。在分离发生后,J L模式计算的分离结构与实验差别很大,没有计算出二次分离线;S L模式的计算结果较好,计算得到的表面摩擦力线图谱中主要分离结构与实验结果较接近,分离区的表面压力系数分布也与实验值定性相符;在S L模式计算结果中,壁面附近出现小尺度的分离再附结构.     

二维喷管分离流线性/非线性湍流模式的研究

采用三个线性和七个非线性涡粘性湍流模式，对二维非对称喷管中的不可压缩分离流进行了研究。通过将数值计算结果和实验结果进行比较，对有关的湍流模式进行了评估和分析。计算结果表明，非线性模式的模化系数与平均流动应变不变量以及旋转不变量有关，反映了湍流的各向异性，比线性模式优越得多。另外，所有线性／非线性涡粘性湍流模式不能够捕捉流动的非定常分离过程。     

带有横喷控制的导弹流场数值模拟

从包含多种组分的N-S方程出发，考虑两方程湍流模型，采用NND2M差分格式，对带有横向喷流的双锥旋称体高超声速绕流场进行数值模拟，计算结果与已有的试验数据进行对比，符合较好。在此基础上对带有横向喷流控制系统、型尾翼布局、高超声速飞行的导弹外流场进行数值模拟，研究了迎角、多个喷口、热喷流效应和湍流模型对气动力特性的影响；计算表明在无尾翼情况下有／喷流的气动力差别较小，喷流影响随迎角变化不敏感；对带有尾翼的气动布局，喷口位于背风区时喷流影响较小，喷口位于迎风面时气动力变化较大，压心明显前移；多喷口产生的附加推力和力矩不等于每个单喷口线性相加；湍流模型和热喷流效应引起流场结构改变，但是对总的气动力影响不大。     

浦东地区近地强风特性观测研究

2000年8月当“派比安”和“杰拉华”台风接近上海时，浦东气象局的10m高茺极值风速记录超过24m/s，笔者用超声风速仪在离地20m高处采集到近20h的三维强风样本。经过实测风速数据的分析，得到了平均风速和风向、分析结果表明：近地强风的湍流度和阵风因子较高，湍流积分长度约在80m左右，水平（纵向横向）湍流功率谱密度函数与imiu谱基本 一致，但垂直流端功率谱与Panofsky谱相差较大。     

用预处理方法数值模拟街区峡谷内的污染物扩散

通过耦合求解二维可压缩Navier-Stokes(NS)方程和k-ε湍流模型方程及平流扩散方程，数值模拟了街区谷（Canyon)内的低速湍流流场及大气污染物在其中的扩散。为了加快计算的收敛，对控制方程进行了预处理。预处理后的NS方程的求解采用LU－SGS隐式方法和JST中心差分格式，数据模拟结果表明预处理方法求解低速流场可加快敛，同时就街峡欲内的污染物扩散状况得出了有意义的结论。