Development of a grid-supported marker particle scheme for interfacetracking

A methodology is presented to simulate the growth and interaction of unstable fronts. Such fronts are found to be important in instabilities arising in several natural and industrial processes, such as solidification, spray dynamics, and bubble growth. The numerical simulation of such phenomena is challenging on account of the highly distorted moving boundary at which, often, curvature-dependent boundary conditions need be applied in each phase. Herein is presented a numerical technique to capture highly distorted interfaces. The interface is represented employing marker particles. Joining successive markers with circular arcs yields values of curvatures and normals on the interface. The markers are followed over an underlying cartesian grid and new marker particles are generated at each time step by an intersection procedure. The issue of mergers of interfaces is also attacked and the use of cells permits the simulation of merger-breakup processes. Thus, the method presented here, unlike previous marker particle based schemes, can handle particle depletion/accumulation and mergerbreakup issues with good accuracy. Results are presented, employing velocity functions modeled to mimic the actual instability phenomena, to demonstrate the accuracy and capability of the scheme developed. * Research Assistant ** Professor, Member AIAA Copyright @ 1993 by authors. Published by the American Institute of Aeronautics and Astronautics with permission.

[1]  Lyle H. Ungar,et al.  Cellular interface morphologies in directional solidification. The one-sided model , 1984 .

[2]  Martin E. Glicksman,et al.  Interaction of Flows with the Crystal-Melt Interface , 1986 .

[3]  Martin E. Glicksman,et al.  FUNDAMENTALS OF DENDRITIC SOLIDIFICATION—I. STEADY-STATE TIP GROWTH , 1981 .

[4]  Hideaki Miyata,et al.  Finite-difference simulation of breaking waves , 1986 .

[5]  Stephen H. Davis,et al.  Hydrodynamic interactions in directional solidification , 1990, Journal of Fluid Mechanics.

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

[7]  Bart J. Daly,et al.  A technique for including surface tension effects in hydrodynamic calculations , 1969 .

[8]  P. Tabeling,et al.  An experimental study of the Saffman-Taylor instability , 1987, Journal of Fluid Mechanics.

[9]  G. Taylor,et al.  The penetration of a fluid into a porous medium or Hele-Shaw cell containing a more viscous liquid , 1958, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[10]  Rosenberger,et al.  Morphological evolution of growing crystals: A Monte Carlo simulation. , 1988, Physical review. A, General physics.

[11]  L. Schwartz,et al.  A boundary-integral method for two-phase displacement in Hele-Shaw cells , 1986, Journal of Fluid Mechanics.

[12]  Liang Random-walk simulations of flow in Hele Shaw cells. , 1986, Physical review. A, General physics.

[13]  G. Goldbeck-Wood,et al.  Numerical simulation of dendritic growth. , 1988, Physical review. A, General physics.

[14]  P. Y. Liang,et al.  Numerical method for calculation of surface tension flows in arbitrary grids , 1991 .

[15]  P. Liang A numerical method for calculation of surface tension flows in arbitrary grids , 1988 .

[16]  David H. Sharp,et al.  Front Tracking Applied to Rayleigh–Taylor Instability , 1986 .

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

[18]  Jerzy M. Floryan,et al.  Numerical Methods for Viscous Flows With Moving Boundaries , 1989 .

[19]  Wei Shyy,et al.  An interface tracking method applied to morphological evolution during phase change , 1992 .

[20]  J. Langer Instabilities and pattern formation in crystal growth , 1980 .

[21]  Joe F. Thompson,et al.  Numerical grid generation , 1985 .

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

[23]  The development of solidification microstructures in the presence of lateral constraints , 1990 .

[24]  S. H. Tirmizi,et al.  Effect of natural convection on growth velocity and morphology of dendritic ice crystals , 1987 .

[25]  S. Forth,et al.  Coupled convective and morphological instability in a simple model of the solidification of a binary alloy, including a shear flow , 1992, Journal of Fluid Mechanics.

[26]  Nasser Ashgriz,et al.  A computational method for determining curvatures , 1989 .

[27]  B. J. Daly Numerical Study of Two Fluid Rayleigh‐Taylor Instability , 1967 .

[28]  J. Strain A boundary integral approach to unstable solidification , 1989 .

[29]  Alexandre J. Chorin,et al.  Curvature and solidification , 1985 .

[30]  Elaine S. Oran,et al.  Numerical Simulation of Reactive Flow , 1987 .

[31]  Gretar Tryggvason,et al.  Vortex dynamics of passive and active interfaces , 1984 .