确定大气边界层顶高度的新方法及数值实验

提出了一种确定大气边界层顶高度的数值微分新方法,该方法使用了正则化技术,把对弯角廓线求导数的数值微分问题转化为求目标泛函极小值的问题,采用双参数模型函数方法来选择正则化参数,最后利用最大梯度法确定边界层顶高度.首先通过两个数值实验验证了新方法的有效性,实验结果显示,随着掩星资料噪音的增多,由差分法和结合L曲线方案的数值微分方法得到的边界层顶高度波动增大,而通过双参数模型函数方法得到的高度很稳定,这说明新方法能够很好地过滤噪音,从而保留廓线的主要信息.随后基于2007-2011年1,4,7,10月的COSMIC弯角数据,利用新方法分析了全球海洋大气边界层顶高度的季节特征,并与用掩星资料自带的大气边界层顶高度数据zbalmax得到的季节分布进行对比.结果表明,两者的季节分布特征十分一致:海温相对周围海域高的区域,边界层顶高度较高,反之,边界层顶高度较低;在暖流经过的海域,边界层顶高度较高,在寒流经过的海域,边界层顶的高度相对较低.     

壁面局部吹吸边界层感受性的数值研究

目前理论、实验以及数值模拟主要研究自由来流中的小扰动与壁面局部粗糙相互作用激发边界层感受性问题. 但是, 针对自由来流湍流与壁面局部吹吸相互作用诱导边界层感受性的相关报道甚少. 本文采用直接数值模拟和快速傅里叶变换的方法, 数值研究了二维平板壁面具有局部吹吸的边界层感受性问题. 结果发现, 在二维边界层内能找到一组被激发产生的Tollmien-Schlichting(T-S)波波包的包络序列以及从波包中能够分离出一组稳定的、中性的和不稳定的T-S波, 证明了二维边界层内感受性现象的存在性. 经数值计算获得了T-S波波包传播的群速度; 并建立了自由来流湍流强度、壁面局部吹吸强度和长度与二维边界层感受性之间的关系, 获得了与Dietz感受性实验相类似的结论. 另外, 还发现在自由来流湍流与壁面局部吹、吸相互作用下能诱导二维边界层内产生相位相反的T-S波. 依据这一理论机理来优化设计局部吹吸装置, 不但能促使层流向湍流转捩的提前, 也可以延迟转捩过程的发生, 达到控制湍流运动的目的.     

来流边界层对三角翼诱导涡特性的影响

三角翼是涡发生器的主要结构型式,用于涡发生器的三角翼常处于翼面边界层内。本文采用数值求解雷诺平均NS方程的方法,以平板上一个三角翼模型为研究对象,研究不同来流边界层速度型(均匀分布、湍流边界层速度型和层流边界层速度型)和5种边界层厚度条件下,三角翼诱导涡沿流向的最大涡量分布规律及衰减规律。结果表明,湍流边界层下三角翼诱导涡最大涡量值要大于层流边界层,且边界层高度越低,诱导涡最大涡量值越大。     

转捩边界层中流向条纹的新特性

基于边界层转捩后阶段的高精度直接数值模拟结果,发现流向条纹结构的低速条纹的演化过程中存在不连续现象,以及随高速条纹的发展会出现称之为"高速斑"的新特性. 通过详细剖析边界层转捩过程中的复杂涡系结构以及上喷下扫流动现象,证实流向高低速条纹新特性与流场涡系结构的演化过程紧密相关. 关键词： 流向条纹 边界层 转捩 直接数值模拟     

可压缩流动稳定性方程求解与转捩预测

讨论了Ｒａｙｌｅｉｇｈ逆迭代法与边界层渐近匹配合的方法，并应用于求解高速可压缩边界怪流动稳定性和预报转捩点的数值计算方法。给出了平板边界层第一模式的稳定性分析和转捩区数值预报结果，与已有的实验结果进行了比较。     

平面叶栅边界层的转捩试验及数值模拟

基于温敏漆(TSP)技术,以跨音速平面叶栅为研究对象,开展平面叶栅边界层转捩位置的试验以及数值模拟研究。试验结果表明：利用TSP测量技术可以得到平面叶栅表面的温度分布,进而判断边界层转捩位置。通过CFD(Computational Fluid Dynamics)数值方法完成验证工作。     

用STC格式求解二维激波—边界层相互作用问题

将空间—时间守恒(STC)格式应用于求解N-S方程,并对激波—边界层相互作用问题进行了计算。结果表明,该方法可捕获激波与边界层相互作用的各种现象,显示了优良的数值模拟性能。     

超声速平板圆台突起物绕流实验和数值模拟研究

高速飞行器表面不可避免的存在突起物并形成复杂流场, 从而引起飞行器气动特性和热载荷的变化; 同时, 突起物是流动控制的重要方法之一, 合适的突起物形状及安装位置对于改善冲压发动机进气道性能有重要意义. 本文采用基于纳米粒子的平面激光散射技术(NPLS)研究了马赫3.0来流边界层为层流的平板上三个不同高度圆台突起物绕流流场, 主要关注了突起物后方的尾迹边界层, 并采用高精度的显式五阶精度加权紧致非线性格式(WCNS-E-5)离散求解Navier-Stokes方程模拟了该流场. 获得了超声速圆台绕流精细流场结构, 观察到突起物后方尾迹区域边界层发展的过程. 结合实验和数值模拟结果可以发现, 当圆台高度接近或者小于当地边界层厚度时, 突起物对边界层的扰动非常弱, 圆台后方尾迹边界层能够维持较长距离的层流状态, 在边界层转捩阶段也有清晰的发卡涡结构出现; 反之, 边界层受到的扰动明显增大, 在突起物后方很快发展为湍流; 风洞噪声对本文研究圆台引起的边界层扰动有一定影响, 实验获得的边界层转捩位置要比数值结果靠前. 基于NPLS流场图像, 采用间歇性方法分析了圆台突起物后方边界层的特性, 对于高度大于边界层厚度的圆台其间歇性曲线较为接近并且更加饱满, 边界层的脉动也更为强烈.     

钝锥绕流流动稳定性分析与转捩预报

研究了超音速钝锥绕流的稳定性和转捩点预报的数值计算方法,首先采用Euler方程求解钝锥绕流基本流场,用所得到的物面压力分布作为粘性边界层的外缘压力分布,给出了基本流场的初值;然后应用反迭代法与边界层渐近匹配的方法求解了钝锥边界层的稳定性方程,得到了钝锥边界层转捩数据.该方法可提高计算精度,节约计算时间.     

转捩边界层中的上喷下扫运动

采用高精度直接数值模拟方法,研究平板可压缩流边界层转捩后期过程中的上喷和下扫运动.模拟结果清楚地显示了边界层转捩过程中的上喷下扫运动的复杂现象,展示与典型涡系及尖峰结构的紧密联系.详细剖析多层次的上喷下扫运动的形成机制,数值结果证实了实验的新发现,并进行相应的机理分析.     

基于涡流发生器的风洞试验段附面层控制数值模拟

沿试验段侧壁发展的附面层是影响飞行器半模型实验数据精准度的主要因素之一.利用数值模拟方法验证了涡流发生器减小附面层影响的可行性,重点分析了安装角度、结构尺寸、安装位置及个数等设计参数对附面层内速度分布的影响规律,对涡流发生器尾涡强度以及沿流向的发展规律进行了初步探讨.结果表明,涡流发生器产生的尾涡能够有效改善附面层内的速度分布,进而减小附面层厚度,降低附面层影响;涡流发生器的后缘应略高于当地附面层厚度,安装角度、位置、个数等参数必须合理设计以减小涡流发生器对试验段主气流的影响.基于计算结果初步设计了可用于2.4 m跨声速风洞半模试验段的涡流发生器,在亚声速范围内能够减小模型区侧壁附面层厚度66%左右,对核心流Mach数影响小于0.003,为涡流发生器的实际应用提供了依据.      

叶顶间隙大小和壁面相对运动对低速轴流压气机孤立转子性能的影响

本文针对一低速轴流压气机孤立转子,利用数值方法分析了不同叶顶间隙大小、机匣和叶顶之间相对运动对其总体性能的影响。 结果表明存在着最佳顶部间隙,原因是顶部区域吸力面边界层引起的损失由于叶顶间隙流的作用而减少,虽然存在泄漏损失,但总损失却小于无顶部间隙时的损失。壁面相对静止和相对运动时的顶部区域二次流和损失分布有着明显的区别,因此某些没有考虑壁面相对运动的平面叶栅实验或计算结果在用于压气机转子设计时,需要经过一定的修正。     

The influence of boundary layer on sound propagation in optical fibers

Physical models of two-layered and three-layered fiber were built. These models were used to calculate acoustic properties of optical fibers. Acoustic properties of fibers with boundary layer and without boundary layer calculated from these model were compared. The propagation of acoustic cylindrical symmetric waves in optical fibers was considered. This problem was treated analytically and numerically in continual approximation. The possibility of simultaneous propagation of two cylindrical symmetric waves in optical fibers with a boundary layer was shown (in a three-layered model). The physical phenomenon of the presence of two waves in the fiber with boundary layer is proposed for study of the boundary layer. Formulas for the calculation of the properties of the boundary layer from the acoustic experiment, verified by numerical calculations, are represented.     

Spatial resolution analysis of planar PIV measurements to characterise vortices in turbulent flows

The effects of spatial resolution of planar particle image velocimetry (PIV) on vortex size, swirling strength, circulation and population density characterisation are analysed using a series of experimental and numerical databases. The databases comprise a PIV database of an adverse-pressure-gradient turbulent boundary layer (APG TBL), a PIV database of a zero-pressure-gradient (ZPG) TBL in streamwise-wall-normal planes and streamwise-wall-normal slices of a direct numerical simulation (DNS) of a ZPG TBL. The effects of interrogation window and mesh sizes on the vortex parameters are analysed in the outer region of these flows using different qualitative and quantitative approaches. The quantitative analysis mainly capitalises on the possibility of mimicking the PIV data-sets with the DNS one. These approaches allow us to not only isolate the effects of mesh size and the interrogation window size but also to deduce the combined effects of other measurement errors in PIV. Typical values of mesh size and interrogation window size (0.01–0.03 of the boundary layer thickness) and typical levels of measurement uncertainties have significant effects on the vortex parameters. Moreover, each PIV error source affects the vortex parameters in different and frequently opposite manners. Hence, an optimal selection of measurement parameters such as the interrogation window size is indispensable in order to minimise the effects of spatial resolution and other measurement errors on the vortex parameters. Guidelines are presented in the Conclusions section of this paper. Finally, it is found that all the vortex parameters, when averaged across the outer region, are reasonably comparable in the ZPG and APG TBLs despite the fact that these are very different flows.     

Extensive characterisation of a high Reynolds number decelerating boundary layer using advanced optical metrology

Over the last years, the observation of large-scale structures in turbulent boundary layer flows has stimulated intense experimental and numerical investigations. Nevertheless, partly due to the lack of comprehensive experimental data at sufficiently high Reynolds number, our understanding of turbulence near walls, especially in decelerating situations, is still quite limited. The aim of the present contribution is to combine the equipment and skills of several teams to perform a detailed characterisation of a large-scale turbulent boundary layer under adverse pressure gradient. Extensive particle image velocimetry (PIV) measurements are performed, including a set-up with 16 sCMOS cameras allowing the characterisation of the boundary layer on 3.5 m, stereo PIV and high resolution near wall measurements. In this paper, detailed statistics are presented and discussed, boundary conditions are carefully characterised, making this experiment a challenging test case for numerical simulation.     

Diffraction at an elliptical cylinder with a strongly prolate cross section

The problem of diffraction of a high-frequency plane wave by an infinite cylinder, whose cross section is a strongly prolate ellipse, was considered. The field asymptotics in the boundary layer near the surface were obtained. These asymptotics contain a parameter that characterizes the degree of ellipse elongation, which is equal to the ratio of the transverse wave size squared to the longitudinal size. The approximation accuracy and the applicability domain of the asymptotic formulas were investigated by comparing them to the numerical data obtained in the pdetool environment of the MatLab package.     

Two hardening mechanisms in single crystal thin films studied by discrete dislocation plasticity

Two-dimensional discrete dislocation plasticity simulations of the evolution of thermal stress in single crystal thin films on a rigid substrate are used to study size effects. The relation between the residual stress and the dislocation structure in the films after cooling is analyzed using dislocation dynamics. A boundary layer characterized by a high stress gradient and a high dislocation density is found close to the impenetrable film-substrate interface. There is a material-dependent threshold film thickness above which the dislocation density together with the boundary layer thickness and stress state are independent of film thickness. In such films the stress outside the boundary layer is on average very low, so that the film-thickness-independent boundary layer is responsible for the size effect. A larger size effect is found for films thinner than the threshold thickness. The origin of this size effect stems from nucleation activity being hindered by the geometrical constraint of the small film thickness, so that by decreasing film thickness, the dislocation density decreases while the stress in the film increases. The size dependence is only described by a Hall–Petch type relation for films thicker than the threshold value.     

High-order non-uniform grid schemes for numerical simulation of hypersonic boundary-layer stability and transition

The direct numerical simulation of receptivity, instability and transition of hypersonic boundary layers requires high-order accurate schemes because lower-order schemes do not have an adequate accuracy level to compute the large range of time and length scales in such flow fields. The main limiting factor in the application of high-order schemes to practical boundary-layer flow problems is the numerical instability of high-order boundary closure schemes on the wall. This paper presents a family of high-order non-uniform grid finite difference schemes with stable boundary closures for the direct numerical simulation of hypersonic boundary-layer transition. By using an appropriate grid stretching, and clustering grid points near the boundary, high-order schemes with stable boundary closures can be obtained. The order of the schemes ranges from first-order at the lowest, to the global spectral collocation method at the highest. The accuracy and stability of the new high-order numerical schemes is tested by numerical simulations of the linear wave equation and two-dimensional incompressible flat plate boundary layer flows. The high-order non-uniform-grid schemes (up to the 11th-order) are subsequently applied for the simulation of the receptivity of a hypersonic boundary layer to free stream disturbances over a blunt leading edge. The steady and unsteady results show that the new high-order schemes are stable and are able to produce high accuracy for computations of the nonlinear two-dimensional Navier–Stokes equations for the wall bounded supersonic flow.