Library Subscription: Guest
Begell Digital Portal Begell Digital Library eBooks Journals References & Proceedings Research Collections
Multiphase Science and Technology
SJR: 0.124 SNIP: 0.222 CiteScore™: 0.26

ISSN Print: 0276-1459
ISSN Online: 1943-6181

Multiphase Science and Technology

DOI: 10.1615/MultScienTechn.v21.i1-2.130
pages 169-183


Srdjan Sasic
Department of Applied Mechanics, Chalmers, 412 96 Göteborg, Sweden
Berend G. M. van Wachem
Formerly at the Department of Mechanical Engineering, Imperial College London, SW7 2AZ, UK; Lehrstuhl für Technische Thermodynamik, Otto-von-Guericke-Universität Magdeburg, 39106 Germany


An immersed boundary method for three-dimensional, time-dependent flows is presented in this work and applied to simulating the behavior of an individual fiber in various flow regimes. The fiber is placed in a periodic box and has either a fixed position or is allowed to move freely (including translation and rotation) through the domain. The immersed boundary method is used to match the fluid velocity with the velocity of the interface of the fiber by mirroring the velocity field along the normal of the local triangulated immersed boundary segment to guarantee that the fluid accurately takes into account the presence of an immersed body. As a result of the procedure, there is a fictitious velocity field inside the immersed boundary, mirroring the boundary layer. The method applied is second-order accurate for the drag on the fiber and is intended to be used for fully resolving the flow field around arbitrary moving bodies immersed in a fluid. The immersed boundary method is employed on a selection of different fiber shapes, aiming at predicting the behavior of real fibers in realistic flow situations. A grid refinement study and study of the influence of the size of the periodic box used for the simulations are carried out. It is shown by grid refinement that the simulations performed here are truly direct numerical simulation (DNS). The force exerted by the fluid on a fiber is directly calculated by integrating the pressure and viscous forces over the objects immersed. The resulting coarse-grained drag and lift force functions can be employed by calculations of fluid-fiber flows on a larger scale (e.g., Eulerian-Eulerian simulations of air-fiber flows).


  1. Balay, S., Buschelman, K., Ejikhout, V., Gropp, W. D., Kaushnik, D., Kneple, M. G., and Smith, B. F., PETSc User Manual. DOI: 10.2172/464159

  2. Chhabra, R. P., Agarwal, L., and Sinha, N. H., Drag on non-spherical particles: An evaluation of available methods. DOI: 10.1016/S0032-5910(98)00178-8

  3. Clift, R. and Gauvin, W. H., Proceedings CHEMECA '70.

  4. Ganser, G. H., A rational approach to drag prediction of spherical and non-spherical particles.

  5. GTS, Triangulated Surface Library.

  6. Hu, H. H., Patankar, N. A., and Zhu, M. Y., Direct numerical simulation of fluid-solid systems using arbitrary Lagrangian-Eulerian technique. DOI: 10.1006/jcph.2000.6592

  7. Lai, M. C. and Peskin, C. S., An immersed boundary method with formal second-order accuracy and reduced numerical viscosity. DOI: 10.1006/jcph.2000.6483

  8. Loth, E., Drag of non-spherical solid particles of regular and irregular shape. DOI: 10.1016/j.powtec.2007.06.001

  9. Mark, A. and van Wachem, B., Derivation and validation of a novel implicit second order immersed boundary method. DOI: 10.1016/

  10. Mohd-Yusof, H., Development of Immersed Boundary Methods for Complex Geometries.

  11. Shepard, D., A two-dimensional interpolation function for irregularly spaced data. DOI: 10.1145/800186.810616

  12. van Doormaal, J. P. and Raithby, G. D., Enhancements of the simple method for predicting incompressible fluid flows. DOI: 10.1080/01495728408961817

  13. van Wachem, B., Multiflow - A Coupled Multiphase Flow Solver.