# Numerical simulation of 2D electrothermal flow using boundary element method

Qinlong Ren, Cholik Chan, Alberto L. Arvayo

Research output: Chapter in Book/Report/Conference proceedingConference contribution

2 Citations (Scopus)

### Abstract

Microfluidics and its applications to Lab-on-a-Chip have attracted a lot of attention. Because of the small length scale, the flow is characterized by a low Re number. The governing equations become linear. Boundary element method (BEM) is a very good option for simulating the fluid flow with high accuracy. In this paper, we present a 2D numerical modeling of the electrothermal flow using BEM. In electrothermal flow the volumetric force is caused by electric field and temperature gradient. The physics is mathematically modeled by (i) Laplace equation for the electrical potential, (ii) Poisson equation for the heat conduction caused by Joule heating, (iii) continuity and Stokes equation for the low Reynolds number flow. We begin by solving the electrical potential and electric field. The heat conduction is caused by the Joule heating as the heat generation term. Superposition principle is used to solve for the temperature field. The Coulomb and dielectric forces are generated by the electrical field and temperature gradient of the system. We analyze the Stokes flow problem by superposition of fundamental solution for free-space velocity caused by body force and BEM for the corresponding homogeneous Stokes equation. It is well known that a singularity integral arises when the source point approaches the field point. To overcome this problem, we solve the free-space velocity analytically. For the BEM part, we also calculate all the integral terms analytically. With this effort, our solution is more accurate. In addition, we improve the robustness of the matrix system by combining the velocity integral equation with the traction integral equation. Our purpose is to design a pump for the microfluidics system. Since the system is a long channel, the flow is fully developed in the area far away from the electrodes. With this assumption, the velocity profile is parabolic at the inlet and outlet of the channel. So we can get appropriate boundary conditions for the BEM part of Stokes equation. Consequently, we can simulate the electrothermal flow in an open channel. In this paper, we will present the formulation and implementation of BEM to model electrothermal flow. Results of electrical potential, temperature field, Joule heating, electrothermal force, and velocity field will be presented.

Original language English (US) ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013 American Society of Mechanical Engineers (ASME) https://doi.org/10.1115/MNHMT2013-22075 Published - 2013 ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013 - Hong Kong, ChinaDuration: Dec 11 2013 → Dec 14 2013

### Other

Other ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013 China Hong Kong 12/11/13 → 12/14/13

### Fingerprint

Boundary element method
Joule heating
Computer simulation
Microfluidics
Heat conduction
Integral equations
Temperature distribution
Electric fields
Lab-on-a-chip
Laplace equation
Poisson equation
Heat generation
Linear equations
Flow of fluids
Reynolds number
Physics
Boundary conditions
Pumps
Electrodes

### ASJC Scopus subject areas

• Fluid Flow and Transfer Processes

### Cite this

Ren, Q., Chan, C., & Arvayo, A. L. (2013). Numerical simulation of 2D electrothermal flow using boundary element method. In ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013 American Society of Mechanical Engineers (ASME). https://doi.org/10.1115/MNHMT2013-22075

Numerical simulation of 2D electrothermal flow using boundary element method. / Ren, Qinlong; Chan, Cholik; Arvayo, Alberto L.

ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013. American Society of Mechanical Engineers (ASME), 2013.

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Ren, Q, Chan, C & Arvayo, AL 2013, Numerical simulation of 2D electrothermal flow using boundary element method. in ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013. American Society of Mechanical Engineers (ASME), ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013, Hong Kong, China, 12/11/13. https://doi.org/10.1115/MNHMT2013-22075
Ren Q, Chan C, Arvayo AL. Numerical simulation of 2D electrothermal flow using boundary element method. In ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013. American Society of Mechanical Engineers (ASME). 2013 https://doi.org/10.1115/MNHMT2013-22075
Ren, Qinlong ; Chan, Cholik ; Arvayo, Alberto L. / Numerical simulation of 2D electrothermal flow using boundary element method. ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013. American Society of Mechanical Engineers (ASME), 2013.
@inproceedings{1696a5f0ec9c49d3afa3d19a4ffca5c0,
title = "Numerical simulation of 2D electrothermal flow using boundary element method",
abstract = "Microfluidics and its applications to Lab-on-a-Chip have attracted a lot of attention. Because of the small length scale, the flow is characterized by a low Re number. The governing equations become linear. Boundary element method (BEM) is a very good option for simulating the fluid flow with high accuracy. In this paper, we present a 2D numerical modeling of the electrothermal flow using BEM. In electrothermal flow the volumetric force is caused by electric field and temperature gradient. The physics is mathematically modeled by (i) Laplace equation for the electrical potential, (ii) Poisson equation for the heat conduction caused by Joule heating, (iii) continuity and Stokes equation for the low Reynolds number flow. We begin by solving the electrical potential and electric field. The heat conduction is caused by the Joule heating as the heat generation term. Superposition principle is used to solve for the temperature field. The Coulomb and dielectric forces are generated by the electrical field and temperature gradient of the system. We analyze the Stokes flow problem by superposition of fundamental solution for free-space velocity caused by body force and BEM for the corresponding homogeneous Stokes equation. It is well known that a singularity integral arises when the source point approaches the field point. To overcome this problem, we solve the free-space velocity analytically. For the BEM part, we also calculate all the integral terms analytically. With this effort, our solution is more accurate. In addition, we improve the robustness of the matrix system by combining the velocity integral equation with the traction integral equation. Our purpose is to design a pump for the microfluidics system. Since the system is a long channel, the flow is fully developed in the area far away from the electrodes. With this assumption, the velocity profile is parabolic at the inlet and outlet of the channel. So we can get appropriate boundary conditions for the BEM part of Stokes equation. Consequently, we can simulate the electrothermal flow in an open channel. In this paper, we will present the formulation and implementation of BEM to model electrothermal flow. Results of electrical potential, temperature field, Joule heating, electrothermal force, and velocity field will be presented.",
author = "Qinlong Ren and Cholik Chan and Arvayo, {Alberto L.}",
year = "2013",
doi = "10.1115/MNHMT2013-22075",
language = "English (US)",
booktitle = "ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013",
publisher = "American Society of Mechanical Engineers (ASME)",

}

TY - GEN

T1 - Numerical simulation of 2D electrothermal flow using boundary element method

AU - Ren, Qinlong

AU - Chan, Cholik

AU - Arvayo, Alberto L.

PY - 2013

Y1 - 2013

N2 - Microfluidics and its applications to Lab-on-a-Chip have attracted a lot of attention. Because of the small length scale, the flow is characterized by a low Re number. The governing equations become linear. Boundary element method (BEM) is a very good option for simulating the fluid flow with high accuracy. In this paper, we present a 2D numerical modeling of the electrothermal flow using BEM. In electrothermal flow the volumetric force is caused by electric field and temperature gradient. The physics is mathematically modeled by (i) Laplace equation for the electrical potential, (ii) Poisson equation for the heat conduction caused by Joule heating, (iii) continuity and Stokes equation for the low Reynolds number flow. We begin by solving the electrical potential and electric field. The heat conduction is caused by the Joule heating as the heat generation term. Superposition principle is used to solve for the temperature field. The Coulomb and dielectric forces are generated by the electrical field and temperature gradient of the system. We analyze the Stokes flow problem by superposition of fundamental solution for free-space velocity caused by body force and BEM for the corresponding homogeneous Stokes equation. It is well known that a singularity integral arises when the source point approaches the field point. To overcome this problem, we solve the free-space velocity analytically. For the BEM part, we also calculate all the integral terms analytically. With this effort, our solution is more accurate. In addition, we improve the robustness of the matrix system by combining the velocity integral equation with the traction integral equation. Our purpose is to design a pump for the microfluidics system. Since the system is a long channel, the flow is fully developed in the area far away from the electrodes. With this assumption, the velocity profile is parabolic at the inlet and outlet of the channel. So we can get appropriate boundary conditions for the BEM part of Stokes equation. Consequently, we can simulate the electrothermal flow in an open channel. In this paper, we will present the formulation and implementation of BEM to model electrothermal flow. Results of electrical potential, temperature field, Joule heating, electrothermal force, and velocity field will be presented.

AB - Microfluidics and its applications to Lab-on-a-Chip have attracted a lot of attention. Because of the small length scale, the flow is characterized by a low Re number. The governing equations become linear. Boundary element method (BEM) is a very good option for simulating the fluid flow with high accuracy. In this paper, we present a 2D numerical modeling of the electrothermal flow using BEM. In electrothermal flow the volumetric force is caused by electric field and temperature gradient. The physics is mathematically modeled by (i) Laplace equation for the electrical potential, (ii) Poisson equation for the heat conduction caused by Joule heating, (iii) continuity and Stokes equation for the low Reynolds number flow. We begin by solving the electrical potential and electric field. The heat conduction is caused by the Joule heating as the heat generation term. Superposition principle is used to solve for the temperature field. The Coulomb and dielectric forces are generated by the electrical field and temperature gradient of the system. We analyze the Stokes flow problem by superposition of fundamental solution for free-space velocity caused by body force and BEM for the corresponding homogeneous Stokes equation. It is well known that a singularity integral arises when the source point approaches the field point. To overcome this problem, we solve the free-space velocity analytically. For the BEM part, we also calculate all the integral terms analytically. With this effort, our solution is more accurate. In addition, we improve the robustness of the matrix system by combining the velocity integral equation with the traction integral equation. Our purpose is to design a pump for the microfluidics system. Since the system is a long channel, the flow is fully developed in the area far away from the electrodes. With this assumption, the velocity profile is parabolic at the inlet and outlet of the channel. So we can get appropriate boundary conditions for the BEM part of Stokes equation. Consequently, we can simulate the electrothermal flow in an open channel. In this paper, we will present the formulation and implementation of BEM to model electrothermal flow. Results of electrical potential, temperature field, Joule heating, electrothermal force, and velocity field will be presented.

UR - http://www.scopus.com/inward/record.url?scp=84901754977&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84901754977&partnerID=8YFLogxK

U2 - 10.1115/MNHMT2013-22075

DO - 10.1115/MNHMT2013-22075

M3 - Conference contribution

AN - SCOPUS:84901754977

BT - ASME 2013 4th International Conference on Micro/Nanoscale Heat and Mass Transfer, MNHMT 2013

PB - American Society of Mechanical Engineers (ASME)

ER -