The saturation equation - an analytical expression for partial saturation during wicking flow in paper microfluidic channels

Satvik Verma,Bhushan J. Toley

DOI: https://doi.org/10.26434/chemrxiv-2024-cznzf

2024-01-08

Abstract:The design and fabrication of paper based microfluidic devices is critically dependent on modelling fluid flow through porous paper membranes. A commonly observed phenomenon is partial saturation, i.e., regions of the paper membrane not being filled completely due to pores of different sizes. The most comprehensive model till date of partial saturation during wicking flow in paper is the Richard’s equation. However, the nonlinear nature of this equation and the requirement for numerical solvers for its solution make it largely inaccessible to the paper microfluidics and lateral flow assay community. Moreover, the parameters used in Richard’s equation often need to be tweaked (or assumed) to match experimental data. This necessitates the need for a simple and appropriate model of partial saturation in paper membranes, easily usable by the wider research community. In the current work, we present an approach to model paper membranes as a bundle of parallel capillaries whose radii follow a two parameter log normal distribution. Application of the Washburn equation to the bundle provides a distribution of fluid fronts, which can be used to calculate saturation. Using this approach, we developed the first analytical expression for spatiotemporal variation of saturation in 1D wicking flow. Experimentally obtained data for spatiotemporal saturation for four different paper materials were fit to this analytical model to obtain parameters for each material. Results obtained from this analytical model match both experimental data as well as numerical results obtained from the Richards equation model well. The availability of an analytical expression for partial saturation in wicking flow promises to significantly increase access to such modelling among the wider research community.

Chemistry

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the partial saturation phenomenon during the capillary - action - flow process in paper - based microfluidic channels. Specifically, when the liquid flows in the paper membrane, not all pores will be completely filled, which leads to the partial saturation phenomenon. This phenomenon is crucial for the design and manufacturing of paper - based microfluidic devices, because partial saturation will affect the performance of the devices, especially in rapid - diagnosis tests and point - of - care testing applications.

Existing models such as the Lucas Washburn equation and Darcy's law cannot accurately describe the partial saturation phenomenon, and although the more complex Richard equation can simulate this phenomenon, it requires the use of a numerical solver (such as COMSOL), which makes it difficult for many researchers to use these tools. Therefore, researchers need a simple and easy - to - use model to describe the partial saturation phenomenon in the paper membrane.

To solve this problem, this paper proposes a new method, that is, by modeling the paper membrane as a bundle of parallel capillaries, the radii of these capillaries follow a two - parameter log - normal distribution. Through the application of this model, the researchers developed the "saturation equation" - an explicit analytical expression for calculating the change of saturation with space and time during one - dimensional capillary - action - flow. This model can not only accurately predict experimental data, but also is consistent with the numerical results obtained using the Richard equation, thus providing important theoretical support for the design of paper - based microfluidic devices.

### Formula Summary
1. **Lucas Washburn Equation**:
\[
L^2=\frac{\gamma r_{\text{pore}}t}{2\nu}
\]
where \( r_{\text{pore}} \) is the effective radius of the pore, \( \gamma \) is the surface tension of the liquid - air interface (including the contact - angle - dependent term), \( \nu \) is the dynamic viscosity of the liquid, and \( L \) is the distance that the liquid travels in the paper within time \( t \).

2. **Log - normal Distribution**:
\[
f(r)=\frac{1}{r\sigma\sqrt{2\pi}}\exp\left(-\frac{(\ln(r)-\mu)^2}{2\sigma^2}\right)
\]
where \( \mu \) is the log - median and \( \sigma \) is the log - standard deviation.

3. **Saturation Equation**:
\[
\theta(x,t)=\frac{1}{2}\text{erfc}\left(\frac{\ln(10^6x^2)-\mu-\ln(A^2t)-\frac{5}{2}\sigma^2}{\sqrt{2}\sigma}\right)
\]
where \( A = \sqrt{\frac{\gamma}{2\nu}} \), and \( \text{erfc} \) is the complementary error function.

Through these formulas, the researchers successfully provided a simple and accurate model for describing the partial saturation phenomenon in paper - based microfluidic devices.