Skip to main content
  • Other Publications
    • Philosophical Transactions B
    • Proceedings B
    • Biology Letters
    • Open Biology
    • Philosophical Transactions A
    • Proceedings A
    • Royal Society Open Science
    • Interface
    • Interface Focus
    • Notes and Records
    • Biographical Memoirs

Advanced

  • Home
  • Content
    • Latest issue
    • All content
    • Subject collections
    • Videos
  • Information for
    • Authors
    • Reviewers
    • Readers
    • Institutions
  • About us
    • About the journal
    • Editorial board
    • Author benefits
    • Policies
    • Citation metrics
    • Publication times
    • Open access
  • Sign up
    • Subscribe
    • eTOC alerts
    • Keyword alerts
    • RSS feeds
    • Newsletters
    • Request a free trial
  • Submit
You have accessRestricted access

Optimal concentrations in transport systems

Kaare H. Jensen, Wonjung Kim, N. Michele Holbrook, John W. M. Bush
Published 17 April 2013.DOI: 10.1098/rsif.2013.0138
Kaare H. Jensen
Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wonjung Kim
Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USADepartment of Mechanical Engineering, Sogang University, 1 Shinsu-dong, Seoul 121-742, Republic of Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
N. Michele Holbrook
Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John W. M. Bush
Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Many biological and man-made systems rely on transport systems for the distribution of material, for example matter and energy. Material transfer in these systems is determined by the flow rate and the concentration of material. While the most concentrated solutions offer the greatest potential in terms of material transfer, impedance typically increases with concentration, thus making them the most difficult to transport. We develop a general framework for describing systems for which impedance increases with concentration, and consider material flow in four different natural systems: blood flow in vertebrates, sugar transport in vascular plants and two modes of nectar drinking in birds and insects. The model provides a simple method for determining the optimum concentration copt in these systems. The model further suggests that the impedance at the optimum concentration μopt may be expressed in terms of the impedance of the pure (c = 0) carrier medium μ0 as μopt∼2αμ0, where the power α is prescribed by the specific flow constraints, for example constant pressure for blood flow (α = 1) or constant work rate for certain nectar-drinking insects (α = 6). Comparing the model predictions with experimental data from more than 100 animal and plant species, we find that the simple model rationalizes the observed concentrations and impedances. The model provides a universal framework for studying flows impeded by concentration, and yields insight into optimization in engineered systems, such as traffic flow.

1. Introduction

Transport systems are ubiquitous in nature and technology. Whether biological such as the vascular systems of plants and animals or engineered such as man-made pipes, roads, electrical grids and the Internet, they serve to move matter, energy or information from one place to another. Owing to the cost of constructing and maintaining redundant channels, it is advantageous for biological transport systems to distribute matter efficiently [1,2]. Oxygen transport in vertebrates [3,4], sugar transport in plants [5] and drinking strategies of many animals [6,7] are known to be optimized for efficient transport of energy and material. Engineered systems must likewise be cost-effective and able to provide efficient transport under a variety of conditions; for example, considerable resources are spent annually to ease traffic congestion.

In our examination of transport systems, we consider material flow in four different natural systems: blood flow in vertebrates, sugar transport in vascular plants, and two modes of nectar drinking in birds and insects. A common feature of these and other transport systems is that the flow impedance depends on concentration. While the most concentrated solutions offer the greatest potential in terms of material transfer, the increase of impedance with concentration also makes them the most difficult to transport. Additionally, most transport systems are subject to a set of limiting constraints. For example, nectar feeders are typically constrained by a constant work rate, which in turn is a function of flow impedance and hence concentration [7–9]. These transport systems may thus be characterized in terms of an optimization problem subject to appropriate constraints. This approach has been used to rationalize observed concentrations in a wide range of natural systems [4,8–10]. For example, the observed volume fraction of erythrocytes (red blood cells), typically approximately 40–50% in humans, has been shown to maximize oxygen transport [4,10]. With the widespread use of bioinspired design in the development of novel engineered systems, it seems likely that man-made transport systems such as roads or the electrical grid may benefit from improved understanding of natural transport systems.

We here develop a general framework for determining the concentration that maximizes material transfer in transport systems. By drawing on a number of natural examples—either new or drawn from the biology literature—we show how these can be treated within a single framework that provides new insight into the efficiency of transport systems. We compare our model predictions with experimental data from more than 100 animal and plant species collected from the literature. Finally, we show that similar optimization criteria may be applied to engineered systems, and consider traffic flow in the context of our new framework.

2. General formulation

We consider systems in which the material transfer rate (material flow) J can be expressed as the product of a volumetric flow rate (volume flow) Q and a concentration of material cEmbedded Image 2.1We express the volume flow as Q = Xf/μ, where X is a constant geometric factor, f quantifies the mechanism driving the flow and μ characterizes the impedance. The material flow in equation (2.1) can be expressed asEmbedded Image 2.2where f and μ can depend on c. We now seek the optimal concentration c = copt that maximizes J in equation (2.2) subject to a set of constraints; constant driving force f = f0 or constant work rate W ∼ Qf ∼ f2/μ. We thus consider constraints of the form f ∼ μγ, where γ = 0 corresponds to constant driving force and Embedded Image to constant work rate when f depends on μ. Other values of γ are possible, if there is a direct coupling between the driving force and impedance; for example, Embedded Image for bees that use viscous dipping to drink nectar (see §3.4).

Although we may have limited knowledge of the exact functional form of the concentration-dependent material flow J(c) in equation (2.2), two general statements can be made. First, we expect J to be proportional to the concentration c at low concentrations and to approach zero with the concentration, i.e. J(c) ∝ c, when c ≪ copt. Second, we are concerned with situations where the system impedance increases with concentration, and J increases monotonically up to a maximum value J(copt). To describe the system, we thus propose the first order governing equationEmbedded Image 2.3where J* = J(c)/J(copt) and c* = c/copt is the normalized material flow and concentration, respectively. A and B are parameters that are determined by the boundary conditions: J* (0) = 0, J* (1) = 1 and Embedded Image. This leads toEmbedded Image 2.4While equation (2.4) does not reveal the absolute value of the optimum concentration copt, it does contain information concerning the impedance at the optimum concentration μ (copt). By expressing the constraint as f ∝ μγ, we find from equation (2.2) that the normalized flux J* =J(c)/J(copt) can be written asEmbedded Image 2.5By using equation (2.4), the normalized impedance μ* = μ(c)/μ(copt) may be expressed asEmbedded Image 2.6It follows that the impedance at the optimum concentration isEmbedded Image 2.7where the power α = 1/(1 − γ) = log2(μ(copt)/μ(0)) and μ (0) is the impedance at zero concentration.

Equations (2.4), (2.6) and (2.7) provide a general framework for analysing optimization of concentration-impeded material flow in biological and engineered systems. To test the quantitative predictions of the theory, we proceed in §3 by considering a series of biological examples where the flux J can be optimized along the lines outlined above. In §4, we apply our model to traffic flow. Finally, in §5, we consider universal properties of concentration-impeded material transport systems.

3. Biological transport systems

3.1. Nectar drinking from a tube

Perhaps the simplest situation in which we may apply equation (2.2) is drinking from a cylindrical tube. Many insects and birds such as butterflies and hummingbirds [7] feed on floral nectar, an aqueous solution of sugars, through tubes formed from proboscises or tongues. Quick energy ingestion is advantageous for nectar feeders owing to the threat of predation. While the sweetest nectar offers the greatest energetic rewards, the increase of viscosity with sugar concentration also makes the sweetest nectar the most difficult to transport [11]. An optimal concentration may thus be sought for maximizing energy uptake rate.

Two different suction mechanisms are typically used by nectar feeders: active suction and capillary suction [7–9,11]. Active suction feeders such as butterflies use muscle contraction to suck nectar through their roughly cylindrical proboscises. In the limit of low Reynolds number Hagen–Poiseuille flow, the nectar mass flow rate Js can be expressed asEmbedded Image 3.1where a is the radius and l is the length of the proboscis, Embedded Image is the w/w sugar concentration, η is the viscosity (see appendix A), ρ is the density of the nectar solution and Δp is the pressure difference generated by muscular contraction. The manner in which biological constraints determined the dependence of the pressure Δp on nectar viscosity has been treated elsewhere [8,9]. Active suction feeders are typically constrained by constant work rate W = QΔp = πa4/(8ηl)Δp2, so the pressure Δp = (8Wl/(πa4))1/2η1/2 depends on viscosity and, hence, concentration. Comparing the nectar flow rate in equation (3.1) with the general expression in (2.2), we find that the impedance corresponds to the viscosity of the sugar solution μ = η, the concentration to Embedded Image, the geometric factor to X = πa4/(8l), and the driving mechanism to the pressure f = Δp. We can thus express the constraint as f = (8Wl/(πa4))1/2μ1/2 (table 1).

View this table:
  • View inline
  • View popup
Table 1.

Parameters describing the material flow J = Qc = Xfc/μ (see equation (2.2)) for each of the systems considered. See appendices A and B for details on the viscosity η of blood, nectar and phloem sap (BHN, Bando, Hasebe and Nakayama).

Capillary suction feeders such as hummingbirds use surface tension to draw nectar along their tongue, during repeated cycles of tongue insertion and retraction [7,9]. If the duration of a cyclic motion is the sum of the nectar loading time T and unloading time T0, the average nectar mass flow rate Js can be expressed as Embedded Image where πa2l(T) is the time-dependent nectar volume extracted during the loading. In the loading phase, the volumetric flow rate is given by πa2(dl/dt) = πa4Δp/(8ηl), where Δp = 2σ/a is the capillary pressure. The solution with initial condition l(0) = 0 is given by l(t) = (aσt/(2η))1/2, which depends explicitly on viscosity and hence concentration. The average nectar mass flow rate Js can be expressed asEmbedded Image 3.2By comparing the nectar flow rate in equation (3.1) with the general expression in equation (2.2), we find that X = πa3, f = (σηT/(2a(T + T0)2))1/2 and μ = η. If T/(T+T0)2 is assumed to be independent of viscosity [12], we find a relation between driving force and impedance: f ∝ μ1/2 (table 1).

For both active and capillary suction, we find that f ∝ μ1/2 (i.e. Embedded Image) and the optimal concentration Embedded Image can thus be found by maximizing c/μ1/2. For nectar sugar solutions, we thus predict that Embedded Image w/w and Embedded Image (i.e. α = 2 in (2.7)). This is in good agreement with experimental data on 16 butterfly and hummingbird species (figure 1a and table 2), where optimal concentrations in the range 30–45% are reported.

View this table:
  • View inline
  • View popup
Table 2.

Comparison between theoretical predictions (T) and experimental observations (E) of the optimum concentration copt, the optimum viscosity μopt and the exponent α. Concentration units are % w/w for nectar drinking and sugar transport in plants, % v/v for blood flow and % vehicle density/max. vehicle density for traffic flow. The experimental data are available in the electronic supplementary material.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Optimal concentrations in biological transport systems. (a) Drinking from a tube. Histogram showing the distribution of observed sugar concentrations that maximize nectar uptake for 16 bird and insect species that use muscular contractions or surface tension to feed through cylindrical tubes [9,13]. Normalized sugar mass flow Js/Js,max (solid line, equations (3.1) and (3.2)) and nectar viscosity μ/μ0 (dashed line, data from [14]) are plotted as a function of nectar sugar concentration Embedded Image. Mass flow is predicted to be maximum when c̄opt = 35%, in good agreement with the observed average nectar concentration (37%). (b) Blood flow. Histogram showing the distribution of observed red blood cell concentrations (haematocrit) from 57 vertebrate species [4]. Normalized oxygen flow Jr/Jr,max (solid line, equation (3.3)) and blood viscosity μ/μ0 (dashed line, see appendix B) are plotted as a function of haematocrit Embedded Image. Flow is predicted to be maximum when c̃opt = 39%, in good agreement with the observed average haematocrit (40%). (c) Sugar transport in plants. Histogram showing the distribution of observed sugar concentrations from 28 plant species that use active sugar loading [15]. Normalized sugar flow Jp/Jp,max (solid line, equation (3.4)) and sap viscosity μ/μ0 (dashed line, data from [14]) are plotted as a function of phloem sugar concentration Embedded Image. Mass flow is predicted to be at a maximum when c̄opt = 24%, in good agreement with the observed average sugar concentration (22%). (d) Nectar drinking by viscous dipping. Histogram showing the distribution of observed sugar concentrations that maximize nectar uptake for six insect species that use viscous dipping [9,13]. Normalized sugar mass flow Jv/Jv,max (solid line) and nectar viscosity μ/μ0 (dashed line, data from [14]) are plotted as a function of nectar sugar concentration Embedded Image. Mass flow is predicted to be at a maximum when c̄opt = 57%, in good agreement with the observed average nectar concentration (55%). In (a)–(d), the numbers given above the bins (coloured bars) indicate the percentage of species in the bin. The experimental data are available in the electronic supplementary material.

3.2. Blood flow in vertebrates

Another biological flow problem that can be analysed within our framework is oxygen and nutrient transport within the cardiovascular system of vertebrate animals. Here, red blood cells transport oxygen between the lungs and distal parts of the organism. The cells are suspended in blood plasma, which primarily consists of water [16]. Red blood cells typically measure 10 µm in diameter [17], and the blood's bulk viscosity increases with the haematocrit Embedded Image, the volume concentration of red blood cells (see appendix B). While blood with the highest haematocrit is the most oxygen rich, the increase of viscosity with the haematocrit also makes such blood the most difficult to transport. Accordingly, an optimal haematocrit may be sought for maximizing oxygen transport.

In the limit of low Reynolds number Hagen–Poiseuille flow, the red blood cell volume flow rate Jr in vessels larger than 1 mm can be expressed asEmbedded Image 3.3where a and l are the radius and length of the blood vessel, η is the blood viscosity (see appendix B) and Δp is the pressure difference generated by the heart. The diastolic pressure is of the order of 10 kPa for most animals [18], but the pressure difference associated with flow in large vessels is small; most of the pressure drop in blood flow occurs in vessels with size comparable to red blood cells. The dependence of blood pressure on the haematocrit is, however, negligible [4], i.e. Δp does not depend on Embedded Image. Although blood viscosity η generally depends on the shear rate, this dependence is weak for typical blood conditions, specifically in vessels with diameters larger than 1 mm and shear rates greater than 50 s−1 [16]. Comparing the blood flow rate in equation (3.3) with the general expression in equation (2.2), we find that Embedded Image, μ = η, X = πa4/(8l) and f = Δp. We express the constraint as that of constant pressure f = f0 = Δp (table 1).

For blood flow, we thus find that Xf ∝ μ0 (i.e. γ = 0) and the optimum concentration Embedded Image can be found by maximizing c/μ. We thus predict that Embedded Image v/v and Embedded Image (i.e. α = 1, in (2.7)), in good agreement with experimental data from 57 species observed throughout the animal kingdom (figure 1b and table 2). We attribute the significant variation in the observed concentrations in part to the complex interactions between red blood cells and flow in smaller vessels, which we do not consider in our model. An important feature of the red blood cell volume flow rate Jr is also that it varies by less than 20 per cent over the range of concentrations from c = 20–60%, suggesting that concentrations in this interval are acceptable given additional biological constraints. For example, we note that, for some diving mammals (e.g. Weddell seals and whales), oxygen storage in the blood may also be an important factor, resulting in a higher haematocrit (up to 63%; [4]). It is also likely that lack of thermoregulation may explain why poikilothermic animals (e.g. the rainbow trout) have a lower haematocrit value (23%) than the average, probably as a result of thermally induced variations in blood viscosity [19].

3.3. Sugar transport in plants

Plants, like animals, rely on vascular systems for distribution of energy and nutrients. Energy distribution in plants takes place in the phloem vascular system. Here, an aqueous solution of sugars, amino acids, proteins, ions and signalling molecules flows through a series of narrow elongated cylindrical cells, known as sieve tube elements, that lie end-to-end, forming a microfluidic distribution system spanning the entire length of the plant. The flow is driven by differences in chemical potential between distal parts of the plant [20]. While phloem sap with high sugar concentration has the greatest potential for energy transfer, the increase of viscosity with sugar concentration makes it the most difficult to transport. Accordingly, an optimal concentration may again be sought for maximizing energy flow.

Assuming low Reynolds number Hagen–Poiseuille flow, the phloem sugar mass flow rate Jp can be expressed asEmbedded Image 3.4where a is the radius of the phloem sieve tube (a ≃ 10 µm), l is the length of the plant, Embedded Image is the sugar concentration, η is the phloem sap viscosity and Δp is the pressure difference driving the flow. By comparing the sugar flow rate (3.4) with the general expression in (2.2), we find that Embedded Image, Embedded Image, Embedded Image and f = Δp. We express the constraint as that of constant pressure f = f0 = Δ p (table 1).

For sugar transport in plants, we thus find that Xf ∝ μ0 (i.e. γ = 0), and the optimum concentration Embedded Image can thus be found by maximizing c/μ. We find that Embedded Image w/w and Embedded Image (i.e. α = 1, in (2.7)), in good agreement with experimental data (figure 1c and table 2). While sugar concentrations observed in plants generally span a wide range, this analysis provides a rationale for the observation that plants that use active sugar loading (data shown in figure 1c) typically have a higher sugar concentration than plants that use passive loading [15]. Active loaders expend metabolic energy to increase the sugar concentration in the phloem [21]. The process is driven by membrane transporters and sugar polymerization and occurs against a sugar concentration gradient. However, in passive loading species, sugars move into the phloem without the use of metabolic energy by travelling down a concentration gradient from sites of carbohydrate synthesis and/or storage to the phloem [15]. We also note that plants with the highest sugar concentrations are crop plants, for example potato (50%) and maize (40%), suggesting that selection for high crop yield tends to lead to increased sugar concentration in the phloem sap [15].

3.4. Drinking by viscous dipping

So far, we have limited our attention to transport in closed channels. However, it is straightforward to extend the problem to situations where free surfaces are involved. Most bees whose tongues are solid rather than hollow use a drinking style termed ‘viscous dipping’ in which the fluid is entrained by the tongue surface. The average nectar volume entrained can be expressed by Q ∼ 2πaeu, where a is the tongue radius, e is the thickness of the nectar layer on the tongue and u is the tongue extraction speed. Based on Landau–Levich–Derjaguin theory when the Reynolds number Re ≪ 1 and Bond number Bo ≪ 1, the nectar film thickness is given by e ∼ aCa2/3, where Ca = ηu/σ ≪ 1 is the ratio of viscous to capillary forces [22]. Because the fluid is entrained on the tongue by viscous forces, we define the driving force and geometric factor as f = ηue/a2 and X = 2πa3. The movement of the tongue in the fluid requires power W ∼ ηu2l to overcome the viscous drag, where l is the immersed tongue length. Assuming a constant work rate W for a given creature leads to the constraint on velocity u ∼ (W/(ηl))1/2, which in turn leads to f ∝ μ 5/6 ([9]; table 1).

For viscous dipping, we find that f ∝ μ 5/6 (i.e. Embedded Image) and the optimum concentration Embedded Image can thus be found by maximizing c/μ1/6. We find that copt = 57% w/w and μ (copt) = 64 μ0 (i.e. α = 6, cf. equation (2.7)). This is in reasonable agreement with experimental data on six bees species (figure 1d and table 2), where optimal concentrations in the range 50–60 % are found. This may explain why the nectar concentration of flowers pollinated by bees is generally higher than that of those pollinated by tube-feeding butterflies and hummingbirds [9].

4. Applications to engineered transport systems: traffic flow

We have thus far seen many qualitative similarities between different biological flows. Although the detailed physiological and physical mechanisms are different, provided increased concentration leads to greater impedance, we can rationalize the optimal concentrations. An interesting question naturally arises. In which engineered systems might one expect to observe similar phenomena? It appears likely that most efficient communication and transport systems will exhibit similar features. Nevertheless, we limit our discussion to traffic congestion on highways.

A measure of the efficiency of a given section of road is the vehicle flow Jv, the number of vehicles passing a given point per unit time [23–26]. Designers of road networks strive to maximize the vehicle flow that can be expressed as Jv = ρv, where v is the speed of the individual vehicle and ρ is the number of vehicles per unit length of roadway. Generally, car speed v = v(ρ) is a decreasing function of density ρ. At very low densities, where inter-vehicle interaction is negligible, however, the speed approaches the speed limit vmax and the vehicle flux is proportional to density Jv ≃ ρvmax. At higher vehicle densities, interaction between adjacent cars leads to flow impedance and a significant reduction in the speed of individual vehicles, causing congestion and a net decrease in the flux Jv. The vehicle interactions initially take the form of synchronized flow, a form of congested traffic in which each driver attempts to maintain a safe distance from the neighbouring cars. As the density increases, wide moving jams form, that is, stop-and-go traffic in which the vehicle flux approaches zero [23]. From these considerations, one anticipates an optimal vehicle density ρopt that maximizes the vehicle flux.

To estimate ρopt, we require v(ρ), which can be either found empirically or deduced from vehicle interaction models. One of the simplest models that leads to a reasonable expression for v(ρ) was proposed by Greenberg [27], who treated traffic flow as a one-dimensional flow of an ideal compressible gas. He assumed (i) that the local speed is a function of density only v = v(ρ(x,t)), (ii) that vehicles are conserved Embedded Image, (iii) that vehicle flow satisfies the Euler equation Embedded Image and (iv) that traffic ‘pressure’ is proportional to density p = 𝒞2ρ. This leads to the relation v(ρ) = 𝒞 ln(ρmax/ρ), where ρmax is the density at which traffic stops owing to congestion. The vehicle flow rate Jv = 𝒞ρ ln(ρmax/ρ) is at a maximum when ρ = ρopt = ρmax/e, and the constant 𝒞 = v(ρopt) is the vehicle speed at the optimal concentration. Because vehicles typically occupy 7.5 m in a totally congested flow [23], we estimate that ρmax ≃ 133 vehicles per kilometre. Greenberg's model overestimates the optimal density, predicting ρopt = ρmax/e ∼ 50 vehicles per kilometre, whereas the true value is known to be ∼ 20 vehicles per kilometre. Nevertheless, the vehicle flow rate Jv is qualitatively consistent with empirical traffic data (figure 2). The data are plotted as a function of vehicle concentration c = ρ/ρmax in figure 2 along with Greenberg's flow rate Jv, deduced using ρmax = 133 vehicles per kilometre.

A shortcoming of Greenberg's theoretical model is that the vehicle speed v diverges when the car density is very low. To ensure that v(ρ/ρmax → 0) = vmax and to account for other aspects of traffic flows, numerous other models have been proposed [23–26,28]. For example, Bando, Hasebe and Nakayama (BHN) [28] suggested a traffic model in which the vehicle speed depends on the distance from the car in front, Δx. This leads to v = vmaxtanh(Δx/s), where s is a fixed length scale determined by the road conditions. With a minimum vehicle distance L = 7.5 m, we can express the density in terms of Δx as ρ = 1/(L+Δx). This leads to v(ρ) = vmaxtanh[(1/ρ − L)/s], in which case the flow rate Jv = vρ is optimized when ρ = 0.21 and ρmax = 28 vehicles per kilometre. With vmax = 120 km h−1 and s = 60 m, the BHN model provides a better quantitative fit to the empirical data than does Greenberg's model (figure 2).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Optimal vehicle concentration for maximizing traffic flow. Grey dots show measured vehicle flow rate Jv plotted as a function of vehicle concentration c = ρ/ρopt, where ρopt = 133 vehicles per kilometre. The flow rate is normalized by 1483 vehicles per hour, which corresponds to Jv(ρopt) = Jv,max in Bando, Hasebe and Nakayama's (BHN) model [28]. Histograms show the states occupied by the system in the morning (green, 06.00–08.00) and evening (blue, 16.00–18.00) rush-hour traffic. The data were collected by the Minnesota Department of Transportation from a sensor on the westbound direction of I-94 (Minneapolis, MN, USA) on Fridays (7, 14, 21, 28) in September 2012 [29]. The predicted vehicle transport rate Jv/Jv,max (thick solid black line, BHN model; thin solid red line, Greenberg's model) and traffic impedance μ/μ0 (dashed line, BHN model) are plotted as a function of vehicle concentration c. The experimental data are available in the electronic supplementary material.

Comparing the Greenberg and BHN models of traffic flow with the formulation introduced in equation (2.2), we see that traffic flow can be treated in the same general framework whereEmbedded Image 4.1for Greenberg's model, andEmbedded Image 4.2for the BHN model. In both cases, N is the number of lanes.

Comparing traffic flow with the biological transport problems considered above, we find that the normalized flux and impedance curves follow the same pattern (figure 2). While traffic flow can be treated in the same framework as biological flows, it is important to note that the congested highway (figure 2, data recorded from 16.00 to 18.00) is very far from being optimized. This is presumably the result of two main effects. First, the individual vehicle operator attempts to minimize his or her own travel time, which does not necessarily optimize the overall vehicle flow Jv. Second, traffic flows are intrinsically time dependent, which leads to the formation of travelling density waves and shocks [23–26].

5. Universal properties of transport systems

To compare characteristics of the particular biological and man-made transport systems considered in §3 and §4 with the general formulation in equations (2.4) and (2.6), normalized material flow and impedance curves are plotted in figure 3. Despite the complex dependence of impedance on concentration (see appendices A and B), both the material flow J* and impedance μ* are adequately approximated by the simple forms given in equations (2.4) and (2.6). From (2.6), it follows that the impedance at the optimum concentration is μopt = 2αμ0, where μ0 is the impedance of the pure carrier medium (with c = 0) and the power α is determined by the flow constraints. In the cases of vascular transport in plants and animals, the power α = 1, because there is no coupling between the constant driving pressure (f) or the vascular geometry (X) and the impedance (μ). This suggests that the optimum in material flow should occur when the blood or phloem sap is twice as viscous as water, i.e. ηopt = 2η0, in good agreement with observed values (table 2).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Universal properties of biological and engineered flows. (a) Normalized flow rate J* = J(c)/J(copt) plotted as a function of normalized concentration c* = c/copt. The solid thick black line shows the prediction of equation (2.4). (b) Normalized impedance μ* = μ(c)/μ(copt) plotted as a function of normalized concentration c*. The solid and dashed thick black lines show the predictions of equation (2.6). The inset indicates the dependence of (μ*)1−γ on c*.

In transport systems that are constrained, for example by constant work rate, α will generally be greater than unity, because of the coupling between flow and impedance. The impedance at the optimum concentration μopt = 2αμ0 can thus be significantly greater than that of the carrier medium. This is most clearly seen in the case of viscous dipping (§3.4), where the observed nectar viscosity is up to 50 times greater than that of water, roughly consistent with the value (26 = 64) predicted by our simple model (table 2).

These observations suggest that this general framework may also provide the rationale for the viscosities found in other biological transport systems where efficient transport is favoured. Examples of systems with constant forcing include mammals that drink whole milk (observed viscosity: η ∼ 2η0; [30]), and the macro-alga Chara where streaming distributes the content of the cell cytosol (observed viscosity: ∼3η0 [31]). Although detailed studies of these systems are left for future consideration, we note that both are roughly consistent with the predictions of our general theory with α = 1.

Comparing traffic flow with the biological transport problems considered, we find that the normalized flux and impedance curves follow the same pattern (figure 3 and table 2). Since the speed limit vmax, which is fixed on a given road section, corresponds to the flow-driving mechanism in the BHN model, traffic flow is analogous to vascular transport in animals and plants that operate at constant pressure. Our model thus indicates that the flow constraint does not couple to impedance, Xf ∝ μ0 (γ = 0, α = 1), and hence that the optimal impedance is μopt = 2 μ0. This is in rough accord with the BHN model that yields μopt = 1.9 μ0.

6. Discussion and conclusion

We have seen many qualitative and quantitative similarities between different natural and engineered transport systems. Although the detailed transport mechanisms are different, key common features have allowed us to develop a general framework. Provided impedance increases with concentration, our model provides a means of rationalizing the optimal concentrations. Collecting data from more than 100 plant and animal species, we have observed that optimization of material flow appears to be a universal feature of biological transport systems. This deduction provides a rationale for the observation that the simple model introduced in §2 collapses flow and impedance curves for all the systems considered (figure 3), suggesting a universal component to all natural transport systems.

Finally, we have shown that an interesting analogy can be made between biological systems and self-driven systems such as traffic flows. Here, we find that the impedance analogy is still valid, but that the system is far from optimized owing to conflicting interests between individuals and the collective. The consideration of other man-made transport systems, such as the electrical grid or the Internet, is left for future consideration.

Acknowledgements

The authors wish thank Maciej Zwieniecki, Ruben Rosales, Jessica Savage, Nick Carroll, Kenneth Ho and David Weitz. This work was supported by the NSF (grant nos 1021779 and DMS-0907955) and the Materials Research Science and Engineering Center (MRSEC; grant no. DMR-0820484) at Harvard University.

Appendix A. Viscosity and density of nectar and phloem sap

Phloem sap and flower nectar consist of an aqueous solution of sugars, amino acids, proteins and other nutrients. Sugars, of which sucrose, fructose and glucose are the most abundant types, constitute about 90 per cent of the total solute mass [32]. To approximate the viscosity η and density ρ of phloem sap and nectar, we thus used data from sucrose solutions of concentration Embedded Image obtained from Hainsworth [12]. Least squares fits to sucrose data yield the approximate expressions for viscosity Embedded Image and Embedded Image and density Embedded Image. We note that viscosity and density data from other sugar types (glucose and fructose) are well approximated by the fit, suggesting that the major determinant of viscosity is the mass fraction Embedded Image, and not the type of sugar.

Appendix B. Viscosity of blood

Vertebrate blood is composed of blood cells suspended in blood plasma, a liquid that consists mostly of water. The viscosity of blood η depends primarily on the volume concentration Embedded Image (haematocrit) of red blood cells, and on temperature [4,19]. As demonstrated by Saitô [33] and Stark & Schuster [4], blood viscosity is well described by the function Embedded Image, which for blood vessels with diameters larger than 1 mm is consistent with empirical data with less than 5 per cent error for 0 < Embedded Image ([34,35]; electronic supplementary material, table S6).

  • Received February 12, 2013.
  • Accepted March 25, 2013.
  • © 2013 The Author(s) Published by the Royal Society. All rights reserved.

References

  1. ↵
    1. LaBarbera M
    . 1990 Principles of design of fluid transport systems in zoology. Science 249, 992–1000. doi:10.1126/science.2396104 (doi:10.1126/science.2396104)
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Vogel S
    . 2004 Living in a physical world. J Biosci. 29, 391–397. doi:10.1007/BF02712110 (doi:10.1007/BF02712110)
    OpenUrlCrossRefPubMedWeb of Science
  3. ↵
    1. Murray CD
    . 1926 The physiological principle of minimum work. I. The vascular system and the cost of blood volume. Proc. Natl Acad. Sci. USA 12, 207–214. doi:10.1073/pnas.12.3.207 (doi:10.1073/pnas.12.3.207)
    OpenUrlFREE Full Text
  4. ↵
    1. Stark H,
    2. Schuster S
    . 2012 Comparison of various approaches to calculating the optimal hematocrit in vertebrates. J. Appl. Physiol. 113, 355–367. doi:10.1152/japplphysiol.00369.2012 (doi:10.1152/japplphysiol.00369.2012)
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Jensen KH,
    2. Lee J,
    3. Bohr T,
    4. Bruus H,
    5. Holbrook NM,
    6. Zwieniecki MA
    . 2011 Optimality of the Münch mechanism for translocation of sugars in plants. J. R. Soc. Interface 8, 1155–1165. doi:10.1098/rsif.2010.0578 (doi:10.1098/rsif.2010.0578)
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Kim W,
    2. Bush JWM
    . 2012 Natural drinking strategies. J. Fluid Mech. 705, 7–25. doi:10.1017/jfm.2012.122 (doi:10.1017/jfm.2012.122)
    OpenUrlCrossRefWeb of Science
  7. ↵
    1. Kim W,
    2. Peaudecerf F,
    3. Baldwin MW,
    4. Bush JWM
    . 2012 The hummingbird's tongue: a self-assembling capillary syphon. Proc. R. Soc. B 279, 4990–4996. doi:10.1098/rspb.2012.1837 (doi:10.1098/rspb.2012.1837)
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Pivnick K,
    2. McNeil J
    . 1985 Effects of nectar concentration on butterfly feeding: measured feeding rates for Thymelicus lineola (Lepidoptera: Hesperiidae) and a general feeding model for adult Lepidoptera. Oecologia 66, 226–237.
    OpenUrlWeb of Science
  9. ↵
    1. Kim W,
    2. Gilet T,
    3. Bush JWM
    . 2011 Optimal concentrations in nectar feeding. Proc. Natl Acad. Sci. USA 108, 16 618–16 621. doi:10.1073/pnas.1108642108 (doi:10.1073/pnas.1108642108)
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Birchard GF
    . 1997 Optimal hematocrit: theory, regulation and implications. Integr. Comp. Biol. 37, 65–72. doi:10.1093/icb/37.1.65 (doi:10.1093/icb/37.1.65)
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Kingsolver J,
    2. Daniel T
    . 1983 Mechanical determinants of nectar feeding strategy in hummingbirds: energetics, tongue morphology, and licking behavior. Oecologia 60, 214–226. doi:10.1007/BF00379523 (doi:10.1007/BF00379523)
    OpenUrlCrossRefWeb of Science
  12. ↵
    1. Hainsworth F
    . 1973 On the tongue of a hummingbird: its role in the rate and energetics of feeding. Comp. Biochem. Physiol. 46, 65–78. doi:10.1016/0300-9629(73)90559-8 (doi:10.1016/0300-9629(73)90559-8)
    OpenUrlCrossRefPubMedWeb of Science
  13. ↵
    1. Nicolson SW,
    2. Nep M,
    3. Pacinir E
    1. Nicolson SW
    . 2007 Nectar consumers. In Nectaries and nectar (eds Nicolson SW, Nep M, Pacinir E), pp. 289–342. Dordrecht, The Netherlands: Springer.
  14. ↵
    1. Haynes WM
    (ed.). 2012 CRC handbook of chemistry and physics, 93rd edn. Boca Raton, FL: CRC Press.
  15. ↵
    1. Jensen KH,
    2. Savage JA,
    3. Holbrook NM
    . 2013 Optimal concentration for sugar transport in plants. J. R. Soc. Interface 10, 20130055. doi:10.1098/rsif.2013.0055 (doi:10.1098/rsif.2013.0055)
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Robertson AM,
    2. Sequeira A,
    3. Kameneva MV
    . 2008 Hemodynamical flows. In Oberwolfach seminars, vol. 37. Basel, Switzerland: Birkhäuser Basel. doi:10.1007/978-3-7643-7806-6 (doi:10.1007/978-3-7643-7806-6)
    OpenUrlCrossRef
  17. ↵
    1. Gregory T
    . 2000 Nucleotypic effects without nuclei: genome size and erythrocyte size in mammals. Genome 901, 895–901. doi:10.1139/g00-069 (doi:10.1139/g00-069)
    OpenUrlCrossRef
  18. ↵
    1. Seymour R,
    2. Blaylock A
    . 2000 The principle of Laplace and scaling of ventricular wall stress and blood pressure in mammals and birds. Physiol. Biochem. Zool. 73, 389–405. doi:10.1086/317741 (doi:10.1086/317741)
    OpenUrlCrossRefPubMedWeb of Science
  19. ↵
    1. Eckmann DM,
    2. Bowers S,
    3. Stecker M,
    4. Cheung AT
    . 2000 Hematocrit, volume expander, temperature, and shear rate effects on blood viscosity. Anesth. Analg. 91, 539–545. doi:10.1213/00000539-200009000-00007 (doi:10.1213/00000539-200009000-00007)
    OpenUrlCrossRefPubMedWeb of Science
  20. ↵
    1. Jensen KH,
    2. Liesche J,
    3. Bohr T,
    4. Schulz A
    . 2012 Universality of phloem transport in seed plants. Plant Cell Environ. 35, 1065–1076. doi:10.1111/j.1365-3040.2011.02472.x (doi:10.1111/j.1365-3040.2011.02472.x)
    OpenUrlCrossRefPubMedWeb of Science
  21. ↵
    1. Rennie EA,
    2. Turgeon R
    . 2009 A comprehensive picture of phloem loading strategies. Proc. Natl Acad. Sci. USA 106, 14 162–14 167. doi:10.1073/pnas.0902279106 (doi:10.1073/pnas.0902279106)
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. de Gennes P-G,
    2. Brochard-Wyart F,
    3. Quere D
    . 2003 Capillarity and wetting phenomena: drops, bubbles, pearls, waves (Google eBook). New York, NY: Springer.
  23. ↵
    1. Schadschneider A
    . 2002 Traffic flow: a statistical physics point of view. Phys. A Stat. Mech. Appl. 313, 153–187.
    OpenUrlCrossRef
    1. Seibold B,
    2. Flynn MR,
    3. Kasimov AR,
    4. Rosales RR
    . 2012 Constructing set-valued fundamental diagrams from jamiton solutions in second order traffic models. (http://arxiv.org/abs/1204.5510)
    1. Helbing D
    . 2001 Traffic and related self-driven many-particle systems. Rev. Mod. Phys. 73, 1067–1141. doi:10.1103/RevModPhys.73.1067 (doi:10.1103/RevModPhys.73.1067)
    OpenUrlCrossRefWeb of Science
  24. ↵
    1. Flynn MR,
    2. Kasimov AR,
    3. Nave J-C,
    4. Rosales RR,
    5. Seibold B
    . 2009 Self-sustained nonlinear waves in traffic flow. Phys. Rev. E 79, 56113. doi:10.1103/PhysRevE.79.056113 (doi:10.1103/PhysRevE.79.056113)
    OpenUrlCrossRef
  25. ↵
    1. Greenberg H
    . 1959 An analysis of traffic flow. Oper. Res. 7, 79–85. doi:10.1287/opre.7.1.79 (doi:10.1287/opre.7.1.79)
    OpenUrlCrossRefWeb of Science
  26. ↵
    1. Bando M,
    2. Hasebe K,
    3. Nakayama A
    . 1995 Dynamical model of traffic congestion and numerical simulation. Phys. Rev. E 51, 1035–1042. doi:10.1103/PhysRevE.51.1035 (doi:10.1103/PhysRevE.51.1035)
    OpenUrlCrossRefWeb of Science
  27. ↵
    Minnesota Department of Transportation. 2013 Mn/DOT traffic data. See http://data.dot.state.mn.us/datatools.
  28. ↵
    1. McSweeney P,
    2. Fox PF
    1. McCarthy O,
    2. Singh H
    . 2009 Physico-chemical properties of milk. In Advanced dairy chemistry (eds McSweeney P, Fox PF), pp. 691–758. New York, NY: Springer.
  29. ↵
    1. Scherp P,
    2. Hasenstein KH
    . 2007 Anisotropic viscosity of the Chara (Characeae) rhizoid cytoplasm. Am. J. Bot. 94, 1930–1934. doi:10.3732/ajb.94.12.1930 (doi:10.3732/ajb.94.12.1930)
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Wardlaw IF,
    2. Passioura JB
    1. Pate J
    . 1976 Nutrients and metabolites of fluids recovered from xylem and phloem: significance in relation to long-distance transport in plants. In Transport and transfer processes in plants (eds Wardlaw IF, Passioura JB), pp. 253–281. New York, NY: Academic Press.
  31. ↵
    1. Saitô N
    . 1950 Concentration dependence of the viscosity of high polymer solutions. I. J. Phys. Soc. Jpn. 5, 4–8. doi:10.1143/JPSJ.5.4 (doi:10.1143/JPSJ.5.4)
    OpenUrlCrossRef
  32. ↵
    1. Pries A,
    2. Neuhaus D,
    3. Gaehtgens P
    . 1992 Blood viscosity in tube flow: dependence on diameter and hematocrit. Am. J. Physiol. 263, 1770–1778.
    OpenUrl
  33. ↵
    1. Green H
    . 1950 Circulatory system: physical principles. Med. Phys. 2, 228–251.
    OpenUrl
View Abstract
PreviousNext
Back to top
PreviousNext
6 June 2013
Volume 10, issue 83
Journal of The Royal Society Interface: 10 (83)
  • Table of Contents
  • About the Cover
  • Index by author
  • Ed Board (PDF)

Keywords

transport systems
optimal concentration
vascular transport
drinking strategies
traffic flow
Share
Optimal concentrations in transport systems
Kaare H. Jensen, Wonjung Kim, N. Michele Holbrook, John W. M. Bush
J. R. Soc. Interface 2013 10 20130138; DOI: 10.1098/rsif.2013.0138. Published 17 April 2013
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Connotea logo Facebook logo Google logo Mendeley logo
Email

Thank you for your interest in spreading the word on Journal of The Royal Society Interface.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Optimal concentrations in transport systems
(Your Name) has sent you a message from Journal of The Royal Society Interface
(Your Name) thought you would like to see the Journal of The Royal Society Interface web site.
Print
Manage alerts

Please log in to add an alert for this article.

Sign In to Email Alerts with your Email Address
Citation tools

Optimal concentrations in transport systems

Kaare H. Jensen, Wonjung Kim, N. Michele Holbrook, John W. M. Bush
J. R. Soc. Interface 2013 10 20130138; DOI: 10.1098/rsif.2013.0138. Published 17 April 2013

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Download

Article reuse

  • Article
    • Abstract
    • 1. Introduction
    • 2. General formulation
    • 3. Biological transport systems
    • 4. Applications to engineered transport systems: traffic flow
    • 5. Universal properties of transport systems
    • 6. Discussion and conclusion
    • Acknowledgements
    • Appendix A. Viscosity and density of nectar and phloem sap
    • Appendix B. Viscosity of blood
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

See related subject areas:

  • mathematical physics
  • biophysics
  • biomechanics

Related articles

Cited by

Interface Focus To read similar articles, check out our sister journal

Open biology

  • Interface
    • About this journal
    • Contact information
    • Purchasing information
    • Submit
    • Author benefits
    • Open access membership
    • Recommend to your library
    • FAQ
    • Help

Royal society publishing

  • ROYAL SOCIETY PUBLISHING
    • Our journals
    • Open access
    • Publishing policies
    • Conferences
    • Podcasts
    • News
    • Blog
    • Manage your account
    • Terms & conditions
    • Cookies

The royal society

  • THE ROYAL SOCIETY
    • About us
    • Contact us
    • Fellows
    • Events
    • Grants, schemes & awards
    • Topics & policy
    • Collections
    • Venue hire
1742-5662

Copyright © 2018 The Royal Society