We examine two dimensional properties of vortex shedding past elliptical cylinders through numerical simulations. Specifically, we investigate the vortex formation length in the Reynolds number regime 10 to 100 for elliptical bodies of aspect ratio in the range 0.4 to 1.4. Our computations reveal that in the steady flow regime, the change in the vortex length follows a linear profile with respect to the Reynolds number, while in the unsteady regime, the time averaged vortex length decreases in an exponential manner with increasing Reynolds number. The transition in profile is used to identify the critical Reynolds number which marks the bifurcation of the Karman vortex from steady symmetric to the unsteady, asymmetric configuration. Additionally, relationships between the vortex length and aspect ratio are also explored. The work presented here is an example of a module that can be used in a project based learning course on computational fluid dynamics.