Carbon dioxide bubble removal at the anode of a direct methanol fuel cell (DMFC) is an important technique especially for applications in the portable power sources. This paper presents numerical investigations of the two-phase flow, CO2 bubbles in a liquid methanol solution, in the anode microchannels from the aspect of microfluidics using a thermal lattice Boltzmann model (TLBM). The main purpose is to derive an efficient and effective computational scheme to deal with this technical problem. It is then examined by a commercially available software using Navier-Stokes plus volume of fluid (VOF) method. The latter approach is normally employed by most researchers. A simplified microchannel simulation domain with the dimension of 1.5 mu m in height (or width) and 16.0 mu m in length has been setup for both cases to mimic the actual flow path of a CO2 bubble inside an anodic diffusion layer in the DMFC. This paper compares both numerical schemes and results under the same operation conditions from the viewpoint of fuel cell engineering. (C) 2009 Elsevier B.V. All rights reserved.