This book deals with a boundary element formulation employing a penalty function technique for viscous flow and thermal convection problems. The self-adaptive coordinate transformation technique is employed for computing efficiently and accurately quasi-integral boundary integrals. Some numerical solutions are shown and compared with previous published results to demonstrate the feasibility of the proposed formulation.