A method for automatically generating hex-dominant meshes for CFD applications is presented in this paper. Two important regions of the mesh for any CFD simulation are the interior mesh and the boundary layer mesh. The interior mesh needs to be fine in the critical flow regions to ensure accurate solutions. The proposed method uses directionality and sizing control functions while generating the mesh that offers a flexible control over mesh sizing and local mesh refinement. Boundary layer mesh is important to accurately model the physics of boundary layer near the geometry walls. In the proposed method, boundary elements in the mesh are split into multiple divisions with first division having smallest thickness to ensure it lies inside the physical boundary layer. Rest of the layers expand linearly in thickness until they match in value with the interior mesh size. Hex-dominant meshes were created for two different flow geometries. Flow simulations were performed in open source code OpenFOAM and predictions were compared with available experimental data.