Interaction of blood flow and leaflet behavior in a bileaflet mechanical heart valve was investigated using computational analysis. Blood flows of a Newtonian fluid and a non-Newtonian fluid with Carreau model were modeled as pulsatile, laminar, and incompressible. A finite volume computational fluid dynamics code and a finite element structure dynamics code were used concurrently to solve the flow and structure equations, respectively, where the two equations were strongly coupled. Physiologic ventricular and aortic pressure waveforms were used as flow boundary conditions. Flow fields, leaflet behaviors, and shear stresses with time were obtained for Newtonian and non-Newtonian fluid cases. At the fully opened phase three jets through the leaflets were found and large vortices were present in the sinus area. At the very final stage of the closing phase, the angular velocity of the leaflet was enormously large. Large shear stress was found on leaflet tips and in the orifice region between two leaflets at the final stage of closing phase. This method using fluid-structure interaction turned out to be a useful tool to analyze the different designs of existing and future bileaflet valves.