A quadrilateral mesh generation scheme using recursive domain bi divisions is enhanced in this paper. Robust mesh templates are designed

to provide a unified way for sub domain meshing. The templates pose no restrictions on the numbers of boundary nodes (e.g. six or eight) of sub domains

and can mesh the ‘big’ regular sub domains bounded linearly in many cases

thus to highly improve the timing and storage performance of the algorithm. Coupling the new sub domain mesh generation procedure with an automatic domain decomposition procedure

we obtain a new automatic unstructured quadrilateral mesh generation algorithm. Finally

mesh examples and simulation results are presented to demonstrate the capabilities of the new algorithm.