紧致插值曲线方法在流向强迫振荡圆柱绕流中的应用

利用紧致插值曲线(constrained interpolation profile method in Zhejiang University, CIP-ZJU) 数学模型, 对低科勒冈-卡朋特(Keulegan–Carpenter) 数KC 静止流体中振荡圆柱以及雷诺数Re = 200 时流向强迫振荡圆柱绕流进行了数值模拟. 模型在直角坐标系统下建立, 采用紧致插值曲线方法作为流场的基本求解器离散了纳维-斯托克斯方程, 基于多相流的理论实现流固耦合同步求解, 利用浸入边界方法处理固体边界. 模拟结果与现有文献结果进行比较, 二者吻合情况较好, 验证了此方法对于计算复杂流动问题的可靠性.      

NUMERICAL STUDY OF THE BLOCKAGE EFFECTS ON VISCOUS FLOW PAST A CIRCULAR CYLINDER

In various numerical solutions of flow around bluff bodies the unbounded physical domain is replaced by a restricted computational one whose extent depends on the size of the computational grid network. The truncation of the solution domain in the cross-flow direction reduces the computer time required for the solution, but introduces numerical blockage effects which influence considerably the values of the various flow parameters. In the present paper the finite element solution of steady and unsteady flow around a circular cylinder at Re=106 is presented for blockage ratios of 0·05, 0·15 and 0·25. A boundary condition was tested for which the streamfunction values at the outer boundaries were those of the irrotational solution around a circular cylinder. The size of the standing vortices decreases with the blockage ratio when the flow is steady, while the spacing of the vortices decreases in both directions with increasing blockage ratio when the wake becomes unsteady. The hydrodynamic forces on the cylinder and the Strouhal number are magnified as the blockage ratio increases. The application of the streamfunction values derived from the irrotational solution at the outer boundaries reduced blockage effects only at high blockage ratio.     

均匀来流中旋转圆柱黏性绕流的数值研究

从不可压非定常Ｎ－Ｓ方程出发，首次数值求解了均匀来流中圆柱作周向旋转振荡的黏性绕流问题。探讨了旋转角速度振幅、振荡频率及Ｒｅ数等因素对流场结构及其非定常演化过程的影响，并根据计算结果，给出了在旋转振动频率。速度振幅平面内流场涡结构的分区图。     

均匀来流中旋转圆柱黏性绕流的数值研究

从不可压非定常Ｎ－Ｓ方程出发，首次数值求解了均匀来流中圆柱作周向旋转振荡的黏性绕流问题。探讨了旋转角速度振幅、振荡频率及Ｒｅ数等因素对流场结构及其非定常演化过程的影响，并根据计算结果，给出了在旋转振动频率。速度振幅平面内流场涡结构的分区图。     

APPLICATION OF THE LES METHOD TO THE OSCILLATING FLOW PAST A CIRCULAR CYLINDER

The oscillating flow of an incompressible fluid over a circular cylinder is investigated by solving the vorticity/stream-function version of the two-dimensional Navier Stokes equations using a finite-difference Large Eddy Simulation (LES) method. Two different subgrid scale (SGS) models are tested. They are the classical Smagorinsky model and Yoshizawa's Two-Scale Direct-Interaction Approximation (TSDIA) model, both of which require the input of constants in the specification of the subgrid scale viscosity. In addition, a solution was obtained with no subgrid scale model for comparison purposes. Yoshizawa's model is based on the combination of Kraichman's Direct Interaction Analysis (1964) with the separation of scales of mean (grid scale) and fluctuating (subgrid scale) fields. This TSDIA model has the length scale as a function of space and time so that it incorporates the effects of local time rate of change and advection of the resolvable field. Both models are applied to the classical problem of an oscillating two-dimensional flow over a circular cylinder at various KC values atβ =1035.Calculated lift and in-line force coefficients compare favorably with available experimental data for KC≤2 with the best solutions developed from the Yoshizawa model. The solution for KC>2 demonstrates that the two-dimensional model is inadequate to describe the flow accurately.     

旋转振动圆柱绕流周期解和Floquet稳定性

对低雷诺数旋转振动圆柱绕流问题运用低维Galerkin方法将N－S方程约化为一组非线性常微分方程组。运用打靶法数值求解了这组方程的周期解,并用Tloquet理论对周期解的稳定性进行了分析,确定了流动失稳的机制。     

A NUMERICAL SIMULATION OF VORTEX SHEDDING FROM AN OSCILLATING CIRCULAR CYLINDER

Vortex shedding from an oscillating circular cylinder is studied by numerical solutions of the two-dimensional unsteady Navier–Stokes equations. A physically consistent method is used for the reconstruction of velocity fluxes which arise from discrete equations for the mass and momentum balances. This method ensures a second-order accuracy. Two phenomena are investigated and, in both cases, the cylinder oscillation is forced. The first is the flow induced by the harmonic in-line oscillation of cylinder in water at rest. The Reynolds number is equal to 100 and the Keulegan–Carpenter number is equal to 5. A comparison of phase-averaged velocity vectors between measurements and predictions is presented. Applying the widely used model of Morison to the computed in-line force history, the drag and the inertia coefficients are calculated and compared for different grid levels. Using these to reproduce the force functions, deviations from those originally computed are revealed. The second problem is an investigation of a transversely oscillating cylinder in a uniform flow at fixed Reynolds number equal to 185. The cylinder oscillation frequency ranges between 0·80 and 1·20 of the natural vortex-shedding frequency, and the oscillation amplitude is 20% of the cylinder diameter. As the frequency of excitation of the cylinder increases relative to the inherent vortex formation frequency, the initially formed concentration of vorticity moves closer to the cylinder until a limiting position is reached. At this point, the vorticity concentration abruptly switches to the opposite side of the cylinder. This process induces distinct changes of the topology of the corresponding streamline patterns.     

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

基于格子Boltzmann方法(LBM)对不可压横向振荡圆柱绕流问题进行了数值研究. 与传统的求解宏观的N-S方程的数值方法不同, LBM求解此类问题不需要采用动网格, 而且不需要对网格进行特殊处理, 从而节约了计算成本. 结果显示, 当振荡频率增加到相应的静止圆柱绕流的自然涡脱落频率附近时, 圆柱后最新形成的集中涡距离柱体越来越近, 直到达到一个极限位置. 随后, 集中涡突然转向圆柱体另一侧脱落. 当振荡频率接近于静止圆柱的自然涡脱落频率时, 发生频率同步的现象. 随着振荡频率远离自然涡脱落频率, 同步现象消失. 在几种次谐振荡和超谐振荡下, 尾流区的涡脱落频率仍为相应的静止圆柱绕流的自然涡脱落频率.      

弹性支撑圆柱绕流稳定性分析

基于CFD 技术,采用系统辨识方法,建立了亚临界雷诺数(Re < 47) 下绕圆柱流动的非定常气动力模型(reduced order model, ROM). 耦合结构运动方程和降阶气动力模型,建立了弹性支撑圆柱绕流的稳定性分析模型. 算例分析了亚临界雷诺数下,结构固有频率、质量比等参数以及支撑方式对弹性系统稳定性的影响. 对于单自由度横向支撑圆柱,当结构固有频率趋近流动最不稳定模态频率时,弹性系统会在一定频率范围内失稳,这种现象最低可在Re~20 时出现. 旋转自由度的释放能够进一步降低系统的稳定性,可将临界雷诺数进一步降低至18 左右. ROM 方法不仅具有很高的效率,而且清晰地指出了弹性系统失稳的根本原因:流动模态和结构模态耦合作用导致结构模态失稳所致. 因此,失稳状态下系统振荡频率锁定于结构固有频率. 基于ROM 技术预测的失稳边界与直接CFD/CSD 仿真结果吻合,证明了该方法的正确性和精度.      

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

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

VORTEX SHEDDING RESONANCE FROM A ROTATIONALLY OSCILLATING CYLINDER

Vortex shedding resonance of a circular cylinder wake to a forced rotational oscillation has been investigated experimentally by measuring the velocity fluctuations in the wake, pressure distributions over the cylinder surface, and visualizing the flow field with respect to cylinder oscillations. The vortex shedding resonance occurs near the natural shedding frequency at small amplitude of cylinder oscillations, while the peak resonance frequency shifts to a lower value with an increase in oscillation amplitude. The drag and lift forces acting on the cylinder at fixed forcing Strouhal number indicate that the phase lag of fluid forces to the cylinder oscillations increases with an increase in oscillation amplitude, supporting the variation of resonance frequency with oscillation amplitude. The comparative study of the measured pressure distributions and the simultaneous flow visualizations with respect to cylinder rotation shows the mechanisms of phase lag, which is due to the strengthened vortex formation and the modification of the surface pressure distributions.     

流向强迫振荡圆柱绕流的涡脱落模态分析

通过求解采用ALE方法描述的运动坐标系Navier-Stokes方程组,分析均匀来流下雷诺 数为150的静止和流向振荡的圆柱绕流. 主要研究了强迫振荡频率和较大振幅比 (A/D=0.3-1.2)对圆柱升力、阻力变化特性以及涡脱落模态的影响. 研究表 明,流向振荡圆柱绕流存在多种涡脱落模态,如对称S以及反对称A-I, A-III, A-IV等多种形式;比较研究结果,拓展了各模态下对应的锁定区域,并将其分为5个 子区;A-I模态中圆柱受力较以前所知更复杂;通过分析计算结果,发现最大加速度 比Af_{c}^{2}/Df_{s0}^{2}可能是涡脱落模态(尤其是对称S模态)最有效的控制参数.     

NUMERICAL STUDY OF HIGH-REYNOLDS-NUMBER FLOW PAST A BLUFF OBJECT

A semi-explicit finite difference scheme is proposed to study unsteady two-dimensional, incompressible flow past a bluff object at high Reynolds number. The bluff object comes from a class of elliptical cylinders in which the aspect ratio and the angle of attack are two controlled parameters. Associated with the streamfunction–vorticity formulation, the interior vorticity, streamfunction and wall vorticity are updated in turn for each time step. The streamfunction and wall vorticity are solved by means of a multigrid method and a projection method respectively. In regard to the vorticity transport equation, implicitness is merely associated with the diffusion operator, which can be made semi-explicit via approximate factorization. Low-diffusive upwinding is devised to handle the convection part. Numerical results are reported for Reynolds numbers up to 40,000. Comparisons with other numerical or physical experiments are included.     

均匀来流中横向振动圆柱近迹涡结构的数值模拟

从涡量流函数形式的N-S方程出发,在不同的振动频率、振幅及Re数下数值研究了均匀来流中横向振动圆柱粘性统流的涡脱泻现象。着重探讨了近迹复杂的涡结构及其非定常演化过程,以及它们对物体受力特性的影响,并首次成功地模拟了近年来实验研究中所发现的一些重要的流动现象,如相位“开关”现象(phase switch phenomena)及复杂涡结构(complex vortex structure),等等。通过数值模拟,不仅能够再现实验研究中所发现的一些重要的流动现象,还可进一步预示某些新的流动现象,使数值计算起到与实验研究相辅相成的作用。     

流向磁场作用下圆柱绕流的直接数值模拟

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.      

圆柱绕流的非线性动力学

本文利用非线性动力学的概念和方法来研究粘性不可压缩流体的二维圆柱绕流问题．计算了定常流失稳以及出现混沌的临界雷诺数，并估计了混沌解的一些统计特征．     

小Kc数振荡流中圆柱绕流流场结构的数值分析

本文综合Chorin的涡方法、网格涡方法和泊松方程快速解法的优点,采用固定离散涡的数目和位置,随时间变化调整离散涡强度的计算方法,对Kc数在1到5之间,β数在10到1665之间圆柱在振荡流动中的流场结构进行了数值模拟。并与U型振荡水槽中的流场涡结构显示照片比较,计算结果表明,数值解在稳定性和收敛性方面都是比较好的。与显示照片相比,其符合程度也是令人满意的。     

DIRECT NUMERICAL SIMULATIONS ON THE TURBULENT FLOW PAST A CONFINED CIRCULAR CYLINDER WITH THE INFLUENCE OF THE STREAMWISE MAGNETIC FIELDS1)

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.     

幂律流体圆柱绕流的格子波尔兹曼模拟

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明:插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.