高雷诺数下多柱绕流特性研究

采用改进的延迟分离涡方法数值模拟了高雷诺数下的柱体绕流,包括单圆柱绕流、单方柱绕流、串列双圆柱绕流和串列双方柱绕流,研究了不同雷诺数下圆柱绕流与方柱绕流的水动力特性.计算结果与实验数据及其他文献的数值计算结果吻合良好,研究表明,单方柱绕流在2.0×10~3Re1.0×10~7范围内未出现类似于单圆柱绕流的阻力危机现象,其平均阻力系数C_d、升力系数均方根C'_1及斯特劳哈尔数S t维持在一定范围内波动.串列双圆柱绕流与串列双方柱绕流中,均选取L/D=2.0,2.5,3.0,3.5和4.0这五中间距比进行计算.串列双圆柱绕流中,当Re=2.2×10~4时,在3.0L/D3.5内存在一临界间距比(L_c/D)使得L_c/D前后上下游圆柱的升阻力系数发生跳跃性变化,且当L/DL_c/D时,下游圆柱的阻力系数为负数.而当Re=3.0×10~6时,则不存在临界间距比,且下游圆柱的阻力系数始终为正数.串列双方柱绕流在Re=1.6×10~4和Re=1.0×10~6两种工况下的临界间距比分别处于3.0L/D3.5和3.5L/D4.0区间内,且当L/DL_c/D时,两个雷诺数下的下游方柱阻力系数均为负数.     

大湍流度高雷诺数时并列双圆柱的平均和脉动压力分布

本文通过风洞实验研究了来流湍流度,Iu=10％雷诺数分别为Re=1.95×10~9和Re=6.5×10~5时单个圆柱和不同间距比下并列双圆柱的平均和脉动压力分布。结果表明:在Re=1.95×10~5时单个圆柱的平均压力分布类似于低湍流度高超临界雷诺数时的压力分布;当雷诺数增大至6.5×10~5时,绕圆柱表面流动的分离点前移和背压绝对值提高,总的阻力系数随之增加。并列双圆柱的间距比变化对圆柱表面压力分布影响很大,在极小间距比(N/d=1.05)时,双圆柱间的缝隙流使附近柱面产生高达-5的压力系数峰值(Re=6.5×10~5),同时脉动压力也大为增加;在较小间距比时(1.5     

相同雷诺数下串列双圆柱绕流的仿真分析

为研究均匀水流场中串列排布的柱群之间的干涉影响,本文以三维串列双圆柱为例,通过计算流体力学(CFD)软件FLUENT15.0中双方程k-ε模型,分析模拟了双圆柱所受平均阻力、平均升力、后柱周向压力、斯特劳哈尔数等水动力特性。结果表明:在雷诺数为Re=2×10~4的串列双圆柱绕流中,两圆柱中心间距L与圆柱直径D的比值为L/D=4时,后柱受前柱绕流尾流影响大,明显高于单圆柱绕流的平均阻力系数,后柱的周向压力值也随前柱尾流的摆动呈现显著的不对称性;当L/D=8时,前柱绕流尾流对后柱影响逐渐减弱;当L/D=12时,两圆柱之间的相互干扰几乎可以忽略,可以看作是相互独立的单圆柱绕流。最后,计算的斯特劳哈尔数与单圆柱绕流对应的斯特劳哈尔数相近且仿真数值在计算数值范围之内,验证了整个仿真分析的准确性,也进一步说明了双圆柱绕流的柱群的干涉影响。双圆柱间距越大,前、后柱之间的干涉影响越弱。     

串列双圆柱体绕流特性与互扰效应研究

基于大涡模拟(LES)方法对亚临界雷诺数(Re=3900)下三维串列双圆柱体绕流问题进行了数值计算。首先,通过求解单圆柱算例来验证计算模型及参数的正确性。然后,着重分析了不同间距比对双圆柱体的流体力系数的影响,并阐述了双圆柱体流场特性变化及其互扰效应内在机理。研究表明:雷诺数Re=3900时,串列双圆柱体绕流临界间距比在3.9～4.0之间;随着间距比的增加,双圆柱体临近流场中二次涡团形成的区域与三维涡结构均会发生变化,导致其结构表面所受的流体力系数在时间与空间上变化的规律性逐渐减弱;达到临界间距比时,流体力系数的变化会呈现出较强的规律性。     

基于格子Boltzmann方法的低雷诺数四圆柱绕流流动模式与受力特性模拟

利用格子Boltzmann方法模拟了雷诺数为100时,均匀来流条件下的二维菱形排布的四柱绕流现象,得到了不同柱间距比下的绕流流动模式及阻力变化规律。结果表明:圆柱互扰效应与柱间距比有关,当L/D≤1.2时为单钝体模式,圆柱互扰效应以临近效应为主;当1.2相似文献   

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

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

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

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

虚拟边界法研究正交双圆柱及串列双圆球绕流

把Goldstein等人提出的虚拟边界法推广到三维情况，研究了 Re=150时不同间距下正交双圆柱绕流，和Re=250时不同间距下串列双 圆球绕流流场. 对于正交双圆柱绕流，当间距比大于3，下游圆柱对上游圆柱尾流的影响只 限定在下游圆柱的尾流所扫过的范围之内；当间距比小于等于3，下游圆柱对上游圆柱尾流 的影响扩大，下游圆柱尾流扫过区上下出现两排三维流向二次涡结构. 对于串列圆球绕流， 研究发现，在小间距比(L/D≈ 1.5)的情况下，由于上下游圆球尾流区的相互抑 制消除了压力不稳定性，整个流场呈现稳 态轴对称特征；间距比为2.0时，周向压力梯度诱发出流体的周向输运，流场呈现稳态非对 称性，但流场中存在特定的对称面；间距比增大到2.5后，绕流场开始周期振荡，原有的对 称面依旧存在；在间距比3.5时下游圆球下表面的涡结构强度有所减弱，导致占优频率发生 交替；间距比增至7.0时，整个流场恢复稳态特征，两圆球尾部同时出现双线涡，这时流场 对称面的位置发生了变动.     

NUMERICAL INVESTIGATION ON FLOW-INDUCED VIBRATION OF TWO CYLINDERS IN TANDEM ARRANGEMENTS AND ITS COUPLING MECHANISMS

对雷诺数Re= 100 条件下串列双圆柱的流致振动进行了数值模拟. 圆柱的质量比m^*均为2.0,间距比L/D 为2.0 5.0. 考虑两种工况:(a) 上游圆柱固定,下游圆柱可沿横流向自由振动;(b) 上、下游圆柱均可沿横流向自由振动. 结果表明:无论上游圆柱静止或者振动,下游圆柱横向振幅明显大于单圆柱的. 工况(b) 的下游圆柱最大振幅要大于工况(a) 的,这是由于两圆柱均振动时,两圆柱之间耦合作用增强,上游圆柱的尾流和下游圆柱的振动之间“相互调节” 作用显著. 对工况(b) 的下游圆柱振动和间隙流之间的作用机制进行了详细的研究,发现当上游圆柱脱落的自由剪切层重新附着于下游圆柱上并且完全从间隙之间通过时,下游圆柱的振幅最大.     

串列双圆柱同步测压实验结果的相关分析

本文在雷诺数2×104下,同步测量了12个不同间距下串列双圆柱的表面压力分布,积分得到脉动升、阻力的时间历程,并对前、后柱之间的脉动升、阻力以及脉动升阻力和圆柱表面的脉动压力进行了相关分析.在本实验中,串列双圆柱流态切换的临界间距比在3.5～4.0之间,在临界间距前后,相关曲线的形态差别很明显.当串列双圆柱之间的距离小于临界间距时,前后柱之间的影响显著,造成脉动升/阻力和脉动压力之间的相关程度下降.串列双圆柱之间的距离大于临界距离时,前、后柱的脉动升力之间的相位差随间距线性增加.分离点的脉动压力和脉动升力之间以及背压和脉动阻力之间都有较好的相关性,可以用分离点的脉动压力和背压来分别表征脉动升力和脉动阻力.     

流向振荡圆柱绕流的格子Boltzmann方法模拟

用一种新近发展起来的格子Boltzmann方法(LBM)在相对较小的雷诺数(Re \le 200)条件下模拟了不可压缩的流向振荡圆柱绕流问题, 考查了涡脱落模态和升阻力特性. 通过模拟, 在近尾流区发现了实验研究中已经发现的对称/反对称的涡脱落模态, 包括有些传统数值方法未发现的模态. 研究了频率锁定区域的范围及其与振幅的关系, 发现振幅越大, 发生锁定的频率区域越宽. 此外还对升阻力进行了定量意义的模拟,研究了振荡频率和振幅与升阻力的关系.      

Resolution sensitivity of momentum‐exchange and immersed boundary methods for solid–fluid interaction in the lattice Boltzmann method

In the lattice Boltzmann method (LBM), the mechanism of fluid–solid interaction can be effectively captured by appropriately enforcing the no‐slip conditions in shear direction, and bounce‐back of the non‐equilibrium distribution portion in the normal direction at fluid–solid interfaces. Among various solid–fluid interaction schemes being proposed for LBM in recent decades, two simple fluid–solid interaction methods—the momentum exchange algorithm (MEA) and the immersed boundary scheme (IBS)—were developed based on the above concept. In this paper, MEA and IBS are implemented in a D2Q9 LBGK system and applied to measure the wall correction factors of drag force upon a stationary circular particle midway in the Poiseuille channel flow at very low Reynolds number and drag coefficients at low to moderate Reynolds numbers. MEA and IBS are also employed to compare the fluid‐induced torque over the cylinder in the Taylor–Couette flow, and the steady velocity of a particle settling under the influence of gravity inside a tube. The above experiments show that IBS seems to be more accurate and less demanding on lattice resolution. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical characterization of high-temperature filament interaction with blunt cylinder at Mach 3

Energy deposition is a robust technique for various high-speed flow control applications including drag reduction. A numerical study of energy deposition is performed for perfect gas flow approaching a blunt cylinder at Mach 3. The energy deposition is simulated by a high-temperature filament injected instantaneously in front of the cylinder. The effect of important dimensionless parameters is studied to characterize the drag modification. The results indicate a saturation effect on maximum drag reduction at higher magnitudes of energy deposition. The computations reveal that the discharge location of the filament does not significantly impact the drag. A phenomenological examination of the interaction is performed. The effectiveness and efficiency of the filament on drag reduction are investigated. A one-dimensional analytical approach is studied to describe the numerical results.     

Flow around a cylinder at the beginning of the critical region

Results of experiments with a turbulent flow around a transversely aligned circular cylinder located at identical distances from the walls of a rectangular channel are reported. Data on averaged velocity fields around the cylinder are obtained by means of particle image velocimetry (PIV). Based on these fields, the near wake behind the cylinder is studied, and the kinematic characteristics for flow regimes with and without cavitation are compared. Based on the vector fields of averaged velocity, the angles of separation of the boundary layer from the cylinder surface in the considered flow regimes are determined. The drag coefficients of the cylinder for different flow regimes are calculated. It is demonstrated that the vortex region behind the cylinder and the drag coefficient of the cylinder increase in the case with cavitation. It is also shown that vortex shedding from the cylinder may be irregular, despite the fact that this process is quasi-periodic for most of the time.     

Application of the Lattice-Boltzmann Method for Particle-laden Flows: Point-particles and Fully Resolved Particles

The Lattice-Boltzmann-Method (LBM) is a powerful and robust approach for calculating fluid flows over or through complex geometries. This method was further developed for allowing the calculation of several problems relevant to dispersed particle-laden flows. For that purpose two approaches have been developed. The first approach concerns the coupling of the LBM with a classical Lagrangian procedure where the particles are considered as point-masses and hence the particles and the flow around them are numerically not resolved. As an example of use, the flow through a single pore representing a single element of a filter medium was considered and the deposition of nano-scale particles was simulated. The temporal evolution of the deposit structures is visualised and both the filtration efficiency and the pressure drop are simulated and compared with measurements. In the second developed LBM-approach, the particles are fully resolved by the numerical grid whereby the flow around particles is also captured and it is possible to effectively calculate forces on complex particles from the bounce-back boundary condition. As a case study the flow around spherical agglomerates consisting of poly-sized spherical primary particles with sintering contact is examined. Using local grid refinement and curved wall boundary condition, accurate simulations of the drag coefficient of such complex particles were performed. Especially the effect of porosity on the drag was analysed. Moreover, the flow about very porous fractal flocks, generated by a random process, was simulated for different flock size and fractal dimension. The drag coefficients resulting from LBM simulations were compared to theoretical results for Stokes flow. Finally, scenarios with moving particles were considered. First, the sedimentation of a single particle towards a plane wall was simulated and compared with measurements for validation. Secondly, the temporal sedimentation of a cluster of 13 particles was studied. Here, the primary particles were allowed to stick together and form agglomerates. This research will be the basis for further analysing agglomerate formation in laminar and turbulent flows.     

不同壁面绕流特性的格子Boltzmann模拟研究

为了探讨不同壁面的绕流特性,针对粘性流场中,不同壁面诱导的涡脱落现象以及升阻力系数等流场特性进行了格子Boltzmann数值研究。利用基于分子动理论的格子Boltzmann方法(LBM)求解Navier-Stokes方程,实现对流体运动的描述,针对不同的壁面条件,分别采用不同的格子Boltzmann流-固壁面处理方法。采用Half-way反弹边界条件来处理平直壁面,而曲壁面则采用LBM与有限差分法相结合的形式进行处理,计入了壁面与标准网格不重合对结果造成的影响。开发相应的计算程序,计算结果与已发表文献结果吻合良好,验证了数值模型的正确性。同时,探讨了进出口边界与钝体中心的距离对结果的影响。对比分析了不同壁面的绕流模型中升阻力系数、斯托罗哈数和涡量云图等,并进一步研究了雷诺数条件的影响。结果表明,不同壁面的绕流特性具有明显差异,且同时受雷诺数的显著影响;一般地,平直壁面对于来流作出的响应更迅速。     

Viscoplastic flow around a cylinder kept between parallel plates

Numerical simulations have been undertaken for the creeping pressure-driven flow of a Bingham plastic past a cylinder kept between parallel plates. Different gap/cylinder diameter ratios have been studied ranging from 2:1 to 50:1. The Bingham constitutive equation is used with an appropriate modification proposed by Papanastasiou, which applies everywhere in the flow field in both the yielded and practically unyielded regions. The emphasis is on determining the extent and shape of yielded/unyielded regions along with the drag coefficient for a wide range of Bingham numbers. The present results extend previous simulations for creeping flow of a cylinder in an infinite medium and provide calculations of the drag coefficient around a cylinder in the case of wall effects.     

Lattice Boltzmann method for simulating particle–fluid interactions

The lattice Boltzmann method （LBM） has gained increasing popularity in the last two decades as an alternative numerical approach for solving fluid flow problems. One of the most active research areas in the LBM is its application in particle-fluid systems, where the advantage of the LBM in efficiency and parallel scalability has made it superior to many other direct numerical simulation （DNS） techniques. This article intends to provide a brief review of the application of the LBM in particle-fluid systems. The numerical techniques in the LBM pertaining to simulations of particles are discussed, with emphasis on the advanced treatment for boundary conditions on the particle-fluid interface. Other numerical issues, such as the effect of the internal fluid, are also briefly described. Additionally, recent efforts in using the LBM to obtain closures for particle-fluid drag force are also reviewed.     

Active control of flow past an elliptic cylinder using an artificial neural network trained by deep reinforcement learning

The active control of flow past an elliptical cylinder using the deep reinforcement learning (DRL) method is conducted. The axis ratio of the elliptical cylinder $\Gamma$ varies from 1.2 to 2.0, and four angles of attack $\alpha=0^\circ, 15^\circ, 30^\circ$, and $45^\circ$ are taken into consideration for a fixed Reynolds number $Re=100$. The mass flow rates of two synthetic jets imposed on different positions of the cylinder $\theta_1$ and $\theta_2$ are trained to control the flow. The optimal jet placement that achieves the highest drag reduction is determined for each case. For a low axis ratio ellipse, i.e., $\Gamma=1.2$, the controlled results at $\alpha=0^\circ$ are similar to those for a circular cylinder with control jets applied at $\theta_1=90^\circ$ and $\theta_2=270^\circ$. It is found that either applying the jets asymmetrically or increasing the angle of attack can achieve a higher drag reduction rate, which, however, is accompanied by increased fluctuation. The control jets elongate the vortex shedding, and reduce the pressure drop. Meanwhile, the flow topology is modified at a high angle of attack. For an ellipse with a relatively higher axis ratio, i.e., $\Gamma\ge1.6$, the drag reduction is achieved for all the angles of attack studied. The larger the angle of attack is, the higher the drag reduction ratio is. The increased fluctuation in the drag coefficient under control is encountered, regardless of the position of the control jets. The control jets modify the flow topology by inducing an external vortex near the wall, causing the drag reduction. The results suggest that the DRL can learn an active control strategy for the present configuration.