The classical model for studying one-phase Hele-Shaw flows is based on a highly nonlinear moving boundary problem with the fluid velocity related to pressure gradients via a Darcy-type law. In a standard configuration with the Hele-Shaw cell made up of two flat stationary plates, the pressure is harmonic. Therefore, conformal mapping techniques and boundary integral methods can be readily applied to study the key interfacial dynamics, including the Saffman–Taylor instability and viscous fingering patterns. As well as providing a brief review of these key issues, we present a flexible numerical scheme for studying both the standard and nonstandard Hele-Shaw flows. Our method consists of using a modified finite-difference stencil in conjunction with the level-set method to solve the governing equation for pressure on complicated domains and track the location of the moving boundary. Simulations show that our method is capable of reproducing the distinctive morphological features of the Saffman–Taylor instability on a uniform computational grid. By making straightforward adjustments, we show how our scheme can easily be adapted to solve for a wide variety of nonstandard configurations, including cases where the gap between the plates is linearly tapered, the plates are separated in time, and the entire Hele-Shaw cell is rotated at a given angular velocity.

A new Semi-Lagrangian scheme is proposed to discretize the surface convection-diffusion equation. The other involved equations including the the level-set convection equation, the re-initialization equation and the extension equation are also solved by S-L schemes. The S-L method removes both the CFL condition and the stiffness caused by the surface Laplacian, allowing larger time step than the Eulerian method. The method is extended to the block-structured adaptive mesh. Numerical examples are given to demonstrate the efficiency of the S-L method.

Semi-Lagrangian (S-L) methods have no CFL stability constraint, and are more stable than the Eulerian methods. In the literature, the S-L method for the level-set re-initialization equation was complicated, which may be unnecessary. Since the re-initialization procedure is auxiliary, we propose to use the first-order S-L scheme coupled with a projection technique to improve the accuracy at the grid points just adjacent to the interface. Standard second-order S-L method is used for evolving the level-set convection equation. The implementation is simple, including on the block-structured adaptive mesh. The efficiency of the S-L method is demonstrated by extensive numerical examples including passive convection of interfaces with corners/kinks/large deformation under given velocity fields, a geometrical flow with topological changes, simulations of bubble/ droplet dynamics in incompressible two-phase flows. In terms of accuracy it is comparable to the other existing methods.

We propose and analyze a constrained level-set method for semi-automatic image segmentation. Our level-set model with constraints on the level-set function enables us to specify which parts of the image lie inside respectively outside the segmented objects. Such a-priori information can be expressed in terms of upper and lower constraints prescribed for the level-set function. Constraints have the same conceptual meaning as initial seeds of the popular graph-cuts based methods for image segmentation. A numerical approximation scheme is based on the complementary-finite volumes method combined with the Projected successive over-relaxation method adopted for solving constrained linear complementarity problems. The advantage of the constrained level-set method is demonstrated on several artificial images as well as on cardiac MRI data.

In this paper, a numerical method is presented for simulating the 3D interfacial flows with insoluble surfactant. The numerical scheme consists of a 3D immersed interface method (IIM) for solving Stokes equations with jumps across the interface and a 3D level-set method for solving the surfactant convection-diffusion equation along a moving and deforming interface. The 3D IIM Poisson solver modifies the one in the literature by assuming that the jump conditions of the solution and the flux are implicitly given at the grid points in a small neighborhood of the interface. This assumption is convenient in conjunction with the level-set techniques. It allows standard Lagrangian interpolation for quantities at the projection points on the interface. The interface jump relations are re-derived accordingly. A novel rotational procedure is given to generate smooth local coordinate systems and make effective interpolation. Numerical examples demonstrate that the IIM Poisson solver and the Stokes solver achieve second-order accuracy. A 3D drop with insoluble surfactant under shear flow is investigated numerically by studying the influences of different physical parameters on the drop deformation.

In this paper, we numerically investigate the effects of surfactant on drop-drop interactions in a 2D shear flow using a coupled level-set and immersed interface approach proposed in (Xu et al., J. Comput. Phys., 212 (2006), 590-616). We find that surfactant plays a critical and nontrivial role in drop-drop interactions. In particular, we find that the minimum distance between the drops is a non-monotone function of the surfactant coverage and Capillary number. This non-monotonic behavior, which does not occur for clean drops, is found to be due to the presence of Marangoni forces along the drop interfaces. This suggests that there are non-monotonic conditions for coalescence of surfactant-laden drops, as observed in recent experiments of Leal and co-workers. Although our study is two-dimensional, we believe that drop-drop interactions in three-dimensional flows should be qualitatively similar as the Maragoni forces in the near contact region in 3D should have a similar effect.

The level set method has become widely used in shape optimization where it allows a popular implementation of the steepest descent method. Once coupled with a ersatz material approximation [Allaire et al., J. Comput. Phys.194 (2004) 363–393], a single mesh is only used leading to very efficient and cheap numerical schemes in optimization of structures. However, it has some limitations and cannot be applied in every situation. This work aims at exploring such a limitation. We estimate the systematic error committed by using the ersatz material approximation and, on a model case, explain that they amplifies instabilities by a second order analysis of the objective function.