PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Circ Res. Author manuscript; available in PMC 2010 November 11.
Published in final edited form as:
Circ Res. 1966 August; 19(2): 332–346.
PMCID: PMC2978656
NIHMSID: NIHMS233693

Plasma Indicator Dispersion in Arteries of the Human Leg

Abstract

Indicator-dilution curves were recorded from the femoral and dorsalis pedis arteries of five normal men after injections of indocyanine green into the superior vena cava or thoracic aorta. By considering the femoral curves as inputs to a mathematically linear system and the dorsalis pedis curves as outputs, transfer functions (the distribution of transit times) for the arterial segment between these sites were obtained in terms of a four-parameter model, the lagged normal density curve, over a sixfold range of flow rates. The parameters of the spread (dispersion) of 57 transfer functions were proportional to the mean transit time. The mean difference between transit time and appearance time was 0.30 t; the square root of the variances was 0.18 t. These linear relationships suggest that flow rate has no significant influence on dispersion and that, since no transition from laminar to turbulent flow was apparent, arterial flow characteristics were not significantly changed over a wide range of flow rates. The secondary implication is that the rate of spatial longitudinal spreading of indicator with distance traveled is primarily a function of the geometry of the arterial system, not of the rate of flow, and, therefore, that the spatial distribution at any instant is a function of this rate and of the distance traveled through the system.

Keywords: arterial blood flow, cardiac output velocity profile, vasodilatation by adenosine triphosphate infusion circulation model, arterial transfer functions distribution of transit times through an arterial segment, indocyanine green indicator dilution, unanesthetized man

Recorded concentration-time curves show the effects of dispersion of indicator at the injection site and in the sampling system. The latter effect can be lessened by the use of sampling systems of high dynamic response, but the former effect is unavoidable. However, circulatory transport through a given segment of the circulation can be characterized experimentally. After injection of indicator at a point upstream from the segment, simultaneous indicator concentration-time curves may be recorded, with identical sampling-recording systems, from the entrance and from the exit of the segment. From these curves, the transfer function of the segment (the frequency distribution of transit times through it) can be estimated.

In this study, the transfer functions were defined in terms of a model—the lagged normal density curve.1 The flow of blood through the leg was increased by the infusion, at varied constant rates, of adenosine triphosphate in order to assess the influence of flow rate on these transfer functions.

Methods and Computations

Five normal men were studied while they were awake in the supine position. The details of catheter sizes and locations, the methods of sampling, the calibration of the densitometers, and the recording and preliminary reduction of data have been described.1 A transfer function for the arterial segment lying between the femoral artery and the dorsalis pedis artery was obtained by considering the concentration-time curve recorded at the femoral artery as the input to a linear system and that recorded at the dorsalis pedis artery as the output.

MATHEMATICAL APPROACH TO THE USE OF THE CONVOLUTION INTEGRAL

The approach used in this study is essentially that of Stephenson:2 the arterial system is considered subject to the types of analysis applicable to linear filters with lumped characteristics. It is assumed that the transfer function has its domain along only the time axis. This implies the secondary assumptions that each sampling site acts as a point and that the concentration at that point represents the amount of indicator per volume of fluid passing that point (that is, the concentration is volume-averaged)3, 4 or is the mean flow concentration.5 The transfer function of a linear system is a unique expression, h(t), describing the gain and distortion of any signal passing through the system. h(t) is the response of the network to a unit impulse as input. Given h(t) and any input, f(t), the output or system response, g(t), has a unique description:

g(t)=0th(tλ)f(λ)dλ
(1)

or

g=f*h=h*f
(1a)

in which λ is a dummy variable of integration and * is an abbreviation for the convolution. Equation 1 is the convolution integral, limited to systems which do not respond unless there is an input and in which f(t) = 0 for t ≤ 0. The present experimental situation is represented in figure 1. Points 1, 2, and 3 represent the injection site and the sampling sites in the femoral artery and dorsalis pedis artery, respectively. Points 4 and 5 represent the densitometers whose outputs are the recorded concentration-time curves, Ci (t) and Co(t), respectively. L1; the vascular system between the injection site and the femoral sampling site, has a transfer function, h1(t), which is the response to a unit impulse, δ(t). Similarly, L2, the arterial segment between the two sampling sites, has a unit impulse response, h2(t). It is the purpose of these experiments to determine h2(t) under various circumstances. This cannot be done directly because a unit impulse cannot be introduced into the femoral artery (point 2) nor can the concentration, g(t), at the dorsalis pedis artery (point 3) be obtained directly. The form of the injection cannot be characterized with current methods; therefore, the input function, e(t), is unknown and the impulse response, h1(t), of system L1 cannot be determined. f(t) and g(t) cannot be recorded directly but are distorted by the sampling systems, LS1 and LS2.

FIGURE 1
Block diagram of the experimental situation. Point 1 (left) is the injection site, 2 is the tip of the femoral sampling system, 3 is the tip of the sampling system in the dorsalis pedis artery, and 4 and 5 are the recorded outputs of the densitometers ...

The recorded curves are the arterial curves distorted by the sampling systems:

Ci(t)=f(t)*S1(t)
(2)

and

Co(t)=g(t)*S2(t)
(3)

in which S1(t) and S2(t) are the sampling-system responses to a unit impulse. In these experiments the sampling systems have been designed to have nearly identical responses to a step change of dye concentration at the input1 and, therefore, their time-domain transfer functions, S1(t) and S2(t), are very nearly identical:

S1(t)=S2(t)=S(t).
(4)

This simplifies the analysis considerably because then Co(t) = g(t) * S(t). Substituting for g(t) by using equation la,

Co(t)=h2(t)*f(t)*S(t)
(5)

and then from equation 2,

Co(t)=h2(t)*Ci(t)
(6)

or more properly:

Co(t)=+lh2(λ)Ci(tλ)dλ.
(7)

In the experimental situation the injection of indicator at time zero resulted in the concentrations C1 and Co, and therefore the lower limit of integration may be changed to ignore the zero concentrations preceding the injection:

Co(t)=0th2(λ)Ci(tλ)dλ.
(8)

The experimental determination of h2(t), hence-forth called simply h2(t), depends on the veracity of equation 8, which depends on S1(t) and S2(t) being equal and on the assumption of mathematical linearity. A consistent difference between S1(t) and S2(t) will result in a systematic error in the transfer function, but it was thought that the 0.25-second difference in the mean transit times of S1(t) and S2(t) present in this study was sufficiently small relative to an average mean transit time of 9 seconds for the transfer functions of these experiments that it could be ignored.

The determination of h2(t), is entirely dependent on the applicability of the superposition theorem. The requirements are that the system be linear and stationary.5, 6 (1) Stationarity means constant flow. (2) The indicator must travel in exactly the same fashion as the substance which it is considered to label. (In this experiment, plasma protein or perhaps simply plasma is labelled by indocyanine green but erythrocytes are not.) (3) There must be no loss or gain of indicator or its carrier during passage through the system between the femoral and dorsalis pedis arteries. (4) The sampling of the bloodstream for indicator concentration should be not only from all of the cross section of the vessel at the tip of the sampling system but also from each portion of this cross section in proportion to the flow across it (flow-averaged sample).5 (5) The system must also be linear such that, when

f1(t)g1(t)

and

f2(t)g2(t)

then

a1f1(tt1)+a2f2(tt2)a1g1(tt1)+a2g2(tt2).

Experimentally, t1 or t2 cannot be zero and therefore a test of linearity using the latter equation is dependent on the persistence of steady flow over a period of at least a few minutes.

While it is obvious that none of these conditions can be completely fulfilled, the deviations from ideal are not so great as to vitiate the practical application of superposition theory.7

COMPUTATION OF THE CONVOLUTION INTEGRAL

Digital computation of integrals of continuous functions is done by means of numerical approximations. The accuracy of a numerical integration is increased by decreasing the size of the intervals chosen. For these experiments, the interval, Δt, was 0.25 second which was 1/35 to 1/100 of the passage time of the primary curves. The functions f(t), h(t), and g(t) are considered as sequences fi, hi, and gi where i is an integer subscript valued 1 to n. fi is the amplitude of a pseudo-impulse and substitutes for the area of f(t) during a period, Δt, whose midpoint is ti.

fi=tiΔt/2ti+Δt/2f(t)dt=f¯(ti)Δt.

The convolution integral (equation 1) was calculated in the recursive form:

g1=f1h1g2=f2h1+f1h2g3=f3h1+f2h2+f1h3

gi=f=1f=ifjhij+1.
(9)

This Euler integration was found to be as accurate for these experiments as a Simpson’s rule or a Runge-Kutta integration and was very much faster.

ESTIMATION OF THE TRANSFER FUNCTION

The transfer function appropriate to each pair of recorded curves, Ci(t) and Co(t), was determined in terms of a suitably defined lagged normal density curve1 whose equation is:

h(t)=0t1σ2πe1/2(λtoσ)21τe(tλ)τdλ,
(10)

in which tc is the median of a normal density curve of standard deviation, σ, and τ is the time constant of a superimposed exponential (first order) lag. This model, which is the same one as that used by Nicholes and co-workers,8 was generated using a Runge-Kutta method. Values of h(t) for t < 0 were ignored and, in practice, values were less than 0.5% of the peak value until t > 0.5 (τ + tc). The curve was arbitrarily truncated on its downslope following the peak when values were less than 1% of the peak. The first, second, and third moments of h(t) are, respectively, τ + tc, σ2 + τ2, and 2τ8.

It should be emphasized that this model is used to provide a concise numerical description of the transfer function. It is not unique, for several other models could be used, and there is no direct means of interpreting its parameters in physical terms, such as a velocity profile or other blood-flow characteristics. It is simply a stochastic, descriptive model which has been chosen on the basis of its mathematical similarity to equations for turbulent diffusion and the random walk and because, empirically, it does provide a sufficiently accurate description of the time-domain transfer function.

A trial transfer function, h(t), was generated by using arbitrary values for σ, τ, and tc and had values at intervals of 0.25 second. A trial output curve, C′o(t), was computed by the convolution of Ci(t) and h(t), by equation 9; C′o(t) was compared with Co(t) by determining the coefficient of variation. If the coefficient of variation was higher than a prescribed value, usually 0.005, then h(t) was changed by adjusting the parameters of the lagged normal density curve, the convolution was repeated, and the new C′o(t) was compared with Co(t). This process was repeated either until the coefficient was less than the prescribed value or until 10 trials were completed, at which point the best result was printed out and plotted. The adjustment of the parameters was made automatically by the computer: if C′o(t) was too steep on the upslope, σ was increased; if C′o(t) was too steep on the downslope, τ was increased; if the peak was too early, tc was increased. Interaction between the effects of these adjustments was taken into account (for example, an increase in τ also decreases the upslope and delays the peak) but there still tended to be a small oscillation about the best solution. This was seen when the prescribed goal for the coefficient of variation was set at an impossibly low value such as 0.002.

Results

USE OF THE LAGGED NORMAL DENSITY CURVE AS THE TRANSFER FUNCTION

The transfer function was determined in terms of a lagged normal density curve for 29 pairs of curves recorded after injection of indicator into the superior vena cava and for 28 pairs after injection into the aorta. Typical results are shown in figure 2. In each panel the computed downstream curve, C′o(t), represented by the + signs, is superimposed on the continuous line representing the dye curve recorded at the dorsalis pedis artery, Co(t). The parameters of the transfer function are listed in each panel. The coefficients of variation between C′o(t) and Co(t) were significantly (P < 0.001) less for pairs of curves recorded after injection of dye into the superior vena cava (mean coefficient of variation = 0.0053, sem = 0.0004) than for those recorded after aortic injection (0.0094 ± 0.0008). The average for the 57 transfer functions was 0.0070 ± 0.0005 with a range from 0.002 to 0.020. In only four instances was the coefficient of variation greater than 0.011.

FIGURE 2
Use of a lagged normal density curve as the transfer function relating pairs of time-concentration curves. The pair of curves (continuous lines) in each panel were recorded simultaneously from the femoral and dorsalis pedis arteries. Femoral curves were ...

Under conditions of stationarity of flow, the transfer function of the segment during the period of recirculation should be the same as that for the dye passing the sampling sites for the first time. The better-than-average result in the left upper panel of figure 2 shows this most clearly. Following aortic injection, the recirculation peaks were smaller because smaller amounts of dye were injected and because a larger proportion of the indicator goes into the vasculature of the lower part of the body, from which recirculation is much slower.9 Adjustment of the parameters of the transfer function was made in accordance with comparisons of the primary portions of curves Co(t) and C′o(t). As a result the later parts of the curves, representing recirculated indicator, did not match as well. This may be due to variability in the physiologic state (change of flow or distribution of flow) or perhaps variation in sensitivity or drift in base line of the densitometers during the recording.

The ratios of areas of the femoral curves to those at the dorsalis pedis artery averaged 0.96 (sem = 0.011). The areas should be identical and the variation is an expression of the experimental error. Differences could be due to lack of stationarity, instrumental errors, or, least likely, to a difference in the error due to time-averaged sampling instead of volume-averaged sampling at the two sites.10 No significant change in sensitivity of the densitometers to known concentrations of dye in blood from the subjects was apparent on tests done periodically during each experiment. Lack of representative mixing was probably a minor factor because the variation in area was slightly greater following injections into the aorta. The systematic difference was due primarily to the error in extrapolation of the femoral curves (see fig. 5 of reference 1).

Lack of stationarity is a probable source of much of the variation noted in this study. Positive evidence of a change in flow is seen in figure 3. These data were recorded while the subject’s physiologic status was apparently undisturbed and unchanging (no changes in heart rate, blood pressure, or respiratory rate). A sudden injection of dye was made into the thoracic aorta, and the upper pair of curves was recorded. The transfer function had the shape represented by the open circles in the left upper panel. The computed output curve, C′o(t), is superimposed on the recorded curve, Co(t), in the right upper panel (coefficient of variation = 0.006). Fifty seconds after the first injection, two injections of dye were made three seconds apart into the aorta, and the curves shown in the lower panels were recorded. As above, the transfer function, h(t) (labelled “dispersion process”), and a theoretical output curve, C′o(t), were computed, and C′o(t) was compared with Co(t) in the lower right panel. The coefficient of variation was 0.010. Lack of stationarity is shown by the difference in mean transit time: τ + tc equalled 14.5 seconds initially (upper panel) but was 13.5 seconds, less than a minute later (lower panel). If the change in flow occurred in the few seconds between these recordings, then both transfer functions may be valid. If a gradual change was occurring, then the error in each will be small. If the change occurred during the recording of one of the primary peaks the error would be quite significant. Differences of 0.5 to 1.0 second in t were found fairly frequently even when the determinations were made at intervals as short as 50 seconds. This reflects presumably a continual variation in total blood flow in the leg or of distribution of flow in the leg vessels.

FIGURE 3
Test of constancy of the transfer function at normal femoral flow rate. Left panels: curves recorded from the femoral artery (asterisks) and the calculated transfer function (zeros) whose parameters are listed and whose amplitude was 1. Right panels: ...

TEST OF THE MATHEMATICAL LINEARITY OF THE ARTERIAL SEGMENT

The data of figure 4 were obtained in a manner similar to those in figure 3, except that adenosine triphosphate was being infused into the right common iliac artery at a rate of 1.0 mg/min, producing a marked decrease in mean transit time and in dispersion. The mean transit time between sampling sites was the same (3.1 sec) after a single injection as it was, less than a minute later, after a pair of injections 4 sec apart. If this is taken as evidence for stationarity, then the curves may be used to test whether the transfer function is the same for input curves of widely differing form. The transfer function for the curves of the upper panel is described by a lagged normal density curve having σ = 0.20 sec, τ = 0.26 sec, and tc = 2.8 sec. That for the doubly peaked curves of the lower panels had σ = 0.19 sec, τ = 0.31 sec, and tc = 2.75 sec. The coefficients of variation were 0.007 and 0.011. The conclusion is that, when stationarity is present, linearity was demonstrated by showing that the transfer function was independent of the shape of the recorded concentration-time curves. During infusion of adenosine triphosphate, stationarity was usually observed and five pairs of computed transfer functions obtained in the above manner were similarly independent of the form of the dye dilution curves. During the control state there were no consecutive pairs of transfer functions obtained at precisely constant flow or mean transit time and therefore linearity could not be tested.

FIGURE 4
Test of the constancy of the transfer function at a high femoral flow rate (during infusion of adenosine triphosphate into the common iliac artery at a rate of 1.0 mg/min). The situation is similar to that of figure 3. Lower panels: two injections 4 sec ...

RELATIONSHIPS BETWEEN PARAMETERS OF THE TRANSFER FUNCTION

The transfer function was determined only in terms of the parametric model, the lagged normal density curve. Figure 5 (left panel) shows that the square root of the variance, (σ2 + τ2)½, is linearly proportional to the mean transit time, τ + tc. That the shape was rather variable is shown by the relationships of σ and τ to τ + tc (fig. 5, middle and right panels). The least-squares regression lines, calculated on the basis that error exists in values of ordinates and abscissas, have approximately the same slopes as observed1 for the same relationships for the model fitted to the recorded dye curves. The large variation is a reflection of the fact that the shape of the transfer function is a caricature of the difference in the shapes of Ci(t) and Co(t). The parameters σ and τ were unrelated (τ = 0.06).

FIGURE 5
Parameters of the model simulating the transfer function. Left panel dispersion, as estimated by the square root of the variance of the transfer function, (σ2 + τ2)½, is linearly related to the mean transit time, τ + t ...

The theoretical values for t and π2½ (the square root of the variance) were used in figure 5. The actual calculated moments gave slightly smaller values because the curve ordinates were computed at finite intervals (0.25 sec) and h(t) was terminated when its value on the downslope was less than 1.5% of the peak value. As calculated, τ + to = 0.23 + 1.008 t with a standard deviation of 0.033 and a correlation coefficient, r, of 0.999; similarly, (σ2 + τ2)½ = − 0.08 + 1.13π2½ (sd = 0.051, τ = 0.998). These relationships should fall on the line of identity; such errors are more obvious with the higher moments.

The regression lines for the appearance time, ta and tta versus t for h(t) show relationships (fig. 6) similar to those obtained for the recorded curves. These lines might be expected to pass through the origin because there should be no influence from dispersion at the injection site or upstream from the arterial segment under study. That they do not, suggests that there may be some dependence of shape on the mean transit time. tta is proportionately larger when t is smaller, suggesting that h(t) is more symmetric when t is small. The same tendency is shown in the relationship between tta and π2½ (fig. 6, right panel) and in the interrelationships between σ, τ, tta, and π2½. For example, σ = −0.07 + 2.24(tta) (sd = 0.24, r = 0.73); this is very close to the ratio, σ/(tta) = 0.26. Thus, the relationship σ = 0.45 + 0.19 π2½ (sd = 0.33, r = 0.38) must and does show a positive ordinate intercept, while values of τ were relatively small when t and π2½ were small (fig. 5, right panel). τ = 0.51 + 1.23 π2½ (sd = 0.30, r = 0.94) and τ = −0.92 + 0.82(tta) (sd = 0.60, r = 0.72).

FIGURE 6
Dispersion of the transfer function. The relationships of ta and (t — ta) to the mean transit time between the sampling sites are similar to those relating the corresponding parameters of the recorded curves, but the scatter is ...

Table 1 reflects the same tendency. The variations in σ/t and in τ/t were so great that no significant differences due to injection site were seen. The average ratios are similar to those obtained previously1 for recorded dye curves but are smaller following aortic injections, π2½/t was significantly (P < 0.02) less with aortic injections. Similarly, the average (tta)/ t was smaller and ta/t was larger for aortic injections.

TABLE 1
Ratios of Parameters of the Transfer Function

COMPARISON OF TRANSFER FUNCTION WITH RECORDED CURVES

The shape and variance of the transfer function may be compared with the recorded curves. Figure 7 illustrates the relationships of the first and second moments of the transfer function to those of the recorded curves. They should follow the line of identity. The systematic deviation of the second moments from the relationship, π2h = π2DP — π2F (the subscripts DP and F indicate the downstream [dorsalis pedis] and upstream [femoral] curves, respectively), is due in part to integration error and to truncation error in the calculation. The scatter observed may be expected in a comparison of an error-containing parameter, (σ2 + τ2)½, with the relatively small difference between two large and some-what erroneous values, π2DP and π2F.

FIGURE 7
Relationships between moments of the transfer function and moments of the recorded curves. The systematic deviation of the relationship between the variances, the second moments, is partly due to error in computing π2.

Discussion

LAGGED NORMAL DENSITY CURVE AS A MODEL

The purpose of these experiments was to compute a transfer function which was independent of the influences of dispersion at the site of injection and of distortion by the sampling system. Since these factors unavoidably affect the recorded concentration-time curve, use of this curve as the transfer function, as Zierler11 advocated, serves only as a rough approximation.

In order to ascertain the transfer function without using transforms, h(t) may be assumed and adjusted by some trial and error process until a suitable h(t) is found. Tedious as this is, such a method has certain advantages because it avoids the division of the frequency-domain transform of the recorded output curve by that of the input curve. Such a division is akin to differentiation, and when the data are not smooth but have random or systematic variation superimposed on the ideal Co(t) and Ci(t), the solution tends to have wide fluctuations and may oscillate. Stephenson2 predicted this, and it has also been observed by Paynter12 and by Parrish and co-workers.13 The use of the iterative computation of the output function in the time domain (equation 9) is akin to integration or the use of linear filters and, because h(t) is smooth, it results in C′o(t) being smoother than Ci(t) or Co(t).

A model is used because it provides convenience in defining and adjusting h(t). The lagged normal density curve is thought to be easier to adjust than other models whose forms are not visualized so readily from the parameters. It is a better model for the transfer function than it is for recorded curves. For the 57 determinations, the average coefficient of variation between Co(t) and C′o(t) was 0.007, which appears to be almost an order of magnitude better than the 0.056 found when comparing the model to the recorded curves.1 The improvement is not so great as these numbers would imply because the calculation comparing C′o(t) to Co(t) uses the data points beginning at the appearance time of Co(t) and ending at the last data point during the recirculation; this produces a relative increase in the denominator and reduces the coefficient of variation.

Various probability models could very likely be used as a transfer function in the manner demonstrated here. This has not been attempted in this study, but the cascaded delay lines used by Parrish and co-workers13 form transfer functions very similar to a normal density curve. The response of a single delay line to a step input is a symmetric sigmoid curve which is remarkably similar to a normal distribution curve.* Theirs is a complex model but, since each delay line may be roughly approximated by a normal density curve, a cascade of n such units in series may be described by:

h(t)=i=1i=nAiσ2πie1/2(tiΔtσi)2
(11)

in which i is the position of the delay line in the series, Ai is the coefficient for each delay line and i=1i=nAi=1.0, Δt is the mean delay time of one unit, and σ is the square root of the variance of h(t) for one unit. Thus, the system is describable by the n + 2 parameters, Ci, Δt, and σ. If the nonzero coefficients are all equal, then h(t) will be symmetric. Skewing to the right was produced by having higher coefficient settings for the longer delay lines than for the shorter. If only one coefficient were used, then the model would be virtually identical to the lagged normal density curve. The virtuosity of the Parrish model could be greatly extended and the number of parameters reduced by the inclusion of exponential lags and, even in its present state, could encompass the situation in which there are parallel systems with different mean transit times. This, of course, can also be done with the lagged normal density curve by using the sum of two curves as a transfer function. Recirculation to the heart by separate systemic pathways can be described only in such terms.7,8,14,15

MATHEMATICAL LINEARITY OF THE VASCULAR SYSTEM

The most formidable requirement for linearity is stationarity.5 Zierler16 has said, “If phasic, not necessarily regular, alterations in distribution of transit times fluctuate rapidly about some central value, and if the periods of the phases are brief compared to the time required for evolution of the sudden-injection indicator concentration-time curve, then the violation of stationarity may not be important.” The same condition reduces the error induced by time-averaged sampling of the flowing blood. It is then not surprising that the coefficients of variation were smaller for the broad curves following injections into the superior vena cava, because the passage time of the aortic curves may be as short as 5 to 7 sec.

Variation in vasomotor activity and in peripheral flow, as observed by Burton17 and Allwood and Burry,18 is a more significant source of error. Figure 3 illustrates that this does not necessarily prevent the calculation of transfer functions, but their validity must be doubted because it is not known when the change of flow occurred. In certain instances the flow must have increased during the period the indicator was passing the distal sampling site because Co(t) was narrower than Ci(t) and the calculation of a transfer function was ridiculous. For such reasons, similarity of area is a good criterion for selection of pairs of curves from which to determine the transfer function.

Nicholes and co-workers8 described a test for linearity, with which they obtained a satisfactory result. Injections of 0.5, 1.0, and 2.0 times the usual amount of dye were made into the superior vena cava, and the blood was sampled from the pulmonary artery. The resulting concentration-time curves had the same shape, and their areas were in proportion to the dye dose. This test confirms both the laws of conservation of material and the basis for the use of dye techniques to measure flow, but it is not a critical test of linearity because the shape of the primary portion of the dye curve is so ubiquitous. Transfer functions for indicator recirculating to the pulmonary artery were calculated from the primary curve (Ci) and the “first recirculation” curve (Co) but were inconstant, which probably reflects the lack of stationarity of flow in the peripheral vascular system.

The test shown in figure 4 is a good test of linearity over a specific vascular segment but, since the two injections (lower panels) were 4 sec apart and the transfer function was short, the two peaks remained separate throughout the segment. The shape of the curves changed very little from those in the upper panels and therefore the criticism must be made that there is not a great difference in the form of Ci(t) in the upper and lower panels. The assumption of linearity has been justified by experiments in which transfer functions were obtained simultaneously for three segments of the aorta in dogs.7

INTERRELATIONSHIPS BETWEEN PARAMETERS OF THE TRANSFER FUNCTION

The curves recorded at the two sampling sites have been shown1 to be basically similar in shape. The transfer function should therefore also have the same shape. The output curve, C′o(t), computed by means of the convolution integral, must always have a degree of skewness or kurtosis intermediate between those of the input curve Ci(t) and the transfer function h(t), and it will be more spread out than either. For example, the lagged normal density curve is the convolution of a normal density curve and a single exponential curve; it is neither as skewed as the exponential nor as symmetric as the Gaussian curve. Therefore, when Ci(t) and Co(t) differ in shape, the transfer function, h(t), must be more skewed than Co(t) if Co(t) is more skewed than Ci(t). Or, if σ/t of the model fitting Co(t) is greater than is σ/t of the model fitting Ci(t), then σ/t of h(t) must be greater than τ/t of Co(t) and τ/t of h(t) must be less than τ/t of Co(t). The statement regarding τ follows from the fact that the variance of Ci(t) plus the variance of h(t) should equal the variance of Co(t). In other words,

π2f+π2h=π2g
(12)

or

(σf2+τf2)+(σh2+τh2)=(σg2+τg2)
(13)

coupled with the empirical observation that π2½ = 0.18 t, which holds approximately for both the recorded curves and the transfer function (fig. 5 and table 1). This phenomenon of the transfer function exaggerating the differences between Ci(t) and Co(t) is similar to the exaggeration of irregularities produced by differentiating a function. This is the main reason for the poor correlation between σ or τ and t (fig. 5), and it results in the apparently random relationship between σ and τ.

The relationships of ta and (tta) to t (fig. 6) show much less scatter than do those of σ or τ to t. The difference between the average transit time and the shortest transit time is closely related to the variance of the transfer function (right panel) and a quick and approximate estimate of the variance could be obtained by the equation π2½ = 0.55 (tta). This value is a rough average for the relationships of figure 5.

The observation that π2½/t and (tta)/t had smaller average values (table 1) following aortic injection might suggest that the dispersion in the arterial system of the leg is influenced by the site of injection. Such a thought is contrary to physiologic and mathematical expectations, and the deviation is probably due to a small systematic experimental error.

RELATIONSHIPS OF PARAMETERS OF THE TRANSFER FUNCTION TO PARAMETERS OF THE RECORDED CURVES

In table 2 are listed the average values of the slopes of the relationships of various parameters to the mean transit times of the recorded curves1 (left column) and the slopes obtained for the parameters of the transfer function (right column). These slopes and the ratios of the transfer function parameters to t (next to last column, table 1) should be influenced only by dispersion within the vascular segment between the sampling sites. The similarity in the values indicates that a rather good prediction of the transfer function could be made from the parameters of the recorded curves.

TABLE 2
Comparison of Slopes Relating Parameters to Mean Transit Times for Recorded Curves and for Those of the Transfer Function

The relationship between the transfer function mean transit time and the difference in mean transit times between the curves recorded from the densitometers is obviously a necessary one and should fall on the line of identity if the sampling systems are identical (fig. 7, left panel). Some variation is expected because the three mean transit times necessary to plot each point are calculated independently. Greater error is found in the estimates of variance principally because of errors in the tails of the curves, and the variation from the line of identity (equation 12) is larger (fig. 7, right panel). Computation error is greater for the higher moments and produces a systematic underestimation.

SPATIAL AND TEMPORAL DISTRIBUTION OF INDICATOR PARTICLES

Because of the proportionality of spread to mean transit time and because of the small influence of flow on this relationship, certain statements can be made about C(x, t), the distribution of the indicator particles with respect to distance x along the system. In general, if any two of the following are known, then the third may be calculated: (1) C(x1, t) at a point in the system x1; (2) C(x1, t) the spatial distribution of concentration with respect to distance at a particular time t1; (3) the rate of increase of spread of C(x) with time. This is somewhat oversimplified since details such as cross-sectional versus flow sampling must be considered also.

One way of looking at these phenomena is to consider a system in which there is constant proportionality of distribution of flow. Such a system might have turbulent or laminar flow but the simplest example is a system (fig. 8) consisting of a number of parallel pathways, each having piston flow and each carrying a constant proportion of the total flow. (Piston flow defines the velocity to be the same at all points in a cross section of the tube.) Indicator injected at point A in such a system will be dispersed, if the pathways have varied mean transit times, at a given total flow rate, Q. It is obvious that, under these conditions, when an impulse input is made at A, the spatial distribution, C(x), will be the same around B no matter what the flow rate is. The temporal distribution at the point B, C(t), is a function of the spatial distribution and of the flow rate. Let the curve at the lower left comer of the figure represent h1(t) when the flow is stationary at Q1 ml/sec. The amount of indicator and the time for it to travel from A to B through pathway D is given by the rectangle at time, tD1.

FIGURE 8
Diagrammatic representation of the effect of change of flow rate in a generalized flow system. When the rate is doubled, the dilution curve is halved in area and the transit time through any particular path, D or E, between A and B is also halved. In ...

In a geometrically stable system the flow rate and mean transit time are related reciprocally:

V=Q1t¯1=Q2t¯2
(14)

in which V is the Stewart-Hamilton or the mean transit time volume of the system and need not be defined anatomically. Equation 14 is applicable to the portion of the circulatory network (including branches) in which there is constant proportionality of distribution of flow. As a result, the transfer functions, h(t), at different flow rates are related by:

h2(t)=t¯1t¯2h1(t¯1t¯2t)=Q2Q1h1(Q2Q1t).
(15)

The scale factor, t1/t2 = Q2/Q1, is necessary to keep the transfer function at unit area. When the concentration-time curves are related in the same fashion:

C2(t)=C1(t¯1t¯2t)=C1(Q2Q1t)
(16)

then the peak heights, C1p and C2p, are the same and the areas are inversely related to the flow, Q. Similarly, C(tD2) = C(tD1) and C(tE2) = C(tE1) while tD2 = Q1tD1/Q2, and tE2 = Q1tE1/Q2. This may be viewed merely as normalization of the time scale to t/t.

The model illustrated in figure 8 is a generalized streamline flow model. The parameters of the spread of C(t), viz., π2½, tta, σ, τ, and any others, are in constant proportion to t. These relationships will apply to any constant-volume streamline flow model over any range of flows so long as the requirement for constant proportionality of flow in each stream is fulfilled. They, and also equations 14, 15, and 16, apply equally well to any turbulent or disturbed flow system of constant volume in which the mean velocity profile is unchanged by flow changes. Because the systemic pressure was constant in each subject despite large changes in femoral artery flow, it is not surprising that the arterial volume was quite constant. It is unlikely that the highly distensible venous system would have so little change in volume over such a range of flows.

The restriction of geometric stability and constant proportionality of distribution may be relaxed to some extent in an open-ended system such as that between the femoral artery and the dorsalis pedis artery. Stationarity need only be maintained in the main arterial stream between these two points. Thus, the distribution of flow to each primary branch must be constant, but in secondary branches the flow need not be clearly defined. The boundaries of the volume, V, may be changed at will within the second and higher order branches of the system.

It is important to realize that streamline flow is not necessarily implied by the observations of this study. The generalized flow model has been discussed because it is simple and because streaming flow is probable in the vascular system, even though true laminar flow is not, because of the presence of erythrocytes. Any type of turbulent flow whose spatial dispersion is related linearly or nonlinearly to the distance traveled and not to the flow will give similar results. The six-fold range of mean transit times of h(t)′s indicates a similar range of flows, the highest of which would result in Reynolds numbers of about 2,000 in the femoral artery if the flow were not pulsatile. It seems unlikely that laminar flow could persist at such flows even if the artery were a uniform cylinder. The similarity of the relationship of spread to mean transit time at low flows to that at high flows, where turbulence is almost certainly present, suggests that flow is usually turbulent. Turbulence may be expected at relatively low flows of nonhomogeneous fluids driven by a pulsatile head of pressure through elastic, branched, tapering, curved tubes. On the other hand, “complete” turbulence resulting in constant concentration over the cross section of the artery and in longitudinal dispersion of indicator seems little more probable than Newtonian laminar flow. Most likely, a compromise model embodying longitudinal and lateral dispersion superimposed on some rather blunt velocity profile will turn out to be the most realistic.

ROLE OF CIRCULATORY DISPERSION IN STUDIES ON DIFFUSIBLE INDICATORS

After injection of a bolus containing diffusible and nondiffusible indicators, simultaneous concentration-time curves for the several substances may be obtained by sampling from the venous outflow of an organ. Such data can be used to estimate the volume of distribution of each indicator and rates of transfer of diffusible indicators from the bloodstream into the tissue.1921 Combination of this knowledge with information on the mode of transport within various vascular segments should produce relatively complete descriptions of the movement of a substance from one tissue to another, and complex models of the whole body or of single organ systems could be formed without some of the simplifying assumptions made by Jacquez and co-workers. 2223

Acknowledgment

Dr. Earl H. Wood not only made this work possible but also aided in the initiation and execution of the experiments and in the revision of the manuscript The author also thanks Dr. Homer R. Warner, Department of Biophysics and Bio-engineering, Latter Day Saints Hospital, University of Utah, for his role in starting the project. Mrs. Jean Frank and her co-workers assisted in the preparation of the illustrations. Dr. K. K. Nicholes helped in developing the initial programs.

Supported in part by Research Grants HE-04664 and FR-00007 from the National Institutes of Health, U. S. Public Health Service.

Footnotes

Reprints: Information about reprints can be found online at http://www.lww.com/reprints

Dr. Bassingthwaighte is a recipient of a Research Career Development Award of the U. S. Public Health Service.

*For graphic illustrations, see “Time domain synthesis based on a passive tardigrade module.” In Lightning Empiricist: 11 (3) July, 1963. Boston, G. A. Philbrick Researches, Inc.

References

1. Bassingthwaighte JB, Ackerman FH, Wood EH. Applications of the lagged normal density curve as a model for arterial dilution curves. Circulation Res. 1966;18:398. [PMC free article] [PubMed]
2. Stephenson JL. Theory of the measurement of blood flow by the dilution of an indicator. Bull. Math. Biophys. 1948;10:117. [PubMed]
3. Visscher MB, Johnson JA. The Fick principle: Analysis of potential errors in its conventional application. J. Appl. Physiol. 1953;5:635. [PubMed]
4. Stow RW. Systematic errors in flow determinations by the Fick method. Minn. Med. 1954;37:30. [PubMed]
5. González-Fernández JM. Theory of the measurement of the dispersion of an indicator in indicator-dilution studies. Circulation Res. 1962;10:409. [PubMed]
6. Meier P, Zierler KL. On the theory of the indicator-dilution method for the measurement of blood flow and volume. J. Appl. Physiol. 1954;6:731. [PubMed]
7. Bassingthwaighte JB, Ackerman FH, Essex HE. Mathematical linearity of the dog’s arterial system. Federation Proc. 1965;24:335.
8. Nicholes KK, Warner HR, Wood EH. A study of dispersion of an indicator in the circulation. Ann. N. Y. Acad. Sci. 1964;115:721. [PubMed]
9. Nicholes KK, Warner HR, Toronto AF, Bassingthwaighte JB. Effect of exercise on the distribution of transit times of an indicator through various routes in the systemic circulation. Federation Proc. 1963;22:523. (abstr.)
10. Bassingthwaighte JB, Edwards AWT, Wood EH. Areas of dye-dilution curves sampled simultaneously from central and peripheral sites. J. Appl. Physiol. 1962;17:91. [PubMed]
11. Zierler KL. Circulation times and the theory of indicator-dilution methods for determining blood flow and volume. In: Hamilton WF, Dow P, editors. Handbook of Physiology. sec. 2. vol. 1. Washington, D. C.: 1962. pp. 585–615. Circulation, Am. Phys. Soc.
12. Paynter HM. Methods and results from M.I.T. studies in unsteady flow. Boston Soc. Civil Engr. J. 1952;39:120.
13. Parrish Daniel, Hayden DT, Garrett Wayne, Huff RL. Analog computer analysis of flow characteristics and volume of the pulmonary vascular bed. Circulation Res. 1959;7:746. [PubMed]
14. Sheppard CW. An electromathematical theory of circulatory mixing transients. Proceedings of the First National Biophysics Conference; Yale University Press; 1959. pp. 476–492.
15. Sato T, Takahashi K, Dohi Y. The IV World Congress of Cardiology. Mexico City: 1962. Sep, Critical study of the Hamilton-Stewart principle for hemodynamic analysis; pp. 312–313.
16. Zierler KL. Theoretical basis of indicator-dilution methods for measuring flow and volume. Circulation Res. 1962;10:393.
17. Burton AC. The range and variability of the blood flow in the human fingers and the vasomotor regulation of body temperature. Am. J. Physiol. 1939;127:437.
18. Allwood MJ, Burry HS. The effect of local temperature on blood flow in the human foot. J. Physiol. 1954;124:345. [PubMed]
19. Chinard FP, Enns T, Nolan MF. Indicator-dilution studies with “diffusible” indicators. Circulation Res. 1962;10:473. [PubMed]
20. Goresky CA. A linear method for determining liver sinusoidal and extra-vascular volumes. Am. J. Physiol. 1963;204:626. [PubMed]
21. Zierler KL. Theory of use of indicators to measure blood flow and extracellular volume and calculation of transcapillary movement of tracers. Circulation Res. 1963;12:464.
22. Jacquez JA, Bellman R, Kalaba R. Some mathematical aspects of chemotherapy. II. The distribution of a drug in the body. Bull. Math. Biophys. 1960;22:309.
23. Bellman R, Jacquez JA, Kalaba R. Some mathematical aspects of chemotherapy. I. One-organ models. Bull. Math. Biophys. 1960;22:181.