Experimental investigation of a hydrofoil designed via hydrostructural optimization

In the last decade, there has been an increased interest in the use of multidisciplinary optimization techniques for the design of aerospace, maritime, and wind engineering systems. However, validation of numerically optimized results using experimental measurements has been scarce. In this paper, numerical predictions are compared with experimental measurements of the hydrodynamic forces, deformations, and cavitation performance for a baseline NACA 0009 hydrofoil and an optimized hydrofoil. Both hydrofoils are made of solid aluminum, and are cantilevered at the root. One of the hydrofoils is optimized using a high-fidelity hydrostructural solver combined with a gradient-based optimizer, as detailed by Garg et al. (2017). The numerical predictions agree well with experimental measurements for both the baseline NACA 0009 and the optimized hydrofoils. For the optimized hydrofoil, the mean differences between the predicted and measured values for mean lift, drag coefficient, and moment coefficients, are 2.9%, 5.1%, and 3.0%, respectively. For the non-dimensional tip bending deflection, the mean difference is 3.4%. Although the optimized hydrofoil is significantly thicker to withstand higher loads than the baseline, it yields an overall measured increase in the lift-to-drag ratio of 29% for lift coefficients ranging from $-0.15$ to 0.75 and exhibits significantly delayed cavitation inception compared to the baseline. The improvement in hydroelastic and cavitation performance is attributed to the changes in the distribution of camber, twist, thickness, and the leading edge radius of the optimized hydrofoil. The results validate the analysis and optimization of the high-fidelity hydrostructural design optimization approach, and opens up new possibilities for the design of high-performance hydrofoils, marine propellers, and turbines.     

An improved progressive preconditioning method for steady non‐cavitating and sheet‐cavitating flows

An improved progressive preconditioning method for analyzing steady inviscid and laminar flows around fully wetted and sheet‐cavitating hydrofoils is presented. The preconditioning matrix is adapted automatically from the pressure and/or velocity flow‐field by a power‐law relation. The cavitating calculations are based on a single fluid approach. In this approach, the liquid/vapour mixture is treated as a homogeneous fluid whose density is controlled by a barotropic state law. This physical model is integrated with a numerical resolution derived from the cell‐centered Jameson's finite volume algorithm. The stabilization is achieved via the second‐and fourth‐order artificial dissipation scheme. Explicit four‐step Runge–Kutta time integration is applied to achieve the steady‐state condition. Results presented in the paper focus on the pressure distribution on hydrofoils wall, velocity profiles, lift and drag forces, length of sheet cavitation, and effect of the power‐law preconditioning method on convergence speed. The results show satisfactory agreement with numerical and experimental works of others. The scheme has a progressive effect on the convergence speed. The results indicate that using the power‐law preconditioner improves the convergence rate, significantly. Copyright © 2010 John Wiley & Sons, Ltd.     

圆柱绕流的电磁控制

流体绕过非流线物体时，在物体尾部形成涡街，使其表面周期性变化的阻力和升力增加，从而导致物体振荡，产生噪音。本文通过实验和计算，研究圆柱绕流的电磁控制，阐述浸于弱电介质溶液中，表面包覆电磁激活板的圆柱，在电磁力作用下的流体控制原理，讨论电磁力的消涡、减阻和减振过程。     

Frequency effects on lift and drag for flow past an oscillating cylinder

A transversely oscillating cylinder in a uniform flow is modeled to investigate frequency effects of flow-induced wake on lift and drag of the cylinder. Specifically, verified unsteady fluid dynamic simulations using an immersed-boundary method in a fixed Cartesian grid predict the flow structure around the cylinder and reveal how the integration of surface pressure and shear distributions provides lift and drag on the oscillating cylinder. In this study, frequency ranges to be considered are both near and away from the natural frequency of wake vortex shedding. Subsequently, the effects of frequency lock-in, superposition and demultiplication on lift and drag are discussed based on the spectral analysis of time histories of lift and drag.     

Influences of Lorentz force on the hydrofoil lift

In this paper, Lorentz forces are proved to be able to suppress separation in flows over hydrofoils. Furthermore, a differential equation of pressure distributions on the hydrofoil surface is derived, from which it is found that BVF （boundary vortex flux） σ is a suitable criterion for describing the lift coefficient variations during the electromagnetic control process. According to our numerical results, the periodic variations of lift for a hydrofoil at an attack angle of 17 ° are analyzed and its inherent mechanism is discussed in detail with the concept of BVE On the other hand, the effects of Lorentz force on the hydrofoil＇s lift are investigated both experimentally and numerically for different magnitudes and locations.     

A NUMERICAL STUDY OF OPTIMAL DRAG REDUCTION FOR TURBULENT TRANSONIC PROJECTILES USING A PASSIVE CONTROL

Abstract

The purpose of this research is to numerically study a drag reduction method—passive control of shock/boundary layer interaction, which is applied to the boattail portion of a secant-ogive-cylinder-boattail projectile in turbulent transonic flows. The flow pattern and the components of aerodynamic drag computed from numerical data are analyzed. The effectiveness of this method is studied by varying the values of parameters such as porosity distribution, maximum porosity factor and size of porous region. The conditions for optimal drag reduction are investigated and reported. The present results show that the use of this passive control method can not only reduce the boattail drag but also the base drag, and results in an additional 8% total drag reduction compared to that without the passive control technique. This passive control method can be an effective approach for the design of high-performance projectiles in the transonic regime.     

Reduced-order modeling of three-dimensional unsteady partial cavity flows

An efficient reduced-order modeling to analyze three-dimensional unsteady partial cavity flows is proposed. The proposed approach is based on the boundary element method along with the potential flow assumption. To this end, a novel non-iterative method based on the flow eigenmodes of three-dimensional partial cavity flows is applied. Eigenanalysis and reduced-order modeling for unsteady flows over a three-dimensional hydrofoil with various sections are performed. The results obtained from the present analysis are compared with those reported in the literature to verify the strength of the proposed approach. In order to examine the performance of the introduced algorithm for unsteady cavitating flows, various simulations for several reduced frequencies, hydrofoil geometries and different cavitation numbers are also investigated. Comparison between the obtained results using the novel and conventional methods indicates that the present algorithm works very well with sufficient accuracy. Moreover, it is shown that the proposed method is computationally more efficient than the conventional ones for unsteady sheet cavitation analysis on three-dimensional hydrofoils.     

On the origin of drag increase in varying-phase opposition control

This study investigates the physics underlying the drag increase in a low Reynolds number turbulent channel flow due to varying-phase opposition control by means of direct numerical simulation and modal analysis. The drag increase occurs for an extended region of the parameter space and we consider a controller with a positive phase shift in Fourier domain between sensor measurement and actuator response as a representative example for this regime. Analyses of instantaneous flow fields as well as spatial power spectra show that the structure of drag-increased flows is remarkably different from that of drag-reduced and canonical flows. In particular, the near-wall region is dominated by structures of short streamwise and large spanwise extent. Isolation of a representative control scale shows that these energetic structures can be characterized as spanwise rollers, which induce strong ejection and sweep motions and lead to drag increase. The presence of rollers, and therefore drag increase, in the full nonlinear system correlates well with the presence of an amplified eigenvalue in the eigenspectrum of the linearized Navier–Stokes operator. It is further shown that the scales responsible for drag increase at positive phase shifts are inactive at negative phase shifts and do not contribute to drag reduction. These scales can therefore be excluded from a controller aimed at drag reduction, which relaxes the spatial resolution requirements on the control hardware. The eigenspectrum may be used as a computationally cheap tool to identify such detrimental scales during an early design stage.     

高超声速钝前缘乘波构型优化设计研究

乘波体的高升阻比优势使其在高超声速飞行器设计中极具应用前景. 在实际工程应用中, 为了满足防热要求, 乘波体前缘必须进行钝化处理, 前缘钝化对乘波体气动性能会产生显著影响. 因此, 原始尖前缘最优乘波体并不能保证钝化后仍为最优. 针对这一问题, 首先研究了前缘钝化对不同构型升阻特性的影响程度和作用机理. 结果表明: 前缘钝化会造成乘波体升力小幅度降低, 阻力大幅增加, 升阻比显著降低; 其中钝前缘本身的波阻在阻力增量中起主导作用, 而钝前缘本身的摩阻增加量与物面的摩阻降低量十分接近. 基于上述结果, 提出了一种高效评估钝前缘乘波体气动力的方法, 并结合遗传算法, 开展了直接考虑前缘钝化影响的乘波构型优化设计研究, 获得了钝前缘最优构型. 通过CFD数值模拟对最优构型的气动力特性进行评估, 结果表明: 在不同飞行高度、不同升力和不同钝化半径约束下, 相比尖前缘最优构型, 钝前缘最优构型宽度变窄, 相同纵向位置处的后掠角增大, 且升阻比显著提升. 在M_∞ = 15, H = 50 km, CL = 0.3约束条件下, 钝化半径R = 10 mm的钝前缘最优构型设计点升阻比相比尖前缘最优构型提升量可达9.32%.      

Determination of the lift of a partially cavitated hydrofoil

New methods for the velocity circulation around partially cavitated hydrofoils are proposed for both short and long cavities. The calculated results are compared with available experimental data. The calculated lift and cavitation number dependences remain bounded as the cavity end approaches the trailing edge of the hydrofoil.St. Petersburg. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 71–74, November–December, 1994.     

Influence of cavitation on the hydroelastic stability of hydrofoils

This work numerically examines the effect of turbulent and cavitating flow on the hydroelastic response and stability of a hydrofoil. A cantilevered, rectangular, chordwise rigid hydrofoil is modeled as a 2-degrees-of-freedom structure for its spanwise bending and torsional flexibilities. The fluid flow is modeled with the incompressible, Unsteady Reynolds Averaged Navier–Stokes equations using an eddy-viscosity turbulence closure model that is corrected for the presence of cavitation, and with a transport equation based cavitation model. The results show that, in general, massive cavitation tends to: (i) reduce the mean lift, (ii) increase the mean drag, (iii) lower the mean deformations, and (iv) delay static divergence, while unsteady sheet/cloud cavitation promotes flow induced vibrations. Such vibrations and load fluctuations could be as large as (and even greater than) the mean values for cases with unsteady cavitation, so dynamic and viscous fluid–structure models are needed to simulate flexible hydrofoils in cavitating flows. In general, the flow induced vibrations, and hence the drag force, are higher with decreasing stiffness. For small leading edge partial cavitation, increasing foil flexibility increases the maximum cavity length and reduces the cavity shedding frequency; however, the influence of foil flexibility is limited for cases where the maximum cavity length is near or beyond the foil trailing edge, because of the relocation of the center of pressure at the elastic axis, near the mid-chord. The results show that the mean deformations are generally limited by stall, and by the quasi-steady linear theory predictions at the fully-wetted and supercavitating limits. Furthermore, frequency focusing can occur when the cavity shedding frequency is near the fundamental system resonance frequencies, and broadening of the frequency spectrum can occur due to excitation of the sub-harmonics and/or modulation induced by the fluctuating cavities, if the cavity shedding frequency is away from the fundamental system resonance frequencies.     

Reynolds-averaged modeling of polymer drag reduction in turbulent flows

A novel eddy viscosity model for predicting friction drag reduction induced by polymers in turbulent wall-bounded flows is presented. The approach is based on the elliptic relaxation model modified to account for the modified Reynolds-stress equilibrium established by the presence of elastic polymer chains in the fluid. The increased wall damping of the turbulent fluctuations is obtained by modifying the pressure–strain redistribution term. Polymer solutions are represented using the Finite Extensibility Non-linear elastic FENE-P dumbbell model; only one transport equation for the elongation of the polymer chains is considered. The model reproduces the level of drag reduction observed over a wide range of rheological parameters. In addition, both the mean velocity and the turbulent fluctuations are predicted with good accuracy. The approach is computationally attractive because of its limited increase in computational cost in comparison with its Newtonian counterpart.     

定常吹气对倾转旋翼机向下载荷的影响

为进一步提高倾转旋翼机悬停状态下的有效载重,开展了定常吹气流动控制对向下载荷的影响研究。首先应用延迟脱体涡模拟(DDES)方法对翼型-90°迎角下非定常大范围分离流动结构进行了数值分析;然后分别开展了前缘吹气、后缘吹气降载措施研究,揭示了吹气降载的机理,并对不同吹气口位置和吹气动量系数的影响进行了定量分析,最后开展了前、后缘同时吹气作用下降载数值模拟研究。计算结果表明:前缘最佳吹气位置在翼型的前缘点,而后缘吹气最佳位置位于襟翼弦长的15%处;前缘吹气的降载效果要优于后缘吹气,而且吹气动量系数对向下载荷的影响较小;相对于初始未施加流动控制构型,阻力系数减小量可达到32.72%。     

圆柱附加旋转整流罩表面压力分布及气动力特性

对圆柱附加固定整流罩的已有研究表明,它在降低升阻力和抑制涡激振动方面有优良的效果。但固定整流罩具有方向敏感性,当来流方向改变后效果会受到显著影响,甚至起到增加升阻力和加剧涡激振动的反作用。本文给圆柱附加了圆弧直径为40mm,形状夹角α分别为30°、45°、60°、75°和90°五种尺寸的旋转整流罩,并进行了风洞实验。其中整流罩可以自由地围绕圆柱轴线旋转。实验结果表明:旋转整流罩在流体力产生的力矩作用下,旋转至一个偏离尾流中心线固定角度的动态平衡位置,而平衡位置偏转角δ随着形状夹角α的增大而增大。附加旋转整流罩后,相对单圆柱能够提高尾迹区域压力,并能使时均阻力和脉动升力分别在α=30°和α=75°时获得最大43.5%和67.0%的降低。此外,对于小α(α≤60°)情况,漩涡脱落频率明显高于单圆柱情况,而对于大α(α≥75°)情况,则与单圆柱情况相接近。所有旋转整流罩升力主频的幅值较之单圆柱有了很大程度的降低,可见旋转整流罩在抑制漩涡脱落方面有很好的效果。     

Turbulent Duct Flow Controlled with Spanwise Wall Oscillations

The spanwise oscillation of channel walls is known to substantially reduce the skin-friction drag in turbulent channel flows. In order to understand the limitations of this flow control approach when applied in ducts, direct numerical simulations of controlled turbulent duct flows with an aspect ratio of A R = 3 are performed. In contrast to channel flows, the spanwise extension of the duct is limited. Therefore, the spanwise wall oscillation either directly interacts with the duct side walls or its spatial extent is limited to a certain region of the duct. The present results show that this spanwise limitation of the oscillating region strongly diminishes the drag reduction potential of the control technique. We propose a simple model that allows estimating the achievable drag reduction rates in duct flows as a function of the width of the duct and the spanwise extent of the controlled region.     

Investigation of vortical flow over bluff bodies with base cavities

Based on our previous research about drag reduction in term of the base cavity length using two dimensional simulations, this paper describes a numerical study of a bluff body of which the number of base cavities is successively increased and the cavity geometries are also modified to assume different shapes. Here we attempt to find an effective configuration to reduce the drag by increasing the number of base cavities. The numerical simulations examining varied number of base cavities reveal the presence of different strength of vortices in the wake zone which is the reason why the drag coefficients are distinctly different for different cases. In the case with double and triple rectangular cavities, we use the pressure contours snapshots at successive time instants to describe the wake evolution. We further investigate the effect of variable base cavity shapes for a constant cavity length at an identical time instant. A total of two different geometries of base cavities are discussed here: the rectangular and the sinusoidal cavities with sharp and rounded trailing edges, respectively. The numerical results reveal that the former is an effective drag reduction configuration which can produce a significant base pressure recovery corresponding to the strength of the vortices shown in the pressure contour figures. While the latter shows no obvious reduction in drag coefficient and a similar intensity of vortex in the wake zone compared with the unmodified case. Reductions in drag are observed for all the investigated cavity configurations, and additionally it is found that the magnitude of the reduction bears a direct relationship with the number of the cavities up to a certain minimum value.     

高升阻比乘波构型优化设计

在M∞ =6, 30km高空条件下,以升阻比为目标函数,进行了锥形流乘波体的黏性优化设计,讨论 了影响乘波体升阻比的因素,并对优化结果进行了数值验证. 结果表明：对于升阻比最大的 黏性优化乘波体,存在最优圆锥角使得源自该基本流场的乘波体升阻比最大;摩阻和波阻处 于同一量级;体积率、细长比和展长比都随着基本流场圆锥角的增大而增大.     

Migration of settling particles in a horizontal viscous flow through a vertical slot with porous walls

Particle migration in a horizontal flow of dilute suspension through a vertical slot with porous walls is studied using the two-continua approach. The lateral migration is induced by two opposite effects: an inertial lift force due to particle settling and directed toward the slot centre-line, and a drag due to leak-off entraining particles toward the walls. An expression for the inertial lift on a settling particle in a horizontal channel flow found recently is generalized to the case of a low leak-off velocity. The evolution of an initial uniform particle concentration profile is studied within the full Lagrangian approach. Four migration regimes are found differing by the direction of particle migration and numbers of equilibrium positions. Conditions of the regime change and a critical value of dimensionless leak-off velocity for particle deposition on the walls are obtained analytically. Suspension flows with zones where the particle concentration is zero or increases infinitely, are studied numerically.     

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

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

电磁极宽度对圆柱绕流场影响的数值分析

在雷诺数Re=200的情况,利用Maxwell方程直接数值计算表面包覆电极与磁极圆柱体产生的电磁力分布,将其加入到动量方程中,在各种电磁力作用参数和电磁极宽度的组合下,对表面覆盖电磁极圆柱体在弱电解质中的绕流场结构及其升阻力特性进行了数值模拟与分析.结果表明,当电磁极宽度较小时,圆柱体绕流场的分离点越容易接近后驻点,而电磁力对总阻力的影响并不明显,但对压差和摩擦阻力均有明显影响.当电磁极宽度较大时,圆柱体尾部区域越容易产生射流现象,而且总阻力随电磁力作用参数和电磁极宽度增大而减小.在电磁力尚不足以完全抑制周期性涡脱落的情况下,升力幅值随电磁力作用参数增大而减小,但随电磁极宽度则先减小后略有增加,升力脉动频率则均随电磁力作用参数和电磁极宽度增大而增加.研究表明,电磁力可以有效地改善圆柱体绕流场结构,达到减小圆柱体阻力并抑制其脉动升力之目的,因此是圆柱型结构的一种有效流动控制手段.