Bubble damping of non-stationary oscillatory flow stabilization in microfluidic systems

TL;DR

Oscillatory flows in microfluidics are prone to instability due to non-stationary forcing. The authors derive a first-principles model that couples air compressibility (via $\mathcal{C}=V_A/K$) with viscous resistance (via $\mathcal{R}$) to predict the transmitted liquid-front motion under oscillation; the model reveals that the dimensionless product $\tau\omega$ governs amplitude reduction, phase shift, and transient drift. The solution yields $V_{front}(t)$ and $Q(t)$ expressions, with amplitude ratio $V_{front}/V_0 = 1/\sqrt{1+(\tau\omega)^2}$ and phase $\phi = \arctan(\tau\omega)$, and includes a decaying transient term. The theory is validated across a range of air volumes, frequencies, and tube lengths using a custom oscillatory syringe-pump setup, establishing a predictive design framework. Practically, it turns trapped air into a tunable design element for stabilizing oscillatory microfluidic flows.

Abstract

The inherent instability of oscillatory flows presents a significant challenge in microfluidics, impairing performance in different applications from particle detachemnt to organs-on-a-chip. Trapped air inside a microfluidic system passively dampens these fluctuations because of the compressible nature of air. However, a foundational theoretical model that describes this effect has remained elusive. Here, a first-principles model that fully characterizes the effects of a trapped air volume in oscillatory microfluidic flow is derived. The model identifies a dimensionless product as the governing parameter, unifying the interplay between air compressibility and fluidic resistance. It precisely predicts the volume displacement dynamics of the liquid front, which compared with the original flow, it presents amplitude reduction, phase shift, and transient drift. The theoretical framework was validated with different experiments across a broad range of conditions. This work transforms trapped air from a source of unpredictability into a powerful, predictable element for tailoring oscillatory flow stability, providing a rigorous design tool for microfluidic systems.

Figures (6)

• Figure 1: Schematic of the oscillatory flow system with a trapped air volume. The piston undergoes sinusoidal motion, compressing and expanding the air volume, which in turn drives the liquid through a narrow tube with significant fluidic resistance.
• Figure 2: Custom oscillatory syringe pump. (A) Schematic diagram of the crank mechanism converting rotary motor motion to sinusoidal plunger displacement. (B) Photograph of the manually fabricated prototype showing the motor, reduction gear, connecting rods, and syringe assembly.
• Figure 3: Schematic of the experimental setup for validating the damping model. The system consists of the oscillatory pump connected to a flexible tube, a glass observation tube containing the liquid-air interface, and a narrow capillary tube that provides dominant fluidic resistance, ending in a reservoir maintained at atmospheric pressure.
• Figure 4: Experimental time-series of the volume displaced by the piston, $V_p(t)$ (blue, Eq. \ref{['eq:V_piston']}), and the liquid front, $V(t)$ (orange, Eq. \ref{['eq:V_solution_final']}), for three systems with varying $\tau\omega$. Experimental data for the liquid front are shown as points. System parameters: top: $\omega=0.37\ \mathrm{s}^{-1}$, $\mathcal{R}=2.94\times 10^{12}\ \mathrm{Pa\ s/m}^3$; middle: $\omega=1.37\ \mathrm{s}^{-1}$, $\mathcal{R}=1.96\times 10^{12}\ \mathrm{Pa\ s/m}^3$; bottom: $\omega=8.37\ \mathrm{s}^{-1}$, $\mathcal{R}=1.47\times 10^{12}\ \mathrm{Pa\ s/m}^3$; all with $\mathcal{C}=5.94\times 10^{-12}\ \mathrm{m}^3/\mathrm{Pa}$.
• Figure 5: Amplitude relation between the plunger motion and the fluid front as a function of the $\tau \omega$ product for each experiment. The theoretical amplitude ratio obtained in equation \ref{['eq:amplitude_ratio']} is plotted as a solid orange line.