A front-tracking method with Catmull-Clark subdivision surfaces for studying liquid capsules enclosed by thin shells in shear flow

  • Authors:
  • Duc-Vinh Le;Sum Thai Wong

  • Affiliations:
  • Institute of High Performance Computing, A*STAR, 1 Fusionopolis Way, #16-16 Connexis, Singapore 138632, Singapore;Institute of High Performance Computing, A*STAR, 1 Fusionopolis Way, #16-16 Connexis, Singapore 138632, Singapore

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

Quantified Score

Hi-index 31.45

Visualization

Abstract

This paper presents a front-tracking method for studying the large deformation of a liquid capsule enclosed by a thin shell in a shear flow. The interaction between the fluid and the shell body is accomplished through an implicit immersed boundary method. An improved thin-shell model for computing the forces acting on the shell middle surface during the deformation is described in surface curvilinear coordinates and within the framework of the principle of virtual displacements. This thin-shell model takes full account of in-plane tensions and bending moments developing due to the shell thickness and a preferred three-dimensional membrane structure. The approximation of the shell middle surface is performed through the use of the Catmull-Clark subdivision surfaces. The resulting limit surface is C^2-continuous everywhere except at a small number of extraordinary nodes where it retains C^1 continuity. The smoothness of the limit surface significantly improves the ability of our method in simulating capsules enclosed by hyperelastic thin shells with different shapes and physical properties. The present numerical technique has been validated by several examples including an inflation of a spherical shell and deformations of spherical, ellipsoidal and biconcave capsules in the shear flow. In addition, different types of motion such as tank-treading, swinging, tumbling and transition from tumbling to swinging have been studied over a range of shear rates, viscosity ratios and bending modulus.