Vol. 83, Issue 6, 1799-1813, December 1997
Coupled vs. uncoupled pericardial constraint: effects on
cardiac chamber interactions
Masao
Takata,
Yasuhiko
Harasawa,
Sadek
Beloucif, and
James L.
Robotham
Pathophysiology Research, National Children's Medical Research
Center, Tokyo 154; Division of Cardiology, Kyushu University
Hospital, Fukuoka 812-82, Japan; and Department of Anesthesiology
and Critical Care Medicine, The Johns Hopkins Medical Institutions, Baltimore, Maryland 21205
ABSTRACT
INTRODUCTION
GLOSSARY
METHODS
RESULTS
DISCUSSION
APPENDIX
ACKNOWLEDGEMENTS
FOOTNOTES
REFERENCES
ABSTRACT
Takata, Masao, Yasuhiko Harasawa, Sadek Beloucif, and James
L. Robotham. Coupled vs. uncoupled pericardial constraint: effects
on cardiac chamber interactions. J. Appl.
Physiol. 83(6): 1799-1813, 1997.
The effects of
pericardial constraint on cardiac chamber interactions were evaluated
by mathematical model analyses based on a novel concept of coupled vs.
uncoupled pericardial constraint. We hypothesized that the nature of
pericardial constraint can be classified as a "coupled"
constraint exerted by uniform liquid pressure or an "uncoupled"
constraint exerted by regional surface pressure. The numerical solution
of the model of atrioventricular interaction produced the
characteristic waveforms in venous flows and right atrial/ventricular
pressures in classical pericardial diseases. Coupled constraint
accounted for the patterns in cardiac tamponade; uncoupled constraint
accounted for those in constrictive pericarditis. Analytic solution of
the model of ventricular interdependence demonstrated that coupled
constraint (tamponade) produced greater gains in ventricular
interdependence, increasing the occurrence of pulsus paradoxus, whereas
uncoupled constraint (constriction) produced a greater effective right
ventricular elastance, increasing the likelihood of Kussmaul's sign.
Thus the concept of coupled vs. uncoupled constraint may offer a
coherent framework to understand the characteristic steady-state and
respiratory-induced hemodynamic events in multiple forms of pericardial
diseases.
atrioventricular interaction; constrictive pericarditis; model; tamponade; venous return; ventricular interdependence
INTRODUCTION
CARDIAC CHAMBERS are physically linked to each other,
hydrodynamically connected by blood flows, and situated in a space
limited by the juxtacardiac structures. Filling and emptying of one
cardiac chamber are thus influenced by changes in volume or pressure of other chambers. The pericardium has been shown to modulate such cardiac
chamber interactions (13): the vertical interaction between an atrium
and the corresponding ventricle within a cardiac cycle, i.e.,
atrioventricular interaction (4), and the horizontal interaction
between the two ventricles, i.e., ventricular interdependence (6, 21).
However, the mechanism whereby the pericardium in disease states
influences such cardiac chamber interactions has not been well
characterized.
The differences in hemodynamic profiles between cardiac tamponade and
constrictive pericarditis have fascinated clinicians and physiologists
for many years. There are classic observations of differences in
vascular pressure waveforms in pericardial diseases, i.e., a prominent
x-descent in atrial pressure with
tamponade (18, 23) vs. a large
y-descent in atrial pressure and a
"dip-and-plateau" pattern in ventricular pressure with
constriction (19, 23). Recent studies demonstrated that venous flow
waveforms, normally biphasic with one peak during ventricular systole
and the other during ventricular diastole (2, 25, 27), are also altered in pericardial diseases. With tamponade, we found in dogs that venous
flows become mainly systolic with an almost absent diastolic flow (4),
consistent with Doppler flow velocity studies in humans (3, 5, 8, 11).
With constriction, an enhanced diastolic venous flow with a relatively
small systolic component was observed by Doppler studies in humans (5,
8, 9, 11). These findings suggest that the manner in which the
pericardium constrains the heart may be quite different between
tamponade and constriction, despite similarly increased pericardial
pressure. The differences in venous pressure and flow waveforms may be
related to the unique nature of pericardial constraint in each
pathology and its influence on atrioventricular interaction.
In addition to those steady-state pressure and flow changes, it has
been well documented in the literature that respiratory-induced hemodynamic signs are manifested differently in cardiac tamponade and
constrictive pericarditis. Pulsus paradoxus, i.e., an inspiratory decrease in arterial pressure associated with a decrease in left ventricular stroke volume, is observed much more often in tamponade than in constriction (22). Conversely, Kussmaul's sign, i.e., a
paradoxical inspiratory increase in right atrial pressure, is occasionally associated with constriction but rarely with tamponade (18, 19). Respiration produces large changes in right heart volume due
to an inspiratory increase in systemic venous return (25, 27), which
then decrease left heart filling via ventricular interdependence,
leading to pulsus paradoxus (7, 17). Thus it seems a reasonable
assumption that the status of ventricular interdependence is different
in nature or degree between tamponade and constriction. The specific
manner of pericardial constraint and its effects on ventricular
interdependence possibly could explain the differences in incidence of
Kussmaul's sign in the two forms of pericardial diseases.
The goals of the present study are to better understand the
characteristic hemodynamic signs in pericardial diseases and, more
generally, to develop a conceptual framework characterizing the
influence of the pericardial constraint on cardiac chamber interactions. We hypothesized that the nature of the pericardial constraint may be classified as one of two types:
1) a "coupled" constraint
exerted by uniform liquid pressure over the heart (e.g., cardiac
tamponade), which couples changes in all cardiac chamber volumes, and
2) an "uncoupled" constraint
exerted by regional surface pressure over the heart (e.g., constrictive
pericarditis), which independently restricts each chamber volume. On
the basis of this hypothesis, the effects of coupled vs. uncoupled
pericardial constraint on cardiac chamber interactions were studied by
use of mathematical model analyses. Flow-mediated atrioventricular interaction was studied with a numerical model and computer simulation, whereas pressure-mediated ventricular interdependence was evaluated with an analytic model and symbolic mathematical analysis. The results
suggest that the construct of coupled vs. uncoupled pericardial constraint may offer a physiological rationale to interpret the differences in the hemodynamic profiles of cardiac tamponade and constrictive pericarditis and provide insights necessary to develop better diagnostic and therapeutic strategies.
GLOSSARY
Model of atrioventricular interaction
| Cp |
Pulmonary arterial capacitance
|
| Cv |
Systemic venous capacitance
|
| Epe |
Pericardial elastance over RA and RV for coupled constraint
|
| Epera |
Pericardial elastance over RA for uncoupled constraint
|
| Eperv |
Pericardial elastance over RV for uncoupled constraint
|
| Era |
Time-varying elastance of RA
|
| Eramax |
Maximum elastance of RA
|
| Eramin |
Minimum elastance of RA
|
| Erv |
Time-varying elastance of RV
|
| Ervmax |
Maximum elastance of RV
|
| Ervmin |
Minimum elastance of RV
|
| Lv |
Systemic venous inertance
|
| Pd |
Downstream pressure source
|
| Pp |
Pressure at the pulmonary arterial capacitance
|
| Ppe |
Pericardial pressure over RA and RV for coupled constraint
|
| Ppera |
Pericardial pressure over RA for uncoupled constraint
|
| Pperv |
Pericardial pressure over RV for uncoupled constraint
|
| Pra |
RA pressure
|
| Pratm |
Transmural RA pressure
|
| Prv |
RV pressure
|
| Prvtm |
Transmural RV pressure
|
| Pu |
Upstream pressure source
|
| Pv |
Pressure at the systemic venous capacitance
|
| PRI |
P-R interval
|
p |
Pulmonary arterial flow
|
t |
Tricuspid flow
|
v |
Combined vena caval flow
|
| Rpd |
Distal pulmonary arterial resistance
|
| Rpp |
Proximal pulmonary arterial resistance
|
| Rt |
Transtricuspid valve resistance
|
| Rvd |
Distal systemic venous resistance
|
| Rvp |
Proximal systemic venous resistance
|
| RA |
Right atrium
|
| RV |
Right ventricle
|
| T |
Cardiac cycle time
|
| Tmaxra |
Contraction time of RA
|
| Tmaxrv |
Contraction time of RV
|
| Vra |
Volume of RA above its unstressed volume
|
| Vrv |
Volume of RV above its unstressed volume
|
| Vtotal |
Total right heart volume above its unstressed volume
|
Model of ventricular interdependence
| Elvf |
Elastance of LV free wall
|
| Epe |
Pericardial elastance over LV and RV for coupled constraint
|
| Epelv |
Pericardial elastance over LV for uncoupled constraint
|
| Eperv |
Pericardial elastance over RV for uncoupled constraint
|
| Epetotal |
Elastance of the total pericardium surrounding both ventricles
|
| Erveff |
Effective RV elastance with interdependence
|
| Ervf |
Elastance of RV free wall
|
| Es |
Elastance of the septum
|
| GP |
Right-to-left pressure interdependence gain
|
| GV |
Right-to-left volume interdependence gain
|
| LV |
Left ventricle
|
| Plv |
LV pressure
|
| Ppe |
Pericardial pressure over LV and RV for coupled constraint
|
| Ppelv |
Pericardial pressure over LV for uncoupled constraint
|
| Pperv |
Pericardial pressure over RV for uncoupled constraint
|
| Prv |
RV pressure
|
| Vlv |
Volume of LV
|
| Vlvf |
Volume of LV when transseptal pressure is zero
|
| Vrv |
Volume of RV
|
| Vrvf |
Volume of RV when transseptal pressure is zero
|
| Vs |
Volume contribution of the septal shift to either ventricle when
transseptal pressure is not zero
|
METHODS
Concept of Coupled vs. Uncoupled Pericardial Constraint
To characterize the nature of the pericardial constraint, we have
utilized a classic concept in pulmonary mechanics, i.e., surface vs.
liquid pressure (1, 16, 24, 26). Surface pressure is a force per unit
area acting on a surface and equals the sum of any liquid pressure and
any deformational forces produced locally by the apposition of two
surfaces. When the pericardial pressure is a pure liquid pressure
produced by a liquid column in the pericardial space, the pressure over
all cardiac chambers is homogeneous. Any change in the pressure over
one chamber will affect the pressure over all other chambers
effectively instantaneously. In contrast, when there is no liquid
column or it is not of sufficient volume, the heart surface is
macroscopically in contact with the pericardium. The resultant
pericardial surface pressure is essentially of regional nature, such
that the pressure on each chamber may not be equal in a short period of
time. Changes in pericardial pressure over a particular area of the
heart are determined mainly by local events occurring within that area.
On the basis of this reasoning, we classified the nature of pericardial
constraint as 1) a coupled
constraint exerted by uniform pericardial liquid pressure over all four
cardiac chambers or 2) an uncoupled
constraint produced by independent and different local pericardial
surface pressures over each cardiac chamber. The essence of this
concept is that a coupled constraint restricts the volumes of four
cardiac chambers together, whereas an uncoupled constraint
independently restricts the volume of each cardiac chamber. The present
study evaluates how such differences in the nature of pericardial
constraint can be predicted to modify the status of cardiac chamber
interactions.
Numerical Approach: Model of Atrioventricular Interaction
Analyses of vertical cardiac chamber interactions require consideration
of direct flow interaction between an atrium and the corresponding
ventricle through the atrioventricular valve and time-varying elastic
properties of the chambers within a cardiac cycle. To evaluate the
effects of coupled vs. uncoupled pericardial constraint on
atrioventricular interaction, we developed an open-loop mathematical
model of the right heart circulation. Numerical solution of the model
was performed by computer simulation, and the results were compared
with the experimental data in the literature regarding steady-state
venous pressure and flow waveforms in pericardial diseases.
Model description.
The model is a lumped, linear parameter, open-loop model including all
elements of right heart circulation (Fig.
1). The right atrium (RA) and ventricle
(RV) were characterized as time-varying elastances
(Era and
Erv) (10, 14, 20), which relate
the transmural pressure (relative to pericardial pressure) to the stressed volume of each chamber. The preloading and afterloading systems to the right heart were modeled on the basis of the impedance, rather than steady-state, characteristics of each vasculature. The
systemic venous system was modeled by an upstream pressure source
(Pu) and a four-element venous
impedance network, including a capacitance
(Cv), an inertance
(Lv), and two resistances
(Rvd and
Rvp) (20). The pulmonary
arterial system consists of a three-element windkessel impedance
network, including a capacitance
(Cp) and two resistances
(Rpp and
Rpd), and a downstream pressure
source (Pd). Two one-way valves,
represented by diodes in Fig. 1, were interposed between
Era and
Erv (tricuspid valve) and between
Erv and the pulmonary arteries
(pulmonic valve). A resistance placed between
Era and
Erv represents transtricuspid
valve resistance (Rt). The
driving force for flow in the model is provided by two mechanisms: a
pressure gradient between Pu and
Pd and active periodic increases
in Era and
Erv with the two competent valves.
Fig. 1.
Electrical analog of numerical model of atrioventricular interaction.
See Glossary for definition of
abbreviations. A coupled pericardial constraint was modeled by addition
of a single external elastance
(Epe) over
Era and
Erv; an uncoupled pericardial
constraint was modeled by addition of 2 different external
elastances
(Epera and
Eperv) on
Era and
Erv, respectively.
[View Larger Version of this Image (25K GIF file)]
A coupled pericardial constraint was modeled by adding a single
external elastance (Epe) over
Era and
Erv (Fig. 1). In this situation, a
single charge pressure (Ppe) of
Epe influences
Era and
Erv, simulating a uniform
pericardial liquid pressure that restricts the RA and RV volumes
simultaneously. An uncoupled pericardial constraint was modeled by
adding two different external elastances (Epera and
Eperv) on
Era and
Erv (Fig. 1). The charge pressure of Epera
(Ppera) acts
only on the RA; the charge pressure of
Eperv (Pperv)
influences only the RV. This simulates different regional pericardial
surface pressures that individually restrict the RA and RV volumes.
The behavior of this model can be characterized by sets of differential
state equations as described below. A detailed mathematical description
is provided in the APPENDIX
(Mathematical description of the numerical
model).
With the coupled constraint
|
(1a)
|
|
(1b)
|
|
(1c)
|
|
(1d)
|
|
(1e)
|
|
(1f )
|
With the uncoupled constraint
|
(2a)
|
|
(2b)
|
|
(2c)
|
|
(2d)
|
|
(2e)
|
|
(2f )
|
|
(2g)
|
where
tricuspid flow
(
t) and
pulmonary flow
(
p) are
given by the following equations that define the status of the
tricuspid and pulmonic valves.
When Pra
Prv
|
(3a)
|
and
Pra < Prv
|
(3b)
|
When Prv
Pp
|
(3c)
|
and
Prv < Pp
|
(3d)
|
Numerical solution of the model.
Numerical solutions of the above equations can be performed on a
computer if all the parameters and initial conditions of the state
variables are given. We incorporated experimental data available in the
literature whenever possible to assign specific values for the cardiac
and vascular parameters (Table 1). Details of the justification of these values are included in the
APPENDIX (Parameter
justification and sensitivity analysis for the numerical model). Time-varying properties of
Era and
Erv were modeled by use of cosine
functions (10, 14, 20) to periodically increase from minimum diastolic
values to maximum end-systolic values (Fig. 2).
|
Table 1.
Baseline parameter values used in simulation
|
| Abbreviations |
Values |
|
| Systemic
venous impedance |
Rvd
(mmHg · s · ml 1) |
4
|
Rvp
(mmHg · s · ml 1) |
0.14
|
| Cv (ml/mmHg) |
6.1 |
Lv
(mmHg · s2 · ml 1)
|
0.007 |
| Pulmonary arterial
impedance |
Rpp
(mmHg · s · ml 1) |
0.2
|
Rpd
(mmHg · s · ml 1) |
0.8
|
| Cp (ml/mmHg) |
1.2
|
| Right heart parameters
|
Rt
(mmHg · s · ml 1) |
0.025
|
| Eramax (mmHg/ml) |
1.2
|
| Eramin (mmHg/ml) |
0.6
|
| Tmaxra (s) |
0.12
|
| Ervmax (mmHg/ml) |
4
|
| Ervmin (mmHg/ml) |
0.4
|
| Tmaxrv (s) |
0.18 |
| T
(s) |
0.6 |
| PRI (s) |
0.12 |
|
|
See Glossary for definition of abbreviations.
|
|
Fig. 2.
Time-varying elastances of right atrium and ventricle as functions of
time defined in numerical model.
Era and
Erv were modeled by use of
modified cosine functions.
[View Larger Version of this Image (13K GIF file)]
When 0
tra
2Tmaxra
|
(4a)
|
and
2Tmaxra
tra
T
|
(4b)
|
When 0
trv
2Tmaxrv
|
(4c)
|
and
2Tmaxrv
trv
T
|
(4d)
|
where
tra and
trv represent a
given time after the start of a contraction of the RA and RV until the
next contraction and T is the cardiac
cycle time (reciprocal of heart rate). The heart rate was set at 100 beats/min with a delay time of 0.12 s between the RA and RV contraction
(P-R interval).
Equations 1-4 were programmed and
solved on a Macintosh computer (Apple Computer, Cupertino, CA) by use
of a general purpose simulation software (Extend, Imagine That, San
Jose, CA). The modified Euler method was used for integration with a
simulation step of 0.002 s (500 Hz). On each step, the valve status was
checked and modified according to Eq.
3. A simulation run of at least 45,000 steps (90 s) was
performed to achieve a steady state.
Simulation protocols.
The control condition with no pericardial constraint was simulated by
assigning an effectively nil value (e.g., 0.001 mmHg/ml) for the
pericardial elastances. The initial values for the state variables and
the values of Pu and
Pd were adjusted to best
approximate the normal pressure and flow patterns with an end-diastolic
RV pressure of ~5 mmHg, observed in our canine experiments (4, 25,
27). The pericardial elastances
(Epe in coupled constraint; Epera and
Eperv in
uncoupled constraint) were then incrementally increased from the level
comparable to
Eramin or
Ervmin to the level 10 times larger, i.e., from 0.5, 1.0, and 2.5 to 5.0 mmHg/ml. Pu and
Pd were adjusted such that the
end-diastolic RV pressure increased in 2- to 3-mmHg increments along
with increases in the pericardial elastances (Table
2). Changes in all hemodynamic variables
were calculated, including combined vena caval flow (
v),
tricuspid flow
(
t), volumes
of the RA and RV above their unstressed volumes
(Vra and
Vrv), and total right heart
volume above its unstressed volume
(Vtotal = Vra + Vrv). Systolic-diastolic distribution of
v was analyzed
by calculating the ratio of diastolic to systolic venous inflow volumes
into the RA (D/S

v ratio), as used in our previous study of tamponade in vivo (4). The
contribution of RA contraction to RV filling was evaluated by the ratio
of the two peaks in
t, i.e., the
ratio of the peak in atrial A wave to that in rapid filling E wave (A/E
t ratio). Systole and diastole were defined with reference to the
atrioventricular volume transfer: diastole as the time when the
tricuspid valve is open (i.e.,
t > 0), and
systole as including the periods of isovolumic contraction and
relaxation. Systole and diastole refer to ventricular events unless
specified as atrial.
|
Table 2.
Changes in D/S  v
ratio and A/E t ratio with increases in
pericardial constraint
|
|
Pu, mmHg |
Pd, mmHg
|
Prv,* mmHg |
Cardiac Output |
D/S
 v Ratio |
A/E
t Ratio |
|
| Control condition (no
elastances) |
50 |
5 |
4.5 |
100 |
0.92 |
0.72 |
| Coupled
constraint |
| Epe = 0.5 mmHg/ml |
40 |
5
|
7.0 |
73 |
0.62 |
0.59 |
| Epe = 1.0 mmHg/ml |
35 |
13 |
10.2 |
57 |
0.22 |
0.58
|
| Epe = 2.5 mmHg/ml |
30 |
14 |
12.5 |
39
|
0.02 |
0.52 |
| Epe = 5.0 mmHg/ml |
27
|
17 |
14.3 |
28 |
0.04 |
0.55 |
| Uncoupled constraint
|
| Epera = Eperv = 0.5 mmHg/ml |
40 |
5 |
6.5 |
75 |
1.44
|
0.56
|
| Epera = Eperv = 1.0 mmHg/ml |
35 |
13 |
9.6 |
58 |
1.53
|
0.49
|
| Epera = Eperv = 2.5 mmHg/ml |
30 |
14 |
13.0 |
40 |
3.51
|
0.38
|
| Epera = Eperv = 5.0 mmHg/ml |
27 |
17 |
16.4 |
25 |
11.15
|
0.07 |
|
D/S  v ratio, ratio of diastolic
to systolic venous inflow volumes into RA; A/E
t ratio, ratio of peak in A wave to peak in E wave in t.
*
Values at
end-diastole.
Expressed as relative values (control
value = 100%). With increases in pericardial elastances, values
for Pu and Pd were adjusted by approximately
equal increments to produce desired changes in RV end-diastolic
pressure.
|
|
The model contains a large number of parameters, which may
substantially influence the simulation results. We therefore performed an extensive analysis regarding the relative impact of changes in
cardiac and vascular parameter values (other than pericardial elastances) on the simulation results, particularly the steady-state venous pressure and flow waveforms. Details of the sensitivity analysis
for the parameters are included in the
APPENDIX (Parameter justification and sensitivity analysis for the numerical
model).
Analytic Approach: Model of Ventricular Interdependence
Horizontal cardiac chamber interactions can be considered as
pressure-mediated phenomena, because direct flow interaction is not
present between the two side-by-side chambers. Therefore, we
constructed a simple analytic model of ventricular interdependence in
which the effects of coupled vs. uncoupled pericardial constraint on
ventricular interdependence can be directly assessed by symbolic mathematical analysis. The results were extrapolated to interpret manifestations of the two classical respiratory-induced hemodynamic signs in pericardial diseases.
Model description.
The model was based on a volume elastance model of ventricular
interdependence described by Maughan et al. (15) (Fig.
3). The left ventricle (LV) and RV were
assumed to consist of three volume elastances: LV free wall
(Elvf), RV free wall
(Ervf), and septum
(Es). Thus the volume of the LV
or RV (Vlv or
Vrv) has free wall
(Vlvf or
Vrvf) and septal
(Vs) components.
Vlvf and Vrvf are defined as the volumes of
the LV and RV when the transseptal pressure, i.e., LV pressure
(Plv) minus RV pressure
(Prv), is zero with the septum
in an unstressed neutral position.
Vs is defined as the volume
contribution of the septal shift to either ventricle when the
transseptal pressure is not zero. In this model, all
myocardium-mediated interactions, including transseptal as well
as transcommon fiber interactions, are conceptually
integrated into the interaction through the septal elastance
(15). Coupled and uncoupled pericardial constraints were modeled as
additional volume elastances (Fig. 3), similar to the model of
atrioventricular interaction. With a coupled constraint,
Elvf and
Ervf shared a single external
elastance (Epe) and a uniform
pericardial liquid pressure
(Ppe). With an uncoupled
constraint, two different external elastances
(Epelv and
Eperv) and
regional pericardial surface pressures
(Ppelv and
Pperv) were
added over Elvf and
Ervf independently.
Fig. 3.
Schematic illustration of analytic model of ventricular
interdependence. With a coupled pericardial constraint, both ventricles were assumed to share a single external elastance
(Epe) and a uniform pericardial
liquid pressure (Ppe). With an
uncoupled constraint, 2 different external elastances
(Epelv and
Eperv) and
regional pericardial surface pressures
(Ppelv and
Pperv) were
added over each ventricle independently.
[View Larger Version of this Image (33K GIF file)]
Analytic solution of the model.
By analytically solving the model, it is possible to directly
characterize the "status" of ventricular interdependence as functions of ventricular and pericardial elastances. Details of the
mathematical solution are included in the
APPENDIX
(Mathematical description of the analytic
model). Briefly, the following three interdependence
parameters were derived: right-to-left volume interdependence gain
(GV)
|
(5a)
|
right-to-left
pressure interdependence gain
(GP)
|
(5b)
|
and
effective RV elastance with interdependence
(Erveff)
|
(5c)
|
GV
and GP represent the degree of
ventricular interdependence in right-to-left direction, equivalent to
"cross-talk gains" defined by Maughan et al. (15). As
GV or
GP becomes larger, a given
increase in right heart volume or pressure will produce a greater
increase in left heart pressure.
Erveff is
different from the conventional simple RV elastance, because it
includes influences of the left heart. This index should be viewed as
the effective elastance of the combined RV and LV seen from the
systemic venous port, under conditions when ventricular interdependence is present. As
Erveff becomes
larger, a given increase in right heart volume will produce a greater
increase in right heart pressure. These interdependence parameters are
useful to theoretically predict how frequently pulsus paradoxus or
Kussmaul's sign may be manifest with increases in coupled vs.
uncoupled constraint (see Pericardial Constraint and
Respiratory-Induced Hemodynamic Signs).
Protocols for quantitative assessment.
To quantitatively compare the effects of coupled with uncoupled
constraint on these interdependence parameters, a set of values was
first assigned for ventricular elastances to fix the influences of
myocardium-mediated interdependence. Pericardial elastances (Epe in coupled constraint;
Epelv and
Eperv in
uncoupled constraint) were then incrementally increased, and changes in
the interdependence parameters were numerically calculated and compared
between the two constraint conditions.
Elvf and
Ervf were set at 0.5 mmHg/ml on
the basis of the studies with an isolated canine heart preparation (14,
20). Maughan et al. (15) showed that
Es is 7-15 times larger than the free wall elastances, and relative stiffness of the septum to free
walls plays an important role in determining the cross-talk gains.
Therefore, baseline Es was set at
5.0 mmHg/ml, and situations with lesser
(Es = 2.5 mmHg/ml) or greater
(Es = 10 mmHg/ml) values of
Es were also studied. Under each
of the three situations, Epe or
Epelv and
Eperv were
increased, with the rate of increase adjusted for appropriate
comparison between the coupled and uncoupled constraint conditions.
Epe represents the behavior of the
total pericardium, whereas
Epelv or
Eperv reflects
only each regional portion of the pericardium. If the elastance of the
total pericardium surrounding both ventricles
(Epetotal)
were the same in the two constraint conditions, we would assume that
the elastic constraining capability of the pericardium would be similar
but the manner of constraint would be different. Thus
Epe or
Epelv and
Eperv were
increased to achieve the same value of
Epetotal.
Because
Epetotal can be given as
|
(6a)
|
|
(6b)
|
Epelv
and Eperv were
increased two times more than Epe,
with the assumption that
Epelv and
Eperv equal
each other in Eq. 6b.
RESULTS
Results of Numerical Approach
Under the control condition with no pericardial constraint, the
simulated flows and pressures well approximated the characteristics of
normal venous flows and pressures (4, 25, 27) (Fig. 4A).
v exhibited a
biphasic pattern, with the systolic component larger than the diastolic
component.
t
showed a large rapid-filling E wave with a smaller atrial A wave.
Pra had the systolic
x- and diastolic
y-descents, and
Prv showed little change during
diastole. The simulated Vra,
Vrv, and
Vtotal were also consistent with
expected changes in those volumes during a cardiac cycle (Fig.
5A).
Vra increased at middiastole,
corresponding with the diastolic component of
v.
Vrv showed an increase during
diastole, corresponding with the A wave in
t.
Fig. 4.
Typical simulation traces in numerical model showing changes in
v,
t,
Pra, and
Prv with increases in pericardial
constraint. A: control condition with
no pericardial elastance; B: increased coupled constraint; C: increased
uncoupled constraint. S, ventricular systole; D, ventricular diastole;
E, rapid filling wave in tricuspid flow; A, atrial contraction
component in tricuspid flow; x,
x-descent in right atrial pressure;
y,
y-descent in right atrial pressure.
[View Larger Version of this Image (26K GIF file)]
Fig. 5.
Typical simulation traces in numerical model showing changes in
Vra,
Vrv, and
Vtotal during a cardiac cycle.
A: control condition with no
pericardial constraint; B: increased
coupled constraint; C: increased
uncoupled constraint. S, ventricular systole; D, ventricular
diastole.
[View Larger Version of this Image (24K GIF file)]
With an increased coupled constraint (Figs.
4B and
5B),
v became a
predominantly systolic flow with only a small diastolic component.
Pra showed a prominent
x-descent with a small
y-descent. Consistent with the change
in
v,
Vra and
Vrv changed almost reciprocally
during a cardiac cycle. Vrv
retained its biphasic pattern during diastole with an increasing phase
due to the RA contraction. Vtotal
was minimum at midsystole but almost constant from end systole to end
diastole. In contrast, with an increased uncoupled constraint (Figs.
4C and
5C),
v became a
predominantly diastolic flow. The A wave in
t was markedly
attenuated. The y-descent in
Pra was prominent with a small
x-descent, and
Prv exhibited a steep transient
decrease at early diastole followed by a relatively unchanged portion
in late diastole, i.e., a dip-and-plateau pattern.
Vra showed only a small change
during a cardiac cycle. Vrv became
monophasic during diastole with a minimal increase due to the RA
contraction.
Figure 6 illustrates relationships of
Prv, transmural
Prv (relative to pericardial
pressure), and pericardial pressure
(Ppe or
Pperv) to the
increased coupled or uncoupled constraint. Under both conditions, most
of Prv during diastole is
attributable to the pericardial pressure, with a small contribution
from the transmural Prv. With the
increased coupled constraint, Ppe
showed a transient decrease at midsystole but had returned to its
end-diastolic level at end systole, i.e., a pattern similar to the
changes seen in Vtotal. With the
increased uncoupled constraint,
Pperv showed its maximum at end diastole and minimum at end systole with a waveform
similar to Vrv. Only the uncoupled
constraint produced a "dip" in
Prv in early diastole.
Fig. 6.
Simulation traces in numerical model showing relationships among
Prv, transmural
Prv, and pericardial pressure
during a cardiac cycle. S, ventricular systole; D, ventricular
diastole.
[View Larger Version of this Image (28K GIF file)]
Table 2 summarizes changes in the D/S

v
and A/E
t
ratios with increases in pericardial constraint. Under the control
condition, the D/S

v ratio was 0.92 with an A/E
t ratio of 0.72. With increases in coupled constraint, D/S

v
ratio decreased; i.e., the systolic flow became dominant. The A/E
t ratio showed a
small decrease but remained >0.5. With an
Epe of >2.5 mmHg/ml, the D/S

v
ratio was nearly zero, indicating almost absent diastolic flow in
v. With
increases in uncoupled constraint, the D/S

v
ratio increased and the A/E
t ratio
decreased. With
Epera and
Eperv of 5.0 mmHg/ml, the D/S

v
ratio exceeded 10, indicating that the systolic component was less than
one-tenth of the diastolic component in
v. The A/E
t ratio was only
0.07, consistent with the markedly attenuated A wave in
t.
Results of Analytic Approach
The analytic solution of the model yielded the following formulas for
the interdependence parameters (see
APPENDIX,
Mathematical description of the analytic
model).
With no pericardial constraint
|
(7a)
|
|
(7b)
|
|
(7c)
|
With the coupled
constraint
|
(8a)
|
|
(8b)
|
|
(8c)
|
With the uncoupled constraint
|
(9a)
|
|
(9b)
|
|
(9c)
|
Thus,
GV,
GP, and
Erveff would
increase as the pericardial elastances increase with the coupled or
uncoupled constraint. However, as evident in Eqs.
8 and 9, the manner in
which the pericardial elastances affect these parameters was quite
different between the two conditions.
Epe increases
GV,
GP, and
Erveff in a
totally different fashion from any of ventricular elastances, whereas Epelv and
Eperv increase
them in a manner essentially similar to the free wall elastances
(Elvf and
Ervf).
Figure 7 illustrates quantitative
differences in the interdependence parameters between the coupled and
uncoupled constraints. At a given level of
Epetotal,
GV and
GP were higher with the coupled than with the uncoupled constraint (Fig. 7,
A and
B). On the other hand,
Erveff was
higher with the uncoupled than with the coupled constraint (Fig.
7C). The differences in
GV,
GP, and
Erveff between the two constraint conditions became larger as
Es increased.
Fig. 7.
Quantitative assessment of changes in interdependence parameters with
increases in
Epetotal.
A: changes in
GV;
B: changes in
GP;
C: changes in
Erveff with
interdependence. Differences in
GV,
GP, and
Erveff between
coupled and uncoupled constraint conditions became larger as
Es increased from 2.5 to 5 to 10 mmHg/ml.
[View Larger Version of this Image (39K GIF file)]
DISCUSSION
On the basis of the novel concept of "coupled vs. uncoupled"
pericardial constraint, we evaluated the effects of pericardial constraint on cardiac chamber interactions by use of mathematical model
analyses. The numerical model of atrioventricular interaction well
approximated the steady-state venous flow and pressure waveforms observed in pericardial diseases. Increased coupled constraint accounted for the patterns in cardiac tamponade, and increased uncoupled constraint accounted for those in constrictive pericarditis. On the other hand, the analytic model enabled quantitative comparisons of the status of ventricular interdependence between the two constraint conditions. Increased coupled constraint (tamponade) produced greater
interdependence gains, which should lead to manifestation of a pulsus
paradoxus (17), whereas increased uncoupled constraint (constriction)
was associated with a greater effective RV elastance and, hence, the
increased likelihood of a Kussmaul's sign (25). These findings provide
a basis for the pathogenesis of the characteristic steady-state and
respiratory-induced hemodynamic signs seen in tamponade and
constriction. Thus the construct of coupled vs. uncoupled pericardial
constraint may offer a useful conceptual framework to understand the
pathophysiology in various forms of pericardial diseases.
Pericardial Constraint and Atrioventricular Interaction
The numerical model demonstrated that, with increases in coupled
constraint, the normal biphasic patterns in
v and
Pra were replaced by a
predominantly systolic
v with a
prominent x-descent in
Pra. As the degree of coupled
constraint increased, the ratio of diastolic to systolic venous inflow
volumes to the RA (D/S 
v
ratio) was reduced, such that most of venous return occurred during
systole. These changes are consistent with the characteristic flow and
pressure patterns in cardiac tamponade observed by us (4) and others
(3, 5, 8, 11). With increases in uncoupled constraint, the simulated
v became mainly
diastolic with a prominent y-descent
in Pra, consistent with the
reported findings of constrictive pericarditis (5, 8, 9, 11).
A question then arises as to why such different flow and pressure
patterns take place with the increased coupled (tamponade) or uncoupled
constraint (constrictive pericarditis). It is important to appreciate
that, as the degree of tamponade or constriction increases, the
pericardial pressure increases and approaches the intraluminal
Pra with a decrease in the
transmural Pra (18, 19, 23). In
severe forms of both pericardial diseases, the intraluminal
Pra, i.e., the downstream pressure
for systemic venous return, will be determined almost entirely by the
pericardial pressure, which is the function of pericardial elastic
recoil with minimal if any contribution of the intrinsic time-varying elastic properties of the underlying cardiac chambers.
With the increased coupled constraint, "elevated" and
"uniform" Ppe all over the
heart restricts the sum of the volumes of the RA and RV.
Ppe is a function of the total
right heart volume (Vtotal).
During systole, ventricular ejection moves blood out of the
intrapericardial space, reducing
Vtotal and
Ppe. The decrease in
Ppe will lead to a decrease in the
intraluminal Pra
(x-descent), enhancing venous return
and producing a systolic antegrade
v. If sufficient
volume can move into the RA during systole to replace the ejected
stroke volume, the decreased Ppe
and Pra will again increase to the
level of the upstream venous pressure, and
v should become
zero at end systole. During diastole, the RA emptying and RV filling
via the tricuspid valve do not directly influence Vtotal.
Ppe is not affected by any
intrapericardial volume transfer with the coupled constraint. Thus
Ppe and
Pra will be almost unchanged during diastole, producing an almost absent diastolic
v and diminished
y-descent in
Pra.
With the increased uncoupled constraint, elevated and
"regional"
Ppera and
Pperv
individually restrict the RA and RV.
Ppera is a
function of Vra, whereas
Pperv is
determined by Vrv. During systole, ventricular ejection decreases
Vrv and
Pperv but does not affect
Ppera. Thus the
intraluminal Pra shows little
change during systole (diminished
x-descent) and the systolic
v is small. During diastole, the RA emptying (the decrease in
Vra) into the RV produces a
decrease in
Ppera. The
resultant decrease in the intraluminal
Pra
(y-descent) will enhance venous return, producing a diastolic antegrade
v.
The simulation produced an early diastolic dip in
Prv only with the uncoupled
constraint, but not with the coupled constraint (Fig. 5). This finding
is consistent with the clinical observations that a dip-and-plateau
pattern or a square root sign in
Prv is a characteristic phenomenon
of constrictive pericarditis (18, 19, 23). An insight as to the
pathogenesis of this sign may be derived from the fact that the
pericardial pressure over the RV is a function of
Vtotal with the coupled
constraint, whereas it is determined only by
Vrv with the uncoupled constraint.
In tamponade, Ppe decreases
initially during systole but returns to the previous end-diastolic
level at end systole, along with the corresponding change in
Vtotal. In constrictive
pericarditis, however,
Pperv continues
to decrease during systole and reaches its minimum value at end
systole, along with the change in
Vrv. When the tricuspid valve is
opened, the decreased
Pperv will
start to contribute to generation of a pressure gradient for
venous return. In early diastole,
Pperv increases
as the RV fills, whereas the transmural
Prv decreases because the RV
continues to relax. As a result of these changes, the intraluminal
Prv will exhibit a dip in early
diastole. Because the combined elastance of the RV and pericardium is
large, this dip in Prv will soon
be lost as the RV continues to fill rapidly and then abruptly cease
filling, resulting in a plateau in
Prv in late diastole.
Pericardial Constraint and Atrial Function
Atrial function in relation to ventricular filling has been
conceptually classified in three forms: booster, reservoir, and conduit
functions (10). The atrium boosts ventricular filling during atrial
contraction (atrial kick), acts as a compliant chamber to pool blood
during ventricular systole and supply the ventricle with this blood
during ventricular diastole, and serves as a low-resistance conduit
between the peripheral venous system and the ventricle during
ventricular diastole. The simulation results, particularly the changes
in cardiac chamber volumes (Fig. 5), suggest that these atrial
functions are substantially modulated by the increased pericardial
constraint, even if the atrial elastic properties per se are not
changed.
With the increased coupled constraint, the reciprocal changes in
Vra and
Vrv during a cardiac cycle imply
that the RA fills with a volume nearly equal to stroke volume during
the RV ejection and empties it into the RV during the succeeding
diastole. Because
v is almost nil
during diastole, direct filling of the RV from the venous system should
be minimal during diastole. Thus, with cardiac tamponade, the atrium
would not function well as a passive conduit for venous return during
diastole, whereas it would serve as an efficient reservoir for
ventricular filling by utilizing a decrease in
Vtotal during ventricular
ejection. Atrial function as a booster can still be preserved, as
evidenced by the presence of the A wave in
t and the late
diastolic increase in Vrv.
With the increased uncoupled constraint,
Vra exhibited only small changes
during a cardiac cycle. Venous return occurred mainly during diastole,
with the RA serving effectively as a passive conduit for blood flow
from the venous system to the RV. The decreased A wave in
t and monophasic
diastolic increase in Vrv suggest that the contribution of atrial contraction to ventricular filling was
diminished. Because the difference between the maximum and minimum
elastances of the RA is small compared with that of the RV, a large
external elastance
(Epera) added
over the RA would produce a combined RA-pericardium elastance that
varies minimally during a cardiac cycle. Any influence of atrial
contraction and relaxation on systolic venous flow would be masked by
the increased uncoupled constraint. Filling of the right heart from the
venous system would occur only during diastole when the tricuspid valve is opened and the elastances of the RV are connected parallel to the
elastances of the RA. Thus, in constrictive pericarditis, the atrium
appears to serve mainly as a conduit, and atrial function as a booster
or a reservoir would be markedly attenuated, despite an unchanged
atrial contractility.
Our results suggest that echocardiographic observation of
atrioventricular volume changes within a cardiac cycle may provide an
alternative diagnostic strategy for early detection of tamponade or
constrictive physiology. Marked coupling between atrial and ventricular
volume changes (i.e., reciprocal changes within a cardiac cycle) in
addition to dominant systolic venous flow patterns should be useful
signs of cardiac tamponade, whereas minimal changes in atrial volume
and predominant diastolic venous flow should suggest a constrictive
pericarditis.
Pericardial Constraint and Ventricular Interdependence
The analytic model characterized the effects of pericardial constraint
on ventricular interdependence by defining three interdependence parameters as functions of pericardial elastances. The results demonstrated that the coupled and uncoupled constraints enhanced the
interdependence gains (GV and
GP) as well as the effective RV
elastance
(Erveff), but
the degree of enhancement was quantitatively different between the two
conditions. At a given level of total pericardial constraint,
GV and
GP increased more with the coupled than with the uncoupled constraint, whereas
Erveff
increased more with the uncoupled than with the coupled
constraint.
Our model highlighted similarities and differences among the
interdependence mechanisms mediated by the myocardium per se, coupled pericardium, and uncoupled pericardium. Under conditions with
no pericardial constraint, Eq. 7 clarified how the myocardial factors modulate the interdependence
parameters. First, GV and GP were enhanced by
Es and
Elvf and
Ervf. Although the septal displacement is the sole element to produce myocardium-mediated interactions in the model, the free walls are still able to enhance them by augmenting the transseptal interaction (15). When the septum is
shifted leftward, a stiffer LV free wall will augment the degree of
constraint of a given LV volume, creating a greater increase in LV
pressure. Second,
Erveff was
influenced not only by the RV-related elastances
(Ervf and
Es) but also by the LV-related elastance (Elvf). Thus
factors determining
Erveff are not
necessarily limited to the pure RV factors and will also be
influenced by interdependence.
With the increased coupled constraint,
Ppe restricts the LV and RV
volumes together, increasing GV
and GP. Because the RV and LV are
exposed to the coupled pericardial elastance,
Erveff would
also increase. It is important to note that the enhancement of
interdependence produced by the coupled pericardium does not require
involvement of the transseptal interaction. Equation
8 demonstrated that
Epe affected the interdependence
parameters in a fashion totally different and independent from the
septal or free wall elastances. This implies that the coupled
pericardium provides a unique interdependence mechanism in addition to
the myocardium-mediated interactions.
With the increased uncoupled constraint,
Ppelv and
Pperv impose
local constraining forces on the LV and RV. As a result, the effective
stiffness of each ventricle and, hence,
Erveff would increase. More interestingly, despite the uncoupled nature of the
constraint, the intrapericardial volume coupling between the two
ventricles was also enhanced, and
GV and
GP increased. The clue to
understanding this apparently paradoxical finding is given in
Eq. 9, in which
Epelv and
Eperv affected
the interdependence parameters essentially in a manner similar to
Elvf and
Ervf. Thus the uncoupled regional
pericardium would behave as additional free walls, increasing the
effective stiffness of the LV and RV free walls and augmenting the
transseptal interaction. The interdependence mechanism mediated by the
uncoupled pericardium can therefore be considered as an augmentation of
the myocardium-mediated interactions.
The differences in the interdependence mechanism between the coupled
and uncoupled pericardium may be further clarified by calculating
changes in GV,
GP, and
Erveff when
Es approaches infinity in
Eqs. 7-9.
As Es
, with no
pericardial constraint
|
(10a)
|
with
the coupled constraint
|
(10b)
|
and
with the uncoupled constraint
|
(10c)
|
In
this situation, the septum becomes extremely rigid, so that the
myocardium-mediated interactions become negligible, as represented by
zero values of GV and
GP in Eq. 10a. With the coupled constraint, however,
GV and
GP did not become zero, still
being under the influence of Epe.
In contrast, GV and
GP approached zero with the
uncoupled constraint.
Erveff became
equal to Ervf + Epe with the coupled constraint,
whereas it was equal to Ervf + Eperv with the
uncoupled constraint.
Understanding of these equations can provide an intuitive explanation
as to why quantitative differences were observed in the degree of
interdependence between the coupled and uncoupled constraint conditions
(Fig. 7). As the degree of the pericardial constraint increases, the
coupled pericardium produces a complete intrapericardial volume
coupling between the two ventricles, whereas the uncoupled pericardium
only provides a partial volume coupling by enhancing the already
present transseptal interaction.
GV and GP should therefore be larger with
the coupled than with the uncoupled constraint. However, the RV is
connected to the elastance of the total pericardium
(Epe) with the coupled
constraint, whereas it is connected only to the elastance of the
right-sided pericardium (Eperv) with
the uncoupled constraint. The connection of the RV to
the other half of the uncoupled pericardium
(Epelv) is
indirect, mediated by the transseptal interaction.
Eperv is
substantially larger than Epe,
because it reflects only a regional portion of the pericardium. Thus it
is likely that
Erveff
increases more with the uncoupled than with the coupled constraint at a
given level of pericardial stiffness. The differences between the two constraint conditions should increase as the influence of
transseptal interaction is attenuated (i.e.,
Es increases), as graphically shown in Fig. 7.
Pericardial Constraint and Respiratory-Induced Hemodynamic Signs
The interdependence parameters defined in the model can also be used to
predict the likelihood that the respiratory-induced hemodynamic signs
will be manifest under the coupled or uncoupled constraint conditions.
As the interdependence gains (GV
and GP) increase, the
inspiratory increase in RV volume or pressure would produce a greater
rise in transmural LV diastolic pressure (relative to pleural pressure)
at a given LV volume. In other words, the rise in effective LV
diastolic elastance during inspiration would be greater with higher
values of GV and
GP. This results in larger decreases in LV filling from the pulmonary circulation, thereby increasing the degree or likelihood of pulsus paradoxus (7, 17). Thus
our findings provide a basis for why a pulsus paradoxus is manifest not
only in cardiac tamponade but also in constrictive pericarditis and why
it should be observed to a greater extent or more frequently with
tamponade (coupled) than with constriction (uncoupled).
Although Kussmaul's sign is accepted as a useful clinical sign for
pericardial pathology (18, 19), its pathogenesis had not been well
explained. We recently demonstrated in canine experiments that a
Kussmaul's sign should only occur during inspiration with an active
diaphragmatic descent (25). When an inspiratory increase in systemic
venous return is mainly attributed to a large increase in abdominal
pressure, the rise in transmural RA pressure may exceed the fall in
pleural pressure, leading to manifestation of a Kussmaul's sign. With
a greater
Erveff under
conditions of uncoupled pericardial constraint, the inspiratory
increase in right heart volume would produce greater increases in
transmural RV diastolic or RA pressure relative to pleural pressure.
Thus the results are consistent with the classic observation that a Kussmaul's sign is observed relatively frequently in constrictive pericarditis (uncoupled) but rarely in cardiac tamponade (coupled). In
an intuitive sense, with constriction only the elastances of the right
heart and right-sided pericardium accept the enhanced venous return,
whereas with tamponade the left-sided pericardium would also
participate in buffering the effects of the increased right heart
volume, resulting in a decreased likelihood of the manifestation of a
Kussmaul's sign at a similar level of total pericardial constraint.
Critique
The present study uses two separate mathematical models to analyze the
effects of pericardial constraint on cardiac chamber interactions.
Because atrioventricular interaction within a cardiac cycle is a
dynamic flow-mediated phenomenon in a relatively high frequency range,
we constructed an open-loop numerical model of the right heart
circulation based on the actual experimental data of vascular
impedances. On the other hand, ventricular interdependence has been
successfully analyzed by several previous studies as a
pressure-mediated phenomenon without an element of time taken into
account (12, 14, 21). We thus utilized a simple volume elastance model
that can be directly solved by symbolic mathematical manipulation.
It is theoretically possible to model the overall cardiovascular
system, including the four cardiac chambers as well as
systemic/pulmonary arterial and venous beds, as a single closed circuit
on a beat-to-beat basis. In this situation, however, the model must
simultaneously fulfill two difficult and sometimes conflicting
requirements. First, each vascular bed needs to be characterized as
input impedance to simulate instantaneous venous and arterial pressure
waveforms with reasonable precision. Second, the steady-state
characteristics of each vascular bed must also be considered to
simulate the effects of volume redistribution among the vascular beds
in a closed-circuit circulation. Unfortunately, no model has been
successful in perfectly reconciling the beat-to-beat (high-frequency)
and closed-circuit (low-frequency or steady-state) characteristics into
a single mathematical model mainly because of the lack of necessary
experimental data of vascular properties from these two viewpoints.
Thus our approach is appropriate, although not perfect, given the
present level of experimental knowledge of circulatory parameters.
A potential problem may exist in our model because of the assumption of
linear elements. The pressure-volume relationships of the pericardium
and diastolic ventricular elastances are known to be highly nonlinear.
Although it was possible to arbitrarily choose certain nonlinear
mathematical functions for these elements, there are few experimental
data on which we could reliably define such functions under the normal
and pericardial disease conditions. Other limitations in our numerical
model may include lack of inertances of the tricuspid and pulmonic
valves and lack of the effects of descent of the base of the heart
produced by ventricular contraction. All these factors may cause some
problems in simulating accurate flow waveforms across the tricuspid
valve or may underestimate venous flow during early ventricular
systole. Nevertheless, the overall good agreement between the model and
reality suggests that our simple model is valid as a first
approximation, and the principles developed on the basis of the coupled
vs. uncoupled constraint appeared useful in providing a rationale to
interpret hemodynamic events in pericardial diseases.
ACKNOWLEDGEMENTS
The authors thank Drs. W. L. Maughan, W. Mitzner, K. Sunagawa, and
K. Miyasaka for insightful and encouraging discussions during the
development of this work.
FOOTNOTES
This work is supported in part by National Heart, Lung, and Blood
Institute Grant RO-1-39138-04 and Ministry of Health and Welfare (Japan) Grant H8-PK5-02.
Address for reprint requests: M. Takata, Cardiovascular Research
Center, Massachusetts General Hospital
East, 149 13th St., 4th Fl.,
Charlestown, MA 02129.
Received 4 June 1997; accepted in final form 31 July 1997.
APPENDIX
Mathematical description of the numerical model.
The flows across Rvd and
Rpd
(
1 and
2) are given
by
|
(A1a)
|
|
(A1b)
|
The
pressure-flow relationships at the compliances, inertance, and cardiac
chamber elastances in the model can be expressed as
|
(A2a)
|
|
(A2b)
|
|
(A2c)
|
|
(A2d)
|
|
(A2e)
|
where
Pratm and
Prvtm represent
the transmural pressures of the RA and RV, respectively. With
the coupled constraint,
Pratm, Prvtm, and the
pressure-flow relationship at Epe
should be given by
|
(A3a)
|
|
(A3b)
|
|
(A3c)
|
In
contrast, with the uncoupled constraint,
Pratm,
Prvtm, and the
pressure-flow relationships at
Epera and
Eperv will be
|
(A4a)
|
|
(A4b)
|
|
(A4c)
|
|
(A4d)
|
Solving
Eqs. A1-A3 for six state
variables or solving Eqs. A1, A2, and
A4 for seven state variables, we can
derive six sets (coupled constraint) or seven sets (uncoupled
constraint) of linear first-order differential equations, i.e.,
Eq. 1 or
2.
Parameter justification and sensitivity analysis for the numerical
model.
The parameter values used in the numerical model (Table 1) are adjusted
to approximate cardiovascular properties for an ~15-kg dog.
Characterization of the vascular systems as input impedance, particularly the systemic venous impedance, is essential to accurately simulate dynamic venous pressure and flow waveforms in this model. Sagawa et al. (20) measured hydraulic admittance of the systemic veins
in dogs by means of a two-port analysis. To the best of our knowledge,
this is the only study in the literature that experimentally characterized the dynamic properties of the venous system. Thus the parameters for the systemic veins in the model were chosen to
approximate the actual impedance values obtained in their study. The
parameters for the pulmonary arteries were based on a traditional three-element windkessel impedance model taken from the canine experiments by Westerhof et al. (28). Right heart parameters were
derived from the well-established studies using isolated canine hearts:
RA parameters from the study by Lau et al. (10) and RV parameters from
the study by Maughan et al. (14). The ratio of maximum to minimum
elastances was set at 2 for the RA and 10 for the RV. Tricuspid
resistances were estimated in view of a normal pressure gradient
between the RA and RV.
Because the model contains a large number of parameters, we also
performed an extensive analysis regarding the relative impact of
changes in such parameters on the simulation results of interest, i.e.,
steady-state venous pressure and flow waveforms. Table
3 summarizes the changes in the D/S
v ratio as
each cardiac or vascular parameter was changed (decreased or increased
by factors of 2) from its baseline value. It was apparent that changes
in vascular parameters, once the baseline values are appropriately
estimated from the experimental data, had relatively little effect on
the D/S
v ratio.
Changes in Pu and
Pd did not substantially affect
the D/S
v ratio. In
contrast, changes in cardiac parameters, particularly the RA
parameters, had significant effects on the D/S
v ratio.
However, even the twofold changes in the RA elastance parameters did
not produce such large changes in the D/S
v ratio as
observed with changes in pericardial elastances (e.g., D/S
v ratio
changed from 0.92 to
0.04 with
Epe of 5 mmHg/ml to 11.15 with
Epera and
Eperv of 5 mmHg/ml). These results suggest that the intrinsic time-varying cardiac properties, particularly the atrial parameters, are important, but the
effects of coupled vs. uncoupled pericardial constraint, if present and
increased to a substantial level, may dominate in determining the
status of atrioventricular interaction and venous pressure and flow
waveforms.
|
Table 3.
Changes in D/S  v ratio with
changes in each model parameter
|
|
Changes in D/S
 v Ratio
|
| 0.5 × Baseline
|
Baseline |
2 × Baseline |
|
| Systemic venous
impedance |
| Rvd |
0.97 |
0.92 |
0.82
|
| Rvp |
0.86 |
0.92 |
0.93 |
| Cv
|
0.86 |
0.92 |
0.93 |
| Lv |
0.88 |
0.92
|
0.99 |
| Pulmonary arterial impedance |
| Rpp
|
0.96 |
0.92 |
0.80 |
| Rpd |
1.03 |
0.92
|
0.74 |
| Cp |
0.88 |
0.92 |
0.93 |
| Right
heart parameters |
| Rt |
0.95 |
0.92 |
0.87
|
| Eramax |
1.71 |
0.92 |
0.60
|
| Eramin |
0.56 |
0.92 |
2.93
|
| Ervmax |
0.64 |
0.92 |
1.06
|
| Ervmin |
1.01 |
0.92 |
0.83
|
| Upstream and downstream pressure sources |
| Pu
|
0.81 |
0.92 |
0.98 |
| Pd |
0.98 |
0.92
|
0.81 |
| Pu and Pd* |
0.92 |
0.92
|
0.92 |
|
|
*
Pu and Pd were simultaneously
changed.
|
|
Mathematical description of the analytic model.
The behavior of the model can be analyzed by expressing ventricular
pressures (Plv and
Prv) as functions of ventricular
volumes (Vlv and
Vrv).
Vlv and
Vrv can be written as
|
(A5a)
|
|
(A5b)
|
The
pressure-volume relationship at
Elvf,
Ervf, or
Es can be given as
|
(A6a)
|
|
(A6b)
|
|
(A6c)
|
Plvtm
and Prvtm
represent the transmural pressures and
Vlv 0 and
Vrv 0 are the unstressed
volumes of the LV and RV, respectively. Plvtm and
Prvtm under
conditions with no pericardial constraint should be given by
|
(A7a)
|
|
(A7b)
|
with
the coupled constraint by
|
(A8a)
|
|
(A8b)
|
|
(A8c)
|
and
with the uncoupled constraint by
|
(A9a)
|
|
(A9b)
|
|
(A9c)
|
|
(A9d)
|
where Vpe 0 is the unstressed
volume of the coupled total pericardium and
Vpelv 0 and
Vperv 0 are
the unstressed volumes of the uncoupled LV and RV pericardia,
respectively. Solving Eqs. A5-A9
and eliminating Ppe,
Ppelv, and
Pperv, we can
derive two equations expressing
Plv and
Prv as functions of
Vlv and
Vrv in the following format
|
(A10a)
|
|
(A10b)
|
Differentiating these equations with Vrv at a
constant Vlv should yield
|
(A11a)
|
|
(A11b)
|
|
(A11c)
|
Once
the concrete forms of the coefficients of Eq.
A10 are specified for each of the three constraint
conditions, the actual formulas of
GV,
GP, and
Erveff can be
derived as functions of the ventricular and pericardial elastances.
Under conditions with no pericardial constraint, solving
Eqs. A5-A7 and differentiating yield Eq. 7. Similarly,
Eqs. A5, A6, and
A8 for the coupled constraint and
Eqs. A5, A6, and
A9 for the uncoupled constraint yield
Eqs. 8 and 9, respectively.
REFERENCES
| 1.
|
Agostoni, E.
Mechanics of the pleural space.
Physiol. Rev.
52:
57-128,
1972[Free Full Text].
|
| 2.
|
Appleton, C. P.,
L. K. Hatle,
and
R. L. Popp.
Superior vena cava and hepatic vein Doppler echocardiography in healthy adults.
J. Am. Coll. Cardiol.
10:
1032-1039,
1987[Abstract].
|
| 3.
|
Appleton, C. P.,
L. K. Hatle,
and
R. L. Popp.
Cardiac tamponade and pericardial effusion: respiratory variation in transvalvular flow velocities studied by Doppler echocardiography.
J. Am. Coll. Cardiol.
11:
1020-1030,
1988[Abstract].
|
| 4.
|
Beloucif, S.,
M. Takata,
M. Shimada,
and
J. L. Robotham.
Influence of pericardial constraint on atrioventricular interactions.
Am. J. Physiol.
263 ((Heart Circ. Physiol. 32):
H125-H134,
1992[Abstract/Free Full Text].
|
| 5.
|
Borganelli, M.,
and
B. F. Byrd.
Doppler echocardiography in pericardial disease.
In: Cardiology Clinics: Doppler Echocardiography, edited by N. B. Schiller. Philadelphia, PA: Saunders, 1990, p. 333-348.
|
| 6.
|
Bove, A. A.,
and
W. P. Santamore.
Ventricular interdependence.
Prog. Cardiovasc. Dis.
23:
365-388,
1981[Medline].
|
| 7.
|
Brinker, J. A.,
J. L. Weiss,
D. L. Lappe,
J. L. Rabson,
W. R. Summer,
S. Permutt,
and
M. L. Weisfeldt.
Leftward septal displacement during right ventricular loading in man.
Circulation
61:
623-633,
1980.
|
| 8.
|
Byrd, B. F.,
and
R. W. Linden.
Superior vena cava Doppler flow velocity patterns in pericardial disease.
Am. J. Cardiol.
65:
1464-1470,
1990[Medline].
|
| 9.
|
Hatle, L. K.,
C. P. Appleton,
and
R. L. Popp.
Differentiation of constrictive pericarditis and restrictive cardiomyopathy by Doppler echocardiography.
Circulation
79:
357-370,
1989[Abstract/Free Full Text].
|
| 10.
|
Lau, V. K.,
K. Sagawa,
and
H. Suga.
Instantaneous pressure-volume relationship of right atrium during isovolumic contraction in canine heart.
Am. J. Physiol.
236 ((Heart Circ. Physiol. 5):
H672-H679,
1979[Abstract/Free Full Text].
|
| 11.
|
Linden, R.,
and
B. F. Byrd.
Superior vena cava Doppler: a non-invasive method for the diagnosis of pericardial disease.
Int. J. Cardiol.
16:
145-153,
1987[Medline].
|
| 12.
|
Little, W. C.,
F. R. Badke,
and
R. A. O'Rourke.
Effect of right ventricular pressure on the end-diastolic left ventricular pressure-volume relationship before and after chronic right ventricular pressure overload in dogs without pericardia.
Circ. Res.
54:
719-730,
1984[Abstract/Free Full Text].
|
| 13.
|
Maruyama, Y.,
K. Ashikawa,
S. Isoyama,
H. Kanatsuka,
E. Ino-Oka,
and
T. Takishima.
Mechanical interactions between four heart chambers with and without the pericardium in canine hearts.
Circ. Res.
50:
86-100,
1982[Abstract/Free Full Text].
|
| 14.
|
Maughan, W. L.,
A. A. Shoukas,
K. Sagawa,
and
M. L. Weisfeldt.
Instantaneous pressure-volume relationship of the canine right ventricle.
Circ. Res.
44:
309-315,
1979[Free Full Text].
|
| 15.
|
Maughan, W. L.,
K. Sunagawa,
and
K. Sagawa.
Ventricular systolic interdependence: volume elastance model in isolated canine hearts.
Am. J. Physiol.
253 ((Heart Circ. Physiol. 22):
H1381-H1390,
1987[Abstract/Free Full Text].
|
| 16.
|
Permutt, S.,
and
P. Caldini.
Tissue pressures and fluid dynamics of the lungs.
Federation Proc.
35:
1876-1880,
1976[Medline].
|
| 17.
|
Peters, J.,
C. Fraser,
R. S. Stuart,
W. Baumgartner,
and
J. L. Robotham.
Negative intrathoracic pressure independently decreases both left ventricular inflow and outflow.
Am. J. Physiol.
257 ((Heart Circ. Physiol. 26):
H120-H131,
1989[Abstract/Free Full Text].
|
| 18.
|
Reddy, P. S.
Hemodynamics of cardiac tamponade in man.
In: Pericardial Disease, edited by P. S. Reddy,
D. F. Leon,
and J. A. Shaver. New York: Raven, 1982, p. 161-185.
|
| 19.
|
Reddy, P. S.
Hemodynamics of constrictive pericarditis.
In: Pericardial Disease, edited by P. S. Reddy,
D. F. Leon,
and J. A. Shaver. New York: Raven, 1982, p. 275-297.
|
| 20.
|
Sagawa, K.,
W. L. Maughan,
H. Suga,
and
K. Sunagawa.
Cardiac Contraction and the Pressure-Volume Relationship. New York: Oxford University Press, 1988.
|
| 21.
|
Santamore, W. P.,
T. Shaffer,
and
L. Papa.
Theoretical model of ventricular interdependence: pericardial effects.
Am. J. Physiol.
259 ((Heart Circ. Physiol. 28):
H181-H189,
1990[Abstract/Free Full Text].
|
| 22.
|
Shabetai, R.
The pathophysiology of pulsus paradoxus in cardiac tamponade.
In: Pericardial Disease, edited by P. S. Reddy,
D. F. Leon,
and J. A. Shaver. New York: Raven, 1982, p. 215-230.
|
| 23.
|
Shabetai, R.,
N. O. Fowler,
and
W. G. Guntheroth.
The hemodynamics of cardiac tamponade and constrictive pericarditis.
Am. J. Cardiol.
26:
480-489,
1970[Medline].
|
| 24.
|
Smiseth, O. A.,
M. A. Frais,
I. Kingma,
E. R. Smith,
and
J. V. Tyberg.
Assessment of pericardial constraint in dogs.
Circulation
71:
158-164,
1985[Abstract/Free Full Text].
|
| 25.
|
Takata, M.,
S. Beloucif,
M. Shimada,
and
J. L. Robotham.
Superior and inferior vena caval flows during respiration: pathogenesis of Kussmaul's sign.
Am. J. Physiol.
262 ((Heart Circ. Physiol. 31):
H763-H770,
1992[Abstract/Free Full Text].
|
| 26.
|
Takata, M.,
W. Mitzner,
and
J. L. Robotham.
Influence of the pericardium on ventricular loading during respiration.
J. Appl. Physiol.
68:
1640-1650,
1990[Abstract/Free Full Text].
|
| 27.
|
Takata, M.,
and
J. L. Robotham.
Effects of inspiratory diaphragmatic descent on inferior vena caval venous return.
J. Appl. Physiol.
72:
597-607,
1992[Abstract/Free Full Text].
|
| 28.
|
Westerhof, N.,
G. Elzinga,
and
P. Sipkema.
An artificial arterial system for pumping hearts.
J. Appl. Physiol.
31:
776-781,
1971[Free Full Text].
|
0161-7567/97 $5.00
Copyright © 1997 the American Physiological Society