Control Of Vortex Shedding From A Square Cylinder

Small circular, square, and thin-strip cross-sectional elements are used to suppress vortex shedding from a square cylinder at Reynolds numbers in the range of 1.12 x 10(4)-1.02 x 10(5). The axes of the element and cylinder are parallel. The element's size, position, and angle of attack are varied. Measurements of the fluctuating surface pressures and wake velocities, together with smoke flow visualization, show that vortex shedding from both sides of the cylinder is suppressed and the mean drag and fluctuating lift on the cylinder is reduced if the element is installed in an effective zone downstream of the cylinder. The effective zone of the circular element is shown to be much smaller than those of the other elements. The effects of Reynolds number and blockage ratio are investigated. A phenomenon of monoside vortex shedding is observed. The role of the element's bluffness is investigated and the suppression mechanism is discussed.

A Note On The Numerical Simulations Of Flow Past A Wavy Square-Section Cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Karman vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.

MILU-CG method and the numerical study on the flow around a rotating circular cylinder

A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

Streamline topology and dilute particle dynamics in a Karman vortex street flow

Three types of streamline topology in a Karman vortex street flow are shown under the variation of spatial parameters. For the motion of dilute particles in the Karman vortex street flow, there exist a route of bifurcation to a chaotic orbit and more attractors in a bifurcation diagram for the proportion of particle density to fluid density. Along with the increase of spatial parameters in the flow field, the bifurcation process is suspended, as well as more and more attractors emerge. In the motion of dilute particles, a drag term and gravity term dominate and result in the bifurcation phenomenon.

Numerical estimate of the stability curve for the flow past a rotating cylinder

It is demonstrated that the primary instability of the wake of a two-dimensional circular cylinder rotating with constant angular velocity can be qualitatively well described by the Landau equation. The coefficients of the Landau equation are determined by means of numerical simulations for the Navier-Stokes equations. The critical Reynolds numbers, which depend on the angular velocity of the cylinder, are evaluated correctly by linear regression. (C) 2004 American Institute of Physics.

Variation of added mass and its application to the calculation of amplitude response for a circular cylinder

In the present study, analyzed are the variation of added mass for a circular cylinder in the lock-in ( synchronization) range of vortex-induced vibration (VIV) and the relationship between added mass and natural frequency. A theoretical minimum value of the added mass coefficient for a circular cylinder at lock-in is given. Developed are semi-empirical formulas for the added mass of a circular cylinder at lock-in as a function of flow speed and mass ratio. A comparison between experiments and numerical simulations shows that the semi-empirical formulas describing the variation of the added mass for a circular cylinder at lock-in are better than the ideal added mass. In addition, computation models such as the wake oscillator model using the present formulas can predict the amplitude response of a circular cylinder at lock-in more accurately than those using the ideal added mass.

Compact finite difference - Fourier spectral method for three-dimensional incompressible Navier-Stokes equations

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.

Numerical study of bluff body flow structures

Our recent progress in numerical studies of bluff body flow structures and a new method for the numerical analysis of near wake flow field for high Reynolds number flow are introduced. The paper consists of three parts. In part one, the evolution of wake vortex structure and variation of forces on a flat plate in harmonic oscillatory flows and in in-line steady-harmonic combined flows are presented by an improved discrete vortex method, as the Keulegan-Carpenter number (KC) varies from 2 to 40 and ratios of U-m to U-0 are of O(10(-1)), O(10) and O(10), respectively. In part 2, a domain decomposition hybrid method, combining the finite-difference and vortex methods for numerical simulation of unsteady viscous separated flow around a bluff body, is introduced. By the new method, some high resolution numerical visualization on near wake evolution behind a circular cylinder at Re = 10(2), 10(3) and 3 x 10(3) are shown. In part 3, the mechanism and the dynamic process for the three-dimensional evolution of the Karman vortex and vortex filaments in braid regions as well as the early features of turbulent structure in the wake behind a circular cylinder are presented numerically by the vortex dynamics method.

The wave model of mesothermal plasma near wakes and Korteweg-de Vries equation

The stationary two-dimensional (x, z) near wakes behind a flat-based projectile which moves at a constant mesothermal speed (V∞) along a z-axis in a rarefied, fully ionized, plasma is studied using the wave model previously proposed by one of the authors (VCL). One-fluid theory is used to depict the free expansion of ambient plasma into the vacuum produced behind a fast-moving projectile. This nonstationary, one-dimensional (x, t) flow which is approximated by the K-dV equation can be transformed, through substitution, t=z/V∞, into a stationary two-dimensional (x, z) near wake flow seen by an observer moving with the body velocity (V∞). The initial value problem of the K-dV equation in (x, t) variables is solved by a specially devised numerical method. Comparisons of the present numerical solution for the asymptotically small and large times with available analytical solutions are made and found in satisfactory agreements.

钝体尾流控制机理及方法研究进展

首先从涡脱落生成理论出发对钝体尾流控制方法进行了分类,并简单介绍r国内尾流控制研究情况.之后介绍r我们用窄条或小方柱取代小圆柱后,对Strykowsky和Sreenivasan控制方法的改进及其在高雷诺数下对圆柱和方柱尾流涡脱落的有效抑制情况,并探讨了控制件钝度对抑制效果的影响.第3部分用实验数据对各个涡脱落生成模型做了分析与检验,指出控制件方法的机理与改变钝体分离位置、减小钝体背压吸力、改变流动的腱向相关性、防止钝体两侧剪切层相互作用等无关,而与钝体近尾流速度剖面的局部修正及其稳定性的改变有关.最后简单介绍了控制件方法今后研究工作展望及其工程应用前景.

湍流等离子体尾迹雷达散射截面的计算及其影响因素分析

针对一阶畸变波Born近似模型,深入分析了湍流内外尺度和电子数密度脉动值对雷达散射截面的影响,并且给出了关于湍流外尺度的一个经验的,能被工程上较好使用的公式.在以上分析的基础上,计算了几种高程条件下再入小钝头锥体等离子尾迹的雷达散射截面,与已有实验结果进行了对比分析.分析和计算结果表明,湍流的外尺度和局部电子数密度值对雷达散射截面值影响较大;湍流内尺度变化的影响不大.

再入飞行器湍流尾迹流场研究

再入飞行器湍流尾迹流场状况，直接系统到飞行器的雷达散射特性。对再入飞行器湍流尾迹等离子体场理论模型，试图通过湍流模式理论来表达，即使用k-ε-g模型方程来封闭平均化的全Navier-Stokes方程，从而准确获得流动平均场和脉动场住处。使用的N－S平均议程由质量加权平均过程产生，湍流模型方程也经过可压缩性能修正。真实气体效应重点考察空气处于局部热化学平衡状态。流动控制方程运用一个二阶TVD格式有限体积法求解。以一典型小钝锥体零攻角再入飞行为例，计算了在两高程（H = 40 km和H = 30 km）条件下的高超声速湍流尾迹流场。获得的尾迹流场参数与流动物理状况符合，并且湍流脉动参数与已有相应的实验结果定性一致，初步证实该方法合理。

再入尾迹湍流对雷达散射截面影响分析

通过探讨高超声速再入湍流尾迹等离子体场中电磁波的散射机制，推导出在工程上描述湍流亚密等离子体雷达散射的一阶畸变波Born近似模型，分析了该模型在充分发展湍流尾迹等离子体场中的适用性，完成了适用于三维尾迹等离子体场的程序设计。以已有的湍流尾迹等离子体流场数据为基础，分析了再入尾迹湍流等离子体流动对雷达散射截面的影响。选取考察的几个有代表性的因素为：湍流模型、转捩过程、湍流尺度、电子组分脉动初始条件等。由结果可以看到，湍流转捩过程和湍流尺度对雷达散射截面值影响不大；电子组份脉动强度初始值影响较明显；湍流模型在特定条件下影响亦不大。

Experimental Study on Flow Induced Vibration of a Cylinder with Two Degrees of Freedom Near a Rigid Wall

The flow-induced vibration of a cylinder with two degrees of freedom near a rigid wall under the action of steady flow is investigated experimentally. The vibration amplitude and frequency of the cylinder and the vortex shedding frequency at the wake flow region of the cylinder are measured. The influence of gap-to-diameter ratio upon the amplitude response is analyzed. The experimental results indicate that when the reduced velocity (Vr) is in the range of 1.2 < Vr < 2.6, only streamwise vibration with small amplitude occurs, whose frequency is quite close to its natural frequency in the still water. When the reduced velocity Vr > 3.4, both the streamwise and transverse vibrations of the cylinder occur. In this range, the amplitudes of transverse vibration are much larger than those of streamwise vibrations, and the amplitudes of the streamwise vibration also get larger than those at the range of 1.2 < Vr < 2.6. At the range of Vr > 3.4, the frequency of streamwise vibration undergoes a jump at certain values of Vr, at which the streamwise vibrating frequency is twice as much as the transverse one. However, when the streamwise vibration does not experience a jump, its frequency is the same as that of the transverse vibration. The maximum values of second streamwise and transverse vibration amplitudes increase with increasing gap-to-diameter ratios.

再入尾迹电磁特性的湍流效应研究

本文旨意在于通过探讨高超声速再入尾迹中的湍流等离子体与电磁波相互作用的机理，以及建立能反映此机理的应用性理论模型，从而提供一套可进行目标特性分析的方法，以便为工程部门的突防技术服务。本题目在再入气动物理现象研究中具有重要意义。综合分析指出，地面雷达观测到的非相干散射信号主要来源于再入尾迹的亚密湍流区产生的体积散射。因此，电磁散射特性分析主要针对尾迹亚密湍流等离子体。并且，这里所有的分析都是根据在工程应用中最成熟的一阶畸变波Born近似理论模型。再入尾迹电磁特性的湍流效应研究，着眼点就在于湍流等离子体场的研究。对湍流等离子体场理论模型，本文试图通过模式理论来表达，即求解平均化的全Navier-Stokes方程及其封闭方程k-ε-g模型，从而准确获得流动平均场和脉动场信息。这种表达方式较以前有了较大改进。注意到高超声速流动具有强烈可压缩性的特点，故使用的N-S平均方程由质量加权平均过程产生，湍流模型方程也经过可压缩性修正。方程的离散求解方法，都是运用带矢通量分裂的二阶TVD格式的有限体积法。再入尾迹湍流场的初始条件由近尾迹（底部）流动经N-S方程求解给定，初始值更加准确可靠。尾迹从层流到湍流的转捩过程采用相对成熟的半经验公式确定。飞行器的高超声速再入过程必然导致它周围的空气温度升高，使得流动表现出真实气体效应。对重点考察的湍流流动而言，真实气体效应主要表现为气体处于热化学平衡状态。就工程部门面临的实际问题，把一阶畸变波Born近似的解算方法做些改进，使其能够处理的范围从轴对称尾迹扩展到三维湍流等离子体场是必要的。这为深入的理论分析提供了有力的保障。在能够准确模拟湍流流动的刻划雷达散射截面的基础上，考察亚密湍流等离子体对电磁散射的影响。通过选择的几个有代表性的因素进行讨论，初步结果表明：湍流转捩方式、湍流尺度对尾迹雷达散射截面值计算影响不大，而电子组份脉动能初始值影响较明显，且在特定条件下湍流模型的影响亦不大。但由于湍流模型涉及脉动初始值，其影响需进一步确定。同时，一些今后开展继续此项研究工作的有益建议也提了出来。