Reduced-order model for laminar vortex-induced vibration of a rigid circular cylinder with an internal nonlinear absorber

The nonlinear interaction of a laminar flow and a sprung rigid circular cylinder results in vortex-induced vibration (VIV) of the cylinder. Passive suppression of the VIV by attaching an internal nonlinear vibration absorber that acts, in essence, as a nonlinear energy sink (NES) to the cylinder has been observed in finite-element computations involving thousands of degrees of freedom (DOF). A single-DOF self-excited oscillator is developed to approximate the limit-cycle oscillation (LCO) of the cylinder undergoing VIV. This self-excited oscillator models the interaction of the flow and the cylinder. Then, a two-DOF reduced-order model for the system with the internal NES is constructed by coupling the single-DOF NES to the single-DOF self-excited oscillator. Hence, the complicated high-dimensional system of flow-cylinder-NES involving thousands of DOF is reduced to a two-DOF model. The two targeted energy transfer mechanisms responsible for passive VIV suppression that are observed in the finite-element computations are fully reproduced using the two-DOF reduced-order model. This reduction of the dynamics to an easily tractable low-dimensional reduced-order model facilitates the approximate analysis of the underlying dynamics. Moreover, the underlying assumptions of the order reduction, and the parameter ranges of validity of the reduced-order model are formulated and systematically studied.     

Analysis of a conducting fluid in a thin annulus with rotating insulated walls under radial magnetic effect

This work considers an electrically conducting fluid filled between two concentric cylindrical walls relatively close to each other. A theoretical solution for the steady Taylor–Couette flow between these two electrically insulated rotating cylinders under the influence of a radial magnetic field is provided in this work. By solving the appropriate set of governing equations simultaneously, the profiles of fluid tangential velocity component and induced magnetic field were obtained as complicated functions involving the modified Bessel functions of the first and second kinds of the first-order in terms of radial coordinates and Hartmann number. A computational study was also performed to validate the present theoretical solution. The analytical and computational results are identical when Ha = 1 while these results only slightly deviate from each other as Ha increases. Current results show that, the presence of the external magnetic field causes the flow close to the slower cylinder to accelerate while that close to the faster cylinder to decelerate. This has clearly implied the fact that an external magnetic field tends to make the velocity distribution across the inner and outer cylinders more uniform.     

Effect of wall alignment in a very short rotating annulus

This paper reports numerical results of the study of effects of cylinders wall alignment in a small aspect ratio Taylor–Couette system. The investigation concerns bifurcations of steady vortical structures when the cylindrical walls defining the gap are not perfectly parallel. The imperfection is introduced by opening the outer fixed cylinder with a certain angle with regard to the vertical to form a tapered very short liquid column and keeping the inner rotating cylinder wall vertical. The numerical results obtained for the velocity components have revealed that bifurcation from a particular mode to another one occurs at a range of specific values of the inclination angle of the outer cylinder. The band width of the angle at which bifurcation occurred depended on the Reynolds number Re and was found to become narrower as Re increased. It is shown that geometrically broken symmetry can yield flow symmetry for specific combinations of geometrical and dynamical parameters.     

Streamwise oscillations of a cylinder beneath a free surface: Free surface effects on vortex formation modes

A computational study of a viscous incompressible two-fluid model with an oscillating cylinder is investigated at a Reynolds number of 200 and at a dimensionless displacement amplitude of A=0.13 and for the dimensionless forcing cylinder oscillation frequency-to-natural vortex shedding frequency ratios, f/f₀=1.5,2.5,3.5. Specifically, two-dimensional flow past a circular cylinder subject to forced in-line oscillations beneath a free surface is considered. The method is based on a finite volume discretization of the two-dimensional continuity and unsteady Navier-Stokes equations (when a solid body is present) on a fixed Cartesian grid. Two-fluid model based on improved volume-of-fluid method is used to discretize the free surface interface. The study focuses on the laminar asymmetric flow structure in the near wake region and lock-on phenomena at a Froude number of 0.2 and for the dimensionless cylinder submergence depths, h=0.25, 0.5 and 0.75. The equivorticity patterns and pressure distribution contours are used for the numerical flow visualization. The code validations in special cases show good comparisons with previous numerical results.     

Expansions at small Reynolds numbers for the flow past a porous circular cylinder

The purpose of this article is to use the method of matched asymptotic expansions (MMAE) in order to study the two-dimensional steady low Reynolds number flow of a viscous incompressible fluid past a porous circular cylinder. We assume that the flow inside the porous body is described by the continuity and Brinkman equations, and the velocity and boundary traction fields are continuous across the interface between the fluid and porous media. Formal expansions for the corresponding stream functions are used. We show that the force exerted by the exterior flow on the porous cylinder admits an asymptotic expansion with respect to low Reynolds numbers, whose terms depend on the characteristics of the porous cylinder. In addition, by considering Darcy's law for the flow inside the porous circular cylinder, an asymptotic formula for the force on the cylinder is obtained. Also, a porous circular cylinder with a rigid core inside is considered with Brinkman equation inside the porous region. Stress jump condition is used at the porous–liquid interface together with the continuity of velocity components and continuity of normal stress. Some particular cases, which refer to the low Reynolds number flow past a solid circular cylinder, have also been investigated.     

An automated curve fairing algorithm for cubic B-spline curves

The fairing of contours is an important part of the computerised production of curved objects. A number of different fairing strategies have been proposed. In a recent paper we have introduced an extension of Kjellander's algorithm for fairing parametric B-splines, which can be applied to a wide range of two- and three-dimensional curves. In this paper we describe developments towards a fully automated fairing procedure based on our new algorithm. Like that of Kjellander, our algorithm fairs by an iterative process. The key problems are to decide which points need to be faired and how many times to iterate. Sapidis (1992) has proposed a curve fairness indicator to locate points to be faired and a criterion for termination of fairing. However, we have found that for interpolating curves with great variation in curvature Sapidis' criterion tends to concentrate fairing on regions with large curvature. Therefore we have developed a new scale-independent curve fairness indicator which does not suffer from this drawback. A number of examples of faired curves are presented.     

气-液横向流动下悬臂柱体结构涡激振动机理研究北大核心CSCD

针对潜艇侦查望远镜举升水面时产生的涡激振动现象,该文建立了悬臂柱体结构受气-液两种不同横向流作用下的涡激振动理论模型.研究了两种流体不同的分布比和密度比对柱体结构涡激振动行为的影响规律.基于Galerkin法和Runge-Kutta法,得到了柱体结构涡激振动响应的数值结果.研究表明,柱体结构的涡激振动锁频区随着流体分布比的增大而增大,自由端最大幅值随着流体分布比的增大先增大后减小.当分布比为0.5附近时,振幅出现极大值,该极大值随着流体密度比的减小呈现明显的增大趋势.此外,柱体的动力学行为随着流体分布比的变化呈现出周期和多周期等振动模式.该研究可为潜艇侦查望远镜结构的设计与分析提供理论指导意义.     

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode’“dual vortex pair mode” downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f _m in addition to the vortex shedding frequency f _s is identified. Frequency coupling between f _m and f _s occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.     

On the spectrum for the artificial compressible system

Stability of stationary solutions of the incompressible Navier–Stokes system and the corresponding artificial compressible system is considered. Both systems have the same sets of stationary solutions and the incompressible system is obtained from the artificial compressible one in the zero limit of the artificial Mach number ? which is a singular limit. It is proved that if a stationary solution of the incompressible system is asymptotically stable and the velocity field of the stationary solution satisfies an energy-type stability criterion by variational method with admissible functions being only potential flow parts of velocity fields, then it is also stable as a solution of the artificial compressible one for sufficiently small ?. The result is applied to the Taylor problem.     

hp-Finite element for free vibration analysis of the orthotropic triangular and rectangular plates

The hp-version of the finite element method based on a triangular p-element is applied to free vibration of the orthotropic triangular and rectangular plates. The element's hierarchical shape functions, expressed in terms of shifted Legendre orthogonal polynomials, is developed for orthotropic plate analysis by taking into account shear deformation, rotary inertia, and other kinematics effects. Numerical results of frequency calculations are found for the free vibration of the orthotropic triangular and rectangular plates with the effect of the fiber orientation and plate boundary conditions. The results are very well compared to those presented in the literature.     

Steady-state flow of an incompressible viscoelastic polymer fluid between two coaxial cylinders

A boundary value problem for a quasi-linear equation determining the velocity profile of a flow of a polymer fluid in a pipe formed by two coaxial cylinders is considered. On the basis of methods of approximation without saturation, a computational algorithm of increased accuracy is developed, making it possible to solve the problem in a wide range of parameters, including record-low values of r ₀, the radius of the inner cylinder.     

A new triangular spectral element method I: implementation and analysis on a triangle

This paper serves as our first effort to develop a new triangular spectral element method (TSEM) on unstructured meshes, using the rectangle–triangle mapping proposed in the conference note (Li et al. 2011). Here, we provide some new insights into the originality and distinctive features of the mapping, and show that this transform only induces a logarithmic singularity, which allows us to devise a fast, stable and accurate numerical algorithm for its removal. Consequently, any triangular element can be treated as efficiently as a quadrilateral element, which affords a great flexibility in handling complex computational domains. Benefited from the fact that the image of the mapping includes the polynomial space as a subset, we are able to obtain optimal L ²- and H ¹-estimates of approximation by the proposed basis functions on triangle. The implementation details and some numerical examples are provided to validate the efficiency and accuracy of the proposed method. All these will pave the way for developing an unstructured TSEM based on, e.g., the hybridizable discontinuous Galerkin formulation.     

The scattering of electromagnetic wave at oblique incidence in a homogeneous chiral medium

We consider the scattering of a time-harmonic electromagnetic wave by a perfectly and imperfectly conducting infinite cylinder at oblique incidence respectively. We assume that the cylinder is embedded in a homogeneous chiral medium and the cylinder is parallel to the z axis. Since the x components and y components of electric field and magnetic field can be expressed in terms of their z components, we can derive from Maxwell's equations and corresponding boundary conditions that the scattering problem is modeled as a boundary value problem for the z components of electric field and magnetic field. By using Rellich's lemma and variational approach, the uniqueness and the existence of solutions are justified.     

A finite element method for the axisymmetric three‐field Stokes system with a discontinuous pressure

A quadrilateral based velocity‐pressure‐extrastress tensor mixed finite element method for solving the three‐field Stokes system in the axisymmetric case is studied. The method derived from Fortin's Q₂ − P₁ velocity‐pressure element is to be used in connection with the standard Galerkin formulation. This makes it particularly suitable for the numerical simulation of viscoelastic flow. It is proven to be second‐order convergent in the natural weighted Sobolev norms, for the system under consideration. The crucial result that the method is uniformly stable is proven for the case of rectangular meshes. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 739–763, 1999     

Simulation and study of stratified flows around finite bodies

The flows past a sphere and a square cylinder of diameter d moving horizontally at the velocity U in a linearly density-stratified viscous incompressible fluid are studied. The flows are described by the Navier–Stokes equations in the Boussinesq approximation. Variations in the spatial vortex structure of the flows are analyzed in detail in a wide range of dimensionless parameters (such as the Reynolds number Re = Ud/ν and the internal Froude number Fr = U/(Nd), where ν is the kinematic viscosity and N is the buoyancy frequency) by applying mathematical simulation (on supercomputers of Joint Supercomputer Center of the Russian Academy of Sciences) and three-dimensional flow visualization. At 0.005 < Fr < 100, the classification of flow regimes for the sphere (for 1 < Re < 500) and for the cylinder (for 1 < Re < 200) is improved. At Fr = 0 (i.e., at U = 0), the problem of diffusion-induced flow past a sphere leading to the formation of horizontal density layers near the sphere’s upper and lower poles is considered. At Fr = 0.1 and Re = 50, the formation of a steady flow past a square cylinder with wavy hanging density layers in the wake is studied in detail.     

Extracting energy from Vortex-Induced Vibrations: A parametric study

Here, Vortex-Induced Vibrations (VIVs) of a circular cylinder are analyzed as a potential source for energy harvesting. To this end, VIV is described by a one-degree-of-freedom model where fluid forces are introduced from experimental data from forced vibration tests. The influence of some influencing parameters, like the mass ratio m^∗ or the mechanical damping ζ in the energy conversion factor is investigated. The analysis reveals that: (i) the maximum efficiency η_M is principally influenced by the mass-damping parameter m^∗ζ and there is an optimum value of m^∗ζ where η_M presents a maximum; (ii) the range of reduced velocities with significant efficiency is mainly governed by m^∗, and (iii) it seems that encouraging high efficiency values can be achieved for high Reynolds numbers.     

The periodic motions of a non-autonomous Hamiltonian system with two degrees of freedom at parametric resonance of the fundamental type

The motion of an almost autonomous Hamiltonian system with two degrees of freedom, 2π-periodic in time, is considered. It is assumed that the origin is an equilibrium position of the system, the linearized unperturbed system is stable, and its characteristic exponents ±iωj (j = 1,2) are pure imaginary. In addition, it is assumed that the number 2ω₁ is approximately an integer, that is, the system exhibits parametric resonance of the fundamental type. Using Poincaré's theory of periodic motion and KAM-theory, it is shown that 4π-periodic motions of the system exist in a fairly small neighbourhood of the origin, and their bifurcation and stability are investigated. As applications, periodic motions are constructed in cases of parametric resonance of the fundamental type in the following problems: the plane elliptical restricted three-body problem near triangular libration points, and the problem of the motion of a dynamically symmetrical artificial satellite near its cylindrical precession in an elliptical orbit of small eccentricity.     

Stabilization of traffic flow in optimal velocity model via delayed-feedback control

Traffic jams may occur due to various reasons, such as traffic accidents, lane reductions and on-ramps. In order to suppress the traffic congestion in an optimal velocity traffic model without any driver’s delay taken into account, a delayed-feedback control of both displacement and velocity differences is proposed in this study. By using the delay-independent stability criteria and the H_∞-norm, the delayed-feedback control can be determined to stabilize the unstable traffic flow and suppress the traffic jam. The numerical case studies are given to demonstrate and verify the new control method. Furthermore, a comparison is made between the new control method and the method proposed by Konishi et al. [K. Konishi, M. Hirai, H. Kokame, Decentralized delayed-feedback control of an optimal velocity traffic model, Eur. Phys. J. B 15 (2000) 715–722]. The results show that the new control method makes the traffic flow more stable and improves the control performance.     

On the accuracy of the viscous form in simulations of incompressible flow problems

We present a numerical study of drag/lift and flux estimates using two forms of Navier‐Stokes equations (NSE) that are equivalent in the continuum formulation but not in the discrete finite element formulation. The two investigated forms of the NSE differ in the viscous term that is represented in one form by νΔ u with ν being the viscosity and 2ν?·?^S u in the other form where ?^S represents the deformation tensor. The study consists of numerical analysis of the two forms and computations of drag/lift, pressure drop on the cylinder problem and computations of flux for the Poiseuille flow. The main objective is to provide a clear comparison of the reference values for the maximal drag and lift coefficient at the cylinder and for the pressure difference between the front and the back of the cylinder at the final time for the two forms of NSEs. Our computational results of the reference values do not differ significantly between the two forms, but the differences are there. For the Poiseuille flow, the differences in the flux computations were much smaller, and this agreed with the computationally obtained results of the divergence of the velocity field. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 523–541, 2012     

The effect of actuator parameters on the critical flow velocity of a fluidic amplifier

The effect of actuator parameters on the critical flow velocity of the fluidic amplifier in liquid-jet hammers has been investigated numerically and experimentally. The flow in the fluidic amplifier and actuators coupled with the rigid body movement of the impacting body has been simulated using a commercial CFD software package, Fluent. The flow is modeled by the RNG-based κ–ε turbulence model and the incompressible Navier–Stokes equations. Dynamic layering method and a user-defined function written in C programming language are used to update the mesh in the simulations. The results show that, increasing the piston diameter decreases rapidly the critical flow velocity as the piston diameter is less than a certain value. The critical flow velocity increases sharply as the piston rod diameter is greater than a certain value, and increases nearly linearly with mass of the impacting body, and is independent on stroke length of the impacting body. The findings of the numerical investigations agree well with corresponding experimental results.