Numerical Estimation of Fluidelastic Instability in Tube Arrays

2010 ◽  
Vol 132 (4) ◽  
Author(s):  
Marwan Hassan ◽  
Andrew Gerber ◽  
Hossin Omar

This study investigates unsteady flow in tube bundles and fluid forces, which can lead to large tube vibration amplitudes, in particular, amplitudes associated with fluidelastic instability (FEI). The fluidelastic forces are approximated by the coupling of the unsteady flow model (UFM) with computational fluid dynamics (CFD). The CFD model employed solves the Reynolds averaged Navier–Stokes equations for unsteady turbulent flow and is cast in an arbitrary Lagrangian–Eulerian form to handle any motion associated with tubes. The CFD solution provides time domain forces that are used to calculate added damping and stiffness coefficients employed with the UFM. The investigation demonstrates that the UFM utilized in conjunction with CFD is a viable approach for calculating the stability map for a given tube array. The FEI was predicted for in-line square and normal triangle tube arrays over a mass damping parameter range of 0.1– 100. The effect of the P/d ratio and the Reynolds number on the FEI threshold was also investigated.

Author(s):  
H. Omar ◽  
M. Hassan ◽  
A. Gerber

This study investigates the unsteady flow and the resulting fluidelastic forces in a tube bundle. Numerical simulations are presented for normal triangle tube arrays with pitch-to-diameter (P/d) ratios of 1.35, 1.75, and 2.5 utilizing a 2-dimensional model. In this model a single tube was forced to oscillate within an otherwise rigid array. Fluid forces acting on the oscillating tube and the surrounding tubes were estimated. The predicted forces were utilized to calculate fluid force coefficients for all tubes. The numerical model solves the Reynolds-Average Navier-Stokes (RANS) equations for unsteady turbulent flow, and is cast in an Arbitrary Lagrangian-Eulerian (ALE) form to handle mesh the motion associated with a moving boundary. The fluidelastic instability (FEI) was predicted for both single and fully flexible tube arrays over a mass damping parameter (MDP) range of 0.1 to 200. The effect of the P/d ratio and the Reynolds number on the FEI threshold was investigated in this work.


1991 ◽  
Vol 113 (3) ◽  
pp. 409-417 ◽  
Author(s):  
A. El Chebair ◽  
A. K. Misra

This paper investigates theoretically for the first time the dynamical behavior and stability of a simply supported shell located coaxially in a rigid cylindrical conduit. The fluid flow is incompressible and the fluid forces consist of two parts: (i) steady viscous forces which represent the effects of upstream pressurization of the flow; (ii) unsteady forces which could be inviscid or viscous. The inviscid forces were derived by linearized potential flow theory, while the viscous ones were derived by means of the Navier-Stokes equations. Shell motion is described by the modified Flu¨gge’s shell equations. The Fourier transform technique is employed to formulate the problem. First, the system is subjected only to the unsteady inviscid forces. It is found that increasing either the internal or the annular flow velocity induces buckling, followed by coupled mode flutter. When both steady viscous and unsteady inviscid forces are applied, for internal flow, the system becomes stabilized; while for annular flow, the system loses stability at much lower velocities. Second, the system is only subjected to the unsteady viscous forces. Calculations are only performed for the internal flow case. The results are compared to those of inviscid theory. It is found that the effects of unsteady viscous forces on the stability of the system are very close to those of unsteady inviscid forces.


Author(s):  
Salim El Bouzidi ◽  
Marwan Hassan

Fluidelastic instability is generally regarded as the most severe type of flow excitation mechanism. When this mechanism prevails, it could cause serious damage to tube arrays in a very short period of time. This mechanism is characterized by a critical flow velocity beyond which the tubes undergo unstable oscillations. Recently, a number of experimental investigations showed that it is possible to have instability in the streamwise direction; previously, it was believed that fluidelastic instability was only a concern in the direction transverse to the flow. The purpose of this study is to characterize the flow in the channels surrounding a vibrating tube in a normal triangular bundle with P/d = 1.5. The tube is oscillating in the streamwise direction with a constant amplitude. Numerical simulations were conducted by solving the unsteady Reynolds Averaged Navier-Stokes equations (uRANS) cast in Arbitrary Lagrangian-Eulerian (ALE) form. The unsteady flow perturbation is estimated along the flow channel. The pressure perturbation is used to compute the streamwise unsteady force coefficients in the context of Chen’s model. The perturbation phase and decay are extracted and utilized in the framework of the Lever & Weaver model to study the stability of tube bundles due to tube motion in the streamwise direction.


2021 ◽  
Vol 11 (11) ◽  
pp. 4748
Author(s):  
Monika Balázsová ◽  
Miloslav Feistauer ◽  
Jaromír Horáček ◽  
Adam Kosík

This study deals with the development of an accurate, efficient and robust method for the numerical solution of the interaction of compressible flow and nonlinear dynamic elasticity. This problem requires the reliable solution of flow in time-dependent domains and the solution of deformations of elastic bodies formed by several materials with complicated geometry depending on time. In this paper, the fluid–structure interaction (FSI) problem is solved numerically by the space-time discontinuous Galerkin method (STDGM). In the case of compressible flow, we use the compressible Navier–Stokes equations formulated by the arbitrary Lagrangian–Eulerian (ALE) method. The elasticity problem uses the non-stationary formulation of the dynamic system using the St. Venant–Kirchhoff and neo-Hookean models. The STDGM for the nonlinear elasticity is tested on the Hron–Turek benchmark. The main novelty of the study is the numerical simulation of the nonlinear vocal fold vibrations excited by the compressible airflow coming from the trachea to the simplified model of the vocal tract. The computations show that the nonlinear elasticity model of the vocal folds is needed in order to obtain substantially higher accuracy of the computed vocal folds deformation than for the linear elasticity model. Moreover, the numerical simulations showed that the differences between the two considered nonlinear material models are very small.


2013 ◽  
Vol 136 (2) ◽  
Author(s):  
Graham Ashcroft ◽  
Christian Frey ◽  
Kathrin Heitkamp ◽  
Christian Weckmüller

This is the first part of a series of two papers on unsteady computational fluid dynamics (CFD) methods for the numerical simulation of aerodynamic noise generation and propagation. In this part, the stability, accuracy, and efficiency of implicit Runge–Kutta schemes for the temporal integration of the compressible Navier–Stokes equations are investigated in the context of a CFD code for turbomachinery applications. Using two model academic problems, the properties of two explicit first stage, singly diagonally implicit Runge–Kutta (ESDIRK) schemes of second- and third-order accuracy are quantified and compared with more conventional second-order multistep methods. Finally, to assess the ESDIRK schemes in the context of an industrially relevant configuration, the schemes are applied to predict the tonal noise generation and transmission in a modern high bypass ratio fan stage and comparisons with the corresponding experimental data are provided.


1976 ◽  
Vol 78 (2) ◽  
pp. 355-383 ◽  
Author(s):  
H. Fasel

The stability of incompressible boundary-layer flows on a semi-infinite flat plate and the growth of disturbances in such flows are investigated by numerical integration of the complete Navier–;Stokes equations for laminar two-dimensional flows. Forced time-dependent disturbances are introduced into the flow field and the reaction of the flow to such disturbances is studied by directly solving the Navier–Stokes equations using a finite-difference method. An implicit finitedifference scheme was developed for the calculation of the extremely unsteady flow fields which arose from the forced time-dependent disturbances. The problem of the numerical stability of the method called for special attention in order to avoid possible distortions of the results caused by the interaction of unstable numerical oscillations with physically meaningful perturbations. A demonstration of the suitability of the numerical method for the investigation of stability and the initial growth of disturbances is presented for small periodic perturbations. For this particular case the numerical results can be compared with linear stability theory and experimental measurements. In this paper a number of numerical calculations for small periodic disturbances are discussed in detail. The results are generally in fairly close agreement with linear stability theory or experimental measurements.


2000 ◽  
Author(s):  
B. V. Rathish Kumar ◽  
T. Yamaguchi ◽  
H. Liu ◽  
R. Himeno

Abstract Unsteady flow dynamics in a doubly constricted vessel is analyzed by using a time accurate Finite Volume solution of three dimensional incompressible Navier-Stokes equations. Computational experiments are carried out for various values of Reynolds number in order to assess the criticality of multiple mild constrictions in series and also to bring out the subtle 3D features like vortex formation. Studies reveal that pressure drop across a series of mild constrictions can get physiologically critical. Further this pressure drop is found to be sensitive to the spacing between the constrictions and also to the oscillatory nature of the inflow profile.


Author(s):  
Elisabeth Longatte

This work is concerned with the modelling of the interaction of a fluid with a rigid or a flexible elastic cylinder in the presence of axial or cross-flow. A partitioned procedure is involved to perform the computation of the fully-coupled fluid solid system. The fluid flow is governed by the incompressible Navier-Stokes equations and modeled by using a fractional step scheme combined with a co-located finite volume method for space discretisation. The motion of the fluid domain is accounted for by a moving mesh strategy through an Arbitrary Lagrangian-Eulerian (ALE) formulation. Solid dyncamics is modeled by a finite element method in the linear elasticity framework and a fixed point method is used for the fluid solid system computation. In the present work two examples are presented to show the method robustness and efficiency.


2013 ◽  
Vol 2013 ◽  
pp. 1-6
Author(s):  
Wen-Juan Wang ◽  
Yan Jia

We study the stability issue of the generalized 3D Navier-Stokes equations. It is shown that if the weak solutionuof the Navier-Stokes equations lies in the regular class∇u∈Lp(0,∞;Bq,∞0(ℝ3)),(2α/p)+(3/q)=2α,2<q<∞,0<α<1, then every weak solutionv(x,t)of the perturbed system converges asymptotically tou(x,t)asvt-utL2→0,t→∞.


Author(s):  
Fayçal Hammami ◽  
Nader Ben Cheikh ◽  
Brahim Ben Beya

This paper deals with the numerical study of bifurcations in a two-sided lid driven cavity flow. The flow is generated by moving the upper wall to the right while moving the left wall downwards. Numerical simulations are performed by solving the unsteady two dimensional Navier-Stokes equations using the finite volume method and multigrid acceleration. In this problem, the ratio of the height to the width of the cavity are ranged from H/L = 0.25 to 1.5. The code for this cavity is presented using rectangular cavity with the grids 144 × 36, 144 × 72, 144 × 104, 144 × 136, 144 × 176 and 144 × 216. Numerous comparisons with the results available in the literature are given. Very good agreements are found between current numerical results and published numerical results. Various velocity ratios ranged in 0.01≤ α ≤ 0.99 at a fixed aspect ratios (A = 0.5, 0.75, 1.25 and 1.5) were considered. It is observed that the transition to the unsteady regime follows the classical scheme of a Hopf bifurcation. The stability analysis depending on the aspect ratio, velocity ratios α and the Reynolds number when transition phenomenon occurs is considered in this paper.


Sign in / Sign up

Export Citation Format

Share Document