In this paper, a numerical simulation of steady laminar and turbulent flow in a two-dimensional model for the total artificial heart is presented. A trileaflet polyurethane valve was simulated at the outflow orifice while the inflow orifice had a trileaflet or a flap valve. The finite analytic numerical method was employed to obtain solutions to the governing equations in the Cartesian coordinates. The closure for turbulence model was achieved by employing the k-epsilon-E model. The SIMPLER algorithm was used to solve the problem in primitive variables. The numerical solutions of the simulated model show that regions of relative stasis and trapped vortices were smaller within the ventricular chamber with the flap valve at the inflow orifice than that with the trileaflet valve. The predicted Reynolds stresses distal to the inflow valve within the ventricular chamber were also found to be smaller with the flap valve than with the trileaflet valve. These results also suggest a correlation between high turbulent stresses and the presence of thrombus in the vicinity of the valves in the total artificial hearts. The computed velocity vectors and turbulent stresses were comparable with previously reported in vitro measurements in artificial heart chambers. Analysis of the numerical solutions suggests that geometries similar to the flap valve (or a tilting disk valve) results in a better flow dynamics within the total artificial heart chamber compared to a trileaflet valve.