颈动脉粥样硬化斑块模型的谱元法计算-.pdfVIP

第32卷 第4期 工 程 数 学 学 报 Vol. 32 No. 4 2015年08月 CHINESE JOURNAL OF ENGINEERING MATHEMATICS Aug. 2015 doi: 10.3969/j.issn.1005-3085.2015.04.008 文章编号: 1005-3085(2015)04-0546-11 颈动脉粥样硬化斑块模型的谱元法计算∗ 宋方应, 刘晓玲, 许传炬 (厦门大学数学科学学院，厦门 361005) 摘 要: 本文考虑人体颈动脉粥样硬化斑块的力学建模及其数值计算．颈动脉粥样硬化斑块的 动力学行为由控制血液流动的不可压缩Navier-Stokes 方程和控制动脉形变的弹性方程 耦合描述，因此完整的硬化斑块模型是一组三维流体/固体耦合方程，其数学理论和数 值计算都面临极大的困难．但在一定的假设条件下，该问题可简化为一个弹性力学方 程．本文考虑粥样硬化斑块简化模型的数值计算问题．我们构造了一个针对模型简化 方程的有效算法，即基于空间方向上的谱元离散和时间方向的Newmark 格式的计算方 法．分析显示前者对光滑解具有指数收敛性，而后者具有二阶收敛精度并且在特定的 参数条件下无条件稳定．论文不仅讨论了连续问题以及半离散问题的稳定性，而且给 出了全离散格式的最优误差估计．最后通过数值模拟验证了所提算法的有效性和理论 分析结果，并对一些具有实际背景参数的模型问题进行了模拟． 关键词: 动脉粥样硬化斑块；弹性方程；谱元法；数值分析 分类号: (2000) 74F10; 65M70; 65M12 中图分类号: O302; O34; O241 文献标识码: A 引引引言言言 心脑血管疾病是威胁人们健康的一大疾病，当前主要是通过医学图像(例如：核磁共 振成像(MRI)、X 射线扫描(CT)、超声) 等手段进行诊断．医学上假设冠心病、脑中风等 疾病主要是由于血管阻塞引起的脑部供氧不足引起的，而血管粥样硬化斑块的形成是造成 血管阻塞的主要原因．最近，计算方法的发展使得通过数值模拟研究病变过程成为可能． 数值模拟正成为医生对这类疾病进行诊断的有效工具． 人体血管粥样硬化过程的动力学行为可由控制血液流动的不可压缩Navier-Stokes 方 程和控制动脉形变的弹性方程耦合描述，但完整的硬化斑块模型是一组三维流体/固体耦 合方程，数值计算非常困难．一些研究发现，在一定的假设条件下斑块的力学模型可简 [1] [1] 化为一个弹性力学方程 ，还出现了一些求解该简化模型的数值方法，如无网格法 、有 [2,3] 限差分法和有限体积 等．考虑到该问题有复杂的计算区域和相对较好的正则性，本文 提出一个具有高阶精度特征的空间谱元法和较好稳定性的时间Newmark 推进格式．一方 面，谱方法与低阶方法相比可用较少的节点数达到更高的精度；另一方面，谱元法的区域 剖分技术又允许灵活处理复杂的计算区域；还有，谱元法可扩展性强更容易实现大规模并 收稿日期: 2014-03-17. 作者简介: 宋方应(1985年9月生)，男，博士. 研究方向：微分方程数值解及并行计算. 基金项目: 国家自然科学基金 11421110001;. 第4期 宋方应，等：颈动脉粥样硬化斑块模型的谱元法计算 547 行计算．本文可视为硬化斑块力学建模及其数值计算研究的第一步，旨在为将来实现三维 流固耦合模型的高精度模拟提供基础算法．本文的主要结构如下，在接下来的一节中我们 将给出颈动脉粥样动脉硬化的二维弹性动力学模型；在第三节我们讨论该模型的数值方法 及理论分析，主要包括问题的稳定性分析以及问题的离散形式和误差分析；最后给出一些 数值算例来验证分析所得的结论，并提供一个具有实际背景参数的模型问题的模拟