scholarly journals A finite element immersed boundary method for fluid flow around rigid objects

2011 ◽  
Vol 65 (7) ◽  
pp. 856-875 ◽  
Author(s):  
F. Ilinca ◽  
J.-F. Hétu
Keyword(s):  
Finite Element ◽  
Fluid Flow ◽  
Immersed Boundary Method ◽  
Immersed Boundary ◽  
Rigid Objects ◽  
Boundary Method

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.

COMPUTATIONAL SIMULATION OF FLUID FLOW OVER A TRIANGULAR STEP USING IMMERSED BOUNDARY METHOD

2013 ◽  
Vol 10 (04) ◽  
pp. 1350016 ◽  
Author(s):  
C. A. SALEEL ◽  
A. SHAIJA ◽  
S. JAYARAJ
Keyword(s):  
Fluid Flow ◽  
Immersed Boundary Method ◽  
Reynolds Numbers ◽  
Stream Line ◽  
Immersed Boundary ◽  
Complex Geometries ◽  
Backward Facing Step ◽  
Fractional Step Method ◽  
Boundary Method ◽  
Fluid Solid Interaction

Handling of complex geometries with fluid–solid interaction has been one of the exigent issues in computational fluid dynamics (CFD) because most engineering problems have complex geometries with fluid–solid interaction for the purpose. Two different approaches have been developed for the same hitherto: (i) The unstructured grid method and (ii) the immersed boundary method (IBM). This paper details the IBM for the numerical investigation of two-dimensional laminar flow over a backward facing step and various geometrically configured triangular steps in hydro-dynamically developing regions (entrance region) as well in the hydro-dynamically developed regions through a channel at different Reynolds numbers. The present numerical method is rooted in a finite volume approach on a staggered grid in concert with a fractional step method. Geometrical obstructions are treated as an immersed boundary (IB), both momentum forcing and mass source terms are applied on the obstruction to satisfy the no-slip boundary condition and also to satisfy the continuity for the mesh containing the immersed boundary. Initially, numerically obtained velocity profiles and stream line plots for fluid flow over backward facing step is depicted to show its excellent agreement with the published results in various literatures. There after profiles and plots in the channel with triangular steps are also being unveiled with in depth elucidation. Results are presented for different Reynolds numbers.

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.

Numerical Simulation of Flow in a Wavy Wall Microchannel Using Immersed Boundary Method

2020 ◽  
Vol 13 (2) ◽  
pp. 118-125
Author(s):  
Mithun Kanchan ◽  
Ranjith Maniyeri
Keyword(s):  
Reynolds Number ◽  
Fluid Flow ◽  
Immersed Boundary Method ◽  
Stokes Equations ◽  
Flow Characteristics ◽  
Immersed Boundary ◽  
Solid Boundary ◽  
Wavy Wall ◽  
Dirac Delta ◽  
Boundary Method

Background: Fluid flow in microchannels is restricted to low Reynolds number regimes and hence inducing chaotic mixing in such devices is a major challenge. Over the years, the Immersed Boundary Method (IBM) has proved its ability in handling complex fluid-structure interaction problems. Objectives: Inspired by recent patents in microchannel mixing devices, we study passive mixing effects by performing two-dimensional numerical simulations of wavy wall in channel flow using IBM. Methods: The continuity and Navier-Stokes equations governing the flow are solved by fractional step based finite volume method on a staggered Cartesian grid system. Fluid variables are described by Eulerian coordinates and solid boundary by Lagrangian coordinates. A four-point Dirac delta function is used to couple both the coordinate variables. A momentum forcing term is added to the governing equation in order to impose the no-slip boundary condition between the wavy wall and fluid interface. Results: Parametric study is carried out to analyze the fluid flow characteristics by varying amplitude and wavelength of wavy wall configurations for different Reynolds number. Conclusion: Configurations of wavy wall microchannels having a higher amplitude and lower wavelengths show optimum results for mixing applications.

Sign in / Sign up

Export Citation Format

Share Document