Cardiovascular physiology — the heart

Intuition [Beginner]

The heart is a muscular pump that drives blood through two circuits. The pulmonary circuit sends deoxygenated blood to the lungs; the systemic circuit sends oxygenated blood to the rest of the body. Each circuit is served by one side of the heart: the right side pumps to the lungs, the left side pumps to the body.

A single heartbeat — one cardiac cycle — consists of two main phases. During systole, the heart muscle contracts and ejects blood. During diastole, the heart relaxes and fills with blood. This cycle repeats roughly 70 times per minute at rest, producing the pulse you feel at your wrist.

The amount of blood the heart pumps each minute is called cardiac output. It depends on two things: how fast the heart beats (heart rate) and how much blood it squeezes out per beat (stroke volume). A typical resting adult has a heart rate of 72 beats per minute and a stroke volume of 70 mL, giving a cardiac output of about 5 litres per minute. That is roughly the entire blood volume of an average adult recycled once every minute.

Visual [Beginner]

The cardiac cycle can be visualised as a loop on a pressure-volume diagram. The horizontal axis shows the volume of the left ventricle. The vertical axis shows the pressure inside it.

Starting at the bottom-right corner (low pressure, high volume — end of filling), the heart contracts isovolumically (pressure rises, volume stays constant), then ejects blood through the aortic valve (pressure stays high, volume drops), then relaxes isovolumically (pressure drops, volume stays constant), then fills again (volume rises, pressure stays low). The enclosed area on the diagram equals the mechanical work done by the heart each beat.

The electrocardiogram (ECG) records the electrical activity that triggers each contraction. Three main deflections mark the cycle: the P wave (atrial depolarisation), the QRS complex (ventricular depolarisation), and the T wave (ventricular repolarisation). Each wave corresponds to a specific mechanical event.

Worked example [Beginner]

Calculate cardiac output for a resting adult.

Given: heart rate = 72 beats/min, stroke volume = 70 mL/beat.

Step 1. Multiply heart rate by stroke volume:

$$\text{CO}=\text{HR}\times \text{SV}=72\times 70=5040\text{ mL/min}.$$

Step 2. Convert to litres per minute:

$$5040\text{ mL/min}\approx 5.0\text{ L/min}.$$

This value — about 5 litres per minute — matches the total blood volume of a typical 70-kg adult. During intense exercise, cardiac output can rise to 20-25 L/min through increases in both heart rate (up to 180-200 bpm) and stroke volume (up to 120-150 mL).

Check your understanding [Beginner]

Formal definition [Intermediate+]

The cardiac cycle is a periodic process with period $T=1/\text{HR}$ (where HR is in beats per second). Each cycle consists of four phases applied to each ventricle:

1. Isovolumic contraction. All four heart valves are closed. Ventricular pressure rises from its end-diastolic value to the aortic (or pulmonary) pressure. Volume is constant at the end-diastolic volume (EDV).

2. Ventricular ejection (systole). The semilunar valve (aortic or pulmonary) opens. Ventricular pressure exceeds arterial pressure, driving blood out. The volume drops from EDV toward the end-systolic volume (ESV). The stroke volume is:

$$\text{SV}=\text{EDV}-\text{ESV}.$$

3. Isovolumic relaxation. All valves close again. Ventricular pressure drops. Volume stays at ESV.

4. Ventricular filling (diastole). The atrioventricular valve (mitral or tricuspid) opens. Blood flows passively from the atrium into the ventricle, then an atrial contraction (atrial kick) tops off the volume to EDV.

Cardiac output is the volume of blood pumped by one ventricle per unit time:

$$\text{CO}=\text{HR}\times \text{SV}=\text{HR}\times (\text{EDV}-\text{ESV}).$$

Ejection fraction quantifies contractile efficiency:

$$\text{EF}=\frac{\text{SV}}{\text{EDV}}=\frac{\text{EDV}-\text{ESV}}{\text{EDV}}.$$

Normal ejection fraction is 55-70%. Values below 40% indicate systolic heart failure.

Mean arterial pressure (MAP) relates cardiac output to systemic vascular resistance (SVR):

$$\text{MAP}=\text{CO}\times \text{SVR}+\text{CVP},$$

where CVP is central venous pressure (approximately 0 mmHg at the level of the heart, often neglected). In practical terms, MAP is estimated from cuff measurements:

$$\text{MAP}\approx \text{diastolic}+\frac{1}{3}(\text{systolic}-\text{diastolic}).$$

The Frank-Starling law

The Frank-Starling mechanism states that, within physiological limits, the stroke volume of the ventricle increases in proportion to its end-diastolic volume. That is, the heart pumps whatever volume of blood returns to it. The mechanism arises because greater stretch of cardiac myocytes increases the sensitivity of myofilaments to calcium, producing more forceful contraction.

Quantitatively, the relationship between ventricular end-diastolic volume and stroke volume (or stroke work) is plotted as a ventricular function curve. Increased sympathetic tone or inotropic agents shift this curve upward (greater SV at the same EDV); decreased contractility (heart failure) shifts it downward.

Haemodynamics

Blood flow through a vessel is governed by a fluid-mechanics analogue of Ohm's law:

$$Q=\frac{\Delta P}{R},$$

where $Q$ is flow rate, $\Delta P$ is the pressure difference across the vessel segment, and $R$ is vascular resistance. For a cylindrical vessel of radius $r$ and length $L$ carrying a Newtonian fluid of viscosity $\eta$, the Poiseuille equation gives:

$$R=\frac{8\eta L}{\pi r^4}.$$

The $r^4$ dependence means that small changes in vessel radius produce large changes in resistance. Halving the radius increases resistance 16-fold. This is the primary mechanism by which arteriolar vasoconstriction and vasodilation regulate blood pressure and flow distribution.

Key theorem with proof [Intermediate+]

Theorem (Frank-Starling equilibrium). Under steady-state conditions, the left and right ventricles must have equal cardiac outputs over time. If transiently the right ventricular output exceeds the left, pulmonary blood volume rises, increasing left ventricular preload, which — via the Frank-Starling mechanism — increases left ventricular output until equilibrium is restored.

Proof. Let $Q_L$ and $Q_R$ denote the cardiac outputs of the left and right ventricles respectively. The total blood volume $V_{\text{total}}$ is partitioned among four compartments: systemic veins ($V_s$), systemic arteries, pulmonary veins ($V_p$), and pulmonary arteries. Conservation of blood volume requires:

$$\frac{d{V}_s}{dt}=Q_L-Q_R.$$

If $Q_R>Q_L$, then $d{V}_s/dt<0$ and systemic venous volume falls, while pulmonary venous volume rises by the same amount. The increase in pulmonary venous volume increases left atrial pressure, which increases left ventricular end-diastolic volume (preload). By the Frank-Starling mechanism, increased preload increases $Q_L$. Simultaneously, decreased systemic venous return reduces right ventricular preload, reducing $Q_R$. The system converges to $Q_L=Q_R$. The converse argument applies if $Q_L>Q_R$. The equilibrium is stable because the Frank-Starling curve is monotonically increasing. $\square$

This theorem explains why cardiac output is self-regulating: the two sides of the heart automatically match their outputs without neural intervention. Clinical corollary: left ventricular failure causes pulmonary congestion (blood accumulates upstream of the failing pump) rather than systemic hypotension as the immediate problem.

Blood pressure regulation

Short-term blood pressure regulation occurs via the baroreceptor reflex. Stretch receptors in the carotid sinus and aortic arch sense arterial pressure. Increased pressure increases baroreceptor firing rate, which signals the medullary cardiovascular centre to increase parasympathetic (vagal) tone and decrease sympathetic tone. This reduces heart rate, contractility, and vascular resistance, lowering blood pressure back toward the set point. The reflex operates on a timescale of seconds to minutes.

Long-term regulation involves the kidneys and the renin-angiotensin-aldosterone system (RAAS), covered in unit 18.08.01 pending.

Exercises [Intermediate+]

Cardiac mechanics — Frank-Starling, length-tension, and the pressure-volume loop [Master]

The mechanical behaviour of the ventricle as a pump is encoded in two coupled constitutive laws — the end-systolic pressure-volume relationship (ESPVR) and the end-diastolic pressure-volume relationship (EDPVR) — together with a time-varying elastance function $E(t)$ that interpolates between them across the cardiac cycle. The whole apparatus is the systems-physiology lift of the actin-myosin cross-bridge cycle of 18.04.02 pending to a chambered, geometrically organised pump.

The Frank-Starling law in its modern form is an empirical curve on the (EDV, SV) plane. The ESPVR provides its mechanistic foundation: it is the locus of end-systolic pressure-volume points obtained by varying preload at constant inotropy. Over the physiological range, the ESPVR is approximately linear,

$$P_{\mathrm{es}}(V)=E_{\mathrm{es}}\cdot (V-V_0),$$

where $E_{\mathrm{es}}$ is the end-systolic elastance (a contractility index, $mmHg/mL$) and $V_0$ is the volume-axis intercept (the unstressed volume of the ventricle at zero pressure). The EDPVR is approximately exponential, $P_{\mathrm{ed}}(V)=\alpha (e^{\beta (V-V_0)}-1)$, reflecting the nonlinear stiffness of the relaxed myocardium and the pericardial constraint at high volumes. Each pressure-volume loop is bounded above by the ESPVR (the upper-left corner is the end-systolic point) and below by the EDPVR (the lower-right corner is the end-diastolic point).

The mechanistic origin of the Frank-Starling law lies in length-dependent activation of the sarcomere. At longer sarcomere length (within the physiological range of $1.9-2.3\mu \mathrm{m}$ for ventricular myocardium), three effects combine. First, the filament-overlap length-tension relation of 18.04.02 pending places the sarcomere on the ascending limb of the Gordon-Huxley-Julian curve, so more myosin heads can engage actin. Second, myofilament calcium sensitivity increases with stretch: at longer sarcomere length, the same intracellular $[{\mathrm{Ca}}^{2+}]$ transient produces a higher fraction of activated cross-bridges. Third, the lattice spacing between thick and thin filaments shrinks at longer sarcomere length (the myocyte preserves volume as it lengthens), bringing myosin heads closer to actin and increasing attachment probability. The combined effect is a steep stretch-dependent rise in force at a given calcium transient — the cellular substrate of the chamber-level Frank-Starling law.

A complementary chamber-level description uses the time-varying elastance model of Suga and Sagawa (1973). The model treats the ventricle as a one-port elastic chamber whose elastance $E(t)=P(t)/(V(t)-V_0)$ rises during systole from $E_{\min }$ (end-diastolic stiffness) to $E_{\max }=E_{\mathrm{es}}$ (end-systolic elastance) and returns to $E_{\min }$ during diastole. The time course of $E(t)$ is approximately load-independent across a range of physiological loading conditions — a remarkable empirical regularity that elevates $E(t)$ from a fitting parameter to a state function of the myocardium itself. Inotropic interventions (sympathetic activation, calcium sensitiser drugs) shift $E_{\max }$ upward; preload and afterload changes leave $E(t)$ essentially unchanged but move the operating PV loop within the bounded region.

The Hill force-velocity relation of 18.04.02 pending is the muscle-fibre constitutive law underlying the chamber-level PV loop. For an isotonic afterloaded contraction, the rate of shortening $v$ of the wall at any instant is determined by the instantaneous wall tension $T$ through the Hill curve $(T+a)(v+b)=(T_0+a)b$, with $T_0$ the isometric tension at the current sarcomere length. The chamber-level analogue is the dependence of ejection velocity $dV/dt$ on the difference between ventricular and aortic pressure during ejection. At zero afterload (ideal limit, not physiologically achievable), the ventricle would eject at $v_{\max }$ in a single rapid stroke; at infinite afterload (aortic valve refusing to open), the contraction would be purely isovolumic and the entire chamber-level force would equal the isometric chamber force at end-diastolic volume.

Stroke work is the area enclosed by the PV loop:

$$W_{\mathrm{stroke}}=\oint PdV.$$

A normal resting ventricle does $\sim 1\mathrm{J}$ per beat of mechanical work, giving a cardiac mechanical-power output of $\sim 1-1.5\mathrm{W}$ at rest. Pressure-volume area (PVA) — the sum of stroke work and the triangular "potential energy" region between the ESPVR and the isovolumic relaxation line — correlates linearly with myocardial oxygen consumption $V{O}_2$ per beat (Suga 1979). The slope of $V{O}_2$ vs PVA is independent of preload, afterload, and heart rate, varying only with contractile state and basal metabolism. This PVA-$V{O}_2$ relation is the operational definition of cardiac mechanical efficiency, $\eta =W_{\mathrm{stroke}}/(\Delta H_{{\mathrm{O}}_2}\cdot V{O}_2)$, which in resting myocardium runs around $20-25\%$ — comparable to skeletal muscle.

The link to molecular cross-bridge mechanics is direct. The chamber elastance $E(t)$ is a population-averaged moment of the cross-bridge distribution: more attached cross-bridges per unit cross-section means higher wall tension at given length, hence higher elastance. The Suga-Sagawa load-independence of $E(t)$ corresponds, at the molecular level, to the relatively slow timescale of sarcomere length change compared to cross-bridge cycle time — the ensemble equilibrates fast enough that the population-averaged force tracks the instantaneous geometry. The Frank-Starling steepness is set by the calcium-sensitivity slope of the cross-bridge activation curve.

A concrete edge case clarifies the framework. In systolic heart failure (dilated cardiomyopathy), $E_{\mathrm{es}}$ falls — the ESPVR rotates clockwise (less steep slope) and shifts rightward (larger $V_0$). The same end-diastolic volume yields a smaller stroke volume, so ejection fraction falls below the $\sim 55\%$ normal threshold. Compensatory ventricular dilation increases EDV further, partially recovering stroke volume at the cost of operating at higher wall stress (Laplace's law $T=Pr/(2h)$ for a thin-walled sphere — wall stress rises with chamber radius). The chronic remodeling that follows alters the EDPVR (less compliant chamber, steeper passive curve at lower volumes) and the inotropic state ($E_{\mathrm{es}}$ continues to fall as remodeling progresses). In diastolic heart failure (heart failure with preserved ejection fraction), $E_{\mathrm{es}}$ is approximately normal but the EDPVR is steeper — the ventricle cannot fill adequately at normal venous pressures. The PV-loop framework partitions cardiac dysfunction cleanly into these two categories, and each category has distinct pharmacological targets.

Vascular hemodynamics — Poiseuille, Windkessel, pulse waves, and Bernoulli at valves [Master]

The arterial tree is a branching network of compliant elastic tubes through which pulsatile flow propagates. Three classical results — Poiseuille resistance, the Windkessel model, and the Moens-Korteweg pulse-wave velocity — together with Bernoulli's principle at valves and stenoses, encode most of the relevant fluid mechanics.

Poiseuille flow. For steady laminar flow of a Newtonian fluid of viscosity $\eta$ through a long cylindrical tube of radius $r$ and length $L$, the volumetric flow rate is

$$Q=\frac{\pi r^4}{8\eta L}\Delta P,$$

which rearranges to a Hagen-Poiseuille resistance $R=8\eta L/(\pi r^4)$. The $r^4$ dependence is the central fact of vascular control: a $20\%$ reduction in radius approximately doubles the resistance, and the systemic peripheral resistance is dominated by the arterioles (radius $\sim 30-50\mu \mathrm{m}$) precisely because the calibre at that level is the most actively regulated. Resistance contributions from large arteries (radius $\sim 1\mathrm{cm}$) are negligible by comparison even though their length is much greater — the $r^4$ dominates the $L$.

Poiseuille flow is the asymptotic small-Reynolds-number limit of Navier-Stokes in a straight pipe; it requires (i) Newtonian fluid behaviour, (ii) no slip at the wall, (iii) fully developed parabolic velocity profile (no entrance effects), (iv) steady flow, (v) rigid pipe. None of these holds exactly in vivo: blood is a non-Newtonian suspension whose effective viscosity decreases at high shear (shear-thinning) and increases dramatically below a critical shear rate (Fåhraeus-Lindqvist effect changes effective viscosity in microvessels by a factor of two or more); flow is pulsatile, so the Womersley profile (a complex generalisation of the parabolic profile) is the correct steady-state for sinusoidal driving; and vessel walls are compliant, not rigid. Despite these caveats, the Poiseuille law remains the operational definition of vascular resistance in clinical hemodynamics and the right first-order estimate in most situations.

Two-element Windkessel model. The arterial tree is treated as a single compliant reservoir of compliance $C$ in parallel with a peripheral resistance $R$. Aortic pressure $P(t)$ obeys

$$C\frac{dP}{dt}+\frac{P-P_{\mathrm{venous}}}{R}=Q_{\mathrm{in}}(t),$$

where $Q_{\mathrm{in}}(t)$ is the pulsatile flow from the ventricle. The model has two distinct phases per cycle. During ejection, $Q_{\mathrm{in}}>0$ and the compliance absorbs roughly half the stroke volume while peripheral resistance bleeds off the other half; pressure rises to a systolic peak. During diastole, $Q_{\mathrm{in}}=0$, the compliance discharges through the resistance, and pressure decays exponentially $P(t)=P_{\mathrm{es}}{e}^{-t/RC}$ with time constant $\tau =RC$. Typical resting values: $R\approx 1,000\mathrm{dyn}\cdot s/c{\mathrm{m}}^5$, $C\approx 1,500\mathrm{c}{\mathrm{m}}^3/dyn\cdot 10^{-6}$, giving $\tau \sim 1.5\mathrm{s}$ — comfortably longer than diastolic duration ($\sim 0.5\mathrm{s}$ at rest), so pressure barely decays to its diastolic value before the next ejection.

The diastolic decay $\tau =RC$ is exploited clinically: $\tau$ can be fit from the diastolic limb of the aortic pressure waveform, $R$ is computed from $\mathrm{MAP}/\mathrm{CO}$, and $C$ is then estimated as $\tau /R$. Arterial stiffening with age and atherosclerosis reduces $C$, shortening $\tau$, widening pulse pressure (the systolic-diastolic difference), and increasing the cardiac afterload at the same mean pressure. The widening of pulse pressure with age is one of the most reproducible findings in cardiovascular epidemiology and is mechanistically captured by the Windkessel framework.

The three-element Windkessel of Westerhof (1969) adds a characteristic impedance $Z_0$ in series with the input, where $Z_0$ matches the acoustic impedance of the proximal aorta. This addition reduces the unphysical infinite high-frequency input impedance of the two-element model and captures the early-systolic pressure-flow phase relationship. The four-element Windkessel further adds an inductive element $L$ to represent blood inertia, improving fit at higher frequencies. These models trade anatomical realism for tractability; for whole-system simulation, transmission-line models that integrate the Navier-Stokes equations along a one-dimensional approximation of the vascular tree are used (e.g., Sherwin-Parker-Olufsen frameworks).

Pulse wave velocity and Moens-Korteweg. Because the arterial wall is elastic, a pressure perturbation propagates as a wave rather than instantaneously. The Moens-Korteweg equation for the wave speed in a thin-walled elastic tube of wall thickness $h$, lumen radius $r$, wall elastic modulus $E$, and blood density $\rho$ is

$$c=\sqrt{\frac{Eh}{2r\rho }}.$$

For young healthy human aorta, $c\approx 5m/s$; in elderly or atherosclerotic aorta, $c$ can reach $15m/s$ or higher. The increase is dominated by the rise in effective modulus $E$ as collagen replaces elastin in the wall and as wall thickening alters the geometry. Clinically, carotid-femoral pulse wave velocity (cfPWV) is measured by timing the foot of the pulse waveform between two arterial sites a known distance apart; cfPWV is now one of the standard non-invasive indices of arterial stiffness and an independent predictor of cardiovascular mortality.

The Moens-Korteweg equation derives from a linearised one-dimensional fluid-structure interaction: conservation of mass and momentum in the lumen coupled to a linear wall law $A(P)=A_0(1+(P-P_0)/(Eh/(2r)))$ where $A$ is cross-sectional area. The hyperbolic system for $(A,u)$ admits travelling-wave solutions whose speed is given by the formula above. The same derivation, applied to a stiffer compressible-tube model, recovers the water-hammer formula that pulse pressure $\Delta P=\rho c\Delta u$ where $\Delta u$ is the velocity perturbation — this is the relationship that lets a clinician estimate stroke volume from the foot-to-peak velocity change measured by Doppler.

Bernoulli at valves. In the valve plane of the aortic or pulmonary outlet, blood transitions from the high-volume low-velocity ventricular cavity to a much narrower valve orifice. The flow accelerates as it converges; by mass conservation, $u_{\mathrm{valve}}=u_{\mathrm{cavity}}\cdot A_{\mathrm{cavity}}/A_{\mathrm{valve}}$. Bernoulli's principle (steady, inviscid, incompressible, along a streamline) gives

$$\frac{1}{2}\rho u_{\mathrm{cavity}}^2+P_{\mathrm{cavity}}=\frac{1}{2}\rho u_{\mathrm{valve}}^2+P_{\mathrm{valve}}.$$

Neglecting the cavity velocity (much smaller than valve velocity at normal stroke volume), the pressure gradient across the valve is

$$\Delta P=P_{\mathrm{cavity}}-P_{\mathrm{valve}}\approx \frac{1}{2}\rho u_{\mathrm{valve}}^2.$$

Substituting $\rho \approx 1.06\cdot 10^{3}kg/{\mathrm{m}}^3$ for blood and converting to clinical units yields the famous simplified Bernoulli equation of echocardiography:

$$\Delta P(\mathrm{mmHg})\approx 4u_{\mathrm{valve}}^2(m/s{)}^2.$$

A peak aortic-jet velocity of $u=4m/s$ measured by continuous-wave Doppler implies a peak transvalvular gradient of $\sim 64\mathrm{mmHg}$, indicating severe aortic stenosis. The same calculation applied to mitral regurgitation, tricuspid regurgitation (allowing non-invasive estimation of pulmonary artery pressure), and septal-defect shunts has made the simplified Bernoulli equation perhaps the single most-used quantitative result in clinical cardiology.

The full Bernoulli equation includes a viscous-loss term and a convective-acceleration term, both of which are negligible for short valvular jets but become important for tubular stenoses and across long restrictions. At very high Reynolds number the flow downstream of the valve becomes turbulent, dissipates the kinetic energy of the jet as heat, and recovers only a fraction of the dynamic pressure — a separate clinical phenomenon (pressure recovery) that complicates the interpretation of Doppler gradients in some anatomical configurations.

Cardiac excitation-contraction coupling — calcium-induced calcium release [Master]

Cardiac muscle inherits the actin-myosin cross-bridge cycle of 18.04.02 pending but couples it to the cardiac action potential of 18.02.02 through a distinctive amplifying mechanism: calcium-induced calcium release (CICR), established as a quantitative mechanism by Fabiato (1983, 1985). The CICR step is the chemical bridge between the membrane-electrical event of the action potential and the mechanical event of contraction, and it sets cardiac contractility apart from the mechanical-coupling architecture of skeletal muscle.

The cardiac action potential of a ventricular myocyte is a long-plateau waveform of duration $\sim 200-300\mathrm{ms}$. The phase-2 plateau is sustained by a steady L-type calcium current ($I_{\mathrm{Ca},\mathrm{L}}$) flowing through Ca${}_v$1.2 channels (the cardiac dihydropyridine receptor, DHPR). The trans-sarcolemmal calcium influx during the plateau is modest, raising intracellular calcium concentration by perhaps $0.1-0.2\mu \mathrm{M}$ on its own — insufficient to activate full contraction. The amplifying step is the ryanodine receptor RyR2 on the junctional sarcoplasmic reticulum (jSR), which sits in the dyadic cleft within $\sim 15\mathrm{nm}$ of the DHPR cluster on the T-tubular membrane. The localised calcium signal in the cleft (a "calcium spark" if from one RyR2 cluster, or a coordinated transient if from many clusters firing in synchrony) triggers RyR2 opening through calcium-binding to a cytoplasmic activation site. Each RyR2 cluster then dumps the much larger SR-luminal calcium store ($\sim 1\mathrm{mM}$ free Ca${}^{2+}$ within the SR, $\sim 10,000$-fold over cytoplasmic) into the cleft and surrounding cytoplasm.

The amplification factor (RyR2-released calcium / trigger calcium) is $\sim 10$ under physiological conditions. Activation is graded — the amplitude of the calcium transient scales smoothly with the trigger calcium, rather than being all-or-none — a property crucial for the heart's ability to modulate contractility beat-to-beat. The graded character is itself an emergent property of the dyadic-cleft architecture: each RyR2 cluster fires stochastically with probability dependent on local trigger calcium, and the global transient is the sum over thousands of independent clusters per cell. The stochastic-but-averaged nature is the cardiac analogue of the law-of-large-numbers smoothing that gives skeletal muscle its smooth force despite molecular discreteness.

The released calcium binds cardiac troponin C (cTnC) at the EF-hand calcium-binding sites on the thin filament, with affinity tuned for the physiological calcium transient range. The cTnC isoform has lower calcium sensitivity than skeletal TnC and a single regulatory calcium-binding site (compared to two in skeletal TnC), giving cardiac myofilaments a less steep activation curve and a wider operating range. Calcium-bound cTnC pulls the inhibitory subunit cTnI off its actin contacts, the troponin-tropomyosin complex rotates azimuthally on the actin filament, and the myosin binding sites become available. The cross-bridge cycle of 18.04.02 pending then proceeds — phosphate release coupled to lever-arm rotation, ADP release, ATP-driven detachment — at rates set by cardiac myosin heavy-chain isoforms ($\beta$-MHC dominant in adult human ventricle, with slower ATPase and longer attached-state lifetime than the $\alpha$-MHC predominant in small mammals).

Calcium removal terminates the transient. Three pathways act in parallel: SERCA2a (the sarco/endoplasmic-reticulum calcium-ATPase) pumps calcium back into the SR (recycles it for the next beat), the sodium-calcium exchanger NCX1 extrudes one calcium for three sodium across the sarcolemma (a net loss to be balanced by next-beat L-type influx), and mitochondrial calcium uptake through the MCU contributes a small steady-state component. The balance between SERCA2a recycling and NCX1 extrusion is dynamically regulated and is a key determinant of heart-failure pathophysiology — SERCA2a expression falls in failing myocardium, slowing relaxation (lusitropic failure) and depleting the SR store between beats.

Three layers of regulation tune the CICR mechanism. Phospholamban (PLB) is a transmembrane SR protein that inhibits SERCA2a when dephosphorylated; phosphorylation of PLB by protein kinase A (downstream of ${\beta }_1$-adrenergic activation) relieves the inhibition, accelerates calcium uptake into the SR, and produces the positive lusitropic effect of sympathetic stimulation. Cardiac troponin I (cTnI) phosphorylation at Ser23/24 (also PKA-mediated) decreases myofilament calcium sensitivity, producing faster relaxation by accelerating calcium dissociation from cTnC. CaMKII phosphorylation of RyR2 alters its open probability and is implicated in the calcium-handling abnormalities of heart failure and arrhythmogenesis. Together these regulatory inputs let the cardiac cycle accelerate and amplify under sympathetic drive — an essential adaptation to exercise.

The cardiac cycle's integration of CICR with the cross-bridge cycle bridges three levels of organisation. At the molecular level, the cycle is a Markov chain on myosin chemo-mechanical states driven by ATP hydrolysis (as in 18.04.02 pending). At the cellular level, the cycle is a periodic calcium transient coupled to a periodic actin-myosin activation envelope, with the period set by the spontaneous depolarisation of pacemaker cells (sinoatrial node, treated in 18.02.02). At the chamber level, the cycle is the time-varying elastance $E(t)$ of the Suga-Sagawa framework, with $E_{\max }$ tracking the peak of the calcium transient and the relaxation phase tracking the SERCA-mediated calcium clearance. CICR is the chemical step that knits these three levels together: it transforms a brief membrane-voltage perturbation (the action-potential plateau) into a graded, amplified, recycled calcium transient that drives the cross-bridge population to its activated state.

The system-level bridge to 18.02.02 is that the L-type calcium current — the trigger for CICR — is shaped by the kinetics of the cardiac action potential plateau, which in turn depends on the balance between depolarising L-type calcium current and repolarising potassium currents ($I_{\mathrm{Kr}}$, $I_{\mathrm{Ks}}$, $I_{K1}$). Pharmacological agents that prolong the action-potential plateau (class-III antiarrhythmics) prolong the calcium transient, but also risk early afterdepolarisations and torsades-de-pointes arrhythmia; agents that block L-type calcium channels (verapamil, diltiazem) attenuate both the contractile transient and the chronotropic effect of sympathetic stimulation. The pharmacology of cardiovascular drugs is largely the pharmacology of the CICR coupling itself.

Cardiovascular regulation — autonomic control, baroreflex, RAAS, and exercise [Master]

The cardiovascular system maintains tissue perfusion under perturbations that vary by factor of five or more in metabolic demand (rest to vigorous exercise) and by similar factors in posture, temperature, and intravascular volume. The regulatory architecture is layered across three timescales — neural (seconds), hormonal (minutes-to-hours), and renal-volume (hours-to-days) — each acting on common effector limbs (heart rate, contractility, vasomotor tone, blood volume) but tuned to a different perturbation spectrum.

Autonomic control of cardiac function. Sympathetic and parasympathetic divisions of the autonomic nervous system innervate the heart and vasculature with opposing actions. Sympathetic fibres release norepinephrine at ${\beta }_1$-adrenergic receptors on cardiac myocytes; receptor activation couples through $G_s$ to adenylyl cyclase, raising intracellular cAMP, activating protein kinase A (PKA), and triggering phosphorylation of L-type calcium channels (greater $I_{\mathrm{Ca},\mathrm{L}}$ per beat), phospholamban (faster SR calcium reuptake), cardiac troponin I (lower myofilament calcium sensitivity, faster relaxation), and RyR2 (modulated open probability). The integrated effect is positive chronotropy (faster heart rate), positive inotropy (greater contractility, steeper ESPVR), positive lusitropy (faster relaxation, allowing adequate filling at higher heart rates), and positive dromotropy (faster atrioventricular conduction).

Parasympathetic (vagal) fibres release acetylcholine at ${\mathrm{M}}_2$ muscarinic receptors, primarily on the sinoatrial node, atrial myocardium, and AV node. Receptor activation couples through $G_i$ to inhibit adenylyl cyclase (opposing sympathetic effects) and through $G_{\beta \gamma }$ directly to open GIRK potassium channels (immediate hyperpolarisation of pacemaker cells). The integrated effect is negative chronotropy, negative dromotropy, and a smaller negative inotropic effect on atrial muscle. Ventricular vagal innervation is sparse; the resting heart rate of $\sim 70\mathrm{bpm}$ in the presence of an intrinsic pacemaker frequency of $\sim 100-120\mathrm{bpm}$ reflects continuous resting vagal tone.

Peripheral vasomotor control. Sympathetic noradrenergic fibres innervate arterioles throughout the body. ${\alpha }_1$-adrenergic receptors on vascular smooth muscle couple through $G_q$ to phospholipase C, producing $\mathrm{I}{\mathrm{P}}_3$ and diacylglycerol, releasing SR calcium, activating myosin light-chain kinase, phosphorylating regulatory myosin light chains, and producing smooth-muscle contraction (vasoconstriction). Selected vascular beds carry ${\beta }_2$-adrenergic receptors (skeletal muscle, coronary, hepatic) that mediate vasodilation at low circulating epinephrine — providing the regional redistribution of cardiac output toward exercising muscle without overall hypotension. Resting sympathetic tone provides the baseline vasomotor tone; withdrawal of sympathetic tone (e.g., orthostatic vasovagal episodes, spinal cord injury) produces vasodilation and venous pooling.

The baroreceptor reflex. Stretch-sensitive afferents in the carotid sinus (cranial nerve IX) and aortic arch (cranial nerve X) fire at a rate that increases with vessel-wall stretch (proportional to mean arterial pressure with an added pulse-pressure-dependent component). Afferent input projects to the nucleus tractus solitarius (NTS) in the medulla, which projects in turn to the nucleus ambiguus (vagal motor output to the heart) and to the rostral ventrolateral medulla (sympathetic premotor output to the spinal cord). The reflex arc is fast: a pressure perturbation produces a corrective autonomic response within $\sim 1\mathrm{s}$. The reflex's operating range is $\sim 60-160\mathrm{mmHg}$ mean arterial pressure, with the steepest gain around the resting set point.

The reflex's signature failure mode is resetting. Sustained elevation of mean arterial pressure (over hours to days) shifts the baroreceptor pressure-firing-rate curve rightward, so the new elevated pressure becomes the new operating point. The reflex retains its short-term buffering capacity around the elevated set point but no longer corrects the chronic elevation. This resetting is the reason the baroreflex does not, by itself, defend the long-term arterial pressure against chronic hypertension; long-term defence requires the renal-volume axis (below) acting on the underlying blood volume rather than on its acute distribution.

A mathematical model of the baroreflex as a negative-feedback controller has the form

$$\dot{u}=-k_p(P-P_{\mathrm{set}})-k_i{\int }_0^t(P-P_{\mathrm{set}})d{t}^{\prime },$$

with $u$ the autonomic-output index (a one-dimensional summary of the sympathetic-parasympathetic balance), $P$ the measured mean arterial pressure, and $P_{\mathrm{set}}$ the medullary set point. The proportional gain $k_p$ controls the fast response (sub-second), the integral term $k_i$ provides the slow adaptation that produces the resetting. The set point $P_{\mathrm{set}}$ is itself slowly modifiable by higher centres (hypothalamic, cortical) — the basis of stress-induced and emotional cardiovascular responses.

Hormonal regulation — RAAS and ANP. The renin-angiotensin-aldosterone system (RAAS) is the dominant slow regulator. Renal juxtaglomerular cells release renin in response to (i) decreased renal perfusion pressure, (ii) decreased sodium delivery to the macula densa, or (iii) ${\beta }_1$-adrenergic stimulation. Renin cleaves angiotensinogen (a liver-produced precursor) to angiotensin I; angiotensin-converting enzyme (ACE, primarily pulmonary endothelial) cleaves angiotensin I to angiotensin II, the active mediator. Angiotensin II has multiple actions: direct vasoconstriction through AT${}_1$ receptors on vascular smooth muscle, stimulation of aldosterone release from the adrenal cortex (aldosterone promotes renal sodium retention, expanding extracellular fluid volume), promotion of antidiuretic hormone (ADH/vasopressin) release from the posterior pituitary (further water retention), and central effects on thirst and sympathetic outflow.

Atrial natriuretic peptide (ANP) is the counter-regulatory hormone — released by atrial myocytes in response to atrial stretch (a proxy for blood-volume overload). ANP promotes renal sodium excretion (natriuresis), vasodilation, and reduction in renin and aldosterone secretion. Together, RAAS and ANP form a slow negative-feedback loop on blood volume that operates on the timescale of hours-to-days. Pharmacological interruption of the RAAS axis (ACE inhibitors, angiotensin-receptor blockers, mineralocorticoid-receptor antagonists) is the most consequential therapeutic intervention in modern cardiovascular medicine, with mortality reductions across heart failure, hypertension, post-myocardial-infarction remodelling, and chronic kidney disease.

Exercise physiology. Dynamic exercise (running, cycling) presents the cardiovascular system with its largest physiological challenge. Oxygen consumption $V{O}_2$ rises from $\sim 0.25L/min$ at rest to $\sim 3L/min$ in moderately fit individuals and up to $\sim 6L/min$ in elite endurance athletes. The cardiovascular response involves: (i) sympathetic activation and parasympathetic withdrawal, raising heart rate from $\sim 70$ to $180-200\mathrm{bpm}$; (ii) increased contractility, raising stroke volume from $\sim 70$ to $120-150\mathrm{mL}$; (iii) increased cardiac output from $\sim 5$ to $\sim 25L/min$; (iv) regional redistribution of flow toward exercising muscle through ${\beta }_2$-mediated and metabolite-mediated (adenosine, ${\mathrm{H}}^{+}$, lactate, $K^{+}$) vasodilation of muscle arterioles and ${\alpha }_1$-mediated vasoconstriction of splanchnic and renal beds; (v) the muscle pump of contracting skeletal muscle squeezing veins, augmenting venous return; (vi) widening of arteriovenous oxygen difference as exercising muscle extracts more O${}_2$ per unit blood. The Fick principle $V{O}_2=\mathrm{CO}\cdot (a-v){\mathrm{O}}_2$ partitions the rise in oxygen consumption between cardiac-output increase and AV-difference widening; both contribute roughly equally to maximal exercise capacity.

Long-term endurance training induces cardiac hypertrophy with chamber dilation (an adaptive phenotype distinct from pathological pressure-overload hypertrophy), increasing maximum stroke volume and shifting the resting operating point — trained athletes have resting heart rates of $40-50\mathrm{bpm}$ at the same cardiac output as an untrained individual at $70\mathrm{bpm}$. The remodeling involves coordinated changes in myocyte size, capillary density, mitochondrial content, and SERCA2a expression — a multi-tissue, multi-protein response to the chronic demand signal of repeated exercise sessions.

Synthesis. Cardiovascular function builds toward 18.02.02 through the calcium-handling apparatus that translates the cardiac action potential of pacemaker and working myocardium into the calcium transient driving cross-bridge activation, and appears again in 18.04.02 pending as the same actin-myosin chemistry that powers muscle contraction generally — the cardiac chamber inheriting the molecular machinery and reshaping it via length-dependent activation and CICR amplification. The foundational reason cardiac output is self-regulating is that the heart is an elastic chamber whose end-systolic elastance and end-diastolic compliance together fix the operating PV loop; this is exactly the framework that lets the Frank-Starling law emerge as a constitutive consequence of length-dependent calcium sensitivity rather than as an independent regulatory principle. Putting these together with the Windkessel arterial model, the Poiseuille resistance, and the Moens-Korteweg pulse-wave speed, one obtains a closed lumped-parameter model of the cardiovascular system on a $\sim 10$-dimensional state space (chamber volumes, arterial pressures, sympathetic and vagal tones, blood-volume hormones), and the bridge is to the kidney-volume axis via RAAS / ANP that closes the long-term loop. The central insight is that fluctuation buffering at every timescale — beat-to-beat by the Frank-Starling law, breath-to-breath by the baroreflex, day-to-day by the RAAS axis — identifies the cardiovascular system as a hierarchically nested set of negative-feedback controllers on the same regulated variable (mean arterial pressure, with mean cardiac output as its rate-limited proxy), each controller dominant at its characteristic timescale and the pattern recurring across organ systems as the general design principle of redundant homeostasis.

Full proof set [Master]

Proposition 1 (Stability of the Frank-Starling fixed point for the coupled ventricles). Let $Q_L(V_p)=E_L\cdot (V_p-V_0^L)$ be the left-ventricular output as a monotonically increasing function of pulmonary venous volume $V_p$ (the left-side preload proxy), and $Q_R(V_s)=E_R\cdot (V_s-V_0^R)$ the right-ventricular output as a monotonically increasing function of systemic venous volume $V_s$. Let total blood volume satisfy $V_s+V_p=V_{\mathrm{tot}}$ (a fixed conservation law on the timescale of seconds, neglecting arterial volumes for clarity). Then the system $d{V}_p/dt=Q_R-Q_L$, $V_s=V_{\mathrm{tot}}-V_p$ has a unique stable fixed point at which $Q_L=Q_R$.

Proof. Substituting the conservation law,

$$\frac{d{V}_p}{dt}=Q_R(V_{\mathrm{tot}}-V_p)-Q_L(V_p)=E_R(V_{\mathrm{tot}}-V_p-V_0^R)-E_L(V_p-V_0^L).$$

The right-hand side is affine in $V_p$ with slope $-(E_L+E_R)<0$. There is a unique zero (the fixed point) at

$$V_p^{\ast }=\frac{E_R(V_{\mathrm{tot}}-V_0^R)+E_L{V}_0^L}{E_L+E_R},$$

and at this fixed point $Q_L=Q_R$ by construction. The negative slope of the right-hand side at the fixed point gives linear stability: a small perturbation $\delta V_p$ relaxes exponentially with rate $E_L+E_R$. Nonlinearly, the right-hand side remains monotonically decreasing in $V_p$ as long as both $Q_L$ and $Q_R$ are monotonically increasing in their respective preloads (the Frank-Starling assumption), so the fixed point is globally attracting on the physiological range. $\square$

Remark. The fixed point's existence and stability do not require neural or hormonal regulation. The mechanism is purely mechanical: the more compliant side of the heart receives more volume and ejects more, until the two outputs match. This is the formalisation of the clinical observation that the heart "pumps what it receives" — and the reason an autonomously beating, denervated, transplanted heart can sustain normal cardiac output at rest.

Proposition 2 (Diastolic decay of aortic pressure in the two-element Windkessel). In the two-element Windkessel model, during diastole ($Q_{\mathrm{in}}=0$), aortic pressure decays exponentially from its end-systolic value to its asymptotic value $P_{\mathrm{venous}}$ with time constant $\tau =RC$, where $R$ is peripheral resistance and $C$ arterial compliance.

Proof. Setting $Q_{\mathrm{in}}(t)=0$ in the Windkessel ODE,

$$C\frac{dP}{dt}+\frac{P-P_{\mathrm{venous}}}{R}=0.$$

Substituting $p=P-P_{\mathrm{venous}}$,

$$\frac{dp}{dt}=-\frac{1}{RC}p,p(t)=p(0)e^{-t/RC}.$$

The end-systolic boundary condition $p(0)=P_{\mathrm{es}}-P_{\mathrm{venous}}$ gives

$$P(t)=P_{\mathrm{venous}}+(P_{\mathrm{es}}-P_{\mathrm{venous}})e^{-t/RC}.$$

The time constant $\tau =RC$ is independent of stroke volume and initial conditions, so $\tau$ can be measured from the diastolic decay slope independently of beat-to-beat variation. With $P_{\mathrm{venous}}\approx 0$ (the standard clinical approximation), the formula reduces to the textbook $P(t)=P_{\mathrm{es}}{e}^{-t/RC}$. $\square$

Proposition 3 (Pressure-flow relationship in steady Poiseuille flow). Steady laminar flow of a Newtonian fluid of viscosity $\eta$ through a long cylindrical tube of radius $r$, length $L$, with constant pressure drop $\Delta P$ end-to-end, has volumetric flow rate $Q=\pi r^4\Delta P/(8\eta L)$, hence resistance $R=8\eta L/(\pi r^4)$.

Proof. Inside the tube, the axisymmetric steady-state Navier-Stokes equation reduces (under the assumptions: no swirl, no axial variation other than imposed pressure gradient $-dP/dz=\Delta P/L$, no slip at the wall $r^{\prime }=r$) to

$$\eta \frac{1}{r^{\prime }}\frac{d}{d{r}^{\prime }}\left(r^{\prime }\frac{du}{d{r}^{\prime }}\right)=-\frac{\Delta P}{L},$$

with $u(r)=0$ (no slip) and $u(0)$ finite. Integrating twice,

$$u(r^{\prime })=\frac{\Delta P}{4\eta L}(r^2-r^{\prime 2}).$$

The velocity profile is parabolic with maximum $u(0)=\Delta P{r}^2/(4\eta L)$ at the centre. The volumetric flow rate is

$$Q={\int }_0^{r}2\pi r^{\prime }u(r^{\prime })d{r}^{\prime }=\frac{\pi \Delta P}{2\eta L}{\int }_0^r(r^2{r}^{\prime }-r^{\prime 3})d{r}^{\prime }=\frac{\pi \Delta P}{2\eta L}\frac{r^4}{4}=\frac{\pi r^4\Delta P}{8\eta L}.$$

Resistance $R=\Delta P/Q=8\eta L/(\pi r^4)$. The fourth-power radius dependence follows directly from the integration of the parabolic profile and the area element $r^{\prime }d{r}^{\prime }$ in cylindrical coordinates. $\square$

Remark. The $r^4$ law is a robust mathematical consequence of axisymmetric Newtonian flow with no slip; it does not depend on the specific viscosity value, only on viscosity being constant. Departures from $r^4$ in vivo arise from non-Newtonian effects (small-vessel viscosity decrease via the Fåhraeus-Lindqvist effect), pulsatility (Womersley correction), and entrance effects in short vessel segments.

Connections [Master]

• Cardiac action potentials and pacemaker physiology 18.02.02. The membrane-electrical event that initiates each heartbeat. The L-type calcium current shaping the action-potential plateau is the trigger for the calcium-induced calcium release of cardiac excitation-contraction coupling treated in the Master subsection above. The sinoatrial-node pacemaker mechanism sets the heart rate that the autonomic and hormonal regulators modulate.

• Muscle contraction — the actin-myosin cycle 18.04.02 pending. The molecular substrate of cardiac contraction. Cardiac myocytes deploy the same Lymn-Taylor cross-bridge cycle, with cardiac-specific isoforms (cTnC, cTnI, cTnT, $\beta$-myosin heavy chain) tuning the kinetics and calcium sensitivity. Length-dependent activation of the sarcomere — the molecular origin of the Frank-Starling law — is one of the cardiac-specific refinements of the general cross-bridge machinery.

• Skeletal muscle physiology 18.04.01. The companion muscle physiology that shares the actin-myosin core with cardiac muscle but differs in regulation (mechanical DHPR-RyR1 coupling vs. cardiac CICR), in syncytial organisation, in fatigability, and in fibre-type composition. The comparison clarifies which aspects of cardiac mechanics are inherited from the general cross-bridge cycle and which are cardiac specialisations.

• Respiratory physiology and gas exchange 18.03.01 pending. Cardiac output delivers blood to the pulmonary capillary bed where gas exchange occurs. Ventilation-perfusion matching is a downstream organ-system integration; the haemodynamic quantities (CO, pulmonary arterial pressure, pulmonary vascular resistance) defined here are the inputs to that downstream system.

• Renal physiology and homeostasis 18.08.01 pending. Closes the long-term blood-pressure feedback loop. The renin-angiotensin-aldosterone system originates in the kidney and acts on the cardiovascular effectors discussed in the Master regulation subsection. Renal pressure-natriuresis is the slow homeostatic mechanism that defends the long-term arterial-pressure set point.

• Endocrine hormones and regulation 18.07.01. Adrenal catecholamines, thyroid hormone, antidiuretic hormone, and the natriuretic peptide family all act on cardiovascular function. The sympathetic-adrenal-medullary axis discussed in the regulation subsection is one branch of this broader endocrine integration.

• Cellular neuroscience — resting membrane potential and ion channels 17.09.01. The cellular biophysics of excitable membranes underlies the cardiac action potential. The Nernst and Goldman-Hodgkin-Katz equations developed in the cellular-neuroscience strand apply directly to cardiac myocytes; cardiac-specific ion channels (Ca${}_v$1.2, $\mathrm{HCN}$ for pacemaking, the cardiac potassium-channel families) are the cardiac specialisations of the general ion-channel framework.

• Oxidative phosphorylation and ATP synthesis 17.04.02 pending. The mitochondrial ATP supply that fuels the cross-bridge cycle and the SERCA2a calcium pump. Cardiac myocytes are remarkable for their density of mitochondria ($\sim 35\%$ of cell volume in working myocardium) and their near-exclusive dependence on oxidative metabolism — a single cardiac myocyte may consume $\sim 30$ ATP per second per cross-bridge under load, with negligible glycolytic capacity to fall back on.

Historical & philosophical context [Master]

William Harvey's Exercitatio Anatomica de Motu Cordis et Sanguinis in Animalibus (1628) established that blood circulates in a closed loop driven by the heart, overturning the Galenic doctrine that blood was consumed at the periphery and produced anew in the liver. Harvey's argument was quantitative: he estimated cardiac output (heart capacity times beat frequency) and observed that the daily output exceeded by orders of magnitude what the liver could possibly manufacture, forcing the conclusion that the same blood must return. The capillary anastomoses connecting arteries to veins, predicted by Harvey but not visualised in his lifetime, were demonstrated by Marcello Malpighi in 1661 using the newly invented microscope on frog lung tissue. Harvey's quantitative-circulatory framework remains the conceptual foundation of cardiovascular physiology to the present.

Otto Frank's Zur Dynamik des Herzmuskels in Zeitschrift für Biologie 32 (1895) ^{[Frank 1895]} demonstrated the length-tension relationship in isolated frog heart muscle: greater initial stretch of the ventricular wall produced greater force of contraction. Frank's pressure-volume diagrams from that paper are the earliest published PV loops and contain the structural ingredients of the modern Suga-Sagawa framework — the looping trajectory, the ESPVR upper boundary, the EDPVR lower boundary. Ernest Starling at University College London extended Frank's work to the mammalian heart-lung preparation and articulated the relationship as a general law in The Linacre Lecture on the Law of the Heart ^{[Starling 1918]}. Starling's framing — that the heart's output is determined by its filling, and that the heart adjusts its output to its input without external regulation — became the operational principle now called the Frank-Starling law. The relationship anchors clinical hemodynamics and remains one of the most experimentally robust regulatory principles in cardiovascular physiology.

Carl Wiggers at Western Reserve University in American Journal of Physiology 56 (1921) ^{[Wiggers 1921]} produced the canonical time-aligned diagrams of left ventricular pressure, aortic pressure, left atrial pressure, ventricular volume, electrocardiogram, and phonocardiogram across one cardiac cycle. The Wiggers diagram is the standard pedagogical depiction of the cardiac cycle and the format in which most subsequent measurements have been displayed. Willem Einthoven's 1903 invention of the string galvanometer and the resulting body-surface electrocardiogram earned the 1924 Nobel Prize and gave clinicians their first non-invasive window into cardiac electrical activity.

The Windkessel concept originated with Stephen Hales in Statical Essays: Containing Haemastaticks (1733), where he observed that the elastic large arteries smooth the pulsatile output of the heart into a more continuous flow at the periphery, analogous to the air chamber (Windkessel) used in fire-fighting hand pumps of the period to convert intermittent pumping into continuous spray. The mathematical formalisation as a two-element electrical-analogue model is due to Otto Frank's 1899 monograph Die Grundform des arteriellen Pulses; the three-element extension is due to Nicolaas Westerhof in 1969. The Moens-Korteweg pulse-wave-velocity equation was derived independently by Adriaan Moens (Leiden, 1878) and Diederik Korteweg (Amsterdam, 1878) ^{[Moens 1878]} within months of each other from a one-dimensional fluid-structure interaction model of a thin-walled elastic tube — the same Korteweg who later co-discovered the Korteweg-de Vries equation for shallow-water waves.

The calcium-induced calcium release mechanism that integrates cardiac action potentials with myofilament activation was established quantitatively in skinned cardiac myocytes by Alexandre Fabiato in Journal of General Physiology 85 (1985) ^{[Fabiato 1985]}, following two decades of work distinguishing cardiac from skeletal-muscle excitation-contraction coupling. The mechanism's molecular substrate (the dyadic cleft, RyR2 cluster architecture, SERCA2a / phospholamban regulation) was elucidated over the subsequent two decades through structural and physiological work in the Bers, Lederer, Cannell, and Marks laboratories among many others. Hiroyuki Suga and Kiichi Sagawa's time-varying elastance framework, formalised in a series of papers in Circulation Research in the 1970s (Suga 1971 Am. J. Physiol.; Suga 1979 Circ. Res. on PVA-$V{O}_2$), elevated the PV loop from a descriptive diagram to a quantitative chamber-mechanics framework. Modern computational cardiology builds on the Hodgkin-Huxley-style ionic-current models of cardiac action potentials initiated by Denis Noble in 1962 ^{[Noble 1962]} — the first detailed mathematical model of any cardiac cell — and now extended to whole-organ electromechanical simulations.

Bibliography [Master]

@book{HarveyDeMotuCordis1628,
  author = {Harvey, William},
  title = {Exercitatio Anatomica de Motu Cordis et Sanguinis in Animalibus},
  publisher = {Guilielmi Fitzeri},
  address = {Frankfurt},
  year = {1628},
}

@article{Frank1895,
  author = {Frank, Otto},
  title = {Zur Dynamik des Herzmuskels},
  journal = {Zeitschrift f\"ur Biologie},
  volume = {32},
  pages = {370--437},
  year = {1895},
}

@book{Starling1918,
  author = {Starling, Ernest Henry},
  title = {The Linacre Lecture on the Law of the Heart},
  publisher = {Longmans, Green \& Co.},
  address = {London},
  year = {1918},
}

@article{Wiggers1921,
  author = {Wiggers, Carl J.},
  title = {Studies on the consecutive phases of the cardiac cycle},
  journal = {American Journal of Physiology},
  volume = {56},
  pages = {415--438},
  year = {1921},
}

@book{Hales1733,
  author = {Hales, Stephen},
  title = {Statical Essays: Containing Haemastaticks},
  publisher = {W. Innys},
  address = {London},
  year = {1733},
}

@article{Moens1878,
  author = {Moens, Adriaan I.},
  title = {Die Pulskurve},
  publisher = {Brill},
  address = {Leiden},
  year = {1878},
}

@article{Fabiato1985,
  author = {Fabiato, Alexandre},
  title = {Time and calcium dependence of activation and inactivation of calcium-induced release of calcium from the sarcoplasmic reticulum of a skinned canine cardiac Purkinje cell},
  journal = {Journal of General Physiology},
  volume = {85},
  pages = {247--289},
  year = {1985},
}

@article{Suga1979,
  author = {Suga, Hiroyuki},
  title = {Total mechanical energy of a ventricle model and cardiac oxygen consumption},
  journal = {American Journal of Physiology},
  volume = {236},
  pages = {H498--H505},
  year = {1979},
}

@article{Noble1962,
  author = {Noble, Denis},
  title = {A modification of the Hodgkin-Huxley equations applicable to Purkinje fibre action and pacemaker potentials},
  journal = {Journal of Physiology},
  volume = {160},
  pages = {317--352},
  year = {1962},
}

@article{Westerhof1969,
  author = {Westerhof, Nicolaas and Bosman, Frans and De Vries, Christiaan J. and Noordergraaf, Abraham},
  title = {Analog studies of the human systemic arterial tree},
  journal = {Journal of Biomechanics},
  volume = {2},
  pages = {121--143},
  year = {1969},
}

@book{Katz2022,
  author = {Katz, Arnold M.},
  title = {Physiology of the Heart},
  edition = {6th},
  publisher = {Lippincott Williams \& Wilkins},
  year = {2022},
}

@book{BoronBoulpaep2017,
  author = {Boron, Walter F. and Boulpaep, Emile L.},
  title = {Medical Physiology},
  edition = {3rd},
  publisher = {Elsevier},
  year = {2017},
}

@book{GuytonHall2020,
  author = {Guyton, Arthur C. and Hall, John E.},
  title = {Textbook of Medical Physiology},
  edition = {14th},
  publisher = {Elsevier},
  year = {2020},
}

@book{Bers2001,
  author = {Bers, Donald M.},
  title = {Excitation-Contraction Coupling and Cardiac Contractile Force},
  edition = {2nd},
  publisher = {Springer},
  year = {2001},
}

@book{NicholsORourke2011,
  author = {Nichols, Wilmer W. and O'Rourke, Michael F. and Vlachopoulos, Charalambos},
  title = {McDonald's Blood Flow in Arteries: Theoretical, Experimental and Clinical Principles},
  edition = {6th},
  publisher = {CRC Press},
  year = {2011},
}