圆管湍流减阻电磁力控制的直接数值模拟

利用直接数值模拟研究了圆管湍流减阻的电磁力控制问题. 在圆管表面布置一层交替排 列的条形电极和磁极，施加电压后在导电的流体介质中产生周向随时间变化的电磁力，通过 改变电磁力的大小和变化周期等参数可以使阻力大幅减少. 研究减阻率和控制参数的 关系，并研究了电磁力减阻的机理. 电磁力控制存在最优周期，阻力最大降 低24.2{\%}; 电磁力控制抑制了流向和法向速度的脉动，而周向脉动速度随振荡周 期的增大而增大.     

流体边界层上电磁力的控制效应研究

利用作用于流体边界层上的电磁体积力改变流体边界层的结构，研究电磁力对流场的控制 作用效果. 电极与磁极交替分布的电磁场激活板包覆在圆柱体表面置于流动的电解质溶液 中，产生的电磁力沿圆柱体表面分布，可以改变流体边界层的结构，从而实现对流场的控制. 用电磁屏蔽和时域控制的方法调整电磁力的时空分布参数，圆柱绕流分离点可以在前驻点和 后驻点之间变动，产生不同的控制效果. 流体边界层上的电磁力能连续控制圆柱绕流、尾流 涡街的形态. 正向电磁力具有较好的消涡、减震和减阻控制效应. 反向电磁力具有明显的增 涡控制效应，具有较强的制动控制效应，此时圆柱体表面涡量分布的对称性和稳定性被破坏.     

圆柱表面包覆电磁场消涡与增涡实验研究

利用电磁体积力改变流体边界层的结构，用作用于流体边界层上的电磁力进行消涡与增涡控制。交替分布的电极和磁极包覆在圆柱体的表面置于电介质溶液中，简单调整电磁力的时间与空间分布可以方便地控制圆柱绕流的形态。通过理论分析和数值模拟确定了实验控制的关键参数，实现了电磁力消涡和增涡的连续控制，电磁力作用下的圆柱绕流的分离点可以在前驻点和后驻点之间变动。     

并列弹性双圆柱在均匀流场中的流体激励力的研究

实验研究了并列弹性双圆柱体在均匀流场中的水动力特性，其中包括柱表面稳态及脉动压力系数，稳态阻力及升力系数，脉动阻力和升力系数，结果表明：Ｖ＿ｒ＞Ｖ＿（ｒｃ）时，柱的振动对其流体激励力有明显影响。     

带有新型涡激振动抑制罩的圆柱体的水动力特性

通过模型实验和数值模拟计算,研究了带有涡激振动抑制罩的圆截面柱体的水动力特性.模型实验主要测试了柱体上附加谐波型和类圆锥型涡激振动抑制罩的单摆结构在不同流速下发生涡激振动的性质;数值模拟则针对谐波型和圆锥型扰动,在雷诺数Re为102到105范围内,研究其水动力参数,如阻力、升力和涡脱落频率等,随扰动波长和波动强度的变化.模型实验结果表明,在直圆柱开始发生共振的流速下,带抑制罩的柱体的振幅显著降低,而在更高流速下则显著增大.数值模拟结果表明,谐波型和圆锥型扰动具有相似的水动力特性;且在不同Re时,阻力、升力和涡脱落频率具有相似的变化规律;随波动强度的增大,阻力一般逐渐增大,升力则在多数情况下先减小而后增大,而涡脱落频率一般逐渐减小.      

利用O型环抑制圆柱绕法流涡脱落的数值研究

圆柱绕流涡脱落诱发较大的振动和声,如何有效地抑制值得关注.利用大涡模拟技术求解了Navier-Stokes方程,得到了涡脱落频率,升力脉动幅值及平均阻力系数.计算表明二维模拟不能体现流动基本特征,三维计算与实验吻合较好.为了抑制涡脱落,在直径为D的圆柱表面装入间距为1D,直径为0.0167D的O型环.通过升力、速度谱分析以及柱向横截面流场分析可知,在光滑圆柱外表面加入O型环能诱发流体边界层分离,有效地抑制涡脱落现象,升力脉动和观测点速度脉动幅值几乎完全消失,阻力系数也略微降低,适合在实际工程中采用.     

正方形排列四圆柱体绕流特性及互扰效应研究

基于计算流体动力学理论,运用大涡模拟方法对雷诺数Re=3900三维正方形排列四圆柱体结构群的绕流问题进行数值计算,主要分析来流攻角与间距比两个参数对四圆柱体结构群流体参数及流场模态的影响。结果表明:来流攻角与间距比均对四圆柱体结构群绕流特性有较强的影响;来流攻角θ=0°、22.5°、45°下,临界间距比分别为3.5、4.0、3.0;间距比的变化会导致下游圆柱表面压力系数分布发生改变;另一方面,间距比较小时,四圆柱体结构之间的互扰作用均以临近效应为主;随间距比增大,上游圆柱尾流对下游圆柱有显著影响,其互扰作用会转变尾激效应。     

基于浸入边界-格子Boltzmann通量求解法的椭圆柱流动特性分析

基于浸入边界-格子Boltzmann通量求解法,开展了雷诺数Re=100不同几何参数下单椭圆柱及串列双椭圆柱绕流流场与受力特性对比研究。结果表明,随长短轴比值的增加,单椭圆柱绕流阻力系数先减小后缓慢上升,最大升力系数则随长短轴比值的增大而减小;尾迹流动状态从周期性脱落涡到稳定对称涡。间距是影响串列圆柱及椭圆柱流场流动状态的主要因素,间距较小时,串列圆柱绕流呈周期性脱落涡状态,而椭圆柱则为稳定流动;随着间距增加,上下游圆柱及椭圆柱尾迹均出现卡门涡街现象,且串列椭圆柱临界间距大于串列圆柱。串列椭圆柱阻力的变化规律与圆柱的基本相同,上游平均阻力大于下游阻力;上游椭圆柱阻力随着间距的变大先减小,下游随间距的变大而增加,当间距达到临界间距时上下游阻力跃升,随后出现小幅度波动再逐渐增加,并趋近于相同长短轴比值下单柱体绕流的阻力。     

非结构隐式方法在空心弹丸流场模拟中的应用

采用非结构网格的LU-SGS隐式算法计算三维Euler方程,数值模拟了不同马赫数以及不同攻角下某空心弹丸绕流流场,分析了流场的波系结构及其升阻力特性,计算结果表明空心弹丸的阻力系数比同口径的普通弹丸的阻力系数大约小30％,空心弹丸的阻力系数以及升力系数随攻角的变化规律与普通弹丸一致．     

基于浸入边界-格子Boltzmann通量求解法的椭圆柱流动特性分析

基于浸入边界-格子Boltzmann通量求解法,开展了雷诺数Re=100不同几何参数下单椭圆柱及串列双椭圆柱绕流流场与受力特性对比研究。结果表明,随长短轴比值的增加,单椭圆柱绕流阻力系数先减小后缓慢上升,最大升力系数则随长短轴比值的增大而减小;尾迹流动状态从周期性脱落涡到稳定对称涡。间距是影响串列圆柱及椭圆柱流场流动状态的主要因素,间距较小时,串列圆柱绕流呈周期性脱落涡状态,而椭圆柱则为稳定流动;随着间距增加,上下游圆柱及椭圆柱尾迹均出现卡门涡街现象,且串列椭圆柱临界间距大于串列圆柱。串列椭圆柱阻力的变化规律与圆柱的基本相同,上游平均阻力大于下游阻力;上游椭圆柱阻力随着间距的变大先减小,下游随间距的变大而增加,当间距达到临界间距时上下游阻力跃升,随后出现小幅度波动再逐渐增加,并趋近于相同长短轴比值下单柱体绕流的阻力。     

Electromagnetic control of seawater flow around circular cylinders

We investigate numerically the electromagnetic control of seawater flows over an infinitely long circular cylinder. Stripes of electrodes and magnets, wrapped around the cylinder surface, produce a tangential body force (Lorentz force) that stabilizes the flow. This mechanism delays flow separation, reduces drag and lift, and finally suppresses the von Kármán vortex street. Results from two-dimensional simulations of the Navier–Stokes equations in a range 10<Re<300 and Lorentz force calculations are presented. Emphasis is placed on the disclosure of physical phenomena as well as a quantitative detection of the flow field and forces. It is shown that the drag strongly depends on the geometry of the electromagnetic actuator and on its location at the cylinder surface. The effect of flow control increases with larger Reynolds numbers, since the boundary layer thickness and the penetration depth of the Lorentz force are closely connected.     

电磁力控制翼型绕流分离的增升减阻效率研究

电磁力可有效对流体流动进行控制,增升减阻,抑制流动分离,制约其推广应用的瓶颈为控制效率问题.为提高其控制效率,基于翼型绕流的电磁力控制,对电磁力增升减阻的控制效率问题进行数值研究. 根据能量守恒定律,推导电磁力控制能耗的比,基于升力和阻力计算节省能量. 定义电磁力的控制效率为能量节省与电磁力控制所需能耗的比值,研究不同工况下电磁力增升减阻的控制效率. 发现在控制开始阶段,电磁力能量主要消耗在增加边界层流体的动能上,电磁力控制效率非常低,但电磁力控制效率会随着电磁力工作时间的增长而增加;电磁力控制效率随着来流速度的增加呈指数下降;通过增加电磁力激活板的输入能量可增强电磁力的控制效果,但无法明显增加其控制效率.      

包覆电磁场激活板的圆柱尾迹的数值研究

利用电磁场作用于电介质溶液的Lorentz力可以控制溶液的流动.对置于弱电介质溶液中,包覆电磁场激活板的圆柱周围的绕流进行数值研究,讨论了电磁场激活板的安装位置,激活方式等对圆柱尾迹的影响.     

Vortex shedding and aerodynamic forces on a circular cylinder in linear shear flow at subcritical Reynolds number

Vortex shedding and aerodynamic forces on a circular cylinder in a linear shear flow with its axis normal to the plane of the velocity shear profile at subcritical Reynolds number are investigated experimentally. The shear parameter β, which is based on the velocity gradient, cylinder diameter and upstream mean velocity at the center plane of the cylinder, varies from 0 to 0.27. The Strouhal number has no significant variation with the shear parameter. The time-mean base pressure increases and the fluctuating component of the base pressure decreases significantly with increasing shear parameter. Vortex shedding is suppressed by the velocity shear. Dislocation of the stagnation point takes place and this influences the pressure distribution around the cylinder together with the velocity shear. A mean lift force arises in the shear flow due to asymmetry of the pressure distribution, and it acts from the high velocity side to the low velocity side. In addition, the lift coefficient increases and the drag coefficient decreases with increasing shear parameter.     

Interference effect of two equal‐sized square cylinders in tandem arrangement: With planar shear flow

Numerical simulations have been performed for flow past two equal‐sized square cylinders in tandem arrangement subjected to incoming planar shear flow. Effect of L/d ratio and the shear parameter has been studied. The range of L/d ratio (ratio of center‐to‐center distance (L) to cylinder width (d)) is varied from 2 to 7 and the non‐dimensional shear parameter (K) is varied from 0.0 to 0.4 in steps of 0.1. For all the cases the Reynolds number (Re) based on centerline velocity and cylinder width is fixed at 100. The results are compared with that of isolated square cylinder with uniform flow. Strouhal number decreases with increasing shear parameter. There are more than one shedding frequency at high shear parameters and L/d ratios. The mean drag coefficient is decreased with shear parameter and lesser than that of the single cylinder. The root mean square (RMS) value of both lift and drag coefficients is higher for the downstream cylinder for all values of shear parameter. With increasing L/d ratio, for both lift and drag, the RMS value increases and then decreases for upstream cylinder, whereas it continuously increases for the downstream cylinder. The stagnation point is moved towards the top leading edge with increasing shear. The critical L/d ratio, which is defined as the distance between two cylinders, beyond which the vortex shedding from the upstream cylinder occurs, decreases with increasing shear parameter. Copyright © 2007 John Wiley & Sons, Ltd.     

Instability of cylinder wake under open-loop active control

Instability of a wake controlled by a streamwise Lorentz force is investigated through a Floquet stability analysis. The streamwise Lorentz force, which is a two-dimensional control input created by an electromagnetic actuator located on the cylinder surface,adjusts the base flow to affect the three-dimensional wake instability and achieve wake stabilization and transition delay. The instability mode at a Reynolds number Re = 300 can be transformed from B to A with N = 1.0, where N is an interaction number representing the strength of the Lorentz force relative to the inertial force in the fluid. The wake flow is Floquet stable when N increases to 1.3. The spanwise perturbation wavelengths are 3.926 D and 0.822 D in the modes A and B, respectively, where D is the cylinder diameter. In addition, the oscillating amplitudes of drag and lift are reduced with the increase in the interaction number. Particle tracing is used to explore the essential physical mechanism for mode transformation. The path lines show that suppression of flow separation hinders the fluid deformation and rotation, leading to the decrease in elliptic and hyperbolic instability regions, which is the material cause of mode transformation.All of the results indicate that wake stabilization and transition delay can be achieved under open-loop active control via the streamwise Lorentz force.     

Numerical simulation of Reynolds number effects on velocity shear flow around a circular cylinder

Three-dimensional Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed to investigate the shear effects on flow around a circular cylinder at Reynolds numbers of Re=60–1000. The shear parameter, β, which is based on the velocity gradient, cylinder diameter and upstream mean velocity at the center plane of the cylinder, varies from 0 to 0.30. Variations of Strouhal number, drag and lift coefficients, and unsteady wake structures with shear parameter are studied, along with their dependence on Reynolds number. The presented simulation provides detailed information for the flow field around a circular cylinder in shear flow. This study shows that the Strouhal number exhibits no significant variation with shear parameter. The stagnation point moves to the high-velocity side almost linearly with shear parameter, and this result mainly influences the aerodynamic forces acting on a circular cylinder in shear flow. Both the Reynolds number and shear parameter influence the movement of the stagnation point and separation point. Mode A wake instability is suppressed into parallel vortex shedding mode at a certain shear parameter. The lift force increases with increasing shear parameter and acts from the high-velocity side to the low-velocity side. In addition, a simple method to estimate the lift force coefficient in shear flow is provided.     

双柱体绕流中高度变化对其阻力的影响

对三维槽道内双柱体可压绕流进行了大涡模拟,揭示了绕流过程中柱体表面分离涡的生成、扩散与相互作用过程,并且数值模拟了两个柱体高度(阻塞比)的变化对整个流场的影响以及两柱体阻力系数的变化情况。结果表明,当两个柱体高度同时增加时,上游柱体阻力系数的大小以及下游柱体阻力系数的振幅都急剧变大,这是工程领域中所不期望的。而当仅增加下游柱体高度时,上游柱体阻力系数会略有降低,下游柱体阻力系数虽有增加,但仍小于前种情况的上游柱体阻力系数,且其振幅相对较小,因而有利于改善两柱体总体受力情况。     

Numerical study of the suppression mechanism of vortex-induced vibration by symmetric Lorentz forces

In this paper, the electro-magnetic control of vortex-induced vibration (VIV) of a circular cylinder is investigated numerically based on the stream function–vorticity equations in the exponential–polar coordinates attached on the moving cylinder for Re=150. The effects of the instantaneous wake geometries and the corresponding cylinder motion on the hydrodynamic forces for one entire period of vortex shedding are discussed using a drag–lift phase diagram. The drag–lift diagram is composed of the upper and lower closed curves due to the contributions of the vortex shedding but is magnified, translated and turned under the action of the cylinder motion. The Lorentz force for controlling the vibration cylinder is classified into the field Lorentz force and the wall Lorentz force. The symmetric field Lorentz force will symmetrize the flow passing over the cylinder and decreases the lift oscillation, which, in turn, suppresses the VIV, whereas the wall Lorentz force has no effect on the lift. The cylinder vibration increases as the work performed by the lift dominates the energy transfer. Otherwise, the cylinder vibration decreases. If the net transferred energy per motion is equal to zero, the cylinder will vibrate steadily or be fixed.     

Numerical investigation of fluid flow past circular cylinder with multiple control rods at low Reynolds number

Laminar flow past a circular cylinder with multiple small-diameter control rods is numerically investigated in this study. The effects of rod-to-cylinder spacing ratio, rod and cylinder diameter ratio, cylinder Reynolds number, number of control rods and angle of attack on the hydrodynamics of the main circular cylinder are investigated. Four different flow regimes are identified based on the mechanism of lift and drag reduction. The range of rod-to-cylinder spacing ratio where significant force suppression can be achieved is found to become narrower as the Reynolds number increases in the laminar regime, but is insensitive to the diameter ratio. The numerical results for the case with six identical small control rods at Re=200 show that the lift fluctuation on the main cylinder can be suppressed significantly for a large range of spacing ratio and various diameter ratios, while the drag reduction on the main cylinder is also achieved simultaneously. The six-control-rod arrangement has shown better performance in flow control than the arrangements with less control rods, especially in terms of force reduction at various angles of attack.