Abstract: The aerodynamic performance of a wind turbine largely determines the safety and the efficiency. An efficient and accurate numerical simulation for the aerodynamics is a challenging task owing to its complexity. This paper presents a wind turbine modeling method based on immersion boundary method, including the whole process of the wind turbine modeling, the grid dispersion and the numerical simulation. The homotopy transform is used to generate a smooth blade model; the affine transform is used to deal with the taper and the twist of the airfoil. The numerical accuracy of the algorithm is tested for the lift and the drag of the two-dimensional airfoil. Based on the numerical results, the Richardson linear extrapolation is proposed to improve the numerical prediction of the drag substantially with only a mild increase of the computational complexity. The influence of the arch curvature and the thickness on the lift and drag forces of the two-dimensional airfoil is investigated. The power coefficients of a single wind turbine (including the tower) with different tip speed ratios are studied, along with the aerodynamic interaction between the tower and the blades. Finally, the aerodynamic interference of two wind turbines in-tandem with various separation distances in the wind field is studied. The feasibility of the integrated framework proposed here is verified and it can be used for future research projects of the automatic optimization and the selection of airfoils under prescribed constraints.