0
Research Papers

Algebraic Dimensional Reduction for Microfluidic Simulation

[+] Author and Article Information
Josh Danczyk

Department of Mechanical Engineering, University of Wisconsin Madison, 2059 Mechanical Engineering Building, 1513 University Avenue, Madison, WI 53706-1572jrdanczyk@wisc.edu

Krishnan Suresh

Department of Mechanical Engineering, University of Wisconsin Madison, 2059 Mechanical Engineering Building, 1513 University Avenue, Madison, WI 53706-1572suresh@engr.wisc.edu

J. Comput. Inf. Sci. Eng 9(3), 031001 (Aug 04, 2009) (6 pages) doi:10.1115/1.3184590 History: Received November 02, 2007; Revised October 16, 2008; Published August 04, 2009

Microfluidic devices exhibit a high-aspect ratio in that their channel-widths are much smaller than their overall lengths. High-aspect geometry leads to an unduly large finite element mesh, making the (otherwise popular) finite element method (FEM) a poor choice for modeling microfluidic devices. An alternate computational strategy is to exploit well-known analytical solutions for fluid flow over the narrow channels of a device, and then either (a) assume the same analytical solutions for the cross-flow regions, or (b) exploit these solutions to set-up artificial boundary conditions over the cross-flow regions. Such simplified models are computationally far superior to brute-force FEM, but do not support the generality or flexibility of FEM. In this paper, we propose a third strategy for exploiting the analytical solutions: (c) directly incorporate them into standard FE-based analysis via algebraic reduction techniques. The advantages of the proposed strategy are (1) designers can use standard computer-aided design/computer-aided engineering (CAD/CAE) environments to model, analyze, and postprocess microfluidic simulation; (2) well-established dual-weighted residuals can be used to estimate modeling errors; and (3), if desired, one can eliminate the dependency on analytical solutions over selected regions, and instead revert to brute-force FEM. The simplicity and generality of the proposed method is inherited from the model reduction process, so are its theoretical properties, while simultaneously its computational efficiency is inherited from the use of analytical solutions.

FIGURES IN THIS ARTICLE
<>
Copyright © 2009 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.

References

Figures

Grahic Jump Location
Figure 1

A step-channel microfluidic device

Grahic Jump Location
Figure 2

A FE mesh over the step-channel

Grahic Jump Location
Figure 3

A H-cell microfluidic device

Grahic Jump Location
Figure 4

The computed pressure field over the H-cell

Grahic Jump Location
Figure 5

Computed pressure over a line segment 5≤x≤119, y=0 μm

Grahic Jump Location
Figure 6

Accuracy versus computational effort

Grahic Jump Location
Figure 7

Computed concentration over two line segments 5≤x≤119, y=±1 μm

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In