A Finite Element Approach to the Immersed Boundary Method

2009 ◽  
pp. 271-298 ◽  
Author(s):  
D. Boffi ◽  
L. Gastaldi ◽  
L. Heltai
Keyword(s):  
Finite Element ◽  
Immersed Boundary Method ◽  
Immersed Boundary ◽  
Finite Element Approach ◽  
Boundary Method ◽  
Element Approach

FINITE ELEMENT APPROACH TO IMMERSED BOUNDARY METHOD WITH DIFFERENT FLUID AND SOLID DENSITIES

2011 ◽  
Vol 21 (12) ◽  
pp. 2523-2550 ◽  
Author(s):  
DANIELE BOFFI ◽  
NICOLA CAVALLINI ◽  
LUCIA GASTALDI
Keyword(s):  
Finite Element ◽  
Immersed Boundary Method ◽  
Stokes Equations ◽  
Numerical Procedure ◽  
Biological Tissues ◽  
Immersed Boundary ◽  
Fluid Structure ◽  
Finite Element Approach ◽  
Boundary Method ◽  
Element Approach

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.

A Divergence-Free Immersed Boundary Method and Its Finite Element Applications

2020 ◽  
Vol 36 (6) ◽  
pp. 901-914
Author(s):  
Chuan Zhou ◽  
Jianhua Li ◽  
Huaan Wang ◽  
Kailong Mu ◽  
Lanhao Zhao
Keyword(s):  
Finite Element ◽  
Immersed Boundary Method ◽  
Body Force ◽  
Iterative Scheme ◽  
Current Method ◽  
Slip Condition ◽  
Immersed Boundary ◽  
Boundary Method ◽  
Divergence Free ◽  
Characteristic Based Split

ABSTRACTIn order to maintain the no-slip condition and the divergence-free property simultaneously, an iterative scheme of immersed boundary method in the finite element framework is presented. In this method, the Characteristic-based Split scheme is employed to solve the momentum equations and the formulation for the pressure and the extra body force is derived according to the no-slip condition. The extra body force is divided into two divisions, one is in relation to the pressure and the other is irrelevant. Two corresponding independent iterations are set to solve the two sections. The novelty of this method lies in that the correction of the velocity increment is included in the calculation of the extra body force which is relevant to the pressure and the update of the force is incorporated into the iteration of the pressure. Hence, the divergence-free properties and no-slip conditions are ensured concurrently. In addition, the current method is validated with well-known benchmarks.

scholarly journals A Revisit of Implicit Monolithic Algorithms for Compressible Solids Immersed Inside a Compressible Liquid

Fluids ◽  
2021 ◽  
Vol 6 (8) ◽  
pp. 273
Author(s):  
Sheldon Wang
Keyword(s):  
Finite Element ◽  
Finite Element Methods ◽  
Immersed Boundary Method ◽  
Compressible Liquid ◽  
Working Fluid ◽  
Immersed Boundary ◽  
Constant Density ◽  
Liquid Density ◽  
Boundary Method ◽  
Matrix Free

With the development of mature Computational Fluid Dynamics (CFD) tools for fluids (air and liquid) and Finite Element Methods (FEM) for solids and structures, many approaches have been proposed to tackle the so-called Fluid–Structure Interaction or Fluid–Solid Interaction (FSI) problems. Traditional partitioned iterations are often used to link available FEM codes with CFD codes in the study of FSI systems. Although these procedures are convenient, fluid mesh adjustments according to the motion and finite deformation of immersed solids or structures can be challenging or even prohibitive. Moreover, complex dynamic behaviors of coupled FSI systems are often lost in these iterative processes. In this paper, the author would like to review the so-called monolithic approaches for the solution of coupled FSI systems as a whole in the context of the immersed boundary method. In particular, the focus is on the implicit monolithic algorithm for compressible solids immersed inside a compressible liquid. Notice here the main focus of this paper is on liquid or more precisely liquid phase of water as working fluid. Using the word liquid, the author would like to emphasize the consideration of the compressibility of the fluid and the assumption of constant density and temperature. It is a common practice to assume that the pressure variations are not strong enough to alter the liquid density in any significant fashion for acoustic fluid–solid interactions problems. Although the algorithm presented in this paper is not directly applicable to aerodynamics in which the density change is significant along with its relationship with the pressure and the temperature, the author did revisit his earlier work on merging immersed boundary method concepts with a fully-fledged compressible aerodynamic code based on high-order compact scheme and energy conservative form of governing equations. In the proposed algorithm, on top of a uniform background (ghost) mesh, a fully implicit immersed method is implemented with mixed finite element methods for compressible liquid as well as immersed compressible solids with a matrix-free Newton–Krylov iterative solution scheme. In this monolithic approach, with the simple modulo function, the immersed solid or structure points can be easily located and thus the displacement projections and force distributions stipulated in the immersed boundary method can be effectively and efficiently implemented. This feature coupled with the key concept of the immersed boundary method helps to avoid topologically challenging mesh adjustments and to incorporate parallel processing commands such as Message Passing Interface (MPI) and further vectorization of the numerical operation. Once these high-performance procedures are implemented coupled with the monolithic implicit matrix-free Newton–Krylov iterative scheme with immersed methods, effective and efficient reduced order modeling techniques can then be employed to explore phase and parametric spaces. The in-house developed programs are at the moment two-dimensional. Furthermore, based on the same approach implemented in one-dimensional test example with one continuum immersed in another continuum, such monolithic implicit matrix-free Newton–Krylov iterative approach can be extended for the study of composites with deformable aggregates and matrix.

Sign in / Sign up

Export Citation Format

Share Document