This work is intended to investigate both T-shaped micromixers and interdigital micromixers by using CFD simulations. In the first part of this work, the influences of design and operation parameters upon mixing, heat transfer, chemical reaction, and residence time distribution are investigated. The simulation results show that the convective flow in T-shaped micromixers is effective in increasing heat transfer to realize uniform temperature distribution. In addition, the effects of convective flow are highly dependent on design and operation parameters. In the second part of this work, a systematic design method is proposed. Namely, response surface method (RSM) is introduced into the CFD-based optimal design of micromixers. RSM is one of effective statistical analysis approaches and is often used in the product design works based on experimental data. The validity of the proposed design method is assessed through case studies. The combination of CFD and RMS promotes the design efficiency of various types of micromixer.