Numerical simulations of two-dimensional foam by the immersed boundary method

  • Authors:
  • Yongsam Kim;Ming-Chih Lai;Charles S. Peskin

  • Affiliations:
  • Department of Mathematics, Chung-Ang University, Dongjakgu Heukseokdong, Seoul 156-756, Republic of Korea;Center of Mathematical Modeling and Scientific Computing and Department of Applied Mathematics, National Chiao Tung University, 1001, Ta Hsueh Road, Hsinchu 300, Taiwan;Courant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New York, NY 10012, USA

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2010

Quantified Score

Hi-index 31.46

Visualization

Abstract

In this paper, we present an immersed boundary (IB) method to simulate a dry foam, i.e., a foam in which most of the volume is attributed to its gas phase. Dry foam dynamics involves the interaction between a gas and a collection of thin liquid-film internal boundaries that partition the gas into discrete cells or bubbles. The liquid-film boundaries are flexible, contract under the influence of surface tension, and are permeable to the gas, which moves across them by diffusion at a rate proportional to the local pressure difference across the boundary. Such problems are conventionally studied by assuming that the pressure is uniform within each bubble. Here, we introduce instead an IB method that takes into account the non-equilibrium fluid mechanics of the gas. To model gas diffusion across the internal liquid-film boundaries, we allow normal slip between the boundary and the gas at a velocity proportional to the (normal) force generated by the boundary surface tension. We implement this method in the two-dimensional case, and test it by verifying the von Neumann relation, which governs the coarsening of a two-dimensional dry foam. The method is further validated by a convergence study, which confirms its first-order accuracy.