A coupled quadrilateral grid level set projection method applied to ink jet simulation

A coupled level set-projection method on quadrilateral grids is developed for piezoelectric ink jet simulations. The model is based on the Navier-Stokes equations for incompressible two-phase flows in the presence of surface tension and density jump across the interface separating ink and air, coupled to an electric circuit model which describes the driving mechanism behind the process, and a macroscopic contact model which describes the air-ink-wall dynamics. We simulate the axisymmetric flow on quadrilateral grids using a combination of second-order finite difference projection methods to solve the fluid equations and level set methods to track the air/ink interface. To improve the mass conservation performance of the coupled level set method, a bicubic interpolation is combined with the Fast Marching Method for level set re-initialization on quadrilateral grids. The numerical method is used to analyze the motion of the interface, droplet pinch off, formation of satellites, effect of nozzle geometry on droplet size and motion, and the dynamics for droplet landing. The simulations are faithful to the dimensions and physics of a particular class of inkjet devices.

[1]  John B. Bell,et al.  A Numerical Method for the Incompressible Navier-Stokes Equations Based on an Approximate Projection , 1996, SIAM J. Sci. Comput..

[2]  A. Chorin Numerical solution of the Navier-Stokes equations , 1968 .

[3]  P. Colella A Direct Eulerian MUSCL Scheme for Gas Dynamics , 1985 .

[4]  Phillip Colella,et al.  A Projection Method for Incompressible Viscous Flow on Moving Quadrilateral Grids , 2001 .

[5]  J. Tsitsiklis Efficient algorithms for globally optimal trajectories , 1995, IEEE Trans. Autom. Control..

[6]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

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

[8]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[9]  Peter Smereka,et al.  Axisymmetric free boundary problems , 1997, Journal of Fluid Mechanics.

[10]  John B. Bell,et al.  Approximate Projection Methods: Part I. Inviscid Analysis , 2000, SIAM J. Sci. Comput..

[11]  Michael Griebel,et al.  Numerical Simulation in Fluid Dynamics: A Practical Introduction , 1997 .

[12]  John B. Bell,et al.  A projection method for viscous incom-pressible flow on quadrilateral grids , 1961 .

[13]  J A Sethian,et al.  Computing geodesic paths on manifolds. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[14]  E. Puckett,et al.  A High-Order Projection Method for Tracking Fluid Interfaces in Variable Density Incompressible Flows , 1997 .

[15]  P. Colella,et al.  An Adaptive Level Set Approach for Incompressible Two-Phase Flows , 1997 .

[16]  J. Sethian,et al.  Fast methods for the Eikonal and related Hamilton- Jacobi equations on unstructured meshes. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[17]  John B. Bell,et al.  A Cartesian Grid Projection Method for the Incompressible Euler Equations in Complex Geometries , 1994, SIAM J. Sci. Comput..

[18]  J. Sethian,et al.  Projection methods coupled to level set interface techniques , 1992 .

[19]  Alexander Vladimirsky,et al.  Ordered Upwind Methods for Static Hamilton-Jacobi Equations: Theory and Algorithms , 2003, SIAM J. Numer. Anal..

[20]  J. Sethian AN ANALYSIS OF FLAME PROPAGATION , 1982 .

[21]  S. Osher,et al.  A Level Set Formulation of Eulerian Interface Capturing Methods for Incompressible Fluid Flows , 1996 .

[22]  J. Sethian Curvature and the evolution of fronts , 1985 .

[23]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[24]  P. Colella Multidimensional upwind methods for hyperbolic conservation laws , 1990 .

[25]  J. Bell,et al.  A Second-Order Projection Method for Variable- Density Flows* , 1992 .

[26]  S. Sakai,et al.  Dynamics of Piezoelectric Inkjet Printing Systems. , 2000 .

[27]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[28]  J. Sethian,et al.  A coupled level set projection method applied to ink jet simulation , 2003 .

[29]  David L. Chopp,et al.  Some Improvements of the Fast Marching Method , 2001, SIAM J. Sci. Comput..