垂荡双气室振荡水柱波能装置水动力特性研究

振荡水柱(OWC)波能转换装置因其结构简单、便于安装维护等特点, 被公认为最具应用前景的波能转换技术. 本研究以垂荡式双气室OWC波能转换装置为研究对象, 借助开源代码平台OpenFOAM及基于interFoam求解器开发的造/消波工具箱waves2Foam, 采用流体体积法(VOF)捕捉自由面和六自由度(6DOF)动网格求解器模拟垂荡运动响应, 数值研究在不同入射规则波作用下, 前后气室相对宽度、弹簧弹性系数对装置捕能宽度比及水动力特性的影响规律. 通过与已有的固定情况下的双气室OWC装置结果对比, 并通过对比自由衰减运动响应验证动网格技术, 揭示了本研究中数值模型的合理性和有效性. 计算结果表明, 较宽的后气室结构布置有利于双气室振荡水柱装置在垂荡状态下的波能提取; 前后气室宽度比为1/2时, 垂荡式双气室OWC装置在测试波频段具有最优的捕能宽度比; 相较于固定状态, 垂荡装置的后气室在中高波频段有着更高的捕能宽度比; 装置前后气室内水柱与OWC装置垂荡运动间存在的相位差使得气室内水面相对振幅和相对压强在测试波频段存在多峰值现象, 进一步发现弹装置通过垂向弹簧进行相位控制, 可显著拓宽高效频谱带, 实现较大的捕能宽度比.      

Hardware-in-the-loop simulation of wave energy converters based on dielectric elastomer generators

Dielectric elastomer generators (DEGs) are soft electrostatic generators based on low-cost electroactive polymer materials. These devices have attracted the attention of the marine energy community as a promising solution to implement economically viable wave energy converters (WECs). This paper introduces a hardware-in-the-loop (HIL) simulation framework for a class of WECs that combines the concept of the oscillating water columns (OWCs) with the DEGs. The proposed HIL system replicates in a laboratory environment the realistic operating conditions of an OWC/DEG plant, while drastically reducing the experimental burden compared to wave tank or sea tests. The HIL simulator is driven by a closed-loop real-time hydrodynamic model that is based on a novel coupling criterion which allows rendering a realistic dynamic response for a diversity of scenarios, including large scale DEG plants, whose dimensions and topologies are largely different from those available in the HIL setup. A case study is also introduced, which simulates the application of DEGs on an OWC plant installed in a mild real sea laboratory test-site. Comparisons with available real sea-test data demonstrated the ability of the HIL setup to effectively replicate a realistic operating scenario. The insights gathered on the promising performance of the analysed OWC/DEG systems pave the way to pursue further sea trials in the future.

     

振动翼型和襟翼的绕流特性研究

用数值模拟手段详细地研究了振动翼型和襟翼的绕流问题,数值模拟的出发方程为Euler和N-S方程,格式为Bcam-Warming格式的改进型。数值实验主要针对流场的二大特性进行的,即振动对激波的影响和振动对分离的抑制作用,结果表明:(1)随翼型或襟翼的振动激波强度和位置也相应地变化但这一变化滞后于攻角的变化;(2)振幅加大激波强度的变化和激波运动范围也加大;(3)振动频率越高对激波的影响反而较低频时要小;(4)流动条件的不同可使升力回线的走向发生变化;(5)振动对分离有明显的抑制作用。     

The study of sound wave propagation in rarefied gases using unified gas-kinetic scheme

Sound wave propagation in rarefied monatomic gases is simulated using a newly developed unified gaskinetic scheme (UGKS). The numerical calculations are carried out for a wide range of wave oscillating frequencies. The corresponding rarefaction parameter is defined as the ratio of sound wave frequency to the intermolecular particle collision frequency. The simulation covers the flow regime from the continuum to free molecule one. The treatment of the oscillating wall boundary condition and the methods for evaluating the absorption coefficient and sound wave speed are presented in detail. The simulation results from the UGKS are compared to the Navier-Stokes solutions, the direct simulation Monte Carlo (DSMC) simulation, and experimental measurements. Good agreement with the experimental data has been obtained in the whole flow regimes for the corresponding Knudsen number from 0.08 to 32. The current study clearly demonstrates the capability of the UGKS method in capturing the sound wave propagation and its usefulness for the rarefied flow study.     

Cellular-level near-wall unsteadiness of high-hematocrit erythrocyte flow using confocal μPIV

In hemodynamics, the inherent intermittency of two-phase cellular-level flow has received little attention. Unsteadiness is reported and quantified for the first time in the literature using a combination of fluorescent dye labeling, time-resolved scanning confocal microscopy, and micro-particle image velocimetry (μPIV). The near-wall red blood cell (RBC) motion of physiologic high-hematocrit blood in a rectangular microchannel was investigated under pressure-driven flow. Intermittent flow was associated with (1) the stretching of RBCs as they passed through RBC clusters with twisting motions; (2) external flow through local obstacles; and (3) transitionary rouleaux formations. Velocity profiles are presented for these cases. Unsteady flow clustered in local regions. Extra-cellular fluid flow generated by individual RBCs was examined using submicron fluorescent microspheres. The capabilities of confocal μPIV post-processing were verified using synthetic raw PIV data for validation. Cellular interactions and oscillating velocity profiles are presented, and 3D data are made available for computational model validation.     

受径向振荡激励的黏弹性液滴稳定性分析

当液滴受到外部周期性的径向激励时, 在其表面会形成驻波模式的不稳定, 这就是在球面上的Faraday不稳定问题. 不稳定的表面波的振荡频率根据流体物性参数和所施加激励条件的不同呈现为谐波或是亚谐波模式的振荡. 本文基于线性小扰动理论, 研究了受径向振荡体积力的黏弹性液滴表面波的不稳定性. 振荡的体积力导致动量方程为含有时间周期系数的Mathieu方程, 系统因此变成参数不稳定问题, 采用Floquet理论进行求解. 本模型中将黏弹性的特征处理为与流变模型参数相关的等效黏度, 从而简化了问题的求解. 基于对中性稳定曲线及增长率的分析, 研究了黏弹性参数对液滴稳定性的影响. 结果表明零剪切黏度和应变驰豫时间的增加具有抑制液滴表面波增长的作用, 提高了使液滴表面发生谐波不稳定的激励幅值. 随着振荡幅值的增加, 增长率不稳定的区域减少, 且随着振荡频率的增加, 液滴表面波增长率减小. 通过对增长率的分析可以得出, 应力松弛时间的增加使得增长率增加, 从而促进了液滴表面波的增长.      

Schlieren measurement of axisymmetric internal wave amplitudes

"Synthetic schlieren", which has been used to measure the amplitude of two-dimensional internal wave beams generated from an oscillating cylinder, is adapted to analyze axisymmetric internal waves generated by an oscillating sphere. This nonintrusive technique uses elementary inverse tomographic methods to measure the amplitude of the conical-structured wave beams everywhere in space and time. We compare the results with in situ probe measurements, and we examine the structure of the wave beams generated by a sphere oscillating at different amplitudes.     

Numerical study on the suppression of the vortex-induced vibration of an elastically mounted cylinder by a traveling wave wall

In the present paper, the commercial CFD code “Fluent” was employed to perform 2-D simulations of an entire process that included the flow around a fixed circular cylinder, the oscillating cylinder (vortex-induced vibration, VIV) and the oscillating cylinder subjected to shape control by a traveling wave wall (TWW) method. The study mainly focused on using the TWW control method to suppress the VIV of an elastically supported circular cylinder with two degrees of freedom at a low Reynolds number of 200. The cross flow (CF) and the inline flow (IL) displacements, the centroid motion trajectories and the lift and drag forces of the cylinder that changed with the frequency ratios were analyzed in detail. The results indicate that a series of small-scale vortices will be formed in the troughs of the traveling wave located on the rear part of the circular cylinder; these vortices can effectively control the flow separation from the cylinder surface, eliminate the oscillating wake and suppress the VIV of the cylinder. A TWW starting at the initial time or at some time halfway through the time interval can significantly suppress the CF and IL vibrations of the cylinder and can remarkably decrease the fluctuations of the lift coefficients and the average values of the drag coefficients; however, it will simultaneously dramatically increase the fluctuations of the drag coefficients.     

End reflections in a layered piezoelectric circular cylinder

End reflection phenomenon in a semi-infinitely long layered piezoelectric circular cylinder is constructed with modal data from a spectral decomposition of the differential operator governing its natural vibrations. These modal data consist of all propagating modes and edge vibrations and they constitute the basis for a wave function expansion of the reflection of waves arriving at the traction-free end of the cylinder. Without any other external stimulus, a passive reflection event occurs. This traction-free end condition is enforced at the Gaussian integration points over the end cross-section on the combination of incoming and reflected wave fields. Reflections due to monochromatic incoming axisymmetric (m = 0) and flexural (m = 1) waves are studied and two numerical examples are presented.For an incoming axisymmetric wave, there is a particular frequency that induces an end resonance, which is characterized by high (but finite) amplitudes of end displacements vis-a-vis those of neighboring (i.e., slightly different) frequencies. This phenomenon is illustrated in the two cylinder examples.It is possible to modify the passive reflection event by imposing some voltage distribution over the free end. For an oscillating end voltage that is out-of-phase with the incoming wave, it is possible to extract electrical energy from it, i.e., energy harvesting. Examples of such an oscillating voltage with a particular radial distribution are given, that illustrate the amount of extracted energy as a function of the frequency of the incident monochromatic wave.     

A calculating method of shock wave oscillating frequency due to turbulent shear layer fluctuations in supersonic flow

One of the more severe fluctuating pressure environments encountered in supersonic orhypersonic flows is the shock wave oscillation driven by interaction of a shock wave withboundary layer.The high intensity oscillating shock wave may induce structure resonanceof a high speed vehicle.The research for the shock oscillation used to adopt empirical orsemiempirical methods because the phenomenon is very complex.In this paper atheoretical solution on shock oscillating frequency due to turbulent shear layer fluctuationshas been obtained from basic conservation equations.Moreover,we have attained theregularity of the frequency of oscillating shock varying with incoming flow Mach numbersM_∞and turning angleθ.The calculating results indicate excellent agreement withmeasurements.This paper has supplied a valuable analytical method to study aeroelasticproblems produced by shock wave oscillation.     

基于修正光滑粒子流体动力学算法的低能量耗散数值波浪水池开发

目前, 无网格光滑粒子流体动力学SPH粒子法在波浪与结构物相互作用研究方面得到广泛应用, 但该方法模拟波浪远距离传播时, 常常面临严重的能量耗散问题, 导致波高非物理性降低, 给大范围海域、长时间作用下的波-物耦合作用研究带来一定困难. 对此, 本文采用一种核函数修正算法, 在确保粒子间相互作用对称性的同时, 改进压力梯度离散项的计算精度, 设法解决SPH方法中能量非物理性耗散的难题. 相较于前人减缓能量非物理性衰减的方法, 本文的修正SPH算法避免了自由液面搜索等复杂处理过程, 并能保证动量守恒特性. 数值结果中, 采用振荡液滴、规则波、不规则波等算例, 验证本修正SPH算法的准确性和有效性. 结果表明, 该修正SPH算法能准确模拟振荡液滴形态变化, 且动能保持较好守恒性. 通过数值水池与物理水池两者规则波与不规则波结果的对比分析表明, 基于本文修正SPH算法建立的数值波浪水池具有较好的抗能量衰减效果, 能实现长时间、远距离波浪传播的准确模拟. 此外, 本算法能在低光滑长度系数条件下, 实现精确模拟, 将极大缩减三维SPH模拟的时间, 从而节约计算成本.      

Modal analysis and harmonic response of resonators: An extension of a mapped orthogonal functions technique

A mapped orthogonal functions technique is extended for solving three-dimensional acoustical wave problems with taking into account electrical sources. For validation and illustration purposes, it is applied to the determination of the frequency spectrum of a piezoelectric disc resonator. Formulation is given for the calculation of the electric input admittance along with specializations for open- and short-circuit boundary conditions. In view of validation, normal frequencies and electric input admittance are obtained for specific geometries where a one-dimensional analytical model is working and compared against results from the analytical approach. Illustrative results, dispersion curves for a PZT5A resonator and electromechanical coupling coefficients as a function of the diameter to thickness ratio for the first radial mode of a respectively PZT5A, AlN and SiC resonator are also given.     

Water surface wave radiation generated by multiple cylinders oscillating with identical frequency

ＩｎｔｒｏｄｕｃｔｉｏｎＡｆｕｎｄａｍｅｎｔａｌｓｕｂｊｅｃｔｉｎｔｈｅｓｔｕｄｙｏｆｉｎｔｅｒａｃｔｉｏｎｂｅｔｗｅｅｎｗａｖｅｓａｎｄｂｏｄｙｉｓｔｏｓｏｌｖｅｔｈｅｔｗｏ＿ｄｉｍｅｎｓｉｏｎａｌｒａｄｉａｔｉｏｎｐｒｏｂｌｅｍｇｅｎｅｒａｔｅｄｂｙｃｙｌｉｎｄｅｒｓｏｓｃｉｌｌａｔｉｎｇｉｎｏｒｂｅｌｏｗｔｈｅｆｒｅｅｗａｔｅｒｓｕｒｆａｃｅ,ｗｈｉｃｈｈａｓａｔｔｒａｃｔｅｄｗｉｄｅａｔｔｅｎｔｉｏｎ．ＥｖｅｒｓｉｎｃｅＵｒｓｅｌｌｐｕｂｌｉｓｈｅｄｔｈｅｐａｐｅｒ［１］ｉｎｖｏｌｖｉｎ…     

Flow measurement around a model ship with propeller and rudder

For the design of hull forms with better resistance and propulsive performance, it is essential to understand flow characteristics, such as wave and wake development, around a ship. Experimental data detailing the local flow characteristics are invaluable for the validation of the physical and numerical modeling of computational fluid dynamics (CFD) codes, which are recently gaining attention as efficient tools for hull form evaluation. This paper describes velocity and wave profiles measured in the towing tank for the KRISO 138,000 m³ LNG carrier model with propeller and rudder. The effects of propeller and rudder on the wake and wave profiles in the stern region are clearly identified. The results contained in this paper can provide an opportunity to explore integrated flow phenomena around a model ship in the self-propelled condition, and can be added to the International Towing Tank Conference benchmark data for CFD validation as the previous KCS and KVLCC cases.     

Numerical study on mixed convection from a constant wall temperature circular cylinder in zero-mean velocity oscillating cooling flows

Fluid flow and heat transfer of mixed convection from a constant wall temperature circular cylinder in zero-mean velocity oscillating cooling flows have been simulated based on the projection method with two dimensional exponential stretched staggered cylindrical meshes. Cycle mean temperature and secondary streaming are obtained by the method of partial sums of the Fourier series. Present numerical results are validated by comparing the heat transfer results of free convection and the secondary streaming of pure oscillating flow over a circular cylinder to published experimental and numerical results. The complete structures of the cycle mean temperature and secondary streaming patterns are provided by numerical simulations over wide ranges of the Reynolds number, the Keulegan–Carpenter number and the Richardson number. Based on turning points of the curves of the overall Nusselt numbers versus Reynolds numbers and the characteristics of the cycle averaged temperature and flow patterns, the heat transfer can be divided into three linear regimes (conduction, laminar convection, and turbulent convection dominated regimes) and two non-linear transition regimes. The effects of wave directions, amplitudes, frequencies, and buoyancy forces on the enhancement of heat transfer are also investigated. The effective ranges of the governing parameters for heat transfer enhancement are identified.     

Coupled fluid structure analysis for wing 445.6 flutter using a fast dynamic mesh technology

The flow around wing 445.6 was modelled using Navier–Stokes equations and S-A model. The wing vibration and flow mesh deformation were computed using a fast dynamic mesh technology proposed by our own group. Wing 445.6 flutter was analysed through a strong coupling between the wing vibration and flow. The reduced flutter velocity was predicted and results are in good agreement with the experimental data. It is found that the subsonic flutter is mainly induced by the flow separation and the transonic and supersonic flutter are mainly caused by the oscillating shock wave and its induced flow separation. The positive aerodynamic work increases due to the oscillating shock wave when the subsonic flow becomes transonic reducing the flutter velocity. While the positive aerodynamic work induced by the oscillating shock wave decreases when the transonic flow becomes supersonic increasing the flutter velocity. That is why the transonic dip exists.     

Waves in a heavy two-layer liquid excited by an oscillating sphere

An approximate solution of an initial-boundary-value problem appropriate for the semiaxist>0 (t is time) is constructed for a system of integrodifferential equations which describes the waves excited in an initially stationary unbounded heavy two-layer fluid by a vertically oscillating sphere located at a distance from the interface that is significantly greater than its radius. The shape of the steady-state wave is found by passing to the limit as time increases indefinitely. The wave resistance experienced by the sphere during the transient process and in the steady-state regime is studied as a function of frequency. Nizhnii Novgorod. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 120–133, March–April, 1998.     

Justification of the Nonlinear Schrödinger Equation for the Evolution of Gravity Driven 2D Surface Water Waves in a Canal of Finite Depth

In 1968 V.E. Zakharov derived the Nonlinear Schrödinger equation for the two-dimensional water wave problem in the absence of surface tension, that is, for the evolution of gravity driven surface water waves, in order to describe slow temporal and spatial modulations of a spatially and temporarily oscillating wave packet. In this paper we give a rigorous proof that the wave packets in the two-dimensional water wave problem in a canal of finite depth can be approximated over a physically relevant timespan by solutions of the Nonlinear Schrödinger equation.     

Reflection and transmission of internal solitary waves across a barrier

Transmission and reflection of an internal solitary wave on a two-layer fluid as it encounters a barrier is considered. An edgelayer theory is used to derive the reduced boundary conditions relevant to the shallow-water equations. It is found that for a moderate height of the barrier the incident wave is not affected except by being phase shifted. On the other hand, for a relatively high barrier the incident wave splits into a transmitted wave and a reflected wave. The reflected wave evolves into a small soliton with decaying ripples, and the transmitted wave evolves into a larger soliton and an oscillating tail.     

Oscillatory line source for water waves in shear flow

The linearized water-wave radiation problem for the oscillating 2D submerged source in an inviscid shear flow with a free surface is investigated analytically. The vorticity is uniform, with zero velocity at the free surface. Then there will be at most two emitted waves, and no Doppler effects. Exact far-field waves are derived, with radiation conditions applied at infinity. An upstream wave will always exist, whereas the downstream wave exists only when the angular frequency of oscillation exceeds the vorticity. The wave radiation problem is solved also for oscillating vortex and dipoles. The amplitudes and energy fluxes are calculated.