The circular disk straddling the interface of a two-phase flow

The motion is determined for a thin circular disk straddling the plane interface of an immiscible two phase creeping flow and moving parallel or perpendicular to the interface. Expressions are derived for the drag coefficient on the disk.     

Reduction in drag and vortex shedding frequency through porous sheath around a circular cylinder

A numerical study on the laminar vortex shedding and wake flow due to a porous‐wrapped solid circular cylinder has been made in this paper. The cylinder is horizontally placed, and is subjected to a uniform cross flow. The aim is to control the vortex shedding and drag force through a thin porous wrapper around a solid cylinder. The flow field is investigated for a wide range of Reynolds number in the laminar regime. The flow in the porous zone is governed by the Darcy–Brinkman–Forchheimer extended model and the Navier–Stokes equations in the fluid region. A control volume approach is adopted for computation of the governing equations along with a second‐order upwind scheme, which is used to discretize the convective terms inside the fluid region. The inclusion of a thin porous wrapper produces a significant reduction in drag and damps the oscillation compared with a solid cylinder. Dependence of Strouhal number and drag coefficient on porous layer thickness at different Reynolds number is analyzed. The dependence of Strouhal number and drag on the permeability of the medium is also examined. Copyright © 2009 John Wiley & Sons, Ltd.     

Influence of curvature on drag reduction by opposition control in turbulent flow along a thin cylinder

Opposition controlled fully developed turbulent flow along a thin cylinder is analyzed by means of direct numerical simulations. The influence of cylinder curvature on the skin-friction drag reduction effect by the classical opposition control (i.e., the radial velocity control) is investigated. The curvature of the cylinder affects the uncontrolled flow statistics; for instance, skin-friction coefficient increases while Reynolds shear stress (RSS) and turbulent intensity decrease. However, the control effect in the case of a small curvature is similar to that in channel flow. When the curvature is large, the maximum drag reduction rate decreased. However, the optimal location of the detection plane is the same as that in a flat plate. Further, the drag reduction effect is achieved even on a high detection plane where the drag increases in the flat plate. Although a difference in the drag reduction effect can be observed with a change in the curvature, its mechanism considered in this analysis based on the transport of the Reynolds stress is similar to that of the flat plate.     

New formula for drag coefficient of cylindrical particles

The drag force on a cylindrical particle is calculated using lattice Boltzmann method.The results show that the drag coefficient of a particle with different orientation angles decreases with increasing Reynolds number.When the principal axis of the particle is parallel to flow,the drag coefficient is much larger than that of others and decreases fastest with increasing Reynolds number,which becomes more obvious with increasing particle aspect ratio.When the principal axis of the particle is inclined to flo...     

Self‐propulsion of oscillating wings in incompressible flow

In this paper, we show that the oscillatory motion of an airfoil (wing) in an ideal fluid can determine the apparition of thrust. In the framework of the linearized perturbation theory, the pressure jump over the oscillating wing is the solution of a two‐dimensional hypersingular integral equation. Using appropriate quadrature formulas, we discretize the oscillatory lifting surface integral equation in order to obtain the jump of the pressure across the surface. Integrating numerically, we obtain the drag coefficient. For some oscillatory motions, if the frequency of the oscillations surpasses a certain value, the drag coefficient becomes negative, i.e. there appears a propulsive force. Copyright © 2007 John Wiley & Sons, Ltd.     

平板大攻角绕流升力和阻力系数的计算

二维平板或二维对称薄翼型大攻角绕流升力和阻力系数与攻角之间存在的函数关系一般用数据表格的形式给出。本文根据垂直平板绕流阻力实验数据和对称薄翼型全攻角绕流实验数据,分析得到了平板大攻角绕流总压力及其升力分量和阻力分量系数的近似计算公式。结果表明：平板总压力系数约等于攻角正弦值的2倍;总压力的阻力分量系数约等于攻角正弦值平方的2倍;升力分量系数约为攻角2倍的正弦值。计算结果与两组试验数据具有较好的一致性。     

Variation of drag coefficients in an interacting drop stream

Variation of the drag coefficient of closely spaced drops in a stream injected into a turbulent flow is studied experimentally. Three different regions are identified. In the first region, close to the injector, drops flow in the wake of each other. The drag coefficient in this region is much smaller (by a factor of 4 to 5) than the standard drag coefficient, and its magnitude is dependent on the drop initial spacing. Shortly downstream of the injection point, the transition region starts, where the drag coefficient increases rapidly approaching the drag coefficient of a single isolated drop. And further down-stream when drops are dispersed significantly, the drag coefficient will behave the same way as a single isolated drop.     

钝后缘风力机翼型的环量控制研究

钝后缘风力机翼型具有结构强度高、对表面污染不敏感等优点,但其较大的阻力系数使得翼型的整体气动特性不够理想. 利用环量控制方法对钝后缘风力机翼型进行了流动控制,以改善钝后缘风力机翼型的气动特性,减弱尾迹区脱体涡强度. 通过对钝后缘风力机翼型环量控制方法进行相关的数值模拟,对比研究了环量控制方法的增升减阻效果, 研究了环量控制下翼型升阻力特性随射流动量系数的变化规律,并对不同射流动量系数下环量控制方法的气动品质因子和控制效率进行了分析. 研究结果表明:环量控制方法能够大幅提升钝后缘风力机翼型的升力系数,同时有效地降低翼型的阻力系数; 翼型的升力系数随射流动量系数的增大而增大,表现出很明显的分离控制阶段和超环量控制阶段的变化规律; 射流能耗的功率系数随射流动量系数的增大而增大,且增长速率逐渐增大;实施环量控制方法后叶片的输出功率同样随射流动量系数增大而增大,但增长速率逐渐降低. 总体来说,环量控制方法可以有效地改善钝后缘风力机翼型的气动特性以及功率输出特性,在大型风力机流动控制中具有很好的应用前景.      

Numerical study of flow-induced oscillations of two rigid plates elastically hinged at the two ends of a stationary plate in a cross-flow

The flow-induced oscillation (FIO) of bluff bodies is commonly encountered in the fluid structure interaction (FSI) problems. In this study, we use an unstructured moving grid strategy and simulate the FIO of two rigid plates, which are elastically hinged at the two ends of a fixed flat plate in a cross-flow. We use a hybrid finite-element-volume (FEV) method in an arbitrary Lagrangian–Eulerian (ALE) framework to study FIO of the two hinged plates. The current simulations are carried out for wide ranges of flow Reynolds number (50–175), spring stiffness coefficient, and the two hinged plates' moment of inertia magnitudes. The influences of these parameters are investigated on the magnitudes of maximum deflection angle, the amplitude of oscillation, the total lift and drag coefficients, and so on. The study is also carried out in the transition period to describe the in-phase and out-of-phase angular oscillations occurring for the two elastically hinged plates with respect to each other. After the transition period, the two hinged plates eventually arrive to a similar periodic oscillation; however, with some phase lags. We find that the achieved phase lag is equal to the phase lag between the two pairs of flow vortices, which are alternatively shed into the flow from the upper and lower hinged plates. Similar to past FIO problems, the current model also exhibits two important lock-in and phase-switch FSI phenomena; however, in angular directions. There is a phase jump of approximately 170° between the aerodynamic lift coefficient and angular oscillations of hinged plates, which nearly occurs in the middle of lock-in region. Indeed, our literature review shows that this is the first time to report the phase-switch phenomenon in angular oscillations of three-element bluff bodies in a FSI problem.     

Refined Formula for Calculating the Hydraulic Drag Coefficient of a Pipeline

It is shown that the hydraulic drag coefficient depends significantly on the ratio of the equivalent hydraulic and geometric roughnesses for developed turbulent pipeline flow in the regime of partial manifestation of roughness. It is proposed to classify roughness types using this ratio. Introducing this parameter into the formula for computing the drag coefficient can considerably extend the domain of applicability of the computational method to different roughness types.     

融合稀疏八叉树与卷积神经网络的汽车风阻系数预测

针对汽车风阻系数预测研究中参数化方法难以准确表征汽车外造型的难题,提出融合稀疏八叉树与卷积神经网络的汽车风阻系数预测方法。将汽车外造型按照八叉树结构离散,使用平均法向量对离散的复杂曲面进行简化,利用卷积神经网络对八叉树形式的汽车外造型进行特征提取,进而对汽车风阻系数进行快速预测。通过改变卷积层数与全连接层数,研究了不同卷积神经网络结构对风阻系数预测精度的影响。与参数化方法相比,本文提出的外造型表示方法能更好地描述模型细节,构建的卷积神经网络结构对风阻系数预测的最小相对误差为1.453%,且计算速度是CFD仿真的1620倍,具有较高的精度及计算效率。     

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

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

Experimental aerodynamic study of a car-type bluff body

The Ahmed body is used as a reference model for fundamental studies of car-type bluff body aerodynamics, in particular focused on the influence of the rear slant angle on the drag coefficient. The objectives of the present work are to obtain reliable drag coefficient comparable to the literature and to explain, based on the nature of the flow, its variation when changing the rear slant angle from 10° to 40°. The drag coefficients measured in both an open and a closed test sections differ by less than 0.5% which proves the reliability and reproducibility of the results. The sensitivity of the drag coefficient to some parameters such as the model roughness or the oncoming boundary layer and the lack of precise information on these parameters in the literature could explain the difference observed with the Ahmed drag coefficient data. The various types of measurement techniques used in the study underline their complementarity. The combination of particle image velocimetry and oil visualization provides a deeper understanding of the flow behaviour around the Ahmed body and a physical interpretation of the drag coefficient evolution.     

带控制片方柱绕流的非定常数值模拟

对高雷诺数下带控制片的方体(在方柱之前放置一小尺度的柱形薄片)非定常绕流进行了数值模拟。数值模拟方法采用RNG重整化群紊流模型，SIMPLE算法，在非定常计算中引入双重时间步方法，将方体和控制片区域作为求解域的一部分作整体求解。结果表明，控制片与方柱之间距离的不同，使得控制片所形成的涡被抑止或产生卡门涡街，从而对柱体侧面边界层的分离发生影响，产生三种典型的流态，柱体时均阻力系数相对于无控制片的情况迅速减小，当控制片偏离柱体轴心时，柱体升力系数迅速增大，偏心至柱体外壁面附近达到极大值。阻力系数达到最小值后将产生跳跃式的增大，而升力系数在达到最大值后也将产生跳跃式的减小。     

Turbulent drag reduction by spanwise wall oscillations

In the present work a technique is numerically investigated, which is aimed at reducing the friction drag in turbulent boundary layers and channel flows. A cyclic spanwise oscillation of the wall with a proper frequency and amplitude is imposed, allowing a reduction of the turbulent drag of up to 40%. The present work is based on the numerical simulation of the Navier-Stokes equations in the simple geometry of a plane channel flow. The frequency of the oscillations is kept fixed at the most efficient value determined in previous studies, while the choice of the best value for the amplitude of the oscillations is evaluated not only in terms of friction reduction, but also by taking into consideration the overall energy balance and the power spent for the motion of the wall. The analysis of turbulence statistics allows to shed some light on the way oscillations interact with wall turbulence, as illustrated by visual inspection of some instantaneous flow fields. Finally, a simple explanation is proposed for this interaction, which leads to a rough estimate of the most efficient value for the frequency of the oscillations.     

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

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

CFD在飞翼标模支撑干扰影响研究中的应用

支撑干扰修正一直是风洞试验的一个重点和难点,特别是对于尾部扁平的飞翼布局飞行器,阻力和俯仰力矩是一个小量,进行风洞试验时,尾部需要局部放大,对阻力和俯仰力矩的影响相对较大,这就使得风洞试验时对支撑影响修正方法提出了更高的要求。本文采用CFD数值模拟手段针对小展弦比飞翼标模开展了网格收敛性和支撑干扰影响研究,分析了尾支撑对飞翼气动特性、压力分布和流场结构的影响,为进一步开展风洞试验的支撑干扰修正提供了依据。在典型工况下,支撑干扰会导致升力系数减小,阻力系数减小,抬头力矩增大;试验模型后体放大和支撑存在会导致飞翼标模上表面压力增大,激波强度变弱,尾部分离区域减小。     

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

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

Optimal control of cylinder wake flow by electro-magnetic force based on adjoint flow field

A nonlinear adjoint-based optimal control approach of cylinder wake flow by using electro-magnetic forcing has been proposed and investigated numerically in the paper. A cost functional representing the balance between the enstrophy Ω² and the interaction parameter N has been developed, and its corresponding adjoint equations have been derived. The sensitivity of the cost functional is found to be a simple function of the adjoint stream function in the adjoint field. Under the action of optimal force, N(t), the flow separation is suppressed successfully, the oscillations of drag and lift disappear and the total drag coefficient decreases dramatically. Furthermore, the global enstrophy is proved to be equal to the total squared strain during the control, which validates the conservation rule of Okubo-Weiss function.     

Slow motion of a slip spheroid along its axis of revolution

A combined analytical and numerical study of the Stokes flow caused by a rigid spheroidal particle translating along its axis of revolution in a viscous fluid is presented. The fluid is allowed to slip at the surface of the particle. The general solution for the stream function in prolate and oblate spheroidal coordinates can be expressed in an infinite-series form of semi-separation of variables. The slip boundary condition incorporating the shear stress at the particle surface is applied to this general solution to determine its unknown coefficients of the leading orders. The solution of these coefficients can be either numerical results obtained from a boundary-collocation method or explicit formulas derived analytically. The drag force exerted on the spheroidal particle by the fluid is evaluated with good convergence behavior for various values of the slip parameter and aspect ratio of the particle. The agreement between our hydrodynamic drag results and the relevant numerical solutions obtained previously using a singularity method is excellent. Although the drag force acting on the translating spheroid normalized by that on a corresponding sphere with equal equatorial radius increases monotonically with an increase in the axial-to-radial aspect ratio for a no-slip spheroid, it decreases monotonically as this aspect ratio increases for a perfect-slip spheroid. The normalized drag force exerted on a spheroid with a given surface slip coefficient in between the no-slip and perfect-slip limits is not a monotonic function of its aspect ratio. For a spheroid with a fixed aspect ratio, its drag force is a monotonically decreasing function of the slip coefficient of the particle.