Modeling of a Windkessel Windkessel physiology remains a relevant yet dated description of important clinical interest. The historic mathematical definition of systole and diastole in the model are obviously not novel but are here elementally staged to four degrees. Reaching five would be original work.
Two-element It is assumed that the ratio of pressure to volume is constant and that outflow from the Windkessel is proportional to the fluid pressure. Volumetric inflow must equal the sum of the volume stored in the capacitive element and volumetric outflow through the resistive element. This relationship is described by a
differential equation: I(t)={P(t)\over R}+C{dP(t)\over dt}
I(t) is volumetric inflow due to the pump (heart) and is measured in volume per unit time, while
P(t) is the pressure with respect to time measured in force per unit area,
C is the ratio of volume to pressure for the Windkessel, and
R is the resistance relating outflow to fluid pressure. This model is identical to the relationship between current,
I(t), and
electrical potential,
P(t), in an electrical circuit equivalent of the two-element Windkessel model. In the blood circulation, the passive elements in the circuit are assumed to represent elements in the
cardiovascular system. The resistor,
R, represents the total peripheral resistance and the capacitor,
C, represents total arterial compliance. During
diastole there is no blood inflow since the aortic (or pulmonary valve) is closed, so the Windkessel can be solved for
P(t) since
I(t) = 0: P(t)=P(t_d)e^{-(t-t_d)\over (RC)} where
td is the time of the start of
diastole and
P(td) is the blood pressure at the start of diastole. This model is only a rough approximation of the arterial circulation; more realistic models incorporate more elements, provide more realistic estimates of the blood pressure waveform and are discussed below.
Three-element The three-element Windkessel improves on the two-element model by incorporating another resistive element to simulate resistance to blood flow due to the characteristic resistance of the aorta (or pulmonary artery). The
differential equation for the 3-element model is: (1+{R_1\over R_2})I(t)+CR_1{dI(t)\over dt}= {P(t)\over R_2}+C{dP(t)\over dt} where
R1 is the characteristic resistance (this is assumed to be equivalent to the characteristic impedance), and the pulmonary artery in a pig
Four-element The three-element model overestimates the compliance and underestimates the characteristic impedance of the circulation. The four-element model includes an
inductor,
L, which has units of mass per length, ({M\over l^4}), into the proximal component of the circuit to account for the
inertia of blood flow. This is neglected in the two- and three- element models. The relevant equation is: (1+{R_1\over R_2})I(t)+(R_1C+{L\over R_2}){dI(t)\over dt}+LC{d^2I(t)\over dt^2}={P(t)\over R_2}+C{dP(t)\over dt} == Applications ==