Deriving Pressure-Flow Relations in Lumped Circulatory Networks

Pressure Differences Drive Flow in a Lumped Model

A circulation model needs a rule connecting vascular pressure to blood flow. For a segment represented by one effective resistance, the rule is simple: flow is proportional to the pressure difference across that segment. In a network, this relation is paired with conservation of flow at each node. Together, the two rules let us calculate both the pressure at an internal junction and the flows through the beds connected to it.

The earlier material on coupled cardiovascular-renal state variables and feedback topology places circulatory variables within a larger physiological system. Here the focus is narrower: a steady, lumped pressure-flow calculation. It does not represent transient storage of blood in compliant vessels or the time course of feedback.

The sign convention matters. A positive \(Q_{ij}\) means flow from \(i\) to \(j\); a negative result means the actual flow is opposite to the direction chosen for writing the equation. For a positive resistance, flow is always directed from higher pressure toward lower pressure.

Why the Relation Holds, and What It Assumes

Resistance is defined in the linear model as the pressure drop required per unit flow: \(\Delta P=RQ\). Rearranging gives \(Q=\Delta P/R\). This is a constitutive relation for a lumped segment: all the detailed effects of vessel geometry and blood properties are summarized by one effective parameter, \(R\).

A familiar physical basis is Poiseuille flow through a rigid cylindrical tube containing a Newtonian fluid under steady laminar conditions. In that idealization, \(R=8\mu L/(\pi r^4)\), where \(\mu\) is viscosity, \(L\) is tube length, and \(r\) is radius. A vascular bed is not a single rigid tube, so its effective resistance is generally a model parameter rather than a direct description of one vessel. The idealized law nevertheless motivates the proportional relation between pressure drop and flow.

The calculation relies on several assumptions. Pressures at each node are treated as spatially uniform; each segment has a specified, positive resistance; and the resistance is constant over the range of pressures and flows being considered. The calculation is steady, so there is no net accumulation of blood at a node. Finally, the network must be connected to pressure conditions that determine the unknown node pressures. If these assumptions are not reasonable, the algebra can still be solved but may not describe the physiology adequately.

Flow Conservation at a Node

A node is a junction where segments meet. In a steady model with no source, sink, or storage at the junction, total flow entering must equal total flow leaving. For a node \(k\), this can be written as the sum of signed flows equal to zero. Each flow term is calculated from the pressure difference across its own segment.

In this expression, each term is defined as flow from node \(k\) toward a connected node \(j\). Terms are positive when pressure at \(k\) is higher and negative when the connected node has the higher pressure. This signed form handles either direction without needing to guess the direction in advance.

Two Vascular Beds Connected to a Shared Node

Consider a source at pressure \(P_s\) connected through an upstream resistance \(R_0\) to a shared arterial node at unknown pressure \(P_n\). Two vascular beds leave that node and drain to a common venous pressure \(P_v\). Let their resistances be \(R_1\) and \(R_2\), and define their flows as \(Q_1\) and \(Q_2\), directed from the node toward the venous side. Define \(Q_0\) as flow from the source to the node.

Because the node has no storage or additional inflow or outflow in this steady model, its balance is \(Q_0=Q_1+Q_2\). Substitution produces one equation for the unknown node pressure. After solving for \(P_n\), calculate each flow separately and verify the balance numerically. This procedure keeps the pressure and flow for each bed explicit rather than assuming a particular flow split.

Derivation of the Node Equation and Practical Checks

The two-bed equation follows directly from conservation rather than from an assumed pressure at the internal node. In the example, the source flow entering the node is \((P_s-P_n)/R_0\). The two flows leaving are \((P_n-P_v)/R_1\) and \((P_n-P_v)/R_2\). Equating entry and exit yields:

For the stated positive resistances, the left side decreases as \(P_n\) increases, while the right side increases as \(P_n\) increases. Therefore, when the source pressure exceeds the venous pressure, these two sides can intersect at no more than one node pressure. In this connected network, solving the linear equation gives that pressure; calculating the flows afterward supplies an independent conservation check.

Useful diagnostics are straightforward. First check units: dividing mmHg by mmHg·min/L must give L/min. Next check signs: a positive flow should correspond to a pressure drop in the stated direction. Finally, calculate the node residual, \(Q_0-Q_1-Q_2\). A residual other than zero beyond rounding suggests an arithmetic, sign, or equation-setup error. A computed negative flow is not automatically an error; it indicates flow opposite to the direction initially defined and should prompt a check that the model structure and pressure boundary conditions allow that direction.

Choosing and Interpreting the Model

Direct segment calculations are appropriate when both endpoint pressures and the resistance are known. If an internal pressure is unknown, the node-balance equations provide the appropriate extension: write a conservation equation for each unknown internal node and solve the simultaneous system. A dynamic model is a different approach, appropriate when pressures and flows change over time and storage or feedback cannot be ignored. Such a model requires additional state variables and equations; the steady calculation alone cannot predict a transient response.

The lumped linear relation is useful for transparent accounting and for checking whether proposed pressures, resistances, and flows are mutually consistent. It has important limitations. Real vessels can change radius, resistance can depend on flow or physiological state, and pressure and flow can vary over the cardiac cycle. A single resistance and a steady pressure may therefore be inadequate for a question about pulsatility, changing vascular tone, or local flow regulation. In those settings, the values calculated here should be treated as model-dependent approximations, not as universal physiological constants.

A common mistake is to calculate each bed flow from the source pressure rather than from the pressure difference across that bed. For example, using \((P_s-P_v)/R_1\) in the two-bed example incorrectly assigns the full source-to-vein pressure drop to bed 1, ignoring the upstream resistance and the shared node pressure. The correction is to identify each segment’s actual endpoint pressures, write its own pressure difference, and enforce conservation at the junction.