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.     

Numerical Simulation of Two-Dimensional Flow over Three Cylinders by LatticeBoltzmann Method

The numerical simulation using the multiple relaxation time lattice Boltzmann method (MRT-LBM) is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these three circular cylinders, one of the three cylinders on which a forced in-line vibrating is used to do this research and attempt to find out the effects of the movingcylinder and the other two rigid cylinders on the wake characteristics and vortex formation. As a benchmark problem to discuss the problem of lift coefficient r.m.s for these cylinders with spacing ratios T/D between other rigid side-by-side cylinders, and the calculation is carried out with two compared cases at Reynolds number of 100, two of the cylinders are rigid and the other one is an in-line vibrated cylinder lying downstream, in addition, forced vibrating amplitude and frequency are A/D= 0.5 and f_ν=0.4 (where A is the forced amplitude, D is the cylinder diameter, and f_ν stands for the vibrating frequency, respectively). The calculated results not onlyindicate that the spacing ratios T/D (T is the center-to-center spacing between the two upstream cylinders) have influence on the wake patterns and the formation of vortex shedding, but also analyze the lift coefficient r.m.s for the three cylinders with the spacing ratios S/D (where S is the center-to-center spacing between the center of upstream two side-by-side cylinders and downstream cylinder).     

Lattice Boltzmann simulation of flow around two,three and four circular cylinders in close proximity

Cross-flows around two, three and four circular cylinders in tandem, side-by-side, isosceles triangle and square arrangements are simulated using the incompressible lattice Boltzmann method with a second-order accurate curved boundary condition at Reynolds number 200 and the cylinder center-to-center transverse or/and longitudinal spacing 1.5D, where D is the identical circular cylinder diameter. The wake patterns, pressure and force distributions on the cylinders and mechanism of flow dynamics are investigated and compared among the four cases. The results also show that flows around the three or four cylinders significantly differ from those of the two cylinders in the tandem and side-by-side arrangements although there are some common features among the four cases due to their similarity of structures, which are interesting, complex and useful for practical applications. This study provides a useful database to validate the simplicity, accuracy and robustness of the Lattice Boltzmann method.

     

On the application of the two-level large-eddy simulation method to turbulent free-shear and wake flows

We present application of the hybrid two-level large-eddy simulation (TLS-LES) method, a multi-scale simulation model, to turbulent free-shear and wake flows at moderately high Reynolds number. The TLS-LES method combines the scale-separation-based two-level simulation (TLS) model with the spatial-filtering-based conventional large-eddy simulation (LES) model in an additive manner using a normalised blending function. The additive blending can be performed in a static or a dynamic manner. We demonstrate that the method, which has been originally developed for wall-bounded flows, can be used to simulate flows in complex configurations without requiring any further adjustments to the model. In this study, three canonical flows are simulated, which are representative of free-shear and wake flows. These cases include a temporally evolving mixing layer, flow past a circular cylinder in a uniform flow and flow past a finite-span airfoil placed in a uniform flow at three different angle of attacks. We analyse the role of static and dynamic blending functions, large-scale grid resolution and the effect of small scales on the instantaneous flow features and turbulence statistics. The results obtained from these cases demonstrate robustness, accuracy and consistency of the multi-scale TLS-LES method and show that the method is suitable for investigation of turbulent flows that encompass features such as massive separation, reattachment, transition to turbulence and unsteady wake, which are challenging to model numerically.     

有变高度圆柱干扰的平板表面流谱拓扑结构

采用分区方法和Roe三阶流通量差分分裂格式数值求解雷诺平均N-S方程,湍流附加粘性系数由修正的Baldwin-Lomax模型计算,数值模拟了高超声速绕变高度圆柱流动,进而根据拓扑学理论给出了有变高度圆柱干扰的平板表面流谱拓扑结构,并结合对称截面流谱进行了简要的分析与探讨,指出了需要进一步深入研究的问题。     

Transition of free-surface flow modes in taylor-couette system

This paper presents our numerical and experimental results of the bifurcation found in Taylor-Couette system with a free surface. The lengths of the two concentric cylinders are finite and their axes are parallel to the direction of the gravitational force. When the end walls of the cylinders are fixed and stationary, numerical and experimental studies have shown that the flow has multiple patterns depending on the cylinder lengths and the Reynolds numbers. Experimental studies on flows with free surfaces also gave various flow modes. Our result shows that the measured and predicted time-dependent displacements of the free surface are in favorable agreement. In case of the cylinder length comparable with the gap width between the cylinders, gradual accelerations of the inner cylinder bring the normal mode flows with one, three and five toroidal vortices. The exchanges of stabilities between these flow modes are summarized in a phase diagram.     

An immersed boundary method based on discrete stream function formulation for two- and three-dimensional incompressible flows

An immersed boundary method is proposed in the framework of discrete stream function formulation for incompressible flows. In order to impose the non-slip boundary condition, the forcing term is determined implicitly by solving a linear system. The number of unknowns of the linear system is the same as that of the Lagrangian points representing the body surface. Thus the extra cost in force calculation is negligible if compared with that in the basic flow solver. In order to handle three-dimensional flows at moderate Reynolds numbers, a parallelized flow solver based on the present method is developed using the domain decomposition strategy. To verify the accuracy of the immersed-boundary method proposed in this work, flow problems of different complexity (decaying vortices, flows over stationary and oscillating cylinders and a stationary sphere, and flow over low-aspect-ratio flat-plate) are simulated and the results are in good agreement with the experimental or computational data in previously published literatures.     

颗粒增强湍流的新模型和气粒流动的湍流变动

用雷诺应力方程模型和极细的网格系对单个颗粒受湍流气体绕流进行了数值模拟,研究了改变颗粒直径和气体相对速度时颗粒增强气体湍流的规律.据此构造了颗粒尾涡增强气体湍流的新模型.将此子模型加入到两相流动模型中,对竖直和水平通道内气粒两相流动进行了数值模拟,和实验结果的对照表明,考虑颗粒尾涡增强气体湍流效应得到的气体湍流脉动速度的模拟结果比不考虑此效应的模拟结果好得多.     

Acoustic testing of the vortex structure produced by an air flow about an array of cylinders

Results of experiments on the scattering of a plane ultrasonic wave from a vortex wake formed in an air flow behind a lattice of vertical cylinders are presented. The lattice is periodic in the direction perpendicular to the oncoming flow. The experiments are performed in a wind tunnel for two values of the Reynolds number, namely, Re = 75 and 500, and for lattices with different numbers of cylinders and with different lattice periods g = (2.5–15)d (where d is the diameter of the cylinders). The measured parameters of the scattered waves are used to estimate the degree of transverse correlation between the vortex wakes formed behind the cylinders for flows with different Reynolds numbers. The results obtained from an analysis of the characteristics of the scattered sound are compared with the results of direct hot-wire anemometer measurements and with the data obtained by other researchers.     

A non-linear fluid force model for vortex-induced vibration of an elastic cylinder

Even under the assumption of a sinusoidal lift and drag force at a single frequency for a stationary cylinder in a cross flow, higher harmonics that represent non-linearity in the fluid-structure interaction process are present. This fact is considered in the formulation of a non-linear fluid force model for a freely vibrating cylinder in a cross flow. The force model is developed based on an iterative process and the modal analysis approach. The fluid force components in the model can be evaluated from measured vibration data with the help of the auto-regressive moving averaging (ARMA) technique. An example is used to illustrate that non-linear (higher order) force components are present at resonance, even for a case with relatively weak fluid-structure interaction. Further analysis reveals that the fluid force components are dependent on structural damping and mass ratio. The non-linear fluid force model is further modified by taking these considerations into account and is used to predict the dynamic characteristics of a freely vibrating cylinder over a range of Reynolds numbers, mass and structural damping ratios. On comparison with measurements obtained from four different experiments and predictions made by previous single-degree-of-freedom model, good agreement is found over a wide range of these parameters.     

旋流燃烧室内同向和反向旋转射流湍流流动的数值模拟

针对发展高效低污染旋流燃烧技术的需要,对旋流燃烧室内两股同轴旋转射流相互作用的湍流旋流流动进行了数值模拟.计算中采用了一种新的代数Reynolds应力模型和QUICK离散格式.在两股射流同向旋转和反向旋转两种条件下,将模拟得到的燃烧室内湍流旋流流动的时均气体轴向速度场、切向速度场和静压场与实验数据进行了比较.     

Numerical simulation and flow visualization using soap film of the self-organized vortex structure in the wake of an array of cylinders

Abstract  

With the aid of computational fluid dynamics (CFD) and simple flow visualization technique using flowing soap-film, we present here the wake structures behind an array of cylinders for Reynolds numbers corresponding to both laminar and turbulent flow regimes. The image results illustrate interesting vortex interactions past these equally spaced cylinders; for low Reynolds number flow, well-organized wake pattern persists and manifests unsteadily to different symmetry states. An increase of free stream flow velocity causes the wake transition, resulting in the formation of asymmetric flow wake with chaotic mixing at the far wake. Observations from both the numerical simulations and soap-film are in good agreement at least qualitatively.     

A mixed-time-scale low-Reynolds-number one-equation turbulence model

In this paper, a new low-Reynolds-number (LRN) one-equation turbulence model for eddy viscosity is proposed. A mixed time scale, representing a combination of three time scales: two time scales made of strain-rate parameter S and vorticity parameter Ω and the turbulent time scale k/?, is introduced into this model. The proposed model is derived from an LRN k?? two-equation model where the mixed time scale has been proved to be very effective for predicting local flows over complex terrains. In the transport equation of the model, the mixed time scale is included in the production and the dissipation terms. The new model is evaluated in channel flows at various Reynolds numbers, boundary layer flows with or without pressure gradient and backward-facing step flows with different expansion ratios and Reynolds numbers. Then the grid convergence of the model is investigated. Finally, the model performance for different values of the weighting constant C_s in the mixed time scale is assessed. The results show that the proposed model reproduces the correct wall-limiting behaviour of turbulent quantities and performs well in the near-wall region of turbulent flows. The model could be expected to be adopted in hybrid Reynolds averaged Navier–Stokes/large eddy simulation methodology for complex wall-bounded flows at high Reynolds numbers.     

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

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

Flow-induced vibrations of long circular cylinders modeled by coupled nonlinear oscillators

The dynamics of long slender cylinders undergoing vortex-induced vibrations (VIV) is studied in this work. Long slender cylinders such as risers or tension legs are widely used in the field of ocean engineering. When the sea current flows past a cylinder, it will be excited due to vortex shedding. A three-dimensional time domain model is formulated to describe the response of the cylinder, in which the in-line (IL) and cross-flow (CF) deflections are coupled. The wake dynamics, including in-line and cross-flow vibrations, is represented using a pair of non-linear oscillators distributed along the cylinder. The wake oscillators are coupled to the dynamics of the long cylinder with the acceleration coupling term. A non-linear fluid force model is accounted for to reflect the relative motion of cylinder to current. The model is validated against the published data from a tank experiment with the free span riser. The comparisons show that some aspects due to VIV of long flexible cylinders can be reproduced by the proposed model, such as vibrating frequency, dominant mode number, occurrence and transition of the standing or traveling waves. In the case study, the simulations show that the IL curvature is not smaller than CF curvature, which indicates that both IL and CF vibrations are important for the structural fatigue damage. Supported by the National Natural Science Foundation of China (Grant No. 10532070), the Knowledge Innovation Program of Chinese Academy of Sciences (Grant No. KJCX2-YW-L07), and the LNM Initial Funding for Young Investigators     

绕圆柱体自由表面磁流体流动和传热的研究

本文对在不同雷诺数下,绕圆柱体的磁流体自由表面流动及传热进行了模拟,分析了磁场对绕流圆柱尾迹和涡分离的影响,获得了两种雷诺数下的电磁力密度、流场和温度场分布。结果表明,磁场不仅影响了流动的形态,而且对湍流有抑制作用,降低了自由表面的更新机制,从而减少了传热能力;在相同的Hartmann数下,相比低雷诺数下的流动换热情况,高雷诺数下的湍流不能被完全抑制,自由表面与尾迹的相互作用也较强,因而自由表面换热也较强。     

低雷诺数下钝体三维尾迹中的涡量符号律

钝体是目前各种工程中广泛应用的一种结构.钝体绕流的尾迹涡动力学也是经典的流体力学研究对象之一.本文通过直接数值模拟,针对低雷诺数下各种钝体结构的不可压缩绕流,当形成三维尾迹时,研究具有特定符号的涡量分布特征.通过分析两类钝体结构,基本的直柱体和受到几何扰动的柱体,总结并得到了更为广泛适用的涡量符号律.通过对比并分析这两类钝体结构,结合理论证明的结果,进一步厘清了对产生涡量符号律的这两类钝体结构之间的内在物理关联,即引起自然失稳的小扰动在惯性力作用下产生的表面涡量只能向下游演化发展,而几何扰动则根据扰动位置,产生的表面涡量可以向扰动上游或下游演化发展.从而可以推测所有钝体结构尾迹中的各种型式的涡脱落模态,从涡量符号律的演化角度来看,实际上是一致的,都是起源于壁面产生特定符号组合规律的∏型涡.     

3D simulation and visualization of unsteady wake flow behind a cylinder

In this paper we focus on 3D simulation of unsteady wake flow behind a circular cylinder. We show that in addition to accurate formulations and sufficiently-refined meshes, efficient computing methods are essential components of an effective simulation strategy. We use the Multi-Domain Method (MDM) we developed recently in computation of two cases. At Reynolds number 300, we demonstrate how the MDM enables us to use highly-refined meshes to capture wake patterns which we otherwise cannot fully represent. At Reynolds number 140, we show that with the MDM we can extend our computations sufficiently downstream, and with sufficient accuracy, to successfully capture the second phase of the Karman vortex street, which has been observed in laboratory experiments, and which has double the spacing between the vortices compared to the first phase.     

Modeling of complex flows and heat transfer

This paper presents the numerical modeling of complex flows and heat transfer. The Finite Analytic method is used to discretize the transport equations. The diagonal Cartesian method is proposed to model fluid flows and heat transfer over complex geometries. A three-dimensional channel flow with conjugate heat transfer is simulated. By the diagonal Cartesian method and 5-point Finite Analytic scheme, a grooved channel flow and flow in a casting bank at different Reynolds numbers are modeled. Simulations by both the diagonal Cartesian method and the traditional saw-tooth Cartesian method indicates the diagonal Cartesian method improves the modeling of flows, due to the more accurate approximation of complex boundaries. Heat transfer in two-dimensional finned compact heat exchanger is also studied. An improved heat exchanger is proposed based on the numerical prediction of heat transfer.     

通道中串列旋转圆柱的格子Boltzmann-虚拟区域方法的模拟

通过结合格子Boltzmann方法(LBM)和虚拟区域(Fictitiou sDomain)思想,建立格子Boltzmann-虚拟区域(LB-DF/FD)方法.采用两套网格系统,欧拉网格用于流体,拉格朗日网格用于固体.原有的LBM在计算运动固体的受力方面存在数据振荡,LB-DF/FD方法改进了此缺陷.为验证该方法,模拟圆柱绕流、圆形颗粒在无限长通道中平动及在无限大流场中转动三种情况,结果与其他数值解及理论解符合得很好.利用该方法模拟低雷诺数下通道中串列旋转圆柱周围的流场,分析圆柱间距(g)及雷诺数(Re)对流场结构的影响.给出Re=0.001,0.1和10下,0.2≤g≤8.0的流线结构、圆柱升力、阻力以及力矩等数值结果.结果表明,g对流场的结构及圆柱的受力有显著影响,Re对圆柱阻力及Stokes单元数目的影响较大.