二维喷管的初始流动

基于可压缩Navier-Stokes方程,采用大涡模拟方法与高精度混合WENO/TCD格式,对Ma=1.4的超声速平面射流初始流场进行了数值研究.数值结果清晰地描述了超声速平面射流初始流场的结构特征,包括主涡环与激波结构以及它们演变过程.因主涡环内存在涡导激波对,激波与涡相互作用加速射流剪切层失稳,使剪切层首次卷起形成小涡的位置出现在涡导激波对后,此与亚声速射流情况不同.小涡串卷起成后,相继与涡导激波对相互作用,使激波出现明显的变形并加速主涡环失稳.      

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

高马赫数层流摩阻数值计算精度

飞行器在高空高速飞行时黏性效应显著,摩阻预示精度对飞行器的关键气动性能意义重大。目前,摩阻预示主要依赖数值计算,但高马赫数层流摩擦阻力计算与试验测量结果仍存在差距。以具有高速飞行器典型部件特征的球锥、三角翼为对象,结合风洞试验摩阻测量结果,使用国家数值风洞（NNW）数值计算软件和自研CFD程序研究了数值计算中影响摩阻计算精度的格式数值耗散及壁面温度边界条件等重要因素。研究结果表明：数值耗散越小,表面摩阻的计算精度越高;在速度较低的边界层近壁区内关闭熵修正,将有助于提高表面摩阻的预示精度。此外,在高马赫数流动问题的数值模拟中,壁面温度条件对表面摩阻计算同样具有重要影响。最后,基于分析结果和工程需要提出了对高精度摩阻数值预示的研究需求。     

超声速尾迹-剪切流的混合增强

超声速条件下燃料和空气之间的高效混合是超然冲压发动机技术上的主要挑战。基于大涡模拟和流动稳定性分析,针对超声速尾迹-剪切流动开展了混合增强方法研究。尾迹的存在改变了混合层流动的速度剖面,对流动稳定性产生了重要影响,使混合层由三维最不稳定变为二维最不稳定,最不稳定扰动波频率和增长率增大。基于流动稳定性结果引入扰动的混合增强方式依然有效,根据稳定性结果设计了波纹隔板。数值结果表明：二维波纹壁引入的扰动未能增长,不具备混合强化效果,而三维波纹壁引入的扰动能够快速增长,具有混合强化效果,且波纹壁参数越接近最不稳定扰动波参数,混合强化效果越明显。     

Numerical investigation of the impact of asymmetric fuel injection on shock train characteristics

Numerical simulations are carried out to investigate the impact of asymmetric fuel injection on shock train characteristics using the commercial-code FLUENT. The asymmetry of fuel injection is examined by changing the fuel flow rates of the upper and lower wall fuel injectors. The numerical approach solves the two-dimensional Reynolds-averaged Navier–Stokes (RANS) equations, supplemented with a k-ω model of turbulence. As a result, different ways of fuel injections will always lead to shock train transitions, with the variations of shock train structure, strength and leading edge position. For symmetric fuel injection, the flowfield of the isolator is quite asymmetric with the boundary layer of the upper wall side developing much stronger than that of the lower wall, which is due to the heterogeneity of the incoming flow. Regarding to asymmetric fuel injection with more of lower wall side, though the pressures in the combustor are nearly the same, the first shock of the shock train converts between ‘Distinct symmetric X type shock’ and ‘Obscure and weaker asymmetric shock’ and the shock train leading edge moves upstream with the increase of the asymmetry level. With regard to asymmetric fuel injection with more of upper wall side, ‘incomplete asymmetric X type shock’ occurs and the shock train structures keep nearly the same with low level of fuel injection asymmetry. Unexpected results like unstart will happen when increasing the level of fuel injection asymmetry. And the isolator will come back to normal state by decreasing the differential of upper and lower wall sides fuel injections.     

Features of the occurrence of the additional stratification on the bottom-side F region over the equatorial location of Trivandrum

The general features of occurrence of an additional layer on the bottom side of F region, referred to as F0.5 layer in the pre noon period, over the magnetic equatorial location of Trivandrum (8.5° N; 77° E; dip lat of 0.5° N) in India during the period from 2004 to 2007 are presented using ionosonde observations. The F0.5 layer has a June (northern summer) solsticial maximum probability of occurrence with secondary maxima during December (northern winter) solstice. The seasonal as well as the day-to-day variability in the occurrence of F0.5 layer as mentioned in this paper seems to be a result of the variations in the amplitude and phases of the tides and gravity waves, and inventory of the metallic ions of meteoric origin. This study brings out an important manifestation of morning time F layer base region dynamics.     

Stability of the seasonal variations in diurnal and semidiurnal components of mid-latitude F2 layer parameters

We report work utilizing 15-min resolution ionospheric data obtained with DPS-4 digisonde in 2003–2011 to study the seasonal variations in amplitudes and phases of the most powerful spectral components of the F2 layer critical frequency (foF2) and peak height (hmF2) fluctuations over Irkutsk (52.5°N, 104.0°E). We show that fluctuations of both parameters contain quasi-harmonic components with periods of T_n = 24/n h (n = 1–7). The number of distinct spectral peaks varies from 3 in summer to 7 in winter. Amplitude and phase characteristics of the diurnal (n = 1) and semidiurnal (n = 2) components is studied using the data sets extracted from the original data sets with band-pass filter. It has been found that the amplitudes of diurnal/semidiurnal foF2 and diurnal hmF2 components are maximum in winter and minimum in summer. Amplitudes of the diurnal components vary gradually; those of the foF2 semidiurnal one, abruptly, thus forming a narrow winter maximum in November–January. The phase (local time of maximum) of the diurnal foF2 component increases gradually by 4–6 h from winter to summer. The phase of the semidiurnal foF2 component is nearly stable in winter/summer and sharply decreases (increases) by 2–3 h near the spring (autumn) equinox. The phase of the diurnal component of hmF2 (local time of minimum) varies slightly between 1130 and 1300 LT; that of the semidiurnal one decreases (increases) by 4–6 h from January to March (from September to November). The results obtained show that the main features of seasonal variations in the diurnal and semidiurnal components of the mid-latitude F2 layer parameters recur consistently during the solar activity growth and decline phases.     

A comparative sporadic-E layer study between two mid-latitude ionospheric stations

Hourly systematic measurements of the highest frequency reflected by the sporadic-E layer (foEs) recorded from January 1976 to June 2009 at the ionospheric stations of Rome (Italy, 41.8°N, 12.5°E) and Gibilmanna (Italy, 37.9°N, 14.0°E) were considered to carry out a comparative study between the sporadic E layer (Es) over Rome and Gibilmanna. Different statistical analysis were performed taking into account foEs observations near the periods of minimum and maximum solar activity. The results reveal that: (1) independently from the solar activity, Es develops concurrently over extended regions in space, instead of being a spatially limited layer which is transported horizontally by neutral winds over a larger area; especially during summer months, when an Es layer is present at Rome, there is a high probability that an Es layer is also present over Gibilmanna, and vice versa; (2) Es layer lifetimes of 1–5 h were found; in particular, Es layers with lifetimes of 5 h both over Gibilmanna and Rome are observed with highest percentages of occurrence in summer ranging between 80% and 90%, independently from the solar activity; (3) latitudinal effects on Es layer occurrence emerge mostly for low solar activity during winter, equinoctial, and summer months, when Es layers are detected more frequently over Gibilmanna rather than Rome; (4) when the presence of an Es layer over Rome and Gibilmanna is not simultaneous, Es layer appearance both over Rome and Gibilmanna confirms to be a locally confined event, because drifting phenomena from Rome to Gibilmanna or vice versa have not been emphasized.     

超声速复杂流动中湍流模型的性能评估

对8°,16°,20°,24°超声速二维压缩拐角进行了数值模拟,系统评估了SA(Spalart-Allmaras), k-ω, SST(Shear Stress Transport)3种工程常用湍流模型对激波/边界层干扰复杂流动的模拟适用能力.对于"弱"干扰、小分离,工程常用湍流模型能够较准确地预测其压力、摩擦力、热流分布,而对于"强"干扰、大分离则会造成较大的偏差.同时,随着分离区的增大,各湍流模型的适用性能差别更加明显,数值方法也有一定的影响.      

热力学非平衡状态下DSMC绝热壁边界的抽样方法

从直接模拟蒙特卡罗仿真DSMC方法的实现途径出发,经过分析推导,给出了一种在考虑分子转动能、振动能等热力学非平衡状态下统计求解绝热壁温的DSMC抽样方法。然后采用锥体外形的算例验证,分析了该方法在实现绝热壁条件时的统计误差。计算结果表明,此种方法能够有效地实现DSMC方法壁面处的绝热壁条件,随着采样步数的增加,其边界条件实现时的统计误差不断降低,该方法收敛后的统计误差可降至0.5%,能够满足工程要求。