Takamasa Onuki

Computer simulation of geometry and hemodynamics of canine pulmonary arteries

Assuming that along the pulmonary arteries, mathematical expectations of branching are equal to matrix (Ei,j; i: parent branch class, j: daughter branch class number), we made computer drafts of patterns of pulmonary arterial branching trees from the published data of canine pulmonary arterial casts with the “Monte Carlo methods” by FORTRAN, calculating from these data, the numbers of branches, size and resistance of the trees. We also analyzed blood flow distribution in a pulmonary branch according to the pathlength from the entrance of the tree. The graphics of the trees were similar to the original arterial casts, and numbers of the branches were nearly the same as the previous reported values. Calculated resistance was reasonable. The rate of blood flow per unit vessel in the short pathlength group was calculated as 1.5 times that of the long pathlength group of the same diameter. We believe that our method of describing the branching pattern matrix (Ei,j) provides a reasonable simulation of complex branching patterns such as pulmonary arteries and a useful means to analyze local hemodynamics.