Anti-diffusion method for interface steepening in two-phase incompressible flow

In this paper, we present a method for obtaining sharp interfaces in two-phase incompressible flows by an anti-diffusion correction, that is applicable in a straight-forward fashion for the improvement of two-phase flow solution schemes typically employed in practical applications. The underlying discretization is based on the volume-of-fluid (VOF) interface-capturing method on unstructured meshes. The key idea is to steepen the interface, independently of the underlying volume-fraction transport equation, by solving a diffusion equation with reverse time, i.e. an anti-diffusion equation, after each advection time step of the volume fraction. As the solution of the anti-diffusion equation requires regularization, a limiter based on the directional derivative is developed for calculating the gradient of the volume fraction. This limiter ensures the boundedness of the volume fraction. In order to control the amount of anti-diffusion introduced by the correction algorithm we propose a suitable stopping criterion for interface steepening. The formulation of the limiter and the algorithm for solving the anti-diffusion equation are applicable to 3-dimensional unstructured meshes. Validation computations are performed for passive advection of an interface, for 2-dimensional and 3-dimensional rising-bubbles, and for a rising drop in a periodically constricted channel. The results demonstrate that sharp interfaces can be recovered reliably. They show that the accuracy is similar to or even better than that of level-set methods using comparable discretizations for the flow and the level-set evolution. Also, we observe a good agreement with experimental results for the rising drop where proper interface evolution requires accurate mass conservation.

[1]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[2]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[3]  W. Rider,et al.  Reconstructing Volume Tracking , 1998 .

[4]  Martin E. Weber,et al.  Bubbles in viscous liquids: shapes, wakes and velocities , 1981, Journal of Fluid Mechanics.

[5]  R. I. Issa,et al.  A Method for Capturing Sharp Fluid Interfaces on Arbitrary Meshes , 1999 .

[6]  Jack R. Edwards,et al.  An investigation of interface-sharpening schemes for multi-phase mixture flows , 2009, J. Comput. Phys..

[7]  Mark Sussman,et al.  A sharp interface method for incompressible two-phase flows , 2007, J. Comput. Phys..

[8]  S. Zaleski,et al.  DIRECT NUMERICAL SIMULATION OF FREE-SURFACE AND INTERFACIAL FLOW , 1999 .

[9]  S. Osher,et al.  Level set methods: an overview and some recent results , 2001 .

[10]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[11]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[12]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[13]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[14]  G. Kreiss,et al.  A conservative level set method for two phase flow II , 2005, Journal of Computational Physics.

[15]  David J. Torres,et al.  The point-set method: front-tracking without connectivity , 2000 .

[16]  Jean-François Remacle,et al.  A stabilized finite element method using a discontinuous level set approach for the computation of bubble dynamics , 2007, J. Comput. Phys..

[17]  F. Xiao,et al.  Numerical simulations of free-interface fluids by a multi-integrated moment method , 2005 .

[18]  Yehoshua Y. Zeevi,et al.  Forward-and-backward diffusion processes for adaptive image enhancement and denoising , 2002, IEEE Trans. Image Process..

[19]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

[20]  S. Zaleski,et al.  Modelling Merging and Fragmentation in Multiphase Flows with SURFER , 1994 .

[21]  A. Borhan,et al.  BUOYANCY-DRIVEN MOTION OF DROPS AND BUBBLES IN A PERIODICALLY CONSTRICTED CAPILLARY , 1996 .

[22]  Metin Muradoglu,et al.  An auxiliary grid method for computations of multiphase flows in complex geometries , 2006, J. Comput. Phys..

[23]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[24]  A. Gosman,et al.  High resolution NVD differencing scheme for arbitrarily unstructured meshes , 1999 .

[25]  Ying Sun,et al.  Sharp interface tracking using the phase-field equation , 2007, J. Comput. Phys..

[26]  Alan S. Willsky,et al.  Image segmentation and edge enhancement with stabilized inverse diffusion equations , 2000, IEEE Trans. Image Process..

[27]  Jean-François Remacle,et al.  Transient adaptivity applied to two-phase incompressible flows , 2008, J. Comput. Phys..

[28]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[29]  William J. Rider,et al.  Simple modifications of monotonicity-preserving limiters , 2001 .

[30]  H. Rusche.  Computational fluid dynamics of dispersed two-phase flows at high phase fractions , 2003 .

[31]  Jay P. Boris,et al.  Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works , 1973 .

[32]  J. G. Hnat,et al.  Spherical cap bubbles and skirt formation , 1976 .

[33]  M. Welk,et al.  Staircasing in semidiscrete stabilised inverse linear diffusion algorithms , 2007 .

[34]  M. Sussman,et al.  A Coupled Level Set and Volume-of-Fluid Method for Computing 3D and Axisymmetric Incompressible Two-Phase Flows , 2000 .

[35]  Chi-Wang Shu,et al.  Anti-diffusive flux corrections for high order finite difference WENO schemes , 2005 .

[36]  Feng Xiao,et al.  An efficient method for capturing free boundaries in multi‐fluid simulations , 2003 .

[37]  Feng Xiao,et al.  A simple algebraic interface capturing scheme using hyperbolic tangent function , 2005 .

[38]  Thomas Brox,et al.  Stabilised Nonlinear Inverse Diffusion for Approximating Hyperbolic PDEs , 2005, Scale-Space.

[39]  L. Rudin,et al.  Feature-oriented image enhancement using shock filters , 1990 .

[40]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[41]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[42]  Jian-Guo Liu,et al.  Generalized monotone schemes, discrete paths of extrema, and discrete entropy conditions , 1999, Math. Comput..