A numerical study of 2D electrothermal flow using boundary element method

The electrothermal flow phenomena can be applied to many microfluidic devices such as lab-on-a-chip. As a result of the small length scale in these devices, the fluid flow is characterized by a low Reynolds number thus allowing the governing equations to become linear. In this paper, a 2D numerical modeling of the electrothermal flow using boundary element method (BEM) is presented. BEM is an advantageous option for simulating the electrothermal flow. In an electrothermal flow, the volumetric body force depends on the electric field and temperature gradient. The physics is mathematically modeled by (i) Laplace’s equation for the electrical potential, (ii) Poisson’s equation for the heat conduction with Joule heating, and (iii) continuity and Stokes equations for the low Reynolds number flow. When using BEM to solve the equations, it is well known that a singular integral arises when the source point approaches the field point. Accurate evaluation of the singular integral is important to obtain an accurate simulation. To this end, all the singular and non-singular integrals are evaluated analytically. Consequently, an accurate algorithm is obtained. The formulation and implementation of BEM to model the electrothermal flow and the resulting electrical potential, temperature field, Joule heating and velocity field are presented in this paper.