平板湍流边界层底层的不稳定波

湍流边界层底层的相干结构及猝发现象对决定边界层的特性有很重要的意义,但直至现在对相干结构的性质及起因仍不清楚。在本文中,对湍流边界层的底层进行了非线性稳定性分析,发现该处确可产生不稳定波,可能用以说明相干结构的起因。     

低湍流度风洞中湍流度对平板边界层转捩影响的试验研究

本文报告了在西北工业大学壁低流度风洞中进行了平板边界层转捩试验研究的简况及初步结果，试验湍流度为０．０２％、０．１％及０．３３％，用恒温热线风速仪测量时均速度型，求得边界层沿流向的位移厚度分布，并用示波器观察速度脉动脉形变化，从而确定起始转捩点和完全转捩点位置。结果表明，转捩的规律性和国外经典结果极为吻合。     

对流边界层湍流特性的数值研究

采用大涡模拟方法研究了存在逆温层的情况下大气对流边界层的湍流特性。实际大气边界层中出现逆温层是较常见的，逆温层会导致大气边界层湍流结构的变化，从而影响大气的湍流扩散和输运特性。本文比较了不同逆温梯度的工况，着重分析了逆温层对边界层中热量逆梯度输运（counter gradient heat transportation，CGHT）的影响。计算结果表明：逆温梯度越大，对流边界层的发展越受到抑制；逆温层高度降低会影响整个对流边界层的温度抬升；逆温梯度越大，垂直速度方差越小；在逆温梯度较大的情况下，其逆梯度输运区域要略微低一些，初步分析认为是由于逆温层对热对流的抑制造成的；对于逆温层高度不同的情况，高度越低的逆温层对逆梯度输运的抑制作用更明显。     

边界层转捩区不稳定流结构的声辐射

基于Ｌａｕｃｈｌｅ的平板边界层转捩区的声辐射理论,通过估算转捩区内各点“上升时间”,同时采用Ｊｏｓｓｅｒａｎｄ的转捩区内间歇因子和湍猝发频率对随机指示器函数的时－空相关分布做修正,从而估算了平板边界层转捩区辐射声功率谱,并对影响其谱峰和谱级的各种因素进行了分析。     

湍流边界层外区相干结构的三维波模型

根据流动稳定性理论,提出了一种三维波模型来描述湍流边界层外区大尺度相干结构。计算所得流线图和等涡量线图较罗纪生,周恒（1993）的二维波模型更符合实验结果。说明该三维模型能够较好地反映湍流边界层外区大尺度相干结构的物理特征。     

竖直平板间自然对流的湍流频谱特性

根据直接数值模拟的计算结果，对竖直平板间湍流自然对流的脉动动能、速度及温度等物理量的时间序列进行频谱分析．结果表明，流动达到充分发展状态后，小尺度到大尺度的能谱很宽，计算的分辨率足够．从能谱分布可以观察到含能区、惯性子区和耗散区的存在，文中对各区的特性进行分析．由于该流动的强各向异性，惯性子区很窄．并讨论了法向位置对脉动动能的影响以及大尺度结构的特性．     

一种适用于超音速边界层的湍流转捩模式

建立一种合理反映扰动模态和可压缩性影响的新型k-ω-γ转捩模式.其主要特点为:(1)假设脉动动能k 由湍流脉动部分和非湍流脉动部分组成,且在模化后者时采用了稳定性分析的结果; (2)在间歇因子γ方程的源项中,构造了具有``自动判断转捩起始位置'功能的函数; (3)通过构造新型的物面法向长度尺度,保证了模式中所有的表达式均由当地变量构成,可以方便地应用于现代CFD程序之中. 该模式在亚音速、超音速和高超音速条件下的边界层流动中进行了验证. 计算结果表明,该模式可应用于较宽马赫数范围内的自然转捩以及旁路转捩过程,所具有的捕捉流动转捩的性能优于国际上的现有模式.      

一种用于分析封闭腔内湍流自然对流换热的k-ε新模型

由于目前用于求解湍流自然对流流动与传热的k-ε模型在应用过程中存在不足,结合高雷诺数k-ε模型需要借助壁面函数法来确定壁面上相关参数值和低雷诺数k-ε模型在近壁区布置更多节点以便获得粘性底层详细信息的特点,重新定义了湍流普朗特数σt的计算式,提出了一种修正的k-ε新模型;利用该模型对封闭方腔内的湍流自然对流流动与传热进行了数值分析。结果表明:与文献中数值模拟结果相比,当108≤Ra≤1014时本文模型所得壁面平均努塞尔特数更接近文献中的实验值,与实验值之间的相对误差在8%以内;壁面的局部努塞尔特数与文献中的实验值吻合得较好。这说明本文模型用于求解封闭腔内湍流流动与传热问题是合适的,比其它湍流模型更能准确地描述封闭腔内湍流自然对流换热中边界层发展与壁面传热特性之间的内在联系。     

湍流边界层中固体颗粒运动的数值模拟

根据Ｌａｇｒａｎｇｅ颗粒运动微分方程及不可压缩湍流边界层中流体的壁面速度分布规律，数值求解了颗粒在湍流边界层中的运动，考虑了Ｓａｆｆｍａｎ升为对颗粒运动的影响，壁面对运动阻力的影响，给出了固体颗粒沉积边壁，在边界层外缘上所需的最小速度和最小入射角，计算结果还表明边界层对固体颗粒撞击边壁的速度和入射角有较大影响，从数值结果可可以发现一个重要现象。     

The unstable modes of natural convection boundary layer

The instability of natural convection boundary layer around a vertical heated flat plate is analyzed theoretically in this paper. The results illustrate that the “loop” in the neutral curve is not a real loop but a twist of the curve in the frequency-wave number-Grashof number space, and there is only one unstable mode at small Prandtl numbers. Specially, when the Prandtl number is large enough two unstable modes will be found in the “loop” region. Along the amplyfying surface intersection the two unstable modes have the same Grashof number, wave number and frequency but different amplifying rates. Their instability characteristics are analyzed and the criterion for determining the existence of the multi-unstable modes is also discussed.     

An experimental research on the instability of natural convection boundary layer around a vertical heated flat plate

Experimental results on the instability of the isothermal naturalconvection boundary layer around a vertical heated flat plate are presented. It is demonstrated that the characteristics of the instability wave in the outer layer is consistent with the calculation of Brewster & Gebhart. After an initial growth of its low frequency components at the downstream side of the turning point of the neutral curve (Gr≈120) its comparatively higher frequency components develop and become turbulent subsequently with a buoyancy subrange in its power spectra. Simultaneously, in the measurement at the inner layer near the wall a viscous instability signal the same as the Tollmien-Schlichting waves in ordinary boundary layer and its subharmonics in a much higher frequency domain is discovered and an inertial subrange can be observed in the spectra atGr≈378.6. The project supported by the National Natural Science Foundation of China (19572004)     

Numerical study of natural convection flow in a vertical slot

The stability of convective motion, generated by a lateral temperature difference across a vertical slot, is studied numerically over a range ofGr=5000 to 1.5 × 10⁵,Pr=0.01 to 10, andA=8,16 and 20. Various cellular flow structures and temperature patterns are illustrated. Several branches of solutions characterized by different numbers of the cells in the flow patterns as well as by both steady and unsteady multicellular patterns are found for low-Prandtl-number fluid in the vertical slot. Meanwhile, the behaviors of the temperature variation and heat transfer are also discussed. The project supported by the National Natural Science Foundation of China (59776011) and by the Returnee from Abroad Funding of Academia Sinica.     

Transient flow rate behaviour in an external natural convection boundary layer

A theoretical approach is proposed to investigate the transient dynamic behaviour of a free convection boundary layer-type flow. The set of continuity, momentum and energy equations are solved with the classical Boussinesq approximation using the Karman–Pohlhausen integral method. Applying a step variation of the uniform heat flux on a vertical wall, the boundary layer thickness and velocity profiles within the viscous layer, streamline patterns and volumetric flow rate are evaluated as a function of time. In addition, corresponding fully analytical asymptotic solutions are derived to be readily used in engineering applications.     

Numerical modeling of laminar-turbulent transition in a boundary layer at a high freestream turbulence level

Disturbances generated by external turbulence in the boundary layer on a flat plate set suddenly in motion are determined by numerically solving the Navier-Stokes equations. The results of direct numerical simulation of isotropic homogenous turbulence are taken as initial conditions. The solution obtained models laminar-turbulent transition in the flat-plate boundary layer at a high freestream turbulence level, time measured from the onset of the motion serving as the longitudinal coordinate. The solution makes it possible to estimate the effect of different factors, such as flow unsteadiness and nonlinearity and the characteristics of the freestream velocity fluctuation spectrum, on laminar-turbulent transition in the boundary layer.     

A direct boundary-layer stability analysis of steady-state cavity convection flow

Natural convection flow in cavities with insulated top and bottom and heated and cooled walls is known to exhibit travelling wave instabilities in the thermal boundary layers that form on the walls. In water (Pr = 7.5) at Rayleigh number Ra = 6 × 10⁸, these waves have been observed at start-up. However no such waves have been observed for the fully developed flow, although it may be assumed that the stability character of the boundary layers is at least approximately the same. The start-up waves are generated by perturbations to the system. In the present paper, an artificial perturbation is applied to the system to determine the stability character of the boundary layers in fully developed flow. It is shown that the thermal boundary layers in the fully developed flow have approximately the same stability character as the start-up flow.     

The effect of Reynolds number on boundary layer turbulence

A new facility for studying high Reynolds number incompressible turbulent boundary layer flows has been constructed. It consists of a moderately sized wind tunnel, completely enclosed by a pressure vessel, which can raise the ambient air pressure in and around the wind tunnel to 8 atmospheres. This results in a Reynolds number range of about 20:1, while maintaining incompressible flow. Results are presented for the zero pressure gradient flat plate boundary layer over a momentum thickness Reynolds number range 1500–15?000. Scaling issues for high Reynolds number non-equilibrium boundary layers are discussed, with data comparing the three-dimensional turbulent boundary layer flow over a swept bump at Reynolds numbers of 3800 and 8600. It is found that successful prediction of these types of flows must include length scales which do not scale on Reynolds number, but are inherent to the geometry of the flow.     

Effect of the freestream turbulence scale and the leading edge shape on boundary layer laminar-turbulent transition

The effect of the leading edge shape and the turbulence scale on laminar-turbulent transition in the flat-plate boundary layer due to grid turbulence is investigated. In the experiments, the turbulence scale was changed by a factor of three and the bluntness radius of the edge by a factor of four, all other factors being fixed. It is shown that on the plate with a sharp edge the fluctuation growth rate and the laminar-turbulent transition point depend nonmonotonically on the turbulence scale. On the blunt plate transition occurs considerably earlier than on the sharp plate.     

Thermal dispersion effects on non-Darcy natural convection over horizontal plate with surface mass flux

Summary The effect of surface mass flux on the non-Darcy natural convection over a horizontal flat plate in a saturated porous medium is studied using similarity solution technique. Forchheimer extension is considered in the flow equations. The suction/injection velocity distribution has been assumed to have power function form Bx ^l , similar to that of the wall temperature distribution Ax ⁿ , where x is the distance from the leading edge. The thermal diffusivity coefficient has been assumed to be the sum of the molecular diffusivity and the dynamic diffusivity due to mechanical dispersion. The dynamic diffusivity is assumed to vary linearly with the velocity component in the x direction, i.e. along the hot wall. For the problem of constant heat flux from the surface (n=1/2), similarity solution is possible when the exponent l takes the value −1/2. Results indicate that the boundary layer thickness decreases whereas the heat transfer rate increases as the mass flux parameter passes from the injection domain to the suction domain. The increase in the thermal dispersion parameter is observed to favor the heat transfer by reducing the boundary layer thickness. The combined effect of thermal dispersion and fluid suction/injection on the heat transfer rate is discussed. Received 7 December 1995; accepted for publication 7 January 1997