The fundamental action of the bubble-driven inertial micropump is investigated. The pump has no moving parts and consists of a thermal resistor placed asymmetrically within a straight channel connecting two reservoirs. Using numerical simulations, the net flow is studied as a function of channel geometry, resistor location, vapor bubble strength, fluid viscosity, and surface tension. Two major regimes of behavior are identified: axial and non-axial. In the axial regime, the drive bubble either remains inside the channel or continues to grow axially when it reaches the reservoir. In the non-axial regime the bubble grows out of the channel and in all three dimensions while inside the reservoir. The net flow in the axial regime is parabolic with respect to the hydraulic diameter of the channel cross-section but in the non-axial regime it is not. From numerical modeling, it is determined that the net flow is maximal when the axial regime crosses over to the non-axial regime. To elucidate the basic physical principles of the pump, a phenomenological one-dimensional model is developed and solved. A linear array of micropumps has been built using silicon-SU8 fabrication technology, and semi-continuous pumping across a 2 mm-wide channel has been demonstrated experimentally. Measured variation of the net flow with fluid viscosity is in excellent agreement with simulation results.

A one-dimensional model of inertial pumping is introduced and solved. The pump is driven by a high-pressure vapor bubble generated by a microheater positioned asymmetrically in a microchannel. The bubble is approximated as a short-term impulse delivered to the two fluidic columns inside the channel. Fluid dynamics is described by a Newton-like equation with a variable mass, but without the mass derivative term. Because of smaller inertia, the short column refills the channel faster and accumulates a larger mechanical momentum. After bubble collapse the total fluid momentum is nonzero, resulting in a net flow. Two different versions of the model are analyzed in detail, analytically and numerically. In the symmetrical model, the pressure at the channel-reservoir connection plane is assumed constant, whereas in the asymmetrical model it is reduced by a Bernoulli term. For low and intermediate vapor bubble pressures, both models predict the existence of an optimal microheater location. The predicted net flow in the asymmetrical model is smaller by a factor of about 2. For unphysically large vapor pressures, the asymmetrical model predicts saturation of the effect, while in the symmetrical model net flow increases indefinitely. Pumping is reduced by nonzero viscosity, but to a different degree depending on the microheater location.

Summary In this study, we present a method to predict the droplet ejection in thermal inkjet printheads including the growth and collapse of a vapor bubble and refill of the firing chamber. The three‐dimensional Navier–Stokes equations are solved using a finite‐volume approach with a fixed Cartesian mesh. The piecewise‐linear interface calculation‐based volume‐of‐fluid method is employed to track and reconstruct the ink–air interface. A geometrical computation based on Lagrangian advection is used to compute the mass flux and advance the interface. A simple and efficient model for the bubble dynamics is employed to model the effect of ink vapor on the adjacent ink liquid. To solve the surface tension‐dominated flow accurately, a hierarchical curvature‐estimation method is proposed to adapt to the local grid resolution. The numerical methods mentioned earlier have been implemented in an internal simulation code, CFD3. The numerical examples presented in the study show good performance of CFD3 in prediction of surface tension‐dominated free‐surface flows, for example, droplet ejection in thermal inkjet printing. Currently, CFD3 is used extensively for printhead development within Hewlett‐Packard. Copyright © 2015 John Wiley & Sons, Ltd.

Bubble-driven inertial pumps are a novel method of moving liquids through microchannels. We combine high-speed imaging, computational fluid dynamics (CFD) simulations and an effective one-dimensional model to study the fundamentals of inertial pumping. For the first time, single-pulse transient flow through Ushaped microchannels is imaged over the entire pump cycle with 4-µs temporal resolution. Observations confirm the fundamental N-shape flow profile predicted earlier by theory and simulations. Experimental flow rates are used to calibrate the CFD and one-dimensional models to extract an effective bubble strength. Then the frequency dependence of inertial pumping is studied both experimentally and numerically. The pump efficiency is found to gradually decrease once the successive pulses start to overlap in time.