American Institute of Mathematical Sciences

March  2021, 14(3): 819-833. doi: 10.3934/dcdss.2020349

Computational study of immersed boundary - lattice Boltzmann method for fluid-structure interaction

 a. Department of Mathematics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, Trojanova 13,120 00 Prague, Czech Republic b. Department of Heat Engineering and Environment Protection, AGH University of Science and Technology, Al. Mickiewicza 30, 30-059 Krakow, Poland

* Corresponding author: eichlpa1@fjfi.cvut.cz

Received  January 2019 Revised  December 2019 Published  May 2020

In this article, we deal with the numerical immersed boundary-lattice Boltzmann method for simulation of the fluid-structure interaction problems in 2D. We consider the interaction of incompressible, Newtonian fluid in an isothermal system with an elastic fiber, which represents an immersed body boundary. First, a short introduction to the lattice Boltzmann and immersed boundary method is presented and the combination of these two methods is briefly discussed. Then, the choice of the smooth approximation of the Dirac delta function and the discretization of the immersed body is discussed. One of the significant drawbacks of immersed boundary method is the penetrative flow through the immersed impermeable boundary. The effect of the immersed body boundary discretization is investigated using two benchmark problems, where an elastic fiber is deformed. The results indicate that the restrictions placed on the discretization in literature are not necessary.

Citation: Pavel Eichler, Radek Fučík, Robert Straka. Computational study of immersed boundary - lattice Boltzmann method for fluid-structure interaction. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 819-833. doi: 10.3934/dcdss.2020349
References:

show all references

References:
The setup of the computational domain in the case of the Kármán vortex flow problem. Two obstacles are immersed in the fluid. The first one is a rigid cylinder with diameter $d$. The last one is an elastic fiber with one fixed end (filled circle) and one free end (empty circle)
The study of the elastic fiber deformation in the Kármán vortex street region. In Figs. 2a and 2b, the free end trajectory is shown for two methods: (A) Forward Euler method (red), (B) The third order Runge-Kutta method (blue), and the results are compared with those in [22] (black circles). In Fig. 2c, the superposition of the elastic fibers is shown
The study of the elastic fiber deformation in the Kármán vortex street region with neglected inertial forces of the elastic fiber. In Figs. 3a and 3b the free end trajectory is shown for two methods: (A) Forward Euler method (red), (B) The third order Runge-Kutta method (blue), and the results are compared with results in [22] (black triangles). In Fig. 3c, the superposition of the elastic fibers is shown
The setup of the computational domain for the cavity flow problem
Evolution of $E_{rkl}$ w.r.t. ${\sigma}/{h}$ for different values of the stiffness coefficient $k_S$ and for different mesh spacings $h$: $h_1 = 1.6 \cdot 10^{-4}\;\mathrm{m}$, $h_2 = 7.84\cdot 10^{-5}\;\mathrm{m}$, and $h_3 = 5.24\cdot 10^{-5}\;\mathrm{m}$. The results computed at the final time $T = 10\;\mathrm{s}$ when the steady state is reached
Schematic visualization of the fluid penetration through the flexible fiber for three values of ${\sigma}/{h}$. In the figures, the velocity vectors (white arrows) and values of the fluid velocity magnitude (color scale) are shown. The results are at the final time $T = 10\;\mathrm{s}$ when the steady state is reached. All results are computed using $k_S = 100\;\mathrm{kg\, s^{-2}}$ and $h = 5.24\cdot 10^{-5}\; \mathrm{m}$
Schematic visualization of the fiber deformation at the final time $T = 10\, \mathrm{s}$ and for $h = 5.24 \cdot 10^{-5}$ m. Different colors of the fibers represent different values of ${\sigma}/{h}$
The setup of the computational domain in the case of elastic bump deformation
Investigation of fluid flow around the bump. In Figs. (A) and (B), the relation between the surface area $\mu\left(\Omega_u\right)$ and ${\sigma}/{h}$ is shown. Figs. (C) and (D) illustrate the evolution of $E_{rkl}$ with respect to ${\sigma}/{h}$. For $U_\infty = 0.5\, \mathrm{m\, s^{-1}}$, the values are time-averaged over the time interval of $(1, 9.5)$ s. For $U_\infty = 0.5\, \mathrm{m\, s^{-1}}$, the values are computed at the final time $T = 10\; \mathrm{s}$
Schematic visualization of the fluid flow around an impermeable flexible fiber in time $t = 9.5\, \mathrm{s}$. In the figures, the velocity vectors (white arrows) and the values of the fluid velocity magnitude (color scale) are shown
 [1] Armin Lechleiter, Tobias Rienmüller. Factorization method for the inverse Stokes problem. Inverse Problems & Imaging, 2013, 7 (4) : 1271-1293. doi: 10.3934/ipi.2013.7.1271 [2] Qiang Guo, Dong Liang. An adaptive wavelet method and its analysis for parabolic equations. Numerical Algebra, Control & Optimization, 2013, 3 (2) : 327-345. doi: 10.3934/naco.2013.3.327 [3] Tao Wu, Yu Lei, Jiao Shi, Maoguo Gong. An evolutionary multiobjective method for low-rank and sparse matrix decomposition. Big Data & Information Analytics, 2017, 2 (1) : 23-37. doi: 10.3934/bdia.2017006 [4] Deren Han, Zehui Jia, Yongzhong Song, David Z. W. Wang. An efficient projection method for nonlinear inverse problems with sparsity constraints. Inverse Problems & Imaging, 2016, 10 (3) : 689-709. doi: 10.3934/ipi.2016017 [5] Boris Kramer, John R. Singler. A POD projection method for large-scale algebraic Riccati equations. Numerical Algebra, Control & Optimization, 2016, 6 (4) : 413-435. doi: 10.3934/naco.2016018 [6] Petra Csomós, Hermann Mena. Fourier-splitting method for solving hyperbolic LQR problems. Numerical Algebra, Control & Optimization, 2018, 8 (1) : 17-46. doi: 10.3934/naco.2018002 [7] Christina Surulescu, Nicolae Surulescu. Modeling and simulation of some cell dispersion problems by a nonparametric method. Mathematical Biosciences & Engineering, 2011, 8 (2) : 263-277. doi: 10.3934/mbe.2011.8.263 [8] Min Li. A three term Polak-Ribière-Polyak conjugate gradient method close to the memoryless BFGS quasi-Newton method. Journal of Industrial & Management Optimization, 2020, 16 (1) : 245-260. doi: 10.3934/jimo.2018149 [9] Manfred Einsiedler, Elon Lindenstrauss. On measures invariant under diagonalizable actions: the Rank-One case and the general Low-Entropy method. Journal of Modern Dynamics, 2008, 2 (1) : 83-128. doi: 10.3934/jmd.2008.2.83 [10] Xiaomao Deng, Xiao-Chuan Cai, Jun Zou. A parallel space-time domain decomposition method for unsteady source inversion problems. Inverse Problems & Imaging, 2015, 9 (4) : 1069-1091. doi: 10.3934/ipi.2015.9.1069 [11] Marion Darbas, Jérémy Heleine, Stephanie Lohrengel. Numerical resolution by the quasi-reversibility method of a data completion problem for Maxwell's equations. Inverse Problems & Imaging, 2020, 14 (6) : 1107-1133. doi: 10.3934/ipi.2020056 [12] Madalina Petcu, Roger Temam. The one dimensional shallow water equations with Dirichlet boundary conditions on the velocity. Discrete & Continuous Dynamical Systems - S, 2011, 4 (1) : 209-222. doi: 10.3934/dcdss.2011.4.209 [13] M. Mahalingam, Parag Ravindran, U. Saravanan, K. R. Rajagopal. Two boundary value problems involving an inhomogeneous viscoelastic solid. Discrete & Continuous Dynamical Systems - S, 2017, 10 (6) : 1351-1373. doi: 10.3934/dcdss.2017072 [14] Feng Luo. A combinatorial curvature flow for compact 3-manifolds with boundary. Electronic Research Announcements, 2005, 11: 12-20. [15] Elvise Berchio, Filippo Gazzola, Dario Pierotti. Nodal solutions to critical growth elliptic problems under Steklov boundary conditions. Communications on Pure & Applied Analysis, 2009, 8 (2) : 533-557. doi: 10.3934/cpaa.2009.8.533 [16] Yizhuo Wang, Shangjiang Guo. A SIS reaction-diffusion model with a free boundary condition and nonhomogeneous coefficients. Discrete & Continuous Dynamical Systems - B, 2019, 24 (4) : 1627-1652. doi: 10.3934/dcdsb.2018223 [17] Mats Gyllenberg, Jifa Jiang, Lei Niu, Ping Yan. On the classification of generalized competitive Atkinson-Allen models via the dynamics on the boundary of the carrying simplex. Discrete & Continuous Dynamical Systems - A, 2018, 38 (2) : 615-650. doi: 10.3934/dcds.2018027 [18] Carlos Fresneda-Portillo, Sergey E. Mikhailov. Analysis of Boundary-Domain Integral Equations to the mixed BVP for a compressible stokes system with variable viscosity. Communications on Pure & Applied Analysis, 2019, 18 (6) : 3059-3088. doi: 10.3934/cpaa.2019137 [19] Teddy Pichard. A moment closure based on a projection on the boundary of the realizability domain: 1D case. Kinetic & Related Models, 2020, 13 (6) : 1243-1280. doi: 10.3934/krm.2020045 [20] Zhiming Guo, Zhi-Chun Yang, Xingfu Zou. Existence and uniqueness of positive solution to a non-local differential equation with homogeneous Dirichlet boundary condition---A non-monotone case. Communications on Pure & Applied Analysis, 2012, 11 (5) : 1825-1838. doi: 10.3934/cpaa.2012.11.1825

2019 Impact Factor: 1.233