基于格子Boltzmann方法和大涡模拟的颈动脉分叉狭窄流动并行计算

doi:10.11772/j.issn.1001-9081.2019081388

《计算机应用》唯一官方网站 ›› 2020, Vol. 40 ›› Issue (2): 404-409.DOI: 10.11772/j.issn.1001-9081.2019081388

• 2019年全国开放式分布与并行计算学术年会(DPCS 2019)论文 • 上一篇下一篇

基于格子Boltzmann方法和大涡模拟的颈动脉分叉狭窄流动并行计算

张毅卓¹, 葛森², 王良军², 谢江², 曹洁³, 张武¹^,²()

^1.上海大学上海市应用数学与力学研究所，上海 200444
^2.上海大学计算机工程与科学学院，上海 200444
^3.上海交通大学附属仁济医院老年科，上海 310101

收稿日期:2019-07-31 修回日期:2019-08-26 接受日期:2019-09-17 发布日期:2019-10-14 出版日期:2020-02-10
通讯作者: 张武
作者简介:张毅卓（1995—），男，湖北襄阳人，硕士研究生，CCF会员，主要研究方向：计算流体力学、高性能计算
葛森（1994—），男，河南开封人，硕士研究生，主要研究方向：高性能计算
王良军（1991—），男，安徽阜阳人，硕士研究生，主要研究方向：高性能计算、计算流体力学
谢江（1971—），女，湖北恩施人，副教授，博士，CCF高级会员，主要研究方向：生物信息学、高性能计算
曹洁（1980—），女，上海人，主治医师，硕士，主要研究方向：老年血脂异常；
基金资助:
国家自然科学基金面上项目(61873156);国家自然科学基金重大研究计划重点项目(91630206)

Parallel computing of bifurcation stenosis flows of carotid artery based on lattice Boltzmann method and large eddy simulation model

Yizhuo ZHANG¹, Sen GE², Liangjun WANG², Jiang XIE², Jie CAO³, Wu ZHANG¹^,²()

^1.Shanghai Institute of Applied Mathematics and Mechanics，Shanghai University，Shanghai 200444，China
^2.School of Computer Engineering and Science，Shanghai University，Shanghai 200444，China
^3.Department of Geriatrics of RenJi Hospital，Shanghai Jiao Tong University，Shanghai 310101，China

Received:2019-07-31 Revised:2019-08-26 Accepted:2019-09-17 Online:2019-10-14 Published:2020-02-10
Contact: Wu ZHANG
About author:ZHANG Yizhuo， born in 1995， M. S. candidate. His research interests include computational fluid mechanics， high performance computing.
GE Sen， born in 1994， M. S. candidate. His research interests include high performance computing.
WANG Liangjun， born in 1991， M. S. candidate. His research interests include high performance computing， computational fluid mechanics.
XIE Jiang， born in 1971， Ph. D.， associate professor. Her research interests include bioinformatics， high performance computing.
CAO Jie， born in 1980， M. S.， attending physician. Her research interests include dyslipidemia in elderly.
Supported by:
the Surface Program of National Natural Science Foundation of China(61873156);the Key Project of Major Research Program of National Natural Science Foundation of China(91630206)

摘要/Abstract

摘要：

颈动脉斑块的形成与复杂的血流动力学因素密切相关，血液流动状况的精确模拟对颈动脉斑块的临床诊断具有重要意义。为了精确模拟脉动流场，在格子Boltzmann方法（LBM）的基础上，添加大涡模拟（LES）模型，建立了LBM-LES颈动脉模拟算法。利用医学图像重构软件，建立颈动脉狭窄真实几何模型，对颈动脉狭窄脉动流动进行了数值模拟，通过计算血液流动速度、壁面剪切应力（WSS）等，得出了有意义的流动结果，验证了LBM-LES对颈动脉狭窄后段血液流动研究的有效性。基于OpenMP编程环境，在高性能集群机全互联胖节点上进行了千万量级网格的并行计算，结果表明LBM-LES颈动脉模拟算法具有较好的并行性能。

关键词: 格子Boltzmann方法, 大涡模拟, 颈动脉狭窄, 壁面剪切力, OpenMP

Abstract:

The formation of carotid artery plaque is closely related to complex hemodynamic factors. The accurate simulation of complex flow conditions is of great significance for the clinical diagnosis of carotid artery plaque. In order to simulate the pulsating flow accurately， Large Eddy Simulation （LES） model was combined with Lattice Boltzmann Method （LBM） to construct a LBM-LES carotid artery simulation algorithm， and a real geometric model of carotid artery stenosis was established through medical image reconstruction software， thus the high-resolution numerical simulation of carotid artery stenosis flows was conducted. By calculating blood flow velocity and Wall Shear Stress （WSS）， some meaningful flow results were obtained， proving the effectiveness of LBM-LES in the study of blood flow in the carotid artery narrow posterior. Based on the OpenMP programming environment， the parallel computation of the grid of ten million magnitude was carried out on the fully interconnected fat node of high-performance cluster machine. The results show that the LBM-LES carotid artery simulation algorithm has good parallel performance.

Key words: Lattice Boltzmann Method (LBM), Large Eddy Simulation (LES), carotid artery stenosis, Wall Share Stress (WSS), OpenMP (Open Multi-Processing)

中图分类号:

TP301.6

张毅卓, 葛森, 王良军, 谢江, 曹洁, 张武. 基于格子Boltzmann方法和大涡模拟的颈动脉分叉狭窄流动并行计算[J]. 计算机应用, 2020, 40(2): 404-409.

Yizhuo ZHANG, Sen GE, Liangjun WANG, Jiang XIE, Jie CAO, Wu ZHANG. Parallel computing of bifurcation stenosis flows of carotid artery based on lattice Boltzmann method and large eddy simulation model[J]. Journal of Computer Applications, 2020, 40(2): 404-409.

图/表 8

图1 D3Q19离散速度模型

Fig. 1 D3Q19 discrete velocity model

图2 本文并行算法流程

Fig. 2 Flowchart of the proposed parallel algorithm

图3 三维血管模型

Fig. 3 Three-dimensional carotid model

表1 脉动血流参数

Tab. 1 Pulsating blood flow parameter

时间点	流速/（m·s^-1）	时间点	流速/（m·s^-1）
VM	0.18	VDN	0.15
VMAX	0.80	VDFH	0.30

图4 流速云图

Fig. 4 Velocity cloud chart

图5 剪切力矢量图

Fig. 5 WSS vector diagram

图6 LBM和LBM-LES模型在不同核数下的运行时间对比

Fig. 6 Running time comparison of LBM and LBM-LES models with different cores

图7 LBM和LBM-LES模型在不同核数下的加速比对比

Fig. 7 Speedup comparison of LBM and LBM-LES models with different cores

参考文献 23

1	金龙，乔爱科.颈动脉易损斑块的生物力学机制和破裂风险评价指标［J］.医用生物力学，2016，31（1）：89-94. 10.3871/j.1004-7220.2016.01.089
	JIN L， QIAO A K. Biomechanical mechanism and quantitative assessment indices for vulnerable carotid plaques［J］. Journal of Medical Biomechanics， 2016， 31（1）： 89-94. 10.3871/j.1004-7220.2016.01.089
2	刘国荣，王大力，张文丽，等.颈动脉易损斑块与缺血性脑卒中复发的相关性研究［J］.中华老年心脑血管病杂志，2012，14（10）：1067-1070. 10.3969/j.issn.1009-0126.2012.10.019
	LIU G R， WANG D L， ZHANG W L， et al. Correlation between carotid artery vulnerable plaques and recurrent is chemic cerebral stroke［J］. Chinese Journal of Geriatric Heart Brain and Vessel Diseases， 2012， 14（10）：1067-1070. 10.3969/j.issn.1009-0126.2012.10.019
3	杜宜纲，刘德杰，沈莹莹，等. 血管壁面剪切应力的测量及其临床研究进展［J］. 中国生物医学工程学报， 2018， 37（5）： 593-605. 10.3969/j.issn.0258-8021.2018.05.011
	DU Y G， LIU D J， SHEN Y Y， et al. Measurements of wall shear stress for blood vessel and its progress in clinical research［J］. Chinese Journal of Biomedical Engineering， 2018， 37（5）： 593-605. 10.3969/j.issn.0258-8021.2018.05.011
4	MALEK A M， ALPER S L， LZUMO S. Hemodynamic shear stress and its role in atherosclerosis［J］. Journal of American Medical Association， 1999， 282（21）： 2035-2042. 10.1001/jama.282.21.2035
5	GOHIL T， MCGREGOR R H P， SZCZERBA D， et al. Simulation of oscillatory flow in an aortic bifurcation using FVM and FEM： a comparative study of implementation strategies［J］. International Journal for Numerical Methods in Fluids， 2011， 66（8）： 1037-1067. 10.1002/fld.2301
6	何雅玲，王勇，李庆. 格子Boltzmann方法的理论及应用［M］. 北京：科学出版社， 2009：1-4. 10.1007/s11434-009-0681-6
	HE Y L， WANG Y， LI Q. Lattice Boltzmann Method： Theory and Applications［M］. Beijing： Science Press， 2009： 1-4. 10.1007/s11434-009-0681-6
7	AIDUN C K， CLAUSEN J R. Lattice-Boltzmann method for complex flows［J］. Annual Review of Fluid Mechanics， 2010， 42： 439-472. 10.1146/annurev-fluid-121108-145519
8	KANG X， TANG W， LIU S. Lattice Boltzmann method for simulating disturbed hemodynamic characteristics of blood flow in stenosed human carotid bifurcation［J］. Journal of Fluids Engineering， 2016， 138（12）： No.121104. 10.1115/1.4033913
9	KARIMPOUR H， JAVDAN E. Simulation of stenosis growth in the carotid artery by lattice-Boltzmann method［J］. Journal of Mechanics in Medicine and Biology， 2014， 14（2）： No.1450016. 10.1142/s021951941450016x
10	MATYKA M， KOZA Z， MIROSLAW Ł. Wall orientation and shear stress in the lattice Boltzmann model［J］. Computers and Fluids， 2013， 73： 115-123. 10.1016/j.compfluid.2012.12.018
11	杨帆，刘树红，唐学林，等.格子Boltzmann亚格子模型的研究［J］.工程热物理学报，2004，25（S1）：43-46. 10.1115/ht-fed2004-56100
	YANG F， LIU S H， TANG X L， et al. Study on subgrid turbulence model for lattice Boltzmann method［J］. Journal of Engineering Thermophysics， 2004， 25（S1）：43-46. 10.1115/ht-fed2004-56100
12	王玲玲. 大涡模拟理论及其应用综述［J］. 河海大学学报（自然科学版）， 2004， 32（3）：261-265. 10.3321/j.issn:1000-1980.2004.03.005
	WANG L L. Large eddy simulation theory and its application： a review［J］. Journal of Hohai University （Natural Sciences）， 2004， 32（3）： 261-265. 10.3321/j.issn:1000-1980.2004.03.005
13	PAL A， ANUPINDI K， DELORME Y， et al. Large eddy simulation of transitional flow in an idealized stenotic blood vessel： evaluation of subgrid scale models［J］. Journal of Biomechanical Engineering， 2014， 136（7）： No.071009. 10.1115/1.4027610
14	LU H， RUTLAND C J. Structural subgrid-scale modeling for large-eddy simulation： a review［J］. Acta Mechanica Sinica， 2016， 32（4）： 567-578. 10.1007/s10409-016-0556-4
15	LANTZ J， GÅRDHAGEN R， KARLSSON M. Quantifying turbulent wall shear stress in a subject specific human aorta using large eddy simulation［J］. Medical Engineering and Physics， 2012， 34（8）： 1139-1148. 10.1016/j.medengphy.2011.12.002
16	VERGARA C， LE VAN D， QUADRIO M， et al. Large eddy simulations of blood dynamics in abdominal aortic aneurysms［J］. Medical Engineering and Physics， 2017，47： 38-46. 10.1016/j.medengphy.2017.06.030
17	刘智翔，宋安平，徐磊，等. 多重网格格子Boltzmann方法的并行算法［J］. 计算机应用， 2014， 34（11）：3065-3068. 10.11772/j.issn.1001-9081.2014.11.3065
	LIU Z X， SONG A P， XU L， et al. Parallel algorithms for multi-grid lattice Boltzmann method［J］. Journal of Computer Applications， 2014， 34（11）： 3065-3068. 10.11772/j.issn.1001-9081.2014.11.3065
18	SABRI K， RABBANI H， GUNAWAN P H. OpenMP performance for benchmark 2D shallow water equations using LBM［J］. Journal of Physics： Conference Series， 2018， 971： No.012033. 10.1088/1742-6596/971/1/012033
19	李清宝，张平. 基于分布/共享内存层次结构的并行程序设计［J］. 计算机应用， 2004， 24（6）：148-150， 158.
	LI Q B， ZHANG P. Parallel programming based on distributed/shared memory hierarchical architecture［J］. Journal of Computer Applications， 2004， 24（6）： 148-150， 158.
20	QIAN Y H， D’HUMIÈRES， LALLEMAND P. Lattice BGK models for Navier-Stokes equation［J］. Europhysics Letters， 1992， 17（6）： 479-484. 10.1209/0295-5075/17/6/001
21	YU H， LUO L， GIRIMAJI S S. LES of turbulent square jet flow using an MRT lattice Boltzmann model［J］. Computers and Fluids， 2006， 35（8/9）： 957-965. 10.1016/j.compfluid.2005.04.009
22	王娇，刘洋，张晓玲. Mimics软件在医学图像三维重建中的应用［J］. 医疗卫生装备， 2015， 36（2）：115-118. 10.7687/J.ISSN1003-8868.2015.02.115
	WANG J， LIU Y， ZHANG X L. Application of Mimics software to 3D reconstruction of medical image［J］. Chinese Medical Equipment Journal， 2015， 36（2）： 115-118. 10.7687/J.ISSN1003-8868.2015.02.115
23	YU D， REN W， SHYY W. A unified boundary treatment in lattice Boltzmann method ［C/OL］// Proceedings of the 41st Aerospace Sciences Meeting & Exhibit. Reno， Nevada： AIAA， 2003 ［2019-02-24］. . 10.2514/6.2003-953

[1]	刘智翔, 刘慧超, 黄冬梅, 周丽萍, 苏诚. 多种任务调度混合的IB-LBM并行优化方法[J]. 《计算机应用》唯一官方网站, 2020, 40(2): 386-391.
[2]	胡园园, 谢江, 张武. 二维不可压缩Navier-Stokes方程的并行谱有限元法求解[J]. 计算机应用, 2017, 37(1): 42-47.
[3]	刘有耀, 杨鹏程. 基于JavaCC的C代码自动并行化的设计与实现[J]. 计算机应用, 2016, 36(9): 2422-2426.
[4]	吴洁璇, 陈振杰, 张云倩, 骈宇哲, 周琛. 多核CPU下的K-means遥感影像分类并行方法[J]. 计算机应用, 2015, 35(5): 1296-1301.
[5]	刘智翔宋安平徐磊郑汉垣张武. 多重网格格子Boltzmann方法的并行算法[J]. 计算机应用, 2014, 34(11): 3065-3068.
[6]	李双李文敬孙环龙林中明. 基于多核机群的人工鱼群并行算法[J]. 计算机应用, 2013, 33(12): 3380-3384.
[7]	兰中周乐励华高云. 基于格子Boltzmann方法的一维Burgers方程的数值模拟[J]. 计算机应用, 2013, 33(09): 2432-2435.
[8]	巫小婷邓家先. 基于OpenMP的压缩感知并行处理算法[J]. 计算机应用, 2012, 32(03): 617-619.
[9]	白明泽程丽豆育升孙世新. 基于OpenMP的分子动力学并行算法的性能分析与优化[J]. 计算机应用, 2012, 32(01): 163-166.
[10]	林江宏林锦贤吕暾. 多核CPU和GPU加速分子动力学模拟[J]. 计算机应用, 2011, 31(03): 843-847.

基于格子Boltzmann方法和大涡模拟的颈动脉分叉狭窄流动并行计算

Parallel computing of bifurcation stenosis flows of carotid artery based on lattice Boltzmann method and large eddy simulation model

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表 8

参考文献 23

相关文章 10

编辑推荐

Metrics