CC BY
31
УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Том 148, кн. 3
Физико-математические пауки
2006
UDK 517.958:57^519.63
NUMERICAL MODELING OF THE DESIGN OF BIFURCATED PROSTHESES USED IN THE TREATMENT OF ABDOMINAL AORTIC ANEURYSM
S. Lapin, S. Canic
Abstract
The main goal of this paper is to address the numerical solution of the problem of fluid-structure interaction between blood and arterial walls. We use one-dimensional effective model derived asymptotically from the Navier Stokes equations describing blood flow and the linear elastic membrane equation for arterial walls. In this study we consider case with bifurcated compliant arteries with inserted prosthesis. We use the derived model to estimate the influence of a shear stress on the developing of at.herogenesis.
1. Introduction
1.1. Background of the problem. The dynamic interaction between a fluid and a structure is a mechanism that is seen in many different physical phenomena. One of the fluid structure interaction problems receiving a lot of attention in recent years is the problem arising in hemodynamics applications. In this paper we consider numerical simulation of the blood flow through branching arteries with an inserted prostheses called stents and stent-grafts [1]. Both the arteries and the stents are assumed to be elastic. The Young's modulus of elasticity for arteries was obtained from [2]. The Young's modulus of elasticity of stents was obtained from the measurements by Iv. Ravi-Cliandar and R. Wang [3. 4]. The underlying problem consists of modeling the flow through the human aorta branching into iliac arteries, that suffers from an aortic abdominal aneurysm (AAA).
Non-surgical treatment of AAA consists of inserting a prosthesis inside the diseased aorta to redirect the flow of blood and lower the pressure to the aneurysmal walls. There are various complications associated with this procedure, they include stent-graft migration, occlusion of the graft limbs and formation of new aneurysms near the anchoring sites [5]. The purpose of the study presented here is to nse mathematical modeling and numerical simulations to understand some of the hemodynamic factors that might be responsible for the complications listed above, and suggest an improved prosthesis design that might minimize some of the problems.
1.2. Mathematical model. We nse a reduced, one dimensional model, together with the special conditions derived via Riemann invariants, to simulate the flow through the branching arteries. The model is derived using an asymptotic reduction of the Navier Stokes equations modeling blood with the linear elastic membrane equation for the wall.
The resulting set of equations is nsed to study the optimal prostheses design by investigating the impact of shear stress rates on the initiation of focal atherogenesis.
.¿J_____
R+n R
Q
x
Fig. 1. Computational domain
2. Effective model derivation
Wo consider the unsteady axisynimetric flow of a Newtonian incompressible fluid in a thin elastic cylinder whose radius is small with respect to its length. We define the ratio £ = R/L, where R is the radius and L is the length of the cylinder. Now, for every fixed £ > 0 we introduce the computational domain (Figure 1)
QE(t) = {x e R3; x = (r cos©,r sin©,x), r<R + ne(x,t), 0 < x < L},
where ne(x, t) is the radial wall displacement.
This domain is filled with fluid modeled by the incompressible Navior Stokes equations. We assume that the flow is axially symmetric, so the fluid velocity is ve(r, x, t) = (vE(r,x,t),v%(r, x, t)); the pressure defined by pe(r,x,t) — pref with pref being constant reference pressure.
Now the problem in an Eulerian framework in cylindrical coordinates in QE(t) x R+ reads as follows
dvFr ~~dt
dr
- dvr
; dx
M
d2v
dr2
1 dvFr r dr
d2v
dx2
dpe dr
0,
P
dK m
dvE dvE
d2v
dr2
1 dvfv r dr
d2v
dx2
dpe dx
0,
(1)
(2)
^ =0,
(3)
dx
dvi dvt
—L H--^
dr dx r
where m the fluid dynamic viscosity, and p is the fluid density.
Further we assume that the lateral wall of the cylinder £e(t) = x (0, L) is elastic and allows only radial displacement. We can describe its motion (in Lagrangian coordinates) using the linear elastic membrane equation. The radial contact force is given by fC]
{r = R + ne(x,t)}x
h(e)E(e) if i]_ _ d2if 1 -a2 R2 PleiR Pw [t~> dt2
(4)
where Fr is the radial component of external forces (coming from the stresses induced by the fluid), h = h(£) is the membrane thickness, pw is the wall volumetric mass, E = E(£) ^s the Young's modulus, 0 < a < 0.5 is the Poisson ratio. Typical values of the parameters are given in Table 1.
The fluid equations are coupled with the membrane equation through the lateral boundary conditions requiring the continuity of velocity and balance of forces. Now we consider the balance of forces, namely we set the radial contact force equal to the radial component of the force exerted by the fluid.
The fluid contact force is given in Enlerian coordinates as
Ff = ((pe — pref )1 — 2MD(v£))ner
Table 1. Parameter values
Parameters Values
e 0.002 - 0.06
Characteristic radius: Ro 0.0025 - 0.012 m
Characteristic length: L 0.065 - 0.2 m
Young's modulus: E 10b - 10b Pa [2]
Wall thickness: h 1 - 2 x m
Blood density: p 1050 kg/m-"
Reference pressure: pTef 13000 Pa = 97.5 mmHg
where D(ve) is the symmetrized gradient of velocity defined by,
Using the Jacobian of the transformation from Eulerian to Lagrangian coordinates we have (pointwise):
-Fr = ((jf -Plel)I - 2pD(vx))ner(l + + (^j-y- on E° x R+. (5)
The initial data are given by.
']E = In = 0 and l'e = 0 on s<?(0) x (6)
We assume that the end-points of the tube are fixed and that the problem is driven by a time-dependent pressure drop between the inlet and outlet boundary Therefore, we have the following inlet-outlet boundary conditions:
vX =0, px + p(vX)2/2 = Pi(t)+ Pref on (6Qe(t) f {x = 0}) x R+, (7)
vX =0, px + p(vX)2/2 = P2(t)+ Pref on (6Qx(t) f{x = L}) x R+, (8) nx = 0 for x = 0, nx = 0 for x = L and Vt G R+. (9)
We will assume that the pressure drop, A(t) = P1(t) — P2(t) G Co°(0,
2.1. Asymptotic reduction of the model. Now, using the asymptotic techniques described in detail in [7, 8] we will obtain the reduced, two-dimensional equations.
In order to represent the equations (1) (3) in nondimensional form we introduce the independent variables r, x and t
• r = Rr,
• x = Lx
1 ~ , , 1 I hE
• t = —t. where w=—t ——-^ .
ux L\ Rp(l — a2)
Further, we introduce the following asymptotic expansions
Ve = V{i>° + ev1 + •••}, with V=x g2)P, (10)
n 1 , - R2( 1 -E{if + e»7 +•••}, with S = —^-'-v,
hE
pV 2{pn + ep1 + •••}.
(ID (12)
Here P denotes a norm that measures the magnitude of inlet and outlet pressure,
the pressure drop, and the averaged pressure in one cardiac cycle [7].
e2
incompressibility condition become
dV%
1
3V£
' '-<777 ' ' Re
dr
Bp) dr
1 d
dr V dr
d ,„„ , d ,„„ (rVr) + -77T(>'V
(rV:r: ) = 0,
(13)
(14)
(15)
(16)
dr dX
where Sh = (Lw£)/^d Re = (pVR2)/(pL).
The leading order linear elastic membrane equation read
-Fr = P) + O(e2).
The asymptotic form of the contact force derived from (5) becomes
((i/ -Plei)I - 2pD(v£))ner = pV V -Plei + 0(e2)) f 1 +
Therefore, we obtain the following leading order relationship between the pressure and the radial displacement
^ = in ' ifp = pfe^ I1" ^
VR
PV
R
R
) ■■■)■ (17)
We assume that ^/R is of order e so the last term on right-hand side of (17) can be ignored and we obtain
P R
P£ - Pvei = -y^i), (18)
which in dimensional variables gives us the Law of Laplace.
Eh
P - Pref
(1 - a2)R R'
(19)
Having the above results in hand, we proceed by introducing the reduced two-dimensional model written in non-dimensional variables. We define the scaled domain
Cl(i) = \ (x, f) G R2|r < 1 + i), 0 < X < 1
R
and the lateral boundary E(f) = {r = 1 + %t](x,i)} x (0,1).
0
0
Now wo formulate the two-dimensional problem in the following way. Find a triple (Vx,Vr,n) satisfying
d£e Sh—# dt
d{'fv
: dx
-
dr
dp dx
1
Re
ld_ r dr
,dvf
V _
dr
d d
dî-r' + dl = P R
(PE -Pref) = -r^fj,
~ / ~ S ~ / \ dn _ Vr(x, 1 + —?7(x,i),i) =
■with the initial and boundary conditions given by
= ^2=0 at £ = 0,
n =
P : p :
(Pi(£) (P2(£)
dt
Pref )/(pV 2) Pref )/(pV 2)
on (dne(i) n{x on (3Qe(t) n{X
0}) x 1}) x
n = 0 for X = 0, n = 0 for X = 1 and V£ G
(20) (21) (22)
(23)
(24)
(25)
(26) (27)
The model defined above is a free-bonndary degenerate hyperbolic system with parabolic regnlarization. This system, though already simplified, is still quite complex. The theoretical analysis and numerical simulation of (20) (27) is a difficult task: that is why we will proceed further in order to obtain a simplified one-dimensional reduced model.
2.2. Derivation of the one-dimensional reduced model. In order to derive a simplified one-dimensional model we express the two-dimensional equations in terms of the averaged quantities across the cross-sectional area.
Let us introduce A = (1 + and m = AU where
U =
A
xr dr,
(28)
i+S'ï
A U2
2vx2rdr.
(29)
Now we integrate the equations (20) (21) from r = 0 to r = 1 + ^ij and then express them in terms of the averaged quantities. Using the no-slip condition at the lateral boundary we obtain
dA s dm ~dt + R~di
0,
(30)
„, dm d ( am Sh—- + — —
dt dx V A
' dx Re
dvx dr
0
A.A.
Fig. 2. Blood flow through abdominal aorta (A.A) and iliac arteries (I.A.(l) and I.A.(2))
The next step is to specify the axial velocity profile vx in terms of the averaged quantities. The typical approximation (corresponding to Poiseuille flow [9]) is
- 7 + 2,.
vx =-U
1 -
(31)
li].
We assume that 7 = 9, which corresponds to the non-Newtonian flow of blood [10, After plugging (31) into (27) the right hand side becomes
-2(7 + 2) m
R^ X
Using A0 to denote the non-stressed area R2 we obtain a one-dimensional system written in dimensional variables
OA dm
= 0,
~dt 4 dx
dm d f am■
~dt ~ dx A
p dx Eh
A
p := p - pref
V(1 - <r2)R yv Ao
We use a = 0.5, so equation (34) reduces to
A
4--1
(32)
(33)
(34)
p = Go
1
(35)
with G0 = (4Eh)/(3R).
3. Conditions at the bifurcation point
We use the system of equations (32) (33) and (35) to model blood flow through the abdominal aorta branching into the iliac arteries (Fig. 2).
One of the important questions here is how to model the flow at the branching point. We use the continuity of pressure and conservation of mass condition to conple the flow exiting the abdominal aorta with the flow entering the iliac arteries [12]. It was shown in [13] that the 1-D approximation to the flow at the branching point for the
fixed geometry provides e1/2 accuracy to the 3-D flow for small to moderate Reynolds numbers.
We calculate the continuity of pressure and conservation of mass condition at the branching point using the concept of Riemann invariants for a hyperbolic system. More precisely, let us write the system (32) (33) and (35) in qnasilinear form:
At
mt
0
1
m2 1 , ,. 2m
Ax
0
-2-
a — 1
(36)
The eigenvalues of (36) are
m
A
'1 Go (A_ 2 p VA0
1/2
A2
A
II Go (A_ 2 p VA0
1/2
(37)
Let us further use the following notation:
/1 G, 2 ~P
One can show that the right eigenvectors corresponding to A1 and A2 are given by
r 1 =
1
A1
r2 =
1
A2
Consider Riemann invariants w and z corresponding to the eigenvalues A1 and A2 respectively. Then w and z satisfy the following equations [14]:
Vw • rXl =0, Vz • rX2 = 0,
Direct calculations give us
m 11 Go f A
iv = — +4W--- —
A \ 2 p \A0
1/2
A V 2 P \A0
1/2
(39)
(40)
(41)
Now. we assume that the bifurcation occurs at a point and that there is no leakage at the bifurcation. Then the outflow from the abdominal aorta (''parent" vessel) should be balanced by the inflow into the iliac arteries ("daughter" vessels) [12]:
mp = mdi + md2, (42)
which describes conservation of mass. We also assume the continuity of pressure
Pp = Pdi = Pd2, (43)
here subscript p corresponds to the "parent"vessel and dl, d2 correspond to the "daughter" vessels.
The Riemann invariants for the "parent" and the two "daughter" vessels read as follows:
-P = ? + cp^4, (44)
a
A
x
z.i^-c^4, (45)
= (46)
From the equations (35) and (42) (46) it is straightforward to obtain
Cd1
Aodi Aocn
Aop zdi — Aop Zd2
Ami ^ ~ ( AQC12
. A0p y ) + - Cd2 I , A0p .
41/4 = _A0p A0p__U7)
P ( /A \5/4 /, \5/4\ '
Thus, using Riemann invariants and conditions (42) and (43) we can calculate Ap, Ad1, mp, md^d md2 in terms of the Riemann invariants wp, zdi and zd2. The "forward" Riemann invariant wp and the "backward" Riemann invariants zd^d zd2 are known from the initial conditions and the inlet and outlet boundary data.
4. Numerical method
In order to solve equations (32) (33) and (35) we rewrite them first in conservation form. We take into account that Ao can depend on x and write the equations in conservation form as follows
d d
mu + l£F = s> (48)
where
U =
A
F (U ) =
A
G o ¥
A_
M
3/2
An
and
S(U )
-2-
- 1 A
0
Gn ¥
A_
3/2
AO
(49)
(50)
For the numerical solution of the problem stated above we apply the two-step Lax Wendroff method [15]. We assume here that the grid is uniform, with Ax denoting the mesh step size and At the time step, then we define U^ to be the approximation of the solution at (mAx, nAt). The method takes the following form
Un+i= un -
^ <m ^ vri
Ai Ax
F (U^)
- F{U£}<;2j) + % (S{U%\%) + S(U
where
Uj
U n + U n
+ 1/2 _ Lj+1/2 T L j-1/2
Un
F(U
n
3 + 1/2
- F(Un_
Ax
for j = m + 1/2 and j = m — 1/2.
The method is stable if the CFL condition
max I Ai, Ao I —— = max Ax
A
±
]ja{a-l) (
+
Gn 2P
A_
1/2
Af Ax
< 1,
is satisfied.
2
5. Numerical experiments
Wo use the numerical method described above to study blood flow through the aneurysmal prostheses (see Fig. 3). We considered two kinds of prostheses:
• Fabric covered "rigid" stent-grafts, such as AneuRx ™.
•
The goal of this study was to understand the influence of the elastic properties of the prosthesis on the stresses exerted by the prostheses to the walls of the native aorta. In particular we were interested in understanding what hemodynamic factors are related to the occlusion of graft limbs, observed in patients and reported in [5. 16].
Many researchers have shown that the wall shear stress has a significant influence on blood coagulation and thrombosis, endothelial cell structure and function.
C.K. Zarins et al.( [17. 18]) proposed a method called Oscillatory Shear Index (OSI) to quantify the degree of oscillation in shear direction. OSI can be calculated using the following formula:
OSI = lAneg|
|Apos 1 + |Aneg 1
where Ane^d Apos are the areas under the shear stress vs. time when the shear stress is negative and positive respectively.
M. Haidokkor. C. White and J.A. Frangos suggested that the spatial and temporal gradients of the shear stress must be separated in order to identify which one is the primary cause of atherosclerotic plaque. In their paper [19]. they presented a study implicating high shear stress rate (or the temporal gradient of shear stress) as the main hemodynamic factor responsible for the initiation of focal atherogenesis.
We studied the Oscillatory Shear Index and the shear stress rate to detect possible causes for graft limb occlusion.
The following two parameters were varied in our study:
•
•
We have used the following data in our computations:
• Young's modulus Ec = 8700 Pa and Er = 217500 Pa for compliant and rigid stents respectively
• Young's modulus of a human aorta Ea = 105 Pa
Oscillatory Shear Index (OSI) Oscillatory Shear Index (OSI)
Fig. 4. Oscillatory Shear Index for rigid st.eut.-graft with limb diameter equal to 12 mm and main body diameter 23 mm (left.) and with limb diameter equal to 16 mm and main body diameter 28.5 mm (right)
Oscillatory Shear Index (OSI)
Fig. 5. Oscillatory Shear Index for compliant stent with limb diameter equal to 14.5 mm and main body diameter 23 mm
First, wo present the results of the numerically computed OSI for two different prostheses flexibility parameters and different diameters of''parent" (main body) and "daughter" (limb) components.
On each of the Fig. 4, 5 the solid curve shows the Oscillatory Shear Index for the main body component from the anchoring site to the bifurcation point and the ''solidstar" curve presents the OSI corresponding to the limb component from the bifurcation point downstream to the distal anchoring site.The Young's modulus of the prostheses shown in Fig. 4 corresponds to the one of rigid stent-graft, the main body diameters are 23 mm and 28.5 mm and the limb diameters are 12 mm and 16 mm. We observed that the smaller the diameter of the prosthesis the higher the OSI. We also tested a flexible prosthesis with a Young's modulus Ec = 8700 Pa (Fig. 5). The diameter of the main body for this case was taken to be 20 mm and the limb diameter is 14.5 mm. The
5 10 15 20 25 30 35 40 45 50 55 60 20 5 10 15 20 25 30 35 40 45 50 55 60
mesh points mesh points
Fig. 6. Shear Stress Rates for rigid stout-graft. with limb diameter equal to 12 mm and main body diameter 23 mm (left.) and with limb diameter equal to 16 mm and main body diameter 28.5 mm (right)
Time derivative of Shear Stress (main body diameter=28.5mm. limbs diameters=14.5mm) Time derivative of Shear Stress (main body diameter=23mm. limbs diameters=14.5mm)
5 10 15 20 25 30 35 40 45 50 55 60 5 10 15 20 25 30 35 40 45 50 55 60
mesh points mesh points
Fig. 7. Shear Stress Rates for the rigid st.eut.-graft. with a limb diameter equal to 14.5 mm and a main body diameter 28.5 mm (left) and wit.li a limb diameter equal to 14.5 mm and a main body diameter 23 mm (right)
amplitude of the Oscillatory Shear Index drastically increased for this case, indicating a much higher probability for occlusion.
Although the OSI showed that the smaller the diameter of the limbs the higher the OSI. it did not capture the patient-observed prosthesis characteristics that show higher probability for thrombosis at the distal site and not typically at the proximal site of the prosthesis. In contrast, the shear stress rate of the bifurcated prostheses did show to be higher at the distal anchoring site thereby conforming better to the complications observed in patients.
Fig. C 9 present the temporal gradient of shear stress at the systolic peak for different prostheses data. Again, each figure has two curves: solid curve corresponds to the shear stress rate of the main body component and "star" curve shows shear stress rate along the limb.
Fig. C shows the shear stress rates corresponding to the two extreme cases for compliant prosthesis: the prosthesis with the smallest limb diameter (12 mm) and the prosthesis with the largest limb diameter (16 mm), respectively. We observe that the smaller
mesh points
Fig. 8. Shear Stress Rates for the compliant, endoprosthesis with a limb diameter equal to 14.5 mm and a main body diameter of 23 mm
Time derivative of Shear Stress (main body diameter=28.5mm, limbs diameters=20mm)
Oscillatory Shear Index (OSI)
Dm6=28.5mm _ D,=20mm
45 50 55
80
70
50
50
Fig. 9. Shear Stress Rates (left.) and Oscillatory Shear Index (right.) for an "optimal" stent, with a limb diameter equal to 20 mm and a main body diameter of 28.5 mm
the diameter of the limbs the larger the shear stress rate, indicating a higher chance for occlusion. The worst performance was obtained for the smallest prostheses (Fig. 6 (left)) and the best performance was for prostheses with a limb diameter equal to 16 mm (Fig. 6 (right)).
We also observed that the size of the diameter of the main body influences the magnitude of the shear stress rate. Fig. 7 shows the comparison between two prostheses with the same limb diameters, but with different sizes for the main body: one with a main body diameter of 23 mm and the second one with a main body diameter of 28.5 mm. We can see an improvement in the shear stress rates for the larger main body prosthesis.
We calculated the shear stress rates for a bifurcated prosthesis with a compliancy Ec = 8700 Pa. As in the previous experiments for the "stiff" stent, we assumed that the diameter of the main modular component ranged from 23 mm to 28.5 mm, and that
the limbs diameter is between 12 mm and 16 mm. Fig. 8 presents the results of the numerical computations which indicate that the shear stress rates are drastically higher then for the case of stiff prosthesis.
The results presented above indicate that the magnitude of the shear stress rate is affected mainly by the size of the limbs and by the prostheses elasticity. In an attempt to obtain an optimal design minimizing shear stress rates we tested different prosthesis structures varying the limb diameter and elasticity parameters. We have found that the shear stress rate is minimized for the prosthesis with the following two characteristic:
• The diameter of the prosthesis limbs should be around \/2/2 of the diameter of the main body component: if Dmb denotes the diameter of the main body component and Di denotes the diameter of the limb component, then we obtained numerically that the lowest shear stress rates in the limbs are observed when D\ = Dmb■ This relationship is a consequence of the conservation of mass principle.
• Variable elasticity: The prosthesis which is stiffer in the central section, where there is no support to the prosthesis walls by the walls of native aorta, and softer at the overlap with the iliac arteries, showed best performance. In this test we nsed the fact that the stiffness of the composit prosthesis/vessel structure is eqnal to the combined stiffness of each structure, [20]. Hence, the smaller the stiffness of the prosthesis in the overlap anchoring region, the smaller the difference between the stiffness of the native vessel and prosthesis.
Fig. 9 (left) shows the behavior of the shear stress rates for the "optimal" prosthesis with characteristics obtained from the suggestions above:
• Young's modulus in the overlap region Eoveriap = 8700 Pa,
•
100000
• diameter of the main body component is Dmb = 28.5 mm,
• diameter of the limbs is Di = 20 mm (which is about Dmb )•
We observed a drastic improvement in the limb shear stress rates with respect to the previous prosthesis structures.
We also calculated the Oscillatory Shear Index for the "optimal prosthesis" described above (see Fig. 9 (right)). The results also showed a much smaller amplitude of OSI compared to the cases given in Fig. (4), (5).
Резюме
С. А. Лапии, С. Чаиич. Численное моделирование дизайна раздваивающихся протезов используемых в лечении аневризмы брюшной аорты.
В статье рассматривается численное решение задачи взаимодействия крови со стенками артерий. Используется одномерная эффективная модель полученная асимптотически из уравнений Навье Стокса для течения крови и уравнения лилейной эластичной мембраны описывающего поведение стенок артерий. В этой работе рассматривается случай раздваивающихся эластичных артерий со вставленным протезом. Построенная модель используется для определения влияния касательного напряжения па появление и развитие атерогепеза.
Literature
1. Canic S., Krajcer Z., Lapin S. Design of Optimal Endoprosthesis Using Mathematical Modeling // Endovascular Today. - 2006. - May. - P. 48-50 (Cover story).
2. Nichols W.W., O'Rourke M.F. McDonald's Blood Flow in Arteries: Theoretical, experimental and clinical principles. - N. Y., London, Sydney, Auckland: Arnold and Oxford University Press Inc., 1998.
3. Wang R., Ravi-Chandar K. Mechanical response of a metallic aortic stent-I. Pressure-diameter relationship // J. Appl. Mech. - 2004. - V. 71. - P. 697-705.
4. Wang R., Ravi-Chandar K. Mechanical response of a metallic aortic stent-II. A beam on elastic foundation model // J. Appl. Mech. - 2004 - V. 71. - P. 706-712.
5. Umscheid T., Stelter W.J. Time-related alterations in shape, position, and structure of self-expanding, modular aortic stent-grafts // J. Endovasc. Surg. - 1999. - V. 6 - P. 17-32.
6. Canic S., Hartley C.J., Rosenstrauch D., Tambaca J., Guidoboni G., Mikelic A. Blood flow in compliant arteries: an effective viscoelastic reduced model, numerics and experimental validation // Ann. Biomed. Eng. - 2006. - V. 34. - P. 575-592.
7. Canic S., Mikelic S.A., Lamponi D., Tambaca J. Self-consistent effective equations modeling blood flow in medium-to-large compliant arteries // SIAM J. Multiscale Analysis and Simulation. - 2005 - V. 3, No 3. - P. 559-596.
8. Canic S., Kim E.H. Mathematical analysis of the quasilinear effects in a hyperbolic model of blood flow through compliant axisymmetric vessels // Math. Methods in Appl. Sciences. - 2003. - V. 26. - P. 1161-1186.
9. White F.M. Viscous Fluid Flow. - N. Y.: McGraw-Hill, 1974.
10. Smith N.P., Pullan A.J., Hunter P.J. The generation of an anatomically accurate geometric coronary model // Ann. Biomed. Eng. - 2000. - V. 28, No 1. - P. 14-25.
11. Smith N.P., Pullan A.J., Hunter P.J. An anatomically based model of transient coronary blood flow in the heart // SIAM J. Appl. Math. - 2002. - V. 62. - P. 990-1018.
12. Olufsen M., Peskin C., Kim W., Pedersen E., Nadim A., Larsen J. Numerical simulation and experimental validation of blood flow in arteries with structured-tree outflow conditions // Ann. Biomed. Eng. - 2000. - V. 28. - P. 1281-1299.
13. Marusic-Paloka E. Fluid flow through a network of thin pipes // Comptes Rendus de l'Academie des Sciences Paris, Serie II, Fascicule B - Mecanique. - 2001. - V. 329, No 2. -P. 103-108.
14. Smoller J. Shock waves and reaction-diffusion equations. - N. Y.: Springer-Verlag, 1994.
15. Leveque R. Numerical Methods for Conservation Laws. - Basel; Boston: Birkhauser, 1992.
16. Parent F.N., Godziachvilli V., Meier G.H., Parker F.M., Carter K., Gayle R.G., De-massi R.J., Gregory R.T. Endograft limb occlusion and stenosis after ANCURE endovascular abdominal aneurysm repair // J. of Vasc. Surg. - 2002. - V. 35, No 4. -P. 686-690.
17. Ku D.N., Giddens D.P., Zarins C.K., Glagov S. Pulsatile flow and atherosclerosis in the human carotid bifurcation: positive correlation between plaque location and low and oscillating shear stress // Atherosclerosis. - 1985. - V. 5. - P. 293-302.
18. Moore J.E., Xu C., Glagov S., Zarins C.K., Ku D.N. Fluid wall shear stress measurements in a model of the human abdominal aorta: oscillatory behavior and relationship to atherosclerosis // Atherosclerosis. - 1994. - V. 110. - P. 225-240.
19. Haidekker M.A., White C.R., Frangos J.A. Analysis of temporal shear stress gradients during the onset phase of flow over a backward-facing step //J. Biomech. Eng. - 2001. -V. 123, No 5. - P. 455-463.
20. Dyet J.F., Watts W.G., Ettles D.E., Nicholson A.A. Mechanical properties of metallic stents: How do these properties influence the choice of stent for specific leisons? // Cardiovasc. Interv. Radiology. - 2000. - V. 23. - P. 47-54.
Поступила в редакцию 15.08.06
Лапин Сергей Александрович ассистент Департамента математики Университета Хьюстона, г. Хьюстон, США.
E-mail: slapinemath.uh.edu
Чанич Сунчица профессор Департамента математики Университета Хьюстона, г. Хьюстон, США.
E-mail: eanieemath.uh.edu
