Convection of liquid with internal heat release in a rotating container

The convection of heat-generating fluid in a rotating horizontal cylinder is experimentally investigated. The threshold of convection excitation, the structure of convective flows and the heat transfer in the cylinder depending on the heat release capacity, liquid viscosity and aspect ratio of the cavity are studied. It is found that the average convection is excited by the thermovibrational mechanism —the gravity force, rotating in the cavity frame, produces the oscillations of non-isothermal fluid relative to the wall, which in turn result in excitation of mean convective flows. It is shown that the structure of convective flows depends on the dimensionless velocity of rotation. At relatively low rotation velocity the convection develops in the form of a periodic system of vortices regularly distributed along the cylinder axis. The threshold of excitation (critical value of vibration parameter) of three-dimensional vortex structures grows with rotation velocity. Above some definite rotation velocity the convection develops as two-dimensional rolls parallel to the axis of rotation. The threshold of two-dimensional structures excitation does not depend on the rotation velocity. Besides the structure of thermal convective flows the analysis of the relatively weak currents generated by the inertial waves below the threshold of convection is performed.     

失谐叶栅的受迫振动响应特性分析

采用基于计算流体力学(CFD)方法的降阶气动力模型并耦合结构运动方程,实现了存在外激励载荷时失谐叶栅受迫振动响应的快速分析。针对典型的跨声速叶栅,通过求解其位移响应幅值较系统地研究了失谐方式、失谐强度和叶片质量比对失谐叶栅受迫振动响应幅值的影响。研究表明文中刚度失谐形式可以改善叶栅振动的稳定性,同时导致系统受迫振动响应局部化程度的增加,并且受迫响应的最大振幅放大因子随失谐强度增加或者质量比降低存在先增大后减小的一个峰值,不同失谐形式则对这个峰值的大小有着明显的影响。由于该方法可高效地分析失谐叶栅受迫振动各参数对模态局部化的影响,在工程上有一定的应用价值。     

基于IDDES框架的γ-Reθ转捩模型

针对高超声速流动的γ-Re_θ转捩模型在模拟强制转捩时存在捕获边界层内扰动不足的缺点，将RANS（Reynolds Averaged Navier-Stokes）框架改为IDDES（Improved Delayed Detached Eddy Simulation）框架，既能像基于RANS框架的转捩模型一样模拟复杂构型的自然转捩，又能发挥IDDES能够捕获更多脉动信息的优点，较为准确地模拟粗糙颗粒诱导强制转捩。通过对一系列简单构型的自然转捩及来流马赫数为6条件下平板上单个粗糙颗粒诱导强制转捩的模拟表明，模型既能体现γ-Re_θ转捩模型的优点，在自然转捩模拟中具有较强的鲁棒性，能够反映雷诺数等因素对转捩位置的影响规律；也能体现IDDES方法的优点，能够捕捉粗糙颗粒诱导的扰动及涡结构，从而较为准确地刻画出强制转捩的整个流程。     

小振幅振动翼型升力迟滞环变向的成因分析

在低雷诺数条件下研究了SD7037翼型小振幅强迫振动时，振动轴位置和减缩频率对翼型气动特性的影响。通过升力系数随时间变化曲线与迎角随时间变化曲线的对比，发现引起升力系数随迎角变化的迟滞环变向的原因主要是升力系数随时间变化曲线的相位发生变化，并从数学上证明了这个相位是关于减缩频率的函数。而且，当升力系数随时间变化曲线的相位变大时，会使得对应的升力系数随迎角变化的迟滞环曲线从逆时针向顺时针方向变化，这中间必然会经历一个"直线"过程，就好像迟滞现象消失了一样。进一步通过对比流场，发现振动轴位置和减缩频率对此相位的影响机制不同，振动轴位置主要改变翼型的有效迎角，而减缩频率的增大则影响了周围流场结构，使得附加质量带来的反作用力变大，进而提高升力，振动轴向前移动和减缩频率的增大都会增大升力系数曲线的相位。     

多组分液滴自然对流厚交换层蒸发模型及其检验

为发展多组分液体燃料高温蒸发模型,首先以零扩散和无限扩散概念为基础,拓展考虑自然对流的厚交换层高温蒸发模型到多组分液体燃料,提出多组分NC-TEL模型。其次,采用挂滴法对正庚烷-乙醇、正癸烷-乙醇、RP-3航空煤油-乙醇三种混合燃料的单液滴在高温静止和强迫对流条件下的蒸发特性进行实验研究。实验结果显示:混合液滴蒸发速率随温度升高而显著增大,温度越高组分构成比例对液滴蒸发率的影响越明显;本文实验条件下,对流环境对于液滴蒸发的促进作用并不明显。最后,用实验数据检验蒸发模型。模型对比结果显示:总体上,NC-TEL模型优于R-M模型,高温段预测精度平均提升了8%～35%;低温段,零扩散NC-TEL模型与实验结果吻合程度较好,而无限扩散NC-TEL模型与实验结果相比误差略大;高温段,对于正庚烷-乙醇混合燃料液滴,NC-TEL模型预测较为准确,而对于正癸烷/RP-3航空煤油-乙醇混合燃料液滴,NC-TEL模型预测值则偏低,可能的原因是微爆现象和Marangoni现象。     

航天器对流换热地面降压模拟的临界压力准则

航天器密封舱在地面模拟试验中常通过降压法抑制自然对流的影响,但降压比的选择 往往缺乏定量准则,这使得自然对流的抑制效果得不到保证。使用数值模拟方法确定临 界压力比,以几种典型航天器密封舱内的流动换热情况为例,分析了不同Gr／Re 2 数和压力比下密封舱内的流动换热情况。得到了舱内气体温度分布和对流换热系数。通 过比较空间情况和地面情况的计算结果,分析了自然对流给流动换热带来的影响,给出了判 断临界压力的准则式,并给出了临界Gr／Re 2数。结果表明：临界压力 的准则关系和密封舱的形状和内部结构无关。
     

叶片安装角对叶盘结构受迫响应特性的影响

利用简化的叶盘结构有限元模型和基于极值原理的改进Monte Carlo模拟技术,详细分析了叶片安装角失谐对叶盘结构受迫响应特性的影响.分别对6种不同安装角误差的叶盘结构的受迫响应进行了数值计算,得到了各自对应的响应放大因子概率密度曲线,通过曲线比较,分析了叶片安装角误差对叶盘结构受迫响应特性的影响规律.研究表明:安装角误差虽对叶盘结构的固有频率影响较小,但对叶片受迫响应影响是很大的,并且受迫响应幅值放大因子随安装角失谐程度增加而增加.      

环形液池浮力-热毛细对流表面振荡现象的临界温差

对于上部开口的环形液池,加载径向温度梯度将使气液接触面上表面张力分布不均匀,耦合于地面的重力作用,将会驱动薄层流体形成浮力-热毛细对流.当径向温度差△T达到某一临界值△Tc时,一定厚度的液层表面就会出现具有一定规律的振荡现象.实验采用高精度激光位移传感系统,测量了液池表面某点位移随液层厚度h和径向温差△T的变化曲线,讨论了浮力-热毛细对流表面振荡的临界温差及两种驱动力的耦合问题,进行了表面畸变参数的无量纲分析,考虑了起始环境温度对临界温差的影响,推导了振荡相对较强时液层的厚度范围.     

DEVELOPMENT OF HYBRID RANS / LES CODE AND ITS VALIDATION

A computational code based on the hybrid RANSLES approach is developed. The hybrid approach com- bines the delayed detached-eddy simulation（DDES） with an improved RANS-LES hybrid model aiming at wall modeling in LES（WMLES）. In the code, the convective flux is solved using the fourth-order skew-symmetric scheme so as to diminish the negative effect of numerical dissipation. The Spalart-Allmaras（S-A） model is applied as a suhgrid scale（SGS） model. To validate the developed code, homogeneous isotropic turbulence and turbulent channel flow are simulated and the results are compared with experimental data and DNS results. The results of the isotropic turbulence show that the fourth-order skew-symmetric scheme is adequate enough and the model works well coupling with the convective scheme. The results of the turbulent channel flow agree well with the DNS data, the predicted velocity profiles at Reynolds number from 178 to 2 700 match well with the Reichardt＇s law, and the organized vortical structures are well captured.     

基于RANS-LES混合方法的翼型大迎角非定常分离流动研究

使用雷诺平均Navier-Stokes方程-大涡模拟（RANS-LES）混合方法中的延迟分离涡模拟（DDES）方法,模拟了NACA 0015翼型在大迎角下的静态绕流和强迫振荡运动并与实验值进行了比较。在大迎角静态翼型大分离流动模拟中,DDES方法捕获了非定常RANS计算未能获得的机翼背风面的涡脱落现象。在所采用的RANS和DDES模型中,基于剪切应力输运（SST）湍流模型的SST-DDES混合方法给出的时均压力系数分布与实验吻合得最好。在大迎角强迫振荡翼型绕流模拟中,DDES方法得到的非定常气动载荷与实验值吻合得很好,正确地反映了最大迎角处阻力和俯仰力矩的阶跃性突变,而非定常RANS计算则给出了完全错误的趋势。