In this dissertation, we propose conforming finite element methods that yield divergence-free velocity approximations for the steady Stokes problem on cubical and quadrilateral meshes. In the first part, we construct the finite element spaces for the two-dimensional problem on rectangular grids. Then in the second part, we extend these spaces to n-dimensional spaces. We use discrete differential forms and smooth de Rham complexes to verify the stability and the conformity of the proposed methods, and the solenoidality of the velocity approximations. In the third part, we shift our focus from a dimensionwise extension to a meshwise improvement by introducing macro elements on general shape-regular quadrilateral meshes. By utilizing a smooth de Rham complex, we prove that the macro finite element method yields divergence-free velocity solutions, and with the construction of a Fortin operator, we validate the stability of the method. To improve the pressure approximation properties, we compute a post-processed pressure solution locally. In addition, we describe the implementation process of the (velocity) macro finite element. We show that the methods developed in this dissertation yield optimal convergence rates and present numerical experiments which are supportive of the theoretical results. Moreover, we provide experimental results of our method for the Navier-Stokes equations and show that the convergence rates are preserved
