News & Updates
Search Research Content
Resource Finder at Kennedy Krieger Institute
A free resource that provides access to information and support for individuals and families living with developmental disabilities.
New Finite Difference Methods Based on IIM for Inextensible Interfaces in Incompressible Flows.
|Title||New Finite Difference Methods Based on IIM for Inextensible Interfaces in Incompressible Flows.|
|Publication Type||Journal Article|
|Year of Publication||2011|
|Authors||Li Z, Lai M-C|
|Journal||East Asian journal on applied mathematics|
|Date Published||2011 Jan 1|
In this paper, new finite difference methods based on the augmented immersed interface method (IIM) are proposed for simulating an inextensible moving interface in an incompressible two-dimensional flow. The mathematical models arise from studying the deformation of red blood cells in mathematical biology. The governing equations are incompressible Stokes or Navier-Stokes equations with an unknown surface tension, which should be determined in such a way that the surface divergence of the velocity is zero along the interface. Thus, the area enclosed by the interface and the total length of the interface should be conserved during the evolution process. Because of the nonlinear and coupling nature of the problem, direct discretization by applying the immersed boundary or immersed interface method yields complex nonlinear systems to be solved. In our new methods, we treat the unknown surface tension as an augmented variable so that the augmented IIM can be applied. Since finding the unknown surface tension is essentially an inverse problem that is sensitive to perturbations, our regularization strategy is to introduce a controlled tangential force along the interface, which leads to a least squares problem. For Stokes equations, the forward solver at one time level involves solving three Poisson equations with an interface. For Navier-Stokes equations, we propose a modified projection method that can enforce the pressure jump condition corresponding directly to the unknown surface tension. Several numerical experiments show good agreement with other results in the literature and reveal some interesting phenomena.
|Alternate Journal||East Asian J Applied Math|