Vortex-induced vibrations of pipes conveying fluid in the subcritical and supercritical regimes

In this paper, the vortex-induced vibrations of a hinged–hinged pipe conveying fluid are examined, by considering the internal fluid velocities ranging from the subcritical to the supercritical regions. The nonlinear coupled equations of motion are discretized by employing a four-mode Galerkin method. Based on numerical simulations, diagrams of the displacement amplitude versus the external fluid reduced velocity are constructed for pipes transporting subcritical and supercritical fluid flows. It is shown that when the internal fluid velocity is in the subcritical region, the pipe is always vibrating periodically around the pre-buckling configuration and that with increasing external fluid reduced velocity the peak amplitude of the pipe increases first and then decreases, with jumping phenomenon between the upper and lower response branches. When the internal fluid velocity is in the supercritical region, however, the pipe displays various dynamical behaviors around the post-buckling configuration such as inverse period-doubling bifurcations, periodic and chaotic motions. Moreover, the bifurcation diagrams for vibration amplitude of the pipe with varying internal fluid velocities are constructed for each of the lowest four modes of the pipe in the lock-in conditions. The results show that there is a significant difference between the vibrations of the pipe around the pre-buckling configuration and those around the post-buckling configuration.     

Multi-modes approach to modelling of vortex-induced vibration

In this work the fluid–structure interactions are considered by investigating a straight but slender pipe interacting with uniform water flow. Two configurations are studied, namely vertically and horizontally positioned pipes, which are modelled as an Euler–Bernoulli beam with flexural stiffness. Both pretension and length-wise mass distribution are considered. The structure is assumed to be moving only in the direction normal to flow (cross-flow motion) hence its in-line motion is neglected. The external fluid force acting on the structure is the result of the action of sectional vortex-induced drag and lift forces. Only mean drag force is considered, with time varying lift force modelled using a non-linear oscillator equation of the Van der Pol type. The obtained coupled system of non-linear partial differential equations is simplified employing Galerkin-type discretisation. The resulting ordinary differential equations are solved numerically providing multi-mode approximations of cross-flow displacement and non-dimensional lift coefficient. The comparison between the responses of vertical and horizontal structures shows that, as expected, due to a balancing between pretension and weight, in general a higher amplitude of vibration is observed for the vertical configuration than in the same location along the pipe for the horizontal configuration in the lower part of the structure. However, lower amplitudes are obtained in the upper part of the pipe. The horizontal configuration solutions are identical in symmetrical locations along the pipe due to constant pretension. The influence of the wake equation coefficients and the fluid force coefficients on the response amplitudes has been also considered together with the length of the pipe and pretension level, and the appropriate response curves are included. Finally, for the higher mode approximations it has been shown that the vibrations level at lower frequencies is predicted reasonably well by retaining only a small subset of modes.     

Stability analysis of a fractional viscoelastic plate strip in supersonic flow under axial loading

The stability of a viscoelastic plate strip, subjected to an axial load with the Kelvin–Voigt fractional order constitutive relationship is studied. Based on the classical plate theory, the structural formulation of the plate is obtained by using the Newton’s second law and the aerodynamic force due to the fluid flow is evaluated by piston theory. The Galerkin method is employed to discretize the equation of motion into a set of ordinary differential equations. To determine the stability margin of plate the obtained set of ordinary differential equations are solved using the Laplace transform method. The effects of variation of the governing parameters such as axial force, retardation time, fractional order and boundary conditions on the stability margin of fractional viscoelastic panel are investigated and finally some conclusions are outlined.     

Compliance effects on the torsional flow of a viscoelastic fluid

The effects of transducer compliance on transient stress measurements in torsional flows of a viscoelastic fluid are investigated theoretically. The analysis is based on the torsional flow of an upper-convected Maxwell fluid between a rotating and ‘stationary’ disk, which is allowed to twist and displace axially as a result of the stresses exerted on the disk by the fluid. An approximate analytical solution to the governing equations is obtained using a standard perturbation method. Results of the analysis are used to examine how the fluid velocity is altered by the motion of the stationary disk and to gain insight on how transient stress measurements are affected by transducer compliance. The analysis shows that compliance effects increase with applied shear rate and that the effects of torsional and axial compliance are coupled in measurements of the shear stress and first normal stress difference.     

Nonlinear dynamics of extensible viscoelastic cantilevered pipes conveying pulsatile flow with an end nozzle

Nonlinear dynamics of an extensible cantilevered pipe conveying pulsating flow is considered in this paper. The fluid flow fluctuates harmonically and exhausts via a nozzle attached to the end of the pipe. Taking into account the extensibility assumption, the coupled nonlinear lateral–longitudinal equations of motion are derived using Hamilton's principle and discretized via Galerkin's method. The adaptive time step Adams algorithm is applied to extract the time response, and then the bifurcation, power spectral density and phase plane maps are plotted for some case studies. Effects of some geometrical parameters such as flow mass, pulsating flow frequency, gravity, nozzle mass and nozzle aspect ratio parameters are studied on the dynamics of such system and the validity of extensibility assumption is investigated and some conclusions are drawn.     

Aeroacoustical coupling in a ducted shallow cavity and fluid/structure effects on a steam line

A pure tone phenomenon has been observed at 460 Hz in a piping steam line. The acoustical energy has been identified to be generated in an open gate valve and to be of cavity noise type. This energy is then transmitted to the main pipe by fluid/structure coupling. The objectives here are to display the mechanism of the flow acoustic coupling in the cavity and in the duct through an aeroacoustical analysis and to understand the way of energy transfer from the fluid to the main pipe through a vibroacoustical analysis. Concerning the first objective, an experimental study by means of 2/7 scale models in air is analysed by means of numerical flow simulation. The flow acoustic phenomena are modelled by computing the Euler equations. Two different computations are carried out: in the first one, a pure Euler modelling is used, in the second one, a boundary layer obtained from experimental data is introduced in the computation in order to have a realistic flow profile upstream the cavity. The boundary layer flow profile appears to be essential to recover the experimentally observed coupling between the shear-layer instability and the acoustical transverse mode of the pipe. The numerical results confirm that the second aerodynamic mode is responsible for the oscillation. While the predicted frequency agrees about 1% with the scale model experiments, the predicted amplitude is approximately 15 dB too low. For the second objective, fluid/structure coupling in the main pipe is studied using two fully coupled methods. The first method consists in a modal analysis of the line using a fluid–structure finite element model. The second one is based on the analysis of dispersion diagrams derived from the local equations of cylindrical shells filled with fluid. The way of energy transfer in transverse acoustical waves coupled with flexion-ovalization deformations of the pipe is highlighted using both methods. The dispersion diagrams allow a fast and accurate analysis. The modal analysis using a finite-element model may complete the first one with quantitative data. The link between the fluid/acoustic and the fluid/structure analysis is then the excitation of the transverse acoustical mode of the duct.     

Dynamic instability of inclined cables under combined wind flow and support motion

In this paper an inclined nearly taut stay, belonging to a cable-stayed bridge, is considered. It is subject to a prescribed motion at one end, caused by traveling vehicles, and embedded in a wind flow blowing simultaneously with rain. The cable is modeled as a non-planar, nonlinear, one-dimensional continuum, possessing torsional and flexural stiffness. The lower end of the cable is assumed to undergo a vertical sinusoidal motion of given amplitude and frequency. The wind flow is assumed uniform in space and constant in time, acting on the cable along which flows a rain rivulet. The imposed motion is responsible for both external and parametric excitations, while the wind flow produces aeroelastic instability. The relevant equations of motion are discretized via the Galerkin method, by taking one in-plane and one out-of-plane symmetric modes as trial functions. The two resulting second-order, non-homogeneous, time-periodic, ordinary differential equations are coupled and contain quadratic and cubic nonlinearities, both in the displacements and velocities. They are tackled by the Multiple Scale perturbation method, which leads to first-order amplitude-phase modulation equations, governing the slow dynamics of the cable. The wind speed, the amplitude of the support motion and the internal and external frequency detunings are set as control parameters. Numerical path-following techniques provide bifurcation diagrams as functions of the control parameters, able to highlight the interactions between in-plane and out-of-plane motions, as well as the effects of the simultaneous presence of the three sources of excitation.     

粘弹流体轴对称流动的格子Boltzmann方法

基于BGK碰撞模型，通过在迁移方程中引入作用力项，建立了粘弹流体的轴对称格子Boltzmann模型．通过Chapman-Enskog展开，获得了准确的柱坐标下轴对称宏观流动方程．采用双分布函数对运动方程和本构方程进行迭代求解，模拟分析了粘弹流体管道流动，获得了流场中的速度和构型张量的分布，通过与解析解进行比较，验证了模型的准确性．研究了作为粘弹流体流动基准问题的收敛流动，对涡旋位置进行了定量分析，将回转长度的计算结果与有限体积法进行了比较，两种数值结果十分吻合．研究结果表明，模型能够准确表征粘弹流体的轴对称流动，具有较广阔的应用前景．     

Migration of a solid particle in the vicinity of a plane fluid–fluid interface

The slow migration of a small and solid particle in the vicinity of a gas–liquid, fluid–fluid or solid–fluid plane boundary when subject to a gravity or an external flow field is addressed. By contrast with previous works, the advocated approach holds for arbitrarily shaped particles and arbitrary external Stokes flow fields complying with the conditions on the boundary. It appeals to a few theoretically established and numerically solved boundary-integral equations on the particle’s surface. This integral formulation of the problem allows us to provide asymptotic approximations for a distant boundary and also, implementing a boundary element technique, accurate numerical results for arbitrary locations of the boundary. The results obtained for spheroids, both settling or immersed in external pure shear and straining flows, reveal that the rigid-body motion experienced by a particle deeply depends upon its shape and also upon the boundary location and properties.     

Computing axisymmetric,laminar internal flows

A simple computational scheme is developed to compute laminar flows inside axisymmetric ducts. It is based on the Keller box method where the equations are approximated at the centre of the downstream face of each computational box. The coupling between the pressure gradient and the velocities for internal flow has been observed to introduce stability problems for the Keller box method that are not present for external, boundary layer flow problems. The difference scheme for the velocities is coupled to an iterative scheme to solve for the pressure gradient at each axial step. Example results for developing flow in a pipe and in a 2° conical diffuser are presented.     

Deformation of a fluid-filled compliant cylinder in a uniform flow

We investigated surface compliance effects of a fluid-filled object in flow on its shape and internal flow through numerical simulation. A two-dimensional compliant cylinder containing fluid in a flow is a simple model of a cell, e.g. an erythrocyte, leukocyte or platelet. The thin membrane of the cylinder consisted of a network of mass-spring-damper (MSD) systems, representing its mechanical characteristics. We assumed that the stiffness and damping coefficients were those of latex gum. The two-dimensional flow inside and outside the membrane was obtained by solving the two-dimensional Navier–Stokes equations by using the finite element scheme at Re=400, based on the external flow velocity and diameter of an initial circular cylinder. The deformation of the membrane was calculated by solving the equation of motion for an MSD system by using the fourth-order Runge-Kutta method. The compliant cylinder deformed more if its stiffness was smaller than that of latex gum. The initial circular section of the cylinder became oval, with a flat front and a convex rear. The aspect ratio of the lateral to streamwise axis length of the oval became larger than unity, and increased with decreasing stiffness. The drag coefficient of the oval cylinder became larger than that of the circular cylinder, and increased with decreasing stiffness. The partial vibration at the rear, caused by shedding vortices, induced oscillating internal flows between two antinodes of the vibrating membrane. Since the object with smaller stiffness had higher ductility, velocity fluctuations of the external flow influenced the internal flow of the compliant object through deformation of the membrane.     

A perturbation solution for non-linear vibration of uniformly curved pipes conveying fluid

The first-order non-linear interactions between the pipe structure and the flowing fluid are considered to formulate the governing equations of motion for the in-plane vibration of a circular-arc pipe containing flowing fluid. The forces and moments induced in a pipe element by the flowing fluid are analyzed as functions of the instantaneous local curvature of the pipe. The flow field is assumed to be one-dimensional, incompressible and of uniform flow, and to remain independent of pipe motion. For a fixed-end circular-arc pipe with arbitrary arc angle, the non-linear governing equations are solved by the method of multiple scales in conjunction with the Bubnov-Galerkin method. The non-linear solutions indicate that the vibrational behavior of the system can differ substantially from that predicted by a linear analysis.     

Thrust augmentation of flapping airfoils in low Reynolds number flow using a flexible membrane

The unsteady aerodynamic thrust and aeroelastic response of a two-dimensional membrane airfoil under prescribed harmonic motion are investigated computationally with a high-order Navier–Stokes solver coupled to a nonlinear membrane structural model. The effects of membrane prestress and elasticity are examined parametrically for selected plunge and pitch–plunge motions at a chord-based Reynolds number of 2500. The importance of inertial membrane loads resulting from the prescribed flapping is also assessed for pure plunging motions. This study compares the period-averaged aerodynamic loads of flexible versus rigid membrane airfoils and highlights the vortex structures and salient fluid–membrane interactions that enable more efficient flapping thrust production in low Reynolds number flows.     

A numerical study of fully developed laminar flows in pipes with two planar curvatures

The spectral element method is applied on unstructured tetrahedral elements to solve the Navier–Stokes equations for fully developed laminar flow in pipes with two planar curvatures. Specific implementations of the spectral element method to double curved pipes and parallelization are described. Previous studies on flows in pipes focused on constant curvature or torsion geometries, as well as pipes with varying curvature. This study focuses on the periodic variation of both the curvature as well as torsion by analysing a pipe having two planar curvatures. The effects of the three parameters defining the pipe are studied to isolate the curvature and torsion effect on the magnitude and angle of the secondary flow. Furthermore, the geometric effects on the wall shear stress are studied, as it is an important fluid flow property, especially in blood flows. Copyright © 2006 John Wiley & Sons, Ltd.     

DYNAMICS OF A SUBMERGED AND INCLINED CONCENTRIC PIPE SYSTEM WITH INTERNAL AND EXTERNAL FLOWS

In this paper, we formulate a mathematical model to study the dynamics of submerged and inclined concentric pipes with different lengths. The governing equations of motion for the inner pipe are derived under small deformation assumptions and with the consideration of gravitational forces, turbulent boundary layer thickness of external flow, fluid frictional forces, and inertia effects. We obtain discretized dynamical equations using spatial finite-difference schemes and calculate the resonant frequencies of a particular pipe system design. In addition, by varying the operating conditions, we identify a few critical parameters pertaining to the proper design of such pipe systems.     

3D free‐surface flow computation using a RANSE/Fourier–Kochin coupling

A coupling method for numerical calculations of steady free‐surface flows around a body is presented. The fluid domain in the neighbourhood of the hull is divided into two overlapping zones. Viscous effects are taken in account near the hull using Reynolds‐averaged Navier–Stokes equations (RANSE), whereas potential flow provides the flow away from the hull. In the internal domain, RANSE are solved by a fully coupled velocity, pressure and free‐surface elevation method. In the external domain, potential‐flow theory with linearized free‐surface condition is used to provide boundary conditions to the RANSE solver. The Fourier–Kochin method based on the Fourier–Kochin formulation, which defines the velocity field in a potential‐flow region in terms of the velocity distribution at a boundary surface, is used for that purpose. Moreover, the free‐surface Green function satisfying this linearized free‐surface condition is used. Calculations have been successfully performed for steady ship‐waves past a serie 60 and then have demonstrated abilities of the present coupling algorithm. Copyright © 2003 John Wiley & Sons, Ltd.     

FLOW-INDUCED INTERNAL RESONANCES AND MODE EXCHANGE IN HORIZONTAL CANTILEVERED PIPE CONVEYING FLUID （Ⅱ）

Based on the nonlinear mathematical model of motion of a horizontally can-tilevered rigid pipe conveying fluid, the 3:1 internal resonance induced by the minimum critical velocity is studied in details. With the detuning parameters of internal and primary resonances and the amplitude of the external disturbing excitation varying, the flow in the neighborhood of the critical flow velocity yields that some nonlinearly dynamical behaviors occur in the system such as mode exchange, saddle-node, Hopf and co-dimension 2 bifurcations. Correspondingly, the periodic motion losses its stability by jumping or flutter, and more complicated motions occur in the pipe under consideration. The good agreement between the analytical analysis and the numerical simulation for several parameters ensures the validity and accuracy of the present analysis.     

A VORTICITY–STREAMFUNCTION FORMULATION FOR STEADY INCOMPRESSIBLE TWO-DIMENSIONAL FLOWS

A vorticity–streamfunction formulation for incompressible planar viscous flows is presented. The standard kinematic field equations are discretized using centred finite difference schemes and solved in a coupled way via a Newton-like linearization scheme. The linearized system of partial differential equations is handled through the restarting linear GMRES algorithm, preconditioned by means of an incomplete LU approximate factorization. The proposed solution technique constitutes a fast and robust algorithm for treating laminar flows at high Reynolds numbers. The pressure field is obtained at a subsequent step by solving a convection– diffusion equation in terms of the stagnation pressure, which presents certain advantages compared with the widely used static pressure Poisson equation. Results are shown for a wide variety of applications including internal and external flows.     

Energy-optimal time-dependent regimes of viscous incompressible fluid flow

The hydrodynamic equations of a viscous incompressible fluid are modified for axisymmetric flows in a pipe of time-varying radius. A new exact time-dependent solution of these equations which generalizes the well-known classic steady-state Hagen–Poiseuille solution for flow in a pipe of constant radius (independent of time) is obtained. It is shown that the law of time variation in the pipe radius can be determined from the condition of the minimum work done to pump a given fluid volume through such a pipe during the radius variation cycle period. A generalization of the optimal branching pipeline in which, instead of the Poiseuille law, its modification based on the use of the exact solution corresponding to the time-dependent M-shaped regime is employed is suggested. It is shown that the hydraulic resistance can be reduced over a certain range of the parameters of the time-dependent flow regime as compared with the steady-state pipe flow regime. The conclusion obtained can be used for the development of the hydrodynamic basis for simulating the optimal hydrodynamic blood flow regime.