A finite difference level set-projection method on rectangular grid is developed for piezoelectric ink jet simulation. 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 using a combination of second order projection methods to solve the fluid equations and level set methods to track the air/ink interface. The numerical method can be used to analyze the motion of the interface, breakoff and formation of satellites, and effect of nozzle geometry on droplet size and motion. We focus on close comparison of our numerical ink jet simulation with experimental data.

A coupled finite difference algorithm on rectangular grids is developed for viscoelastic ink ejection simulations. The ink is modeled by the Oldroyd-B viscoelastic fluid model. The coupled algorithm seamlessly incorporates several things: (1) a coupled level set-projection method for incompressible immiscible two-phase fluid flows; (2) a higher-order Godunov type algorithm for the convection terms in the momentum and level set equations; (3) a simple first-order upwind algorithm for the convection term in the viscoelastic stress equations; (4) central difference approximations for viscosity, surface tension, and upper-convected derivative terms; and (5) an equivalent circuit model to calculate the inflow pressure (or flow rate) from dynamic voltage.

Two-dimensional first-order governing equations for electroded piezoelectric crystal plates with general symmetry and thickness-graded material properties are deduced from the three-dimensional equations of linear piezoelectricity by Mindlin's general procedure of series expansion. Mechanical displacements and thickness-graded material properties, i.e., the elastic stiffnesses, piezoelectric coefficients, dielectric permittivities, and mass density, are expanded in powers of the thickness coordinate, while electric potential is expanded in a special series in order to accommodate the specified electric potentials at electroded faces of the plate. The effects of graded material properties on the piezoelectrically induced stresses or deformations by the applied surface potentials are clearly exhibited in these newly derived equations which reduce to Mindlin's first-order equations of elastic anisotropic plates when the material properties are homogeneous. Closed form solutions are obtained from the three-dimensional equations of piezoelectricity and from the present two-dimensional equations for both homogeneous plates and bimorphs of piezoelectric ceramics. Dispersion curves for homogeneous plates and bimorphs and resonance frequencies for bimorph strips with finite width are computed from the solutions of three-dimensional and two-dimensional equations. Comparison of the results shows that predictions from the two-dimensional equations are very close to those from the three-dimensional equations.

No abstract

No abstract