A series of Reynolds Average Navier-Stokes simulations are performed for a supersonic, wall-bounded, turbulent corner flow. These are compared to a high order, implicit, LargeEddy Simulation for the same geometry and flow conditions. Since the LES results were obtained for a low Reynolds number, near the boundary of the validity for the turbulence models, the comparison of results with the two approaches was qualitative. Inclusion of the Quadratic Constitutive Relation in the Reynolds Average Navier-Stokes simulation results in significant improvement in qualitative agreement with the Large Eddy Simulation, specifically the presence of secondary flow (a counter rotating vortex pair). The range of valid values for the constant in the Quadratic Constitutive Relation formulation is explored. Additionally, the effects on this range from different turbulence models and momentum thickness Reynolds number for the corner are also examined. These effects of the Quadratic Constitutive Relation are explored using both the finite-difference fluid solver OVERFLOW, and the the finite-volume fluid solver US3D. The results indicate that the Quadratic Constitutive Relation term directly affects the strength of the vortex pair in the secondary flow and that its influence appears directly dependent on all the aforementioned parameters.