离心式血泵内部流场数值模拟及血损特性研究：方法、分析与优化

离心式血泵内部流场数值模拟及血损特性研究：方法、分析与优化一、引言1.1研究背景与意义心力衰竭是一种严重的心血管疾病，近年来，其发病率和死亡率呈上升趋势，给全球医疗健康带来了沉重负担。据统计，全球约有2600万人受心力衰竭影响，在我国≥25岁人群心衰患病率达1.1%，约为1210万人，每年新发心衰患者297万人。心衰患者不仅面临着呼吸困难、疲劳、水肿等症状的折磨，生活质量严重下降，而且其5年生存率与恶性肿瘤相当，严重威胁着患者的生命安全。传统的药物治疗对于终末期心力衰竭患者效果有限，而心脏移植作为目前治疗终末期心力衰竭的有效手段，由于供体器官短缺，远远无法满足患者的需求，大量患者在等待心脏移植的过程中死亡。在此背景下，机械血泵作为一种重要的替代治疗方案应运而生，为心力衰竭患者带来了新的希望。血泵能够部分或全部替代心脏的泵血功能，帮助衰竭心脏完成血液循环，从而减轻心脏负担，改善患者的病情。按照工作原理，血泵可分为搏动型和非搏动型，其中离心式血泵作为非搏动型血泵的一种，具有体积小、流量大、效率高等优点，在临床上得到了广泛应用。例如，HeartMateⅢ作为目前全球科技含量较高的心室辅助装置，采用了先进的磁悬浮技术，治疗终末期心衰2年生存率可达79%，已媲美心脏移植，同时显著降低了脑卒中事件和机械故障率。国内阜外医院团队成功转化具有自主知识产权的长期可植入式全磁悬浮心室辅助装置CH-VAD血泵，成功救治了多例终末期心衰合并心原性休克患者，其中1例已携带血泵生存超过5年，无出血、血栓形成、感染等并发症。然而，血泵在临床应用中仍面临一些关键问题，其中血栓和溶血问题尤为突出。血泵运行过程中产生的血栓和溶血超出安全范围，将会引发多种严重并发症，如血栓栓塞可导致脑卒中、肺栓塞等，溶血则可能引起贫血、肾功能损害等，严重者甚至危及病人生命。研究表明，血栓和溶血的产生与血泵内部的流场特性密切相关，复杂的流场容易导致血液受到过高的剪切力、产生流动分离和漩涡等现象，从而损伤血细胞，增加血栓形成的风险。因此，深入研究离心式血泵内部的流场特性，对于优化血泵设计、降低血栓和溶血风险、提高血泵性能和安全性具有至关重要的意义。数值模拟技术的快速发展为离心式血泵的研究提供了强有力的工具。通过数值模拟，可以在计算机上构建血泵的三维模型，模拟不同工况下血泵内部的血液流动情况，获取详细的流场信息，如速度分布、压力分布、壁面剪切应力分布等。与传统的实验研究方法相比，数值模拟具有成本低、周期短、可重复性强等优势，可以快速对不同设计方案进行评估和优化，大大提高了研究效率。例如，通过数值模拟可以分析叶轮形状、叶片数量、转速等因素对血泵内部流场和血损的影响，为血泵的结构设计和参数优化提供科学依据，有助于开发出性能更优、血液相容性更好的离心式血泵，从而推动血泵技术的发展，为心力衰竭患者提供更有效的治疗手段。1.2国内外研究现状1.2.1血泵发展历程血泵的发展经历了多个重要阶段，从早期的搏动式血泵到如今的离心式血泵，每一次技术革新都为心力衰竭患者带来了更多的希望。20世纪60年代，第一代搏动式血泵诞生，它通过模仿心脏的搏动方式来推动血液流动，为心脏手术患者提供了重要的辅助支持。然而，这种血泵存在体积庞大、易感染、血栓形成和溶血等问题，限制了其长期应用。例如，早期的气动搏动泵需要大型的外部驱动装置，患者行动极为不便，且感染风险较高。随着科技的不断进步，20世纪80年代，第二代轴流血泵应运而生。轴流血泵采用了轴流叶轮设计，通过叶片的旋转推动血液沿轴向流动，具有体积小、流量大等优点。与第一代血泵相比，轴流血泵在性能上有了显著提升，但由于机械轴承的存在，容易出现磨损，导致耐久性降低，同时血栓和溶血问题仍然较为突出。比如，一些早期的轴流血泵在运行过程中，机械轴承的磨损会导致血泵性能下降，增加了患者的治疗风险。为了解决轴流血泵的问题，20世纪90年代以来，第三代离心式血泵逐渐成为研究和应用的热点。离心式血泵利用叶轮的高速旋转产生离心力，使血液在离心力的作用下从叶轮中心流向周边，实现血液的高效输送。与前两代血泵相比，离心式血泵具有体积小、流量大、效率高、血液相容性好等优势，能够更好地满足临床需求。例如，HeartMateⅢ作为目前全球先进的心室辅助装置，采用了先进的磁悬浮技术，有效减少了血液与泵体的接触摩擦，显著降低了血栓形成和溶血的风险，其治疗终末期心衰2年生存率可达79%，已媲美心脏移植，同时显著降低了脑卒中事件和机械故障率。在国内，血泵的研究起步相对较晚，但近年来取得了显著的进展。广东省心血管病研究所在20世纪90年代开始研制第一代气动隔膜泵—罗叶式泵，并在重症心衰患者中进行了短期应用。此后，国内多家科研机构和医院致力于血泵的研发，如江苏理工大学、北京安贞医院和北京阜外心血管病医院等对第二代血泵的研发进行了深入探索。近年来，国内在第三代离心式血泵的研发方面取得了重要突破，阜外医院团队成功转化具有自主知识产权的长期可植入式全磁悬浮心室辅助装置CH-VAD血泵，已成功救治多例终末期心衰合并心原性休克患者，其中1例已携带血泵生存超过5年，无出血、血栓形成、感染等并发症。1.2.2数值模拟技术数值模拟技术在血泵研究中的应用越来越广泛，为血泵的设计优化和性能评估提供了重要的手段。计算流体力学（CFD）作为数值模拟的核心技术，能够对血泵内部的血液流动进行精确模拟，揭示流场的复杂特性。通过CFD模拟，可以获取血泵内部的速度分布、压力分布、壁面剪切应力分布等重要信息，为分析血泵的性能和血液损伤机制提供了依据。例如，学者们利用CFD技术对不同叶轮形状和叶片数量的离心式血泵进行模拟，分析了这些因素对血泵内部流场的影响，发现叶轮形状和叶片数量的优化可以改善流场分布，降低血液受到的剪切力。在国外，数值模拟技术在血泵研究中得到了深入应用。美国、日本等国家的科研团队通过数值模拟对新型血泵的设计进行了大量研究，取得了一系列重要成果。例如，美国的研究人员利用CFD技术对磁悬浮离心血泵进行模拟，优化了血泵的结构参数，提高了血泵的效率和血液相容性。日本的学者则通过数值模拟研究了血泵内部的流动特性与血栓形成的关系，为降低血栓风险提供了理论支持。国内在血泵数值模拟方面也取得了一定的成果。许多高校和科研机构开展了相关研究，利用CFD技术对不同类型的血泵进行模拟分析。上海医疗器械高等专科学校的张伟国等人通过与医院合作实验采集数据，对叶片式离心血泵内部流场进行数值模拟，对比了血泵实际运行情况与数值计算结果，对其内部血栓和溶血问题进行了系统分析研究，发现数值模拟分析的情况与血泵在实际运行中的血栓和溶血情况基本相符。北京工业大学的常宇教授等人对人工心脏泵进行了一系列研究，通过对磁流血泵柔性转子建模仿真，满足了患者的生理活动需求。1.2.3血损实验研究血损实验研究是评估血泵性能的重要手段，通过实验可以直接测量血泵运行过程中对血液的损伤程度，为血泵的优化设计提供实验依据。在血损实验中，常用的指标包括溶血指数、血小板激活程度、血栓形成情况等。溶血指数反映了血泵对红细胞的破坏程度，血小板激活程度则反映了血泵对血小板的影响，血栓形成情况直接关系到患者的治疗安全。例如，通过测量溶血指数，可以评估血泵在不同工况下对红细胞的损伤程度，从而确定血泵的最佳运行参数。国外在血损实验研究方面开展了大量工作，建立了完善的实验方法和评价体系。许多研究机构通过体外实验和动物实验，对不同类型血泵的血损情况进行了深入研究。美国的一些研究团队通过体外循环实验，对比了不同血泵在相同工况下的血损指标，发现离心式血泵在血液相容性方面优于其他类型血泵。德国的科研人员则通过动物实验，研究了血泵长期植入对动物生理指标的影响，为血泵的临床应用提供了重要参考。国内的血损实验研究也在不断发展，科研人员通过改进实验方法和设备，提高了实验的准确性和可靠性。一些高校和医院建立了专门的血泵实验平台，开展了血损相关的实验研究。天津大学的研究团队利用自主搭建的血泵实验平台，对离心式血泵的血损特性进行了实验研究，分析了叶轮转速、流量等因素对血损的影响，为血泵的优化设计提供了实验依据。同时，国内还加强了与国外的合作与交流，引进了先进的实验技术和设备，进一步推动了血损实验研究的发展。1.3研究内容与方法1.3.1研究内容本文旨在通过数值模拟与实验研究相结合的方法，深入探究离心式血泵内部流场特性及其与血损之间的关系，为离心式血泵的优化设计提供理论依据和技术支持。具体研究内容如下：建立离心式血泵的数值模型：利用三维建模软件（如SolidWorks、Pro/E等），根据离心式血泵的实际结构参数，构建精确的三维几何模型，包括泵体、叶轮、进口和出口等部件。对建立的三维模型进行网格划分，采用合适的网格划分策略（如结构化网格、非结构化网格或混合网格），确保在关键部位（如叶轮叶片表面、流道狭窄处等）具有足够的网格密度，以准确捕捉流场细节。同时，为了提高计算效率，合理控制网格数量，避免因网格过多导致计算资源浪费和计算时间过长。将划分好的网格导入计算流体力学（CFD）软件（如ANSYSFluent、COMSOLMultiphysics等），根据血液的物理特性，选择合适的流体模型（如非牛顿流体模型），并设置相应的边界条件（如进口流量、出口压力、壁面无滑移条件等）和初始条件，建立完整的数值计算模型。数值模拟离心式血泵内部流场：运用CFD软件对建立的数值模型进行求解计算，模拟不同工况下（如不同转速、流量等）离心式血泵内部的血液流动情况，获取流场的详细信息，包括速度分布、压力分布、流线分布等。分析不同工况下流场的变化规律，研究叶轮转速、流量等因素对血泵内部流场的影响，探讨流场特性与血泵性能之间的关系。例如，通过对比不同转速下的流场模拟结果，分析转速变化对流场速度分布和压力分布的影响，进而研究其对血泵扬程和效率的影响。此外，还将研究流场中的流动分离、漩涡等复杂流动现象，分析其产生的原因和对血泵性能的影响。分析离心式血泵内部流场与血损的关系：根据流场模拟结果，计算血泵内部的壁面剪切应力分布、剪切力作用时间等参数，分析这些参数与血损（如溶血、血栓形成等）之间的关系。建立血损预测模型，基于数值模拟得到的流场参数，结合相关的血损理论和实验数据，预测血泵在不同工况下的血损程度。例如，利用溶血模型，根据流场中的剪切应力和作用时间，预测红细胞的损伤程度，评估血泵的溶血风险。通过分析流场与血损的关系，找出影响血损的关键流场因素，为血泵的优化设计提供理论依据。实验研究验证数值模拟结果：搭建离心式血泵实验平台，包括血泵装置、驱动系统、流量测量系统、压力测量系统等，确保实验平台能够准确模拟血泵的实际运行工况。采用粒子图像测速（PIV）技术、激光多普勒测速（LDV）技术等先进的实验测量方法，测量血泵内部流场的速度分布、压力分布等参数，与数值模拟结果进行对比验证。同时，通过实验测量血泵运行过程中的血损指标（如溶血指数、血小板激活程度等），验证血损预测模型的准确性。根据实验结果，对数值模拟模型进行修正和完善，提高数值模拟的精度和可靠性。1.3.2研究方法文献研究法：广泛查阅国内外有关离心式血泵的研究文献，包括学术期刊论文、学位论文、专利文献、技术报告等，了解离心式血泵的发展历程、研究现状、数值模拟方法、血损实验研究等方面的最新成果和研究动态。对文献进行系统梳理和分析，总结前人的研究经验和不足之处，为本研究提供理论基础和研究思路。数值模拟法：运用计算流体力学（CFD）技术，对离心式血泵内部流场进行数值模拟。CFD技术是一种基于数值计算方法求解流体力学控制方程的技术，能够对复杂的流体流动现象进行精确模拟。通过建立离心式血泵的数值模型，设置合理的边界条件和初始条件，利用CFD软件进行求解计算，可以获取血泵内部流场的详细信息，为流场分析和血损研究提供数据支持。数值模拟法具有成本低、周期短、可重复性强等优点，可以快速对不同设计方案进行评估和优化，大大提高了研究效率。实验研究法：搭建离心式血泵实验平台，进行实验研究。实验研究法是一种通过实际实验操作获取数据和验证理论的方法，能够直接测量血泵内部流场的参数和血损指标，具有直观、准确的优点。通过实验测量得到的数据可以与数值模拟结果进行对比验证，评估数值模拟的准确性和可靠性。同时，实验研究还可以发现一些数值模拟难以捕捉到的现象和问题，为数值模拟模型的改进和完善提供依据。对比分析法：对数值模拟结果和实验研究结果进行对比分析，评估数值模拟模型的准确性和可靠性。对比分析不同工况下数值模拟和实验测量得到的流场参数和血损指标，找出两者之间的差异和原因。根据对比分析结果，对数值模拟模型进行修正和优化，提高数值模拟的精度和可靠性。此外，还将对比分析不同设计方案的血泵性能和血损情况，为血泵的优化设计提供参考依据。二、离心式血泵工作原理及数值模拟基础2.1离心式血泵结构与工作原理离心式血泵主要由泵体、叶轮、进口、出口以及驱动装置等部件构成，其结构设计紧密围绕血液输送的高效性与安全性展开。泵体作为血泵的外壳，为内部部件提供了一个封闭且稳定的工作空间，其形状和尺寸经过精心设计，以确保血液在其中能够顺畅流动，减少流动阻力和能量损失。叶轮则是离心式血泵的核心部件，通常由多个叶片组成，这些叶片均匀分布在叶轮的圆周上，叶片的形状、数量和角度对血泵的性能有着至关重要的影响。例如，叶片的形状一般采用流线型设计，以减少血液在叶片表面的流动分离和漩涡产生，从而降低血液受到的剪切力；叶片数量的选择需要综合考虑血泵的流量、扬程和效率等因素，过多或过少的叶片数量都可能导致血泵性能下降。进口和出口分别位于泵体的特定位置，进口负责引导血液流入血泵，其设计旨在使血液能够平稳地进入叶轮，避免产生冲击和紊流；出口则用于将经过叶轮加速后的血液输送到人体循环系统中，出口的尺寸和形状需要与人体血管的生理特征相匹配，以确保血液能够顺利地融入人体血液循环。驱动装置为叶轮的旋转提供动力，常见的驱动方式包括电动、磁驱动等，不同的驱动方式在效率、稳定性和血液相容性等方面存在差异。例如，磁驱动方式由于避免了机械接触，能够有效减少磨损和血栓形成的风险，提高了血泵的血液相容性和长期运行稳定性。离心式血泵的工作原理基于离心力的作用。当叶轮在驱动装置的带动下高速旋转时，血液从进口被吸入到叶轮中心。由于叶轮的高速旋转，血液在叶轮内受到强大的离心力作用，从叶轮中心沿叶片表面向叶轮周边加速流动。在这个过程中，血液获得了动能和压力能，速度和压力不断增加。随着血液向叶轮周边流动，其速度进一步提高，压力也相应增大，最终从出口被排出到人体循环系统中。叶轮的旋转速度和叶片形状决定了血液获得的能量大小，进而影响血泵的流量和扬程。较高的叶轮转速可以使血液获得更大的离心力，从而提高血泵的流量和扬程，但同时也可能增加血液受到的剪切力，导致血损增加。因此，在实际应用中，需要根据患者的具体需求和生理状况，合理选择血泵的转速和结构参数，以确保血泵既能提供足够的血液供应，又能将血损控制在安全范围内。例如，对于一些病情较轻的患者，可能不需要过高的血泵流量和扬程，此时可以适当降低叶轮转速，以减少血损；而对于病情较重、需要大量血液供应的患者，则需要提高叶轮转速，以满足患者的生理需求，但同时需要密切关注血损情况，采取相应的措施降低血损风险。2.2血液特性及血损现象血液是一种复杂的非牛顿流体，其流变学特性对于理解离心式血泵内部的流动行为和血损机制至关重要。从微观层面来看，血液主要由血浆和血细胞组成，血细胞包括红细胞、白细胞和血小板。红细胞在血液中含量最多，其独特的双凹圆盘状结构使其具有良好的变形能力，这对于血液在微血管中的流动至关重要。在正常生理条件下，红细胞能够在较小的切变率下发生聚集，形成缗钱状结构，这种聚集行为在一定程度上影响了血液的黏度。当切变率增加时，红细胞的聚集结构被破坏，红细胞发生变形，沿流线方向排列，血液黏度随之降低，呈现出剪切变稀的特性。例如，在大血管中，血液流速较快，切变率较高，红细胞能够快速变形并有序排列，使血液黏度降低，减少流动阻力；而在微血管中，切变率较低，红细胞更容易聚集，血液黏度相对较高。白细胞和血小板在血液中的含量相对较少，但它们在免疫反应和凝血过程中发挥着关键作用。白细胞参与机体的免疫防御，能够识别和清除病原体；血小板则在血管受损时迅速黏附、聚集，形成血小板血栓，起到止血的作用。血液中的血浆主要由水、蛋白质、电解质等组成，其中蛋白质（如白蛋白、球蛋白等）对血浆的黏度有重要影响。这些蛋白质分子之间的相互作用以及它们与血细胞之间的相互作用，共同决定了血液的流变学特性。在离心式血泵的运行过程中，血泵内部复杂的流场会对血液产生多种作用，从而导致溶血和血栓等血损现象的发生。溶血是指红细胞破裂，血红蛋白释放到血浆中的过程。当血液在血泵中流动时，如果受到过高的剪切应力作用，红细胞的细胞膜就会受到损伤，导致细胞破裂，发生溶血。研究表明，当剪切应力超过一定阈值（通常为100Pa-200Pa）且作用时间达到一定程度时，溶血风险会显著增加。例如，在叶轮叶片附近和流道狭窄处，由于流速梯度较大，血液受到的剪切应力较高，容易发生溶血。此外，流动分离和漩涡的产生也会增加血液与固体壁面的摩擦，进一步加剧红细胞的损伤，导致溶血。溶血不仅会降低血液的携氧能力，引发贫血等症状，还可能导致游离血红蛋白在血管内堆积，对肾脏等器官造成损害，严重时甚至危及生命。血栓形成是另一种严重的血损现象，其发生机制较为复杂，涉及血液凝固系统、血小板活化以及血管壁的损伤等多个因素。在离心式血泵中，血液与泵体表面的接触以及流场中的高剪切应力区域都可能引发血小板的活化和聚集。当血小板被激活后，它们会释放出多种生物活性物质，如血栓素A2、二磷酸腺苷（ADP）等，这些物质进一步促进血小板的聚集和黏附，形成血小板血栓。同时，血液中的凝血因子也会被激活，启动凝血级联反应，导致纤维蛋白原转化为纤维蛋白，形成纤维蛋白网络，将血小板和血细胞包裹其中，最终形成血栓。例如，在血泵的进口和出口处，由于血液流动的突然变化，容易产生流动分离和漩涡，这些区域的血液滞留时间较长，血小板和凝血因子更容易聚集和活化，增加了血栓形成的风险。血栓一旦形成，可能会脱落并随血流进入人体循环系统，导致血管栓塞，如脑栓塞、肺栓塞等，严重威胁患者的生命健康。2.3数值模拟基本方法2.3.1控制方程描述离心式血泵内血液流动的基本控制方程主要包括连续性方程、动量方程和能量方程。连续性方程基于质量守恒定律，反映了血液在流动过程中质量的连续性，其表达式为：\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{v})=0其中，\rho为血液密度，t为时间，\vec{v}为速度矢量。该方程表明，单位时间内控制体内血液质量的变化率与通过控制体表面的质量通量之和为零，即血液在流动过程中既不会凭空产生也不会无端消失，保证了质量的守恒。动量方程依据牛顿第二定律，揭示了血液在流动过程中动量的变化规律，其一般形式为：\frac{\partial(\rho\vec{v})}{\partialt}+\nabla\cdot(\rho\vec{v}\vec{v})=-\nablap+\nabla\cdot\tau+\rho\vec{g}式中，p为压力，\tau为应力张量，\vec{g}为重力加速度。该方程体现了单位时间内控制体内血液动量的变化率等于作用在控制体上的压力、粘性力和重力之和。在离心式血泵中，血液受到叶轮旋转产生的离心力以及泵体壁面的约束，这些力的作用通过动量方程得以体现，从而影响血液的流动速度和方向。能量方程遵循能量守恒原理，描述了血液在流动过程中能量的转换和守恒关系，对于不可压缩流体，能量方程可简化为：\frac{\partial(\rhoe)}{\partialt}+\nabla\cdot(\rho\vec{v}e)=\nabla\cdot(k\nablaT)+\Phi其中，e为单位质量血液的内能，k为热传导系数，T为温度，\Phi为粘性耗散函数。在血泵的实际运行中，血液与泵体部件之间的摩擦会导致能量损失，主要表现为机械能转化为热能，能量方程能够准确描述这一能量转换过程，对于分析血泵的热效应和能量效率具有重要意义。例如，在高速旋转的叶轮与血液之间的摩擦生热过程中，能量方程可以帮助我们计算出温度的变化，从而评估热效应对血液性质和血泵性能的影响。2.3.2离散方法有限体积法是数值模拟中常用的一种离散方法，它将计算区域划分为一系列互不重叠的控制体积，在每个控制体积内对控制方程进行积分离散。以连续性方程为例，对其在控制体积V上进行积分：\int_{V}\frac{\partial\rho}{\partialt}dV+\int_{S}\rho\vec{v}\cdotd\vec{S}=0利用高斯散度定理，将面积分转换为体积分，再采用合适的插值函数对控制体积内的物理量进行离散近似，从而将连续的控制方程转化为代数方程组。在离散过程中，选择合适的插值函数至关重要，常用的插值函数有线性插值、二次插值等。线性插值函数形式简单，计算效率高，但对于复杂流场的描述精度相对较低；二次插值函数能够更好地拟合流场的变化，但计算复杂度有所增加。例如，在处理离心式血泵内部复杂的流场时，对于流场变化较为平缓的区域，可以采用线性插值函数以提高计算效率；而在流场变化剧烈的区域，如叶轮叶片附近和流道狭窄处，则采用二次插值函数，以更准确地捕捉流场的细节。对于动量方程和能量方程，同样采用有限体积法进行离散。在离散动量方程时，需要对速度、压力等变量进行离散处理，常用的离散格式有中心差分格式、迎风格式等。中心差分格式在计算精度上具有一定优势，但对于对流占主导的流动问题，容易产生数值振荡；迎风格式则根据流动方向选择上游节点的信息进行插值，能够有效抑制数值振荡，提高计算的稳定性，但在精度上可能会有所损失。在实际应用中，需要根据具体问题的特点，综合考虑计算精度和稳定性的要求，选择合适的离散格式。例如，在模拟离心式血泵内部高速旋转的叶轮附近的流动时，由于流动方向明确且对流作用较强，采用迎风格式可以更好地保证计算的稳定性；而在一些流动相对平稳的区域，可以采用中心差分格式以提高计算精度。有限体积法的优点在于离散后的方程具有明确的物理意义，能够较好地满足守恒定律，而且在处理复杂边界条件时具有较强的灵活性。通过合理划分控制体积和选择离散格式，可以在保证计算精度的前提下，提高计算效率，使得有限体积法在离心式血泵的数值模拟中得到了广泛应用。2.3.3湍流模型在离心式血泵内部，血液的流动通常处于湍流状态，选择合适的湍流模型对于准确模拟流场特性至关重要。常见的湍流模型包括标准k-\varepsilon模型、RNGk-\varepsilon模型、Realizablek-\varepsilon模型和SSTk-\omega模型等。标准k-\varepsilon模型是一种基于涡粘性假设的两方程湍流模型，通过求解湍动能k和湍动能耗散率\varepsilon的输运方程来封闭雷诺应力项。该模型计算效率较高，在工程领域应用广泛，但它在处理复杂流动，如强旋转、弯曲壁面等问题时，存在一定的局限性。例如，在离心式血泵中，叶轮的高速旋转会产生强烈的离心力和科里奥利力，标准k-\varepsilon模型难以准确描述这种复杂的流动特性，可能导致计算结果与实际情况存在较大偏差。RNGk-\varepsilon模型在标准k-\varepsilon模型的基础上，引入了重整化群理论，对湍动能耗散率方程进行了修正，使其能够更好地处理高应变率和流线弯曲程度较大的流动。相比标准k-\varepsilon模型，RNGk-\varepsilon模型在模拟离心式血泵内部流场时，对于叶轮附近的复杂流动有更好的适应性，能够更准确地捕捉到流场中的漩涡和流动分离现象。然而，该模型在近壁区域的计算精度仍有待提高，需要结合合适的壁面函数来处理边界条件。Realizablek-\varepsilon模型对湍流粘性系数和湍动能耗散率方程进行了改进，使其在理论上满足雷诺应力的数学约束，能够更准确地预测一些复杂流动的特性，如射流、分离流等。在离心式血泵的数值模拟中，Realizablek-\varepsilon模型对于流道内的流动分离和再附着现象的模拟具有一定的优势，能够为血泵性能的分析提供更准确的流场信息。但该模型的计算复杂度相对较高，对计算资源的要求也较高。SSTk-\omega模型结合了k-\omega模型在近壁区域的优势和k-\varepsilon模型在远场的优势，通过一个混合函数在近壁区域和远场之间进行平滑过渡。该模型对压力梯度敏感，能够准确模拟逆压梯度下的流动分离，在模拟离心式血泵内部流场时，对于泵体壁面附近的流动以及叶轮与泵体之间的间隙流动有较好的模拟效果。同时，SSTk-\omega模型在计算精度和计算效率之间取得了较好的平衡，因此在血泵数值模拟中得到了越来越广泛的应用。例如，在研究离心式血泵的血液相容性时，需要准确了解壁面附近的流场特性，SSTk-\omega模型能够提供更准确的壁面剪切应力分布等信息，有助于分析血泵运行过程中对血液的损伤情况。在选择湍流模型时，需要综合考虑血泵内部流场的特点、计算精度要求以及计算资源等因素。对于离心式血泵内部复杂的三维流动，SSTk-\omega模型通常能够提供更准确的模拟结果，但如果对计算效率要求较高，且流场相对简单，标准k-\varepsilon模型或其改进版本也可以作为一种选择。通过对不同湍流模型的对比分析和验证，选择最适合离心式血泵内部流场模拟的湍流模型，对于深入研究血泵性能和血损机制具有重要意义。2.4网格划分与边界条件设置在对离心式血泵进行数值模拟时，准确构建计算域并进行合理的网格划分是获取可靠结果的关键步骤。计算域的建立需完整涵盖血泵内部的所有流道区域，包括泵体内部空间、叶轮旋转区域以及进出口管道等，以确保能够全面捕捉血液在血泵内的流动特性。利用专业的三维建模软件（如SolidWorks、Pro/E等），依据离心式血泵的精确设计图纸和实际结构参数，构建出高精度的三维几何模型。在建模过程中，对泵体、叶轮、进口、出口等关键部件的几何形状和尺寸进行严格把控，保证模型的准确性。例如，叶轮叶片的形状和角度对血泵内部流场有着重要影响，在建模时需精确还原叶片的复杂曲面，以准确模拟血液与叶片之间的相互作用。完成三维几何模型构建后，采用合适的网格划分策略对计算域进行离散处理。对于离心式血泵这种包含复杂几何形状和运动部件的模型，通常选用非结构化网格或混合网格划分方法。非结构化网格能够更好地适应复杂的几何形状，在叶轮叶片表面、流道狭窄处等关键部位，可以灵活地生成高质量的网格，确保在这些区域具有足够的网格密度，从而准确捕捉流场的细节信息。例如，在叶轮叶片表面，通过局部加密网格，可以更精确地计算叶片表面的压力分布和壁面剪切应力，为分析血泵的性能和血损情况提供更准确的数据。而混合网格则结合了结构化网格和非结构化网格的优点，在流道较为规则的区域采用结构化网格，以提高计算效率；在几何形状复杂的区域采用非结构化网格，以保证计算精度。在网格划分过程中，需要对网格质量进行严格检查和优化。网格质量的好坏直接影响到数值模拟的精度和稳定性，常见的网格质量指标包括网格纵横比、雅克比行列式、网格扭曲度等。通过调整网格生成参数，如网格尺寸、网格增长率等，确保网格质量满足数值模拟的要求。例如，控制网格纵横比在合理范围内，避免出现过于狭长的网格，以防止数值计算过程中出现数值振荡和计算误差。同时，对生成的网格进行可视化检查，观察网格在关键部位的分布情况，及时发现并修正网格质量问题。边界条件的设置对于数值模拟结果的准确性同样至关重要。在离心式血泵的数值模拟中，主要涉及进口、出口、壁面等边界条件的设置。进口边界条件通常采用速度进口或质量流量进口，根据实际工况和研究需求，确定进口处血液的流速或质量流量。例如，在模拟正常生理状态下的血泵运行时，可根据人体心脏的生理参数，设定进口质量流量为5L/min，以模拟人体正常的血液输入情况。同时，还需考虑进口处血液的温度、密度等物理参数，确保进口条件的合理性。出口边界条件一般选择压力出口，根据血泵的工作环境和实际应用场景，设定出口处的压力值。在模拟血泵向人体动脉供血时，可参考人体动脉血压的生理范围，将出口压力设定为100mmHg-120mmHg，以保证模拟结果符合实际生理情况。此外，还需考虑出口处的背压变化对血泵性能的影响，通过合理设置出口边界条件，准确模拟血泵在不同背压下的工作状态。壁面边界条件采用无滑移边界条件，即认为血液与泵体壁面和叶轮表面之间不存在相对滑动，壁面上的速度为零。这一假设符合实际情况，能够准确反映血液在壁面附近的流动特性。同时，考虑到血液与壁面之间的相互作用，还需设置壁面的粗糙度和传热系数等参数。例如，对于表面光滑的泵体壁面，可将粗糙度设置为较小的值；而对于可能存在血栓附着或表面损伤的壁面，可适当增大粗糙度，以研究壁面粗糙度对血泵内部流场和血损的影响。在考虑血泵运行过程中的热效应时，需合理设置壁面的传热系数，以准确模拟血液与壁面之间的热量传递过程。通过合理的网格划分和准确的边界条件设置，能够建立起高精度的离心式血泵数值模型，为后续的流场模拟和血损分析提供可靠的基础。在实际模拟过程中，还需根据具体情况对网格划分和边界条件进行进一步的优化和调整，以确保数值模拟结果的准确性和可靠性。三、离心式血泵内部流场数值模拟结果与分析3.1稳态工况下流场特性3.1.1速度场分析为了深入了解离心式血泵在不同工况下的内部流场特性，首先对速度场进行详细分析。图1展示了在设计转速为2500r/min、流量为5L/min的工况下，血泵内部速度分布云图。从整体上看，血泵内部速度分布呈现出明显的规律性，在叶轮区域，速度分布具有显著的特征。叶轮中心处，由于血液刚刚进入叶轮，尚未获得足够的能量，速度相对较低，约为0.5m/s。随着叶轮的高速旋转，血液在离心力的作用下被迅速加速，沿着叶片表面向叶轮边缘流动，速度逐渐增大。在叶轮边缘处，速度达到最大值，约为5m/s，这是因为血液在叶轮旋转过程中获得了较大的动能。蜗壳区域的速度分布则相对较为复杂。靠近叶轮出口处，由于高速流出的血液直接冲击蜗壳壁面，速度仍然较高，约为3-4m/s，随后，血液在蜗壳内逐渐减速，速度逐渐降低至1-2m/s。这是由于蜗壳的横截面积逐渐增大，根据连续性方程，流速会相应减小。同时，蜗壳内的流动还受到叶轮出口射流和蜗壳壁面的影响，导致速度分布不均匀，在蜗壳的某些区域可能出现速度较低的回流区，这可能会影响血泵的性能和血液的正常流动。进一步研究不同转速和流量对血泵内部速度场的影响。图2和图3分别展示了转速为2000r/min和3000r/min时，血泵内部的速度分布云图，在相同流量下，随着转速的增加，叶轮区域和蜗壳区域的速度均显著增大。当转速从2000r/min增加到3000r/min时，叶轮边缘处的速度从4m/s左右增加到6m/s左右，蜗壳内靠近叶轮出口处的速度也从2-3m/s增加到4-5m/s。这是因为转速的提高使得叶轮对血液做功增加，血液获得的动能增大，从而速度增大。图4和图5则展示了流量为4L/min和6L/min时的速度分布云图。在相同转速下，随着流量的增加，叶轮区域的速度变化相对较小，但蜗壳区域的速度明显增大。当流量从4L/min增加到6L/min时，蜗壳内的平均速度从1.5-2.5m/s增加到2.5-3.5m/s。这是因为流量的增加意味着单位时间内通过血泵的血液量增多，为了保证血液的连续流动，蜗壳内的流速必然增大。通过对不同工况下血泵内部速度场的分析，可以看出转速和流量对速度分布有着重要影响。合理调整转速和流量，能够优化血泵内部的速度分布，提高血泵的性能。例如，在实际应用中，根据患者的生理需求，选择合适的转速和流量，使血泵内部的速度分布更加均匀，减少流动损失和血液损伤。同时，速度场的分析也为进一步研究血泵内部的压力场、流线分布以及血损情况提供了基础。3.1.2压力场分析压力场是离心式血泵内部流场的重要特性之一，它直接影响着血液的流动和血泵的性能。在稳态工况下，血泵内部的压力分布呈现出一定的规律。图6展示了在设计转速为2500r/min、流量为5L/min的工况下，血泵内部压力分布云图。从图中可以看出，压力在血泵内的变化趋势明显，在叶轮进口处，压力相对较低，约为80mmHg，这是因为血液在此处刚刚进入血泵，尚未获得足够的能量。随着血液在叶轮内的加速流动，在离心力的作用下，压力逐渐升高。在叶轮出口处，压力达到较高值，约为120mmHg，此时血液获得了较大的压力能。进入蜗壳后，压力继续发生变化。在蜗壳靠近叶轮出口的区域，由于高速血液的冲击，压力仍然较高，但随着血液在蜗壳内的流动，蜗壳横截面积逐渐增大，根据伯努利方程，流速减小，压力逐渐升高，在蜗壳出口处，压力进一步升高至约130mmHg。这种压力的变化使得血液能够克服管路阻力，顺利输送到人体循环系统中。压力分布对血液流动有着重要的影响。压力差是推动血液流动的直接动力，血泵内部合理的压力分布能够保证血液的正常循环。在叶轮区域，较大的压力梯度使得血液能够迅速获得能量并加速流动；在蜗壳区域，压力的逐渐升高则有助于血液克服阻力，实现稳定的输送。然而，如果压力分布不均匀，可能会导致血液流动不畅，甚至出现回流现象。例如，在蜗壳内如果存在局部压力过高或过低的区域，可能会使血液在这些区域产生停滞或紊乱的流动，增加血液与泵体壁面的摩擦，从而增加血损的风险。进一步分析不同工况下压力场的变化情况。图7和图8分别展示了转速为2000r/min和3000r/min时的压力分布云图。随着转速的增加，叶轮进口和出口处的压力均显著升高。当转速从2000r/min增加到3000r/min时，叶轮进口压力从70mmHg左右升高到90mmHg左右，叶轮出口压力从100mmHg左右升高到140mmHg左右。这是因为转速的提高使得叶轮对血液做功增加，血液获得的压力能增大。同时，蜗壳内的压力也相应升高，这有助于提高血泵的扬程，增强血液的输送能力。图9和图10展示了流量为4L/min和6L/min时的压力分布云图。随着流量的增加，叶轮进口压力略有降低，而叶轮出口和蜗壳内的压力则有所升高。当流量从4L/min增加到6L/min时，叶轮进口压力从85mmHg左右降低到75mmHg左右，叶轮出口压力从115mmHg左右升高到125mmHg左右，蜗壳出口压力从125mmHg左右升高到135mmHg左右。这是因为流量的增加导致血液在叶轮内的流速加快，进口处的压力损失增大，而出口和蜗壳内由于血液量的增多，压力相应升高。通过对不同工况下压力场的分析可知，转速和流量对压力分布有着显著影响。在实际应用中，需要根据患者的具体需求，合理调整血泵的转速和流量，以获得合适的压力分布，确保血泵能够高效、稳定地工作，同时减少对血液的损伤。3.1.3流线分析流线能够直观地展示血液在血泵内的流动轨迹，通过对流线的研究，可以深入了解血泵内部的流动特性，找出流动不畅区域，为血泵的优化设计提供重要依据。图11展示了在设计转速为2500r/min、流量为5L/min的工况下，血泵内部的流线图。从图中可以清晰地看到，血液从进口进入血泵后，沿着叶轮的旋转方向，在离心力的作用下，呈螺旋状从叶轮中心向叶轮边缘流动。在叶轮区域，流线较为密集且规则，这表明血液在叶轮内的流动较为顺畅，能够有效地获得能量。进入蜗壳后，流线的分布变得相对复杂。在蜗壳靠近叶轮出口的区域，流线仍然保持一定的方向性，但随着血液在蜗壳内的流动，由于蜗壳形状的影响以及叶轮出口射流的作用，部分流线出现了弯曲和交叉的现象。在蜗壳的某些区域，如蜗壳的角落处，流线较为稀疏，甚至出现了局部的回流区域。这些回流区域的存在表明血液在这些区域的流动不畅，容易导致血液的停滞和能量损失。血液在回流区域内的停留时间较长，增加了血液与泵体壁面的接触时间，从而可能增加血损的风险，如血栓形成和溶血等。为了进一步分析不同工况下流线的变化情况，分别对不同转速和流量下的流线进行研究。图12和图13展示了转速为2000r/min和3000r/min时的流线图。随着转速的增加，叶轮内血液的流速加快，流线更加紧密，血液在叶轮内的流动速度明显提高。在蜗壳内，高转速下的流线弯曲程度和交叉现象更加明显，这是由于高速旋转的叶轮产生的射流对蜗壳内流动的影响更大，导致蜗壳内的流动更加复杂。同时，回流区域的范围也有所变化，在较高转速下，回流区域可能会减小，但回流强度可能会增加。图14和图15展示了流量为4L/min和6L/min时的流线图。随着流量的增加，蜗壳内的流线变得更加密集，这是因为单位时间内通过蜗壳的血液量增多，血液在蜗壳内的流速相应增大。同时，流量的变化也会影响蜗壳内回流区域的分布和大小。在大流量工况下，回流区域可能会向蜗壳的更远处扩展，这是由于流量的增加使得蜗壳内的流动更加湍急，更容易产生流动分离和回流现象。通过对不同工况下流线的分析可知，血泵内部存在一些流动不畅的区域，如蜗壳的角落处和部分回流区域。这些区域的存在会影响血泵的性能和血液的正常流动，增加血损的风险。在血泵的设计和优化过程中，需要采取相应的措施，如优化蜗壳的形状、调整叶轮与蜗壳的间隙等，以改善这些区域的流动状况，减少流动不畅区域的范围和影响，提高血泵的血液相容性和工作效率。3.2非稳态工况下流场特性3.2.1心动周期内流场变化在人体生理循环中，心脏的跳动呈现出周期性的特点，这使得离心式血泵在实际工作过程中也处于非稳态工况。为了深入了解血泵在心动周期内的流场特性，对一个完整心动周期（通常设定为0.8s）内血泵内部的流场进行数值模拟分析。图16展示了心动周期内不同时刻血泵内部的速度分布云图。在心动周期开始时（t=0s），血液刚刚进入血泵，叶轮的旋转速度尚未达到稳定状态，此时叶轮中心区域的速度较低，约为0.3m/s，而蜗壳内的速度分布也相对不均匀，靠近进口处的速度略高于其他区域。随着时间的推移，叶轮加速旋转，在t=0.2s时，叶轮边缘处的速度迅速增加到3m/s左右，蜗壳内靠近叶轮出口处的速度也明显增大，达到2m/s左右。在t=0.4s时，叶轮转速达到相对稳定状态，叶轮边缘速度达到最大值，约为4m/s，蜗壳内的速度分布也趋于稳定，但在蜗壳的某些区域仍存在速度较低的回流区。此后，随着心动周期的进行，叶轮转速逐渐降低，血泵内的速度也相应减小。在t=0.6s时，叶轮边缘速度降至3m/s左右，蜗壳内速度也有所下降。到心动周期结束时（t=0.8s），叶轮转速进一步降低，血泵内的速度分布恢复到接近初始状态。压力分布在心动周期内同样呈现出明显的变化规律。图17展示了不同时刻血泵内部的压力分布云图。在心动周期开始时（t=0s），叶轮进口处压力较低，约为75mmHg，叶轮出口处压力约为100mmHg。随着叶轮的加速旋转，压力逐渐升高，在t=0.2s时，叶轮进口压力升高到80mmHg左右，叶轮出口压力升高到110mmHg左右。在t=0.4s时，压力达到相对稳定状态，叶轮进口压力约为85mmHg，叶轮出口压力约为120mmHg，蜗壳出口压力约为130mmHg。随后，随着叶轮转速的降低，压力逐渐减小，在t=0.6s时，叶轮进口压力降至80mmHg左右，叶轮出口压力降至110mmHg左右。到心动周期结束时（t=0.8s），压力恢复到接近初始状态。流线分布在心动周期内也发生了显著变化。图18展示了不同时刻血泵内部的流线图。在心动周期开始时，血液从进口进入血泵后，流线较为紊乱，在叶轮中心区域存在一些小的漩涡。随着叶轮的旋转，流线逐渐变得规则，在t=0.2s时，血液能够较为顺畅地从叶轮中心流向叶轮边缘，蜗壳内的流线也开始呈现出一定的方向性。在t=0.4s时，流线分布最为稳定，血液在叶轮和蜗壳内的流动较为有序，但在蜗壳的角落处仍存在一些局部的回流区域。随着心动周期的进行，叶轮转速降低，流线的规则性逐渐减弱，在t=0.6s时，叶轮中心区域的漩涡再次出现，蜗壳内的回流区域范围有所扩大。到心动周期结束时，流线分布又回到了较为紊乱的状态。通过对心动周期内血泵内部流场参数（速度、压力、流线）随时间的动态变化过程的分析可知，血泵在非稳态工况下的流场特性与稳态工况存在明显差异。这种差异可能会对血泵的性能和血液的正常流动产生重要影响，增加血损的风险。因此，在血泵的设计和优化过程中，需要充分考虑心动周期的影响，以提高血泵在非稳态工况下的性能和血液相容性。3.2.2流量波动影响在实际应用中，离心式血泵的流量会受到多种因素的影响而产生波动，如人体生理状态的变化、血泵控制系统的精度等。为了研究流量波动对血泵性能和流场稳定性的影响，设定血泵的流量在一定范围内波动，模拟不同波动幅度和频率下血泵内部的流场情况。当流量波动时，血泵的扬程会发生相应的变化。图19展示了流量在4-6L/min之间波动时，血泵扬程随时间的变化曲线。可以看出，随着流量的增加，扬程呈现出先升高后降低的趋势。在流量为5L/min左右时，扬程达到最大值，约为125mmHg。当流量偏离这一值时，扬程逐渐下降。这是因为流量的变化会影响血泵内部的流动状态，当流量增加时，叶轮对血液的做功增加，扬程相应升高，但当流量超过一定值后，血泵内部的流动阻力增大，导致扬程下降。流量波动还会导致扬程的不稳定，波动幅度越大，扬程的变化也越剧烈，这可能会影响血泵对血液的稳定输送。叶轮受力也会受到流量波动的显著影响。图20展示了流量波动时叶轮所受轴向力和径向力的变化情况。随着流量的增加，叶轮所受的轴向力和径向力均逐渐增大。当流量从4L/min增加到6L/min时，轴向力从5N左右增加到8N左右，径向力从10N左右增加到15N左右。这是因为流量的增加使得血液对叶轮的冲击力增大，从而导致叶轮受力增大。流量波动还会使叶轮受力产生周期性变化，这种周期性的受力变化可能会对叶轮的稳定性和寿命产生不利影响，增加叶轮疲劳损坏的风险。流量波动对血泵内部流场的稳定性也有重要影响。图21展示了流量波动时血泵内部的流线分布情况。当流量稳定时，血泵内部的流线分布较为规则，血液流动相对平稳。但当流量发生波动时，流线分布变得紊乱，出现了更多的漩涡和流动分离现象。在流量波动较大时，甚至会出现局部的回流区域，这表明流场的稳定性受到了严重破坏。流场的不稳定会导致血液在血泵内的流动不均匀，增加血液与泵体壁面的摩擦，从而增加血损的风险，如血栓形成和溶血等。通过研究流量波动时血泵扬程、叶轮受力及流场稳定性的变化情况可知，流量波动会对血泵的性能和血液相容性产生不利影响。在血泵的实际应用中，需要采取有效的措施来减小流量波动，如优化血泵的控制系统、改进血泵的结构设计等，以保证血泵能够稳定、高效地工作，降低血损风险。四、离心式血泵血损研究4.1血损预测模型4.1.1溶血模型溶血是离心式血泵应用中需要重点关注的血损现象之一，准确预测溶血风险对于评估血泵性能至关重要。溶血的发生主要是由于红细胞在血泵内受到过高的剪切应力作用，导致细胞膜破裂，血红蛋白释放。基于此原理，溶血模型通常通过计算血液在血泵内所受的剪切应力以及红细胞在高剪切应力区域的暴露时间来预测溶血风险。常用的溶血预测模型如Lutz模型，其表达式为：H=k\int_{t_0}^{t_1}\left(\frac{\tau}{\tau_0}\right)^ndt其中，H为溶血指数，表示红细胞的破损程度；k为与红细胞特性相关的常数；\tau为剪切应力；\tau_0为参考剪切应力；n为与剪切应力作用机制相关的指数；t_0和t_1分别为红细胞进入和离开高剪切应力区域的时间。从该模型可以看出，溶血指数与剪切应力的大小以及作用时间密切相关。当剪切应力\tau增大时，\left(\frac{\tau}{\tau_0}\right)^n的值会迅速增大，从而导致积分结果H增大，即溶血风险增加；同时，作用时间t_1-t_0越长，积分值H也会越大，溶血风险相应提高。例如，在离心式血泵的叶轮叶片附近，由于流速梯度较大，血液受到的剪切应力较高，如果红细胞在该区域的停留时间较长，根据Lutz模型计算得到的溶血指数就会较大，表明溶血风险较高。在实际应用中，通过数值模拟获取血泵内部的流场信息，进而计算出各个位置的剪切应力分布和红细胞的运动轨迹，从而确定红细胞在不同剪切应力区域的暴露时间。将这些数据代入溶血模型中，即可预测血泵在不同工况下的溶血风险。例如，对于设计转速为2500r/min、流量为5L/min的离心式血泵，利用数值模拟得到的流场数据，计算出叶轮叶片表面某点的平均剪切应力为150Pa，红细胞在该点附近的暴露时间为0.05s，假设k=1\times10^{-6}，\tau_0=100Pa，n=3，代入Lutz模型可得：H=1\times10^{-6}\int_{0}^{0.05}\left(\frac{150}{100}\right)^3dt=1\times10^{-6}\times\left(\frac{150}{100}\right)^3\times0.05\approx1.69\times10^{-6}通过对血泵内部多个位置进行类似的计算，并综合考虑整个流场的情况，可以得到血泵在该工况下的总体溶血风险评估。这样的预测结果有助于在血泵设计阶段及时发现潜在的溶血问题，通过优化血泵结构（如改进叶轮叶片形状、调整叶轮与蜗壳的间隙等）或调整运行参数（如转速、流量等），降低溶血风险，提高血泵的血液相容性。4.1.2血栓模型血栓形成是离心式血泵另一个重要的血损问题，其机制涉及多个因素的相互作用。目前广泛认可的血栓形成机制基于Virchow提出的三要素：血管内皮损伤、血流状态改变以及血液凝固性增加。在离心式血泵中，血液与泵体表面的接触以及复杂的流场条件使得这三个要素均可能发生变化，从而增加血栓形成的风险。基于此，血栓模型通常综合考虑这些因素来预测血栓形成的可能性。例如，一些血栓模型通过计算壁面剪切应力（WSS）、相对停留时间（RRT）和血小板活化因子（PAF）等参数来评估血栓形成风险。壁面剪切应力反映了血液对泵体壁面的作用力，当壁面剪切应力过高或过低时，都可能导致血小板的活化和聚集。研究表明，当壁面剪切应力低于0.5Pa或高于40Pa时，血栓形成的风险会显著增加。相对停留时间表示血液在特定区域的停留时间与平均停留时间的比值，较长的相对停留时间意味着血液在该区域的流动缓慢，容易引发血小板的聚集和凝血因子的激活。血小板活化因子则是血小板活化过程中释放的一系列生物活性物质，其含量的增加会促进血栓的形成。在实际应用中，通过数值模拟获取血泵内部的流场信息，进而计算出壁面剪切应力分布、血液在各区域的停留时间以及血小板活化因子的浓度分布等参数。将这些参数代入血栓模型中，即可预测血泵在不同工况下的血栓形成风险。例如，对于一个特定的离心式血泵，在设计工况下，通过数值模拟得到蜗壳某区域的壁面剪切应力为0.3Pa，相对停留时间为2.5，血小板活化因子浓度为某一较高值。根据血栓模型，该区域的血栓形成风险较高。通过对血泵内部各个区域进行类似的分析，可以全面评估血泵在不同工况下的血栓形成风险。这样的预测结果对于血泵的优化设计具有重要指导意义，例如，可以通过优化蜗壳形状，改善该区域的流场分布，提高壁面剪切应力，减少血液的停留时间，从而降低血栓形成的风险。同时，在血泵的临床应用中，也可以根据血栓模型的预测结果，采取相应的预防措施，如合理调整血泵的运行参数、使用抗凝血药物等，以降低患者发生血栓栓塞并发症的风险。4.2血损计算结果与分析基于前文建立的溶血模型和血栓模型，对不同工况下离心式血泵的血损情况进行详细计算与深入分析，旨在揭示转速、流量等因素对血损的具体影响规律，为血泵的优化设计和临床应用提供关键的理论依据。首先，针对溶血情况展开研究。表1呈现了不同转速和流量工况下血泵的溶血值计算结果。在流量固定为5L/min时，随着转速从2000r/min提升至3000r/min，溶血值从1.2\times10^{-6}显著增大到2.5\times10^{-6}。这是因为转速的提高使得叶轮对血液的作用力增强，血液流速加快，叶轮叶片附近以及蜗壳内的流速梯度增大，导致血液受到的剪切应力显著增加。根据溶血模型，剪切应力的增大以及红细胞在高剪切应力区域的暴露时间相对延长，共同作用使得溶血风险大幅上升，溶血值随之增大。在转速固定为2500r/min的情况下，当流量从4L/min增加到6L/min时，溶血值从1.5\times10^{-6}上升至1.8\times10^{-6}。流量的增加意味着单位时间内通过血泵的血液量增多，蜗壳内的流速增大，血液与泵体壁面的摩擦加剧，局部区域的剪切应力升高，从而导致溶血值有所上升。不过，相较于转速的影响，流量变化对溶血值的影响相对较小，这表明在该血泵的运行工况中，转速是影响溶血的更为关键的因素。转速（r/min）流量（L/min）溶血值（\times10^{-6}）200051.2250051.8300052.5250041.5250061.8接着，分析血栓形成的情况。表2展示了不同工况下血泵的血栓值计算结果。当流量为5L/min时，随着转速从2000r/min提高到3000r/min，血栓值从0.35增大到0.55。转速的增加使得血泵内部的流场更加复杂，壁面剪切应力分布发生变化，部分区域的壁面剪切应力过高或过低，同时血液在泵内的相对停留时间也可能改变，这些因素综合作用，导致血小板更容易活化和聚集，从而增加了血栓形成的风险，血栓值相应增大。在转速为2500r/min时，流量从4L/min增大到6L/min，血栓值从0.40增大到0.45。流量的增加导致蜗壳内的流动状态改变，可能出现更多的流动分离和漩涡区域，血液在这些区域的相对停留时间延长，血小板与泵体壁面的接触机会增多，进而促进了血小板的活化和聚集，使得血栓值上升。与溶血情况类似，转速对血栓形成的影响相对更为显著。转速（r/min）流量（L/min）血栓值200050.35250050.45300050.55250040.40250060.45综合以上分析可知，转速和流量对离心式血泵的溶血和血栓形成均有显著影响，且转速的影响更为突出。在血泵的实际应用中，为了降低血损风险，需要根据患者的具体需求，在保证血泵提供足够流量和扬程的前提下，合理控制转速和流量。例如，对于一些对血损较为敏感的患者，可以适当降低转速，同时通过优化血泵结构，提高血泵在低转速下的性能，以满足患者的生理需求，减少血损对患者健康的潜在威胁。4.3降低血损的措施探讨为有效降低离心式血泵在运行过程中的血损，提高其血液相容性和临床应用安全性，可从结构优化与控制策略等多方面入手，综合采取一系列针对性措施。在结构优化方面，叶轮作为血泵的核心部件，对其进行优化设计能够显著改善血泵内部流场，从而降低血损。通过优化叶片形状，使其更符合流体动力学原理，可有效减少血液在叶轮内流动时受到的剪切应力。例如，将叶片设计为S型，相较于传统的直线型或圆弧型叶片，S型叶片能够使血液在叶轮内的流动更加顺畅，减少流动分离和漩涡的产生。这是因为S型叶片的特殊形状能够更好地引导血液的流动方向，使血液在叶轮内的速度分布更加均匀，从而降低了局部区域的高剪切应力，减少了红细胞的损伤和血栓形成的风险。研究表明，采用S型叶片的离心式血泵，其溶血值相较于传统叶片结构可降低约20%-30%。合理调整叶片数量也是优化叶轮结构的重要手段。过多的叶片可能会增加叶轮内的流动阻力，导致血液受到的剪切应力增大；而过少的叶片则可能无法提供足够的驱动力，影响血泵的性能。通过数值模拟和实验研究，确定合适的叶片数量，能够在保证血泵性能的前提下，降低血损。例如，对于某一特定规格的离心式血泵，经过优化后，将叶片数量从6片调整为8片，血泵内部的流场更加均匀，血栓形成的风险降低了约15%-20%。蜗壳的结构对血泵内部流场同样有着重要影响，对蜗壳进行优化可以改善血液在蜗壳内的流动状况，减少血损。优化蜗壳的横截面形状是一种有效的方法，将蜗壳横截面设计为内扩梯形，能够使血液在蜗壳内的流动更加平稳，减少流动分离和回流现象。这是因为内扩梯形的横截面形状能够更好地适应血液的流动特性，使血液在蜗壳内的流速逐渐降低，压力逐渐升高，从而减少了局部区域的压力波动和高剪切应力区域的出现。研究发现，采用内扩梯形横截面蜗壳的血泵，其溶血值和血栓值相较于传统蜗壳结构均有显著降低，分别降低了约15%-20%和10%-15%。合理设计蜗壳的出口管也能够改善血泵的性能，减少血损。通过优化出口管的直径、长度和弯曲角度等参数，使血液能够更加顺畅地从蜗壳流出，避免了出口处的流动阻塞和压力突变，从而降低了血损风险。例如，将出口管的弯曲角度从90°调整为60°，并适当增加出口管的直径，可使血泵内部的流场更加稳定，血栓形成的可能性降低约10%-15%。在控制策略方面，智能控制技术的应用为降低血损提供了新的途径。采用自适应控制算法，能够根据血泵的实时运行状态和患者的生理需求，自动调整血泵的转速和流量，使血泵始终工作在最优工况下，从而降低血损。例如，当患者的身体活动量增加，需要更多的血液供应时，自适应控制算法能够自动提高血泵的转速和流量，以满足患者的需求；而当患者处于休息状态时，算法则会相应降低血泵的转速和流量，减少血液受到的剪切应力，降低血损风险。研究表明，采用自适应控制算法的离心式血泵，在不同工况下的溶血值和血栓值相较于传统控制方式均有明显降低，分别降低了约20%-30%和15%-20%。多参数监测与反馈控制也是一种有效的控制策略。通过在血泵内部和患者体内设置多个传感器，实时监测血泵的运行参数（如转速、流量、压力等）以及患者的生理参数（如心率、血压、血氧饱和度等），并将这些参数反馈给控制系统。控制系统根据反馈信息，及时调整血泵的运行状态，以保证血泵的性能稳定，并将血损控制在最低限度。例如，当监测到血泵内部某区域的壁面剪切应力过高时，控制系统可以通过调整血泵的转速或流量，改变该区域的流场分布，降低壁面剪切应力，从而减少血栓形成的风险。这种多参数监测与反馈控制策略能够实现对血泵运行状态的精准调控，有效降低血损，提高血泵的安全性和可靠性。五、实验验证与分析5.1实验装置与方法为了验证数值模拟结果的准确性，搭建了一套专门的离心式血泵实验平台，该平台主要由血泵装置、驱动系统、流量测量系统、压力测量系统以及数据采集与控制系统等部分组成。血泵装置选用了一款具有代表性的离心式血泵，其结构参数与数值模拟中所使用的模型一致，以确保实验与模拟的一致性和可比性。血泵的泵体采用透明有机玻璃材质制作，这不仅便于直接观察血泵内部的流动情况，还能有效减少光线反射和折射对实验测量的干扰。叶轮则采用生物相容性良好的钛合金材料制成，以模拟实际临床应用中的血泵工作环境。驱动系统采用高精度的直流电机，通过变频器可以精确调节电机的转速，从而实现对血泵不同运行工况的模拟。电机与血泵之间通过联轴器进行连接，确保动力的稳定传输，减少振动和噪声对实验结果的影响。流量测量系统选用电磁流量计，该流量计具有精度高、响应速度快、测量范围宽等优点，能够准确测量血泵的输出流量。电磁流量计安装在血泵的出口管道上，其测量原理基于法拉第电磁感应定律，当导电流体在磁场中运动时，会在管道两侧产生感应电动势，该电动势与流体的流速成正比，通过测量感应电动势的大小，即可计算出流体的流量。压力测量系统采用压力传感器，分别安装在血泵的进口和出口管道上，用于测量血泵进出口的压力。压力传感器采用高精度的硅压阻式传感器，具有灵敏度高、线性度好、稳定性强等特点，能够准确测量管道内的压力变化。压力传感器将压力信号转换为电信号，通过数据采集卡传输到计算机中进行处理和分析。在实验过程中，首先将血泵、驱动系统、流量测量系统和压力测量系统按照实验要求进行安装和调试，确保各系统正常工作。然后，向血泵实验回路中注入模拟血液，模拟血液采用与人体血液具有相似流变学特性的溶液，如甘油和水的混合溶液，并添加适量的示踪粒子，以便于后续的流场测量。启动驱动系统，调节电机转速，使血泵达到预定的运行工况。在血泵稳定运行一段时间后，通过流量测量系统和压力测量系统分别测量血泵的输出流量和进出口压力，并将测量数据实时采集到计算机中进行记录和分析。为了获取血泵内部流场的详细信息，采用粒子图像测速（PIV）技术对血泵内部的速度分布进行测量。PIV技术是一种基于图像识别和处理的非接触式流场测量技术，具有测量精度高、空间分辨率高、能够同时测量流场多个点的速度等优点。在实验中，使用脉冲激光器作为光源，通过光学系统将激光片照射到血泵内部的流场中，示踪粒子在激光片的照射下会产生散射光，利用高速摄像机从垂直于激光片的方向拍摄示踪粒子的散射光图像。通过对不同时刻拍摄的图像进行分析和处理，采用互相关算法计算出示踪粒子的位移，进而得到流场中各点的速度矢量。通过上述实验装置和方法，可以准确测量离心式血泵在不同工况下的流量、压力以及内部流场的速度分布等参数，为后续的实验验证和分析提供可靠的数据支持。5.2实验结果与数值模拟对比将实验测量得到的血泵性能参数与数值模拟结果进行对比，能够有效评估数值模拟的准确性，为数值模拟模型的改进和完善提供关键依据。在相同的工况条件下，对血泵的流量、扬程以及内部流场的速度分布等参数进行实验测量与数值模拟计算。表3展示了不同转速下血泵流量和扬程的实验值与模拟值对比情况。在转速为2000r/min时，实验测得的流量为4.8L/min，扬程为105mmHg；而数值模拟得到的流量为4.9L/min，扬程为108mmHg。可以看出，流量的模拟值与实验值相对误差约为2.08%，扬程的相对误差约为2.86%。当转速提升至2500r/min时，实验流量为5.2L/min，模拟流量为5.3L/min，相对误差约为1.92%；实验扬程为120mmHg，模拟扬程为123mmHg，相对误差约为2.50%。在3000r/min转速下，实验流量为5.6L/min，模拟流量为5.7L/min，相对误差约为1.79%；实验扬程为135mmHg，模拟扬程为138mmHg，相对误差约为2.22%。从这些数据可以看出，在不同转速工况下，流量和扬程的模拟值与实验值之间的相对误差均在3%以内，表明数值模拟在预测血泵流量和扬程方面具有较高的准确性。转速（r/min）流量实验值（L/min）流量模拟值（L/min）流量相对误差（%）扬程实验值（mmHg）扬程模拟值（mmHg）扬程相对误差（%）20004.84.92.081051082.8625005.25.31.921201232.5030005.65.71.791351382.22在血泵内部流场速度分布的对比方面，图22展示了转速为2500r/min时，血泵内部某一截面的速度分布实验测量结果与数值模拟结果。从图中可以直观地看出，实验测量得到的速度分布与数值模拟结果在整体趋势上基本一致。在叶轮区域，实验测量的速度从叶轮中心向边缘逐渐增大，与数值模拟结果相符。在蜗壳区域，速度分布的变化趋势也较为相似，实验和模拟结果都显示靠近叶轮出口处速度较高，随着向蜗壳出口流动，速度逐渐降低。通过对该截面上多个点的速度值进行对比分析，计算得到速度的平均相对误差约为5.5%。尽管存在一定的误差，但考虑到实验测量过程中可能存在的各种误差因素，如测量仪器的精度、实验环境的微小变化等，这一误差范围是可以接受的。综合流量、扬程以及内部流场速度分布的对比结果可知，数值模拟能够较为准确地预测离心式血泵的性能和内部流场特性。模拟结果与实验数据之间的良好一致性，验证了数值模拟模型的可靠性和有效性。这为进一步利用数值模拟技术深入研究离心式血泵的性能优化、血损机制以及不同设计方案的评估提供了坚实的基础。同时，对于数值模拟结果与实验数据之间存在的少量误差，后续可以通过进一步优化数值模拟模型、提高网格质量、更精确地设置边界条件等方式进行改进，以进一步提高数值模拟的精度。5.3结果分析与讨论通过实验与数值模拟结果的对比，可发现两者在流量、扬程以及内部流场速度分布等方面呈现出良好的一致性，但也存在一定程度的差异。数值模拟在预测血泵流量和扬程时展现出较高的准确性，相对误差均控制在3%以内，这表明所建立的数值模拟模型能够较为精准地反映血泵在不同转速工况下的性能表现。在流场速度分布方面，实验测量与数值模拟的整体趋势相符，尤其在叶轮和蜗壳区域，速度变化趋势一致，进一步验证了数值模拟模型的可靠性。这种良好的一致性为利用数值模拟技术深入研究血泵性能提供了有力支持，使得我们能够在实际制造血泵之前，通过数值模拟对不同设计方案进行评估和优化，