速比对混合层湍流参数影响的实验研究

本文使用PIV术对低高速侧速比为0.25、0.33和0.5时竖直通道内的混合层流动进行实验研究,基于速度差和通道水力直径的雷诺数范围15840～132000.研究发现混合层内湍流参数的分布不仅和雷诺数有关,还和速比有关.混合层内同一横截面上平均雷诺应力的最大值随雷诺数的增大而增大,而在同一横截面上相同雷诺数时雷诺应力的最大值则随速比的增大而减小.在同一横截面上平均涡量随雷诺数的增大而增大,雷诺数相同时平均涡量的最大值随速比的增大而增大.无量纲平均涡量的最大值随混合层的发展按指数规律衰减,速比越大衰减速度越快.     

可压缩混合层流场光学效应分析与实验研究

利用量级分析和风洞实验研究了末制导光学外冷窗口典型流动(可压缩混合层流动)气动光学效应的规律性.理论分析主要针对视线误差(boresight error, BSE)与混合层流场特征参数之间的关系进行了讨论. 研究结果表明: 在可压缩混合层中影响时均BSE的特征参数主要有 对流马赫数、雷诺数、自由流与混合层界面剪切应力、自由流速度比和密度比等因素; 采用细光束穿越混合层流场的风洞试验结果主要证实了时均BSE与对流马赫数之间的关系. 关键词： 气动光学效应 可压缩混合层 对流马赫数     

高雷诺数下非混相Rayleigh-Taylor不稳定性的格子Boltzmann方法模拟

本文采用相场格子Boltzmann方法研究了竖直微通道内中等Atwoods数流体的单模Rayleigh-Taylor不稳定性问题,系统分析了雷诺数对相界面动力学行为以及扰动在各发展阶段演化规律的影响.数值结果表明高雷诺数条件下,不稳定性界面扰动的增长经历了四个不同的发展阶段,包括线性增长阶段、饱和速度阶段、重加速阶段及混沌混合阶段.在线性增长阶段,我们计算获得的气泡与尖钉振幅符合线性稳定性理论,并且线性增长率随着雷诺数的增加而增大.在第二个阶段,我们观察到气泡与尖钉将以恒定的速度增长,获得的尖钉饱和速度略高于Goncharov经典势能模型的解析解[Phys.Rev.Lett.200288134502],这归因于系统中产生了多个尺度的旋涡,而涡之间的相互作用促进了尖钉的增长.随着横向速度和纵向速度的差异扩大,气泡和尖钉界面演化诱导产生的Kelvin–Helmholtz不稳定性逐渐增强,从而流体混合区域出现许多不同层次的涡结构,加速了气泡与尖钉振幅的演化速度,并在演化后期阶段,导致界面发生多层次卷起、剧烈变形、混沌破裂等行为,最终形成了非常复杂的拓扑结构.此外,我们还统计了演化后期气泡与尖钉的无量纲加速度,发现气泡和尖钉的振幅在后期呈现二次增长规律,其增长率系数分别为0.045与0.233.而在低雷诺条件下,重流体在不稳定性后期以尖钉的形式向下运动而轻流体以气泡的形式向上升起.在整个演化过程中,界面变得足够光滑,气泡与尖钉在后期的演化速度接近于常数,未观察到后期的重加速与混沌混合阶段.     

A numerical study of rotation effects on jets effusing from inclined holes into a cross-flow

In this work, the complexity of the flow field arising from the impact of the interaction of coolant jets with a hot cross-flow under rotation conditions was numerically simulated using large eddy simulation with artificial inflow boundary condition. The finite-volume method and the unsteady PISO (Pressure Implicit with Splitting of Operators) algorithm were applied on a non-uniform staggered grid. The simulations were performed for four different values of rotation number (Ro) of 0.0, 0.03021, 0.06042, and 0.12084, a jet Reynolds number of 4700, based on the hole width and the jet exit velocity. The air jet was injected at 30° and 90° in the streamwise direction with a density ratio of 1.04 and a velocity ratio of 0.5. The flow fields of the present study were compared with experimental data in order to validate the reliability of the LES technique. It was shown that the rotation has a strong impact on the jet trajectory behaviour and the film cooling effectiveness. The film trajectory always inclines centrifugally. Under rotating conditions, the film trajectory departs from the centreline to the left boundary. The deflection becomes greater as Ro increases. Furthermore, it was also found that the injection angle has a strong impact on separation and reattachment behaviour as well as the strength of the penetration into the cross-flow. As it increases, the distribution of the film cooling downstream the jet exit is more non-uniform and the film cooling effectiveness level slightly decreases.     

Jet forking driven by pipe tone

The present work deals with an experimental investigation of flow of air through a square-edged circular orifice at the downstream end of a circular duct. Self-excited acoustic oscillations at the natural duct modes are observed for certain flow velocities when the orifice is sufficiently thick. For a specific Reynolds number based on the orifice diameter and the mean jet velocity (9150 < Re < 9850), the jet forks into two trains, with the alternating vortices falling into the same branch of the forked train. Whereas this phenomenon has been reported earlier to have occurred when the density ratio of the jet is less than 0.72, the present results show that it is possible for a jet having the same density as the ambient atmosphere. The jet forking is coincident with jump in the acoustic frequency from one natural acoustic mode to another with comparable amplitudes of both the modes.     

Effects of dimensional wall temperature on velocity-temperature correlations in supersonic turbulent channel flow of thermally perfect gas

Direct numerical simulations of temporally evolving supersonic turbulent channel flows of thermally perfect gas are conducted at Mach number 3.0 and Reynolds number 4800 for various values of the dimensional wall temperature to study the influence of the latter on the velocity-temperature correlations. The results show that in a fully developed turbulent channel flow, as the dimensional wall temperature increases, there is little change in the mean velocity, but the mean temperature decreases. The mean temperature is found to be a quadratic function of the mean velocity, the curvature of which increases with increasing dimensional wall temperature. The concept of "recovery enthalpy" provides a connection between the mean velocity and the mean temperature, and is independent of dimensional wall temperature. The right tails of probability density function of the streamwise velocity fluctuation grows with increasing dimensional wall temperature. The dimensional wall temperature does not have a significant influence on the Reynolds analogy factor or strong Reynolds analogy(SRA). The modifications of SRA by Huang et al. and Zhang et al. provide reasonably good results, which are better than those of the modifications by Cebeci and Smith and by Rubesin.     

Numerical simulations of the flapping of a three-dimensional flexible plate in uniform flow

The self-sustained flapping of a three-dimensional flexible plate in a uniform viscous flow is numerically simulated with a fictitious domain method. The effects of the various non-dimensional control parameters including the Reynolds number, the density ratio, the dimensionless shear modulus, the length–thickness ratio, and the width–length ratio on the flapping of the plate are investigated. The results show that there exist two flapping modes: symmetrical and asymmetrical flapping about the centerline in the spanwise direction. Near the critical point a decrease in the plate width–length ratio, or the increase in the Reynolds number or the reduced velocity (a combination of the density ratio, the dimensionless shear modulus, and the length–thickness ratio) can make symmetric (or nearly symmetric) flapping become asymmetric. It is found that the flapping amplitude is mainly controlled by the density ratio and the dimensionless elastic modulus, while the frequency by the density ratio and the length–thickness ratio. In addition, the flapping amplitude and frequency are affected significantly by the confinement effect of the computational domain, and normally enhanced as the confinement effect becomes stronger. The effects of the plate width and the mass ratio (i.e., the ratio of the length–thickness and density ratios) on the critical reduced velocities are examined. The results indicate that when the fluid–plate mass ratio (or the plate length–thickness ratio) is relatively small there exist two significantly different critical velocities for the flapping instability, depending on the strength of initial plate deformation, a hysteresis phenomenon. No obvious hysteresis can be observed when the fluid–plate mass ratio (or the plate length–thickness ratio) is large.     

驱动模式对具有库源平衡的黏性流体中空间反射时间反演联合对称性的影响

满足空间反射时间反演parity and time-reversal(PT)联合对称性的库源平衡宏观开放系统近几年成为一个研究热点.本文将PT对称性引入到动力学系统,用格子玻尔兹曼方法求解Navier-Stokes方程,发现在二维黏性流体中,如果进口和出口的条件完全等同,在低雷诺数流动中,流场的PT对称函数(ρ)随雷诺数(Re)的增高以ρn～Ren指数增长.用三种不同的速度剖面来驱动流体,计算流场达到稳定状态时的PT对称性.结果发现,进出口平衡的黏性管流中,ρn～Ren的规律在三种驱动模式中出现,表明流场的PT对称性是由流体本身决定的,与驱动模式没有关系,从此论证所得到的指数率的谱适性.     

On the behavior of suspension of drops on an inclined surface

The flow of drops suspended on an inclined surface, are studied by numerical simulations at finite Reynolds numbers. The flow is driven by the acceleration due to gravity, and there is no pressure gradient in the flow direction. The effect of the Reynolds number, the Capillary number and density ratio on the distribution of drops and the fluctuation energy across the channel are investigated. It is found that drops tend to stay away from the channel floor, which is consistent with the behavior observed in the granular flow regime. Drops that are less deformable will stay further away from the channel floor. Also, drops appear at a larger distance from the floor as the Reynolds number increases. Simulations at large density ratios show that results are more compatible with computer simulations of granular flows. The behavior observed here resembles more the granular flow regime when the restitution coefficient is low.     

Wake-boundary layer interaction behind an elliptic cylinder at different Reynolds numbers

In this paper, hot-wire anemometry (HWA) is used to experimentally investigate interactions between a fully developed turbulent boundary layer and wake of an elliptic cylinder where axis ratio (AR) of the cylinder is 2. The elliptic cylinder was located inside and outside a turbulent boundary layer with a thickness (δ) of 0.38B. Furthermore, experiments were conducted at different Reynolds numbers (13,250 and 26,500) based upon the smallest cylinder diameter (B). Mean velocity, turbulence intensity and higher-order central moments of velocity signals (i.e. skewness and flatness) measurements were performed using HWA upon wake-boundary layer interactions on a flat plate. Results showed that profiles of stream-wise mean velocity and turbulence intensity were greatly dependent on gap ratio (G/B) and Reynolds number (Re) in near-wake region. It was also observed that, except for G/B = 0.1, the wake-boundary layer interactions were faster at Reynolds number of 26,500 rather than 13,250. The interactions occurred earlier upon fluctuating the velocity rather than the case where a fixed mean velocity was considered. The results further show that an increase in the gap ratio increases Strouhal number almost independent of δ/B. Behind the cylinder, relatively smaller wake region was obtained at Re = 26,500 rather than Re = 13,250, where the velocity profiles quickly converged to the flat plate boundary layer velocity profiles.     

Numerical Study of Mixed Convective Peristaltic Flow through Vertical Tube with Heat Generation for Moderate Reynolds and Wave Numbers

The current study focuses on the numerical investigation of the mixed convective peristaltic mechanism through a vertical tube for non-zero Reynolds and wave number. In the set of constitutional equations, energy equation contains the term representing heat generation parameter. The problem is formulated by dropping the assumption of lubrication theory that turns the model mathematically into a system of the nonlinear partial differential equations. The results of the long wavelength in a creeping flow are deduced from the present analysis. Thus, the current study explores the neglected features of peristaltic heat flow in the mixed convective model by considering moderate values of Reynolds and wave numbers. The finite element based on Galerkin's weighted residual scheme is applied to solve the governing equations. The computed solution is presented in the form of contours of streamlines and isothermal lines, velocity and temperature profiles for variation of different involved parameters. The investigation shows that the strength of circulation for stream function increases by increasing the wave number and Reynolds number. Symmetric isotherms are reported for small values of time-mean flow. Linear behavior of pressure is noticed by vanishing inertial forces while the increase in pressure is observed by amplifying the Reynolds number.     

Soret效应对具有自由表面的圆柱形浅池内双组分溶液热对流影响的实验研究

为了探究Soret效应对具有自由表面的圆柱形浅液池内双组分溶液热对流过程的影响, 通过实验观察了质量分数为50%的正癸烷/正己烷混合溶液在不同深宽比的液池内流动失稳后的自由表面耗散结构及液池内的温度波动. 结果表明, 双组分溶液流动失稳的临界热毛细Reynolds数小于纯工质的值, 且其随液层深宽比的变化规律与纯工质相同. 当深宽比小于0.0848时, 流动失稳后在自由表面观察到热流体波, 监测点处温度波动主频随热毛细Reynolds数增大而增加; 当深宽比大于0.0848时, 随热毛细Reynolds数的增大, 流动失稳后自由表面依次呈现轮辐状、花苞状、分离-合并-分离交替变化的条纹状结构.     

深径比对邱克拉斯基结构内旋转驱动流动的影响

为了了解深径比对邱克拉斯基(Czochralski)结构内旋转驱动流动的影响,利用有限容积法进行了三维非稳态数值模拟。结果表明:随着液池深径比的增加,流动逐渐加强,当旋转速度超过某一临界值后,流动转变为三维非稳态振荡流动。随着液池深径比的增加,速度波振荡幅度增大,速度波波数和周向传播方向都随之改变;浅液池内坩埚旋转作用占主导地位,速度波传播方向与坩埚旋转方向相同,深液池内晶体旋转大于坩埚旋转对流动的影响,速度波传播方向和晶体旋转方向相同。     

Simultaneous density and velocity measurements in a supersonic turbulent boundary layer

A novel technique for simultaneous measurements of instantaneous whole-field density and velocity fields of supersonic flows has been developed.The density measurement is performed based on the nano-tracer planar laser scattering(NPLS) technique,while the velocity measurement is carried out using particle image velocimetry(PIV).The present experimental technique has been applied to a flat-plate turbulent boundary layer at Mach 3,and the measurement accuracy of the density and velocity are discussed.Based on this new technique,the Reynolds stress distributions were also obtained,demonstrating that this is an effective means for measuring Reynolds stresses under compressible conditions.     

PIV measurements of the flow and turbulent characteristics of a round jet in crossflow

The instantaneous and ensemble averaged flow characteristics of a round jet issuing normally into a crossflow was studied using a flow visualization technique and Particle Image Velocimetry measurements. Experiments were performed at a jet-to-crossflow velocity ratio, 3.3 and two Reynolds numbers, 1,050 and 2,100, based on crossflow velocity and jet diameter. Instantaneous laser tomographic images of the vertical center plane of the crossflow jet show that there exists very different natures in the flow structures of the near field jet due to Reynolds number effect even though the velocity ratio is same. It is found that the shear layer becomes much thicker when the Reynolds number is 2,100 because of the strong entrainment of the inviscid fluid by turbulent interaction between the jet and crossflow. The mean and second order statistics are calculated by ensemble averaging over 1,000 realizations of instantaneous velocity fields. The detail characteristics of mean flow field, streamwise and vertical rms velocity fluctuations, and Reynolds shear stress distributions are presented. The new PIV results are compared with those from previous experimental and LES studies.     

Provide a suitable range to include the thermal creeping effect on slip velocity and temperature jump of an air flow in a nanochannel by lattice Boltzmann method

The thermal creeping effect on slip velocity of air forced convection through a nanochannel is studied for the first time by using a lattice Boltzmann method. The nanochannel side walls are kept hot while the cold inlet air streams along them. The computations are presented for the wide range of Reynolds number, Knudsen number and Eckert number while slip velocity and temperature jump effects are involved. Moreover appropriate validations are performed versus previous works concerned the micro–nanoflows.The achieved results are shown as the velocity and temperature profiles at different cross sections, streamlines and isotherms and also the values of slip velocity and temperature jump along the nanochannel walls. The ability of the lattice Boltzmann method to simulate the thermal creeping effects on hydrodynamic and thermal domains of flow is shown at this study; so that its effects should be involved at lower values of Eckert number and higher values of Reynolds number especially at entrance region where the most temperature gradient exists.     

On the orientational characteristics and transport coefficients of an oblate spheroidal hematite particle in a simple shear flow

We have developed the basic equation of the orientational distribution function of oblate spheroidal hematite particles with rotational Brownian motion in a simple shear flow under an applied magnetic field. An oblate spheroidal hematite particle has an important characteristic in that it is magnetized in a direction normal to the particle axis. Since a dilute dispersion is addressed in the present study, we have taken into account only the friction force (torque) whilst neglecting the hydrodynamic interactions among the particles. This basic equation has been solved numerically in order that we may investigate the dependence of the orientational distribution on the magnetic field strength, shear rate and rotational Brownian motion and the relationship between the orientational distribution and the transport coefficients such as viscosity and diffusion coefficient. We found that if the effect of the magnetic field is more dominant, the particle inclines in such a way that the oblate surface aligns in the magnetic field direction. If the Peclet number increases and the effect of the shear flow becomes more dominant, the particle inclines such that the oblate surface tilts in the shear flow direction. The viscosity due to the magnetic torque is shown to increase as the magnetic field increases, since the magnetic torque due to the applied magnetic field becomes the more dominant effect. Moreover, the viscosity increase is shown to be more significant for a larger aspect ratio or for a more oblate hematite particle. We have applied the analysis to the problem of particle sedimentation under gravity in the presence of a magnetic field applied in the sedimentation direction. The particles are found to sediment with the oblate surface aligning more significantly in the sedimentation direction as the applied magnetic field strength increases.     

载气环境对液态空心玻璃微球运动状态的影响

 研究了载气组份、温度、压力以及微球直径和壁厚对液态空心玻璃微球在炉内运动状态的影响。结果表明：在用干凝胶法制备空心玻璃微球工艺的常用载气组份、温度和压力范围内,载气的组份、温度和压力对相同直径和壁厚的液态玻璃微球在炉内运动速度的影响小于8.3%,但载气组份和压力对液态玻璃微球运动雷诺数和韦伯数的影响很显著。玻璃微球的直径和壁厚是液态玻璃微球运动速度、雷诺数和韦伯数的重要影响因素。提高载气中的氦气含量或降低载气压力可以降低载气对液态玻璃微球的非球形化作用,提高载气温度可以降低其球形化的阻力。     

Density and velocity statistics in variable density turbulent mixing

We experimentally study variable–density mixing of miscible gases in an open-circuit wind tunnel using simultaneous particle image velocimetry and planar laser-induced fluorescence. Experiments of a high Atwood number (0.6) and low Atwood number (0.1) are performed to compare non-Boussinesq cases with the Boussinesq limit. The higher density gas is injected into the wind tunnel co-flow using a round jet configuration, and near-field and far-field measurements are performed to examine mixing in both momentum and buoyancy-dominated regimes. The effects of buoyancy are measurable and important in both large-scale mixing features and in turbulence quantities. The low Atwood number PDFs (probability density functions) show fast and uniform mixing. The high Atwood number PDFs of density have skewness towards the larger densities, indicating less mixing of the heavy fluid due to its inertia. The skewness in the density gradient PDFs at high Atwood number displays strong density local variations that can enhance mixing at molecular scales. Turbulent kinetic energy decreases with streamwise distance from the jet for low Atwood number but increases for high Atwood number due to larger buoyancy and density-driven shear. Over 3000 experimental realisations are used to calculate statistical characteristics of the mixing, including valuable and rarely given data such as Favre-averaged turbulent quantities: mass flux velocity, Reynolds stress, turbulent kinetic energy, and density-specific volume correlation. Buoyancy effects are observed in these quantities and the trends are compared qualitatively with direct numerical simulations.     

Effect of heat generation on transient flow of micropolar fluid in a porous vertical channel

The present work is performed to study the effect of heat generation on fully developed flow and heat transfer of micropolar fluid between two parallel vertical plates. The rigid plates are assumed to exchange heat with an external fluid by convection. The governing equations are solved by using Crank–Nicolson implicit finite difference method. The effects of governing parameters such as transient, heat generation, micropolar parameter, Prandtl number, Biot number, and Reynolds number on the velocity and temperature profiles are discussed. It is found that the presence of heat generation enhances the velocity and temperature of the micropolar fluid at the middle of the channel.