Accelerated Computational Fluid Dynamics Simulations of Microfluidic Devices by Exploiting Higher Levels of Abstraction

Michel Takken,Robert Wille

2024-01-15

Abstract:The design of microfluidic devices is a cumbersome and tedious process that can be significantly improved by simulation. Methods based on Computational Fluid Dynamics (CFD) are considered state-of-the-art, but require extensive compute time - oftentimes limiting the size of microfluidic devices that can be simulated. Simulation methods that abstract the underlying physics on a higher level generally provide results instantly, but the fidelity of these methods is usually worse. In this work, a simulation method that accelerates CFD simulations by exploiting simulation methods on higher levels of abstraction is proposed. Case studies confirm that the proposed method accelerates CFD simulations by multiple factors (often several orders of magnitude) while maintaining the fidelity of CFD simulations.

Computational Engineering, Finance, and Science,Fluid Dynamics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the excessive time - consuming of computational fluid dynamics (CFD) simulation in the design of microfluidic devices. Specifically, although the current CFD methods used for microfluidic device simulation can provide high - precision results, their calculation time is very long, sometimes even taking several days or weeks, even on dedicated workstations. This limits the wide application of CFD simulation in the design of microfluidic devices, especially in cases where the entire device needs to be simulated to understand the interactions between various components.

To overcome this problem, the author proposes a method to accelerate CFD simulation. By using a higher - level abstract simulation method, the calculation time is significantly reduced while maintaining the accuracy comparable to that of CFD simulation. The core idea of this method is to use a high - level abstract simulation method in areas where the flow field is highly ordered, and use the traditional CFD method in areas where the flow field is complex and difficult to predict, and ensure the continuity and consistency between the two methods through an iterative method.

### Main contributions of the paper
1. **Accelerating CFD simulation**: By combining the high - level abstract simulation method and the traditional CFD method, rapid simulation of microfluidic devices is achieved, and the acceleration factor usually reaches several orders of magnitude.
2. **Maintaining accuracy**: Although the calculation time is greatly reduced, this method can still maintain the simulation accuracy comparable to that of the traditional CFD method.
3. **Generality**: This method is not only applicable to continuous - channel microfluidic devices, but can also be extended to other types of microfluidic platforms, such as droplet - based microfluidics and paper - based microfluidics.

### Implementation methods
1. **Region division**: First, divide the network of the microfluidic device into two parts: the low - abstraction - level region (Ωlow) and the high - abstraction - level region (Ωhigh). The low - abstraction - level region contains complex flow field structures, such as intersections and T - joints, which require the use of the CFD method for accurate simulation; the high - abstraction - level region contains simple flow field structures, such as straight channels, and can use the high - level abstract method for rapid simulation.
2. **Simulation methods**: Use the modified nodal analysis (MNA) method in the high - abstraction - level region and the lattice Boltzmann method (LBM) in the low - abstraction - level region.
3. **Communication mechanism**: Pass the pressure and velocity field information between the high - abstraction - level region and the low - abstraction - level region through an iterative method to ensure that the solution of the entire network is continuous and consistent.

### Case studies and evaluation results
The author conducted case studies through four different continuous - channel microfluidic networks to verify the effectiveness of the proposed method. The results show that compared with the traditional CFD method, the calculation time of this method is significantly reduced, and the acceleration factor ranges from dozens of times to more than 1,000 times, while maintaining a relatively high simulation accuracy.

### Formula summary
- **Hagen - Poiseuille law**:
\[
\Delta p = Q\cdot R_H
\]
where \(\Delta p\) is the pressure drop, \(Q\) is the flow rate, and \(R_H\) is the hydraulic resistance.

- **Washburn equation**:
\[
L^2=\frac{\gamma Dt}{4\mu}
\]
where \(L\) is the moving distance of the liquid front, \(\gamma\) is the effective surface tension, \(D\) is the diffusion coefficient, \(t\) is the time, and \(\mu\) is the dynamic viscosity.

- **Navier - Stokes equations**:
- Mass conservation equation:
\[
\nabla\cdot\mathbf{u} = 0
\]
- Momentum conservation equation:
\[
\frac{\partial\mathbf{u}}{\partial t}+\nabla\cdot(\mathbf{u}\mathbf{u}^T)=-\frac{1}{\rho}\nabla p+\nu\nabla^2\mathbf{u}
\]
where \(\mathbf{u}\) is the flow velocity vector, \(\rho\) is the fluid density, \(p\) is the pressure, and \(\nu\) is the kinematic viscosity.

Through these methods and formulas, the author has successfully solved the design of microfluidic devices.