均方误差
建模和模拟在工程
1687 - 5605
1687 - 5591
Hindawi出版公司
691982年
10.1155 / 2008/691982
691982年
研究文章
血流动力学效应的不平等的大脑前动脉流速前沟通动脉分叉:计算流体动力学的研究
劳
托马斯。
1
他
兴
2
Venugopal
但是过犹不及
2
Vinuela
费尔南多
2
Duckwiler
加里
2
华伦天奴
丹尼尔·J。
2
Pietka
Ewa
1
美国放射学
科罗拉多大学健康科学中心
极光
有限公司80045
美国
rosalindfranklin.edu
2
放射科学系的
加州大学
洛杉矶CA 90095
美国
ucla.edu
2008年
11
01
2009年
2008年
12
11
2007年
11
06
2008年
11
12
2008年
2008年
版权©2008
这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。
计算流体动力学技术被用来调查不平等的大脑前动脉的血流动力学影响流率前大脑前动脉(ACA-ACOM)和分岔。血液动力学一直与脑血管疾病的一个主要因素。使用一个理想化的二维对称ACA-ACOM几何模型的流场和壁面切应力(WSS)分叉区域评估为脉动的流入左向右流动比率1:1,2:1,3:1,4:1。不平等的流速通过ACA父动脉导致分歧的高父流流动和转移的撞击点沿着A2-ACOM相邻墙侧ACA。错流穿过ACOM通常不稳定,导致WSS撞击地区从更高的流量增加父母动脉和双振幅峰值WSS的侧分岔地区从当地再循环的影响。这些结果表明,不对称ACA流速增加血流动力学压力在ACA-ACOM分叉区域和显示船的可能因素削弱和动脉瘤形成。
1。介绍
脑底动脉环的血液动力学多项调查已经进行阐明动脉瘤的发病机制的形成。众所周知,动脉瘤的85%以上发生在前循环和周围形成主要在船顶分岔和锋利的弯曲。在这些地方血流动力学压力是最大和结构蛋白装配的肌肉层的血管壁最不发达
1 ]。普遍接受的理论是一个异常反应血管内皮的振荡壁剪切应力导致进步和不稳定的动脉壁退化和顺向动脉瘤的形成
2 ]。计算方法的有效性和可靠性模型的血流动力学过程多次动力计算流研究描述相关的流模式和壁剪切应力脉动的流动在这些敏感地区(
3 - - - - - -
9 ]。
前的结脑动脉(ACA)和前沟通(ACOM)代表一个独特而复杂的血管几何两个紧密不是分岔(见图
1 )。过去的研究表明,一个强大的协会的发展之间存在动脉瘤和颈部尺寸ACOM结与占主导地位的A1父动脉支线(
10 ]。在这项研究中,计算流体力学(CFD)模拟了一个理想化的二维对称模型的ACA和ACOM分叉调查的血流动力学效应通过A1 ACA段流速不均。模拟运行与左到右流入比例固定在1:1,2:1,3:1,4:1。整体流动模式以及壁面切应力(WSS)沿着相邻壁各点之间的A2和ACOM段时空上分析了撞击点附近的动脉瘤预计发生。不平等的左右A1流速时可能发生总双边上游流电阻是不平等的。临床上,这可能出现单边优质狭窄的病变上游血管动脉环。这可能是在中风患者血栓栓塞事件后,成功在上游动脉血管痉挛或存在先天性发育不全的A1-ACA。A1段之间的不对称流量也会发生暂时性的医源性意味着在快速静脉注入造影剂在诊断脑血管造影术。本研究试图描述一些血流动力学压力模式差异可能发生在这种情况下。
图1
前的3 d模型描述脑循环的模拟。
2。方法
2 d ACA和ACOM对称模型构造使用Solidworks [
11 如图)软件
2 。直径被分配一个单位1的长度,和所有其他尺寸缩放比例的相应和从其他研究
12 ,
13 ]。分岔角的A1, A2段前部脑动脉被选为45度,在临床模型的范围。曲率半径的A1, A2的墙壁和ACOM段选择这样平稳过渡。然后模型导入到流利的啮合的策略。finer-graded网是用于边界层沿墙表面增加梯度的分辨率在墙附近。总共196 944数据点啮合后得到。
与45度(a) ACA-ACOM几何分叉角度。(b)给出的脉冲波形(
3 )。全球最大的发生在0.45
t
/
T
将定义的收缩压和舒张压之间的过渡阶段流波形。
(一)
(b)
替换下面的无量纲速度表达式
u
和
v
、长度
x
和
y
,时间
t
、压力
P
,雷诺数
再保险
让我们undimensionalize执政的连续性和n - s方程在笛卡儿坐标系统。我们把流变量与相应的引用变量用下标
裁判
生产用无量纲变量
*
。让
ρ
血液的密度,
μ
动态粘度的血液,
u
裁判
参考速度,
l
裁判
参考长度,
(1)
u
*
=
u
u
裁判
,
v
*
=
v
u
裁判
,
x
*
=
x
l
裁判
,
y
*
=
y
l
裁判
,
t
*
=
t
u
裁判
l
裁判
,
P
*
=
P
ρ
u
裁判
2
,
再保险
=
ρ
u
裁判
l
裁判
μ
。
连续性和n - s方程成为[
14 ]
(2)
∂
u
*
∂
x
*
+
∂
v
*
∂
y
*
=
0,
∂
u
*
∂
t
*
+
u
*
∂
u
*
∂
x
*
+
v
*
∂
u
*
∂
y
*
=
−
∂
p
*
∂
x
*
+
1
再保险
(
∂
2
u
*
∂
x
*
2
+
∂
2
u
*
∂
y
*
2
)
,
∂
v
*
∂
t
*
+
u
*
∂
v
*
∂
x
*
+
v
*
∂
v
*
∂
y
*
=
−
∂
p
*
∂
y
*
+
1
再保险
(
∂
2
v
*
∂
x
*
2
+
∂
2
v
*
∂
y
*
2
)
。
管理流动方程的使用这种形式消除依赖绝对尺寸,使仿真只依赖雷诺氏数。因此,获得的结果可以更广义。不可压缩牛顿流体,在无量纲n - s方程的计算假定,和执行的解决方案是使用OpenFOAM软件
15 ]。最低雷诺数400年设定的允许流体密度等于1,比例最大的速度剖面,和流体粘度设置为0.0025。雷诺氏数允许改变800的跳动的速度剖面的性质2水湾。Two-plug流动入口被分配在父母前脑动脉动脉波形脉冲定义为(
3 )。无滑动和非弹性条件强加于墙壁和流出边界被设定为一个参考压力为0。模拟被允许运行三到四个完整周期每分钟60次实现融合的解决方案。
模拟运行的比左派和右派之间的流入率A1段设置为1:1,2:1,3:1,4:1。这是由减少的时间振幅A1流入段由恰当的分数(0.5、0.33和0.25)。
速度,
U
,这些模拟是由脉冲配置文件中使用
(3)
U
(
t
)
=
8.419
*
{
0.08225
+
2
(
−
0.013077
因为
(
w
t
)
−
0.002868
罪
(
w
t
)
)
+
2
(
0.001125
因为
(
2
w
t
)
−
0.001625
罪
(
2
w
t
)
)
+
2
(
−
0.000172
因为
(
3
w
t
)
+
0.001631
罪
(
3
w
t
)
)
+
2
(
−
0.001000
因为
(
4
w
t
)
)
}
,
在哪里
w
=
2
π
/
T
,
T
=
1
第二。
收缩压和舒张压阶段之间的过渡点的速度波形将定义在全球最大出现在0.45
t
/
T
。
壁剪切应力(WSS)得到以下方程:
(4)
(
W
x
W
y
)
=
μ
(
2
∂
u
∂
x
∂
u
∂
y
+
∂
v
∂
x
∂
u
∂
y
+
∂
v
∂
x
2
∂
v
∂
y
)
(
n
x
n
y
)
,
在哪里
n
x
和
n
y
是
x
和
y
组件的血管壁的法向量。
3所示。结果
相对流入利率通过两个A1父动脉段大大影响整个流型和墙上的血流动力学压力分岔地区。当左右A1流入利率是不平等的,速度越高父流发生分岔和ACOM成为功能吻合。同时流过A1和A2 ACA段通常是层流,不规则中看到ACA-ACOM分叉部分尤其是在减速阶段的脉动。数据
4 ,
5 ,
6 显示较高的流的分岔父动脉流。流经结似乎更稳定在加速度脉冲周期的阶段。
平等左右A1流速,分离和涡脱落发生在ACOM结流绕过吻合到A2段。因此,虽然呈现ACOM死区,创建一个局部不稳定地区附近的撞击导致不规则,低振幅A2-ACOM WSS模式的墙。图
3 显示父的双边斜撞击区A1段到远端A2-ACOM毗邻的墙。WSS情节的人物
8 和
12 证明小振幅WSS和梯度相比不平等ACA流量情况。
速度场与平等的左派和右派ACA流率(
问
R
=
问
l
)在(a)早期舒张收缩和(b)。ACOM存在的死区在整个脉动周期。(b)不稳定存在双边在A2-ACOM撞击地点与流动分离和旋涡在舒张期阶段。
(一)
(b)
(一)速度场在收缩和(b)与ACA流速设定在舒张早期
问
R
=
0.5
问
l
。离开ACA撞击地点沿A2-ACOM毗邻的向右移动墙的情况相比
问
R
=
问
l
分岔的流入A2-ACA和ACOM段。(b)生成一个大漩涡右边分岔网站在舒张。
(一)
(b)
(一)速度场在收缩和(b)与ACA流速设定在舒张早期
问
R
=
0.33
问
l
。再分岔的A1父动脉流入ACOM和A2段。离开ACA撞击网站进一步向右移动沿着A2-ACOM相邻墙相比,情况
问
R
=
0.5
问
l
。流动分离再次看到右边A2-ACOM撞击网站导致当地再循环的影响。
(一)
(b)
(一)速度场在收缩和(b)与ACA流速设定在舒张早期
问
R
=
0.25
问
l
。分岔的A1父动脉流入ACOM和A2段再次看到。离开ACA撞击网站似乎将进一步向右沿着A2-ACOM相邻墙相比,情况
问
R
=
0.33
问
l
。流动分离再次看到右边A2-ACOM撞击网站导致当地再循环的影响。
(一)
(b)
从正确的ACA是减少流入率相对向左,撞击点的左边沿着A2-ACOM A1流走向侧分岔的墙段相邻,是分裂的。撞击点也观察到的振动向侧ACA段在减速阶段的脉动周期。大部分流入A2段而继续分裂的分数进入ACOM段和加入流过ACA的权利。流经ACOM吻合是主体分离,协商的曲率ACOM在减速阶段的脉动周期。流动分离的结果在一个大第二高峰时间WSS情节由于地方再循环的影响。再循环区是最明显的例子
问
R
=
0.5
问
l
在图
4 (b),但在其他情况下高倍镜下可以看到视图(没有显示)。第二个峰值的大小似乎增加的比率从左到右流率增加,发生在脉动周期的早期。因此,不平等的A1流速ACA-ACOM分岔地区施加额外压力的不平衡度左右A1流入率。
WSS的情节与时间和位置数据所示
8 ,
9 ,
10 ,
11 为各种左向右流动比率。沿着左和右位置撞击点是用L1-L8 R1-R8,分别如图所示
7 。从左到右流速不均的情节人物
9 - - - - - -
11 显示双振幅峰值发生在舒张期阶段时不在左右流速相等(见图
8 )。图
12 情节WSS和时间点左右撞击区域特征。显然在这个图是双振幅峰值发生在不平等的父母A1流和之间的时间连续山峰之间的失衡A1流动更大。
图7
标记点的左派和右派撞击点WSS的情节。
WSS和位置和时间
问
R
=
问
l
在左边(a)和(b)的分岔地区。
(一)
(b)
WSS和位置和时间
问
R
=
0.5
问
l
在左边(a)和(b)的分岔地区。
(一)
(b)
WSS和位置和时间
问
R
=
0.33
问
l
在左边(a)和(b)的分岔地区。
(一)
(b)
WSS和位置和时间
问
R
=
0.25
问
l
在左边(a)和(b)的分岔地区。
(一)
(b)
WSS左边(a)和(b)的撞击点。
(一)
(b)
4所示。讨论
结果表明,流动行为和血流动力学强调ACOM分叉高度依赖之间的相对流量比率A1-ACA父动脉。不对称流条件在A1段对称几何形状导致分歧的高父动脉流入A2和ACOM段。通过ACOM从高流错流父动脉通常是不稳定的分岔地区特别是在舒张。使用创建一个周期性不同地区的流动分离和逆转的位置在整个周期振幅变化明显。分叉区域,这些不稳定的发生,是网站的大型WSS的波动大小和预计与动脉瘤的位置相一致。
不平等之间流动左右WSS A1段产生一种变化模式相比,对称的流动条件。同时WSS振幅的增加撞击点和双振幅峰值在对面的分岔地区观察到由于地方再循环的影响。传入的父动脉流的撞击地点也转向的顶点在A2-ACOM分岔地区毗邻的墙,中间的动脉壁肌层是最不发达
1 ]。这将导致更大的血流动力学压力位置,至少能够承受机械变形和可能的机制从循环疲劳应力。数据
3 - - - - - -
6 也证明撞击地点从占主导地位的A1转变向侧ACA与A1流速不均,流到一个更小的撞击区更锐角撞击表面的正常。此外,冲击地区倾向于沿着A2-ACOM振荡相邻墙与侧转向ACA在收缩和舒张期间亦然。相比之下,流速相等的情况说明(见图
3 )传入流冲击更斜角对撞击区域,减少振动。Cebral等人发现动脉瘤较小和撞击区域和扰流模式的变化与断裂有关。这些相关的发现进一步加强动脉瘤增长的不稳定发展协会拥有大量时间梯度WSS。
数据
4 - - - - - -
6 证明大空间WSS梯度在左边撞击的地方越高流父流生长的影响。空间梯度从区域流扰动干扰层流剪切应力是重要的地方发现了调节器的内皮细胞表达各种生长因子(
16 ]。正确的分岔地区看到的情节人物
4 - - - - - -
6 ,
12 证明流动分离地区的不稳定影响流速不均受到双重高峰WSS的演化。虽然WSS明显少的不平等的A1级流量条件相比,相同流量条件下,它已经表明,低振幅振动不稳定流动的剪切应力内皮细胞营业额的重要决定因素,可能比WSS振幅(更重要
17 ,
18 ]。另一个有趣的观察从图
12 是两座山峰之间的频率增加之间的不平衡的A1流入更大。增加的频率大约是总数的4%之间的脉动周期2:1和3:1例,增加到8%之间2:1和4:1研究。这些差异反映了时间在流动分离的脉动周期发生错流通过导航ACOM各种失衡A1段流。因此,不平等A1-ACA流速引起的血流动力学模式的改变压力可能影响内皮细胞的转导除了机械变形应力。
在某些情况下,流动分离已被证明是足以引起血小板活化,导致自由浮动总量(
19 - - - - - -
21 ]。这表明,某些地区的非定常流可能导致血栓栓塞现象。因此,ACA-ACOM结可作为血小板聚集的关键中心可能释放成分,成为栓子的来源。自发的血栓形成和血栓扩展双边留住的未破裂动脉瘤ACOM报告文学和易感几何证明血流动力学现象会导致脑血管意外(
22 ]。
有趣的是,在血流动力学改变压力发生变化,尽管总体减少流量通过前脑循环。不平等的流过留住是情况可能发生的病理临床情况以及医源性干预措施。颈动脉支架植入后颅内出血是一种已知的致命并发症在治疗后时期(
23 ]。事实上,血管内的血流动力学影响的干预措施,例如狭窄的血管成形术和支架植入值得进一步调查。给出的结果可能会提供一些见解的病因血管壁削弱甚至协助ACA-ACOM动脉瘤患者的临床管理。
目前的研究主要集中在不平等的影响家长ACA流速ACOM的血液动力学。通过使用一个2 d对称模型,可以控制和其他几何因素能获得快速的仿真结果。虽然二维模拟可以预期相当类似于三维流动的主要特征,2 d和3 d模拟之间的显著差异也预期出现的复杂的几何模型。本研究中使用的模型采用高度对称,试图减少这种3 d的混杂因素。更现实的特定病人的3 d模型是未来研究的主题。
5。结论
不平等的悸动的流速通过A1-ACA动脉导致分歧和不稳定通过ACOM从高流错流父动脉。有增加WSS撞击现场级相同的一侧作为父A1段和改变位置的最大相关WSS左向右流动比率。增加,预计动脉瘤血流动力学压力似乎发生发展。在侧A2-ACOM相邻墙,WSS的双峰值振幅从当地发生流动分离带来的不稳定影响横流协商ACOM段。这些增加血流动力学压力分岔地区可能是一个因素,血管壁削弱动脉瘤常常发生。
[
]1
李
r·m·k·W。
脑动脉的形态
药理学和治疗
1995年
66年
1
149年
173年
10.1016 / 0163 - 7258 (94)00071 - a
[
]2
戴维斯
p F。
杜威
c F。
Jr。
Bussolari
s R。
戈登
e . J。
Gimbrone
m·A。
Jr。
血流动力学部队对血管内皮功能的影响。体外研究的剪切应力和胞饮牛主动脉细胞
《临床研究杂志》上
1984年
73年
4
1121年
1129年
10.1172 / JCI111298
[
]3
福特
m D。
Stuhne
g·R。
Nikolov
h . N。
虚拟动脉动脉瘤血流动力学可视化和验证计算模型的
IEEE医学成像
2005年
24
12
1586年
1592年
10.1109 / TMI.2005.859204
[
]4
他
X。
Venugopal
P。
Cebral
j . R。
施密特
H。
华伦天奴
d . J。
DValentino@mednet.ucla.edu
再现性的脑血流动力学模拟:一个inter-solver比较
6143年
2006年医学成像:生理学、功能,从医学图像和结构
2006年2月
圣地亚哥,加州,美国
1
12
学报学报
10.1117/12.653970
[
]5
Cebral
j . R。
卡斯特罗
m·A。
伯吉斯
j·E。
Pergolizzi
r S。
谢里丹
m·J。
核
c . M。
描述脑动脉瘤破裂的风险评估通过使用特定的计算血流动力学模型
美国神经放射学杂志》
2005年
26
10
2550年
2559年
[
]6
Ku
d . N。
动脉血液流动
流体力学的年度审查
1997年
29日
399年
434年
10.1146 / annurev.fluid.29.1.399
[
]7
Ku
d . N。
吉登斯
d . P。
Zarins
c K。
Glagov
年代。
跳动的流和动脉粥样硬化在人类颈动脉分叉。斑块之间的正相关位置和低和振荡剪切应力
动脉硬化
1985年
5
3
293年
302年
[
]8
伯杰
美国一个。
saberger@me.berkeley.edu
周素卿
L.-D。
流在狭窄的血管
流体力学的年度审查
2000年
32
347年
382年
10.1146 / annurev.fluid.32.1.347
[
]9
泰勒
c。
休斯
t·j·R。
Zarins
c K。
有限元建模动脉的血流量
计算机在应用力学和工程方法
1998年
158年
1 - 2
155年
196年
10.1016 / s0045 - 7825 (98) 80008 - x
[
]10
夸克
R。
Ohi
T。
Niizuma
H。
铃木
J。
传入的脖子的动脉和前动脉瘤进行沟通
外科神经学
1980年
13
3
221年
223年
[
]11
http://www.solidworks.com
[
]12
Cassot
F。
Vergeur
V。
博须埃
P。
Hillen
B。
Zagzoule
M。
Marc-Vergnes
j。
前沟通动脉直径的影响在颈内动脉疾病:脑血液动力学模型研究
循环
1995年
92年
10
3122年
3131年
[
]13
蒂芬妮
m·A。
mstefani@vortex.ufrgs.br
施耐德
f . L。
马伦
a·c·H。
塞韦里诺
a·G。
Jackowski
答:P。
华莱士
m . C。
大脑前动脉皮质支的解剖变异
临床解剖学
2000年
13
4
231年
236年
10.1002 / 1098 - 2353(2000)十三4 < 231::AID-CA1 > 3.0.CO; 2 t
[
]14
后面;
G。
介绍了流体动力学
2000年
英国剑桥
剑桥大学出版社
[
]15
http://www.opencfd.co.uk/openfoam
[
]16
内格尔
T。
雷斯尼克
N。
杜威
c F。
Jr。
Gimbrone
m·A。
Jr。
gimbrone@bustoff.bwh.harvard.edu
血管内皮细胞对空间在流体剪切应力梯度增强激活的转录因子
动脉硬化、血栓和血管生物学
1999年
19
8
1825年
1834年
[
]17
戴维斯
p F。
Remuzzi
一个。
戈登
e . J。
杜威
c F。
Jr。
Gimbrone
m·A。
Jr。
湍流剪应力诱导血管内皮细胞体外营业额
美国国家科学院院刊》上的美利坚合众国
1986年
83年
7
2114年
2117年
10.1073 / pnas.83.7.2114
[
]18
商店
E。
jsho@stanford.edu
Nanjo
H。
商店
M。
动脉扩张、弯曲度和内膜的增厚,以应对连续暴露在高和低壁面切应力
血管外科杂志》
2004年
39
3
601年
612年
10.1016 / j.jvs.2003.10.058
[
]19
莫顿
w·A。
有土豆的
e . M。
Petschek
h·E。
研究聚合形成的区域分离血液流动
血栓形成等素质Haemorrhagica
1975年
34
3
840年
854年
[
]20.
Einav
年代。
einav@eng.tau.ac.il
Bluestein
D。
血流动力学和血小板在病态的船只运输
纽约科学院上
2004年
1015年
351年
366年
[
]21
el-Masry
o . A。
福伊尔斯坦
我一个。
轮
g F。
实验评估简化模式和分离流在一系列的分支血管对动脉粥样硬化和血栓形成的影响
循环研究
1978年
43
4
608年
618年
[
]22
Brownlee
r D。
Tranmer
我。
舍夫契克
r . J。
Karmy
G。
咖喱
b . J。
自发颅内前沟通动脉动脉瘤的栓塞:一个不寻常的缺血性中风的原因
中风
1995年
26
10
1945年
1949年
[
]23
莫里斯近日
w·F。
Grahovac
年代。
Douen
一个。
颅内出血后颅外颈动脉狭窄的支架和血管成形术
美国神经放射学杂志》
2000年
21
10
1911年
1916年