Fully-Coupled Fluid Solid Computation for Simulation of Flutter in Tubes and Tube Arrays

Author(s):  
E. Longatte

This work is concerned with the modelling of the interaction of a fluid with a rigid or a flexible elastic cylinder in presence of axial or cross-flow. A partitioned procedure is involved to performe 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 descrete or beam elements in the linear elasticity framework and systems are solved through a finite element method. The resulting strongly coupled fluid solid set of non linear equations is solved by means of a partitioned solution procedure. A fixed point method combined with under-relaxation is involved to ensure the optimal convergence of the iterative procedure. In the present work two examples are presented to show the methodology robustness and efficiency. The purpose is to attempt to simulate a fluid structure interaction resulting in the development of a dynamic instability induced by a positive damping generation of the system. Both flutter of a flexible cylinder conveying an internal fluid and fluid-elastic instability of a tube array submitted to an external cross flow are investigated numerically. According to first results the partitioned procedure relies on consistant numerical methods ensuring energy conservation at the interface and describing with a sufficient accuracy the mechanical energy transfer between fluid and solid systems through the interface with a limited numerical diffusion. Therefore it seems to be qualitatively convenient for simulation of flutter. For a quantitative evaluation of the methodology further complementary simulations validating these developments from a physical point of view will be required in order to confirm these first trends.

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.


2000 ◽  
Vol 123 (2) ◽  
pp. 435-438 ◽  
Author(s):  
A. Beskok ◽  
A. R. Srinivasa

The operation of a micro-pump system driven by a magnetoelastic polymeric membrane developed at Texas A&M University is analyzed by numerical simulations. Unsteady, incompressible Navier-Stokes equations in a moving boundary system are solved by a spectral element methodology, employing an Arbitrary Lagrangian Eulerian (ALE) formulation on unstructured meshes. The performance of the micro-pump is evaluated as a function of the Reynolds number and the geometric parameters. The volumetric flowrate is shown to increase as a function of the Reynolds number. The system is simulated by assuming the deformation of the membrane. The required voltage and current are then calculated by a lumped parameter analysis.


Author(s):  
Berend van Wachem

In recent years, the application of coupled solver techniques to solve the Navier-Stokes equations has become increasingly popular. The main reason for this, is the increased robustness originating from the implicit and global treatment of the pressure-velocity coupling. The drawback of a coupled solver are the increase in memory requirement and the increased complexity of implementation. However, fully coupled methods are reported to have an overall favorable computational cost when a suitable pre-conditioner and algorithm for solving the resulting set of linear equations are employed. In solving multi-phase flow problems, the coupled solver approach is even more advantageous than in single-phase, due to the presence of large source terms arising from the coupling of the phases. In this chapter, various strategies for the fully coupled approach are discussed. These strategies include employing artificial compressibility, applying physically consistent cell face interpolation, and applying momentum weighted cell face interpolation. The idea behind the strategies is outlined and their advantages and disadvantages are discussed. The treatment of source terms and volume fraction in coupled methods is also shown. Finally, a number of examples of implementations and calculations are presented.


2020 ◽  
Vol 8 (6) ◽  
pp. 387
Author(s):  
Farid Piran ◽  
Hassan Karampour ◽  
Peter Woodfield

Vortex-induced vibrations (VIV) of hexagonal cylinders at Reynolds number of 1000 and mass ratio of 2 are studied numerically. In the numerical model, the Navier–Stokes equations are solved using finite volume method, and the fluid-structure interaction (FSI) is modelled using Arbitrary Lagrangian Eulerian (ALE) Scheme. The numerical model accounts for the cross-flow vibration of the cylinders, and is validated against published experimental and numerical results. In order to account for different angles of attack, the hexagonal cylinders are studied in the corner and face orientations. The results are compared with the published results of circular and square cylinders. Current results show that within the studied range of reduced velocities (up to 20), unlike circular and square cylinders, no lock-in response is observed in the hexagonal cylinders. The maximum normalized VIV amplitudes of the hexagonal cylinders are 0.45, and are significantly lower than those of circular and square cylinders. Vortex shedding regimes of the corner-oriented hexagons are mostly irregular. However, in the face-oriented hexagons, the shedding modes are more similar to the typical P + S and 2P modes.


Author(s):  
José L. D. Alves ◽  
Carlos E. Silva ◽  
Nestor O. Guevara ◽  
Alvaro L. G. A. Coutinho ◽  
Renato N. Elias ◽  
...  

This work presents the development of EdgeCFD-ALE, a finite element system for complex fluid-structure interactions designed for offshore hydrodynamics. Sloshing of liquids in tanks, wave breaking in ships, offshore platforms motions and green water on decks are important examples of these problems. The software uses edge-based parallel stabilized finite elements for the Navier-Stokes equations and the Volume-Of-Fluid method for the free-surface, both described by an Arbitrary Lagrangian Eulerian (ALE) formulation. Turbulence in is treated by a Smagorinsky model. Mesh updating is accomplished by a parallel edge-based solution of a non-homogeneous scalar diffusion problem in each spatial coordinate. Boundary conditions involve the motion of the immersed body’s surface, i.e., the fluid-structure interface, taken as the Lagrangian portion of the domain in the overall problem. The simulation capabilities of the present software are demonstrated in the solution of two problems, the interaction of two cylinders in tandem and the free fall of a sphere.


2016 ◽  
Vol 7 (2) ◽  
Author(s):  
Žarko Milošević ◽  
Dalibor Nikolić ◽  
Igor Saveljić ◽  
Velibor Isailović ◽  
Thanos Bibas ◽  
...  

Benign paroxysmal positional vertigo (BPPV) is the most common type of vertigo. The symptoms of BPPV typically appear after angular movements of the head. BPPV leads to dizziness, nausea and imbalance. In this study, we examined a model of the semi-circular canal (SCC) with fully 3D three dimensional anatomical data from specific patient. A full Navier-Stokes equations and continuity equations are used for fluid domain with Arbitrary-Lagrangian Eulerian (ALE) formulation for mesh motion of finite element. Fluid-structure interaction for fluid coupling with cupula deformation is used. Particle tracking algorithm is implemented for particle motion. Motion of the otoconia particles which is main cause for BPPV is simulated. Velocity distribution, shear stress and force from endolymph side are presented for patient specific three SCC. We compared our numerical models with experiments with head moving and nystagmus eye tracking. Numerical simulation can give more details and understanding of the pathology of the specific patient in standard clinical diagnostic and therapy procedure for BPPV.


2014 ◽  
Vol 945-949 ◽  
pp. 642-645
Author(s):  
Ren Qiang Xi ◽  
Xue Dong Jiang ◽  
Yun Song He

This work is concerned with the modelling of the interaction of fluid flow with flexible solid structures. The fluid flow considered is governed by the incompressible Navier-Stokes equations and modelled with stabilised low order velocity-pressure finite elements. The governing equation of fluid movement is described by an arbitrary Lagrangian-Eulerian (ALE) strategy. The structure is represented by means of an appropriate standard finite element formulation. A simple data transfer strategy based on a finite element type interpolation of the interface degrees of freedom guarantees kinematic consistency and equilibrium of the stresses along the interface. The resulting strongly coupled set of non-linear equations is solved by means of a partitioned solution procedure, which is based on the Newton-Raphson methodology and incorporates the full linearization of the overall incremental problem.


2016 ◽  
Vol 44 (3) ◽  
pp. 174-190 ◽  
Author(s):  
Mario A. Garcia ◽  
Michael Kaliske ◽  
Jin Wang ◽  
Grama Bhashyam

ABSTRACT Rolling contact is an important aspect in tire design, and reliable numerical simulations are required in order to improve the tire layout, performance, and safety. This includes the consideration of as many significant characteristics of the materials as possible. An example is found in the nonlinear and inelastic properties of the rubber compounds. For numerical simulations of tires, steady state rolling is an efficient alternative to standard transient analyses, and this work makes use of an Arbitrary Lagrangian Eulerian (ALE) formulation for the computation of the inertia contribution. Since the reference configuration is neither attached to the material nor fixed in space, handling history variables of inelastic materials becomes a complex task. A standard viscoelastic material approach is implemented. In the inelastic steady state rolling case, one location in the cross-section depends on all material locations on its circumferential ring. A consistent linearization is formulated taking into account the contribution of all finite elements connected in the hoop direction. As an outcome of this approach, the number of nonzero values in the general stiffness matrix increases, producing a more populated matrix that has to be solved. This implementation is done in the commercial finite element code ANSYS. Numerical results confirm the reliability and capabilities of the linearization for the steady state viscoelastic material formulation. A discussion on the results obtained, important remarks, and an outlook on further research conclude this work.


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.


2017 ◽  
Vol 24 (1) ◽  
pp. 152-167
Author(s):  
Izhak Bucher ◽  
Ran Gabai ◽  
Harel Plat ◽  
Amit Dolev ◽  
Eyal Setter

Vibrations are often represented as a sum of standing waves in space, i.e. normal modes of vibration. While this can be mathematically accurate, the representation as travelling waves can be compact and more appropriate from a physical point of view, in particular when the energy flux along the structure is meaningful. The quantification of travelling waves assists in computing the energy being transferred and propagated along a structure. It can provide local as well as global information about the structure through which the mechanical energy flows. Presented in this paper is a new method to quantify the fraction of mechanical power being transmitted in a vibration cycle at a specific direction in space using measured data. It is shown that the method can detect local defects causing slight non-uniformity of the energy flux. Equivalence is being made with the electrical power factor and electromagnetic standing waves ratio, commonly employed in such cases. Other methods to perform experiment based wave identification in one-dimension are compared with the power flow based identification. Including a signal processing approach that fits an ellipse to the complex amplitude curve and Hilbert transform for obtaining the local phase and amplitude. A new representation of the active and reactive power flow is developed and its relationship to standing waves ratio is demonstrated analytically and experimentally.


Sign in / Sign up

Export Citation Format

Share Document