Numerical study of flexible flapping wings with an immersed boundary method: Fluid–structure–acoustics interaction

2019 ◽  
Vol 90 ◽  
pp. 396-409 ◽  
Author(s):  
Li Wang ◽  
Fang-Bao Tian
Author(s):  
Ranjith Maniyeri ◽  
Yong Kweon Suh ◽  
Sangmo Kang

Numerical modeling of fluid-structure interaction problems are challenging in the field of computational fluid dynamics because of the complex geometries involved and freely moving boundaries. Flapping of an inextensible filament in a uniform fluid flow is such a problem which mimics the swimming of energy harvesting eel fish. Recently, immersed boundary method has found much attention in simulating fluid-structure interaction problems due to its easiness in grid generation and memory and CPU savings. In the present work, we employed an improved version of immersed boundary method proposed by Shin et al. [1] which combines the feedback forcing scheme of the virtual boundary method with Peskin’s regularized delta function approach. A FORTRAN code is developed for the simulation of flexible filament flapping in a uniform fluid flow. The code is validated for the bench mark problem of two-dimensional flow over a circular cylinder. A single filament hanging under gravitational force is simulated using the developed code which is analogous to a rope pendulum and the results are compared with available analytical results. The results are found to be in good agreement. Finally, the interaction of the flapping filament in the uniform fluid flow is studied for different flow and structure parameters. The production of a series of vortex procession obtained in the case of flapping of filament is in good agreement with the previous available experimental and numerical results.


2011 ◽  
Vol 21 (12) ◽  
pp. 2523-2550 ◽  
Author(s):  
DANIELE BOFFI ◽  
NICOLA CAVALLINI ◽  
LUCIA GASTALDI

The Immersed Boundary Method (IBM) has been designed by Peskin for the modeling and the numerical approximation of fluid-structure interaction problems, where flexible structures are immersed in a fluid. In this approach, the Navier–Stokes equations are considered everywhere and the presence of the structure is taken into account by means of a source term which depends on the unknown position of the structure. These equations are coupled with the condition that the structure moves at the same velocity of the underlying fluid. Recently, a finite element version of the IBM has been developed, which offers interesting features for both the analysis of the problem under consideration and the robustness and flexibility of the numerical scheme. Initially, we considered structure and fluid with the same density, as it often happens when dealing with biological tissues. Here we study the case of a structure which can have a density higher than that of the fluid. The higher density of the structure is taken into account as an excess of Lagrangian mass located along the structure, and can be dealt with in a variational way in the finite element approach. The numerical procedure to compute the solution is based on a semi-implicit scheme. In fluid-structure simulations, nonimplicit schemes often produce instabilities when the density of the structure is close to that of the fluid. This is not the case for the IBM approach. In fact, we show that the scheme enjoys the same stability properties as in the case of equal densities.


2020 ◽  
Vol 142 (5) ◽  
Author(s):  
Mithun Kanchan ◽  
Ranjith Maniyeri

Abstract Many microfluidics-based applications involve fluid–structure interaction (FSI) of flexible membranes. Thin flexible membranes are now being widely used for mixing enhancement, particle segregation, flowrate control, drug delivery, etc. The FSI simulations related to these applications are challenging to numerically implement. In this direction, techniques like immersed boundary method (IBM) have been successful. In this study, two-dimensional numerical simulation of flexible membrane fixed at two end points in a rectangular channel subjected to uniform fluid flow is carried out at low Reynolds number using a finite volume based IBM. A staggered Cartesian grid system is used and SIMPLE algorithm is used to solve the governing continuity and Navier–Stokes equations. The developed model is validated using the previous research work and numerical simulations are carried out for different parametric test cases. Different membrane mode shapes are observed due to the complex interplay between the hydrodynamics and structural elastic forces. Since the membrane undergoes deformation with respect to inlet fluid conditions, a variation in flowrate past the flexible structure is confirmed. It is found that, by changing the membrane length, bending rigidity, and its initial position in the channel, flowrate can be controlled. Also, for membranes that are placed at the channel midplane undergoing self-excited oscillations, there exists a critical dimensionless membrane length condition L ≥ 1.0 that governs this behavior. Finally, an artificial neural network (ANN) model is developed that successfully predicts flowrate in the channel for different membrane parameters.


Author(s):  
Long He ◽  
Keyur Joshi ◽  
Danesh Tafti

In this work, we present an approach for solving fluid structure interaction problems by combining a non-linear structure solver with an incompressible fluid solver using immersed boundary method. The implementation of the sharp-interface immersed boundary method with the fluid solver is described. A structure solver with the ability to handle geometric nonlinearly is developed and tested with benchmark cases. The partitioned fluid-structure coupling algorithm with the strategy of enforcing boundary conditions at the fluid/structure interaction is given in detail. The fully coupled FSI approach is tested with the Turek and Hron fluid-structure interaction benchmark case. Both strong coupling and weak coupling algorithms are examined. Predictions from the current approach show good agreement with the results reported by other researchers.


Sign in / Sign up

Export Citation Format

Share Document