
MBW:Optimum Design of Blood Vessel BifurcationFrom MathBioContentsExecutive SummaryIn this article, we consider the optimal angle at which bifurcation in blood vessels occur. By minimizing a cost function, we are able to determine the optimal radius of a blood vessel, given its length and the flow of blood through it. We then use this optimal radius to determine our minimum cost as a function of radius and length. The lengths of the blood vessels, however, depend on where the bifurcation occurs. By finding the location of the optimal bifurcation which keeps the cost function minimized, we are able to write expressions for the bifurcation angles as a function of blood vessel radius. Furthermore, using conservation of mass, we are able to eliminate the need for knowing the radius of all three vessels. With this, we can determine the angle of bifurcation of a blood vessel without knowing its radius; instead, we just need to know the radii of the other two blood vessels. Once the expressions are derived, we can then use them to study special cases. The exemplary special case examined here occurs when the radii of the two daughter blood vessels are equal. In this case, the two blood vessels have the same bifurcation angle of . Biological Phenomenon and PurposeThroughout the human body, arteries from the heart bifurcate many times in order to become capillaries. The question proposed in this article is whether or not there exists some optimal angle at which a bifurcation occurs. We shall consider a bifurcation like that shown at right. The bifurcation occurs at the point B. To the left of B, some flow of blood is entering the vessel AB through A. To the right of B, is divided between the two vessels BC and BD, each now with flow and respectively. The goal now is to find a cost function which can be minimized in order to determine the bifurcation pattern of the blood vessels. Then, by varying the location of B, we will be able to determine the angle of bifurcation as a function of the radius of the blood vessels. Determination of Bifurcation AnglePoiseuille's LawBefore trying to determine the bifurcation angle, we must first understand the flow of blood through a given vessel. This is given by Poiseuille's Law, which states that . In this equation, is the viscosity of blood, is the radius and is the length of the vessel and and represent the pressures at one end of the vessel, respectively. By analyzing Poiseuille's Law, we see that even a small decrease in vessel radius requires a large increase in pressure if flow is to be kept constant. It is possible that this could lead to hypertension. (Fluid flow through a blood vessel is also described in MBW:Womersley Arterial Flow) Cost FunctionThe cost function proposed by Murray (1926) and Rosen (1967) is the sum of the rate at which work is done on the blood and the rate at which energy is used up by the blood vessels by metabolism [1]. The rate at which work is done on the blood is equal to , where is flow and is the change in pressure in the vessel. The rate at which energy is used by metabolism is equal to where K is a proportionality constant. The cost function for blood vessels is then . Using Poiseuille's Law, we can rewrite and obtain . To find the optimal radius for a given length and flow, the cost function is minimized with respect to
and is solved for the optimal :
At the optimal radius, the minimum cost function is equal to . Angle DeterminationBy moving the bifurcation point B, we now aim to find the angles at which the cost function is kept minimized. Our system is shown at right. In order for the cost function of the entire system to be optimal, the cost function of the individual blood vessels must also be optimal. With this requirement, we see that the total cost function (now denoted as ) is given by . In this system, the lengths of the blood vessels are dependent on the location of the bifurcation point. Suppose there is a small change in the location of the bifurcation. The change in the cost function is then . However, we want the cost function to still be minimized after this small change, i.e. we wish for to be equal to zero. This occurs at the optimal location of B. To determine the optimal location, we consider three different shifts in B. Shift Along ABSuppose that the bifurcation point is moved to location B' in the direction of the AB blood vessel, as shown at right. Then, the changes in the lengths of the vessels are
and the change in the cost function is . In order for the cost function to remain unchanged, the following relationship must be true:
Shift Along CBNow suppose that the bifurcation point is moved along the direction of CB, as shown at right. In this case, the changes in the lengths of the blood vessels are
The change in the cost function is . In order for the cost function to remain unchanged, the following relationship must be true:
Shift Along DBFinally, suppose the bifurcation point is moved along the direction of CD, as shown at right. In this case, the changes in the lengths of the blood vessels are
The change in the cost function is . In order for the cost function to remain unchanged, the following relationship must be true:
Expressions for AnglesWe now have three equations and three unknowns. The optimal conditions (1), (2), and (3) that resulted from each of the shifts can be solved for , , and . The result is
If we know the radius of each of the three vessels, we can now determine the angle of the bifurcation. Murray's LawBy the principle of the conservation of mass, the flow of blood in the vessel AB must be equal to the sum of the flow of blood through vessels BC and BD; namely,
If we solve our equation for optimal radius for flow, obtaining
and substitute this into our conservation equation, we obtain Murray's law.
Bifurcation Angle EquationsUsing Murray's Law, we can reduce our expressions for the angles of bifurcation so that only two of the three radii need to be known. The final result is the following
Special Case :In this section, we consider the case in which the daughter blood vessels have the same radius. In this case,
and our expressions for and reduce to
which are clearly equal. This means that, when the radii of the daughter blood vessels are equal, their bifurcation angles are also equal. Furthermore, we can determine exactly what the bifurcation angle is in this case. When the , Murray's law becomes
which can be rearranged as
Substituting this into reduced expression for results in
When we solve for , we find that the bifurcation angle is equal to 37.5 degrees.
Project Categorization(a) Mathematics Used: This project uses minimization and pertubation analysis to determine the optimum angles for blood vessel bifurcation. (b) Type of Model: Blood flow through a single vessel is modeled by Poiseuille's Law. This combined by a cost function proposed by Murray and Rosen yields a minimumcost function which is dependent on the length and radius of the vessel. (c) Biological System Studied This project studies what is the optimum angle for bifurcation of a blood vessel. Doing so helps us to understand Further ReadingBlood Vessel Branching: Beyond the Standard Calculus Problem [1] Blood flow dynamics in microvessel bifurcations [2] Related WorkThe work by Murray is also cited in the 2011 paper "Pulse Wave Propagation in the Arterial Tree" by Vosse and Stergiopulos [3]. In this paper the effect of a pulse through a circulatory system was studied. This is important because our blood does not flow at a constant rate but it pumped from the heart in pulses. Thus by being able to study how a pulse will propagate we can form a more accurate model of our circulatory system. The model uses Murray's law to help calculate the resistance to the blood flow in an arterial network. External Links1. Fung, Y.C., Biomechanics: Circulation, Section 3.3. (1997). 2. Mazumadar, J., An Introduction to Mathematical Physiology and Biology, Chapter 6. (1989). 3. Murray, Cecil D. (1926). "The Physiological Principle of Minimum Work: I. The Vascular System and the Cost of Blood Volume". Proceedings of the National Academy of Sciences of the United States of America 12 (3): 207–214. 4. Murray, Cecil D. (1926). "The Physiological Principle of Minimum Work: II. Oxygen Exchange in Capillaries". Proceedings of the National Academy of Sciences of the United States of America 12 (5): 299–304. 