A Discrete-Forcing Immersed Boundary Method for the Fluid-Structure Interaction of an Elastic Slender Body

Author(s):  
Injae Lee ◽  
Haecheon Choi

In the present study, a new immersed boundary method for the simulation of flow around an elastic slender body is suggested. The present method is based on the discrete-forcing immersed boundary method by Kim et al. (J. Comput. Phys., 2001) and is fully coupled with the elastic slender body motion. The incompressible Navier-Stokes equations are solved in an Eulerian coordinate and the elastic slender body motion is described in a Lagrangian coordinate, respectively. The elastic slender body is modeled as a thin flexible beam and is segmented by finite number of blocks. Each block is then moved by the external and internal forces such as the hydrodynamic, tension, bending, and buoyancy forces. With the proposed method, we simulate several flow problems including flows over a flexible filament, an oscillating insect wing, and a flapping flag. We show that the present method does not impose any severe limitation on the size of computational time step. The results obtained agree very well with those from previous studies.

2017 ◽  
Vol 39 (2) ◽  
pp. 109-119
Author(s):  
Cuong Q. Le ◽  
H. Phan-Duc ◽  
Son H. Nguyen

In this paper, a combination of the Proper Generalized  Decomposition (PGD) with the Immersed Boundary method (IBM) for solving  fluid-filament interaction problem is proposed. In this combination, a  forcing term constructed by the IBM is introduced to Navier-Stokes equations  to handle the influence of the filament on the fluid flow. The PGD is  applied to solve the Poission's equation to find the fluid pressure  distribution for each time step. The numerical results are compared with  those by previous publications to illustrate the robustness and  effectiveness of the proposed method.


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.


Author(s):  
Karim M. Ali ◽  
Mohamed Madbouli ◽  
Hany M. Hamouda ◽  
Amr Guaily

This work introduces an immersed boundary method for two-dimensional simulation of incompressible Navier-Stokes equations. The method uses flow field mapping on the immersed boundary and performs a contour integration to calculate immersed boundary forces. This takes into account the relative location of the immersed boundary inside the background grid elements by using inverse distance weights, and also considers the curvature of the immersed boundary edges. The governing equations of the fluid mechanics are solved using a Galerkin-Least squares finite element formulation. The model is validated against a stationary and a vertically oscillating circular cylinder in a cross flow. The results of the model show acceptable accuracy when compared to experimental and numerical results.


2013 ◽  
Vol 477-478 ◽  
pp. 281-284
Author(s):  
Jie Yang ◽  
Song Ping Wu

An immersed boundary method based on the ghost-cell approach is presented in this paper. The compressible Navier-Stokes equations are discretized using a flux-splitting method for inviscid fluxes and second-order central-difference for the viscous components. High-order accuracy is achieved by using weighted essentially non-oscillatory (WENO) and Runge-Kutta schemes. Boundary conditions are reconstructed by a serial of linear interpolation and inverse distance weighting interpolation of flow variables in fluid domain. Two classic flow problems (flow over a circular cylinder, and a NACA 0012 airfoil) are simulated using the present immersed boundary method, and the predictions show good agreement with previous computational results.


Computation ◽  
2018 ◽  
Vol 6 (3) ◽  
pp. 50
Author(s):  
Jonatas Borges ◽  
Marcos Lourenço ◽  
Elie Padilla ◽  
Christopher Micallef

The immersed boundary method has attracted considerable interest in the last few years. The method is a computational cheap alternative to represent the boundaries of a geometrically complex body, while using a cartesian mesh, by adding a force term in the momentum equation. The advantage of this is that bodies of any arbitrary shape can be added without grid restructuring, a procedure which is often time-consuming. Furthermore, multiple bodies may be simulated, and relative motion of those bodies may be accomplished at reasonable computational cost. The numerical platform in development has a parallel distributed-memory implementation to solve the Navier-Stokes equations. The Finite Volume Method is used in the spatial discretization where the diffusive terms are approximated by the central difference method. The temporal discretization is accomplished using the Adams-Bashforth method. Both temporal and spatial discretizations are second-order accurate. The Velocity-pressure coupling is done using the fractional-step method of two steps. The present work applies the immersed boundary method to simulate a Newtonian laminar flow through a three-dimensional sudden contraction. Results are compared to published literature. Flow patterns upstream and downstream of the contraction region are analysed at various Reynolds number in the range 44 ≤ R e D ≤ 993 for the large tube and 87 ≤ R e D ≤ 1956 for the small tube, considerating a contraction ratio of β = 1 . 97 . Comparison between numerical and experimental velocity profiles has shown good agreement.


Sign in / Sign up

Export Citation Format

Share Document