The Effect of Streamwise Tube Motion on the Unsteady Fluid Forces in a Normal Triangle Tube Array

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.

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.


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.


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.


Author(s):  
David Gross ◽  
Yann Roux ◽  
Benjamin Rousse ◽  
François Pétrié ◽  
Ludovic Assier ◽  
...  

The problem of Vortex-Induced Vibrations (VIV) on spool and jumper geometries is known to present several drawbacks when approached with conventional engineering tools used in the study of VIV on risers. Current recommended practices can lead to over-conservatism that the industry needs to quantify and minimize within notably cost reduction objectives. Within this purpose, the paper will present a brief critical review of the Industry standards and more particularly focus on both experimental and Computational Fluid Dynamic (CFD) approaches. Both qualitative and quantitative comparisons between basin tests and CFD results for a 2D ‘M-shape’ spool model will be detailed. The results presented here are part of a larger experimental and numerical campaign which considered a number of current velocities, heading and geometry configurations. The vibratory response of the model will be investigated for one of the current velocities and compared with the results obtained through recommended practices (e.g. Shear7 and DNV guidelines). The strategy used by the software K-FSI to solve the fluid-structure interaction (FSI) problem is a partitioned coupling solver between fluid solver (FINE™/Marine) and structural solvers (ARA). FINE™/Marine solves the Reynolds-Averaged Navier-Stokes Equations in a conservative way via the finite volume method and can work on structured or unstructured meshes with arbitrary polyhedrons, while ARA is a nonlinear finite element solver with a large displacement formulation. The experiments were conducted in the BGO FIRST facility located in La Seyne sur Mer, France. Particular attention was paid towards the model design, fabrication, instrumentation and characterization, to ensure an excellent agreement between the structural numerical model and the actual physical model. This included the use of a material with low structural damping, the performance of stiffness and decay tests in air and in still water, plus the rationalization of the instrumentation to be able to capture the response with the minimum flow perturbation or interaction due to instrumentation.


1994 ◽  
Vol 61 (3) ◽  
pp. 629-633 ◽  
Author(s):  
S. H. Smith

When a stretching surface is moved quickly, for a short period of time, a pulse is transmitted to the surrounding fluid. Here we describe an exact solution in terms of a similarity variable for the Navier-Stokes equations which represents the effect of this pulse for two-dimensional flow. The unusual feature is that this solution is only valid for a limited range of the Reynolds number; outside this domain unbounded velocities result.


Author(s):  
Jaromi´r Hora´cˇek ◽  
Miloslav Feistauer ◽  
Petr Sva´cˇek

The contribution deals with the numerical simulation of the flutter of an airfoil with three degrees of freedom (3-DOF) for rotation around an elastic axis, oscillation in the vertical direction and rotation of a flap. The finite element (FE) solution of two-dimensional (2-D) incompressible Navier-Stokes equations is coupled with a system of nonlinear ordinary differential equations describing the airfoil vibrations with large amplitudes taking into account the nonlinear mass matrix. The time-dependent computational domain and a moving grid are treated by the Arbitrary Lagrangian-Eulerian (ALE) method and a suitable stabilization of the FE discretization is applied. The developed method was successfully tested by the classical flutter computation of the critical flutter velocity using NASTRAN program considering the linear model of vibrations and the double-lattice aerodynamic theory. The method was applied to the numerical simulations of the post flutter regime in time domain showing Limit Cycle Oscillations (LCO) due to nonlinearities of the flow model and vibrations with large amplitudes. Numerical experiments were performed for the airfoil NACA 0012 respecting the effect of the air space between the flap and the main airfoil.


Author(s):  
Juan B. V. Wanderley ◽  
Gisele H. B. Souza ◽  
Carlos Levi

Numerical simulations of Vortex Induced Vibration have been failing to duplicate accurately experimental data mostly due to the complexity of the physics involved in the real problem. Therefore, a careful and comprehensive investigation on CFD algorithms is still required to indicate the most suitable numerical scheme to handle such a complicate problem. Grid generation, boundary condition implementation, and coupling between the fluid flow governing equations and body motion equation are known to have strong influence on the qualities of the numerical results. This work presents results obtained from a long-term investigation featuring different CFD methods. The investigations enabled the selection of a very effective algorithm that showed an outstanding agreement between experiment and numerical simulation of the VIV phenomenon. Good agreement is obtained in the entire range of reduced velocity covered by the experimental investigations. The successful algorithm discussed here applies the Beam and Warming implicit scheme to solve the two-dimensional slightly compressible Navier–Stokes equations with the K-ε turbulence model to simulate the turbulent flow at the wake of the cylinder.


2009 ◽  
Vol 2009 ◽  
pp. 1-10 ◽  
Author(s):  
Uzu-Kuei Hsu

This research adopts a shock tube 16 meters long and with a 9 cm bore to create a supersonic, high-temperature, and high-pressure flowfield to observe the gasification and ignition of HTPB solid fuel under different environments. Also, full-scale 3D numerical simulation is executed to enhance the comprehension of this complex phenomenon. The CFD (Computational Fluid Dynamics) code is based on the control volume method and the pre-conditioning method for solving the Navier-Stokes equations to simulate the compressible and incompressible coupling problem. In the tests, a HTPB slab is placed in the windowed-test section. Various test conditions generate different supersonic Mach numbers and environmental temperatures. In addition, the incident angles of the HTPB slab were changed relative to the incoming shock wave. Results show that as the Mach number around the slab section exceeded 1.25, the flowfield temperature achieved 1100 K, which is higher than the HTPB gasification temperature (930 K~1090 K). Then, gasification occurred and a short-period ignition could be observed. In particular, when the slab angle was7∘, the phenomenon became more visible. This is due to the flow field temperature increase when the slab angle was at7∘.


Sign in / Sign up

Export Citation Format

Share Document