Computation of flow fields and shear rates in an aortic bifurcation.

A finite volume method in a boundary-fitted coordinate system together with a zonal grid method is employed to compute the flow fields and shear stresses in a two-dimensional aortic bifurcation. Eddy is found distal to plaques during pulsatile flow, whereas permanent eddies are observed only during steady flow. The computed flow fields are consistent with those visualized experimentally by other authors. It is also found that although the time averaged shear rates in a pulsatile flow are similar to those of a steady flow with mean Reynolds number in most regions, they are different in recirculation zones. This result implies that care should be taken if a steady flow shear rate were to be used in modeling shear-dependent physiological processes. The non-Newtonian viscosity has only a minor effect on the flows.