Quantitative model for predicting the imbibition dynamics of viscoelastic fluids in nonuniform microfluidic assays

We develop a mathematical model to quantitatively describe the imbibition dynamics of an elastic non-Newtonian fluid in a conical (nonuniform cross section) microfluidic assay. We consider the simplified Phan-Thien-Tanner viscoelastic model to represent the rheology of the elastic non-Newtonian fluid. Our model accounts for the geometrical features of the fluidic assay, the key parameters affecting the rheological behavior of the fluid, and predicts the imbibition dynamics effectively. By demonstrating the temporal advancement of the filling length in the conical capillary graphically, obtained for pertinent parametric values belonging to their physically permissible range, we report an underlying balance between capillary and viscous forces during imbibition resulting in three distinct regimes of filling. Nonuniformity in the capillary cross section gives rise to an alteration in the viscous force being applied at the contact line (manifested through the alteration in shear rate) during the imbibition process, which upon maintaining a balance with the dominant capillary force results in three different regimes of filling. We believe that the present analysis has a twofold significance. First, this work will enhance the understanding of underlying imbibition dynamics of viscoelastic fluids (most of the biofluids exhibit viscoelastic rheology) in nonuniform fluidic pathways. Second, the developed model is of significant practical relevance for the optimum design of microfluidic assays, primarily used for sample diagnostics in biochemical and biomedical applications.