Flow characteristics around airfoils near transonic buffet onset conditions

In transonic flow, buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation. The phenomenon is common in transonic flow, and it has serious impact on the structural strength and fatigue life of aircraft. In this paper, three typical airfoils: the supercritical OAT15A, the high-speed symmetrical NACA64A010, and the thin, transonic/supersonic NACA64A204 are selected as the research objects. The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method, and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD). Qualitative and quantitative analysis of the shock wave motion, shock wave intensity, shock foot bubble and trailing edge separation, and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions. The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack, while the growth rate of the primary mode increases inversely. It was found that at the onset of buffet, there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle, and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.     

Some features of the day-to-day MLT wind variability in winter 2017–2018 as seen with a European/Siberian meteor radar network

We present results of wind measurements near the mesopause carried out with meteor radars (MRs) at Collm (51°N, 13°E), Obninsk (55°N, 37°E), Kazan (56°N, 49°E), Angarsk (52°N, 104°E) and Anadyr (65°N, 178°E) from October 1, 2017 till March 31, 2018. The Collm and Kazan MRs are SKiYMET radars with vertical transmission and radio echo height finding, while the other radars operate with horizontal transmission and without height finding. We paid particular attention to the meridional wind variability with periods of 4–6 days and 9–11 days. The waves with these periods are seen as spots of the wave activity in the wavelet spectra and include oscillations with different periods and different discrete zonal wavenumbers. These wave packets successively propagate as a group of waves from one site to another one in such a way that they are observed at one site and almost disappear at the previous one. The 4–6 wave group includes planetary-scale oscillations (individual spectral components) which have eastward phase velocities and mostly zonal wavenumbers 2 and 3, and the vertical wavelength exceeds 70 km at middle latitudes. The source of the oscillations is the polar jet instability. The wave group itself propagates westward, and the amplitudes of wind oscillations are approximately 5–6 m/s as obtained from the wind data averaged over the meteor zone. The 9–11 day wave set propagates westward as a group and mainly consists of spectral components which have westward phase velocity and zonal wavenumber 1. Amplitudes of these wind perturbations strongly vary from station to station and can reach, approximately, 8 m/s. The vertical wavenumber is 0.014 km⁻¹ as taken from the Kazan and 0.05 km⁻¹ according to the Collm data. We obtained a global view on the waves by using the AURA MLS geopotential data. We found a good correspondence between wave features obtained from the MR wind measurements and the MLS data. To our knowledge, such a wave propagation of planetary wave in the mesosphere/lower thermosphere (MLT) region has so far not obtained much attention.     

Gravity wave mixing effects on the OH*-layer

Based on an advanced numerical model for excited hydroxyl (OH*) we simulate the effects of gravity waves (GWs) on the OH*-layer in the upper mesosphere. The OH* model takes into account (1) production by the reaction of atomic hydrogen (H) with ozone (O₃), (2) deactivation by atomic oxygen (O), molecular oxygen (O₂), and molecular nitrogen (N₂), (3) spontaneous emission, and (4) loss due to chemical reaction with O. This OH* model is part of a chemistry-transport model (CTM) which is driven by the high-resolution dynamics from the KMCM (Kühlungsborn Mechanistic general Circulation Model) which simulates mid-frequency GWs and their effects on the mean flow in the MLT explicitly. We find that the maximum number density and the height of the OH*-layer peak are strongly determined by the distribution of atomic oxygen and by the temperature. As a results, there are two ways how GWs influence the OH*-layer: (1) through the instantaneous modulation by O and T on short time scales (a few hours), and (2) through vertical mixing of O (days to weeks). The instantaneous variations of the OH*-layer peak altitude due to GWs amount to 5–10 km. Such variations would introduce significant biases in the GW parameters derived from airglow when assuming a constant pressure level of the emission height. Performing a sensitivity experiment we find that on average, the vertical mixing by GWs moves the OH*-layer down by ~2 to 7 km and increases its number density by more than 50%. This effect is strongest at middle and high latitudes during winter where secondary GWs generated in the stratopause region account for large GW amplitudes.     

航空发动机进气道出口截面喘振波的分布

用扰流插板使涡扇发动机进气道内产生不对称分布的总压畸变并使发动机出现喘振。在进气道出口即发动机进口截面测取两种不同转速下发动机的喘振压力波。应用小波分析奇异点理论,采用模极大值原理和高阶Daubechies小波多尺度分析,对喘振压力信号进行处理,得到发动机喘振波在其进口截面的波形分布特征。揭示了在两种转速下,喘振波在截面上的分布相似,其频率和幅值与转速成正比。截面中心部分表现为正脉冲,周围部分为负脉冲。在垂直方向上的幅值较大,而两翼较小。对该分布特征进行了分析,初步的推论为:高压压气机引起喘振,稳态压力畸变最严重的区域和未受到畸变扰动的区域,喘振波最为强烈。     

流向磁场作用下二维磁流体槽道湍流直接模拟

对低磁雷诺数近似下流向磁场作用的二维磁流体槽道湍流进行了直接数值模拟(DNS,Direct Numerical Simulation),给出了 Re=10 000时不同磁相互作用数下二维磁流体槽道湍流的近壁速度分布、湍流脉动速度均方根、雷诺应力等统计量的变化,并与中性流体二维槽道湍流进行了比较.结果表明:流向磁场作用会导致对数区上移;雷诺应力最大值随磁相互作用数呈线性变化;随着磁相互作用数的增大,下壁面平均涡量的时间演化由拟周期性向周期性转变,且脉动周期逐渐增大并当流动层流化后下壁面平均涡量成为常值.      

对流层雷暴源所产生的柱形重力波

已有充分的证据表明, 大气对流层的雷暴迹象是大气中间层重力波活动的显著代表源.在雷暴迹象的上方, 通过火箭已观测到大气中间层出现的热效应, 也已通过雷 达探测到大气平流层出现的上行重力波, 从地面和卫星平台上观察到夜间气辉有序而 成环状的重力波波形. 所有这些实验结果都与位于观测点下方的雷暴活动有紧密联系. 此类雷暴通常主要集中在中国东部沿海以及地球其他沿海海湾地域. 关于此类雷暴对大 气中间层的影响尚未被充分研究和了解. 为能有效地探究其成因, 利用所开发的一个二维计算机数值模型模拟和研究大气对流层的雷暴源所引发的上行重力波, 进而揭示 此类重力波产生的基本物理机理, 及其在大气中间层的能量耗散. 通过模拟研究发现, 雷暴源可以大面积高强度地聚集和释放积雨云的能量, 当这种周期性的对流变化引发大气对流层的不稳定性后, 就会有圆柱体重力波的产生和传播.      

Heating heavy ions in the polar corona by collisionless shocks: A one-dimensional simulation

Recently a new model for explaining the observations of preferential heating of heavy ions in the polar solar corona was proposed ( and ). In that model the ion energization mechanism is the ion reflection off supercritical quasi-perpendicular collisionless shocks in the corona and the subsequent acceleration by the motional electric field E = −V × B/c. The mechanism of heavy ion reflection is based on ion gyration in the magnetic overshoot of the shock. The acceleration due to the motional electric field is perpendicular to the magnetic field, giving rise to large temperature anisotropy with T_⊥ ? T_∥, in agreement with SoHO observations. Such a model is tested here by means of a one dimensional test particle simulation where ions are launched toward electric and magnetic profiles representing the shock transition. We study the dynamics of O⁵⁺, as representative of coronal heavy ions for Alfvénic Mach numbers of 2–4, as appropriate to solar corona. It is found that O⁵⁺ ions are easily reflected and gain more than mass proportional energy with respect to protons.     

磁流体管流准一维模型的数值验证与性能分析

为了快速分析磁流体管流,采用涡轮喷气发动机模型对磁流体管流进行类比分析,引入变截面管流理论和修改过的多变效率表达式,推导了磁流体发电器和加速器的准一维模型,采用二维数值方法验证了准一维模型的可行性,利用该准一维模型对磁流体管流进行了性能分析,结果表明,多变效率的引入使准一维模型更加完整,当负载系数不变时,磁流体装置中的平均马赫数越低,其等熵效率越高。     

一类TVD型组合差分方法及其在磁流体数值计算中的应用

根据太阳风数值模拟的特点，考虑到算法的质量（收敛速度、稳定性、精度等），结合磁流体数值计算的特性，对三维球坐标磁流体动力学（MHD）方程组中的流体部分采用一种修正Lax-Friedrichs差分法而对磁场部分采用MacComack格式，发展了一类快捷的具有TVD特性的组合数值新方法，作为格式的检验，在一维情况下，将其与PPM格式进行了比较，对一维快慢磁流体激波问题得到了与PPM格式精度相同的结果，然后将其诮到定态太阳风的数值模拟上，在不同等离子体β情形下，可得到理想的太阳风定态结构，为今后将此数值模式应用到具有复杂磁场位型或三维直实太阳风暴的数值模拟研究奠定了基础。     

正丁醇/煤油混合物非预混燃烧压力振荡特性

为探索生物混合燃料在燃气轮机燃烧室内的应用,将丁醇掺混入航空煤油中,根据燃烧室的实验结果构建熵波对流模型,定量分析其不同体积分数对燃烧室内压力振荡频率的影响。实验结果表明:在燃烧室进口压力1.98 MPa、温度600 K和油气比0.03下,随着丁醇所占体积分数增大,燃烧室内压力振荡频率逐渐下降,而幅值变化无明显规律。熵波对流模型结果与实验吻合较好,通过分析表明,不同燃料组分热值变化引起的对流时间改变,是燃烧室压力振荡频率变化的主要原因,为控制燃烧室内压力振荡提供了参考。