Special Section Articles

An Improved Direct-Forcing Immersed Boundary Method for Fluid-Structure Interaction Simulations1

[+] Author and Article Information
Xiaojue Zhu

Graduate Research Assistant
LNM, Institute of Mechanics,
Chinese Academy of Sciences,
Beijing 100190, China
e-mail: zhuxiaojue@lnm.imech.ac.cn

Guowei He

Professor and Director
LNM, Institute of Mechanics,
Chinese Academy of Sciences,
Beijing 100190, China
e-mail: hgw@lnm.imech.ac.cn

Xing Zhang

Associate Professor
LNM, Institute of Mechanics,
Chinese Academy of Sciences,
Beijing 100190, China
e-mail: zhangx@lnm.imech.ac.cn

This manuscript is submitted to the special issue on immersed boundary method. Preliminary version of this paper appeared as FEDSM2013-16472 in ASME 2013 Fluids Engineering Division Summer Meeting, Incline Village, NV, July 7–11, 2013.

2Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received May 13, 2013; final manuscript received December 4, 2013; published online February 28, 2014. Assoc. Editor: Zhongquan Charlie Zheng.

J. Fluids Eng 136(4), 040903 (Feb 28, 2014) (9 pages) Paper No: FE-13-1310; doi: 10.1115/1.4026197 History: Received May 13, 2013; Revised December 04, 2013

In the present work, we present an improved version of the direct-forcing immersed boundary (IB) method proposed in Wang and Zhang (2011, “An Immersed Boundary Method Based on Discrete Stream Function Formulation for Two- and Three-Dimensional Incompressible Flows,” J. Comput. Phys., 230(9), pp. 3479–3499). In order to obtain an accurate prediction of local surface force, measures have been taken to suppress the unphysical spatial oscillations in the Lagrangian forcing. A fluid-structure interaction (FSI) solver has been developed by using the improved IB method for the fluid and the finite difference method for the structure. Several flow problems are simulated to validate our method. The testing cases include flows over a stationary cylinder and a stationary flat plate, two-dimensional flow past an inextensible flexible filament, and three-dimensional flow past a flapping flag. The results obtained in the present work agree well with those from the literature.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Peskin, C. S., 2003, “The Immersed Boundary Method,” Acta Numer., 11, pp. 479–517. [CrossRef]
Mittal, R., and Iaccarino, G., 2005, “Immersed Boundary Methods,” Ann. Rev. Fluid Mech., 37, pp. 239–261. [CrossRef]
Balaras, E., and Yang, J., 2005, “Nonboundary Conforming Methods for Large-Eddy Simulations of Biological Flows,” ASME J. Fluids Eng., 127(5), pp. 851–857. [CrossRef]
Zhang, N., and Zheng, Z. C., 2007, “An Improved Direct-Forcing Immersed-Boundary Method for Finite Difference Applications,” J. Comput. Phys., 221(1), pp. 250–268. [CrossRef]
Narasimhan, M., Dong, H., Mittal, R., and Singh, S. N., 2005, “Optimal Yaw Regulation and Trajectory Control of Biorobotic AUV Using Mechanical Fins Based on CFD Parameterization,” ASME J. Fluids Eng., 128(4), pp. 687–698. [CrossRef]
Wang, S., and Zhang, X., 2011, “An Immersed Boundary Method Based on Discrete Stream Function Formulation for Two- and Three-Dimensional Incompressible Flows,” J. Comput. Phys., 230(9), pp. 3479–3499. [CrossRef]
Su, S. W., Lai, M. C., and Lin, C. A., 2007, “An Immersed Boundary Technique for Simulating Complex Flows With Rigid Boundary,” Comput. Fluids, 36(2), pp. 313–324. [CrossRef]
Uhlmann, M., 2005, “An Immersed Boundary Method With Direct Forcing for the Simulation of Particulate Flows,” J. Comput. Phys., 209(2), pp. 448–476. [CrossRef]
Zhu, L., and Peskin, C. S., 2002, “Simulation of a Flapping Flexible Filament in a Flowing Soap Film by the Immersed Boundary Method,” J. Comput. Phys., 179(2), pp. 452–468. [CrossRef]
Connell, B. S. H., and Yue, D. K. P., 2007, “Flapping Dynamics of a Flag in a Uniform Stream,” J. Fluid Mech., 581, pp. 33–67. [CrossRef]
Huang, W.-X., Shin, S. J., and Sung, H. J., 2007, “Simulation of Flexible Filaments in a Uniform Flow by the Immersed Boundary Method,” J. Comput. Phys., 226(2), pp. 2206–2228. [CrossRef]
Huang, W. X., and Sung, H. J., 2010, “Three-Dimensional Simulation of a Flapping Flag in a Uniform Flow,” J. Fluid Mech., 653, pp. 301–336. [CrossRef]
Tornberg, A.-K., and Engquist, B., 2004, “Numerical Approximations of Singular Source Terms in Differential Equations,” J. Comput. Phys., 200(2), pp. 462–488. [CrossRef]
Linnick, M. N., and Fasel, H. F., 2005, “A High-Order Immersed Interface Method for Simulating Unsteady Incompressible Flows on Irregular Domains,” J. Comput. Phys., 204(1), pp. 157–192. [CrossRef]
Taira, K., and Colonius, T., 2007, “The Immersed Boundary Method: A Projection Approach,” J. Comput. Phys., 225(2), pp. 2118–2137. [CrossRef]
Majumdar, S., Iaccarino, G., and Durbin, P., 2001, “RANS Solvers With Adaptive Structured Boundary Non-Conforming Grids,” Annual Research Briefs, Center for Turbulence Research, Stanford University, pp. 353–466.


Grahic Jump Location
Fig. 1

Schematic representation of the Lagrangian coordinate system s on the filament. The length of the filament is L.

Grahic Jump Location
Fig. 2

Schematic representation of the Lagrangian coordinate system s1 and s2 on the flag. The longitudinal coordinate s1 ranges from 0 to L, and the spanwise coordinate s2 ranges from 0 to H, where L and H denote the length and width of the flag, respectively.

Grahic Jump Location
Fig. 3

The distribution of Lagrangian force correction F'x along the cylinder's surface by using Eq. (23) and setting the right hand side of Eq. (23) to −100 at each Lagrangian point. The Lagrangian grid width and Eulerian grid width are both 0.02. θ denotes the angle between the radical direction and the horizontal direction.

Grahic Jump Location
Fig. 4

Two ways of computing forcing component on a staggered mesh: (a) original way used in Ref. [6] and (b) more consistent way of computing forcing component. In (a), the Eulerian forcing (vector) is defined at cell centers, and each forcing component is interpolated (individually) to cell edges via simple average. In (b), each Eulerian forcing component is defined at cell edges, and no extra interpolation is needed.

Grahic Jump Location
Fig. 5

The distribution of Lagrangian intermediate velocity U→* along the cylinder's surface. The Eulerian grid width is 0.02 while in case (a) Δs/h = 1 and in case (b) Δs/h = 2.

Grahic Jump Location
Fig. 6

Distribution of (a) pressure coefficient Cp and (b) skin-friction coefficient Cf on the surface of the cylinder for Re = 40

Grahic Jump Location
Fig. 7

Distributions of pressure coefficient difference ΔCp and skin-friction coefficient (Cf) along the flat-plate surface at Re = 200 and two angles of attack: (a) Cf for α = 0deg, (b) ΔCp for α = 10deg, and (c) Cf for α = 10deg

Grahic Jump Location
Fig. 8

Comparison of the predicted free-end position with the analytical result

Grahic Jump Location
Fig. 9

Vorticity contours in wake of a flapping flexible filament

Grahic Jump Location
Fig. 10

Time history of y-position of the trailing edge: (a) γ=0.0015 and (b) γ = 0.0

Grahic Jump Location
Fig. 11

The instantaneous shape of a flapping flag in the three-dimensional simulation

Grahic Jump Location
Fig. 12

Time history of the transverse displacement of point A in Fig. 2 for Fr = 0.0



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In