Comparison of Fluid Structure Interaction Capabilities in Finite Volume and Finite Element Based Codes

Author(s):  
Qiyue Lu ◽  
Alfonso Santiago ◽  
Seid Koric ◽  
Paula Cordoba

Abstract Fluid-Structure Interaction (FSI) simulations have applications to a wide range of engineering areas. One popular technique to solve FSI problems is the Arbitrary Lagrangian-Eulerian (ALE) method. Both academic and industry communities developed codes to implement the ALE method. One of them is Alya, a Finite Element Method (FEM) based code developed in Barcelona Supercomputing Center (BSC). By analyzing the application on a simplified artery case and compared to another commercial code, which is Finite Volume Method (FVM) based, this paper discusses the mathematical background of the solver for domains, and carries out verification work on Alya’s FSI capability. The results show that while both codes provide comparable FSI results, Alya has exhibited better robustness due to its Subgrid Scale (SGS) technique for stabilization of convective term and the subsequent numerical treatments. Thus this code opens the door for more extensive use of higher fidelity finite element based FSI methods in future.

Author(s):  
Yanbo Huang ◽  
Shanshan Li ◽  
Zhenhai Pan

Abstract Fluid-structure interaction (FSI) is an important fundamental problem with wide scientific and engineering applications. The immersed boundary method has proved to be an effective way to model the interaction between a moving solid and its surrounding fluid. In this study, a novel modeling approach based on the coupled immersed-boundary and finite-volume method is proposed to simulate fluid-structure interaction problems. With this approach, the whole computational domain is treated as fluid and discretized by only one set of Eulerian grids. The computational domain is divided into solid parts and fluid parts. A goal velocity is locally determined in each cell inside the solid part. At the same time, the hydrodynamic force exerted on the solid structure is calculated by integrating along the faces between the solid cells and fluid cells. In this way, the interaction between the solid and fluid is solved explicitly and the costly information transfer between Lagranian grids and Eulerian grids is avoided. The interface is sharply restricted into one single grid width throughout the iterations. The proposed modeling approach is validated by conducting several classic numerical experiments, including flow past static and freely rotatable square cylinders, and sedimentation of an ellipsoid in finite space. Throughout the three numerical experiments, satisfying agreements with literatures have been obtained, which demonstrate that the proposed modeling approach is accurate and robust for simulating FSI problems.


2016 ◽  
Vol 307 ◽  
pp. 670-695 ◽  
Author(s):  
Vito Pasquariello ◽  
Georg Hammerl ◽  
Felix Örley ◽  
Stefan Hickel ◽  
Caroline Danowski ◽  
...  

2019 ◽  
Vol 29 (05) ◽  
pp. 967-994 ◽  
Author(s):  
Alessio Castorrini ◽  
Alessandro Corsini ◽  
Franco Rispoli ◽  
Kenji Takizawa ◽  
Tayfun E. Tezduyar

Computational fluid–structure interaction (FSI) and flow analysis now have a significant role in design and performance evaluation of turbomachinery systems, such as wind turbines, fans, and turbochargers. With increasing scope and fidelity, computational analysis can help improve the design and performance. For example, it can help add a passive morphing attachment (MA) to the blades of an axial fan for the purpose of controlling the blade load and section stall. We present a stabilized Arbitrary Lagrangian–Eulerian (ALE) method for computational FSI analysis of passive morphing in turbomachinery. The main components of the method are the Streamline-Upwind/Petrov–Galerkin (SUPG) and Pressure-Stabilizing/Petrov–Galerkin (PSPG) stabilizations in the ALE framework, mesh moving with Jacobian-based stiffening, and block-iterative FSI coupling. The turbulent-flow nature of the analysis is handled with a Reynolds-Averaged Navier–Stokes (RANS) model and SUPG/PSPG stabilization, supplemented with the “DRDJ” stabilization. As the structure moves, the fluid mechanics mesh moves with the Jacobian-based stiffening method, which reduces the deformation of the smaller elements placed near the solid surfaces. The FSI coupling between the blocks of the fully-discretized equation system representing the fluid mechanics, structural mechanics, and mesh moving equations is handled with the block-iterative coupling method. We present two-dimensional (2D) and three-dimensional (3D) computational FSI studies for an MA added to an axial-fan blade. The results from the 2D study are used in determining the spanwise length of the MA in the 3D study.


Author(s):  
Bhuiyan Shameem Mahmood Ebna Hai ◽  
Markus Bause

Advanced composite materials such as Carbon Fiber Reinforced Plastics (CFRP) are being applied to many aircraft structures in order to improve performance and reduce weight. Most composites have strong, stiff fibers in a matrix which is weaker and less stiff. However, aircraft wings can break due to Fluid-Structure Interaction (FSI) oscillations or material fatigue. This paper focuses on the analysis of a non-linear fluid-structure interaction problem and its solution in the finite element software package DOpElib: the deal.II based optimization library. The principal aim of this research is to explore and understand the behaviour of the fluid-structure interaction during the impact of a deformable material (e.g. an aircraft wing) on air. Here we briefly describe the analysis of incompressible Navier-Stokes and Elastodynamic equations in the arbitrary Lagrangian-Eulerian (ALE) frameworks in order to numerically simulate the FSI effect on a double wedge airfoil. Since analytical solutions are only available in special cases, the equation needs to be solved by numerical methods. This coupled problem is defined in a monolithic framework and fractional-step-θ time stepping scheme are implemented. Spatial discretization is based on a Galerkin finite element scheme. The non-linear system is solved by a Newton method. The implementation using the software library package DOpElib and deal.II serves for the computation of different fluid-structure configurations.


Sign in / Sign up

Export Citation Format

Share Document