常曲率弯道二维流动稳定性与非线性演化规律分析

国内外对弯曲河流的研究从线性特性转移到不稳定性与非线性特性上来, 弯曲河流作为一个不稳定的系统, 主要是受到水流和边界的不稳定性的影响.针对弯曲河流中存在的水流动力的不稳定性与非线性特征,利用弱非线性理论与摄动法进行展开, 在微弯条件下建立了时间模式下的常曲率 弯道二维水流扰动幅值与扰动幅角的非线性演化方程.研究了在不同弯曲度下的扰动发展特征,探讨了扰动流场在弯曲度影响下的时空分布的规律,分析了扰动流场中出现的扰动漩涡的运动过程, 阐述了弯曲度对弯曲边界内水流不稳定性的影响, 具体表现为在微弯的情况下, 弯道的弯曲度增大会提高河湾内水流的稳定性, 扰动振幅与扰动幅角会随弯曲度的增大其衰减速度更快, 并且扰动流速的对称性在弯曲度较大时减弱, 逐渐向弯道的凸岸偏斜, 在中性状态附近 的弯道水流动力具有对流不稳定性和非线性衰减的特征.本文的研究成果为构建水流动力非线性与床面形态几何非线性相互作用的模型提供了思路.      

超声速平板边界层斜波失稳转捩过程研究

以5阶迎风和6阶对称紧致格式混合差分求解三维可压缩滤波Navier-Stokes方程，对Mach 数为4.5, Reynolds数为10000的空间发展平板边界层湍流进行了大涡模拟. 时间推进采用 紧致存储3阶Runge-Kutta方法，亚格子尺度模型为修正Smagorinsky涡黏性模型. 通过在 入口边界叠加一对线性最不稳定第一模态斜波扰动，数值模拟得到了平板层流边界层失稳转 捩直至湍流的演化过程. 对流场转捩过程中瞬时量及统计平均量的分析表明，数值模拟结果 与理论吻合，得到的Y型剪切层、交替\Lambda涡结构以及转捩后期的发卡涡结构的发展 变化与相关文献结果一致，湍流流谱定性合理.     

流体黏性对输流管道运动方程及临界流速的影响

考虑实际流体黏性引起的管内流速非均匀分布,针对层流和两种不同的湍流流态,对理想流体情况下输流管道运动方程中的离心力项进行了修正,得到的修正系数分别为1.333(圆管层流)、1.020(光滑管壁圆管湍流)和1.037～1.055(粗糙管壁圆管湍流).根据修正后的运动方程得到的上述3种情况下的发散失稳临界流速比理想流体流动情况下依次分别低13.4%,1.0%和1.8%～2.6%.流体黏性对输流管道运动方程及临界流速的影响只与流态有关,雷诺数决定流态,而黏性系数通过雷诺数间接起作用.     

剪切流及分离流结构和层流-湍流转捩的研究(成就和问题)

本文介绍了剪切层(边界层,Poiseujlle.流、分离流、尾流)中层流向湍流转捩过程的连续各阶段的实验研究结果,包括:给定流动对外部扰动的敏感性,产生的波(拟序结构)的线性发展,以及引起有序的层流流态破坏的非线性过程.     

亚临界雷诺数下圆柱绕流的大涡模拟

本文应用Smagorinsky涡粘性模式和二阶精度的有限体积法对圆柱绕流的流场进行大涡模拟．求解了非正交曲线坐标系下的N-S方程,对雷诺数为100和20000的工况进行了计算．计算结果与实验及动力涡粘性模式的结果进行了比较,表明计算对于层流及高亚临界雷诺数的湍流流动是合理的     

壁面定常波纹状吹吸槽道流中湍流特性的研究

在非平衡湍流中,如具有周期性边界条件的流动,由于雷诺应力与平均流速的变形率有着 不同的性质,当周期性边界条件发生变化时,雷诺应力和平均流速变形率的相位对边界条件 的响应也不同,但是二者的相位差在相当大的范围内是稳定的. 这一特性加深了对雷诺应力 的认识,并对非平衡湍流中的模式理论及大涡模拟中亚格子雷诺应力模式的建立提出了许多 需要注意的问题. 利用层流模型,把空间周期性边界条件作为某种扰动,研究了扰动 及其非线性项的分布以及相位间的关系,得到了一些有益的结果.     

有限水深中二维湍流边界层的发展

有限水深中湍流边界层主要是指水流在重力作用下绕建筑物流动时的边界层现象。它区别于一般在无限流场中绕流物体的边界层。它的特点是具有自由表面,边界层有可能发展至全部水深,质量力不容忽视,受到外流流速变化的影响。本实验采用激光流速仪量测二维明渠水流沿程各断面的流速分布,根据实验结果,分析研究了有限水深、粗糙壁面条件下,二维湍流边界层的流速分布特征和厚度发展规律,以及孤立粗糙体对此二者的影响。      

湍流涡致振动频率特性

用数值模拟方法对固定圆柱湍流涡脱落频率与弹性圆柱湍流涡致振动频率特性进行了研究,湍流计算模型采用标准κ-ε模型,压力泊松方程提法基于非交错网格系统.研究结果表明:固定圆柱湍流绕流涡脱落频率基本不随雷诺数而变,对于同一固有频率弹性圆柱,涡振频率基本不随雷诺数而变;对于某一固定雷诺数流动涡振频率在一定范围内与系统固有频率有关.     

等离子体激励器诱导射流的湍流特性研究

为了进一步掌握等离子体流动控制机理, 完善等离子体激励器数学模型, 提升等离子体激励器扰动能力, 采用粒子图像测速技术, 在静止空气下开展了介质阻挡放电等离子体激励器诱导射流特性研究. 实验时, 将非对称布局激励器布置在平板模型上, 随后将带有激励器的模型放置在有机玻璃箱内, 从而避免环境气流对测试结果的影响. 基于激励器诱导流场, 分析了激励电压对诱导射流特性的影响, 揭示了较高电压下诱导射流近壁区的拟序结构, 获得了卷起涡、二次涡等拟序结构的演化发展过程, 计算了卷起涡脱落频率, 阐述了卷起涡与启动涡的区别, 初步探索了卷起涡的耗散机制. 结果表明: (1)层流射流不能完全概括等离子体诱导射流特性, 激励电压是影响射流特性的重要参数. 当电压较低时, 诱导射流为层流射流; 当电压较高时, 诱导射流的雷诺数提高, 射流剪切层不稳定, 层流射流逐渐发展为湍流射流. (2)等离子体诱导湍流射流包含着卷起涡、二次涡等拟序结构; 在固定电压下, 这些涡结构存在恒定的卷起频率. (3)当激励电压较高时, 流动不稳定使得卷起涡发生了拉伸、变形, 引起了流场湍动能增大, 从而加速了卷起涡的耗散. 研究结果为全面认识激励器射流特性, 进一步挖掘激励器卷吸掺混能力, 提升激励器控制能力积累基础.      

曲率对机翼边界层二次失稳影响

采用非线性抛物化扰动方程(NPSE)计算了层流机翼NLF0415(2)在几个工况下横流涡的非线性发展, 应用Floquet理论分析了横流涡的二次失稳. 较系统地分析了后掠机翼在多种参数下表面曲率对其流动稳定性的影响. Haynes证明曲率对于横流在线性稳定性计算(LST)和非线性抛物化扰动方程计算中(NPSE)都起着非常明显的稳定作用; 然而, 该文计算结果表明, 曲率在二次失稳计算中影响不大.      

Turbulence modification in vertical upward annular flow passing through a throat section

Experimental studies on the turbulence modification in annular two-phase flow passing through a throat section were carried out. The turbulence modification in multi-phase flow due to the interactions between two-phases is one of the most interesting scientific issues and has attracted research attention. In this study, the gas-phase turbulence modification in annular flow due to the gas–liquid phase interaction is experimentally investigated. The annular flow passing through a throat section is under the transient state due to the changing cross sectional area of the channel and resultantly the superficial velocities of both phases are changed compared with a fully developed flow in a straight pipe. The measurements for the gas-phase turbulence were precisely performed by using a constant temperature hot-wire anemometer, and made clear the turbulence structure such as velocity profiles, fluctuation velocity profiles. The behavior of the interfacial waves in the liquid film flow such as the ripple or disturbance waves was also observed. The measurements for the liquid film thickness by the electrode needle method were also performed to measure the base film thickness, mean film thickness, maximum film thickness and wave height of the ripple or the disturbance waves.     

Nonlinear evolution of secondary instabilities in boundary-layer transition

Methods of nonlinear stability theory are applied to analyze the evolution of disturbances in the three-dimensional stage immediately preceding the breakdown of a laminar boundary layer. A perturbation scheme is used to solve the nonlinear equations and to develop a dynamical model for the interaction of primary and secondary instabilities. The first step solves for the two-dimensional primary wave in the absence of secondary disturbances. Once this finite-amplitude wave is calculated, it is decomposed into a basic-flow component and an interaction component. The basic-flow component acts as a parametric excitation for the three-dimensional secondary wave, while the interaction component captures the resonance between the secondary and primary waves. Results are presented in two principal forms: amplitude growth curves and velocity profiles. Our results agree with experimental data and the few available results of transition simulations and, moreover, reveal the origin of the observed phenomena. The method described establishes the basis for physical transition criteria in a given disturbance environment.This work has been supported by the Air Force Office of Scientific Research under Contract F46920-87-K-0005 and Grant AFOSR-88-0186 (TH) and by an ONT Postdoctoral Fellowship (JDC).     

On nonlinear water waves under a light wind and Landau type equations near the stability threshold

Various mechanisms of nonlinear saturation of water wave growth under the action of a light wind are discussed. The unstable wind may be saturated by nonlinear dissipation due to the energy transfer to the damping harmonics of the wave. Other nonlinear saturation mechanisms: nonlinear frequency shift, self-modulation or self-focusing of a wave packet may be effective in certain wavenumber regions. In case the wind speed is close to the critical one, an equation is derived for the complex wave amplitude. This equation describes all these nonlinear effects in near-critical systems. In the one-dimensional case this is the nonlinear Shrödinger equation with complex coefficients. Its solutions under various conditions are discussed.     

Two-dimensional analytical solution for compound channel flows with vegetated floodplains

This paper presents a two-dimensional analytical solution for compound channel flows with vegetated floodplains. The depth-integrated N-S equation is used for analyzing the steady uniform flow. The effects of the vegetation are considered as the drag force item. The secondary currents are also taken into account in the governing equations, and the preliminary estimation of the secondary current intensity coefficient K is discussed. The predicted results for the straight channels and the apex cross-section of meandering channels agree well with experimental data, which shows that the analytical model presented here can be applied to predict the flow in compound channels with vegetated floodplains.     

Calculation of three-dimensional weakly expanding flow in jets and channels

The results are given of a calculation of laminar flow in a channel of square section and the motion of a turbulent jet from a cruciform nozzle in an ambient flow. To calculate the secondary flows, the field of the transverse velocity is decomposed into irrotational and solenoidal components. The results of the calculation of the flow in the channel are compared with the calculations of other authors and experimental data. To calculate the flow in the turbulent jet, a one-parameter turbulence model is used, and the influence of the inhomogeneity of the distribution of the longitudinal component of the velocity on the components of the Reynolds stress tensor is taken into account. The results of calculation of the flow in the jet behind a cruciform nozzle are compared with experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 36–44, July–August, 1984.     

THE RENORMALIZATION-GROUP ANALYSIS FOR THE STATISTICAL PROPERTIES OF RAPIDLY ROTATING TURBULENCE

利用重正化群方法对强旋转湍流场统计性质予以研究, 通过重正化微扰展开, 对高波数速度分量进行逐 阶平均.计算结果显示当旋转角速度Ω → ∞时, 用以表征高波数速度分量对低波数速度分量影响的重正化黏性将趋于0, 这表明在强旋转条件下科氏力将抑制湍流速度分量之间的非线性相互作用, 从而阻碍湍流的能量级串效应, 当Ω → ∞时湍流的能量级串效应消失, 导致湍流脉动消失, 流动将层流化.理论计算结果还显示对于强旋转湍流, 时域-空域联立Fourier的湍流速度分量存在二维化趋势, 球面平均能谱函数有标度关系E(k) ∝ k^-3.     

Development of three-dimensional disturbances in a boundary layer

In the region of transition from a two-dimensional laminar boundary layer to a turbulent one, three-dimensional flow occurs [1–3]. It has been proposed that this flow is formed as the result of nonlinear interaction of two-dimensional and three-dimensional disturbances predicted by linear hydrodynamic stability theory. Using many simplifications, [4, 5] performed a calculation of this interaction for a free boundary layer and a boundary layer on a wall with a very coarse approximation of the velocity profile. The results showed some argreement with experiment. On the other hand, it is known that disturbances of the Tollmin—Schlichting wave type can be observed at sufficiently high amplitude. This present study will use the method of successive linearization to calculate the primary two- and three-dimensional disturbances, and also the average secondary flow occurring because of nonlinear interaction of the primary disturbances. The method of calculation used is close to that of [4, 5], the disturbance parameters being calculated on the basis of a Blazius velocity profile. A detailed comparison of results with experimental data [1] is made. It developed that at large disturbance amplitude the amplitude growth rate differs from that of linear theory, while the spatial distribution of disturbances agree s well with the distribution given by the natural functions and their nonlinear interaction. In calculating the secondary flow an experimental correction was made to the amplitude growth rate.     

摆动河槽水动力稳定性特征分析

河流形态与水动力结构息息相关,形态约束水动力结构,水动力结构则通过泥沙运动进一步塑造形态,在自然界河流中形成一对辩证互馈关系.天然河流形态形式多样,大致可以分为顺直、微弯、分叉和散乱游荡几种类型,其中微弯及多个弯曲构成的河型为河流动力演化中最重要的一环.多个弯曲构成的河型可用正弦派生曲线来描述,它也是天然河流主槽与水动力结构复杂相互作用的结果.作为探讨这一过程的力学作用机理,构建摆动槽道并研究槽道摆动与其内部流动结构的互馈关系,既是流体力学研究的热点内容,也是目前河流动力过程研究的基础内容.在此重点讨论这一互馈关系前一部分,即水流对摆动边界的响应.文中建立了随体坐标系下摆动河槽与内部水流动力响应理论模型,通过给定摆动弯曲槽道的不同特征参数,研究讨论了正弦派生型摆动边界下的槽道水流动力稳定性特征,明确了弯曲槽道摆动对其内部主流及扰动水流结构的影响,确定弯曲槽道摆动波数、摆动频率对扰动流发展影响的相应参数定量关系,得到了槽道弯曲度和摆动特征对其内部水流不同尺度扰动影响的阈值选择性范围.     

The characteristics of spiral waves in an axially rotating pipe

An experimental study of the instability of a flow in an axially rotating pipe is performed by means of LDV and flow visualization technique. It is found that the axial velocity of the rotating pipe flow fluctuates like a sine wave at first, then its fluctuating pattern assumes a somewhat sawtooth wave form as a spiral wave appears, which is predicted by means of linear and nonlinear stability analysis. At a certain rotation rate, the amplitude of the velocity fluctuations amounts to 30% of the axial velocity. At the down-stream section, another fluctuating component appears in the velocity, which interferes with the initially appearing component, then the fluctuation becomes one with broad-band spectral components. There is a close analogy between this spectral evolution and that of a Taylor-Couette flow. Deformation of the velocity distribution is obtained from the velocity fluctuating pattern and its phase, and the structure of the spiral wave is considered. The strength, azimuthal wavenumber and angular velocity of the spiral wave obtained from the velocity data are confirmed by flow visualization. The change of pressure loss in the rotating pipe is compared with the case without rotation.