粗糙物面引起的超声速边界层转捩现象研究

采用风洞实验和数值模拟方法研究平板表面圆柱形粗糙单元引起的M=3超声速边界层转捩问题。实验采用纳米粒子示踪平面激光散射技术（NPLS）对流动结构进行测量。共考察了1mm、2mm和4mm三个不同的粗糙度条件。采用五阶精度加权紧致非线性格式（WCNS-E-5）对风洞实验进行数值模拟和对比。实验及计算表明：粗糙元对边界层干扰后诱导了尾迹流向涡的形成,流向涡通过抬升机制产生剪切层和流向速度条带等不稳定结构;实验流动图像显示,剪切层不稳定在边界层转捩过程中起重要作用;随着粗糙元高度增加,流动分离范围和转捩区域明显扩大,转捩位置有提前的趋势。     

基于升力线理论的机翼几何扭转设计方法

针对飞机设计对于减小诱导阻力的需求,结合升力线理论中三维机翼诱导阻力的相关研究,提出机翼的几何扭转设计方法.此设计方法以总升力、机翼平面形状、翼型气动特性为输入参数,以机翼的几何扭转为输出.在低速情况下,利用此方法得到的机翼具有椭圆形升力分布,从而具有最小的诱导阻力.采用远场法计算诱导阻力,应用参数化方法和CATIA实现机翼的自动生成,应用网格生成工具ICEM的脚本语言实现气动网格自动生成,应用MATLAB的符号运算功能编制程序.该设计方法适用于低速有粘流中平直翼的几何扭转设计,算例的CFD(Computational Fluid Dynamics)计算结果表明:利用该方法设计的机翼展向升力分布为椭圆,并能有效降低机翼的诱导阻力.     

风洞模型静弹性变形对气动力影响研究

介绍了一种基于模型变形视频测量系统和计算空气动力学研究静弹性变形对气动力影响的方法.利用模型变形视频测量系统获取模型在气动载荷作用下的静弹性变形,驱动模型表面网格运动,得到模型变形后的表面CFD计算网格.CFD计算变形前后网格外形下的气动力,研究模型变形对模型气动特性的影响.对一大展弦比连接机翼的测量与计算结果进行了分析,分析结果表明:模型变形对升力系数影响最大发生在升力线性变化的最大迎角附近,模型变形对阻力系数影响最大发生在失速迎角附近,模型静弹性变形对气动力的最大影响量远远超出风洞测力实验的精度指标,因此开展风洞模型静弹性变形影响研究与修正是十分必要的.     

基于5阶精度格式WCNS-E-5的p-multigrid方法研究

p-multigrid方法的基本思想是:在保证收敛结果为高阶精度的同时,利用低阶精度格式耗散大的特点,来改善高精度有限差分格式在迭代计算时收敛速度慢的弱点.本文基于5阶精度WCNS-E-5差分格式,引入1阶精度迎风格式和3阶精度加权格式,构造了p-multigrid方法,在迭代过程中采用了V循环、W循环、S循环、PreV和FMG循环等不同方式来应用这三种格式,并通过典型算例考察了这些循环方式对收敛速度的影响,初步数值试验表明,采用恰当的循环方式,本文所设计的p-multigrid方法能够加快收敛速度,并保证了最终收敛结果与5阶精度WCNS-E-5差分格式的一致性.     

粘性高超声速激波—激波碰撞的计算和分析

本文用特片变量型的隐式ＮＮＤ格式计算了五种类型的高超声速激波－激波干扰流场，得到了清晰的流场结构图画。激波和剪切层等分辨得很清楚。计算的物面压力与实验甚为吻合。借助于激波极曲线，对第４类激波干扰结构进行了分析，解释了干扰引起物面压力增高的原因。     

受DBD等离子体控制的低速流动数值模拟方法研究

基于正离子、电子和氧负离子的三组分空气DBD等离子体动力学模型和流体动力学方程,利用有限体积方法处理控制方程空间项,采用双时间步方法模拟非定常过程,对DBD空气等离子体动力学模型方程和流体动力学方程松耦合求解,数值模拟了高压交流电源激励的低速空气等离子体流动过程。仔细研究了一个周期中真实时间步长对模拟DBD等离子体形成过程的影响,给出了等离子体的形成历程以及等离子体对低速翼型边界层流动的影响,结果在定性上是合理的。     

基于高阶耗散紧致格式的GMRES方法收敛特性研究

计算效率较低是当前限制高阶精度计算方法应用的重要因素。为了提高高阶精度混合型耗散紧致格式（HDCS）的计算效率,发展了适合多块对接网格的广义最小残值（GMRES）方法,并利用GMRES方法开展了HDCS格式的加速收敛研究。首先研究了GMRES的预处理方法、CFL数和内层迭代步数对HDCS数值模拟收敛特性的影响,计算结果显示：点松弛方法是一种高效的预处理方法;CFL数对计算收敛速度影响较大;GMRES方法存在最优的内层迭代步数。利用GMRES方法完成了NACA 0012翼型绕流、NLR 7301翼型绕流和DLR-F4翼身组合体绕流的数值模拟,并与其他隐式时间推进方法进行了对比,GMRES方法计算更加稳定,并且计算效率相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法可以提高5倍以上。研究结果表明,本文发展的GMRES方法在多块对接网格中具有良好的计算稳定性,计算结果的残差可以收敛到更低的量级,并且可以较大幅度地提高高阶精度数值模拟的计算效率。     

三维定常分离流和涡运动的定性分析研究

本文由两部分组成。第一部分研究和分析了物体表面上的分离形态,指出若分离线从奇点始,该奇点为鞍点;若在奇点终,该奇点为结点。若分离线上有很多奇点,其鞍点和结点是交替分布的。分离线的起始,可能有三种形态:一种称之为闭式分离的鞍点起始;一种为正常点起始;一种为鞍、结点组合形态起始。可以把正常点起始和距离很近的鞍、结点组合形态起始称之为开式分离。在一定条件下,闭式分离的形态,先转变为鞍、结点组合的形态,然后过渡到正常点的开式分离。第二部分研究了旋涡沿其轴线发展过程中,其横截面流线形态的发展,指出旋涡轴线上物理量λ=(1/ρ)(ρw/z)是决定流线形态的重要参量。如果λ>0,截面流线在涡心附近是稳定的螺旋点形态。如果λ<0,为不稳定的螺旋点形态,如果λ=0,流线形态是中心型的,当沿涡轴λ由正变负或由负变正时,相应在涡心附近截面流线由稳定螺旋点形态变为不稳定螺旋点形态,或由不稳定螺旋点形态变为稳定螺旋点形态;并且从λ由正到负的变号点起,截面流线产生稳定的极限环,或从λ由负到正的变号点起,产生不稳定的极限环。旋涡破裂只能发生在截面流线为不稳定螺旋点形态的区域中,在通过涡轴的纵向平面内,破裂点处的流线形态是鞍型的。破裂涡的合拢点处在不稳定的极限环区域内。     

强激光与高超声速球锥流场干扰数值模拟研究

本文采用有限差分技术，数值求解带能量源项与空气11组分化学动力学模型的轴对称粘性流体力学方程组，采用类氢原子理论模型计算空气对激光的吸收系数，建立了模拟强激光与超声速钝体流场干扰的计算软件，得到了高能浙江与超声速球锥流场干扰数值模拟的结果。     

高阶精度有限差分方法几何守恒律研究进展

针对高阶精度有限差分格式的几何守恒律问题,系统梳理了国内外离散几何守恒律问题方面的研究工作,以有限差分格式离散后的自由流保持问题为切入点,综述了近年来课题组在有限差分离散几何守恒律方面的研究工作,包括守恒网格导数算法(CMM)、对称守恒网格导数算法(SCMM)等,并通过若干典型算例验证了几何守恒律的满足对高阶精度有限差...