Computation of hydrodynamic mass and damping coefficients for a cavitating marine propeller flow using a panel method

The present paper deals with the numerical calculation of hydrodynamic mass and damping coefficients under consideration of unsteady sheet cavitation on marine propeller flows. In the first part of the paper, the mathematical and numerical background behind the numerical method is introduced. The numerical calculations carried out in this work are based on a low-order panel method. Panel methods belong to the class of collocation techniques and are applied to obtain a numerical solution of a potential flow based system of boundary integral equations. They are suitable for the present application because of their short computation time which makes them applicable in the design process of marine propellers.Additionally, two different approaches for the determination of hydrodynamic masses and damping are introduced in this work. The hydrodynamic masses and damping are important in studies of the ship motion in seaway and in the analysis of vibrations of a vessel and its appendages. The developed approaches are applied on a propeller flow in heave motion. Hereby, the calculations are performed for a non-rotating and rotating propeller under non-cavitating and cavitating conditions. The results obtained from the simulations are discussed in detail and an outlook is given.     

ANALYSIS OF UNSTEADY CAVITATION FLOW OVER HYDROFOIL BASED ON DYNAMIC MODE DECOMPOSITION$^{\bf 1)}$

空化是船舶和水下航行体推进器中经常发生的一种特殊流动现象,它具有强烈的非定常性,空化的发生往往会影响推进器的水动力性能和效率. 为探究绕水翼非定常空化流场结构,本文基于 Schnerr-Sauer 空化模型和 SST $k$-$\omega $ 湍流模型,开展绕二维水翼非定常空化流动数值预报与流场结构分析. 通过将数值预报的空泡形态演变和压力数据与试验结果对比,验证了建立的数值方法的有效性. 并基于动力学模态分解方法对空化流场的速度场进行模态分解,分析了各个模态的流场特征. 结果表明,第一阶模态对应频率为 0,代表平均流场;第二阶模态对应频率约为空泡脱落频率,揭示了空泡在水翼前缘周期性地生长与脱落,第三阶模态对应频率约为第二阶模态的 2 倍,揭示了两个大尺度旋涡在水翼后方存在融合行为. 第四阶模态对应频率约为第二阶模态的 3 倍,具有更高的频率,表征流场中存在一些小尺度旋涡的融合行为. 最后对不同空化数下的空化流场进行了模态分解分析,发现脱落空泡的旋涡结构随着空化数的减小而增大,第二阶模态频率随着空化数的减小而减小.     

基于动力学模态分解法的绕水翼非定常空化流场演化分析

空化是船舶和水下航行体推进器中经常发生的一种特殊流动现象,它具有强烈的非定常性,空化的发生往往会影响推进器的水动力性能和效率. 为探究绕水翼非定常空化流场结构,本文基于 Schnerr-Sauer 空化模型和 SST $k$-$\omega $ 湍流模型,开展绕二维水翼非定常空化流动数值预报与流场结构分析. 通过将数值预报的空泡形态演变和压力数据与试验结果对比,验证了建立的数值方法的有效性. 并基于动力学模态分解方法对空化流场的速度场进行模态分解,分析了各个模态的流场特征. 结果表明,第一阶模态对应频率为 0,代表平均流场;第二阶模态对应频率约为空泡脱落频率,揭示了空泡在水翼前缘周期性地生长与脱落,第三阶模态对应频率约为第二阶模态的 2 倍,揭示了两个大尺度旋涡在水翼后方存在融合行为. 第四阶模态对应频率约为第二阶模态的 3 倍,具有更高的频率,表征流场中存在一些小尺度旋涡的融合行为. 最后对不同空化数下的空化流场进行了模态分解分析,发现脱落空泡的旋涡结构随着空化数的减小而增大,第二阶模态频率随着空化数的减小而减小.      

Cavity induced vibration of flexible hydrofoils

The objective of this work is to investigate the influence of cavity-induced vibrations on the dynamic response and stability of a NACA66 hydrofoil at 8° angle of attack at Re=750 000 via combined experimental measurements and numerical simulations. The rectangular, cantilevered hydrofoil is assumed to be rigid in the chordwise direction, while the spanwise bending and twisting deformations are represented using a two-degrees-of-freedom structural model. The multiphase flow is modeled with an incompressible, unsteady Reynolds Averaged Navier–Stokes solver with the k–ω Shear Stress Transport (SST) turbulence closure model, while the phase evolutions are modeled with a mass-transport equation based cavitation model. The numerical predictions are compared with experimental measurements across a range of cavitation numbers for a rigid and a flexible hydrofoil with the same undeformed geometries. The results showed that foil flexibility can lead to: (1) focusing – locking – of the frequency content of the vibrations to the nearest sub-harmonics of the foil׳s wetted natural frequencies, and (2) broadening of the frequency content of the vibrations in the unstable cavitation regime, where amplifications are observed in the sub-harmonics of the foil natural frequencies. Cavitation was also observed to cause frequency modulation, as the fluid density, and hence fluid induced (inertial, damping, and disturbing) forces fluctuated with unsteady cavitation.     

Time-dependent hydroelastic analysis of cavitating propulsors

A 3-D potential-based boundary element method (BEM) is coupled with a 3-D finite element method (FEM) for the time-dependent hydroelastic analysis of cavitating propulsors. The BEM is applied to evaluate the moving cavity boundaries and fluctuating pressures, as well as the added mass and hydrodynamic damping matrices. The FEM is applied to analyze the dynamic blade deformations and stresses due to pressure fluctuations and centrifugal forces. The added mass and hydrodynamic damping matrices are superimposed onto the structural mass and damping matrices, respectively, to account for the effect of fluid–structure interaction. The problem is solved in the time-domain using an implicit time integration scheme. An overview of the formulation for both the BEM and FEM is presented, as well as the BEM/FEM coupling algorithm. Numerical and experiment validation studies are shown. The effects of fluid–structure interaction on the propeller performance are discussed.     

Unsteady load on an oscillating Kaplan turbine runner

A Kaplan turbine runner oscillating in turbine waterways is subjected to a varying hydrodynamic load. Numerical simulation of the related unsteady flow is time-consuming and research is very limited. In this study, a simplified method based on unsteady airfoil theory is presented for evaluation of the unsteady load for vibration analyses of the turbine shaft line. The runner is assumed to oscillate as a rigid body in spin and axial heave, and the reaction force is resolved into added masses and dampings. The method is applied on three Kaplan runners at nominal operating conditions. Estimates for added masses and dampings are considered to be of a magnitude significant for shaft line vibration. Moderate variation in the added masses and minor variation in the added dampings is found in the frequency range of interest. Reference results for added masses are derived by solving the boundary value problem for small motions of inviscid fluid using the finite element method. Good correspondence is found in the added mass estimates of the two methods. The unsteady airfoil method is considered accurate enough for design purposes. Experimental results are needed for validation of unsteady load analyses.     

Frequency lock-in during nonlinear vibration of an airfoil coupled with van der Pol Oscillator

The frequency lock-in during the nonlinear vibration of a turbomachinery blade is modeled using a spring-mounted airfoil coupled with a van der Pol Oscillator (VDP) oscillator. The proposed reduced-order model uses the nonlinear VDP oscillator to represent the oscillatory nature of wake dynamics caused by the vortex shedding. The damping term in the VDP oscillator is assumed to be nonlinear. The coupled equations governing the pitch and plunge motion of an airfoil are used to approximate the vibration of a turbomachinery blade. Springs having cubic-order nonlinearity for their stiffnesses are used to mount the airfoil. The unsteady lift acting on the blade is modeled using a self-excited nonlinear wake oscillator. The model for wake dynamics takes into account the influence of blade inertia. The nonlinear coupled three degrees of freedom (dof) aeroelastic system is studied for instability resulting in the frequency lock-in phenomenon. The equations are transformed into non-dimensional form, and then the frequencies of the coupled system are plotted to demonstrate the frequency lock-in. Further, the method of multiple scales is used to derive modulation equations which represent the amplitude and phase of the oscillation. The results obtained using the method of multiple scales are compared with direct numerical solutions to verify the present modeling method. The steady-state amplitudes of the response are plotted against the detuning parameter, which represents the frequency response curve. Further, the sensitivity of non-dimensional parameters such as coupling coefficients, mass ratio, reduced velocity, static unbalance, structural damping coefficient and the ratio of uncoupled pitch and plunge natural frequencies on the frequency response is investigated. The study revealed that parameters such as mass ratio, reduced velocity, structural damping coefficient, and coupling coefficients have a stronger influence in suppressing the amplitude of vibration. Meanwhile, parameters such as the frequency ratio, static unbalance, reduced velocity, and mass ratio significantly affect the range of frequency in which the lock-in phenomenon happens. Further, linear perturbation analysis is done to understand the qualitative effect of the system parameters such as coupling coefficients, mass ratio, frequency ratio, and static unbalance on the range of lock-in.     

风力机叶片非线性挥舞分析

将风力机叶片简化为绕轮毂旋转的变截面Euler-Bernoulli悬臂梁,基于Greenberg公式给出非线性气动力,建立叶片挥舞振动非线性控制方程.由于变截面梁的弯曲刚度和线密度是沿梁轴线变化的函数,无法给出模态函数解析式,论文提出使用假设模态法计算的模态函数,作为基函数对控制方程进行Galerkin截断,通过将挥舞振动分解为静态位移和动态扰动合成,对其进行动态响应分析,同时讨论了叶轮转速、风速和旋转位置对振动特性的影响.研究表明:(1)叶轮转速对叶片挥舞特性影响显著,风速和叶片转角对振动特性影响很小.(2)静态位移随风速增加而增大,大体上成线性关系,气动阻尼随风速增加而减小.(3)风速较低时,非线性挥舞振动表现为衰减振动,随着风速增加,振动由衰减振动演化为周期运动,再由周期运动演化为拟周期运动.     

Dynamics of cavitation–structure interaction

Cavitation–structure interaction has become one of the major issues for most engineering applications. The present work reviews recent progress made toward developing experimental and numerical investigation for unsteady turbulent cavitating flow and cavitation–structure interaction. The goal of our overall efforts is to(1) summarize the progress made in the experimental and numerical modeling and approaches for unsteady cavitating flow and cavitation–structure interaction,(2) discuss the global multiphase structures for different cavitation regimes, with special emphasis on the unsteady development of cloud cavitation and corresponding cavitating flow-induced vibrations,with a high-speed visualization system and a structural vibration measurement system, as well as a simultaneous sampling system,(3) improve the understanding of the hydroelastic response in cavitating flows via combined physical and numerical analysis, with particular emphasis on the interaction between unsteady cavitation development and structural deformations. Issues including unsteady cavitating flow structures and cavitation–structure interaction mechanism are discussed.     

绕水翼空化流动及振动特性的实验研究

空化是发生在水力机械内部的一种水动力现象,其发展具有显著的非定常特性.空化流动中空穴的脱落以及溃灭会诱发结构振动,对水力机械的效率、噪声、安全性等造成影响. 研究空化流动中结构的振动特性具有重要的工程意义. 采用实验的方法研究了绕NACA66 水翼空化流动的空穴形态和水翼振动特性. 实验在一闭式空化水洞中进行. 采用高速摄像技术观测不同空化阶段的空穴形态,应用多普勒激光测振仪测量水翼的振动速度,并通过一套同步系统实现了高速相机和多普勒激光测振仪的同步触发和测量. 采用小波分析方法对不同空化阶段下的空穴形态和水翼振动数据在时域和频域中的特性进行了分析.对云状空化阶段的同步测量结果进行了研究,分析了振动与空穴发展过程的联系. 结果表明,随着空化数的降低,流场经历了无空化、初生空化、片状空化和云状空化4个阶段,水翼的振动强度呈逐渐增大趋势. 在片状空化和云状空化阶段,空穴脱落导致水翼振动,诱发的振动频率与空穴脱落频率相同. 对于云状空化,在附着型空穴生长阶段水翼发生高频小幅度振动,在空穴脉动和断裂脱落期间水翼表现为低频大幅振动.      

Analysis of frequency-domain vibration response of thin plate attached with viscoelastic free layer damping

To predict the vibration response of viscoelastic composite structure, two key issues need to be conducted, one is introducing the constitutive model of viscoelastic material into the analysis model and the other is describing the real damping behavior of viscoelastic composite structure. The emphasis of this study is to obtain the effects of frequency dependence on the vibration response of viscoelastic composite structure and the method of introducing two kinds of damping (viscoelastic material damping and remaining equivalent viscous damping). Vibration response analysis in frequency domain was investigated for viscoelastically damped plate. A cantilever plate attached with the ZN_1 viscoelastic free layer damping (FLD) was chosen to demonstrate the developed method. Frequency-domain response of the composite plate were solved and the obtained results were compared with the experimental values for the purpose of assessing the rationality of the proposed method. In addition, in order to obtain the effects of viscoelastic material parameters on vibration response of viscoelastic composite structure, a detailed parametric analysis was performed. This study shows that the frequency dependent characteristic of viscoelastic material has significant influence on the vibration response in the resonant region and acceptable results can be achieved in the non-resonant region if frequency dependent parameters are substituted by average values of the viscoelastic parameters reasonably in the analysis process.     

某型电动飞机螺旋桨振动特性实验研究

某型电动飞机采用螺旋桨产生拉力,为了防止螺旋桨工作时共振,利用ES-2-150振动实验系统进行了两叶木质螺旋桨和碳纤维螺旋桨的振动特性实验,采用谐振搜索与驻留方法测量出木质螺旋桨的第一阶固有频率为36.07Hz,碳纤维螺旋桨的第一阶固有频率为73.58Hz。螺旋桨爬升状态转频为39Hz,这与木质螺旋桨的第一阶固有频率非常接近,导致木质螺旋桨在爬升状态出现比较严重的振动故障。因此,某型电动飞机最终选择两叶碳纤维螺旋桨作为其拉力产生装置。     

Zur Theorie des Trochoidenpropellers

Übersicht In der vorliegenden Arbeit wird eine Theorie zur Behandlung der Strömung durch Trochoidenpropeller entwickelt. Sowohl die Flügelzirkulation und die Flügelkräfte als auch die Propellerkräfte und Wirkungsgrade werden berechnet. Ein Vergleich der theoretischen Ergebnisse mit Experimenten von van Manen zeigt befriedigende Übereinstimmung.

---

Summary In this paper a theoretical hydrodynamic investigation of trochoidal vertical axis propellers is made. The blade circulation and the blade forces as well as propeller forces and propeller efficiency are calculated by using unsteady lifting line theory. The theoretical results are compared with measurements by van Manen; the agreement is shown to be satisfactory.

     

Effect of trailing edge shape on hydrodynamic damping for a hydrofoil

Flow induced vibration on a hydrofoil may be significantly reduced with a slight modification of the trailing edge without alteration of the hydrodynamic performance. Particularly, the so called Donaldson trailing edge shape gave remarkable results and is being used in a variety of industrial applications. Nevertheless, the physics behind vibration reduction is still not understood. In the present study, we have investigated the hydrodynamic damping of a 2D hydrofoil with Donaldson trailing edge shape. The results are compared with the same hydrofoil with blunt trailing edge. The tests are carried out in EPFL high speed cavitation tunnel and two piezoelectric patches are used for the hydrofoil excitation in non-intrusive way. It was found that the hydrodynamic damping is significantly increased with the Donaldson cut. Besides, as the flow velocity is increased, the hydrodynamic damping is found to remain almost constant up to the hydrofoil resonance and then increases linearly, for both tested trailing edge shapes and for both first bending and torsion modes.     

旋转射弹高速倾斜入水多相流场与弹道数值模拟

基于VOF多相流模型和有限体积法求解水、汽、气多相流动的RANS方程,结合重叠网格技术和six DOF算法对某一型号舰载射弹倾斜入水过程进行数值模拟研究。首先基于该方法研究了射弹旋转效应对射弹运动特性及流体动力特性的影响,然后对不同入水角下倾斜入水过程进行分析,得到不同倾角下旋转射弹入水空泡形态发展规律、弹体运动特征及流体动力特性变化规律。研究结果表明:射弹的旋转有利于弹体在初始对称面内的弹道稳定性,但会降低弹体侧向稳定性,使射弹受到的阻力系数、俯仰力矩系数变小;入水角越小,形成的空泡越不对称,由射弹运动状态的改变引起的空泡形态变化越明显,在超空泡航行阶段,弹体运动较稳定,不同角度下流体动力系数差别很小,当弹体下表面刺破空泡壁沾湿时,弹体运动状态发生较大变化,流体动力系数迅速增大,此时入水角度过小,弹体容易失稳;弹体的沾湿对空泡形态、弹体运动稳定性和流体动力特性有着重要的影响。     

DYNAMIC CHARACTERISTICS OF LAG VIBRATION OF A WIND TURBINE BLADE

Lagwise dynamic characteristics of a wind turbine blade subjected to unsteady aero- dynamic loads are studied in this paper. The partial differential equations governing the coupled longitudinal-transverse vibration of the blade with large bending deflection are obtained by ap- plying Hamilton＇s principle. The modal problem of the coupled vibration is handled by using the method of numerical integration of Green＇s function. Influences of the rotating speed, the pitch angle, the setting angle, and the aerodynamic loads on natural frequencies are discussed. Results show that： （I） Lagwise natural frequencies ascend with the increase of rotating speed; effects of the rotating speed on low-frequencies are dramatic while these effects on high-frequencies become less. （2） Influences of the pitch angle on natural frequencies are little; in the range of the normal rotating speed, the first frequency ascends with the increase of the absolute value of the pitch angle, while it is contrary to the second and third frequencies. （3） Effects of the setting angle on natural frequencies depend on the rotating speed; influences are not significant at low speed, while they are dramatic on the first frequency at high speed. （4） Effects of the aerodynamic loads on natural frequencies are very little; frequencies derived from the model considering aerodynamic loads are smaller than those from the model neglecting aerodynamic loads; relative errors of the results corresponding to two models ascend with the increase of the absolute value of the setting angle.     

Fluid–structure coupling for an oscillating hydrofoil

Fluid-structure investigations in hydraulic machines using coupled simulations are particularly time-consuming. In this study, an alternative method is presented that linearizes the hydrodynamic load of a rigid, oscillating hydrofoil. The hydrofoil, which is surrounded by incompressible, turbulent flow, is modeled with forced and free pitching motions, where the mean incidence angle is 0° with a maximum angle amplitude of 2°. Unsteady simulations of the flow, performed with ANSYS CFX, are presented and validated with experiments which were carried out in the EPFL High-Speed Cavitation Tunnel. First, forced motion is investigated for reduced frequencies ranging from 0.02 to 100. The hydrodynamic load is modeled as a simple combination of inertia, damping and stiffness effects. As expected, the potential flow analysis showed the added moment of inertia is constant, while the fluid damping and the fluid stiffness coefficients depend on the reduced frequency of the oscillation motion. Behavioral patterns were observed and two cases were identified depending on if vortices did or did not develop in the hydrofoil wake. Using the coefficients identified in the forced motion case, the time history of the profile incidence is then predicted analytically for the free motion case and excellent agreement is found for the results from coupled fluid–structure simulations. The model is validated and may be extended to more complex cases, such as blade grids in hydraulic machinery.     

Numerical analysis of unsteady tip leakage vortex cavitation cloud and unstable suction-side-perpendicular cavitating vortices in an axial flow pump

The objective of this work is to simulate and analyze the formations of three-dimensional tip leakage vortex (TLV) cavitation cloud and the periodic collapse of TLV-induced suction-side-perpendicular cavitating vortice (SSPCV). Firstly, the improved SST k–ω turbulence model and the homogeneous cavitation model were validated by comparing the simulation result with the experiment of unsteady cavitation shedding flow around the NACA66-mod hydrofoil, and then the unsteady TLV cloud cavitation and unstable SSPCV in an axial flow pump were predicted using the improved numerical method. The predicted three-dimensional cavitation structures of TLV and SSPCV as well as the collapsing features show a good qualitative agreement with the high speed photography results. Numerical results show that the TLV cavitation cloud in the axial flow pump mainly includes tip clearance cavitation, shear layer cavitation, and TLV cavitation. The unsteady TLV cavitation cloud occurs near the blade trailing edge (TE) where the shapes of sheet cavitation and TLV cavitation fluctuate. The inception of SSPCV is attributed to the tail of the shedding cavitation cloud originally attached on the suction side (SS) surface of blade, and the entrainment affect of the TLV and the influence of the tip leakage flow at the tailing edge contribute to the orientation and development of the SSPCV. The existence of SSPCV was evidently approved to be a universal phenomenon in axial flow pumps. At the part-load flow rate condition, the SSPCV may trigger cavitation instability and suppress the tip cavitation in the neighboring blade. The cavitation cloud on the SS surface of the neighboring blade grows massively, accompanying with a new SSPCV in the neighboring flow passage, and this SSPCV collapses in a relatively short time.     

Control of the corner separation in a compressor cascade by steady and unsteady plasma aerodynamic actuation

This paper reports experimental results on using steady and unsteady plasma aerodynamic actuation to control the corner separation, which forms over the suction surface and end wall corner of a compressor cascade blade passage. Total pressure recovery coefficient distribution was adopted to evaluate the corner separation. Corner separation causes significant total pressure loss even when the angle of attack is 0°. Both steady and unsteady plasma aerodynamic actuations suppress the corner separation effectively. The control effect obtained by the electrode pair at 25% chord length is as effective as that obtained by all four electrode pairs. Increasing the applied voltage improves the control effect while it augments the power requirement. Increasing the Reynolds number or the angle of attack makes the corner separation more difficult to control. The unsteady actuation is much more effective and requires less power due to the coupling between the unsteady actuation and the separated flow. Duty cycle and excitation frequency are key parameters in unsteady plasma flow control. There are thresholds in both the duty cycle and the excitation frequency, above which the control effect saturates. The maximum relative reduction in total pressure loss coefficient achieved is up to 28% at 70% blade span. The obvious difference between steady and unsteady actuation may be that wall jet governs the flow control effect of steady actuation, while much more vortex induced by unsteady actuation is the reason for better control effect.     

Analysis of unstable vortical structure in a propeller wake affected by a simulated hull wake

A two-frame PIV (particle image velocimetry) technique was used to investigate the flow characteristics of a complicated propeller wake influenced by a hull wake. As the propeller is significantly affected by the hull wake of a marine vessel, measurements of the propeller wake under the hull wake are certainly needed for more reliable validation of numerical predictions. Velocity field measurements were conducted in a cavitation tunnel with a simulated hull wake. Generally, the hull wake generated by the hull of a marine ship may cause different loading distributions on the propeller blade in both the upper and the lower propeller planes. The unstable propeller wake caused by the ship’s hull was interpreted in terms of turbulent kinetic energy (T _KE) to obtain useful information for flow modeling. The unstable or unsteady phenomenon in the upper propeller wake was identified by using the proper orthogonal decomposition (POD) method to characterize the coherent flow structure with turbulent kinetic energy. Strong unsteadiness appeared in the second and higher modes, largely affecting the downstream flow characteristics. The first eigenmode can be used to appropriately identify the tip vortex positions even in the unstable downstream region, which are helpful for establishing reliable wake modeling.