基于CFD的重力热管流动与传热特性数值模拟及验证研究

基于CFD的重力热管流动与传热特性数值模拟及验证研究一、引言1.1研究背景与意义随着全球工业化进程的加速，能源问题已成为世界各国关注的焦点。一方面，能源需求持续增长，传统能源储量却日益减少，能源供需矛盾日益突出；另一方面，大量使用传统能源所带来的环境污染问题也愈发严重，如温室气体排放导致的全球气候变暖等。因此，提高能源利用效率、开发和利用清洁能源成为解决能源问题的关键。在能源利用过程中，传热是一个至关重要的环节。高效的传热技术能够实现能源的有效传递和利用，降低能源消耗和成本。重力热管作为一种高效的传热元件，在能源领域展现出了巨大的应用潜力。重力热管是一种依靠重力使工作液体回流的热管，其结构简单，通常由管壳、端盖和工作介质组成，内部被抽成负压状态并充入适量的工作液体。与其他类型的热管相比，重力热管具有诸多独特的优势。首先，它的结构相对简单，不需要复杂的吸液芯结构，这使得其制造工艺更为简便，成本更低。其次，重力热管的传热效率极高，其等效导热系数可比传统金属材料高出数倍甚至数十倍，能够在较小的温差下传递大量的热量。此外，重力热管还具有良好的等温性，能够使热源表面的温度分布更加均匀，有效降低热点温度，提高设备的可靠性和稳定性。在能源领域，重力热管的应用十分广泛。在太阳能热发电系统中，重力热管可用于高效收集和传输太阳能，将太阳能转化为热能，提高系统的发电效率。在余热回收领域，重力热管能够有效地回收工业生产过程中产生的废热，将其转化为有用的能源，实现能源的二次利用，降低能源消耗和生产成本。在新能源汽车的电池热管理系统中，重力热管可用于控制电池的温度，确保电池在适宜的温度范围内工作，提高电池的性能和寿命，保障新能源汽车的安全运行。深入研究重力热管的流动与传热特性具有重要的理论和实际意义。从理论层面来看，重力热管内部的流动和传热过程涉及到复杂的多相流和相变传热现象，目前相关理论尚未完全成熟，存在许多有待深入研究的问题。通过对重力热管流动与传热特性的研究，可以进一步揭示其传热机理，完善多相流和相变传热理论，为热管的优化设计和性能提升提供坚实的理论基础。从实际应用角度出发，准确掌握重力热管的流动与传热特性，有助于优化其结构和工作参数，提高其传热效率和可靠性，从而推动重力热管在能源领域以及其他相关领域的更广泛应用，促进能源的高效利用和可持续发展，为解决当前的能源问题做出贡献。1.2国内外研究现状自20世纪60年代热管技术诞生以来，国内外学者对重力热管展开了大量的研究，研究内容涵盖了实验研究和数值模拟等多个方面。在实验研究方面，早期的研究主要集中在重力热管的基本传热性能测试上。国外学者如[具体学者1]等搭建了重力热管实验平台，通过测量不同工况下热管的温度分布、传热功率等参数，对重力热管的传热性能进行了初步探究，发现了重力热管在不同充液率、热负荷等条件下的传热规律。国内学者[具体学者2]等也进行了类似的实验研究，进一步验证和补充了相关结论，并对实验数据进行了更深入的分析。随着研究的深入，实验研究逐渐向多因素、多参数方向发展。国外学者[具体学者3]研究了工作介质、管径、倾角等因素对重力热管传热性能的综合影响，发现工作介质的热物理性质、管径的大小以及热管的倾角会显著影响热管的传热功率和热阻。国内学者[具体学者4]则通过实验研究了不同加热方式、冷却条件下重力热管的启动特性和稳态传热特性，为重力热管的实际应用提供了更丰富的实验数据。在数值模拟方面，随着计算机技术和计算流体力学（CFD）的发展，数值模拟成为研究重力热管流动与传热特性的重要手段。国外学者[具体学者5]等采用VOF（体积模型）来模拟重力热管中汽-液两相的相互作用阶段，通过建立二维模型，对热管内的流动和传热过程进行了数值模拟，并将模拟结果与实验对比，显示模拟的温度分布与实验测得的壁温分布吻合性良好。国内学者[具体学者6]在前人关于重力热管实验和理论研究的基础上，在二维柱坐标系下建立了描述重力热管内部流动和传热过程的数学模型，采用求解压力耦合方程组的半隐式方法（SIMPLE算法）编写C语言计算机程序进行求解，分析了重力热管工作温度、蒸发段热流密度等参数变化对其流动和传热性能的影响。尽管国内外在重力热管的研究方面取得了一定的成果，但仍存在一些不足和空白。一方面，目前的研究大多集中在单一因素对重力热管性能的影响上，对于多因素耦合作用下重力热管的流动与传热特性研究较少。另一方面，在数值模拟中，部分模型的假设条件与实际情况存在一定差异，导致模拟结果的准确性和可靠性有待进一步提高。此外，对于一些新型重力热管或特殊工况下的重力热管，相关研究还较为匮乏。1.3研究内容与方法本文主要研究内容为构建重力热管的数值模型，利用CFD技术对其内部的流动与传热特性进行模拟分析，并通过实验对模拟结果进行验证。具体包括以下几个方面：建立重力热管模型：根据重力热管的实际结构和工作原理，考虑到热管内部的复杂流动和传热过程，建立合理的三维物理模型。在建模过程中，充分考虑热管的几何参数、工作介质的热物理性质以及边界条件等因素。数值模拟计算：采用CFD软件对建立的模型进行数值模拟计算。选择合适的多相流模型来描述热管内的汽液两相流动，如VOF模型、混合物模型或欧拉模型等，并结合传热学相关理论，对热管内的传热过程进行模拟。通过模拟，得到热管内部的速度场、温度场、压力场等分布情况，以及不同工况下热管的传热性能参数，如传热功率、热阻等。实验研究：搭建重力热管实验平台，制备实验用重力热管，确保其尺寸、材料等参数与数值模型一致。通过实验测量不同工况下热管的温度分布、传热功率等数据，并与数值模拟结果进行对比分析。实验过程中，严格控制实验条件，保证实验数据的准确性和可靠性。结果分析与讨论：对数值模拟结果和实验结果进行深入分析，探讨重力热管的流动与传热特性，研究各因素对热管性能的影响规律。通过对比模拟结果和实验结果，验证数值模型的准确性和可靠性，进一步优化数值模型和计算方法。本文采用CFD技术和实验研究相结合的方法。CFD技术具有成本低、效率高、能够深入揭示内部流动和传热细节等优点，但数值模拟结果的准确性依赖于模型的合理性和假设条件的可靠性。实验研究则能够提供真实可靠的数据，验证数值模拟结果的准确性，但实验过程往往受到实验条件、测试手段等因素的限制。将两者结合，可以充分发挥各自的优势，相互补充和验证，从而更全面、准确地研究重力热管的流动与传热特性。二、重力热管基本原理与理论基础2.1重力热管工作原理重力热管主要由管壳、端盖和工作介质组成，内部被抽成负压状态并充入适量的工作液体，其结构通常可分为蒸发段、绝热段和冷凝段。当重力热管的蒸发段与高温热源接触时，工作介质吸收热量，从液态转变为气态，这一过程称为蒸发。由于蒸发过程中工作介质吸收了大量的汽化潜热，使得蒸发段的温度基本保持在工作介质的沸点温度。产生的蒸汽在微小的压力差作用下，迅速从蒸发段流向冷凝段。由于蒸汽的密度远小于液体，在重力的作用下，蒸汽上升的速度较快，能够快速地将热量传递到冷凝段。在冷凝段，蒸汽与低温冷源接触，将热量释放给冷源，自身重新凝结成液态，这一过程称为冷凝。冷凝过程中，蒸汽放出大量的潜热，使得冷凝段的温度也基本保持在工作介质的沸点温度。凝结后的液体在重力的作用下，沿着管壳内壁回流至蒸发段，再次吸收热量并蒸发，如此循环往复，实现了热量从高温端到低温端的持续传递。绝热段位于蒸发段和冷凝段之间，其作用是减少热量在传递过程中的损失。绝热段通常采用导热系数较低的材料制成，或者对其进行良好的保温处理，以确保热量主要在蒸发段和冷凝段之间传递，而不会在绝热段散失到周围环境中。通过这种独特的工作方式，重力热管能够在较小的温差下实现高效的热量传递，其等效导热系数可比传统金属材料高出数倍甚至数十倍，具有极高的传热效率。同时，由于重力热管内部的工质循环是依靠重力和相变过程实现的，无需外加动力装置，因此具有结构简单、运行可靠、维护方便等优点。2.2相关理论基础2.2.1传热学理论重力热管的传热过程涉及多种传热方式，主要包括导热、对流换热和相变传热。在管壳和工作介质内部，热量传递主要以导热的方式进行。导热是指物体内部或相互接触的物体之间，由于温度差而引起的能量传递现象。根据傅里叶定律，导热的基本公式为：q=-\lambda\frac{\partialT}{\partialx}其中，q为热流密度，\lambda为导热系数，\frac{\partialT}{\partialx}为温度梯度。导热系数\lambda是材料的固有属性，表示材料传导热量的能力，不同材料的导热系数差异较大。对于重力热管的管壳材料，通常选择导热性能良好的金属，如铜、铝等，以减少管壳的热阻，提高热量传递效率；而工作介质的导热系数也会影响热管的传热性能，一般来说，工质的导热系数越高，在蒸发和冷凝过程中热量传递就越迅速。在蒸汽和液体的流动过程中，存在对流换热现象。对流换热是指流体与固体壁面之间，由于温度差和流体的宏观运动而引起的热量传递过程。对流换热的强度与流体的性质、流速、固体壁面的形状和粗糙度等因素有关。对流换热的基本公式为牛顿冷却公式：q=h(T_w-T_f)其中，h为对流换热系数，T_w为固体壁面温度，T_f为流体温度。对流换热系数h反映了对流换热的强弱程度，它不是一个物性参数，而是受多种因素影响的综合系数。在重力热管中，蒸汽在管内的流动速度、液体在管壁上的流动状态等都会对对流换热系数产生影响。例如，蒸汽流速较高时，对流换热系数会增大，有利于热量的传递；而液体在管壁上形成的液膜厚度和流动稳定性也会影响对流换热效果，如果液膜过厚或不稳定，会增加热阻，降低对流换热效率。重力热管中最重要的传热方式是相变传热，即工作介质在蒸发段吸收热量蒸发和在冷凝段放出热量凝结的过程。相变传热过程中，工质吸收或释放大量的潜热，使得热管能够在较小的温差下传递大量的热量。相变传热的热流量可通过下式计算：Q=m\cdotr其中，Q为相变传热量，m为相变的工质质量，r为工质的汽化潜热。汽化潜热是工质的重要热物理性质之一，不同的工质具有不同的汽化潜热，一般来说，汽化潜热越大，在相同质量的工质相变时传递的热量就越多。在重力热管的设计和应用中，选择合适的工作介质以充分利用其相变传热特性，是提高热管传热性能的关键因素之一。2.2.2流体力学理论在重力热管内，流体的流动对传热性能有着重要影响。重力热管内的流体流动主要包括蒸汽的上升流动和液体的回流。蒸汽在蒸发段产生后，由于密度差在管内向上流动。蒸汽的流动状态受到多种因素的影响，如热管的管径、蒸汽的生成速率、管内压力分布等。根据流体力学中的连续性方程和伯努利方程，可以分析蒸汽的流动特性。连续性方程表示在稳态流动中，单位时间内通过管道任意截面的流体质量相等，其表达式为：\rho_1v_1A_1=\rho_2v_2A_2其中，\rho为流体密度，v为流体流速，A为管道横截面积。在重力热管中，蒸汽在流动过程中，由于压力和温度的变化，其密度也会发生改变，通过连续性方程可以分析蒸汽流速与管径、密度之间的关系。例如，当管径不变时，蒸汽生成速率增加，根据连续性方程，蒸汽流速会相应增大。伯努利方程则描述了理想流体在稳态流动过程中的能量守恒关系，对于不可压缩流体，其表达式为：p+\frac{1}{2}\rhov^2+\rhogz=C其中，p为流体压力，\frac{1}{2}\rhov^2为动能项，\rhogz为重力势能项，C为常数。在重力热管中，利用伯努利方程可以分析蒸汽在上升过程中压力、速度和高度之间的关系。例如，蒸汽在上升过程中，高度增加，重力势能增大，如果忽略摩擦阻力，根据伯努利方程，蒸汽的压力和速度会发生相应的变化，以保持总能量守恒。液体在冷凝段凝结后，依靠重力沿管壁回流至蒸发段。液体的回流受到重力、管壁的摩擦力以及蒸汽与液体之间的相互作用力等因素的影响。当液体在管壁上形成液膜向下流动时，液膜的厚度和流动稳定性对传热性能有重要影响。如果液膜过厚，会增加液体回流的阻力，同时也会增大热阻，降低传热效率；而液膜不稳定时，可能会出现波动、破裂等现象，影响液体的正常回流和传热的稳定性。此外，重力热管内的汽液两相流还存在一些特殊的流动现象，如汽液界面的波动、相间的质量和动量传递等。这些现象使得重力热管内的流体力学行为更加复杂，需要综合考虑多种因素来分析和理解流体的流动特性，进而优化重力热管的设计和性能。三、数值模拟方法与模型建立3.1数值模拟软件选择在进行重力热管的数值模拟时，可供选择的CFD软件众多，如FLUENT、ANSYSCFX、OpenFOAM等。这些软件各有特点和优势，在不同的工程领域和研究方向中都有广泛应用。ANSYSCFX是一款功能强大的CFD软件，具有先进的数值算法和高效的求解器，能够处理复杂的几何模型和多物理场耦合问题。它在航空航天、汽车工程等领域应用广泛，对于模拟复杂流场和高精度的计算需求有很好的表现。然而，其软件操作相对复杂，学习成本较高，对于一些初学者来说可能存在一定难度。OpenFOAM是一款开源的CFD软件，具有高度的灵活性和可定制性。用户可以根据自己的需求对求解器进行修改和扩展，这对于一些有特定研究需求和算法开发的科研人员来说具有很大的吸引力。但由于其开源特性，软件的稳定性和技术支持相对较弱，在使用过程中可能需要用户具备较强的编程能力和CFD理论基础。FLUENT软件则凭借其丰富的物理模型、友好的用户界面和广泛的应用领域，成为了本次重力热管数值模拟的首选软件。它拥有多种湍流模型、多相流模型以及传热模型，能够准确地模拟重力热管内复杂的汽液两相流动和传热过程。例如，在模拟重力热管中汽液界面的运动和相变传热时，FLUENT提供的VOF模型、混合物模型等可以很好地描述汽液两相的相互作用和界面变化。其用户界面简洁直观，前处理和后处理功能强大，方便用户进行模型的构建、网格划分以及结果的可视化分析。而且，FLUENT在国内外的科研和工程领域都有大量的应用案例和丰富的文献资料，当遇到问题时，用户可以很容易地找到相关的参考和解决方案。综上所述，考虑到重力热管数值模拟的复杂性和对模型准确性、易用性的要求，选择FLUENT软件能够更好地满足研究需求，为深入分析重力热管的流动与传热特性提供有力支持。3.2模型假设与简化由于重力热管内部的流动和传热过程极其复杂，涉及到汽液两相的相变、界面的波动以及复杂的流固耦合等现象，为了便于建立数学模型并进行数值求解，需要对实际物理模型进行合理的假设与简化。假设热管处于稳态工作状态，管内流体为稳态流动。在实际运行中，重力热管在启动阶段会经历一个非稳态过程，但当达到稳定运行状态后，其内部的流动和传热特性基本保持不变。通过这一假设，可以忽略启动阶段的复杂瞬态变化，简化计算过程，更专注于研究热管在稳定工况下的性能。假设热物理常数不随温度变化而变化。尽管工作介质的热物理性质，如密度、粘度、导热系数等，在实际运行中会随着温度的变化而有所改变，但在一定的温度范围内，这种变化相对较小。为了简化计算，将热物理常数视为定值，在保证计算精度的前提下，大大减少了计算的复杂性。设定蒸发段和冷凝段外表面为等热流密度边界条件，热管两端面为绝热面。这一假设是基于实际应用中，重力热管通常与外界热源和冷源通过特定的换热方式进行热量交换，在一定的研究范围内，可以近似认为蒸发段和冷凝段外表面的热流密度均匀分布。而热管两端面由于通常进行了良好的保温处理，热量散失极少，可将其视为绝热面，以减少计算域的边界条件复杂性，提高计算效率。假定汽、液两相为饱和流体且为连续流动，充液率适当，无液池存在或存在一个高度较小的液池。在重力热管正常工作时，工作介质在蒸发段和冷凝段分别处于饱和蒸汽和饱和液体状态，且流动较为连续。适当的充液率是保证热管正常运行的关键因素之一，当充液率合适时，管内不会出现大量的液体积聚形成液池的情况，或者即使存在液池，其高度也较小，对整体流动和传热的影响可以忽略不计。假设汽液流体流动为层流流动，而且蒸汽为亚音速流动。在一些工况下，重力热管内的汽液流动可能处于层流状态，特别是在低负荷或管径较小的情况下。同时，由于蒸汽在管内的流速相对较低，一般可认为蒸汽为亚音速流动。这一假设符合许多实际应用场景，能够简化流体力学方程的求解，避免了复杂的湍流模型和超音速流动计算。将流动与传热问题简化为二维轴对称问题，即蒸气流速和压力及液体流速和压力仅是轴向位置的函数。对于一些长径比较大且结构对称的重力热管，在忽略径向变化的情况下，将其简化为二维轴对称模型可以大大减少计算量，同时又能较好地反映热管内的主要流动和传热特性。通过这一假设，将三维问题转化为二维问题，降低了计算难度，提高了计算效率，并且在一定程度上能够获得较为准确的模拟结果。通过以上合理的假设与简化，能够在保证一定计算精度的前提下，大大降低重力热管数值模拟的复杂性，为后续的数值计算和分析奠定基础。3.3控制方程与边界条件3.3.1控制方程在重力热管的数值模拟中，需要建立一系列的控制方程来描述管内流体的流动和传热过程。这些控制方程基于质量守恒、动量守恒和能量守恒定律，能够准确地反映热管内的物理现象。连续性方程是质量守恒定律在流体力学中的体现，它表示在单位时间内，流入和流出控制体的质量差等于控制体内质量的变化率。对于重力热管内的汽液两相流，连续性方程可以分别对蒸汽相和液相进行描述。假设蒸汽相的密度为\rho_v，速度矢量为\vec{v}_v，液相的密度为\rho_l，速度矢量为\vec{v}_l，则蒸汽相的连续性方程为：\frac{\partial\rho_v}{\partialt}+\nabla\cdot(\rho_v\vec{v}_v)=0液相的连续性方程为：\frac{\partial\rho_l}{\partialt}+\nabla\cdot(\rho_l\vec{v}_l)=0式中，t为时间，\nabla为哈密顿算子，表示对空间坐标的求导。动量方程是动量守恒定律的数学表达，它描述了作用在流体微元上的各种力与流体微元动量变化之间的关系。对于不可压缩流体，动量方程可以简化为纳维-斯托克斯方程（N-S方程）。在重力热管中，考虑到蒸汽和液体的流动以及重力的作用，动量方程为：\rho\left(\frac{\partial\vec{v}}{\partialt}+(\vec{v}\cdot\nabla)\vec{v}\right)=-\nablap+\nabla\cdot\tau+\rho\vec{g}其中，\rho为流体密度（对于汽液两相流，分别为蒸汽相密度\rho_v和液相密度\rho_l），\vec{v}为速度矢量（分别为蒸汽相速度\vec{v}_v和液相速度\vec{v}_l），p为压力，\tau为粘性应力张量，\vec{g}为重力加速度矢量。粘性应力张量\tau与流体的动力粘度\mu和速度梯度有关，对于牛顿流体，其表达式为\tau=\mu\left(\nabla\vec{v}+(\nabla\vec{v})^T\right)，其中(\nabla\vec{v})^T为速度梯度张量的转置。能量方程是能量守恒定律在传热学中的体现，它描述了单位时间内控制体内能量的变化与流入、流出控制体的能量以及控制体内热源之间的关系。在重力热管中，能量方程考虑了流体的内能、动能以及由于热传导和对流引起的热量传递。假设流体的比内能为e，导热系数为\lambda，比热容为c_p，则能量方程为：\rho\left(\frac{\partiale}{\partialt}+(\vec{v}\cdot\nabla)e\right)=\nabla\cdot(\lambda\nablaT)-p\nabla\cdot\vec{v}+\Phi式中，T为温度，\Phi为粘性耗散项，表示由于流体粘性作用而产生的机械能转化为热能的部分。对于理想气体，比内能e与温度T的关系为e=c_vT，其中c_v为定容比热容，且c_p-c_v=R，R为气体常数。在重力热管的蒸发段和冷凝段，还存在相变传热过程，这需要在能量方程中考虑相变潜热的影响。当蒸汽在冷凝段凝结时，会释放出汽化潜热r，此时能量方程中应增加相变潜热项\rho_lr\frac{\partial\alpha_l}{\partialt}，其中\alpha_l为液相体积分数；当液体在蒸发段蒸发时，会吸收汽化潜热，能量方程中的相变潜热项则为-\rho_vr\frac{\partial\alpha_v}{\partialt}，\alpha_v为蒸汽相体积分数。此外，为了描述重力热管内的汽液两相流，还需要引入相体积分数方程。相体积分数表示某一相在混合流体中所占的体积比例，对于汽液两相流，蒸汽相体积分数\alpha_v和液相体积分数\alpha_l满足\alpha_v+\alpha_l=1。相体积分数方程可以通过追踪汽液界面的运动来求解，常用的方法如VOF模型中，通过求解相体积分数的对流扩散方程来确定汽液界面的位置和形状。这些控制方程相互耦合，共同描述了重力热管内复杂的流动和传热现象。在数值模拟过程中，需要采用合适的数值方法对这些方程进行离散和求解，以获得热管内的速度场、温度场、压力场以及相分布等信息。3.3.2边界条件设置边界条件的合理设置对于准确模拟重力热管的流动与传热特性至关重要，它直接影响到数值计算的结果和收敛性。根据重力热管的实际工作情况，主要需要设置蒸发段、冷凝段和绝热段的边界条件。在蒸发段，假设外表面为等热流密度边界条件，即给定蒸发段外表面单位面积上的热流密度q_{evap}。这一热流密度表示外界热源通过蒸发段管壁传递给工作介质的热量，其大小取决于实际的加热工况。根据傅里叶定律，蒸发段外表面的热流密度与管壁温度梯度的关系为q_{evap}=-\lambda_w\frac{\partialT}{\partialn}，其中\lambda_w为管壁材料的导热系数，\frac{\partialT}{\partialn}为管壁外表面法向温度梯度。在数值模拟中，将q_{evap}作为已知条件施加在蒸发段外表面边界上，用于确定蒸发段内工作介质的温度和热流分布。对于冷凝段，同样假设外表面为等热流密度边界条件，给定冷凝段外表面单位面积上的热流密度q_{cond}，它表示工作介质在冷凝段通过管壁传递给外界冷源的热量。与蒸发段类似，q_{cond}=-\lambda_w\frac{\partialT}{\partialn}，在数值计算中，将q_{cond}施加在冷凝段外表面边界上，以求解冷凝段内的温度场和热流传递情况。需要注意的是，q_{cond}的大小与蒸发段的热流密度q_{evap}以及热管的传热效率有关，在稳态工况下，忽略热管的热损失，有q_{cond}=q_{evap}。绝热段的主要作用是减少热量在传递过程中的损失，因此假设绝热段外表面与外界环境之间无热量交换，即绝热段外表面的热流密度为零，q_{ins}=0。这意味着在数值模拟中，绝热段外表面的温度梯度为零，\frac{\partialT}{\partialn}=0。通过这一边界条件的设置，能够保证热量主要在蒸发段和冷凝段之间传递，而不会在绝热段散失到周围环境中，从而更准确地模拟重力热管的传热过程。在热管的两端面，由于通常进行了良好的保温处理，热量散失极少，因此将两端面设置为绝热边界条件，即两端面的热流密度为零，q_{end}=0，温度梯度\frac{\partialT}{\partialn}=0。这一条件确保了热量在热管内部沿着轴向进行传递，而不会从两端面泄漏，符合实际的物理情况。对于汽液两相流的边界条件，在蒸发段和冷凝段的汽液界面上，需要考虑相界面的质量、动量和能量传递。在相界面上，蒸汽和液体的温度相等，均为工作介质的饱和温度T_{sat}。根据质量守恒，蒸汽相和液相在相界面上的质量通量相等，即\rho_vv_{v,n}=\rho_lv_{l,n}，其中v_{v,n}和v_{l,n}分别为蒸汽相和液相在相界面法向的速度分量。在动量传递方面，相界面上的剪应力连续，即\tau_{v,n}=\tau_{l,n}，其中\tau_{v,n}和\tau_{l,n}分别为蒸汽相和液相在相界面法向的剪应力分量。在能量传递方面，考虑相变潜热的影响，相界面上的热流密度满足q_{interface}=\rho_lrv_{l,n}-\lambda_v\frac{\partialT}{\partialn}_v=\lambda_l\frac{\partialT}{\partialn}_l，其中\lambda_v和\lambda_l分别为蒸汽相和液相的导热系数，\frac{\partialT}{\partialn}_v和\frac{\partialT}{\partialn}_l分别为蒸汽相和液相在相界面法向的温度梯度。通过合理设置上述边界条件，能够准确地模拟重力热管在不同工况下的流动与传热特性，为深入研究重力热管的性能提供可靠的数值计算基础。在实际模拟过程中，还需要根据具体的研究需求和模型特点，对边界条件进行适当的调整和优化，以获得更准确的模拟结果。3.4网格划分与独立性验证3.4.1网格划分方法网格划分是数值模拟中的关键步骤之一，它将连续的计算域离散化为有限个小的单元，这些单元的集合构成了数值计算的网格。合理的网格划分能够提高计算精度和效率，准确地捕捉物理现象的细节。对于重力热管模型，本文采用结构化网格和非结构化网格相结合的方法进行划分。结构化网格具有规则的拓扑结构，网格单元排列整齐，节点分布均匀。在重力热管的主体部分，由于其几何形状较为规则，采用结构化网格可以提高网格质量和计算效率。例如，对于圆柱形的热管管壳，可以使用六面体结构化网格进行划分，通过合理设置网格尺寸和增长率，能够在保证计算精度的前提下，减少网格数量，降低计算成本。在划分结构化网格时，通常采用映射网格划分技术，将计算域映射到规则的坐标系中，然后在该坐标系中生成网格，再将网格映射回实际的计算域。这种方法能够保证网格的质量和正交性，有利于数值计算的收敛和精度的提高。然而，在重力热管的一些复杂区域，如蒸发段和冷凝段与外界热源、冷源的接触部位，以及管壳与端盖的连接处，采用结构化网格划分可能会遇到困难，难以准确地描述这些区域的几何形状和边界条件。此时，采用非结构化网格进行划分更为合适。非结构化网格的单元形状和大小可以根据几何形状和物理场的变化进行灵活调整，能够更好地适应复杂的几何模型。在这些复杂区域，可以使用四面体或混合网格（如四面体与六面体混合）进行划分，通过局部加密网格，能够更准确地捕捉这些区域的流动和传热特性。例如，在蒸发段与热源接触的壁面附近，由于温度梯度较大，需要对网格进行加密，以提高温度场的计算精度；在管壳与端盖的连接处，由于几何形状的突变，采用非结构化网格可以更好地拟合边界形状，减少计算误差。在使用Gambit软件进行网格划分时，首先根据重力热管的几何模型创建计算域，然后对计算域进行分区处理。对于适合结构化网格划分的区域，设置相应的映射网格参数，生成结构化网格；对于复杂区域，使用非结构化网格划分工具，如TGrid或MeshTool，生成非结构化网格。在划分过程中，需要注意网格的质量指标，如网格的长宽比、雅克比行列式等，确保网格质量满足数值计算的要求。同时，为了提高计算效率，可以采用局部网格细化技术，在物理量变化剧烈的区域（如汽液界面、温度梯度较大的区域）进行网格加密，而在物理量变化平缓的区域适当降低网格密度。通过结构化网格和非结构化网格相结合的方法，能够充分发挥两种网格的优势，既保证了计算域主体部分的计算效率，又提高了复杂区域的计算精度，为准确模拟重力热管的流动四、数值模拟结果与分析4.1流动特性分析4.1.1速度分布通过数值模拟，获得了不同位置和工况下重力热管内流体的速度矢量图和速度分布曲线，以深入分析其速度变化规律。图4.1展示了在加热功率为100W、充液率为30%工况下，重力热管蒸发段、绝热段和冷凝段的速度矢量图。从蒸发段速度矢量图（图4.1a）中可以清晰地看到，靠近管壁处的液体速度较低，而管中心区域的蒸汽速度较高。这是因为在蒸发段，液体吸收热量蒸发成为蒸汽，蒸汽在管内形成向上的流动。由于蒸汽的密度远小于液体，且受到浮力和压力差的作用，蒸汽在管中心区域快速上升，而靠近管壁的液体则受到管壁的摩擦力和粘性力的影响，速度相对较慢。随着离蒸发段入口距离的增加，蒸汽的速度逐渐增大，这是由于更多的液体蒸发产生蒸汽，使得蒸汽的流量增加，从而导致速度增大。在绝热段（图4.1b），蒸汽继续保持向上的流动，速度分布相对较为均匀，但整体速度略有下降。这是因为绝热段虽然没有热量的输入和输出，但蒸汽在流动过程中会与管壁发生摩擦，消耗一部分能量，导致速度有所降低。同时，由于蒸汽在绝热段没有进一步的质量增加，其流量保持相对稳定，因此速度变化相对较小。进入冷凝段（图4.1c）后，蒸汽与冷源接触开始凝结成液体，速度迅速减小。在冷凝段的上部，蒸汽速度仍然较高，但随着蒸汽的不断凝结，速度逐渐降低。靠近冷凝段出口处，蒸汽基本全部凝结成液体，液体在重力的作用下沿管壁向下流动，速度相对较低。此时，液体的速度分布呈现出靠近管壁处速度较大，管中心区域速度较小的特点，这与蒸发段的速度分布情况相反。为了更直观地分析速度变化规律，图4.2给出了不同工况下重力热管轴向速度分布曲线。从图中可以看出，在不同加热功率和充液率下，速度分布曲线具有相似的趋势，但具体数值存在差异。随着加热功率的增加，蒸汽的生成速率加快，管内蒸汽的速度明显增大。这是因为加热功率的提高使得蒸发段内液体蒸发更加剧烈，产生更多的蒸汽，在相同的管径条件下，蒸汽流量的增加必然导致速度的增大。而充液率对速度的影响则相对较为复杂，当充液率较低时，管内液体量较少，蒸汽的生成空间相对较大，蒸汽速度相对较高；随着充液率的增加，管内液体量增多，蒸汽的生成空间受到一定限制，且液体对蒸汽的流动会产生一定的阻碍作用，导致蒸汽速度有所下降。但当充液率过高时，液体可能会在管内形成较大的液柱，阻碍蒸汽的正常流动，甚至出现液泛现象，使得蒸汽速度急剧下降，严重影响热管的传热性能。因此，在实际应用中，需要选择合适的充液率，以确保重力热管内流体的正常流动和良好的传热性能。（a）蒸发段（b）绝热段（c）冷凝段4.1.2压力分布通过数值模拟得到了重力热管内的压力云图和压力沿程变化曲线，以深入探讨压力分布与流动状态的关系。图4.3展示了在加热功率为150W、充液率为40%工况下，重力热管的压力云图。从压力云图中可以明显看出，蒸发段的压力相对较高，这是因为在蒸发段，液体吸收热量蒸发成为蒸汽，蒸汽的产生使得管内分子数增加，根据理想气体状态方程pV=nRT（其中p为压力，V为体积，n为物质的量，R为气体常数，T为温度），在体积和温度一定的情况下，分子数的增加会导致压力升高。随着蒸汽从蒸发段流向冷凝段，压力逐渐降低。在绝热段，由于没有热量和质量的交换，蒸汽的压力主要由于流动阻力而逐渐减小。而在冷凝段，蒸汽与冷源接触凝结成液体，蒸汽的体积迅速减小，分子数减少，导致压力急剧下降。为了更直观地分析压力沿程变化规律，图4.4给出了不同工况下重力热管压力沿程变化曲线。从图中可以看出，在不同加热功率和充液率下，压力沿程变化曲线具有相似的趋势，但具体数值存在差异。随着加热功率的增加，蒸发段的压力明显升高。这是因为加热功率的提高使得蒸发过程更加剧烈，单位时间内产生的蒸汽量增加，从而导致蒸发段的压力上升。同时，由于蒸汽在管内流动时需要克服更大的阻力，从蒸发段到冷凝段的压力降也会增大。充液率对压力分布也有显著影响，当充液率增加时，管内液体量增多，液体对蒸汽的流动会产生更大的阻碍作用，使得蒸汽在流动过程中的压力损失增大，从而导致从蒸发段到冷凝段的压力降增大。此外，充液率过高还可能导致液体在管内积聚，形成液柱，进一步增加蒸汽的流动阻力，使压力分布更加不均匀。压力分布与流动状态密切相关。在蒸发段，较高的压力驱动蒸汽向冷凝段流动，压力差是蒸汽流动的主要动力来源。而在冷凝段，较低的压力有利于蒸汽的凝结，促进液体的回流。如果压力分布不合理，如蒸发段压力过低或冷凝段压力过高，可能会导致蒸汽流动不畅，甚至出现倒流现象，严重影响重力热管的传热性能。因此，在设计和运行重力热管时，需要充分考虑压力分布对流动状态的影响，通过合理调整加热功率、充液率等参数，优化压力分布，确保重力热管内流体的正常流动和高效传热。4.1.3流型分析通过数值模拟，成功识别出了不同工况下重力热管内的流型，并结合理论深入分析了流型转变的条件和影响因素。图4.5展示了在不同加热功率和充液率下重力热管内的流型。在低加热功率和较低充液率工况下（图4.5a），重力热管内主要呈现泡状流。此时，蒸汽以小气泡的形式分散在连续的液相中，气泡在浮力的作用下缓慢上升。这是因为在低加热功率下，液体蒸发速率较慢，产生的蒸汽量较少，蒸汽气泡在液体中分散分布，尚未形成连续的蒸汽相。而较低的充液率使得管内液体相对较少，蒸汽气泡有足够的空间在液体中分散存在，不易聚合形成较大的气泡或气柱。随着加热功率的增加或充液率的提高（图4.5b），流型逐渐转变为弹状流。在弹状流中，蒸汽气泡逐渐聚合长大，形成较大的气弹，气弹之间被液柱隔开。加热功率的增加使得液体蒸发速率加快，产生的蒸汽量增多，蒸汽气泡更容易聚合。而充液率的提高则使得管内液体量增加，液体对蒸汽气泡的约束作用增强，促使气泡聚合形成气弹。在弹状流状态下，气弹的快速上升会引起管内流体的剧烈波动，对传热性能产生重要影响。当加热功率进一步增大或充液率进一步提高时（图4.5c），流型可能转变为环状流。在环状流中，蒸汽在管中心形成连续的气相，而液体则在管壁上形成一层连续的液膜。此时，蒸汽的流速较高，对液膜产生较大的剪切力，使得液膜能够稳定地附着在管壁上。加热功率的进一步增大使得蒸汽量大幅增加，蒸汽的流速和流量都达到较高水平，从而能够将液体推向管壁形成液膜。而充液率的进一步提高则保证了有足够的液体形成液膜。环状流状态下，由于蒸汽和液膜之间的接触面积较大，传热效率相对较高，但如果液膜过薄或不稳定，可能会导致干涸现象的发生，影响热管的正常运行。流型转变的条件和影响因素较为复杂。加热功率是影响流型转变的重要因素之一，加热功率的增加会使液体蒸发速率加快，蒸汽的生成量和流速增大，从而促使流型从泡状流向弹状流、环状流转变。充液率也对流型转变有显著影响，充液率的变化会改变管内液体和蒸汽的相对含量和分布情况，进而影响流型。此外，管径、重力加速度等因素也会对流型产生影响。管径的减小会使蒸汽和液体的流动空间受限，更容易导致流型的转变；而重力加速度的变化则会影响蒸汽和液体的受力情况，对气液两相的相对运动和流型产生影响。在实际应用中，需要根据具体的工况条件，合理选择加热功率、充液率等参数，以确保重力热管内形成合适的流型，实现高效的传热性能。（a）泡状流（b）弹状流（c）环状流4.2传热特性分析4.2.1温度分布通过数值模拟得到了重力热管内的温度云图和不同位置的温度分布曲线，以深入分析温度场的形成和变化原因。图4.6展示了在加热功率为200W、充液率为50%工况下，重力热管的温度云图。从温度云图中可以清晰地看到，蒸发段的温度最高，这是因为蒸发段与高温热源接触，工作介质吸收热量，发生相变从液态转变为气态，在这个过程中，热量不断地从热源传递到工作介质中，使得蒸发段的温度升高。随着蒸汽从蒸发段流向冷凝段，温度逐渐降低。在绝热段，由于绝热段采用了良好的保温措施，热量散失极少，温度基本保持不变，蒸汽主要以等温流动的形式通过绝热段。而在冷凝段，蒸汽与低温冷源接触，放出热量重新凝结成液态，热量从工作介质传递到冷源中，导致冷凝段的温度降低。为了更直观地分析温度分布规律，图4.7给出了不同工况下重力热管轴向温度分布曲线。从图中可以看出，在不同加热功率和充液率下，温度分布曲线具有相似的趋势，但具体数值存在差异。随着加热功率的增加，蒸发段和冷凝段的温度都明显升高。这是因为加热功率的提高意味着更多的热量输入到重力热管中，在蒸发段，更多的热量使得工作介质蒸发加剧，蒸汽的温度升高；而在冷凝段，由于蒸汽携带的热量增加，需要释放更多的热量才能凝结成液体，因此冷凝段的温度也相应升高。充液率对温度分布也有一定的影响，当充液率增加时，管内液体量增多，液体的比热容较大，能够吸收更多的热量，从而在一定程度上降低了蒸发段和冷凝段的温度。但充液率过高时，可能会导致液体在管内积聚，影响蒸汽的正常流动和传热，使得温度分布不均匀，甚至出现局部过热或过冷的现象。温度场的形成和变化与重力热管的传热过程密切相关。在蒸发段，热量通过导热和对流的方式从热源传递到工作介质中，使工作介质蒸发，形成高温的蒸汽。蒸汽在管内流动过程中，通过对流换热将热量传递给管壁，管壁再通过导热将热量传递给周围环境。在冷凝段，蒸汽与冷源接触，通过对流换热和相变传热将热量释放给冷源，使蒸汽凝结成液体。因此，了解温度分布规律对于优化重力热管的传热性能具有重要意义。通过合理调整加热功率、充液率等参数，可以优化温度场分布，提高重力热管的传热效率和稳定性。4.2.2热流密度分布通过数值模拟获得了重力热管内的热流密度分布情况，以深入研究热流密度在不同区域的变化规律及影响因素。图4.8展示了在加热功率为250W、充液率为60%工况下，重力热管蒸发段和冷凝段的热流密度分布云图。从蒸发段热流密度分布云图（图4.8a）中可以看出，靠近热源的管壁处热流密度最高，这是因为热源的热量首先通过管壁传递给工作介质，在这个过程中，管壁与工作介质之间存在较大的温度差，根据傅里叶定律q=-\lambda\frac{\partialT}{\partialx}（其中q为热流密度，\lambda为导热系数，\frac{\partialT}{\partialx}为温度梯度），较大的温度差导致热流密度较高。随着离热源距离的增加，热流密度逐渐减小，这是因为热量在工作介质中传递时，会受到工作介质的导热阻力和对流换热的影响，使得热量逐渐分散，温度梯度减小，从而热流密度降低。在蒸发段的中心区域，由于蒸汽的生成和流动，热流密度相对较低，这是因为蒸汽的导热系数相对较小，且蒸汽的流动会带走一部分热量，使得中心区域的热量传递相对较弱。在冷凝段（图4.8b），热流密度分布情况与蒸发段相反。靠近冷源的管壁处热流密度最高，这是因为蒸汽在冷凝段与冷源接触，放出热量凝结成液体，热量从蒸汽传递到管壁再传递到冷源，在这个过程中，管壁与蒸汽之间存在较大的温度差，导致热流密度较高。随着离冷源距离的增加，热流密度逐渐减小，这是因为蒸汽在冷凝过程中，热量逐渐释放，温度逐渐降低，温度梯度减小，从而热流密度降低。在冷凝段的中心区域，由于蒸汽的逐渐凝结，热流密度也相对较低。为了更直观地分析热流密度变化规律，图4.9给出了不同工况下重力热管蒸发段和冷凝段热流密度沿程变化曲线。从图中可以看出，在不同加热功率和充液率下，热流密度沿程变化曲线具有相似的趋势，但具体数值存在差异。随着加热功率的增加，蒸发段和冷凝段的热流密度都明显增大。这是因为加热功率的提高意味着更多的热量输入到重力热管中，在蒸发段，更多的热量需要通过工作介质传递，导致热流密度增大；在冷凝段，由于蒸汽携带的热量增加，需要释放更多的热量，使得热流密度也相应增大。充液率对热流密度的影响较为复杂，当充液率增加时，管内液体量增多，液体的导热性能和对流换热能力会发生变化，从而影响热流密度。在一定范围内，充液率的增加可能会使热流密度增大，这是因为更多的液体可以更好地传递热量；但当充液率过高时，液体可能会在管内积聚，影响蒸汽的正常流动和传热，导致热流密度下降。热流密度在不同区域的变化规律受到多种因素的影响，除了加热功率和充液率外，还与工作介质的热物理性质、管径、管壁材料等因素有关。工作介质的导热系数、汽化潜热等热物理性质会直接影响热量的传递效率和热流密度的大小。管径的大小会影响流体的流动状态和传热面积，进而影响热流密度。管壁材料的导热系数也会对热流密度产生影响，导热系数较高的管壁材料能够更有效地传递热量，降低热阻，提高热流密度。因此，在设计和优化重力热管时，需要综合考虑这些因素，以获得最佳的热流密度分布和传热性能。（a）蒸发段（b）冷凝段4.2.3传热效率分析为了准确评估重力热管的传热性能，定义传热效率\eta为：\eta=\frac{Q_{actual}}{Q_{input}}\times100\%其中，Q_{actual}为重力热管实际传递的热量，可通过冷凝段释放的热量来计算，即Q_{actual}=m\cdotr（m为冷凝段凝结的工质质量，r为工质的汽化潜热）；Q_{input}为输入到蒸发段的热量，可根据加热功率P和时间t计算，即Q_{input}=P\cdott。通过数值模拟分析了不同参数对传热效率的影响，以得出提高传热效率的途径。图4.10展示了不同加热功率和充液率下重力热管的传热效率变化曲线。从图中可以看出，随着加热功率的增加，传热效率呈现先增加后减小的趋势。在较低加热功率下，随着加热功率的提高，蒸发段的温度升高，工作介质蒸发加剧，蒸汽的生成量和流速增大，使得热量传递更加迅速，传热效率提高。然而，当加热功率超过一定值后，继续增加加热功率，传热效率反而下降。这是因为过高的加热五、实验研究与验证5.1实验装置设计与搭建为了对重力热管的流动与传热特性进行实验研究，并验证数值模拟结果的准确性，精心设计并搭建了一套实验装置。该装置主要包括重力热管本体、加热系统、冷却系统和测量系统。在重力热管制作过程中，选用外径为25mm、壁厚为2mm的紫铜管作为管壳材料，紫铜管具有良好的导热性能，能够有效减少管壁的热阻，提高热管的传热效率。管壳的长度根据实验需求确定为1500mm，其中蒸发段长度为500mm，绝热段长度为500mm，冷凝段长度为500mm。端盖采用相同材质的紫铜加工而成，通过氩弧焊将端盖与管壳进行密封焊接，确保热管内部的密封性，防止工作介质泄漏。在完成管壳和端盖的焊接后，对热管进行严格的清洗和除油处理，以去除表面的油污和杂质，保证热管内部的清洁度。然后，将热管抽真空至10^-3Pa以下，以排除管内的空气和其他不凝性气体，为后续充注工作介质创造良好的真空环境。充注工作介质时，选择去离子水作为工作介质，根据实验设计的充液率，使用高精度的微量注射泵将适量的去离子水注入热管内。充液完成后，再次对热管进行密封处理，确保工作介质不会泄漏。加热系统采用电加热方式，使用加热功率可调节的电加热器对重力热管的蒸发段进行加热。电加热器由镍铬合金电阻丝和绝缘耐高温材料制成，电阻丝均匀缠绕在蒸发段管壳外表面，通过调节输入电压来控制加热功率。为了减少热量散失，在蒸发段外表面包裹一层厚度为50mm的硅酸铝纤维保温棉，保温棉的导热系数低，能够有效地阻止热量向周围环境传递，提高加热效率。冷却系统采用循环水冷却方式，通过冷却水箱、循环水泵和冷凝器组成的循环回路对重力热管的冷凝段进行冷却。冷却水箱中储存有一定量的冷却水，循环水泵将冷却水从水箱中抽出，送入冷凝器中进行冷却，然后再将冷却后的水输送到重力热管的冷凝段，吸收热管内蒸汽放出的热量后，返回冷却水箱。在冷凝器中，通过与外部空气进行热交换，将冷却水中的热量散发到大气中，从而保证冷却水的温度始终保持在较低水平。为了精确控制冷却水的流量和温度，在循环水管路上安装了转子流量计和温度传感器，可实时监测和调节冷却水的流量和温度。测量系统主要用于测量重力热管的温度、压力等参数。在重力热管的蒸发段、绝热段和冷凝段的管壳表面均匀布置了10个K型热电偶，用于测量不同位置的管壁温度，热电偶的测量精度为±0.5℃。通过数据采集系统将热电偶测量的温度数据实时采集并传输到计算机中进行记录和分析。在蒸发段和冷凝段的进出口位置分别安装了压力传感器，用于测量蒸汽和液体的压力，压力传感器的测量精度为±0.01MPa。此外，还在加热系统和冷却系统的管路上安装了功率表和流量计，用于测量加热功率和冷却水流量。通过以上设计和搭建，构建了一套完整的重力热管实验装置，能够满足对重力热管流动与传热特性实验研究的需求，为后续的实验测试和结果分析提供了可靠的硬件基础。5.2实验方案与步骤确定实验工况和变量控制方法，以全面研究重力热管在不同条件下的流动与传热特性。实验工况主要包括不同的加热功率、充液率和冷却水温。加热功率设置为50W、100W、150W、200W、250W五个水平，以研究加热功率对热管性能的影响。充液率设置为20%、30%、40%、50%、60%五个水平，通过控制充入热管内的工作介质的量来改变充液率，探究充液率对热管性能的作用。冷却水温设置为20℃、25℃、30℃三个水平，通过调节冷却系统中冷却水的温度，分析冷却水温对热管传热性能的影响。在实验过程中，严格控制其他变量保持不变。重力热管的管径、长度、管壳材料等几何参数在整个实验过程中固定不变。实验环境的大气压力和温度也尽量保持稳定，以减少环境因素对实验结果的干扰。每次实验前，确保实验装置处于正常工作状态，各测量仪器经过校准，测量数据准确可靠。实验操作步骤如下：首先，检查实验装置的各个部分，确保加热系统、冷却系统、测量系统等连接正确，无泄漏现象。开启冷却系统，启动循环水泵，使冷却水在冷却回路中循环流动，并调节冷却水温至设定值。开启加热系统，将加热功率设置为初始值，开始对重力热管的蒸发段进行加热。在加热过程中，密切观察测量系统的数据变化，记录热管各部位的温度、压力以及加热功率、冷却水流量等参数。当热管达到稳定工作状态后，即热管各部位的温度和压力在一段时间内基本保持不变时，每隔30s记录一次测量数据，连续记录10组数据，以保证数据的准确性和可靠性。按照上述步骤，依次改变加热功率、充液率和冷却水温，进行不同工况下的实验测试。数据采集流程为：在实验过程中，K型热电偶将测量的温度信号转换为电信号，通过数据采集线传输到数据采集仪中。压力传感器将测量的压力信号也传输到数据采集仪中。数据采集仪对采集到的温度和压力信号进行放大、滤波等处理后，通过RS485通信接口将数据传输到计算机中。计算机通过专用的数据采集软件对数据进行实时显示、记录和存储。同时，加热系统的功率表和冷却系统的流量计也将测量的数据传输到计算机中，与温度和压力数据一起进行综合分析。通过以上实验方案和步骤，能够系统地研究不同因素对重力热管流动与传热特性的影响，为深入了解重力热管的工作原理和性能提供实验依据。5.3实验结果与数值模拟对比验证将实验测量的温度、压力等数据与数值模拟结果进行对比，以评估模拟的准确性。图5.1展示了在加热功率为150W、充液率为40%、冷却水温为25℃工况下，重力热管轴向温度分布的实验结果与数值模拟结果对比。从图中可以看出，实验测量的温度分布与数值模拟结果总体趋势一致，在蒸发段，温度迅速升高，达到工作介质的沸点温度后基本保持不变；在绝热段，温度略有下降；在冷凝段，温度逐渐降低。然而，仔细观察发现，实验结果与数值模拟结果在某些位置存在一定的差异。在蒸发段靠近入口处，实验测量的温度略高于数值模拟结果，这可