High-order curvilinear hybrid mesh generation for CFD simulations