离心式心脏泵内部流场：定常与非定常数值模拟的深度剖析与对比研究

离心式心脏泵内部流场：定常与非定常数值模拟的深度剖析与对比研究一、引言1.1研究背景与意义心脏病作为全球范围内的高发病症，严重威胁着人类的生命健康。据世界卫生组织（WHO）统计数据显示，每年因心脏病导致的死亡人数在全球总死亡人数中占比颇高。在传统治疗手段中，药物维持治疗成本高昂，且难以从根本上治愈疾病；心脏移植手术虽能显著改善患者病情，但供体心脏的严重短缺极大地限制了其应用。在这样严峻的医疗形势下，人工心脏技术应运而生，为心脏病患者带来了新的希望。离心式心脏泵作为人工心脏的关键部件，凭借其独特的结构和工作原理，在医疗领域中展现出了巨大的应用价值。其结构相对简洁，主要由叶轮、泵壳、进口和出口等部分构成。在工作时，叶轮高速旋转，使血液在离心力的作用下从叶轮中心被甩向四周，从而实现血液的泵送。这种工作方式使得离心式心脏泵具有较高的泵送效率，能够满足人体对血液循环的需求。在实际应用中，离心式心脏泵被广泛应用于体外辅助循环系统，如心室辅助装置（VAD）和体外膜肺氧合（ECMO）机等。在心室辅助装置中，离心式心脏泵可部分或完全替代心脏的泵血功能，帮助心脏功能受损的患者维持血液循环；而在体外膜肺氧合机中，它则负责驱动血液在体外循环中流动，经过氧合后再输回体内，为心肺功能严重衰竭的患者提供有效的生命支持。由于血液的流动特性复杂，且对血细胞的损伤、血栓形成等问题极为敏感，离心式心脏泵内部的流场特性对其性能和临床应用效果有着决定性影响。不合理的流场结构可能导致血液流速分布不均，局部区域出现过高的剪切应力，从而造成血细胞的破坏，引发溶血现象；流场中的流动分离和漩涡也可能促使血栓的形成，一旦血栓脱落进入血液循环，将引发严重的栓塞并发症，危及患者生命。因此，深入研究离心式心脏泵内部的流场特性，对于优化其设计、提高血液相容性和可靠性具有至关重要的意义。数值模拟技术，尤其是计算流体动力学（CFD）方法，为离心式心脏泵的研究提供了强大的工具。通过建立精确的数值模型，CFD能够对心脏泵内部的流场进行详细的模拟和分析，获取传统实验手段难以测量的流场参数，如速度分布、压力分布、涡度等。与实验研究相比，数值模拟具有显著的优势。一方面，它能够大幅降低研究成本，避免了实验所需的昂贵设备和大量耗材；另一方面，数值模拟可以灵活地改变各种参数，如泵的结构参数、运行工况等，快速评估不同设计方案对流场特性的影响，从而为离心式心脏泵的优化设计提供全面、准确的理论依据。综上所述，开展离心式心脏泵内部流场的定常与非定常数值模拟研究，不仅有助于深入理解其内部流动的物理机制，揭示流场特性与泵性能之间的内在联系，还能为离心式心脏泵的设计优化提供关键的技术支持，进而推动人工心脏技术的发展，为广大心脏病患者带来更好的治疗效果和生活质量。1.2国内外研究现状在离心式心脏泵内部流场数值模拟领域，国内外学者开展了大量研究工作，推动着该领域不断发展。国外在该领域起步较早，取得了一系列具有开创性的成果。美国学者在早期利用CFD技术对离心式心脏泵进行模拟研究，通过建立简化的模型，初步分析了泵内的速度和压力分布。例如，[具体学者1]等在[具体年份1]使用商业CFD软件对一款早期的离心式心脏泵进行了定常数值模拟，成功获取了泵内基本的流场参数分布，为后续研究奠定了基础。随着研究的深入，[具体学者2]团队在[具体年份2]针对不同叶轮结构的离心式心脏泵开展了详细的数值模拟，系统地分析了叶轮叶片数、叶片形状等因素对流场特性的影响，发现叶片数的增加在一定程度上能提高泵的扬程，但同时也会增加流动损失，导致局部区域剪切应力增大，对血液相容性产生不利影响。欧洲的研究团队在离心式心脏泵的非定常数值模拟方面成果显著。[具体学者3]等在[具体年份3]运用先进的动网格技术，对离心式心脏泵启动和停止过程中的非定常流场进行了模拟，清晰地捕捉到了瞬态过程中流场的动态变化，揭示了启动和停止瞬间叶轮与流体之间的相互作用机制，以及由此产生的流动不稳定现象，为优化心脏泵的启停控制策略提供了理论依据。德国的研究人员还通过数值模拟与实验相结合的方法，深入研究了离心式心脏泵在不同工况下的性能和流场特性，实验结果与数值模拟结果的良好吻合，验证了数值模拟方法的可靠性，为进一步优化心脏泵的设计提供了有力支持。国内在离心式心脏泵内部流场数值模拟研究方面也取得了长足的进步。近年来，众多高校和科研机构积极投身于该领域的研究。清华大学的研究团队在离心式心脏泵的数值模拟研究中，针对泵内复杂的流动现象，采用高精度的数值计算方法，对泵内的涡结构和湍流特性进行了深入分析。[具体学者4]等在[具体年份4]通过数值模拟发现，离心式心脏泵内部存在多个尺度的涡结构，这些涡结构的产生、发展和相互作用对泵内的能量损失和血液损伤具有重要影响，并基于此提出了通过优化泵壳形状来抑制有害涡结构生成的方法，有效提高了泵的性能和血液相容性。哈尔滨工业大学的学者[具体学者5]等在[具体年份5]针对第三代离心式液力悬浮人工心脏泵，开展了系统的数值模拟研究。通过建立三维流场模型，结合流场分析与溶血数值化分析，深入研究了心脏泵泵体的血流动力学特性以及结构改进对血流动力学性能的影响。研究发现，叶轮入口处增设导流锥可以有效地优化流场，减少流动分离和剪切应力，从而提高人工心脏泵的血液相容性。上海交通大学的研究团队在离心式心脏泵的多物理场耦合数值模拟方面取得了创新性成果。[具体学者6]等在[具体年份6]考虑了血液的非牛顿特性、流固耦合效应以及电磁场对血液流动的影响，建立了多物理场耦合的数值模型，对离心式心脏泵内部的复杂流动进行了全面的模拟和分析。研究结果表明，多物理场耦合作用对泵内流场特性和血液损伤有着不可忽视的影响，为更准确地评估离心式心脏泵的性能和血液相容性提供了新的研究思路。尽管国内外在离心式心脏泵内部流场数值模拟方面已经取得了丰硕的成果，但仍存在一些不足之处。一方面，目前的数值模拟模型在考虑血液的复杂特性，如血液的粘弹性、血细胞的变形性以及血液与泵壁之间的相互作用等方面还不够完善，需要进一步改进和优化数值模型，以更真实地模拟血液在心脏泵内的流动情况；另一方面，对于离心式心脏泵在复杂生理条件下的运行特性，如不同心率、血压以及人体运动状态对泵内流场的影响等方面的研究还相对较少，需要开展更多的研究工作，以全面了解心脏泵在实际应用中的性能表现，为其临床应用提供更坚实的理论基础。1.3研究内容与方法本研究以离心式心脏泵为核心对象，致力于通过先进的数值模拟技术，深入剖析其内部流场的定常与非定常特性，为该领域的发展提供全面而深入的理论依据和实践指导。在研究内容方面，首先，将深入研究离心式心脏泵内部的定常流场特性。运用先进的数值模拟技术，精确获取泵内的速度、压力和涡度分布等关键参数，全面揭示流场的内在特征。通过对不同工况下的定常流场进行模拟分析，深入探究工况变化对这些参数分布的具体影响，从而建立起工况与流场参数之间的内在联系。特别关注叶轮与泵壳之间的间隙流场特性，详细分析间隙大小对流动稳定性、能量损失以及血液损伤的影响机制，为优化泵的结构设计提供关键的数据支持。其次，对离心式心脏泵内部的非定常流场特性展开深入研究。利用高精度的数值计算方法，细致捕捉启动、停止以及变工况过程中流场的动态变化过程，清晰揭示瞬态过程中流场的演变规律和不稳定特性。重点研究叶轮旋转过程中产生的动静干涉现象，深入分析其对泵内压力脉动和流动损失的影响程度，以及由此可能引发的血液损伤风险。同时，探究非定常流场中的涡结构生成、发展和相互作用机制，全面揭示其对泵性能和血液相容性的重要影响，为提升泵的运行稳定性和血液相容性提供理论指导。再者，深入分析离心式心脏泵内部流场特性对泵性能和血液相容性的影响。基于数值模拟获得的流场参数，运用专业的计算方法，准确预测泵的扬程、效率等性能指标，并深入分析流场特性与泵性能之间的内在关联，为优化泵的性能提供科学依据。同时，结合血液损伤模型，系统研究流场中的剪切应力、速度梯度等参数对血细胞的破坏作用，全面评估不同流场条件下的溶血风险和血栓形成倾向，为提高泵的血液相容性提供关键的技术支持。此外，基于数值模拟结果，对离心式心脏泵的结构和运行参数进行优化。通过建立全面的优化模型，综合考虑泵的性能、血液相容性和结构可行性等多方面因素，运用先进的优化算法，对叶轮形状、叶片数、泵壳结构等关键结构参数进行优化设计，同时对转速、流量等运行参数进行优化调整，以实现泵性能和血液相容性的协同提升。通过优化前后的数值模拟对比，全面评估优化效果，验证优化方案的有效性和可行性。在研究方法上，主要采用计算流体动力学（CFD）技术。利用专业的CFD软件，如ANSYSFluent、OpenFOAM等，建立高精度的离心式心脏泵三维流场模型。在模型建立过程中，充分考虑血液的非牛顿特性，采用合适的本构模型来准确描述血液的流变行为，确保模型的准确性和可靠性。对模型进行精细的网格划分，运用先进的网格生成技术，如结构化网格、非结构化网格或混合网格，在保证计算精度的同时，提高计算效率。合理设定边界条件，根据实际运行情况，准确设定进口流速、出口压力等边界条件，确保模拟结果能够真实反映泵内的实际流动情况。同时，结合理论分析方法，深入探讨离心式心脏泵内部流场的基本理论和物理机制。运用流体力学、传热学等相关学科的基本原理，对数值模拟结果进行深入分析和解释，揭示流场特性与泵性能之间的内在联系，为数值模拟提供坚实的理论基础。此外，为了验证数值模拟结果的准确性和可靠性，将开展实验研究。搭建专门的实验平台，采用先进的实验测量技术，如粒子图像测速（PIV）技术、激光多普勒测速（LDV）技术等，对离心式心脏泵内部的流场进行实验测量，获取实际的流场数据。将实验结果与数值模拟结果进行详细对比分析，对数值模型进行验证和修正，进一步提高数值模拟的精度和可信度。二、离心式心脏泵工作原理与结构2.1工作原理离心式心脏泵的工作原理基于离心力的作用，通过叶轮的高速旋转实现血液的高效输送，其工作过程紧密关联着流体力学的基本原理，为维持人体血液循环提供了关键支持。当离心式心脏泵启动时，电机驱动叶轮以一定的角速度绕轴高速旋转。叶轮通常由多个弯曲的叶片组成，这些叶片的形状和布局经过精心设计，以优化血液的流动特性。血液从泵的入口轴向流入叶轮中心区域，此时血液处于相对静止的状态。随着叶轮的快速转动，叶片对血液产生强烈的作用力，使血液在惯性离心力的作用下获得巨大的能量。根据离心力公式F=m\omega^2r（其中F为离心力，m为血液微元的质量，\omega为叶轮的角速度，r为血液微元到叶轮中心的距离），可以清晰地看出，血液在离心力的作用下，会从叶轮中心沿径向迅速被甩向叶轮的外周，速度大幅增加。在血液从叶轮中心向外周运动的过程中，其动能不断增大。当血液到达叶轮外周时，具有较高的流速，一般流速可达数米每秒。此时，血液进入泵壳的蜗形流道。泵壳的设计呈特殊的蜗壳形状，流道的横截面积从叶轮出口到泵的出口逐渐增大。这种渐扩的流道结构具有重要的流体力学意义，它能够使高速流动的血液在流道中逐渐减速。根据伯努利方程p+\frac{1}{2}\rhov^2+\rhogh=C（其中p为压力，\rho为流体密度，v为流速，h为高度，C为常数），在理想情况下，当血液流速降低时，其部分动能会转化为静压能，从而使血液的压力显著升高。这一压力升高的过程至关重要，它确保了血液能够获得足够的压力，以克服人体循环系统中的阻力，顺利地输送到全身各个组织和器官，满足人体正常的生理需求。在整个工作过程中，离心式心脏泵的流量和压力输出与叶轮的转速密切相关。通过调节电机的转速，可以灵活地改变叶轮的旋转速度，进而精确地控制泵的流量和压力。一般来说，随着叶轮转速的增加，离心力增大，血液获得的能量增多，泵的流量和压力也随之增大。这种通过调节转速来控制泵性能的方式，使得离心式心脏泵能够根据患者的不同生理状况和需求，提供个性化的血液泵送支持，具有很强的适应性和灵活性。2.2基本结构离心式心脏泵主要由叶轮、蜗壳、进口管、出口管等部件构成，各部件紧密协作，共同确保心脏泵高效、稳定地工作。叶轮是离心式心脏泵的核心部件，其结构设计直接影响泵的性能和血液流动特性。叶轮通常由轮毂和多个叶片组成，叶片均匀分布在轮毂周围，呈放射状排列。叶片的形状复杂多样，常见的有后弯式、前弯式和径向式等。后弯式叶片在离心式心脏泵中应用较为广泛，其出口角一般小于90°，这种叶片形状能使液体在叶轮出口处获得较为合理的速度和压力分布，有效提高泵的效率，同时减少能量损失和流动损失，降低对血液的损伤风险。叶片的数量也是影响叶轮性能的重要因素，一般来说，叶片数量增多，叶轮对液体的作用力更加均匀，能够提高泵的扬程，但也会增加液体在叶片间的流动阻力，导致能量损失增大，因此需要在设计时进行综合考虑，以确定最佳的叶片数量，通常叶轮的叶片数量在3-8片之间。蜗壳是围绕在叶轮外周的螺旋形流道，其主要作用是收集从叶轮流出的高速液体，并将液体的动能有效地转化为压力能，同时引导液体平稳地流向出口管。蜗壳的截面形状通常为渐扩型，从叶轮出口到蜗壳出口，流道的横截面积逐渐增大。这种渐扩的结构设计符合流体力学原理，当高速流动的液体进入蜗壳后，随着流道面积的增大，液体流速逐渐降低。根据伯努利方程，流速的降低会导致液体压力升高，从而实现了动能到压力能的转化，使液体在蜗壳出口处具有较高的压力，满足人体血液循环的需求。蜗壳的材料一般选用具有良好生物相容性和耐腐蚀性的金属或高分子材料，以确保在长期使用过程中不会对血液产生不良影响。进口管位于离心式心脏泵的入口端，其作用是将血液平稳地引入叶轮中心。进口管的形状和尺寸对血液的进口流速和流动均匀性有着重要影响。为了减少进口处的流动阻力和能量损失，进口管通常设计为光滑的圆形截面，且其内径与叶轮进口直径相匹配，以保证血液能够顺畅地流入叶轮。在进口管的内部，可能会设置一些导流装置，如导流叶片或导流锥等，这些装置能够引导血液以合适的角度和速度进入叶轮，避免血液在进口处发生冲击和漩涡，从而减少对血细胞的损伤。出口管则连接在蜗壳的出口，负责将经过加压的血液输送到人体的血液循环系统中。出口管的内径和长度需要根据实际应用需求进行合理设计，以确保血液能够在一定的压力下顺利输送到目标位置。出口管的材料同样需要具备良好的生物相容性和强度，以承受血液的压力和长期的流动冲击。为了监测和控制血液的输出参数，出口管上可能会安装一些传感器，如压力传感器、流量传感器等，这些传感器能够实时获取血液的压力和流量信息，为心脏泵的运行控制和性能监测提供重要的数据支持。2.3与人体生理需求的适配性离心式心脏泵在工作原理和结构设计上与人体生理需求展现出高度的适配性，这使其能够有效地替代或辅助自然心脏的功能，维持人体正常的血液循环。从工作原理来看，离心式心脏泵的连续泵送特性与人体血液循环的连续性需求相契合。人体的血液循环是一个不间断的过程，心脏需要持续地将血液泵送到全身各个部位，以满足组织和器官对氧气和营养物质的需求，并带走代谢产物。离心式心脏泵通过叶轮的高速旋转，能够持续地将血液从入口吸入并加压后输送到出口，实现了血液的连续流动，确保了人体血液循环的稳定运行。这种连续泵送的方式避免了血液在泵内的停滞，减少了血栓形成的风险，保障了血液在循环系统中的顺畅流动，为人体各器官的正常生理功能提供了可靠的支持。在满足人体不同生理状态下的流量需求方面，离心式心脏泵具有显著的优势。人体在不同的生理状态下，如休息、运动、睡眠等，对血液流量的需求差异较大。例如，在安静休息时，人体的基础代谢率较低，心脏的输出量相对较小，一般成年人的心输出量约为4-6L/min；而在剧烈运动时，肌肉组织对氧气和营养物质的需求急剧增加，心脏需要大幅提高输出量，以满足身体的代谢需求，此时心输出量可增加至15-25L/min甚至更高。离心式心脏泵通过调节叶轮的转速，能够灵活地改变泵的流量输出。通过先进的控制系统，可以根据人体的生理信号，如心率、血压等，实时监测人体的生理状态，并自动调节心脏泵的转速，从而精确地控制泵的流量，使其能够根据人体的实际需求提供相应的血液供应量，满足人体在不同生理状态下的血液循环需求。在压力输出方面，离心式心脏泵同样能够满足人体循环系统的压力要求。人体的血液循环系统存在一定的阻力，心脏需要产生足够的压力来克服这些阻力，推动血液在血管中流动。动脉系统中的血压通常维持在一定的范围内，以保证血液能够顺利地输送到各个器官和组织。离心式心脏泵在工作过程中，通过叶轮对血液施加离心力，使血液获得较高的动能，然后在蜗壳中，血液的动能逐渐转化为压力能，从而使血液在出口处具有足够的压力。根据人体生理需求，离心式心脏泵的出口压力一般能够达到100-150mmHg，这与人体正常动脉血压范围相匹配，能够有效地克服人体循环系统的阻力，确保血液在血管中正常流动，为组织和器官提供充足的血液灌注。从结构设计角度分析，离心式心脏泵的设计充分考虑了血液相容性，以减少对血细胞的损伤。叶轮和泵壳的表面通常经过特殊处理，使其表面光滑，减少血液与壁面之间的摩擦和剪切力，降低血细胞的破坏风险。同时，在叶轮的设计上，通过优化叶片的形状、数量和角度，使血液在叶轮中的流动更加平稳，避免了流动分离和漩涡的产生，进一步减少了对血细胞的损伤。例如，采用流线型的叶片设计，能够引导血液沿着叶片表面平滑流动，减少局部流速的突变和剪切应力的集中，从而降低溶血的可能性。此外，一些离心式心脏泵还在结构中采用了特殊的涂层或材料，以提高血液相容性，减少血栓形成的倾向，确保血液在泵内的流动安全可靠。离心式心脏泵的紧凑结构也使其在实际应用中更具优势。由于其体积小、重量轻，便于植入人体或集成到便携式的医疗设备中，提高了患者的生活质量和活动便利性。在心室辅助装置中，紧凑的离心式心脏泵可以更方便地植入患者体内，减少对患者身体的负担，同时也降低了手术的复杂性和风险。这种紧凑的结构设计使得离心式心脏泵能够更好地适应人体的生理环境，为患者提供长期、有效的心脏辅助支持。三、数值模拟理论基础3.1计算流体力学（CFD）基本理论计算流体力学（ComputationalFluidDynamics，CFD）是一门融合了计算机技术、数值计算方法和流体力学理论的交叉学科，它通过数值计算和图像显示等手段，对包含有流体流动和热传导等相关物理现象的系统进行分析和研究。其核心在于利用计算机强大的计算能力，求解描述流体运动的数学方程，从而获得流体流动的各种参数，如速度、压力、温度等在空间和时间上的分布。CFD的发展历程是一部充满创新与突破的历史，自20世纪中叶诞生以来，它经历了从理论探索到实际应用的巨大跨越。20世纪60年代，随着电子计算机的出现和发展，CFD开始崭露头角。当时，科学家们开始尝试用数值方法求解简单的流体力学方程，尽管计算能力有限，但这些早期的尝试为CFD的发展奠定了基础。例如，在航空航天领域，为了满足飞行器设计对空气动力学性能的精确需求，CFD技术逐渐被应用于模拟飞行器周围的气流，帮助工程师优化飞行器的外形设计。到了70年代，CFD取得了重要进展，有限差分法、有限体积法和有限元法等数值方法逐渐成熟，并被广泛应用于CFD计算中。这些数值方法的出现，使得CFD能够处理更加复杂的流体力学问题。同时，随着计算机性能的不断提升，CFD的计算规模和精度也得到了显著提高。在这一时期，CFD开始在工业领域得到广泛应用，如汽车制造业中用于优化汽车外形以降低风阻，船舶工业中用于模拟船舶的水动力性能等。80年代至90年代，CFD技术迎来了快速发展阶段。计算机技术的飞速进步，特别是并行计算技术的出现，使得CFD能够处理大规模的复杂计算问题。同时，CFD软件的开发也取得了重大突破，商业化的CFD软件如ANSYSFluent、CFX等开始广泛应用于各个领域，这些软件具有友好的用户界面和强大的功能，使得CFD技术更加易于使用和推广。在能源领域，CFD被用于研究燃烧过程，优化燃烧器的设计，提高燃烧效率，减少污染物排放；在环境科学领域，CFD被用于模拟大气和水体的流动，预测污染物的扩散，为环境保护提供科学依据。进入21世纪，CFD技术继续深入发展，在多物理场耦合、高精度算法、复杂几何形状处理等方面取得了一系列重要成果。随着科学研究和工程应用对CFD精度和效率要求的不断提高，研究人员不断开发新的数值算法和模型，以更好地模拟复杂的流体物理现象。例如，在生物医学领域，CFD被用于研究人体血液循环系统和呼吸系统中的流体流动，为疾病的诊断和治疗提供重要的理论支持；在微机电系统（MEMS）领域，CFD被用于研究微尺度下的流体流动特性，为微器件的设计和优化提供技术支持。在流体力学研究中，CFD发挥着不可替代的重要作用，具有传统实验方法难以比拟的优势。首先，CFD能够模拟各种复杂的流动现象，无论是在极端条件下的流动，如高温、高压、高马赫数等，还是在微观尺度下的流动，CFD都能够通过数值计算进行深入研究。例如，在航空发动机的设计中，CFD可以模拟发动机内部高温、高压燃气的流动，帮助工程师优化发动机的燃烧效率和性能，提高发动机的可靠性和耐久性；在微流控芯片的设计中，CFD可以模拟微通道内液体的流动，为芯片的结构优化提供依据，提高芯片的分析精度和效率。其次，CFD具有成本低、周期短的显著优势。与传统的实验研究相比，CFD不需要建造昂贵的实验设备，也不需要消耗大量的实验材料和时间。通过在计算机上进行数值模拟，研究人员可以快速获得大量的流场数据，对不同的设计方案进行评估和优化。这使得CFD在产品研发过程中能够大大缩短研发周期，降低研发成本。例如，在汽车研发过程中，利用CFD技术可以在设计阶段对汽车的空气动力学性能进行模拟分析，快速调整设计方案，减少风洞实验的次数，从而节省大量的研发成本和时间。此外，CFD还能够提供详细的流场信息，这些信息在实验中往往难以测量得到。通过CFD模拟，研究人员可以获得流场中任意位置的速度、压力、温度等参数的分布，以及流场中的涡结构、边界层特性等信息，从而深入了解流体流动的物理机制。在离心泵内部流场的研究中，CFD可以精确地模拟泵内的速度、压力和涡度分布，揭示叶轮与泵壳之间的间隙流场特性，为离心泵的优化设计提供全面而准确的数据支持。CFD也存在一定的局限性。由于数值模拟是基于数学模型和假设条件进行的，模拟结果可能与实际情况存在一定的偏差。在建立CFD模型时，需要对实际问题进行简化和假设，这些简化和假设可能会忽略一些次要因素，但在某些情况下，这些次要因素可能会对模拟结果产生重要影响。此外，CFD计算结果的准确性还受到数值方法的精度、网格质量、边界条件的设定等多种因素的影响。如果数值方法选择不当，网格划分不合理，或者边界条件设定不准确，都可能导致模拟结果的误差增大。因此，在使用CFD技术时，需要对模拟结果进行严格的验证和分析，结合实验数据进行对比和修正，以提高模拟结果的可靠性。3.2控制方程在对离心式心脏泵内部流场进行数值模拟时，控制方程是描述流体运动的核心数学模型，它基于质量守恒、动量守恒和能量守恒等基本物理定律建立，为准确模拟流场特性提供了坚实的理论基础。这些控制方程在流体力学领域中具有广泛的应用，是研究各种流体流动问题的关键工具。通过对控制方程的求解，可以获得流场中速度、压力、温度等重要参数的分布情况，从而深入了解流体的运动规律和物理特性。在离心式心脏泵的研究中，控制方程能够帮助我们揭示泵内血液的流动特性，分析不同结构和工况对泵性能的影响，为泵的优化设计提供理论依据。下面将详细介绍连续性方程、动量方程和能量方程的基本形式及其在离心式心脏泵流场模拟中的应用。3.2.1连续性方程连续性方程是描述流体质量守恒的基本方程，其物理意义在于，在一个封闭的流体系统中，无论流体如何流动，系统内的总质量始终保持不变。这一原理在流体力学中具有基础性的地位，它反映了物质的基本属性，即质量既不会凭空产生，也不会凭空消失。在直角坐标系下，连续性方程的微分形式可表示为：\frac{\partial\rho}{\partialt}+\frac{\partial(\rhou)}{\partialx}+\frac{\partial(\rhov)}{\partialy}+\frac{\partial(\rhow)}{\partialz}=0其中，\rho为流体密度，t为时间，u、v、w分别为流体在x、y、z方向上的速度分量。方程左边第一项\frac{\partial\rho}{\partialt}表示单位时间内流体密度的变化率，它反映了流体密度随时间的动态变化情况；后三项\frac{\partial(\rhou)}{\partialx}+\frac{\partial(\rhov)}{\partialy}+\frac{\partial(\rhow)}{\partialz}则表示通过控制体表面的质量通量，即单位时间内通过单位面积的质量流量，它们描述了流体在空间各个方向上的流动对质量分布的影响。在离心式心脏泵内部流场模拟中，连续性方程起着至关重要的作用。它确保了在数值模拟过程中，进入泵内的血液质量与流出泵的血液质量相等，从而保证了模拟结果的物理合理性。当血液在离心式心脏泵内流动时，由于叶轮的旋转和泵壳的约束，血液的速度和密度在空间上会发生复杂的变化。通过连续性方程，可以准确地描述这些变化，为后续的动量方程和能量方程的求解提供必要的基础。在叶轮区域，血液受到离心力的作用，速度和密度会发生显著变化，连续性方程能够精确地捕捉这些变化，确保质量守恒在整个流场中得到满足。如果连续性方程不成立，就会导致模拟结果出现质量不守恒的错误，如流入泵内的血液质量与流出泵的血液质量不一致，这将严重影响对泵内流场特性的分析和理解。因此，在离心式心脏泵内部流场的数值模拟中，严格满足连续性方程是保证模拟结果准确性和可靠性的关键。3.2.2动量方程动量方程是基于牛顿第二定律推导而来，它深刻地体现了流体的受力和运动变化之间的紧密关系。牛顿第二定律指出，物体的加速度与作用在它上面的合外力成正比，与物体的质量成反比。在流体力学中，动量方程将这一定律应用于流体微元，描述了单位体积流体的动量变化率与作用在该体积上的外力之间的平衡关系。在直角坐标系下，不可压缩粘性流体的动量方程（即Navier-Stokes方程，简称N-S方程）的一般形式为：\rho\left(\frac{\partialu}{\partialt}+u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy}+w\frac{\partialu}{\partialz}\right)=-\frac{\partialp}{\partialx}+\mu\left(\frac{\partial^{2}u}{\partialx^{2}}+\frac{\partial^{2}u}{\partialy^{2}}+\frac{\partial^{2}u}{\partialz^{2}}\right)+\rhof_x\rho\left(\frac{\partialv}{\partialt}+u\frac{\partialv}{\partialx}+v\frac{\partialv}{\partialy}+w\frac{\partialv}{\partialz}\right)=-\frac{\partialp}{\partialy}+\mu\left(\frac{\partial^{2}v}{\partialx^{2}}+\frac{\partial^{2}v}{\partialy^{2}}+\frac{\partial^{2}v}{\partialz^{2}}\right)+\rhof_y\rho\left(\frac{\partialw}{\partialt}+u\frac{\partialw}{\partialx}+v\frac{\partialw}{\partialy}+w\frac{\partialw}{\partialz}\right)=-\frac{\partialp}{\partialz}+\mu\left(\frac{\partial^{2}w}{\partialx^{2}}+\frac{\partial^{2}w}{\partialy^{2}}+\frac{\partial^{2}w}{\partialz^{2}}\right)+\rhof_z其中，p为流体压力，\mu为动力粘度，f_x、f_y、f_z分别为单位质量流体在x、y、z方向上受到的体积力，如重力、电磁力等。方程左边表示单位体积流体的动量变化率，它由两部分组成：\frac{\partialu}{\partialt}、\frac{\partialv}{\partialt}、\frac{\partialw}{\partialt}表示当地加速度，即由于时间变化导致的速度变化率；u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy}+w\frac{\partialu}{\partialz}、u\frac{\partialv}{\partialx}+v\frac{\partialv}{\partialy}+w\frac{\partialv}{\partialz}、u\frac{\partialw}{\partialx}+v\frac{\partialw}{\partialy}+w\frac{\partialw}{\partialz}表示迁移加速度，即由于流体微元在空间位置的变化而引起的速度变化率。方程右边第一项-\frac{\partialp}{\partialx}、-\frac{\partialp}{\partialy}、-\frac{\partialp}{\partialz}表示压力梯度力，它是由于压力分布不均匀而产生的作用力，驱使流体从高压区域流向低压区域；第二项\mu\left(\frac{\partial^{2}u}{\partialx^{2}}+\frac{\partial^{2}u}{\partialy^{2}}+\frac{\partial^{2}u}{\partialz^{2}}\right)、\mu\left(\frac{\partial^{2}v}{\partialx^{2}}+\frac{\partial^{2}v}{\partialy^{2}}+\frac{\partial^{2}v}{\partialz^{2}}\right)、\mu\left(\frac{\partial^{2}w}{\partialx^{2}}+\frac{\partial^{2}w}{\partialy^{2}}+\frac{\partial^{2}w}{\partialz^{2}}\right)表示粘性力，它是由于流体内部各层之间的相对运动而产生的摩擦力，阻碍流体的流动；第三项\rhof_x、\rhof_y、\rhof_z表示体积力，它是作用在整个流体体积上的外力。在离心式心脏泵的工作过程中，动量方程全面地描述了血液在叶轮和泵壳内的受力情况以及由此导致的运动变化。当叶轮高速旋转时，叶片对血液施加切向力，使血液获得切向速度，同时血液在离心力的作用下向叶轮外周运动，速度和方向发生改变。在这个过程中，动量方程中的各项力都发挥着重要作用。压力梯度力推动血液在泵内流动，从叶轮中心流向出口；粘性力则会消耗血液的能量，导致流动损失；如果考虑重力等体积力，它们也会对血液的运动产生一定的影响。通过求解动量方程，可以精确地计算出泵内各点的速度和压力分布，进而深入分析泵的性能和效率。在叶轮出口附近，通过动量方程可以计算出由于叶轮旋转产生的离心力对血液速度和压力的影响，以及粘性力对流动损失的影响，从而为优化叶轮的设计提供理论依据。3.2.3能量方程能量方程是基于能量守恒定律建立的，它在分析流体能量转化过程中具有关键作用。能量守恒定律表明，在一个封闭系统中，能量不会凭空产生或消失，只会从一种形式转化为另一种形式，或者从一个物体转移到另一个物体。在流体力学中，能量方程描述了单位体积流体的能量变化率与外界对流体做功以及流体内部能量转化之间的平衡关系。对于不可压缩流体，忽略重力和热辐射等因素，能量方程的一般形式为：\rhoc_p\left(\frac{\partialT}{\partialt}+u\frac{\partialT}{\partialx}+v\frac{\partialT}{\partialy}+w\frac{\partialT}{\partialz}\right)=k\left(\frac{\partial^{2}T}{\partialx^{2}}+\frac{\partial^{2}T}{\partialy^{2}}+\frac{\partial^{2}T}{\partialz^{2}}\right)+\Phi其中，c_p为流体的定压比热容，T为流体温度，k为热导率，\Phi为粘性耗散项，表示由于流体粘性作用导致的机械能转化为热能的速率。方程左边表示单位体积流体的内能变化率，它由两部分组成：\frac{\partialT}{\partialt}表示当地内能变化率，即由于时间变化导致的内能变化；u\frac{\partialT}{\partialx}+v\frac{\partialT}{\partialy}+w\frac{\partialT}{\partialz}表示迁移内能变化率，即由于流体微元在空间位置的变化而引起的内能变化。方程右边第一项k\left(\frac{\partial^{2}T}{\partialx^{2}}+\frac{\partial^{2}T}{\partialy^{2}}+\frac{\partial^{2}T}{\partialz^{2}}\right)表示热传导项，它描述了由于温度梯度引起的热量传递；第二项\Phi表示粘性耗散项，在流体流动过程中，粘性力会使流体内部的机械能转化为热能，导致流体温度升高，\Phi就是用来衡量这种能量转化的速率。在离心式心脏泵内部流场中，虽然血液温度的变化通常较小，但能量方程对于深入理解泵内的能量转化过程仍然具有重要意义。在血液流动过程中，由于叶轮对血液做功，血液获得机械能，同时粘性力会消耗部分机械能，将其转化为热能。通过能量方程，可以准确地计算出这些能量的转化情况，评估泵的能量利用效率。在叶轮高速旋转时，叶轮对血液做功，使血液的动能增加，同时粘性力在血液与叶轮、泵壳之间产生摩擦，导致部分机械能转化为热能，能量方程能够精确地描述这一能量转化过程。了解泵内的能量转化情况有助于优化泵的设计，减少能量损失，提高泵的效率。通过分析能量方程中的各项参数，可以找出能量损失较大的区域和原因，从而采取相应的措施进行改进，如优化叶轮的形状和表面粗糙度，以减少粘性耗散，提高泵的能量利用效率。3.3湍流模型在离心式心脏泵内部流场的数值模拟中，湍流模型的选择至关重要，它直接影响到模拟结果的准确性和可靠性。湍流是一种高度复杂的三维非稳态、带旋转的不规则流动，其内部存在着各种大小不同的涡结构，这些涡结构的不断产生、发展和相互作用，使得湍流现象的描述和模拟极具挑战性。目前，常用的湍流模型主要包括标准k-ε模型、RNGk-ε模型和Realizablek-ε模型等，它们各自具有独特的特点和适用范围。下面将对这三种湍流模型进行详细介绍和分析。3.3.1标准k-ε模型标准k-ε模型是两方程湍流模型中最具代表性且在工程中应用最为普遍的模式。它通过求解两个单独的输运方程来确定湍流长度和时间尺度，从而描述湍流的特性。其中，k代表紊流脉动动能，k值越大，表明紊流脉动长度和时间尺度越大，反映了湍流中较大尺度涡旋的运动强度；ε表示紊流脉动动量的耗散率，ε值越大，则代表紊流脉动长度和时间尺度越小，体现了小尺度涡旋的能量耗散情况。标准k-ε模型假设湍流粘性和湍动能及耗散率有关，其湍动能k和湍流耗散率ε的输运方程如下：湍动能k方程：\frac{\partial(\rhok)}{\partialt}+\frac{\partial(\rhoku_{i})}{\partialx_{i}}=\frac{\partial}{\partialx_{j}}\left[(\mu+\frac{\mu_{t}}{\sigma_{k}})\frac{\partialk}{\partialx_{j}}\right]+G_{k}+G_{b}-\rho\varepsilon-Y_{M}湍流耗散率ε方程：\frac{\partial(\rho\varepsilon)}{\partialt}+\frac{\partial(\rho\varepsilonu_{i})}{\partialx_{i}}=\frac{\partial}{\partialx_{j}}\left[(\mu+\frac{\mu_{t}}{\sigma_{\varepsilon}})\frac{\partial\varepsilon}{\partialx_{j}}\right]+C_{1\varepsilon}\frac{\varepsilon}{k}(G_{k}+C_{3\varepsilon}G_{b})-C_{2\varepsilon}\rho\frac{\varepsilon^{2}}{k}其中，\rho为流体密度，t为时间，u_{i}为速度分量，\mu为分子粘性系数，\mu_{t}为湍流粘性系数，\sigma_{k}和\sigma_{\varepsilon}分别为k和ε的湍流普朗特数，G_{k}为由于平均速度梯度引起的湍动能生成项，G_{b}为由于浮力引起的湍动能生成项，Y_{M}为可压缩湍流中脉动扩张的贡献，C_{1\varepsilon}、C_{2\varepsilon}和C_{3\varepsilon}为经验常数。该模型的优点在于计算相对简单，对完全湍流的情况具有较好的模拟能力，能够处理较为复杂的紊流，在许多工程领域，如航空航天、汽车工程、能源动力等，都得到了广泛的应用。在航空发动机内部流场的模拟中，标准k-ε模型能够有效地预测气流的流动特性，为发动机的设计和优化提供重要依据；在汽车空气动力学研究中，它可以模拟汽车周围的气流分布，帮助工程师降低风阻，提高汽车的燃油经济性。在离心式心脏泵内部流场模拟中，标准k-ε模型也有一定的应用。它可以对泵内的平均流动特性进行较好的预测，如速度和压力的分布情况。由于该模型在处理近壁面和雷诺数较低的过渡层时存在一定的局限性，对于离心式心脏泵内部流场中一些复杂的流动现象，如叶轮与泵壳间隙处的流动细节、边界层内的流动特性等，模拟结果可能不够准确。在近壁面区域，标准k-ε模型对壁面函数的依赖较强，而壁面函数的准确性在一定程度上影响了模拟结果的精度；在低雷诺数区域，该模型的假设条件与实际情况存在偏差，导致模拟结果与实际流动存在差异。3.3.2RNGk-ε模型RNGk-ε模型是对标准k-ε模型的一种改进，它在标准k-ε模型的基础上，通过重整化群理论（RNG）对湍流进行分析，引入了一些修正项，从而在某些方面提高了模型的性能。与标准k-ε模型相比，RNGk-ε模型在湍动能k方程和湍流耗散率ε方程中增加了一些附加项，以更好地描述湍流的特性。在ε方程中，RNGk-ε模型考虑了平均应变率的影响，对湍流耗散率的计算进行了修正，使其能够更准确地反映湍流的能量耗散过程。此外，RNGk-ε模型还对湍流粘性系数的计算进行了改进，考虑了湍流漩涡的拉伸和旋转等因素，使得模型在处理复杂流动时具有更好的适应性。RNGk-ε模型的优势在于对强旋流、浮力影响以及低雷诺数影响等方面的预测更加准确。在模拟具有强旋流的流动时，如旋风分离器内部的气流运动，RNGk-ε模型能够更精确地捕捉到旋流的特性，包括旋流的强度、方向和分布等，从而为旋风分离器的设计和优化提供更可靠的依据；在处理浮力影响较大的流动，如大气边界层中的热对流现象时，该模型能够充分考虑浮力对湍流的作用，准确地模拟出温度和速度的分布情况。在离心式心脏泵内部流场模拟中，RNGk-ε模型相较于标准k-ε模型具有一定的优势。由于离心式心脏泵内部存在叶轮的高速旋转，会产生较强的旋流，同时血液的流动也会受到一定的浮力影响。RNGk-ε模型能够更好地模拟这些复杂的流动情况，更准确地捕捉到叶轮附近的涡结构和流场的细节信息，从而为分析离心式心脏泵的性能和血液损伤机制提供更全面、准确的数据支持。它可以更精确地预测泵内的压力脉动和流动损失，有助于优化泵的结构设计，提高泵的效率和血液相容性。3.3.3Realizablek-ε模型Realizablek-ε模型是一种基于物理机制的湍流模型，它在湍动能k方程的基础上，对湍流耗散率ε方程进行了重新推导和修正，使其具有更明确的物理意义和更好的通用性。Realizablek-ε模型的一个重要特点是对湍流粘性系数的计算进行了改进，它考虑了湍流的各向异性和旋转效应，使得模型在处理复杂流动时能够更准确地描述湍流的特性。在ε方程中，Realizablek-ε模型引入了与平均应变率和旋转率相关的项，从而能够更好地反映湍流的产生和耗散过程。此外，该模型还对近壁面区域的处理进行了优化，采用了更合适的壁面函数，提高了在近壁面区域的模拟精度。该模型的适用范围广泛，尤其在处理具有复杂几何形状和边界条件的流动问题时表现出色。在叶轮机械内部流场的模拟中，Realizablek-ε模型能够准确地预测叶轮与机壳之间的间隙流动、叶片表面的边界层特性以及尾迹区的流动情况，为叶轮机械的设计和性能优化提供了有力的工具；在汽车发动机的燃烧室内流场模拟中，它可以精确地模拟燃料与空气的混合过程、燃烧过程以及燃烧产物的排放情况，有助于提高发动机的燃烧效率和降低污染物排放。在离心式心脏泵内部流场模拟中，Realizablek-ε模型能够很好地适应泵内复杂的几何结构和流动条件。它可以更准确地模拟血液在叶轮和泵壳内的流动特性，包括速度分布、压力分布和涡度分布等，为研究离心式心脏泵的性能和血液相容性提供了更可靠的依据。通过该模型的模拟，可以更清晰地了解泵内的流动细节，如叶轮入口处的流动分离现象、泵壳内的二次流等，从而为优化泵的结构设计，减少血液损伤和血栓形成提供有效的指导。3.4方程离散化方法在离心式心脏泵内部流场的数值模拟中，方程离散化方法是将描述流体运动的连续控制方程转化为离散的代数方程，以便于计算机进行求解的关键技术。不同的离散化方法具有各自的特点和适用范围，对模拟结果的准确性、计算效率和稳定性都有着重要影响。下面将详细介绍有限体积法和有限差分法这两种常用的方程离散化方法。3.4.1有限体积法有限体积法（FiniteVolumeMethod，FVM）是一种基于守恒原理的数值计算方法，在数值模拟领域应用广泛，尤其在计算流体力学中发挥着重要作用。其基本原理是将计算区域划分为一系列互不重叠的控制体积，然后对每个控制体积内的守恒型控制方程进行积分。以二维连续性方程为例，在笛卡尔坐标系下，其积分形式为：\frac{\partial}{\partialt}\int_{V}\rhodV+\oint_{S}\rho\vec{v}\cdotd\vec{S}=0其中，V为控制体积，S为控制体积的表面，\rho为流体密度，\vec{v}为速度矢量。在有限体积法中，通过将控制体积内的积分转化为对控制体积表面的通量积分，从而得到离散的代数方程。具体来说，将控制体积表面划分为多个面元，对每个面元上的通量进行近似计算，然后将这些面元的通量求和，得到控制体积表面的总通量。通过这种方式，将连续的控制方程离散为一组关于控制体积中心节点上未知变量的代数方程。有限体积法具有诸多优点，使其在离心式心脏泵内部流场模拟中备受青睐。该方法严格满足守恒定律，在离散过程中确保了物理量在每个控制体积内的守恒，这对于模拟离心式心脏泵内部血液的流动至关重要，能够保证模拟结果在物理上的合理性。在处理复杂的几何形状时，有限体积法表现出良好的适应性。它可以根据几何形状的特点，灵活地划分控制体积，无论是规则的还是不规则的几何区域，都能有效地进行离散处理，从而准确地模拟离心式心脏泵内部复杂的流道结构中的流动。有限体积法还能够方便地处理各种边界条件。在离心式心脏泵的模拟中，进口、出口以及壁面等边界条件的准确设定对模拟结果的准确性有着重要影响。有限体积法可以通过在边界上的控制体积表面设置合适的通量条件，轻松地实现对不同边界条件的处理，确保模拟结果能够真实反映实际流动情况。在进口边界，可以根据实际的流量或速度条件，设定进口面元上的通量；在壁面边界，可以采用无滑移边界条件，通过设置壁面处的速度通量为零来实现。在离心式心脏泵内部流场模拟中，有限体积法的应用过程如下：首先，根据离心式心脏泵的几何模型，将其内部流场区域划分为合适的控制体积，通常采用结构化网格或非结构化网格进行划分。对于复杂的叶轮和泵壳结构，可能会采用非结构化网格，以更好地贴合几何形状，提高网格质量。然后，对每个控制体积应用守恒方程，计算控制体积表面的通量，从而得到离散的代数方程组。接着，通过迭代求解这些代数方程组，得到流场中各控制体积中心节点上的物理量，如速度、压力等。在迭代过程中，需要采用合适的数值求解算法，如SIMPLE算法（Semi-ImplicitMethodforPressure-LinkedEquations）及其改进算法，以确保计算的收敛性和稳定性。最后，根据计算得到的离散解，通过插值等方法得到整个流场的物理量分布，从而实现对离心式心脏泵内部流场的数值模拟。3.4.2有限差分法有限差分法（FiniteDifferenceMethod，FDM）是一种较为经典的数值计算方法，它通过将求解区域划分为网格，用差商来近似代替微商，从而将连续的控制方程离散化为代数方程组。以一维非定常对流扩散方程为例，其一般形式为：\frac{\partial\phi}{\partialt}+u\frac{\partial\phi}{\partialx}=\Gamma\frac{\partial^{2}\phi}{\partialx^{2}}其中，\phi为待求解的物理量，u为速度，\Gamma为扩散系数。在有限差分法中，将时间和空间进行离散，用有限个离散点上的函数值来近似表示连续函数。对于时间导数\frac{\partial\phi}{\partialt}，可以采用向前差分、向后差分或中心差分等方法进行近似，如向前差分公式为\frac{\partial\phi}{\partialt}\approx\frac{\phi_{i}^{n+1}-\phi_{i}^{n}}{\Deltat}，其中\phi_{i}^{n}表示在第n个时间步、第i个空间节点上的函数值，\Deltat为时间步长；对于空间导数\frac{\partial\phi}{\partialx}和\frac{\partial^{2}\phi}{\partialx^{2}}，也可以采用类似的差分公式进行近似。通过将这些差分公式代入原方程，就可以得到离散的代数方程组。有限差分法具有计算简单、易于编程实现的特点，在早期的数值模拟中得到了广泛应用。它对于简单几何形状和规则网格的问题能够快速得到数值解，并且计算效率较高。在一些简单的管道流动模拟中，有限差分法可以快速准确地计算出速度和压力分布。该方法在处理边界条件时相对直观，通过在边界节点上直接设定相应的差分公式，就可以实现对边界条件的处理。与有限体积法相比，有限差分法也存在一些局限性。有限差分法的精度对网格的依赖性较强，在复杂几何形状下，为了保证计算精度，往往需要采用非常细密的网格，这会导致计算量急剧增加，计算效率降低。对于离心式心脏泵内部复杂的流道结构，使用有限差分法进行网格划分时，很难保证网格的质量和均匀性，从而影响计算精度。有限差分法在守恒性方面不如有限体积法严格，在处理一些需要严格满足守恒定律的问题时，可能会出现一定的误差。在模拟离心式心脏泵内部血液流动时，守恒性对于保证模拟结果的准确性至关重要，有限差分法在这方面的不足使其在该领域的应用受到一定限制。3.5流场数值计算求解方法在离心式心脏泵内部流场的数值模拟中，选择合适的流场数值计算求解方法至关重要，它直接关系到模拟结果的准确性和计算效率。常用的求解方法包括SIMPLE算法及其改进算法SIMPLEC算法等，这些算法在处理压力和速度耦合问题上各有特点，下面将详细介绍这两种算法。3.5.1SIMPLE算法SIMPLE算法，全称为Semi-ImplicitMethodforPressure-LinkedEquations，即求解压力耦合方程的半隐方法，由Patankar和Spalding于1972年提出。该算法提出后，迅速成为计算不可压流场的主要方法，并随后成功推广到可压缩流场计算中，现已成为一种能够计算任何流速流动的通用数值方法。在求解压力和速度耦合问题时，SIMPLE算法具有独特的步骤。由于动量方程中压力项以源项形式出现且无独立方程，需要设计专门算法来改进压力值。该算法的基本假设是速度场和压力场的假定相互独立，对假定压力场的修正通过已求解速度场的质量守恒条件来实现。具体计算步骤如下：首先，假定一个初始速度分布，记为u_0，v_0，利用此速度分布计算动量离散方程中的系数及常数项。这一步骤为后续的计算提供了初始条件，虽然是假定的速度分布，但对整个计算过程的启动至关重要。接着，假设一个压力场p^*，这个假设的压力场是后续计算的基础，它将在迭代过程中不断得到修正。然后，依次求解动量方程，得到速度场的估计值u^*，v^*。在求解动量方程时，会用到前面计算得到的系数和常数项，以及假设的压力场。之后，求解压力修正方程，得到压力修正值p'。压力修正方程的求解是SIMPLE算法的关键步骤之一，通过求解该方程，可以得到更准确的压力修正值，从而改进速度场。根据p'改进速度场，使速度场更加符合实际情况。利用改进后的速度场求解那些通过源项物性等与速度场耦合的变量，如果这些变量不影响流场，则应在速度场收敛后再求解。利用改进后的速度场重新计算动量离散方程的系数，并利用改进后的压力场作为下一层次迭代计算的初值，重复上述步骤，直到获得收敛的解。在迭代过程中，需要不断调整计算参数，以确保计算的收敛性和稳定性。在离心式心脏泵内部流场模拟中，SIMPLE算法存在一定的局限性。由于在速度修正中忽略了邻近速度修正值的影响，虽然不影响最后收敛的解，但会加重压力修正值p'的负担，导致整个速度场的迭代速度减慢。这在处理离心式心脏泵内部复杂的流场时，可能会增加计算时间和计算成本。为了克服这一问题，通常对p'作亚松弛处理，一般亚松弛因子可取0.8左右；类似地，速度修正值也可考虑亚松弛，一般松弛因子可取0.5左右。但这些参数的选择需要根据具体的模拟情况进行调整，增加了计算的复杂性。3.5.2SIMPLEC算法SIMPLEC算法是对SIMPLE算法的一种改进，其全称是Semi-ImplicitMethodforPressure-LinkedEquations-Consistent，即改进的求解压力耦合方程的半隐方法。该算法主要是为了使速度场的改进进程与压力场的改进进程能够更好地同步进行。SIMPLEC算法对SIMPLE算法的改进主要体现在压力修正方程的推导和速度修正的过程中。在SIMPLE算法中，速度修正时忽略了不同位置速度修正量之间的影响，而SIMPLEC算法通过对压力修正方程的重新推导，考虑了更多的因素，使得速度场和压力场的修正更加合理。在推导压力修正方程时，SIMPLEC算法对一些项进行了近似处理，使得方程的求解更加高效，同时也提高了计算结果的准确性。与SIMPLE算法相比，SIMPLEC算法具有明显的优势。由于其对压力修正方程的改进，使得速度场和压力场的耦合更加紧密，计算收敛速度更快。在处理复杂的离心式心脏泵内部流场时，SIMPLEC算法能够更快地得到收敛解，从而节省计算时间和计算资源。它还可以采用更大的松弛因子，进一步提高计算效率。在一些情况下，SIMPLEC算法的收敛速度比SIMPLE算法快数倍，这对于大规模的数值模拟计算来说，具有重要的意义。在离心式心脏泵内部流场模拟中，SIMPLEC算法能够更准确地捕捉流场中的压力和速度分布，对于分析泵的性能和优化泵的设计具有重要的帮助。它可以更精确地预测叶轮与泵壳之间的间隙流场特性，以及泵内的压力脉动和流动损失，为提高离心式心脏泵的效率和血液相容性提供更可靠的数据支持。四、定常数值模拟4.1建模过程4.1.1软件选择在构建离心式心脏泵模型时，本研究选用了Creo软件。Creo软件由美国参数技术公司（PTC）开发，是一款功能强大的三维计算机辅助设计（CAD）软件。它在机械设计、工业设计、模具设计等众多领域都有广泛应用，尤其在构建复杂几何模型方面表现出色，这使其成为构建离心式心脏泵模型的理想选择。与其他常见建模软件如Pro/E相比，Creo具有显著的优势。Pro/E自1988年问世以来，在三维建模软件领域占据重要地位，其参数化设计功能曾引领设计风潮。随着行业需求的不断升级，PTC公司在2010年将Pro/E更名为Creo，并进行了重大改版，此次改版不仅是名称的变化，更是技术和产品理念的革新。Creo以其开放、灵活的设计环境吸引了众多用户。在参数化设计方面，虽然Pro/E也具备这一功能，但Creo在此基础上进行了优化和拓展，使得参数调整更加便捷高效。通过智能推理、智能捕捉等功能，设计师在调整模型参数时，软件能够自动根据相关约束和规则进行关联调整，大大减少了手动干预的工作量，提高了设计效率。在构建离心式心脏泵模型时，当调整叶轮的叶片角度、数量或泵壳的形状等参数时，Creo能够快速准确地更新整个模型，确保模型的一致性和准确性。在功能模块方面，Pro/E的设计较为一体化，覆盖了参数化功能定义、实体零件及装配建模等众多功能。然而，随着市场需求的多样化，Creo转向了模块化设计，推出了诸多可互操作的应用程序。用户可以根据不同的项目需求进行相应的功能组合，这种灵活性在构建离心式心脏泵模型时尤为重要。在进行离心式心脏泵的设计时，可能需要同时运用到零件设计、装配设计、分析模拟等多个功能模块，Creo的模块化设计允许用户根据实际需求灵活选择和组合这些功能，避免了不必要的功能冗余，提高了设计的针对性和效率。Creo的用户界面设计更符合现代用户的习惯，合理的工具布局加上可定制的工作区，使得用户能够迅速适应。对于初次接触该软件的用户来说，其简洁直观的操作界面能够降低学习成本，提高操作的便捷性。在构建离心式心脏泵模型的过程中，设计师可以根据自己的操作习惯和工作流程，对工作区进行个性化定制，将常用的工具和命令放置在易于访问的位置，进一步提高建模效率。4.1.2模型构建步骤在构建离心式心脏泵三维模型时，使用Creo软件进行建模，具体步骤如下：设计草图绘制：打开Creo软件，进入草绘模块。根据离心式心脏泵的设计要求和参数，利用软件提供的绘图工具，如直线、圆弧、样条曲线等，绘制泵的各个部件的二维草图。在绘制叶轮草图时，精确绘制轮毂的圆形轮廓，确定其直径和中心位置；根据叶片的设计形状，使用样条曲线工具绘制叶片的截面形状，包括叶片的进口角、出口角以及叶片的弯曲程度等关键参数。在绘制泵壳草图时，根据蜗壳的设计，绘制出螺旋形的流道轮廓，确定流道的起始位置、半径变化规律以及出口位置等。在绘制过程中，严格按照设计尺寸进行绘制，并添加必要的几何约束和尺寸约束，确保草图的准确性和规范性。几何约束可以保证草图中线条之间的相对位置关系，如平行、垂直、相切等；尺寸约束则用于确定草图中各个元素的大小和位置，通过精确输入尺寸数值，确保草图与设计要求一致。三维实体建模：完成草图绘制后，退出草绘模块，进入实体建模模块。利用拉伸、旋转、扫描、混合等实体建模工具，将二维草图转化为三维实体模型。对于叶轮，选择绘制好的叶轮草图，使用旋转工具，以轮毂的中心轴为旋转轴，将草图旋转360度，生成叶轮的三维实体模型。在旋转过程中，设置合适的旋转角度和方向，确保模型的正确性。对于泵壳，根据蜗壳的草图，使用扫描工具，沿着螺旋形的轨迹线扫描草图截面，生成蜗壳的三维实体模型。在扫描过程中，准确选择轨迹线和草图截面，调整扫描的参数，如扫描的起始位置、终止位置、扫描的方式等，以获得符合设计要求的蜗壳模型。通过这些建模工具的灵活运用，逐步构建出离心式心脏泵的各个部件的三维实体模型。部件装配：将构建好的各个部件的三维实体模型进行装配，形成完整的离心式心脏泵模型。在装配模块中，使用装配约束工具，如对齐、匹配、插入等，将叶轮、泵壳、进口管、出口管等部件准确地定位和组合在一起。将叶轮的中心轴与泵壳的中心轴进行对齐约束，确保叶轮能够在泵壳内准确地旋转；将进口管和出口管分别与泵壳的进口和出口进行匹配和插入约束，保证管道与泵壳之间的连接紧密且位置准确。在装配过程中，仔细检查各个部件之间的相对位置和装配关系，确保装配后的模型符合设计要求和实际工作情况。通过合理设置装配约束，可以确保模型的装配精度和稳定性，为后续的数值模拟分析提供准确的模型基础。模型检查与修正：完成装配后，对构建好的离心式心脏泵三维模型进行全面检查。利用Creo软件提供的分析工具，检查模型的几何形状、尺寸精度、装配关系等是否符合设计要求。检查叶轮与泵壳之间的间隙是否均匀，进口管和出口管的连接是否顺畅，模型中是否存在干涉或不合理的结构等问题。如果发现模型存在问题，及时进行修正。对于几何形状不符合要求的部分，可以返回草图或实体建模模块进行修改；对于装配关系不正确的部件，重新调整装配约束，确保模型的质量和准确性。通过严格的模型检查和修正，可以提高模型的可靠性，为后续的数值模拟提供高质量的模型，从而保证模拟结果的准确性和可靠性。四、定常数值模拟4.2网格划分4.2.1划分策略（结构化、非结构化网格）在对离心式心脏泵内部流场进行数值模拟时，网格划分是一个至关重要的环节，其质量直接影响到模拟结果的准确性和计算效率。结构化网格和非结构化网格是两种常见的网格划分类型，它们各自具有独特的特点，在离心式心脏泵流场模拟中有着不同的适用性。结构化网格是一种规则的网格划分方式，其节点在空间上按照一定的规律排列，形成整齐的网格结构。在二维情况下，结构化网格通常表现为矩形或四边形网格；在三维情况下，则多为六面体网格。这种网格的优点显著，由于节点排列规则，数据存储和计算相对简单，能够有效地提高计算效率。在进行数值计算时，结构化网格可以利用其规则性，采用高效的数值算法，减少计算量和计算时间。结构化网格的节点编号具有规律性，便于程序的编写和实现，能够降低编程的难度和出错的概率。在一些简单几何形状的流场模拟中，如直管内的流动，结构化网格能够快速生成高质量的网格，准确地模拟流场特性。然而，结构化网格也存在一定的局限性。当面对复杂的几何形状时，如离心式心脏泵的叶轮和泵壳，由于其形状不规则，边界条件复杂，结构化网格的生成难度较大。为了适应复杂的几何形状，可能需要对网格进行大量的局部加密和调整，这不仅会增加网格生成的工作量，还容易导致网格质量下降。在叶轮的叶片与轮毂连接处，以及泵壳的蜗形流道等部位，结构化网格难以精确地贴合几何边界，可能会产生较大的网格误差，影响模拟结果的准确性。非结构化网格则是一种更加灵活的网格划分方式，它不受节点排列规则的限制，可以根据几何形状的特点自由地生成网格。非结构化网格可以由三角形、四面体、多边形等多种形状的单元组成，能够很好地适应复杂的几何形状。在离心式心脏泵的流场模拟中，非结构化网格能够精确地贴合叶轮和泵壳的复杂边界，在边界层和几何形状变化剧烈的区域，通过局部加密网格，可以提高计算精度。在叶轮的叶片表面和泵壳的内壁等关键部位，非结构化网格可以生成细密的网格，准确地捕捉流场的细节信息，从而更准确地模拟流场特性。非结构化网格在处理边界条件时也具有很大的优势，它可以根据边界的形状和要求，灵活地设置网格节点和单元，使边界条件的处理更加准确和方便。在离心式心脏泵的进口和出口边界，非结构化网格可以根据实际的流动情况，合理地布置网格，准确地设定边界条件，提高模拟结果的可靠性。非结构化网格也存在一些缺点，由于其网格单元形状和大小的多样性，数据存储和计算相对复杂，计算效率相对较低。在进行数值计算时，非结构化网格需要处理更多的网格信息，导致计算量增加，计算时间延长。非结构化网格的生成算法相对复杂，需要较高的技术水平和经验，才能生成高质量的网格。在离心式心脏泵内部流场模拟中，需要根据具体的情况选择合适的网格划分策略。对于几何形状相对简单的部分，如进口管和出口管，可以采用结构化网格，以提高计算效率；而对于叶轮和泵壳等几何形状复杂的部分，则采用非结构化网格，以确保网格能够精确地贴合几何边界，提高计算精度。在实际应用中，还可以采用混合网格的方式，将结构化网格和非结构化网格结合起来，充分发挥它们的优势，以获得更好的模拟效果。在叶轮的轮毂部分采用结构化网格，在叶片表面和泵壳内采用非结构化网格，通过合理的过渡和连接，实现两种网格的协同工作，既能保证计算效率，又能提高计算精度。4.2.2网格质量检查与优化在完成离心式心脏泵内部流场的网格划分后，对网格质量进行检查与优化是确保数值模拟结果准确性和可靠性的关键步骤。网格质量直接影响到计算的精度、收敛性和稳定性，因此需要采用一系列指标来检查网格质量，并运用相应的方法进行优化。检查网格质量的指标丰富多样，其中单元形状质量是一个重要指标。理想的单元形状应尽量接近正多边形或正多面体，以保证计算的准确性。在二维网格中，三角形单元的内角应尽量接近60°，四边形单元应尽量接近正方形；在三维网格中，四面体单元的各面应尽量接近正三角形，六面体单元应尽量接近正方体。单元形状质量可以通过计算单元的长宽比、内角比等参数来评估。长宽比过大的单元可能会导致计算精度下降，因为在这种单元中，物理量的变化可能会被过度放大或缩小；内角比不合理的单元则可能会影响计算的稳定性，增加计算过程中出现数值振荡的风险。网格的正交性也是衡量网格质量的重要因素。正交性良好的网格，其单元的边或面之间夹角接近90°，这样在数值计算中能够减少数值误差的产生。在计算流体力学中，当网格的正交性较差时，会导致离散方程中的交叉项增多，从而影响计算的精度和收敛性。在离心式心脏泵内部流场的模拟中，如果网格的正交性不好，可能会使速度和压力的计算结果出现偏差，无法准确反映流场的真实情况。网格的纵横比同样不容忽视，它是指单元在不同方向上的尺寸之比。合理的纵横比能够保证物理量在各个方向上的变化得到准确的捕捉。如果纵横比过大或过小，会导致在某些方向上的计算精度降低，影响模拟结果的可靠性。在离心式心脏泵的叶轮附近，由于流场变化剧烈，需要保证网格在各个方向上的尺寸分布合理，以准确捕捉流场的细节信息。如果网格的纵横比不合理，可能会使叶轮附近的速度和压力分布计算不准确，无法准确评估叶轮的性能。针对检查出的网格质量问题，需要采取相应的优化方法。对于形状质量较差的单元，可以采用网格平滑技术进行优化