A GENERAL FORMULA FOR CALCULATING FORCES ON A 2-D ARBITRARY BODY IN INCOMPRESSIBLE FLOW

In the present paper, a general integral equation is presented to calculate the forces exerted on a two-dimensional (2-D) body of arbitrary shape immersed in unsteady, incompressible flows. By finding the general solutions of a set of Laplace equations with particular boundary conditions, the equation can be simplified to produce a simplified formula for calculating the forces. The simplified formula consists of three parts, representing contributions from different physical phenomena: added mass force and/or inertial force in inviscid flow, the force caused by the deformation of fluid and viscosity and the force caused by the convection of fluid with nonzero circulation. It can be applied to any 2-D arbitrary body in viscous or inviscid, steady or unsteady incompressible flow. As the formula excludes either temporal derivatives of velocity or spatial derivatives of vorticity in the flow field, the numerical errors contained in the numerical solution of velocity and vorticity fields will not be magnified, and therefore the resulting force calculated is more accurate. Most importantly, the formula presents an alternative method for obtaining the added mass of a 2-D body of arbitrary shape accelerating in a fluid. For bodies of simple shape, such as a circle, ellipse and plate, the added masses predicted using the present method are in agreement with that obtained by conventional methods. For bodies of complex shape, the present method only requires the calculation of the first two coefficients of the conformal transformation and cross-sectional area.     

An improved HSFR method for natural vibration analysis of an immersed cylinder pile with a tip mass

Immersed cylinder piles are usually modelled as immersed cantilever cylinder columns carrying a tip mass and rotary moment of inertia. In this paper, the equations of motion of an immersed cylinder pile along transversal modes of vibration are developed. Compressibility of water and structural damping are included in the formulation. Natural frequencies of the immersed pile are obtained from the developed equations using harmonic sweep frequency response analyses. The proposed method is applied to numerical examples, and the results obtained are shown satisfactory when compared to other numerical solutions in the literature, or to finite element solutions and experimental data.     

New analytical solutions to water wave radiation by vertical truncated cylinders through multi-term Galerkin method

New analytical solutions to water wave radiation by vertical truncated circular cylinders are developed based on linear potential flow theory. Two typical cylinder configurations of a surface-piercing cylinder and a submerged floating cylinder are considered. The multi-term Galerkin method is employed in the solution procedure, in which the fluid velocity on the interface between different regions is expanded into a set of basis function involving the Gegenbauer polynomials, and the cube-root singularity of fluid velocity at the side edges of the truncated cylinders is correctly modeled. The present solutions have the merits of very rapid convergence. The results with six-figure accuracy for added mass and radiation damping can be obtained using a few truncated numbers in the basis function for three motions (surge, heave and roll). The calculated results of the present solutions agree well with that by a higher-order boundary element method solution. Calculation examples are presented to investigate the influence of the motion frequency on the added mass and the radiation damping of the truncated cylinders with different geometric parameters. The present solutions can be used as a reliable benchmark for numerical solutions to water wave radiation by complicated structures.

     

STUDY ON THE DYNAMIC BEHAVIOR OF AXIALLY MOVING RECTANGULAR PLATES PARTIALLY SUBMERSED IN FLUID

The natural frequencies, complex modes and critical speeds of an axially moving rectangular plate, which is partially immersed in a fluid and subjected to a pretension, are investigated. The effects of free surface waves, compressibility and viscidity of the fluid are neglected in the analysis. The subsection functions are used to describe the discontinuous characteristics of the system due to partial immersion. The classical thin plate theory is adopted to formulate the equations of motion of a vibrating plate. The velocity potential and Bernoulli's equation are used to describe the fluid pressure acting on the moving plate. The effect of fluid on the vibrations of the plate may be equivalent to the added mass on the plate. The effects of distance ratio, moving speed, immersed-depth ratio, boundary conditions, stiffness ratio and aspect ratio of the plate as well as the fluid-plate density ratios on the free vibrations of the moving plate-fluid system are investigated.     

Improvement of mass source/sink for an immersed boundary method

An improved immersed boundary method using a mass source/sink as well as momentum forcing is developed for simulating flows over or inside complex geometries. The present method is based on the Navier–Stokes solver adopting the fractional step method and a staggered Cartesian grid system. A more accurate formulation of the mass source/sink is derived by considering mass conservation of the virtual cells in the fluid crossed by the immersed boundary. Two flow problems (the decaying vortex problem and uniform flow past a circular cylinder) are used to validate the proposed formulation. The results indicate that the accuracy near the immersed boundary is improved by introducing the accurate mass source/sink. Copyright © 2006 John Wiley & Sons, Ltd.     

ADDED MASS AND OSCILLATION FREQUENCY FOR A CIRCULAR CYLINDER SUBJECTED TO VORTEX-INDUCED VIBRATIONS AND EXTERNAL DISTURBANCE

This paper reports results from an experiment with a lightly damped elastically mounted rigid cylinder subjected to constant flow velocity. The cylinder was allowed to vibrate in the cross-flow direction and was fixed in the flow direction. The Reynolds numbers varied from 10⁴to 6×10⁴. The added mass for the freely vibrating cylinder agreed well with the results found by others in driven cylinder tests. The predicted natural frequency based on the measured added mass was approximately equal to the measured mean oscillation frequency. The added mass calculated from one oscillation cycle to the next varied considerably. The oscillation frequency from one oscillation to the next corresponded to the natural frequency including the added mass for the same cycle. By movement of the attachment point of the elastic members to the external structure a disturbance could be added to the normal vortex-induced vibrations (VIV) response. When an external disturbance was introduced at a frequency other than the VIV frequency, the added mass coefficient was found to be weakly influenced by the external harmonic disturbance.     

Simulation of incompressible viscous flows around moving objects by a variant of immersed boundary‐lattice Boltzmann method

A variant of immersed boundary‐lattice Boltzmann method (IB‐LBM) is presented in this paper to simulate incompressible viscous flows around moving objects. As compared with the conventional IB‐LBM where the force density is computed explicitly by Hook's law or the direct forcing method and the non‐slip condition is only approximately satisfied, in the present work, the force density term is considered as the velocity correction which is determined by enforcing the non‐slip condition at the boundary. The lift and drag forces on the moving object can be easily calculated via the velocity correction on the boundary points. The capability of the present method for moving objects is well demonstrated through its application to simulate flows around a moving circular cylinder, a rotationally oscillating cylinder, and an elliptic flapping wing. Furthermore, the simulation of flows around a flapping flexible airfoil is carried out to exhibit the ability of the present method for implementing the elastic boundary condition. It was found that under certain conditions, the flapping flexible airfoil can generate larger propulsive force than the flapping rigid airfoil. Copyright © 2009 John Wiley & Sons, Ltd.     

Modal analysis of a spinning disk in a dense fluid as a model for high head hydraulic turbines

In high head Francis turbines and pump-turbines in particular, Rotor Stator Interaction (RSI) is an unavoidable source of excitation that needs to be predicted accurately. Precise knowledge of turbine dynamic characteristics, notably the variation of the rotor natural frequencies with rotation speed and added mass of the surrounding water, is essential to assess potential resonance and resulting amplification of vibrations. In these machines, the disk-like structures of the runner crown and band as well as the head cover and bottom ring give rise to the emergence of diametrical modes and a mode split phenomenon for which no efficient prediction method exists to date. Fully coupled Fluid–Structure Interaction (FSI) methods are too computationally expensive; hence, we seek a simplified modeling tool for the design and the expected-life prediction of these turbines.We present the development of both an analytical modal analysis based on the assumed mode approach and potential flow theory, and a modal force Computational Fluid Dynamics (CFD) approach for rotating disks in dense fluid. Both methods accurately predict the natural frequency split as well as the natural frequency drift within 7.9% of the values measured experimentally. The analytical model explains how mode split and drift are respectively caused by linear and quadratic dependence of the added mass with relative circumferential velocity between flexural waves and fluid rotation.     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界,避免生成复杂贴体网格及动网格,极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法,将简化附加体积力以源项形式嵌入动量方程迭代中,通过用户自定义函数对CFD软件FLUENT二次开发,实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证,并探讨了插值函数对计算精度的影响。研究表明,通过引入浸入边界模型,能够提高计算效率,并实现结构网格背景下复杂边界和动边界的高效建模。     

浸没边界-简化热格子Boltzmann方法研究及其应用

针对流固耦合传热问题,本文提出了一种基于浸没边界-简化热格子玻尔兹曼方法(immersed boundary method-simplified thermal lattice Boltzmann method,IB-STLBM)的耦合模型.不同于传统的格子玻尔兹曼方法使用分布函数演化流场和温度场,简化热格子玻尔兹曼方法(simplified thermal lattice Boltzmann method,STLBM)的演化过程不需要依赖分布函数,只涉及平衡态分布函数和非平衡态分布函数,能够直接演化宏观量,极大减小了计算过程中所占用的虚拟内存,简化了边界条件的实现方式,同时具有较高的稳定性.传统的浸没边界法对流场的计算采用欧拉网格,对固体边界采用拉格朗日网格,认为固体边界是对流场产生某种体积力.在应用浸没边界法时,汲取介观的思想,把固体的介入看作是对流场的干扰,打破了固体附近流体介观微团颗粒原始的平衡状态,这种干扰可以看作是在耦合边界上产生的一个非平衡项,可用非平衡态分布函数来表示.基于此,在模型中浸没边界法与简化热格子玻尔兹曼方法更紧密联系在一起,更大程度发挥二者的优点,整个计算过程更加简单直观,符合物理特性.通过对热圆柱绕流和内含热颗粒的封闭方腔自然对流问题的模拟以及对其结果的分析,验证了该算法在求解流固耦合传热问题的有效性和可行性.      

A divergence‐free interpolation scheme for the immersed boundary method

The immersed boundary approach for the modeling of complex geometries in incompressible flows is examined critically from the perspective of satisfying boundary conditions and mass conservation. It is shown that the system of discretized equations for mass and momentum can be inconsistent, if the velocity is used in defining the force density to satisfy the boundary conditions. As a result, the velocity is generally not divergence free and the pressure at locations in the vicinity of the immersed boundary is not physical. However, the use of the pseudo‐velocities in defining the force density, as frequently done when the governing equations are solved using a fractional step or projection method, combined with the use of the specified velocity on the immersed boundary, is shown to result in a consistent set of equations which allows a divergence‐free velocity but, depending on the time step, is shown to have the undesirable effects of inaccurately satisfying the boundary conditions and allowing a significant permeability of the immersed boundary. If the time step is reduced sufficiently, the boundary conditions on the immersed boundary can be satisfied. However, this entails an unacceptable increase in computational expense. Two new methods that satisfy the boundary conditions and allow a divergence‐free velocity while avoiding the increased computational expense are presented and shown to be second‐order accurate in space. The first new method is based on local time step reduction. This method is suitable for problems where the immersed boundary does not move. For these problems, the first new method is shown to be closely related to the second new method. The second new method uses an optimization scheme to minimize the deviation from the interpolation stencil used to represent the immersed boundary while ensuring a divergence‐free velocity. This method performs well for all problems, including those where the immersed boundary moves relative to the grid. Additional results include showing that the force density that is added to satisfy the boundary conditions at the immersed boundary is unbounded as the time step is reduced and that the pressure in the vicinity of the immersed boundary is unphysical, being strongly a function of the time step. A method of computing the total force on an immersed boundary which takes into account the specifics of the numerical solver used in the iterative process and correctly computes the total force irrespective of the residual level is also presented. Copyright © 2007 John Wiley & Sons, Ltd.     

Oscillations of a Circular Cylinder in a Linearly Stratified Fluid

The plane problem of the motions of a three-layer fluid initiated by the oscillations of a circular cylinder is solved in the linear formulation in the Boussinesq approximation. The cylinder is completely immersed in the linearly stratified middle layer, and the upper and lower layers are homogeneous and bounded by rigid horizontal walls. The fluid is assumed to be ideal and incompressible. The added mass and damping coefficients are calculated as functions of the oscillation frequency of the cylinder and the layer thicknesses.     

Experimental investigation of added mass effects on a Francis turbine runner in still water

The mechanical design of hydraulic turbines is conditioned by the dynamic response of the runner that is usually estimated by a computational model. Nevertheless, the runner has complex boundary conditions that are difficult to include in the computational model. One of these boundary conditions is the water in which the runner is submerged. The effect of the added mass and damping of water can modify considerably the natural frequencies of the runner. An experimental investigation in a reduced scale model of a turbine runner, using modal analysis, was carried out. Several impact tests with the runner freely suspended in air and in water were done. The response was measured with accelerometers located in different positions of the runner. From the modal analysis, the natural frequencies, damping ratios, and mode-shapes were determined. The same mode-shapes obtained in air were obtained in water but with lower natural frequencies and higher damping ratios in water. The difference in the natural frequencies is shown to be dependant basically on the added mass effect of the water and not on its added damping. This difference also depends on the geometry of the mode, presenting different values for different mode-shapes. Using nondimensional values, the reduction in the natural frequencies can be extrapolated to other Francis runners presenting similar geometrical characteristics.     

An immersed boundary method for simulation of inviscid compressible flows

In this paper, an immersed boundary method for simulating inviscid compressible flows governed by Euler equations is presented. All the mesh points are classified as interior computed points, immersed boundary points (interior points closest to the solid boundary), and exterior points that are blanked out of computation. The flow variables at an immersed boundary point are determined via the approximate form of solution in the direction normal to the wall boundary. The normal velocity is evaluated by applying the no‐penetration boundary condition, and therefore, the influence of solid wall in the inviscid flow is taken into account. The pressure is computed with the local simplified momentum equation, and the density and the tangential velocity are evaluated by using the constant‐entropy relation and the constant‐total‐enthalpy relation, respectively. With a local coordinate system, the present method has been extended easily to the three‐dimensional case. The present work is the first endeavor to extend the idea of hybrid Cartesian/immersed boundary approach to compressible inviscid flows. The tedious task of handling multi‐valued points can be eliminated, and the overshoot resulting from the extrapolation for the evaluation of flow variables at exterior points can also be avoided. In order to validate the present method, inviscid compressible flows over fixed and moving bodies have been simulated. All the obtained numerical results show good agreement with available data in the literature. Copyright © 2014 John Wiley & Sons, Ltd.     

基于反馈力浸入边界法模拟复杂动边界流动

浸入边界法是模拟流固耦合的重要数值方法之一。本文采用反馈力浸入边界方法,对旋转圆柱和水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟。其中,固体边界采用一系列离散的点近似代替,流体为不可压缩牛顿流体,使用笛卡尔自适应加密网格,利用有限差分法进行求解。固体对流场的作用通过构造适宜的反馈力函数实现。本文首先通过旋转圆柱绕流的计算结果同实验结果进行对比,吻合较好,验证了该计算方法的可靠性。然后针对水电站水力过渡过程中水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟,得到导叶动态绕流后的流场分布特性和涡结构的演化特性。     

Vortex and wake-induced vibrations of a tandem arrangement of two flexible circular cylinders with near wake interference

Results showing the dynamic response of a tandem arrangement of two vertical high aspect ratio (length over diameter) and low mass ratio (mass over mass of displaced fluid) flexible cylinders vibrating at low mode number are presented in this paper. Two circular cylinder models were aligned with the flow, so the downstream or trailing cylinder was immersed in the wake of the leading one. Centre-to-centre distances from 2 to 4 diameters were studied. The models were very similar in design, with external diameters of 16 mm and a total length of 1.5 m. Reynolds numbers up to 12 000 were achieved with reduced velocities, based on the fundamental natural frequency of the downstream cylinder in still water, up to 16. The trailing model had a mass ratio of 1.8 with a combined mass-damping parameter of 0.049, whilst the corresponding figures for the leading cylinder were 1.45 and 0.043, respectively. The dynamic response of the trailing model has been analysed by studying cross-flow and in-line amplitudes, dominant frequencies and modal amplitudes. The dynamic response of the leading one is analysed by means of its cross-flow amplitudes and dominant frequencies and it is also related to the motion of the trailing cylinder by studying the synchronisation between their instantaneous cross-flow motions. Planar digital particle image velocimetry (DPIV) was used to visualise the wake. Different response regimes have been identified based on the type of oscillations exhibited by the cylinders: vortex-induced (VIV), wake-induced (WIV) or combinations of both.     

Application of Differing Forcing Function Models on Simulated Flow Past an Oscillating Cylinder in a Uniform Low Reynolds Number Flow

A Computational Fluid Dynamics (CFD) model is presented for the uniform viscous two dimensional flow past an oscillating cylinder at low Reynolds number. Numerical simulations are made to study the effect of differing forced induced oscillation mechanisms with a large range of cylinder forcing frequencies. In the first case sinusoidal velocity slip boundary conditions are adopted for the cylinder surface to simulate cylinder oscillation. The implication suggests that no modification or additional term need to be added to the Navier-Stokes equations. In the second case this time extra body force terms which are assumed to account for velocity effects due to cylinder movement are included in the Navier-Stokes equations with the imposition of same boundary conditions. Drag and lift coefficients are extracted from present numerical results and other detailed computations of these coefficients are made at a Reynolds number of 80 and an amplitude-to diameter ratio 0.14. The results are found to be in agreement with each other at low force driving frequencies below and near lock-in. However, differences are found at higher frequencies above lock-in. Agreement are also found with experimental results at some frequency ranges.     

Equivalent Cylinder Systems Representing the Soil Matrix in Diffusion-Reaction Models for an Aggregated Soil

The structure of an aggregated soil or a dual-porosity medium is characterized by a function called the generalized surface-area-to-volume ratio, which can be measured as a distance distribution. For an isotropic structure the distance measurements may be done in the digitized image of a cross-section. A diffusion-reaction process in the soil matrix can be described by means of an equivalent cylinder system which preserves the generalized surface-area-to-volume ratio of the structure. The contribution of each cylinder radius is described by a weight function, which is obtained by transforming the measured cross-sectional distance distribution. In order to test the method, a diffusion-reaction model from the literature has been applied to spheres with a lognormal radius distribution and to an irregular, generated structure. Both geometries can be represented by just a few cylinder radii and cylinder weights. The method works best for dual-porosity structures with a convex shape. Deviations occur if parts of the matrix structure have a concave shape.     

ON THE FREE VIBRATION OF A SUBMERGED FGM HOLLOW SPHERE

The free vibration of a functionally graded material hollow sphere submerged in a compressible fluid medium is exactly analyzed. The sphere is assumed to be spherically isotropic with material constants being inhomogeneous along the radial direction. By employing a separation technique as well as the spherical harmonics expansion method, the governing equations are simplified to an uncoupled second-order ordinary differential equation, and a coupled system of two such equations. Solutions to these equations are given when the elastic constants and the mass density are power functions of the radial coordinate. Numerical examples are finally given to show the effect of the material gradient on the natural frequencies. The project is supported by the National Natural Sciences Foundation of China(No. 19872060).