纵向受迫振荡圆柱绕流问题的数值模拟

用有限体积法对平行于均匀来流方向受迫振荡的圆柱绕流问题进行了二维数值模拟.雷诺数选取Re=200、855、4000等几种亚临界雷诺数情况.通过研究不同振幅和振动频率下的流场结构和一些重要流动参数如升阻力系数、Strouhal数等随Re数、KC数、Stokes数的变化关系,验证了实验中观察到的一定条件下发生的"频率锁定"现象,并将涡脱落方式划分为三种主要模态.文中引入网格速度,对常用的处理速度与压力耦合的SIMPLE算法作了适当的补充和修改,以适应随时间变化的网格坐标.     

Experimental investigation on laminar separation control for flow over a two-dimensional bump

The laminar boundary layer separation flow over a two-dimensional bump controlled by synthetic jets is experimentally investigated in a water channel with hydrogen-bubble visualisation and particle image velocimetry (PIV) techniques. The two-dimensional synthetic jet is applied near the separation point. Two Reynolds numbers (Re = 700 and 1120) based on the bump height and free-stream velocity are adopted in this experiment, and seven different excitation frequencies at each Reynolds number are considered, focusing on the separation control as well as the vortex dynamics. The experimental results show that the optimal control can only be achieved within some excitation frequencies at both Reynolds numbers. However, beyond this range, further increasing the excitation frequency leads to an increase in the separation region. The proper orthogonal decomposition (POD) technique and vortex identification by swirling strength (Λ_ci) are applied for the deeper analysis of the separated flow. The reconstructed Λ_ci field by the first four POD modes is used and vortex lock-on phenomenon is observed. It is found that the negative synthetic jet vortex with clockwise rotation draws the separated wake shear layer as it is convected downstream, and then they syncretise together. Thus, the new vortex is induced and shedding downstream periodically.     

Radon变换在渗流力学反问题中的应用

研究了在渗透系数相差不大的渗流场中椭圆型偏微分方程的系数反问题,通过把CT技术中的Radon变换推广到渗流力学中,给出了渗透系数的计算方法和实例。     

Experimental characterization of onset of acoustic instability in a nonpremixed half-dump combustor

This paper reports work on a nonpremixed half-dump combustor, in which methane is injected at the backward-facing step, and mixes and burns with the air flowing past the step in the unsteady recirculation zone. The flow and geometric parameters are widely varied, to gradually change from conditions of low-amplitude noise to excitation of high-amplitude discrete tones. The purpose of the work is to focus on the transition from the former condition to the latter, and to mark the onset of instability. Dimensionless groups such as the Helmholtz and Strouhal numbers are formed based on the observed dominant frequencies, whose variation with the air flow Reynolds number is used to identify the oscillations as those due to the natural acoustic modes or the vortex shedding process. High-speed chemiluminescence imaging reveals shedding of vortical structures in the flame zone. With variation in the conditions, flow-acoustic lock-on and transition from one vortex shedding mode to another is marked by nonlinearity in the corresponding amplitude variations. Such conditions are identified as the onset of instability in terms of the ratio of the flow time scale to the acoustic time scale and mapped against the operating fuel-air equivalence ratio of the combustor.     

Flame-acoustic coupling of combustion instability in a non-premixed backward-facing step combustor: the role of acoustic-Reynolds stress

Combustion instability in a laboratory scale backward-facing step combustor is numerically investigated by carrying out an acoustically coupled incompressible large eddy simulation of turbulent reacting flow for various Reynolds numbers with fuel injection at the step. The problem is mathematically formulated as a decomposition of the full compressible Navier–Stokes equations using multi-scale analysis by recognising the small length scale and large time scale of the flow field relative to a longitudinal mode acoustic field for low mean Mach numbers. The equations are decomposed into those for an incompressible flow with temperature-dependent density to zeroth order and linearised Euler equations for acoustics as a first order compressibility correction. Explicit coupling terms between the two equation sets are identified to be the flow dilatation as a source of acoustic energy and the acoustic Reynolds stress (ARS) as a source of flow momentum. The numerical simulations are able to capture the experimentally observed flow–acoustic lock-on that signifies the onset of combustion instability, marked by a shift in the dominant frequency from an acoustic to a hydrodynamic mode and accompanied by a nonlinear variation of pressure amplitude. Attention is devoted to flow conditions at two Reynolds numbers before and after lock-on to show that, after lock-on, the ARS causes large-scale vortical rollup resulting in the evolution of a compact flame. As compared to acoustically uncoupled simulations at these Reynolds numbers that show an elongated flame with no significant roll up and disturbance in the upstream flow field, the ARS is seen to alter the shear layer dynamics by affecting the flow field upstream of the step as well, when acoustically coupled.     

圆管管翅散热器肋侧主流方向绝对涡通量研究

本文以圆管管翅式散热器为研究对象,采用数值模拟的方法对圆管管翅式散热器肋侧空气通道主流方向绝对涡通量的特性进行了探讨。给出了肋侧空气通道主流方向绝对涡通量无量纲化参数,二次流雷诺数;横向管间距、肋侧空气流动雷诺数对二次流雷诺数及努塞尔数的影响。     

On the influence of the Reynolds number on the intensity of vortex sound flowing around a cylindrical profile

Regularities of the emission of vortex sound (eolian tone) during air flow around stationary and rotating cylindrical profiles have been investigated. The influence of the flow Reynolds number on the intensity of vortex sound emission has been estimated from results of measuring the pressure fluctuation distribution on the surface of stationary cylindrical rods flowed around by air, as well as in the wake behind them. It is shown that the emission intensity depends on the location of the point of flow detachment from the profile surface and the track width near the profile. The ranges of the flow Reynolds numbers where the emission intensity increases with different flow velocities have been determined by analyzing the dependence of the profile lift coefficient on the Reynolds number. An independent way of determining the profile lift coefficient by measuring the vortex sound intensity is proposed. The results explain contradictions between the results of some authors, who experimentally observed different dependences of emission intensity on the flow velocity. The influence of the profile diameter on the vortex sound emission intensity has been investigated. The boundary Reynolds number above which the profile diameter does not affect emission has been established for stationary and rotating cylindrical profiles. It is shown that deposition of rough coatings on the rod surface may reduce the vortex sound emission intensity by affecting the point of flow detachment from the surface.     

The laminar horseshoe vortex upstream of a short-cylinder confined in a channel formed by a pair of parallel plates

A flow visualization experiment was performed in order to characterize the laminar horseshoe vortex system that appears upstream of the junction of a short cylinder and a pair of flat parallel plates. The experiments were performed in a water tunnel and the technique used for flow visualization was laser illumination of seeded particles whose traces were captured using long exposure photography. Geometrical and flow parameters, such as Reynolds number and height-to-diameter ratio of the cylinders, are varied during the experiments and the flow regimes are analyzed as a function of these parameters. The behavior of vortex systems is reported. For low Reynolds number cases, the vortices stay in a fixed position, as the Reynolds number is increased the number of vortices grows and for larger Reynolds numbers the vortex system becomes oscillatory and for further increases it becomes periodic. As for the dimensionless height of the cylinders, the vortex system is weak for short cylinders and increases its strength and number of vortices as the cylinder height-to-diameter ratio is increased. For further increases in height the vortex system do not change, which shows that the flow becomes independent of the height-to-diameter ratio for sufficiently tall cylinders. Information of the frequency of appearance of periodic vortices is also included.     

压缩拐角流动中G?rtler涡特性及热流分布实验研究

对高超声速压缩拐角流动中G?rtler涡特性及热流分布进行了实验研究.开发了温敏漆（temperature sensitive paint,TSP）系统,简要介绍了TSP技术的原理、文章所用的TSP涂料的标定曲线、辅助设备参数、实验过程数据后处理过程,采用基于离散Fourier定律的热流算法.研究在Ma=6低噪声风洞中进行,采用TSP技术,得到压缩拐角斜坡板上的热流分布图像,并对高低热流条带现象做出解释,与G?rtler涡有对应关系.通过改变拐角角度及来流参数,获得了不同拐角和单位Reynolds数条件下的热流分布图像,分析得到压缩拐角斜坡上G?rtler涡特性及热流分布在变参数条件下的变化规律.研究发现:当增加拐角角度或增大单位Reynolds数时,G?rtler涡的波长减小,且涡的起始位置更靠近拐角;随单位Reynolds数增加,斜坡上热流值整体增加,热流峰值位置前移;峰值位置后,热流缓慢减小的区域与G?rtler涡位置相对应.      

Local Heating Effect of Flow Past a Circular Cylinder

Heating effects of air flows past a two-dimensional circular cylinder at low Reynolds numbers and low Mach numbers are investigated by numerical simulation. The cylinder wall is heated partially rather than heated on the whole surface as with previous researches. The heating effects are completely different for various heating locations on the cylinder surface. Heating either windward or leeward side stabilizes the flow and reduces or completely suppresses vortex shedding from the cylinder at supercritical Reynolds numbers, which is consistent with previous results of heating on the whole surface of the cylinder. However, as the lateral sides of the cylinder （perpendicular to the stream-wise direction） are heated, an adverse effect is found for the first time in that the flow is destabilized and vortex shedding can be excited at subcritical Reynolds numbers. As the lateral sides of the cylinder are cooled, the flow is stabilized.     

Experimental investigation of a swirling flow in a cubic container

A flow induced by a rotating disk mounted at the top of a cubic container totally filled with a liquid is studied experimentally. The flow pattern is visualized for Reynolds numbers in the interval 1500–6000, and the velocity variation along the axis of the container, which coincides with the axis of revolution of the disk, is observed by means of Doppler laser anemometry (DLA). As Re grows, the velocity axial component starts fluctuating because of the vortex core precession. The breakdown of the vortex helical structure becomes pronounced at Re>4000 without the formation of the return flow region (vortex breakdown bubble) at the axis. With the Reynolds number and the container height-to-disk radius ratio being the same, the axial flow patterns in standard cylindrical [1] and cubic containers differ radically. In the latter vessel at low Re, the steady flow regime and axisymmetric bubble breakdown of the vortex structure near the axis are absent.     

Control of a coupled map lattice model for vortex shedding in the wake of a cylinder

The flow behind a vibrating flexible cable at low Reynolds numbers can exhibit complex wake structures such as lace-like patterns, vortex dislocations and frequency cells. These structures have been observed in experiments and numerical simulations, and are predicted by a previously developed low-order coupled map lattice (CML). The discrete (in time and space) CML models consist of a series of diffusively coupled circle map oscillators along the cable span. Motivated by a desire to modify the complex wake patterns behind flexible vibrating cables we have studied the addition of control terms into the highly efficient CML models and explored the resulting dynamics. Proportional, adaptive proportional and discontinuous non-linear (DNL) control methods were used to derive the control laws. The first method employed occasional proportional feedback. The adaptive method used spatio-temporal feedback control. The DNL method used a discontinuous feedback linearization procedure, and the controller was designed for the resulting linearized system using eigenvalue assignment. These techniques were applied to a modeled vortex dislocation structure in the wake of a vibrating cable in uniform freestream flow. Parallel shedding patterns were achieved for a range of forcing frequency-forcing amplitude combinations studied to validate the control theory. The adaptive proportional and DNL methods were found to be more effective than the proportional control method due to the incorporation of a spatially varying feedback gain across the cylinder span. The DNL method was found to be the most efficient controller of the low-order CML model. The required control level across the cable span was correlated to the 1/1 lock-on behavior of the temporal circle map.     

三维曲面腔顶盖驱动流的MRT-LBM研究

采用多松弛时间格子玻尔兹曼方法(MRT-LBM)的D3Q15模型分别对长方体腔、圆柱腔、半圆柱腔、旋转双曲面腔、旋转椭球面腔、半球腔以及两种组合腔体的三维顶盖驱动腔流进行数值模拟, 比较分析各腔体内流线分布、流速等值线分布和涡心的发展, 对于典型腔体模拟不同雷诺数下的流动情况。结果表明: 在同一雷诺数下, 曲面边界不仅能消除从边界产生的次涡, 还会导致腔内主涡的分离, 增大中心纵剖面纵向回流速度; “上长方体+下半圆柱”腔内流函数分布与边界贴合度最高。当雷诺数不断增大时, 半圆柱腔内主涡逐渐分离成两个同向涡, “上圆柱+下半球”腔内始终保持着圆柱腔与半球腔内的基本流动特征; 而长方体腔内主涡涡心保持在同一高度, 次涡逐渐增强, “上长方体+下半圆柱”腔内流动愈加规则, 主涡逐渐下沉, 流速等值线分布逐渐趋于中心小、四周大。     

基于间断有限元方法的并列圆柱层流流动特性

基于高阶的间断有限元方法, 数值模拟低马赫数下并列圆柱的可压缩层流流动, 捕捉并列圆柱流场中的漩涡结构, 以便分析并列圆柱尾流的流动特性. 针对二维圆柱的边界形式, 采用曲边三角形单元构造二维圆柱的曲面边界, 以适应高阶离散格式的精度. 在验证方法合理性的基础上, 分析圆柱间距及雷诺数对漩涡脱落及受力特性的影响规律. 研究结果表明: 并列圆柱的间距是影响流场流动特性的一个主要因素, 它会改变圆柱漩涡脱落的形式. 随着圆柱间距的增加, 上下圆柱的平均阻力系数及平均升力系数的绝对值随之显著下降. 雷诺数对于平均阻力系数的影响相对较小. 但随着雷诺数的增加, 上下圆柱的平均升力系数会随之降低, 而漩涡的脱落频率会随之增大.     

气固两相混合层流场双向耦合的数值研究

采用双向耦合模型对有涡配对的二维气固两相混合展数值模拟,在考虑颗粒对流场反作用基础上进一步对颗粒间通过流体的相互作用进行分析。流场用拟谱方法求解,颗粒用颗粒轨道模型跟踪。结果发现,流场中大涡卷起和配对仍居主导地位;颗粒Ｓｔ数为Ｏ（０．１）～Ｏ（１）时,颗粒减弱了流场雷诺应力强度,加快涡量扩散; Ｓｔ数为Ｏ（１）时,颗粒分布极不均匀,主要集中在渴的边缘．     

Numerical analysis of the influence of angle of attack on turbulent flow around a thick goettingen airfoil with vortex cells

On the basis of the solution by multi-block computational technologies of Reynolds equations closed with the aid of the equations of the model of Menter’s shear stresses transfer, an analysis of the flow around a thick airfoil of classical geometry with vortex cells is given at an arrangement of suction from the surface of central bodies placed inside them. The suction velocities, angles of attack, and location of vortex cells on the contour are determined, under which the flow around an airfoil of 35,2% thickness is ensured close to a separation-free flow, for high Reynolds numbers (Re = 10⁵). The integral force characteristics of the Goettingen and EKIP profiles are compared for the distributed and concentrated suction in vortex cells. The work was financially supported by the Russian Foundation for Basic Research (Grants Nos. 06-08-81002 and 05-01-00162) and by the European Union, the program Framework-6 (Project VortexCell2050).     

A study on drag reduction of a rotationally oscillating circular cylinder at low Reynolds number

The flow field around a rotationally oscillating circular cylinder in a uniform flow is studied by using a particle image velocimetry to understand the mechanism of drag reduction and the corresponding suppression of vortex shedding in the cylinder wake at low Reynolds number. Experiments are conducted on the flow around the circular cylinder under rotational oscillation at forcing Strouhal number 1, rotational amplitude 2 and Reynolds number 2,000. It is found from the flow measurement by PIV that the width of the wake is narrowed and the velocity fluctuations are reduced by the rotational oscillation of the cylinder, which results in the drag reduction rate of 30%. The mechanism of drag reduction is studied by phase-averaged PIV measurement, which indicates the formation of periodic small-scale vortices from both sides of the cylinder. It is found from the cross-correlation measurement between the velocity fluctuations that the large-scale structure of vortex shedding is almost removed in the cylinder wake, when the small-scale vortices are generated at the unstable frequency of shear layer by the influence of rotational oscillation.     

A numerical analysis on flow field and heat transfer by interaction between a pair of vortices in rectangular channel flow

This paper is a numerical study of the effect of flow field and heat transfer created by interactions between a pair of vortices generated by a vortex generator in a rectangular channel flow. In order to analyze the vortices produced by the vortex generator, the pseudo-compressibility method is introduced into the Navier–Strokes (NS) equation of a three-dimensional unsteady, incompressible viscous flow. A two-layer k–ε turbulence model is used on the flat plate three-dimensional turbulence boundary to predict the turbulence characteristics of the vortices. The computational results accurately predict the vortex characteristics, which are related to Reynolds stress, turbulent kinetic energy, and flow field. Also, in the prediction of thermal boundary layers, skin friction characteristics, and heat transfers, the present results are reasonably close to the experimental results obtained by other researchers.     

Vortex dynamics of a trapezoidal bluff body placed inside a circular pipe

The influence of Reynolds number and blockage ratio on the vortex dynamics of a trapezoidal bluff body placed inside a circular pipe is studied experimentally and numerically. Low aspect ratio, high blockage ratio, curved end conditions (junction of pipe and bluff body), axisymmetric upstream flow with shear and turbulence are some of the intrinsic features of this class of bluff body flows which have been scarcely addressed in the literature. A large range (200:200,000) of Reynolds number (Re_D) is covered in this study, encompassing all the three pipe flow regimes (laminar, transition and turbulent). Four different flow regimes are defined based on the distinct features of Strouhal number (St)–Re_D relation: steady, laminar irregular, transition and turbulent. The wake in the steady regime is stationary with no oscillations in the shear layer. The laminar regime is termed as irregular owing to irregular vortex shedding. The vortex shedding in this regime is observed to be symmetric. The emergence of separation bubble downstream of the bluff body on either side is another interesting feature of this regime, which is further observed to be symmetric. Two pairs of mean streamwise vortices are noticed in the near-wake regime, which are termed as reverse dipole-type wake topology. Beyond the irregular laminar regime, the Strouhal number falls gradually and vortex shedding becomes more periodic. This regime is named transition and occurs close to the Reynolds number at which transition to turbulence takes place in a fully developed pipe. The turbulent regime is characterised by a nearly constant Strouhal number. Typical Karman-type vortex shedding is noticed in this regime. The convection velocity, wake width formation length and irrecoverable pressure loss are quantified to highlight the influence of blockage ratio. These results will be useful to develop basic understanding of vortex dynamics of confined bluff body flow for several practical applications.     

基于浸入式边界方法的串联双矩形柱绕流数值模拟

基于浸入式边界算法(Virtual Boundary Method)中力源反馈边界的思想,改进其原有内部流体处理方法以减少计算耗费,并结合非等间距网格以便工程应用计算,模拟雷诺数范围内(Re=200–10³)串联双矩形柱绕流,研究表明:Re=200–300时,前柱尾流涡脱处于双剪切层控制阶段;柱间涡街为Karman类涡街,在小间距条件下被抑制,形成涡环;前柱对后柱屏蔽效应体现为后柱阻力系数远小于前柱;临界间距时柱间涡街充分发展,后柱阻力系数等气动参数亦在此发生跃升,但仍小于前柱值;随雷诺数升高,尾流涡街尺寸缩小,临界间距及跃升幅度变小. Re=400时,前柱尾流涡脱进入冲击剪切层控制阶段,阻力系数不再呈现规律性振荡;此后随雷诺数升高,冲击剪切层逐步完善,前柱流动分离使其表面产生更多附着涡,导致尾流旋涡尺寸进一步减小,屏蔽效应消失,涡脱更为剧烈,进而对后柱产生脉动冲击效应;适当间距比条件下此类脉动冲击效应使得后柱阻力系数发生跃升,并略高于前柱. 关键词： 浸入式边界算法 串联双矩形柱 屏蔽效应 临界间距