【技术实现步骤摘要】
本专利技术属于计算流体力学领域，具体涉及一种基于非平衡壁面模型的心脏瓣膜血液流动les/ib计算方法。
技术介绍
1、心血管疾病是当今几种主要的致死性疾病之一。由于疾病或先天问题，自然的心脏瓣膜会变得有缺陷，通过外科手术用人工瓣膜代替人体瓣膜已有五十多年的历史。其中，双叶机械心脏瓣膜(bmhv：bileaflet mechanicalheart valve)是最典型的机械瓣膜，如图1所示。然而，由于会产生非生理性的流动结构，mhv容易导致血栓栓塞等并发症。减弱非生理性流动，降低血栓栓塞风险，是mhv流动研究的重点。bhv的研究注重于瓣膜的有限元动力学分析以及血液与瓣膜之间的相互耦合作用。计算流体力学在心脏瓣膜流动研究中的显著优势是能够高分辨率地模拟三维流场，是推进相关血液动力学特性理解的有效选择。
2、mhv的存在使得流场结构变得更加复杂，一个脉动周期内，会出现流动向紊流转捩以及后来的再层流化等复杂现象。ib(immersed boundary)方法使得运动界面流动的数值模拟可以在固定网格上进行，避免了网格的重新生成以及不同网格之间的数值传递，是目前心脏血液流动模拟中的主流方法。les(large-eddy simulation)模拟中，网格计算节点正比于re13/7，壁面模化技术能够降低心脏瓣膜血液流动les模拟时近壁面附近的网格尺度要求，使得基于ib的湍流模拟数值方法可行。壁面模化技术一般都是针对贴体网格方法提出，ib湍流模拟中的壁面模化已引起了广泛的关注。然而，目前针对心脏瓣膜血液流动的数值模拟，均采用基于平衡
技术实现思路
1、本专利技术针对现有技术中的不足，提供一种基于非平衡壁面模型的心脏瓣膜血液流动les/ib计算方法，实现心脏瓣膜血液流动的复杂流场的准确预测。
2、为实现上述目的，本专利技术采用以下技术方案：
3、一种基于非平衡壁面模型的心脏瓣膜血液流动les/ib计算方法，包括：
4、步骤1、根据给定的血管壁面的形状，生成流场域的背景贴体曲面euler网格，建立全局直角坐标系；
5、步骤2、将机械瓣膜表面进行三角离散，生成瓣膜lagrange网格，将瓣膜lagrange网格各点的初始位置设置为瓣膜瓣叶关闭状态的位置，初始速度设置为零，并浸没在背景贴体曲面euler网格中；
6、步骤3、设置背景贴体曲面euler网格的出口和入口边界条件以及其余边界的边界条件，设置初始血液流场解；
7、步骤4、应用流固强耦合方法，以初始血液流场解、瓣膜lagrange网格各点的初始位置和初始速度为基础，考虑边界条件，迭代求解每一个时间步上的血液流场解、瓣膜位置解和瓣膜速度解。
8、为优化上述技术方案，采取的具体措施还包括：
9、进一步地，步骤3具体为：
10、背景贴体曲面euler网格的右边界为出口，采用对流边界条件，uout是出口的流向平均速度，i取值为1、2、3，u1、u2、u3分别为全局直角坐标下的流向速度、法向速度及展向速度；t和x分别表示时间和全局直角坐标系的流向；背景贴体曲面euler网格左边界为入口，入口血液流量为给定的随时间变化的时间函数，初始血液流场解w0＝(p0，u0)为静止流场的解，p0表示初始时刻血液的压强，u0表示初始时刻血液的流动速度矢量；背景贴体曲面euler网格的其余边界设为无滑移边界条件，背景贴体曲面euler网格节点上u1＝0，u2＝0，u3＝0。
11、进一步地，步骤4具体为：基于第n时间步上血液流场解wn＝(pn，un)、瓣膜lagrange网格各点的位置解xn和速度解vn，n＝0～nmax-1，应用强耦合方法求n+1时间步上血液流场解wn+1、瓣膜lagrange网格各点的位置解xn+1和瓣膜速度解vn+1，nmax为最大时间步，pn表示第n时间步上血液的压强，un表示第n时间步上血液的速度矢量，包括以下步骤：
12、步骤4.1、设置时间步n内的子迭代步序号k＝1，令wk＝wn，xk＝xn，vk＝vn，n＝0～nmax-1；wk为第k子迭代步上血液流场解，xk和vk分别为第k子迭代步瓣膜lagrange网格上各点的位置解和速度解；wn＝(pn，un)，wn表示第n时间步的血液流场解，pn表示第n时间步上血液的压强，un表示n时间步上血液的速度矢量，xn表示第n时间步的瓣膜lagrange网格上各点的位置解，vn表示第n时间步的瓣膜lagrange网格上各点的速度解；nmax表示最大时间步；
13、步骤4.2、根据第k子迭代步瓣膜lagrange网格上各点的位置解xk，将背景贴体曲面euler网格点分为解域内部计算点、外部相关点和外部不参与计算的无关点；
14、步骤4.3、根据第k子迭代步瓣膜lagrange网格上各点的位置解xk及对应的速度解vk，利用边界条件及法向积分的湍流边界层方程，求解非平衡壁面模型，获得当地摩擦速度；
15、步骤4.4、利用当地摩擦速度计算瓣膜lagrange网格上各点的壁面剪切应力；利用壁面剪切应力计算外部相关点上的当地流向速度根据第k子迭代步瓣膜lagrange网格上各点的速度解vk、无穿透边界条件和沿瓣膜壁面法向的线性插值求解当地法向动量方程，获得第k子迭代步外部相关点di上的当地法向速度和压强通过当地坐标系与全局直角坐标系之间的转换关系，将当地流向速度和当地法向速度转换到全局直角坐标系上，得到第k子迭代步外部相关点di在全局直角坐标下的速度矢量获得外部相关点di上的第k子迭代步血液流场解
16、步骤4.5、利用三角线性插值确定当前第k子迭代步参考点fi上的湍流粘性系数然后将第k子迭代步外部相关点di上的湍流粘性系数设为与第k子迭代步参考点fi上的湍流粘性系数相等；
17、步骤4.6、基于外部相关点di上的第k子迭代步血液流场解和湍流粘性系数通过虚拟时间推进获得第k+1迭代步上内部计算点上流动方程的解wk+1；
18、步骤4.7、确定瓣膜lagrange网格各点所在的euler网格单元，利用三角线性插值获得瓣膜lagrange网格点上的流体力；
19、步骤4.8、根据瓣膜lagrange网格点上的流体力求解瓣膜的结构运动方程，得到瓣膜打开角度，所述瓣膜的结构运动方程如下：
20、
21、式中，ired为瓣膜约化的转动惯量，θk+1表示第k+1子迭代步上瓣膜打开角度，t为时间，mi为瓣膜表面所受转动矩的积分；
22、采用亚松弛方法优化当前第k+1子迭代步的瓣膜打开角度θk+1，得到第k+1子迭代步上松弛后的瓣膜打开角度
23、通过松弛后的瓣膜打开角度确定第k+1迭代步上瓣膜lagrange网格各点的位本文档来自技高网...
【技术保护点】
1.一种基于非平衡壁面模型的心脏瓣膜血液流动LES/IB计算方法，其特征在于，包括：
2.如权利要求1所述的基于非平衡壁面模型的心脏瓣膜血液流动LES/IB计算方法，其特征在于，步骤3具体为：
3.如权利要求1所述的基于非平衡壁面模型的心脏瓣膜血液流动LES/IB计算方法，其特征在于，步骤4具体为：基于第n时间步上血液流场解wn＝(pn，un)、瓣膜Lagrange网格各点的位置解Xn和速度解Vn，n＝0～Nmax-1，应用强耦合方法求n+1时间步上血液流场解wn+1、瓣膜Lagrange网格各点的位置解Xn+1和瓣膜速度解Vn+1，Nmax为最大时间步，pn表示第n时间步上血液的压强，un表示第n时间步上血液的速度矢量，包括以下步骤：
4.如权利要求3所述的基于非平衡壁面模型的心脏瓣膜血液流动LES/IB计算方法，其特征在于，步骤4.2具体为：
5.如权利要求3所述的基于非平衡壁面模型的心脏瓣膜血液流动LES/IB计算方法，其特征在于，步骤4.3具体为：
6.如权利要求3所述的基于非平衡壁面模型的心脏瓣膜血液流动LES
...【技术特征摘要】
1.一种基于非平衡壁面模型的心脏瓣膜血液流动les/ib计算方法，其特征在于，包括：
2.如权利要求1所述的基于非平衡壁面模型的心脏瓣膜血液流动les/ib计算方法，其特征在于，步骤3具体为：
3.如权利要求1所述的基于非平衡壁面模型的心脏瓣膜血液流动les/ib计算方法，其特征在于，步骤4具体为：基于第n时间步上血液流场解wn＝(pn，un)、瓣膜lagrange网格各点的位置解xn和速度解vn，n＝0～nmax-1，应用强耦合方法求n+1时间步上血液流场解wn+1、瓣膜lagrange网...
【专利技术属性】
技术研发人员：濮天梅，王静洋，张扬，徐振东，付士凤，祁浩天，
申请(专利权)人：扬州大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。