1. Introduction
The aims of cardiovascular imaging and modeling include
diagnosis and comprehension of patho-physiology, and progress
towards improved interventions for cardiovascular disease.
But the human cardiovascular system is well hidden, extremely
complex and in a continual state of movement and change.
To what extent are we able to image, understand and simulate
the reality of our cardiovascular system? Medical images
and physical or computational models, however sophisticated,
are always far from complete. They represent only limited
aspects of the living system, never with complete accuracy,
and always leaving unseen or un-modeled far more than they
represent. Images and models are inevitably selective. If
blood flow of a single organ or limb, let alone the whole
body, were depicted fully, there would be so much detail
of branched, superimposed streams, micro-streams and counter-streams
that the whole would be incomprehensible in a single image.
Anatomical images are always selective, usually in favor
of larger-scale structures, perhaps depicting one slice
only, or, as in angiography, projecting to a plane the 3-dimensional
distribution of contrast that has arrived in certain vascular
branches. Living tissues are permeated through-and-through
with movement and change, across all scales from the whole
organism down to the inconceivably small scales of molecular
and atomic components.
The limitations of particular imaging or modeling approaches
may not be obvious, especially to those who work with the
most sophisticated equipment and software. The aim of this
paper is to recognize limitations of imaging and modeling
relative to the complexity of the living cardiovascular
system, and to consider how living reality may be questioned
and approached.
Imagination informed and appropriately fluent
can move beyond and between the limitations of separate
images or models. But imagination is liable, on the one
hand, to stray into fantasies and illusions. And, on the
other, it tends to be limited by habits, assumptions and
received interpretations. What is needed is open-minded,
informed imagination, continually taking account of reality
through a range approaches. These may include direct study
of tissues and casts, imaging, critical use of published
information and interactive experiments in the limited but
adaptable context of physical or computational models.
2. Methods available
This section summarizes some approaches to cardiovascular
imaging and modeling, and offers hints towards more inclusive
appreciation of the living cardiovascular system.
2.1. Cardiovascular Imaging
The semi-invasive trans-esophageal approach gives improved
ultrasonic access to more posterior parts of the heart.
In recent years 3-dimensional ultrasound techniques have
been developed involving translation or rotation of a 2-D
scanning plane, but access tends to be restricted, as with
conventional echocardiography.
Cardiovascular Magnetic Resonance Imaging
Magnetic resonance (MRI) is relatively expensive and
complex, but has the advantages of non-invasiveness, safety,
unrestricted access, and unrivaled versatility [Bogaert
et al., 2000; Manning and Pennell, 2001].
All imaging techniques (MRI, x-ray, ultrasound, etc.) require
an energy source, interaction of the energy with body tissue,
and reception of energy coming out again to form an image.
In MRI, the energy source is a radio transmitter, the tissue
interaction involves resonance between radio signal and
oscillations, in the magnetic field, of spinning protons
that constitute the nuclei of hydrogen in body water or
fat. The receiver is a radio aerial, picking up signal that
has been re-emitted by energized protons. What is unique
about magnetic resonance compared with other imaging techniques
is the degree to which interactions at tissue level
nuclear magnetic resonance can be controlled and
manipulated by magnetic gradients. Protons are effectively
played energized by radio pulses and
tuned and re-tuned by magnetic gradients. Different imaging
sequences, which consist of sequences of radio pulse and
magnetic gradient switches, result in different contrast
between tissues. They may also be used to encode velocities
of flow in different directions. This is the basis of the
versatility of magnetic resonance. Magnetic gradients are
applied with respect to x, y and z coordinates, and MRI
is inherently 3-dimensional, allowing resolution of voxels
in 3-dimensional space. Imaging is usually carried out in
2-dimensional planes however, with typical in-plane voxel
dimensions of 1 or 2mm, and a slice thickness of 5 mm or
more. Resolution in the time dimension is achieved by gating
at a series of different delays from the R-wave of the ECG.
Sequential images can be played back as a cine loop which
represents an averaged cardiac cycle, typically acquired
over several heartbeats during a single breath-hold. This
can give very clear images of the heart and vessels, with
wide fields of view, but with a slice thickness of 5mm or
more, spatial resolution is not necessarily as good as it
may appear, and is rarely as good as that provided by computed
x-ray tomography.
One of the great strengths of cardiovascular MRI compared
with other modalities is acquisition of velocity data in
any chosen direction, in or through a chosen plane. This
gives the potential to acquire comprehensive
flow data, i.e. all three directional components of velocity
for points distributed in 3-dimensional space and the dimension
of time [Firmin et al., 1993]. At present, however, constraints
of hardware, software and available time usually limit cine
velocity acquisitions to 1 or 2 directions of velocity in
or through the two dimensions of a plane. Several planes
may be imaged to build up information on a volume [Kilner
et al., 1993 and 2000; Walker et al., 1996]. This type of
velocity data can be correlated with cardiovascular boundaries
detected by other magnetic resonance imaging techniques,
and used as a basis for computational modeling of flow through
heart chambers or large vessels [Saber et al., 2001].
For delineation of 3-dimensional vascular geometry, magnetic
resonance angiography can be performed following injection
in a vein of the contrast agent Gadolinium. This gives fairly
good 3-D spatial resolution, but generally without cardiac
gating, during a breath-hold of about 10 seconds.
Magnetic resonance hardware and software continues to be
developed, and there is great potential for further refinement.
MRI is already the most comprehensive and versatile modality
for cardiovascular imaging and flow measurement, both in
clinical diagnosis and as a basis for cardiovascular modeling.
2.2. Cardiovascular Modeling
The greatest challenges facing those who attempt to model
aspects of cardiovascular flow are probably complexity of
form and complexity of compliance of vascular
walls. Vascular compliance is variable, non-linear and non-isotropic.
Vessels have a mixture of circularly, longitudinally and
irregularly orientated fibers of elastin, plus less-elastic
fibers of collagen which help to limit overstretching
of a vessel as wall tension rises relative to pressure as
diameter increases. In addition, arteries have smooth muscle
fibers that modify elasticity in response to autonomic and
other biochemical signals.
Blood vessels have branching, tree-like structures, dividing
and re-dividing until diameters of less than 100th
of a millimeter are reached in the capillary beds. Increase
in numbers of branches peripherally is so great that, in
spite of extreme diminution in size, summed cross sectional
areas increase about 100-fold from great vessels to micro-vessels.
This goes with decrease in velocities of flow from center
to periphery [Caro et al., 1978].
Centrally, flow through cavities of the heart traces sinuous
paths, with changes of direction at atrial, ventricular
and arterial levels, and quasi-helical twists to the outflow
tracts and great arteries. These turns and twists have significance
for dynamics of flow, especially in the exercising state,
when it would be unrealistic to model flow in one chamber
or great vessel in isolation from those immediately up-
and down-stream [Kilner et al., 1997 and 2000]. In the heart
itself, the blood-muscle boundary is indented and fragmented
by trabeculations consisting of separate but interlinked
muscle bundles that may help to optimize combinations of
contractility and compliance of the muscular wall. The myocardium
as a whole is composed of muscle fibres with different,
more or less oblique orientations in different layers.
Variable viscosity of blood, a non-newtonian fluid, also
contributes to the difficulty of accurate modeling, but
perhaps less than do the complexities of geometry, compliance
and contractility when large-scale cardiovascular flows
are under consideration.
Physical Flow Models
Combinations of form, flow, compliance (and contractility)
found in the cardiovascular system cannot be simulated closely
in models. Form alone can be copied through a process of
casting from post-mortem tissues, if available, for example
using wax or silicone rubber [Kilner et al. 1988]. But although
it is relatively easy to cast the forms of blood spaces
as solid objects, it is much more difficult to re-create
the hollow tubes of vascular branches. Segments of hollow
walls may be cast in silicone rubber around a wax or rubber
cast of the lumen, but their compliance depends on properties
of the rubber and its thickness. Attempts may be made to
include fibers, but I do not think that close approximation
to arterial wall geometry and compliance has ever
been achieved in physical models.
Unrealistic compliance would be relatively unimportant
if flows were continuous, but compliance significantly effects
pressure-flow relationships in pulsatile flow.
For a study of flow through a simple pulsatile cavity,
I have used a pneumatically contained elastic membrane in
circuit of otherwise incompliant polythene tubing [de Leval
et al. 1988]. The model was used only to explore the effect
of pulsatility on flows and pressures in this setting, and
conclusions could only be drawn cautiously in relation to
limitations of the model.
Open-channel flow models are worth exploring even though
they can never give close representations of flow in fully
enclosed blood vessels. It is technically easy to manipulate
and visualize patterns of flow in open channels, which can
be made in malleable materials such as wax or clay. A suspension
of fine metal flakes can be used for flow visualization,
and contours of the free surface may give some insight into
differences of pressure beneath the surface. Much can be
learned about wave propagation, flow separation and flow
instability through interactive modeling in open channels.
Surface gravity waves propagate and interact in ways that
loosely simulate pressure wave propagation in elastic vessels
but the approximation is not close. Also, multi-directional
flows in the half-cylinder of an open channel may loosely
simulate those in one half of a tubular vessel, although
the approximation fails where axes of vessels and branches
have non-planar geometry. Although openchannels fail
to simulate vascular dynamics, their suitability for interactive,
easily viewed experiment makes them invaluable aids to exploring
principles of flow dynamics in relation to variable geometry.
Experience gained is relevant to interpretation of appearances
on flow images and design of other types of flow model.
Computational Fluid Dynamic Modeling
I do not have direct experience of computational flow modeling,
but have collaborated with others in the field [Migliavacca
et al., 1997 and 1999; Saber et al., 2001]. My impression
is that computational simulation of vascular geometry is
as challenging as its physical counterparts, and simulation
of interactions between fluid and compliant boundaries remains
extremely difficult by computation. The great advantage
of computation, however, is that methods are readily stored
and transferred as software, and computation lends itself
to interactive manipulation of selected parameters. Parameters
such as velocity, pressure and shear may be mapped across
the volume of virtual flow. A further, increasingly relevant
advantage is that geometrical information available from
magnetic resonance and other imaging techniques can be used
fairly directly for determination of the geometry of containing
boundaries, although accuracy leaves much to be desired,
mainly due to limitations of spatio-temporal resolution.
We should not be deceived by apparent sophistication of
computational fluid dynamics in attempted simulations of
in vivo flow. The inaccuracies of simulated geometry, compliance
and, in the case of cardiac cavities, contractility remains
significant. As with simplified physical models described
above, I see computational modeling as a valuable means
of exploring and separating out principles of flow-structure
interaction that may be relevant to the living system. Modeling
still falls a long way short of simulating the complexity,
flexibility and interactive responsiveness of living cardiovascular
systems. In life dynamics of flow through the heart and
vessels change radically from resting to exercising states
a fact rarely considered in relation to computational
simulations.
2.3. Complementary Approaches to Cardiovascular Structure
We live in times of ever increasing specialisation. For
this reason, it is important to make efforts to broaden
perspectives and explore alternative approaches. I will
outline two kinds of resource, which complement more familiar
approaches to cardiovascular research. Each allows direct
use of the senses tactile and visual on aspects
of cardiovascular structure, but inevitably not in the living
circulatory system. Then I will list some basic methods
that may foster more comprehensive appreciation of the cardiovascular
system.
Fresh Animal Tissues
A number of years ago, when preparing to work with physical
flow models, I bought and handled fresh animal hearts-and-lungs
from a local meat supplier. I also experimented with incorporation
of fresh vessel segments (caval vein, pulmonary artery and
aorta) in simple pulsatile flow models. These tissues are
of little commercial value to meat traders, but unfortunately,
with precautions against possible spread of infective disease,
fresh animal cardiovascular tissues may be less easy to
obtain now.
Handling fresh tissue gives insight into the varied structure
and elastic properties of veins, atrial and ventricular
cavities, heart valves, pulmonary arteries, aorta and arterial
branches. It is striking how elasticity and strength of
each tissue is adapted to the range of pressures and stresses
it supports in vivo. The non-linear compliance of vascular
walls, both longitudinal and circumferential, can be judged
by stretching with the fingertips. The delicate compliance
and collapsibility of caval veins, normally functioning
with low, sometimes negative transmural pressures, contrasts
with the elastic resilience and circular cross section of
the pulmonary artery and aorta.
Cardiovascular Casts
For a number of decades post mortem casts of cardiovascular
blood spaces have been prepared using injectable, solidifying,
extractable materials such as resin [Tompsett 1970; McMinn
and Hutchings, 1985], silicone rubber [Kilner et al., 1988]
or microfil gel for microvascular casts [Williams
and Warwick, 1980]. There is a superb collection of colored
resin casts of vascular trees of human organs and organ
systems in the anatomical museum of the Royal College of
Surgeons in London. Well-prepared three-dimensional casts
give unique insight into the complexity of cardiovascular
form. Different kinds of cast should be considered in relation
to one another, for example arterial and venous, and casts
that exclude and include very small branches. Larger branches
may only be visible when smaller branches have been removed
from a cast, so dense can be the fine branches towards the
microvessels.
Suggestions Towards an Inclusive, Imaginative Approach
to the Cardiovascular System:
- Practice careful, unbiased, inclusive observation.
As far as possible, observation should begin with direct
sensory perception, unaided by instruments. The living
cardiovascular system is largely inaccessible, however.
Fresh post-mortem tissues and vascular casts provide
non-living substitutes, allowing selected aspects of
structure and form to be studied directly.
- Efforts should be made to withhold interpretation
or explanation while gathering more observations and
information on different aspects of the cardiovascular
system.
- It is worth being alert to the observations of others,
who may notice what you do not.
- Where instruments (imaging equipment, blood pressure
monitors, ECG traces, microscopes, etc) are used, it
is important to be aware of their mode of action, their
limitations, and the ways in which they select or alter
appearances. A microscope or microscopic image, for
example, gives false impressions of massiveness and
distance, which should be countered by efforts to appreciate
the actual smallness of what is seen.
- Attention should be paid to the wider context (organism,
environment) as well as to the details that contribute
to phenomena.
- Attention should be paid to processes and change in
the dimension of time [Hildebrandt 1991; Hidebrandt
et al., 1998], as well as in space.
- Do not necessarily trust what you read or are taught.
Do not even trust what you see, which tends to be biased
by what you expect to see. Repeatedly question whether
another point of view or interpretation is also feasible.
- As you enrich your experience through several observational
and imaging approaches, allow your imagination
informed and appropriately fluid to move among
the phenomena. New kinds of understanding may arise
which can show up the limitations of previously held
interpretations.
- Creative engagement is a potent way of testing understanding.
This is where modeling physical or computational
has its place. Relationships and influences of
one variable on another can be explored, within the
acknowledged limits of a model. Understanding can applied
back to in vivo information, and further investigations
or experiments devised.
- Another kind of creative engagement is drawing, from
image data, the forms and flow paths of the cardiovascular
system. [Kilner et al., 1993; Kilner et al., 2000].
3. Result a Brief Description of the Living
Cardiovascular System
Such a richly varied fluid system cannot be captured in
words or images alone appropriate fluency of thought
is needed to begin to appreciate the dynamic beauty of the
circulatory system.
As you read this, your blood is streaming inwards and outwards
- converging in veins from all parts of your body, turning
through the entwined but separate flow paths of right and
left heart, and branching out again through arteries to
the microscopically fine capillary webs of the lungs and
other organs.
Observation of vascular and micro-vascular casts allows
appreciation of the division and sub-division of vascular
branches, down to the delicate webs of capillaries that
permeate most of the volume of the body. These filament-like
micro-vessels have diameters approaching one tenth of those
of a human hair. They accommodate red blood cells (each
less than one hundredth of a millimeter in diameter) only
in single-file, and with compliant deformation. Capillary
webs in different organs (lung, kidney, liver, muscle etc.)
have varied spatial forms. Capillaries in the lungs, for
example, have the form of a net so dense that there is more
net than hole i.e. a thin membrane of blood effectively
seeps round each alveolus, its bounding endothelial surfaces
linked by tiny columns of endothelium that form the holes
of the net. Capillaries of a renal glomerulus are also densely
packed in a tuft, whereas those of the renal tubules extend
in linear filaments [Williams and Warwick; 1980]. Each organ
or tissue has characteristic capillary morphologies, which
are associated with different metabolic and functional roles.
They vary in compliance, resistance to flow and responsiveness
to biochemical change. The pulmonary vascular bed, for example,
generally has about one fifth of the resistance of all the
systemic vascular beds combined. This is reflected in low
pulmonary artery pressures relative to those of systemic
arteries, despite similar outputs from right and left ventricle.
Each of many billions of capillaries of the body is in
fluid continuity, via arterial branches upstream, and venous
branches downstream, with the entwined streams of the left
and right heart. Arterial and venous branches of each organ
or limb usually run side-by side, so that venous and arterial
trees penetrate the same volume, their branches adjacent
and counter-flowing. Organs may also be penetrated by additional
vascular trees, for example the portal venous and biliary
trees of the liver, the bronchial trees of the lungs, the
tubules of the kidneys, and delicate lymphatic branches
through most tissues of the body.
Components of living systems, more than those of mechanical
or computational systems, are interactive and mutually responsive.
The heart generally thought of rather simply as a
pump propelling blood through vascular resistance
is responsive to and dependent on returning streams in a
number of ways. These include the responsiveness of myocardial
fibres to initial stretch, according to the Frank-Starling
principle [Noble, 1978], and responsiveness to circulating
substances such as adrenaline. Heart rhythm and contractility
are modulated by sympathetic and parasympathetic nerve activity
via various feedback pathways and myocardial contractility
depends on exchanges of nutrients and gases via coronary
blood flow. Flows of left and right ventricles interact
with one another through displacements of the septum, especially
when it is relaxed during ventricular filling.
Thanks to its confluent situation, central to circulatory
branches of the body and lungs, the heart unifies the diversity
of our organism. Blood filling the right atrium with each
beat includes cells that have passed through capillary beds
in all parts of the body. Similarly, blood ejected from
the left ventricle is destined for arterial and microvascular
branches throughout the volume of the body. Different individual
blood cells will take different paths and different periods
of time (from about 0.3 to about 6 seconds) to reach and
pass through capillary beds. Periods of venous return can
be longer and even more variable. This means that the particular
gathering of blood cells in the heart at one moment will
never come together again there will always be thorough
mixing and re-mixing through circulatory dispersions and
confluences.
A single red cell, with a life span of about 3 months,
will trace more than 100,000 systemic and pulmonary journeys,
probably never by exactly the same route twice. It will
gradually trace paths to and from the heart, to and from
all parts of the body and lungs. The circulatory system
is a system of continual transport, mixing and exchange.
In summary, blood circulation and heart maintain unity
in diversity, and continuity through continual change. The
heart is responsive and active, sounding and reliably serving
the diversity of our organism. The seeping of blood through
billions of varied, microscopically fine capillary branches
and, at the same time, the flows and counter-flows in branches
of arteries and veins, and the sinuous streams swirling
through the curvatures of the heart - these are awesome
in their complexity, unity and beauty of form. This fluent
masterwork moves in each one of us, and continues to move,
reliably, but ever changing, day and night, life-long.
Acknowledgements
Approaches to study outlined in this paper owe much to
the natural scientific methods of Johann Wolfgang von Goethe
(1749-1832) [Bortoft 1996; Seamon and Zajonc, 1998]. These
methods were reviewed and developed by Rudolf Steiner (1861-1925)
[Steiner 1968 and 1979], and have since been applied in
biological fields by scientists including Wolfgang Schad
and Jochen Bockemuhl, whose work and teaching I have appreciated
greatly.
References
Bogaert J, Duerinckx AJ, Rademaker FE. . Magnetic Resonance
of the Heart and Great Vessels. Springer Verlag,
Berlin. 2000
Bortoft H. . The Wholeness of Nature Goethes
Way Toward a Science of Conscious Participation in Nature.
Lindisfarne Press, Hudson, NY.1996
Caro CG, Pedley TJ, Schroter RC. . The mechanics of the
circulation. Oxford: Oxford University Press.1978
de Leval M, Kilner PJ, Gewillig M, Bull K. Total cavo-pulmonary
connection: a logical alternative to atrio-pulmonary connection
for complex Fontan patients. J Thorac Cardiovasc Surg,
96: 682-695.1988
Firmin DN, Gatehouse PD, Konrad JP, Yang GZ, Kilner PJ,
Longmore DB. Rapid 7-dimensional imaging of pulsatile flow.
Computers in Cardiology, IEEE computer Society,
London, 353-356.
Hanrath P(Ed). Cardiovascular imaging by utrasound. Kluwer
press.1992
Hildebrandt G, Moser M, Lehofer M. . Chronobiologie und
Chronomedizin. Biologische Rythmen Medizinische Konsequenzen.
Hippokrates Verlag, Stutgart. 1998
Hildebrandt G. . Reactive modifications of the autonomous
time structure in the human organism. J Physiol Pharmacol,
42:5-27.1991
Kilner PJ, Ho SY, Anderson RH. . Cardiovascular cavities
cast in silicone rubber as an adjunct to post-mortem examination
of the heart. Intern J Cardiol, 22: 99-107.1988
Kilner PJ, Yang GZ, Mohiaddin RH, Firmin DN, Longmore DB.
Helical and retrograde secondary flow patterns in the aortic
arch studied by three-directional magnetic resonance velocity
mapping, Circulation. 88:2235-2247.1993
Kilner PJ, Henein MY, Gibson DG. Our tortuous heart in
dynamic mode - an echocardiographic study of mitral flow
and movement in exercising subjects. Heart and Vessels,
12:103-110.1997
Kilner PJ, Yang GZ, Wilkes AJ, Mohiaddin RH, Firmin DN,
Yacoub M. Asymmetric redirection of flow through the heart.
Nature, 404:7759-761.2000
Manning WJ, Pennell DJ. Cardiac Magnetic Resonance
Imaging. Churchill.2001
McMinn RMH, Hutchings RT. A Colour Atlas of Human Anatomy.
Wolfe, London.1985
Migliavacca F, Kilner PJ, Dubini GG, Pennati G, Redaelli
A. Localised peaks of computed flow velocity through asymmetric
stenoses in large vessel models could have implications
for "pressure gradient" estimation. J Vascular
Investigation, 3:17-25.1997
Migliavacca F, Kilner PJ, Pennati G, Dubini G, Pietrabissa
R, Fumero R, de Leval MR. Computational fluid dynamic and
magnetic resonance analyses of flow distribution between
the lungs after total cavopulmonary connection. IEEE
Transactions on Biomedical Engineering., 46(4) 393-399.
1999
Noble MI. The Frank-Starling curve. Clin Sci Mol Med,
54:1-7.1978
Saber NR, Gosman AD, Wood NB, Kilner PJ, Charrier CL, Firmin
DN. Computational flow models of the left ventricle based
on in vivo MRI data: initial experience. Annals of Biomedical
Engineering, 29: 275-283. 2001
Seamon D, Zajonc A (eds). Goethes Way of Science
A Phenomenology of Nature. SUNY Press. 1998
Shapiro LM, Kenny MA. Cardiac ultrasound. Manson press.1998
Steiner R (translated by Wannamaker OD). Theory of Knowledge
implicit in Goethes World Conception. Anthroposophic
Press, New York. 1968
Steiner R (1886) Erkentnistheorie der Goetheschen Weltanschauung.
Rudolf Steiner Verlag, Dornach, Swizerland. 1979.
Tompsett DH. Anatomical Techniques. Livingstone Edinburgh.
1970
Walker PG, Cranney GB, Grimes RY, Delatore J, Rectenwald
J, Pohost GM, et al. Three-dimensional reconstruction of
the flow in a human left heart by using magnetic resonance
phase velocity encoding. Ann Biomed Eng, 24:139-147.
1996
Williams and Warwick (editors). Grays anatomy,
36th edition. Churchill-Livingstone, London,
1400-1401.1980