Aerodynamic force and flow structures of two airfoils in flapping motions

Aerodynamic force and flow structures of two airfoils in a tandem configuration in flapping motions are studied, by solving the Navier-Stokes equations in moving overset grids. Three typical phase differences between the fore- and aftairfoil flapping cycles are considered. It is shown that: (1) in the case of no interaction (single airfoil), the time average of the vertical force coefficient over the downstroke is 2.74, which is about 3 times as large as the maximum steady-state lift coefficient of a dragonfly wing; the time average of the horizontal force coefficient is 1.97, which is also large. The reasons for the large force coefficients are the acceleration at the beginning of a stroke, the delayed stall and the “pitching-up” motion near the end of the stroke. (2) In the cases of two-airfoils, the time-variations of the force and moment coefficients on each airfoil are broadly similar to that of the single airfoil in that the vertical force is mainly produced in downstroke and the horizontal force in upstroke, but very large differences exist due to the interaction. (3) For in-phase stroking, the major differences caused by the interaction are that the vertical force on FA in downstroke is increased and the horizontal force on FA in upstroke decreased. As a result, the magnitude of the resultant force is almost unchanged but it inclines less forward. (4) For counter stroking, the major differences are that the vertical force on AA in downstroke and the horizontal force on FA in upstroke are decreased. As a result, the magnitude of the resultant force is decreased by about 20 percent but its direction is almost unchanged. (5) For 90°-phase-difference stroking, the major differences are that the vertical force on AA in downstroke and the horizontal force on FA in upstroke are decreased greatly and the horizontal force on AA in upstroke increased. As a result, the magnitude of the resultant force is decreased by about 28% and it inclines more forward. (6) Among the three cases of phase angles, inphase flapping produces the largest vertical force (also the largest resultant force); the 90°-phase-difference flapping results in the largest horizontal force, but the smallest resultant force.     

Investigation of flow field of clap and fling motion using immersed boundary coupled lattice Boltzmann method

This paper deals with the investigation of flow field due to clap and fling mechanism using immersed boundary coupled with lattice Boltzmann method. The lattice Boltzmann method (LBM), an alternative to Navier–Stokes solver, is used because of its simplicity and computational efficiency in solving complex moving boundary problems. Benchmark problems are simulated to validate the code, which is then used for simulating flow over two elliptic wing of aspect ratio 5 performing clap and fling flapping motion for different flow parameters such as Reynolds number (Re=75, 100, 150), advance ratio (J=10E−3.10E−2,0.2) and frequency (f=0.05 Hz, 0.25 Hz). Numerical simulation is able to capture typical low Reynolds number unsteady phenomena such as, ׳wake vortex wing interaction׳, ׳Kramer effect׳ and ׳delayed stall׳. The results are both qualitatively and quantitatively consistent with experimental observation. The parametric study involving different combinations of Re, f and J depict distinctly different aerodynamic performances providing physical insights into the flow physics. It is observed that a combination of low f, low J and high Re flow results in better aerodynamic performance. Pronounced lift enhancement via leading edge vortices are obtained in unsteady regime (J<1) compared to quasi-steady regime (J>1). The role of leading edge vortices in enhancing lift are investigated by studying the size and strength of these vortices for different flow conditions. For a given Re, the magnitude of maximum lift coefficient decreases with increasing f irrespective of the value of J; while the same is enhanced with the increasing Re.     

翼型、机翼非定常运动模式下动态数值模拟

基于RANS方程,通过刚性动网格技术实现对翼型和机翼典型运动模式的描述,采用双时间推进方法和Roe空间离散格式对流场求解,构建了一个非定常气动计算平台;以NACA0012翼型为算倒进行了动态数值模拟可信度验证。数值模拟结果与试验数据吻合较好,升力和俯仰力矩的最大计算误差分别为3％和10％,表明了该平台的可靠性。另外,还数值模拟了M6机翼的动态非定常流场,并分析了两种湍流模型对非定常流场激波的捕捉能力。结果表明非定常流动中S-A湍流模型对激波的捕捉较B-L模型更敏感。文中开发的非定常计算平台对进一步解决三维复杂流场的流动问题有很高的工程应用价值。     

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

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

Large eddy simulation and spectrum analysis of flow around two square cylinders arranged side by side

Large eddy simulation cooperated with the second order full extension ETG ( Euler-Taylor-Galerkin ) finite element method was applied to simulate the flow around two square cylinders arranged side by side at a spacing ratio of 1.5. The second order full extension ETG finite element method was developed by Wang and He. By means of Taylor expansion of terms containing time derivative, time derivative is replaced by space derivative. The function of it is equal to introducing an artificial viscosity term. The streamlines of the flow at different moments were obtained. The time history of drag coefficient, lift coefficient and the streamwise velocity on the symmetrical points were presented. Furthermore, the symmetrical problem of the frequency spectrum of flow around two square cylinders arranged side by side were studied by using the spectral analysis technology. The data obtained at the initialstage are excluded in order to avoid the influence of initial condition on the results. The power spectrums of drag coefficient, lift coefficient, the streamwise velocity on the symmetrical points were analyzed respectively. The results show that although the time domain process of dynamic parameters is nonsymmetrical, the frequency domain process of them is symmetrical under the symmetrical boundary conditions.     

Three-dimensional simulations of flow past two circular cylinders in side-by-side arrangements at right and oblique attacks

The flow past two identical circular cylinders in side-by-side arrangements at right and oblique attack angles is numerically investigated by solving the three-dimensional Navier–Stokes equations using the Petrov–Galerkin finite element method. The study is focused on the effect of flow attack angle and gap ratio between the two cylinders on the vortex shedding flow and the hydrodynamic forces of the cylinders. For an oblique flow attack angle, the Reynolds number based on the velocity component perpendicular to the cylinder span is defined as the normal Reynolds number Re_N and that based on the total velocity is defined as the total Reynolds number Re_T. Simulations are conducted for two Reynolds numbers of Re_N=500 and Re_T=500, two flow attack angles of α=0° and 45° and four gap ratios of G/D=0.5, 1, 3 and 5. The biased gap flow for G/D=0.5 and 1 and the flip-flopping bistable gap flow for G/D=1 are observed for both α=0° and 45°. For a constant normal Reynolds number of Re_N=500, the mean drag and lift coefficients at α=0° are very close to those at α=45°. The difference between the root mean square (RMS) lift coefficient at α=0° and that at α=45° is about 20% for large gap ratios of 3 and 5. From small gap ratios of 0.5 and 1, the RMS lift coefficients at α=0° and 45° are similar to each other. The present simulations show that the agreement in the force coefficients between the 0° and 45° flow attack angles for a constant normal Reynolds number is better than that for a constant total Reynolds number. This indicates that the normal Reynolds number should be used in the implementation of the independence principle (i.e., the independence of the force coefficients on the flow attack angle). The effect of Reynolds number on the bistable gap flow is investigated by simulating the flow for Re_N=100–600, α=0° and 45° and G/D=1. Flow for G/D=1 is found to be two-dimensional at Re_N=100 and weak three-dimensional at Re_N=200. While well defined biased flow can be identified for Re_N=300–600, the gap flow for Re_N=100 and 200 changes its biased direction too frequently to allow stable biased flow to develop.