A quadtree adaptive method for simulating fluid flows with moving interfaces

In this paper, a computational method for solving fluid flow problems with moving interfaces is presented. Herein, adaptive quadtree grids are used coupled with the CICSAM [O. Ubbink, Numerical prediction of two fluid systems with sharp interfaces, PhD Thesis, Imperial College of Science, Technology and Medicine, London, 1997] free surface capturing volume of fluid (VoF) method and PLIC reconstruction to interpolate the volume fraction field during refinement and derefinement processes. The combination of high resolution adaptive hierarchical remeshing and CICSAM interface advection is shown to overcome the problems of interface smearing and high CPU intensivity inherent in most VoF schemes. The result is a combination of free surface tracking and free surface capturing in that the interface is effectively tracked by the adapting refinements in the quadtree grid. In this way, a sharp interface is achieved and the advantages of both free surface tracking and capturing are combined. The new method is applied to interface advection examples in translating, rotating and shearing flow fields, and the benefits of using adapting quadtree grids demonstrated.

[1]  Scott D. Phillips,et al.  Computational and experimental analysis of dynamics of drop formation , 1999 .

[2]  Deborah Greaves,et al.  On the use of adaptive hierarchical meshes for numerical simulation of separated flows , 1998 .

[3]  M. Rudman.  INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 24, 671–691 (1997) VOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS , 2022 .

[4]  L. E. Scriven,et al.  Three-Dimensional Stability Analysis of Free Surface Flows , 1999 .

[5]  Tadashi Sakaguchi,et al.  Numerical Analysis of a Single Bubble by VOF Method. , 1993 .

[6]  Elke A. Rundensteiner,et al.  Fast, adaptive summation of point forces in the two-dimensional Poisson equation , 1989 .

[7]  S. Zaleski,et al.  Volume-of-Fluid Interface Tracking with Smoothed Surface Stress Methods for Three-Dimensional Flows , 1999 .

[8]  Messaoud Benantar.  Parallel and adaptive algorithms for elliptic partial differential systems , 1992 .

[9]  Hanan Samet,et al.  Applications of spatial data structures , 1989 .

[10]  Deborah Greaves,et al.  A moving boundary finite element method for fully nonlinear wave simulations , 1997 .

[11]  James M. Hyman,et al.  Numerical methods for tracking interfaces , 1984 .

[12]  P. Woodward,et al.  SLIC (Simple Line Interface Calculation) , 1976 .

[13]  D. P. Young,et al.  A locally refined rectangular grid finite element method: application to computational fluid dynamics and computational physics , 1990 .

[14]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

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

[16]  Deborah Greaves,et al.  Numerical simulation of fluid flows using an unstructured finite volume method with adaptive tri‐tree grids , 2002 .

[17]  Mark S. Shephard,et al.  Automatic three‐dimensional mesh generation by the finite octree technique , 1984 .

[18]  Hanan Samet,et al.  The Design and Analysis of Spatial Data Structures , 1989 .

[19]  Csaba Gáspár,et al.  Fast, adaptive approximation of wind-induced horizontal flow pattens in shallow lakes using quadtree-based multigrid method , 1992 .

[20]  D. D. Zeeuw,et al.  An adaptively refined Cartesian mesh solver for the Euler equations , 1993 .

[21]  Derek M. Causon.  An efficient front-tracking algorithm for multi-component flow calculations with biomedical applications , 1996 .

[22]  Mark Yerry,et al.  A Modified Quadtree Approach To Finite Element Mesh Generation , 1983, IEEE Computer Graphics and Applications.

[23]  O. Ubbink.  Numerical prediction of two fluid systems with sharp interfaces , 1997 .

[24]  Deborah Greaves,et al.  Hierarchical tree-based finite element mesh generation , 1999 .

[25]  Nasser Ashgriz,et al.  FLAIR: fluz line-segment model for advection and interface reconstruction , 1991 .

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

[27]  M. S Carvalho,et al.  Regular ArticleThree-Dimensional Stability Analysis of Free Surface Flows: Application to Forward Deformable Roll Coating☆ , 1999 .

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

[29]  B. P. Leonard,et al.  The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection , 1991 .