Tutorials › Quantitative Systems Pharmacology (QSP) › Arterial Compliance and Two-Element Windkessel Dynamics

Home  /  Tutorials  /  Cardiovascular and Renal Physiology

Cardiovascular and Renal Physiology · Tutorial 363 of 799

Arterial Compliance and Two-Element Windkessel Dynamics

Learn how arterial compliance and peripheral resistance combine to determine pressure decay between heartbeats in the two-element Windkessel model.

Advanced 11 min read

What You'll Learn

  • Derive the pressure equation from arterial volume balance and a linear compliance assumption.
  • Explain why diastolic pressure above the outflow reference decays exponentially.
  • Calculate the Windkessel time constant and interpret its units.
  • Compare pressure decay when compliance or peripheral resistance changes.
  • Estimate the time constant from hypothetical pressure measurements.
  • Recognize the model’s assumptions and common diagnostic warning signs.

Why Arterial Pressure Persists Between Heartbeats

Blood enters the large arteries in pulses, but flow to the tissues does not stop completely between beats. The elastic arteries expand as blood enters and recoil as it leaves. This storage-and-release behavior helps sustain pressure during diastole. The two-element Windkessel model represents that behavior with an arterial compliance and a peripheral resistance.

The model is deliberately compact: it does not describe pressure at every location in the circulation. Instead, it approximates the arterial system as a compliant reservoir draining through a resistance to a downstream reference pressure. Its main dynamic prediction is that, when arterial inflow temporarily stops, pressure above that reference decays exponentially.

Definition: Arterial compliance \(C\) is the change in arterial volume per change in pressure, \(C=dV/dP\). In the two-element Windkessel model, \(C\) is constant, and a peripheral resistance \(R\) relates outflow to the pressure difference across it. The arterial pressure is \(P(t)\), the downstream reference pressure is \(P_{\mathrm{out}}\), arterial inflow is \(Q_{\mathrm{in}}(t)\), and outflow is \(Q_{\mathrm{out}}(t)\).

Here, pressure is measured in mmHg, volume in mL, and flow in mL/s. With these units, \(C\) has units mL/mmHg and \(R\) has units mmHg·s/mL. The model assumes that the arterial compartment can be represented by a single pressure, that compliance is linear and constant over the pressure range of interest, and that outflow is proportional to the pressure difference. It also treats \(P_{\mathrm{out}}\) as fixed over the interval being modeled.

Deriving the Pressure Equation

The starting point is conservation of volume: the rate of change of blood volume stored in the arterial compartment equals inflow minus outflow. The model’s compliance assumption connects that volume change to pressure change, while its resistance assumption specifies outflow.

$$ \frac{dV}{dt}=Q_{\mathrm{in}}-Q_{\mathrm{out}}, \qquad \frac{dV}{dt}=C\frac{dP}{dt}, \qquad Q_{\mathrm{out}}=\frac{P-P_{\mathrm{out}}}{R}. $$

Substituting the compliance and resistance relationships into the volume balance gives the two-element Windkessel pressure equation:

$$ C\frac{dP}{dt} = Q_{\mathrm{in}}(t)-\frac{P-P_{\mathrm{out}}}{R}. $$

This equation says that inflow raises stored arterial volume and pressure, while pressure-driven outflow lowers them. During a period of diastole when arterial inflow is approximated as zero, it becomes a simple decay equation. Define the pressure excess above the outflow reference as \(P-P_{\mathrm{out}}\). Then:

$$ \frac{d(P-P_{\mathrm{out}})}{dt} = -\frac{P-P_{\mathrm{out}}}{RC}. $$

If \(P_{\mathrm{out}}\) is constant, the solution from an initial pressure \(P_0\) at time \(t=0\) is:

$$ P(t)=P_{\mathrm{out}}+(P_0-P_{\mathrm{out}})e^{-t/(RC)}. $$
Key result: The Windkessel time constant is \(\tau=RC\). During zero inflow, the pressure excess \(P-P_{\mathrm{out}}\) falls to \(1/e\), or about 36.8%, of its initial value after one time constant. A larger \(R\), a larger \(C\), or both produce a slower decay.

The units confirm that \(\tau\) is a time: \((\mathrm{mmHg\cdot s/mL})(\mathrm{mL/mmHg})=\mathrm{s}\). The pressure in the exponential formula is the excess above \(P_{\mathrm{out}}\), not the absolute pressure. This distinction matters whenever the downstream reference is not zero.

Worked Example: Calculating Diastolic Pressure Decay

Suppose a hypothetical arterial compartment has \(C=1.6\ \mathrm{mL/mmHg}\), \(R=1.2\ \mathrm{mmHg\cdot s/mL}\), and \(P_{\mathrm{out}}=5\ \mathrm{mmHg}\). At the start of diastole, pressure is \(P_0=100\ \mathrm{mmHg}\). Estimate pressure one second later, assuming inflow is zero during that interval.

The zero-inflow model applies over the interval, and \(R\) and \(C\) are treated as constant. The time constant is:

$$ \tau=RC=(1.2)(1.6)=1.92\ \mathrm{s}. $$

The starting pressure excess is \(100-5=95\ \mathrm{mmHg}\). Substituting into the decay equation:

$$ P(1.00) = 5+95e^{-1.00/1.92} \approx 5+95(0.5940) \approx 61.4\ \mathrm{mmHg}. $$

Thus, under these hypothetical assumptions, pressure is about \(61.4\ \mathrm{mmHg}\) one second into diastole. It remains above \(P_{\mathrm{out}}\), as expected for an exponential decay that approaches the reference pressure rather than crossing it.

How Compliance and Resistance Shape the Decay

The product \(RC\) controls the rate of decay, but the two parameters represent different physiological properties. Resistance determines how readily blood leaves the arterial compartment for a given pressure difference. Compliance determines how much arterial volume change corresponds to a given pressure change. Increasing either one lengthens the time constant, all else equal.

A larger compliance also reduces the pressure change caused by a given volume addition. For a rapid inflow pulse that adds volume \(\Delta V\) before substantial outflow occurs, the linear compliance relation gives \(\Delta P=\Delta V/C\). Therefore, for the same added volume, greater compliance means a smaller pressure increment. This pulse approximation does not specify the full systolic waveform; it isolates the storage effect.

Worked Example: Comparing Compliance and Resistance Changes

Use the initial and downstream pressures from the prior example: \(P_0=100\ \mathrm{mmHg}\) and \(P_{\mathrm{out}}=5\ \mathrm{mmHg}\). Compare pressure after one second for (a) a reduction of compliance to \(0.8\ \mathrm{mL/mmHg}\), with \(R=1.2\ \mathrm{mmHg\cdot s/mL}\), and (b) an increase of resistance to \(2.4\ \mathrm{mmHg\cdot s/mL}\), with \(C=1.6\ \mathrm{mL/mmHg}\). Assume zero inflow and constant parameters in both cases.

(a) Lower compliance. The time constant is \(\tau=(1.2)(0.8)=0.96\ \mathrm{s}\). Therefore:

$$ P(1) = 5+95e^{-1/0.96} \approx 5+95(0.3529) \approx 38.5\ \mathrm{mmHg}. $$

(b) Higher resistance. The time constant is \(\tau=(2.4)(1.6)=3.84\ \mathrm{s}\). Therefore:

$$ P(1) = 5+95e^{-1/3.84} \approx 5+95(0.7707) \approx 78.2\ \mathrm{mmHg}. $$

Both changes alter the time constant, but in opposite directions: lower compliance shortens it and speeds the fall, whereas higher resistance lengthens it and slows the fall. These comparisons hold the other model inputs fixed; they are not, by themselves, claims about the net effect of a physiological intervention that changes several properties at once.

The model also gives a useful relation between average pressure and average flow in a periodic steady state. Integrate the pressure equation over a complete cycle of duration \(T\). If arterial pressure returns to its starting value at the end of the cycle, the integral of \(C\,dP/dt\) is zero. The average inflow must then equal the average outflow, and the model’s linear resistance relation gives:

$$ \overline{Q}_{\mathrm{in}} = \overline{Q}_{\mathrm{out}} = \frac{\overline{P}-P_{\mathrm{out}}}{R}, \qquad \overline{P}=P_{\mathrm{out}}+R\overline{Q}_{\mathrm{in}}. $$

This is a cycle-average consequence of volume balance, not a prediction that pressure is constant during the cycle. The instantaneous pressure still depends on the timing of inflow and on compliance.

Worked Example: Estimating the Time Constant From Decay Data

In a hypothetical measurement, arterial pressure is \(100.0\ \mathrm{mmHg}\) at the start of a zero-inflow interval and \(61.43\ \mathrm{mmHg}\) one second later. The downstream reference pressure is \(5.0\ \mathrm{mmHg}\). Estimate \(\tau\), then estimate \(R\) if independent information gives \(C=1.6\ \mathrm{mL/mmHg}\). Assume the single-compartment exponential model is appropriate over the measured interval.

Rearrange the decay solution using the ratio of pressure excesses:

$$ \frac{P(t)-P_{\mathrm{out}}}{P_0-P_{\mathrm{out}}} =e^{-t/\tau}, \qquad \tau = -\frac{t}{ \ln\left(\dfrac{P(t)-P_{\mathrm{out}}}{P_0-P_{\mathrm{out}}}\right) }. $$

At one second, the ratio is \((61.43-5.0)/(100.0-5.0)=56.43/95.0\approx0.5940\). Thus:

$$ \tau = -\frac{1.00\ \mathrm{s}}{\ln(0.5940)} \approx 1.920\ \mathrm{s}. $$

Since \(\tau=RC\), the resistance estimate is:

$$ R=\frac{\tau}{C} = \frac{1.920\ \mathrm{s}}{1.6\ \mathrm{mL/mmHg}} = 1.20\ \mathrm{mmHg\cdot s/mL}. $$

The decay data identify the product \(RC\); estimating \(R\) separately required an independently specified value of \(C\). Without that information, many resistance-compliance pairs could have the same time constant.

Assumptions, Diagnostics, and Alternatives

The exponential formula is useful only when its assumptions are reasonable over the interval being analyzed. In particular, the interval should have negligible arterial inflow, a stable downstream reference pressure, and approximately constant resistance and compliance. If inflow continues during the interval, the zero-inflow decay formula does not apply. For constant inflow \(Q_{\mathrm{in}}\), the pressure instead approaches \(P_{\mathrm{out}}+RQ_{\mathrm{in}}\), with the same time constant:

$$ P(t) = P_{\mathrm{out}}+RQ_{\mathrm{in}} + \left[P(0)-P_{\mathrm{out}}-RQ_{\mathrm{in}}\right]e^{-t/(RC)}. $$

A useful diagnostic is to plot the logarithm of pressure excess against time during a candidate zero-inflow interval. Under the model, \(\ln(P-P_{\mathrm{out}})\) is a straight line with slope \(-1/(RC)\). Systematic curvature, a changing estimated slope, or an apparent pressure excess that becomes nonpositive can signal that the interval or model assumptions are unsuitable. Measurement noise can also make the transformed plot unstable, especially when pressure is close to the reference.

The two-element model is an approximation, not a detailed account of wave propagation, spatially varying arterial properties, or pressure-dependent compliance. A purely resistive description can relate average pressure difference to average flow but cannot represent storage or diastolic decay. A model with an additional proximal resistance can distinguish more features of pulsatile input and pressure, but adds parameters and is not needed to derive the two-element model’s basic decay relation. Model choice should match the question and the information available.

Common pitfall: A learner may say that increasing compliance makes pressure decay faster because the arteries “stretch more.” That reverses the model’s prediction when resistance is fixed: increasing \(C\) increases \(RC\), so the pressure excess decays more slowly. Compliance can reduce the pressure rise from a fixed incoming volume, but the subsequent decay rate depends on the product \(RC\).

Check Your Understanding

Answer each question, then compare your reasoning with the solutions below.

  1. A hypothetical system has \(R=1.5\ \mathrm{mmHg\cdot s/mL}\) and \(C=2.0\ \mathrm{mL/mmHg}\). What is its time constant, and what fraction of the initial pressure excess remains after one time constant?
  2. Pressure falls from \(90\ \mathrm{mmHg}\) toward a fixed reference of \(10\ \mathrm{mmHg}\) during zero inflow. What pressure does the model predict after one time constant?
  3. If resistance is unchanged and compliance doubles, what happens to the time constant and to the rate of diastolic decay?
  4. Why can a single measured decay time constant not determine both \(R\) and \(C\) without additional information?

Solutions. (1) \(\tau=RC=(1.5)(2.0)=3.0\ \mathrm{s}\). After one time constant, the excess is \(e^{-1}\approx0.368\), or 36.8%, of its initial value. (2) The initial excess is \(90-10=80\ \mathrm{mmHg}\), so after one time constant \(P=10+80e^{-1}\approx39.4\ \mathrm{mmHg}\). (3) Doubling \(C\) doubles \(RC\), so the time constant doubles and the excess pressure decays more slowly. (4) The decay identifies only the product \(RC\). Multiple pairs with the same product produce the same exponential time constant, so one parameter must be known or independently estimated to determine the other.

Previous

Next

Back to Tutorials
Clinical Biostats