Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

The Physics of Coronory Blood Flow - M. Zamir

.pdf
Скачиваний:
157
Добавлен:
10.08.2013
Размер:
5.78 Mб
Скачать

180 6 Composite Pressure-Flow Relations

and opposition to flow consists of only the resistance R, the corresponding flow wave is given by (Eq. 4.8.7)

Q(t) =

P0

 

cos ωt

(6.2.2)

R

 

 

 

 

 

=

P (t)

when P (t), Q(t) are sinusoidal waves

(6.2.3)

R

 

 

 

 

 

 

 

As stated in the introduction to this chapter, P (t) is now being used to denote the pressure drop denoted by Δp in earlier chapters, and P0 is a constant representing the amplitude of the pressure wave and denoted by Δp0 in earlier chapters.

It is clear that in this case the relation between the pressure and flow waves is particularly simple. The flow wave has the same form and the same phase angle as the pressure wave, and the amplitude of the flow wave is given by the amplitude of the pressure wave divided by the resistance R, that is

phase

{Q(t)}

= phase {P (t)}

(6.2.4)

amplitude

{

Q(t)

}

=

amplitude {P (t)}

(6.2.5)

R

 

 

 

 

What is important is that this relationship between pressure and flow is fixed, it does not change within the oscillatory cycle. Equally important, however, this relation is only possible when the pressure wave is a simple sine or cosine function. This is clear from the way these solutions were obtained in Chapter 4.

Thus, the relation between pressure and flow cannot be applied to the composite waves of pressure and flow but it can be applied to their individual harmonics because, as we saw in Chapter 5, these harmonics consist of simple sine and cosine waves. Thus, we decompose the pressure wave into its harmonics by writing, as in Chapter 5,

P (t) =

p

. . .+ p1(t) + p2(t) + + pn(t)

(6.2.6)

where p1(t) . . . pn(t) are the n harmonics of the oscillatory part of p(t) and p is the (constant) average value of P (t) over one cycle, and if each of these parts of P (t) is now treated separately, we obtain the corresponding series of flow rates, using Eq. 6.2.3,

 

 

 

 

 

 

 

 

 

 

 

=

 

p

(6.2.7)

q

 

R

 

 

 

 

q1(t) =

p1(t)

(6.2.8)

 

 

R

 

 

 

 

 

 

 

q2(t) =

p2(t)

(6.2.9)

 

 

R

 

 

 

 

 

 

 

6.3 Example: Cardiac Pressure Wave

181

.

 

 

 

.

 

 

 

.

pn(t)

 

 

qn(t) =

(6.2.10)

R

 

 

 

These components can now be added to give the composite flow wave Q(t) produced by the composite pressure wave P (t), namely

Q(t) =

 

+ q1(t) + q2(t) + . . . + qn(t)

(6.2.11)

q

There are two important points to emphasize:

1. The sum in Eq. 6.2.11 to obtain the total Q(t) is only possible because the equation governing Q(t), namely Eq. 4.7.1, is linear.

2.The simple relation in Eq. 6.2.3 between pressure and flow applies only to the components of the pressure and flow waves but not to the composite waves themselves, that is

Q(t) =

P (t)

when P (t), Q(t) are composite waves (6.2.12)

R

6.3 Example: Cardiac Pressure Wave

As a first example of pressure-flow analysis, consider the cardiac waveform given numerically in Table 6.3.1 and shown in Fig. 6.3.1. The waveform represents cardiac pressure in mm Hg measured in a 20–kg dog. Note that in this case the pressure data include both the mean and the oscillatory parts of the pressure wave, and we use these below to illustrate the results of the previous section.

To obtain the corresponding flow wave using the scheme outlined in the previous section, the pressure waveform is decomposed into its harmonics as was done in Chapter 5, and the results for the first 10 harmonics are shown numerically in Table 6.3.2. Since the only opposition to flow in this case is the resistance R, then the corresponding harmonics of the flow wave are as in Eq. 6.2.10, namely

qn(t) =

pn(t)

n = 1, 2, . . . , 10

(6.3.1)

R

 

 

 

However, the resistance R is not known. Since in this section we are only interested in illustrating the analysis involved, we shall use an estimate of R for this purpose.

In the human cardiovascular system, for a 60–kg man, with a cardiac output of 5 L/min and a mean aortic pressure of 100 mm Hg, an estimate of total resistance in the systemic circulation would be

182 6 Composite Pressure-Flow Relations

 

120

 

 

 

 

 

 

115

 

 

 

 

 

(mmHg)

110

 

 

 

 

 

105

 

 

 

 

 

100

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

pressure

95

 

 

 

 

 

90

 

 

 

 

 

85

 

 

mean = 102.5995

 

 

 

 

 

 

80

 

 

 

 

 

 

75

 

 

 

 

 

 

0

0.2

0.4

0.6

0.8

1

 

 

 

time t (normalized)

 

 

Fig. 6.3.1. A cardiac pressure wave measured in a 20–kg dog and used as an example to illustrate how the corresponding flow wave is obtained from the given pressure waveform.

R (systemic, human)

100

[mm Hg]

 

 

 

 

(6.3.2)

 

 

5 [L/min]

 

 

 

 

 

 

 

 

mm Hg

 

 

 

 

 

=

20

 

 

 

 

 

 

 

 

(6.3.3)

L/min

 

 

 

 

100 × 1, 333

 

 

 

 

 

(6.3.4)

5

×

1, 000/60

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dynes

·

s

 

1, 600

cm5

 

 

(6.3.5)

 

 

If aortic pressure is taken as the input pressure driving flow into the coronary circulation, and if coronary blood flow is estimated at 5% of cardiac output, then an estimate of total resistance in the coronary circulation is

R (coronary, human)

100 [mm Hg]

 

(6.3.6)

0.25

[L/min]

 

 

 

 

 

mm Hg

 

 

= 400

 

 

 

(6.3.7)

L/min

100 × 1, 333

 

(6.3.8)

0.25

×

1, 000/60

 

 

 

6.3 Example: Cardiac Pressure Wave

183

 

 

dynes

·

s

 

32, 000

 

cm5

 

 

(6.3.9)

 

 

The corresponding estimates for a 20–kg dog, taking the mean cardiac output as 2 L/min, the mean aortic pressure as 100 mm Hg, and coronary blood flow again as 5% of cardiac output, we find

 

mm Hg

 

 

 

 

R (systemic, dog) 50

 

 

 

 

 

 

(6.3.10)

L/min

 

 

 

 

 

 

 

dynes

s

 

 

 

4, 000

 

 

·

 

 

(6.3.11)

 

cm5

 

R (coronary, dog) 1, 000

 

mm Hg

 

(6.3.12)

 

L/min

 

 

 

 

 

dynes

·

s

 

80, 000

 

cm5

 

 

(6.3.13)

 

 

Using the estimated value of R for the coronary system of the dog, and the harmonic components of the cardiac pressure wave in Table 6.3.2, where

Table 6.3.1. A numerical description of the cardiac wave shown in Fig. 6.3.1, giving the pressure (P ) at di erent times t within the oscillatory cycle. The oscillatory period has been normalized to 1.0. The pressure data include both the mean and the oscillatory part of the pressure.

t

P

t

P

0.000 96.60 0.500 111.00

0.025 96.21 0.525 108.15

0.050 97.27 0.550 103.65

0.075 95.56 0.575 103.00

0.100 95.34 0.600 103.00

0.125 95.38 0.625 103.00

0.150 93.46 0.650 102.00

0.175 91.92 0.675 101.53

0.200 93.88 0.700 101.00

0.225 100.04 0.725 100.26

0.250 104.50 0.750 101.00

0.275 106.68 0.775 101.00

0.300 108.20 0.800 101.00

0.325 110.00 0.825 102.00

0.350 110.95 0.850 102.00

0.375 112.38 0.875 101.13

0.400 113.80 0.900 100.70

0.425 113.00 0.925 99.86

0.450 113.00 0.950 99.47

0.475 112.93 0.975 98.18

184 6 Composite Pressure-Flow Relations

Table 6.3.2. Fourier coe cients of the first 10 harmonics of the pressure wave in Fig. 6.3.1.

 

n

An

 

 

 

 

Bn

 

 

 

Mn

φn (deg)

 

 

1

-6.5901

0.94298

6.65720

171.8568

 

 

2

1.08200

-4.66400

4.78780

-76.9394

 

 

3

0.74761

0.62007

0.97129

39.6723

 

 

4

0.15931

0.39243

0.42353

67.9050

 

 

5

-0.77719

0.93475

1.21560

129.7415

 

 

6

-0.28180

-0.56377

0.63028 -116.5580

 

 

7

-0.19808

-0.38691

0.43467 -117.1109

 

 

8

0.55090

-0.07722

0.55628

-7.9796

 

 

9

-0.23510

0.22975

0.32872

135.6601

 

 

10

-0.13825

0.18000

0.22696

127.5263

 

 

 

 

 

 

 

2πnt

 

 

 

 

 

 

 

 

pn(t) = Mn cos

 

 

T

− φn

n = 1, 2, . . . , 10

(6.3.14)

we then find

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

p

 

 

 

 

 

 

(6.3.15)

 

 

 

q

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

102.5995

 

 

 

 

(6.3.16)

 

 

1, 000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.1026

[L/min]

(6.3.17)

where p and q are the mean values of p(t) and q(t) over one oscillatory cycle. For the oscillatory parts of the flow wave, using Eq. 6.3.14 and Table 6.3.2, and taking the period T = 1, we find

q1(t) =

p1(t)

 

 

 

(6.3.18)

 

R

 

 

 

 

 

 

 

 

 

 

=

6.6572

× cos (2πt − 171.8568 × π/180)

 

(6.3.19)

 

 

 

 

 

1, 000

 

 

0.0067 × cos (2πt − 171.8568 × π/180)

[L/min]

(6.3.20)

q2(t) =

p2(t)

 

 

 

(6.3.21)

 

R

 

 

 

 

 

 

 

 

 

 

=

4.7878

× cos (4πt + 76.9394 × π/180)

 

(6.3.22)

 

 

 

 

 

1, 000

 

 

0.0048 × cos (4πt + 76.9394 × π/180)

[L/min]

(6.3.23)

.

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

.

p10(t)

 

 

 

 

q10(t) =

 

 

 

(6.3.24)

 

R

 

 

 

 

 

 

 

 

 

 

=

0.22696

 

 

 

 

 

× cos (20πt − 127.5263 × π/180)

(6.3.25)

1, 000

 

6.3 Example: Cardiac Pressure Wave

185

0.0002 × cos (20πt − 127.5263 × π/180) [L/min] (6.3.26)

These components of Q(t) can now be added to give the composite flow wave Q(t) produced by the composite pressure wave P (t), that is

Q(t) =

 

+ q1(t) + q2(t) + + q10(t)

(6.3.27)

q

which is shown graphically together with the pressure wave in Fig. 6.3.2.

P (mmHg)

120

 

 

 

 

 

110

 

 

 

 

 

100

 

 

 

 

 

pressure

 

 

 

 

 

90

 

 

mean = 102.5995

 

 

 

 

 

 

 

 

0.12

 

 

 

 

 

Q (L/min)

0.11

 

 

 

 

 

0.1

 

 

 

 

 

flow

0.09

 

 

mean = 0.1026

 

 

 

 

 

 

0

0.2

0.4

0.6

0.8

1

 

 

 

time t (normalized)

 

 

Fig. 6.3.2. Cardiac pressure wave (top) and corresponding flow wave (bottom) when opposition to flow consists of only resistance R which has been estimated at 1, 000 mm Hg · L/min. We see that in the presence of pure resistance the pressure and flow waves have precisely the same form, the flow wave being only scaled by the value of the resistance R.

It would be convenient to put the pressure and flow waves together, using the same scale, and the results of this section suggest that in order to do so we should plot P (t) and R × Q(t) (instead of Q(t)), as shown in Fig. 6.3.3. This presentation would be useful not only when the opposition to flow consists of pure resistance but also, and particularly, when other elements of the RLC system are present. In such cases, as we shall see, any small change in the form of the flow wave can be detected more easily and can be attributed directly to inertial (L) or capacitive (C) e ects only, because the e ects of resistance (R) have been scaled out. We shall refer to the product R × Q(t) as the “R-scaled” flow.

186 6 Composite Pressure-Flow Relations

 

120

 

 

 

 

 

(mmHg)

115

 

 

 

 

 

110

 

 

 

 

 

flow

 

 

 

 

 

105

 

 

 

 

 

R−scaled

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

pressure,

95

 

 

 

 

 

90

 

 

 

 

 

 

 

 

 

 

 

 

85

 

 

 

 

 

 

0

0.2

0.4

0.6

0.8

1

 

 

 

time t (normalized)

 

 

Fig. 6.3.3. When studying pressure-flow relations it is convenient to plot the pressure and flow waves to the same scale so as to compare their waveforms. This can be achieved as seen here by plotting P (t) (solid curve) and the “R-scaled” flow R × Q(t) (dashed) instead of Q(t). When the opposition to flow consists of only resistance R as it is in this case, the two curves become graphically identical. In this figure they are slightly shifted to make them visibly distinct. The use of R-scaled flow is useful also when other elements of the RLC system are present. In such cases any small change in the form of the flow wave can be detected more easily and can be attributed directly to inertial (L) or capacitive (C) e ects only, because the e ects of resistance (R) have been scaled out.

6.4Composite Pressure-Flow Relations Under General Impedance

When opposition to pulsatile flow consists of more than pure resistance, that is, when either inertial or capacitance e ects or both are involved, the relation between a composite flow wave and the corresponding flow wave is more complicated than was seen in the previous section. For the purpose of considering this relation here, the term “impedance” shall be used in a general sense here to mean any form of opposition to flow beyond that of pure resistance.

It is convenient in this and in subsequent sections to use the concept of complex impedance Z introduced in Section 4.9. It was seen in that section that when the opposition to flow consists of only pure resistance, the complex impedance Z becomes real and equal to R, but when any inertial or capacitive

6.4 Composite Pressure-Flow Relations Under General Impedance

187

e ects are present, that is any “reactive” e ects, Z is a complex quantity whose real and imaginary parts depend on the nature and arrangement of the reactive elements. In these terms, our interest in this section is in pressure-flow relations when Z has both a real and an imaginary part.

Let the composite pressure wave under consideration be denoted by P (t) and the corresponding flow wave be denoted by Q(t). In Chapter 5 we saw that P (t) can always be separated into a steady part p and a purely oscillatory part, such that

P (t) =

 

+ p(t)

(6.4.1)

p

The mean part of the pressure represents the mean value of P (t) over one oscillatory cycle, namely

 

1

0

T

 

p =

P (t)dt

(6.4.2)

T

The oscillation part of the pressure, by definition, has a zero mean. This separation of the composite wave is important because the corresponding flow wave will consist similarly of steady and oscillatory parts, which we shall denote by q and q(t), respectively, and write

Q(t) =

 

+ q(t)

(6.4.3)

q

The relation between the mean components of the pressure and flow waves, p, q, is fundamentally di erent from that between the oscillatory components, p(t), q(t), and the two relations must be dealt with separately. For the steady components, the relation between pressure and flow is simply that established in Section 6.2, namely

 

 

 

 

 

 

 

 

=

 

 

p

(6.4.4)

q

R

 

 

 

While in Section 6.2 the opposition to flow consisted of only pure resistance, this relation remains valid in this section even though reactive elements are assumed to be present. The reason for this is that reactive e ects come into play only when flow is non-steady. And since here we are able to separate the steady and non-steady parts of the flow, as discussed in Section 6.2, Eq. 6.4.4 can be used for the steady parts of the pressure and flow waves.

One of the most important advantages of using complex impedance is that the relation between the oscillatory parts of the pressure and flow waves can then be put in the general form

q(t) =

p(t)

(6.4.5)

Z

 

 

 

However, in this equation the pressure p(t) must be in complex form. And since both p(t) and Z are complex, it follows that q(t) is also complex. Thus, if subscripts r and i are used to denote real and imaginary parts, and we write

188

6 Composite Pressure-Flow Relations

 

 

 

 

p(t) = pr (t) + ipi(t)

 

 

Z = zr + izi

 

 

 

 

q(t) = qr (t) + iqi(t)

then Eq. 6.4.5 can be put in the form

 

 

 

 

qr (t) + iqi(t) =

pr (t) + ipi(t)

 

 

 

 

 

zr + izi

 

 

 

 

 

 

 

 

 

= pr (t)zr + pi(t)zi + i(pi(t)zr − pr (t)zi)

 

 

 

 

 

z2

+ z2

 

 

 

 

 

r

i

and we find

 

 

 

 

 

 

qr (t) =

pr (t)zr + pi(t)zi

 

 

z2

+ z2

 

 

 

 

 

 

r

 

i

 

qi(t) =

pi(t)zr − pr (t)zi

 

 

 

 

z2

+ z2

 

 

 

r

 

i

(6.4.6)

(6.4.7)

(6.4.8)

(6.4.9)

(6.4.10)

(6.4.11)

(6.4.12)

As we saw in earlier sections, the real part of the oscillatory flow rate, namely qr (t), represents the flow rate when the driving pressure is the real part of p(t), namely pr (t), and similarly for the imaginary parts of pressure and flow. But as we see clearly from Eqs. 6.4.11, 12

qr (t) =

pr (t)

 

 

 

 

zr

 

qi(t) =

pi(t)

 

 

 

zi

 

The correct pressure-flow relation is

 

 

 

 

 

qr (t) =

pZ

 

 

 

 

(t)

 

 

 

=

 

r zr + izi

 

 

 

p

(t) + ipi(t)

 

qi(t) =

pZ

 

 

 

 

 

(t)

 

 

 

=

 

rzr

+ iziZ

 

 

 

p

(t) + ipi(t)

 

(6.4.13)

(6.4.14)

(6.4.15)

(6.4.16)

(6.4.17)

(6.4.18)

which yield the results in Eqs. 6.4.11, 12.

It is important to recall that Eq. 6.4.5 and all the above results that follow from it apply only when the driving pressure, pr or pi, is a simple sine or cosine

6.4 Composite Pressure-Flow Relations Under General Impedance

189

wave. Thus Eq. 6.4.5 cannot be applied directly to a composite pressure wave, but it can be applied to its individual harmonics. If there are N harmonics and we denote these by pn(t), then, as we found in Chapter 5, they are given by

pn(t) = Mn cos

2πnt

− φn

n = 1, 2, . . . , N

(6.4.20)

T

and the oscillatory part of the pressure wave is given by

p(t) = p1(t) + p2(t) + . . . + pN (t)

(6.4.21)

where Mn and φn are (real) constants associated with the Fourier series representation of the composite wave.

Eq. 6.4.5 applies only to each of these harmonics individually, and only if each is considered to be the real or the imaginary part of the complex pressure. Thus, for the first harmonic we can introduce a complex pressure

p1(t) = M1 cos

T

− φ1

+ iM1 sin

 

2T − φ1

 

(6.4.22)

 

 

2πt

 

 

 

πt

 

 

and then use Eq. 6.4.5 to get the corresponding harmonic of the flow wave

 

q1

(t) =

p1(t)

 

 

 

(6.4.23)

 

Z

 

 

 

 

 

 

 

 

This can be done with each harmonic, writing, in general

− φn

 

pn(t) = Mn cos

T

− φn + iMn sin

2 T

(6.4.24)

 

2πnt

 

 

 

πnt

 

which can be put in the more compact exponential notation

 

pn(t) = Mnei((2πnt/T )φn ) n = 1, 2, . . . , N

(6.4.25)

The corresponding harmonics of the oscillatory flow wave are then given by

qn(t) =

pn(t)

n = 1, 2, . . . , N

(6.4.26)

Z

 

 

 

 

The oscillatory part of the composite pressure wave is now seen to be the real part of the complex pressure pn(t), which in turn corresponds to the real part of the flow wave, that is, in the notation of Eqs. 6.4.6, 8, where subscripts r and i are used for denoting real and imaginary parts, respectively,

pnr (t) = (pn(t))

(6.4.27)

= Mn cos

2πnt

− φn

n = 1, 2, . . . , N

(6.4.28)

T