超临界二氧化碳燃煤火电冷却壁多物理场耦合：机理、特性与优化策略

超临界二氧化碳燃煤火电冷却壁多物理场耦合：机理、特性与优化策略一、引言1.1研究背景与意义1.1.1研究背景在全球能源需求持续增长以及环境保护意识不断增强的大背景下，能源领域的技术革新成为了至关重要的议题。燃煤发电作为目前全球范围内主要的发电方式之一，在满足能源需求方面发挥着关键作用。然而，传统的燃煤发电技术面临着诸多严峻挑战，如能源利用效率偏低、环境污染问题严重等。据统计，传统燃煤发电技术的能源转换效率通常在30%-40%之间，大量的能源在发电过程中被浪费。同时，其排放的大量污染物，包括二氧化硫、氮氧化物和颗粒物等，对大气环境造成了严重污染，加剧了酸雨、雾霾等环境问题。超临界二氧化碳燃煤发电技术作为一种新兴的高效发电技术，应运而生并逐渐成为研究热点。该技术通过将二氧化碳提升至超临界状态（温度高于31.1℃，压力高于7.38MPa），利用其独特的物理性质，实现了煤炭的高效燃烧和能源的充分利用。与传统的燃煤发电技术相比，超临界二氧化碳燃煤发电技术具有显著优势。其热效率得到了大幅提高，能够达到45%-50%左右，有效减少了能源的浪费。在污染物排放方面，该技术能够更有效地降低二氧化硫、氮氧化物和颗粒物等污染物的排放，对环境保护起到了积极的推动作用。超临界二氧化碳发电系统还具有结构简单、运行灵活、设备体积小等优点，对于提高电力系统的稳定性和可靠性具有重要意义。在超临界二氧化碳燃煤发电系统中，冷却壁作为关键部件之一，承担着保护锅炉炉壁、维持系统稳定运行的重要职责。冷却壁的工作环境极为复杂，处于高温、高压以及强热流的恶劣条件之下。在这样的环境中，冷却壁内部会同时发生多种物理现象，包括热传导、对流换热、流体流动以及化学反应等，这些物理现象相互作用、相互影响，形成了复杂的多物理场耦合效应。这种多物理场耦合效应会对冷却壁的性能和寿命产生显著影响。如果冷却壁的温度分布不均匀，可能会导致热应力集中，从而使冷却壁出现裂纹、变形等损坏情况，严重影响系统的安全运行。冷却壁内的流体流动状态也会影响到换热效率和系统的能耗。因此，深入研究超临界二氧化碳燃煤发电系统冷却壁的多物理场耦合特性，对于保障系统的安全稳定运行、提高发电效率具有至关重要的意义。1.1.2研究意义本研究对于提高超临界二氧化碳燃煤发电系统的发电效率具有重要意义。通过深入研究冷却壁的多物理场耦合特性，可以优化冷却壁的结构设计和运行参数，提高其换热效率，从而减少能源在冷却过程中的损失，使更多的热能能够转化为电能，进而提高整个发电系统的能源转换效率。例如，通过对冷却壁内流体流动和换热过程的精确模拟和分析，可以确定最佳的冷却介质流量和流速，以实现最有效的热量传递，提高发电效率。冷却壁作为超临界二氧化碳燃煤发电系统的关键部件，其性能和寿命直接关系到整个系统的安全稳定运行。通过研究多物理场耦合对冷却壁的影响，可以揭示冷却壁损坏的机理，从而采取针对性的措施来提高冷却壁的可靠性和寿命。例如，通过优化冷却壁的材料选择和结构设计，降低热应力集中，减少裂纹和变形的发生，确保系统能够长期稳定运行，降低设备维护成本和停机时间，提高电力供应的可靠性。超临界二氧化碳燃煤发电技术作为一种新型的发电技术，其发展和应用对于推动能源行业的技术进步具有重要作用。本研究关于冷却壁多物理场耦合的成果，不仅可以为超临界二氧化碳燃煤发电系统的设计、优化和运行提供理论支持和技术指导，还可以为相关领域的研究提供参考和借鉴，促进整个能源行业在高效、清洁发电技术方面的发展，推动能源结构的优化升级，助力实现可持续发展的能源目标。1.2国内外研究现状在国外，对于超临界二氧化碳燃煤发电系统的研究开展较早。美国、欧盟等国家和地区投入了大量资源进行相关技术的研发与探索。在冷却壁多物理场耦合研究方面，美国的一些科研团队通过数值模拟手段，对超临界二氧化碳在冷却壁内的流动与传热特性进行了深入分析。他们利用先进的计算流体力学（CFD）软件，建立了详细的物理模型，模拟了不同工况下冷却壁内的温度分布、速度场以及压力变化等。研究发现，超临界二氧化碳的特殊物性，如在近临界区比热容的急剧变化，会导致冷却壁内的传热出现复杂的非线性行为，对冷却效果产生显著影响。欧盟的相关研究则侧重于实验研究，通过搭建实验平台，对不同结构和材料的冷却壁进行实验测试，获取了大量的实验数据。这些实验研究为数值模拟提供了验证依据，同时也揭示了冷却壁在实际运行过程中面临的问题，如高温腐蚀、热疲劳等，这些问题与多物理场耦合效应密切相关。国内在超临界二氧化碳燃煤发电技术领域的研究虽然起步相对较晚，但近年来发展迅速。众多高校和科研机构积极参与相关研究，取得了一系列重要成果。在冷却壁多物理场耦合研究方面，国内学者从理论分析、数值模拟和实验研究多个角度展开工作。一些高校通过建立理论模型，对冷却壁内的热传导、对流换热以及化学反应等物理过程进行了理论推导和分析，为深入理解多物理场耦合机理提供了理论基础。在数值模拟方面，国内研究团队利用自主开发的软件和商业CFD软件，对冷却壁的多物理场耦合进行了全面的模拟研究。通过模拟，分析了不同运行参数（如压力、温度、流量等）对多物理场耦合特性的影响，为冷却壁的优化设计提供了理论指导。在实验研究方面，国内也搭建了多个实验平台，对冷却壁的传热、流动等性能进行实验测试。例如，通过实验研究了不同冷却介质（如超临界二氧化碳、水等）在冷却壁内的传热性能差异，以及不同结构的冷却壁对多物理场耦合的影响。尽管国内外在超临界二氧化碳燃煤发电系统冷却壁多物理场耦合方面取得了一定的研究成果，但仍存在一些不足之处。目前的研究大多集中在单一物理场的分析或简单的多物理场耦合模拟，对于复杂工况下多物理场之间的强耦合作用及其相互影响机制的研究还不够深入。例如，在实际运行过程中，冷却壁不仅受到高温烟气的热冲击，还可能受到机械振动、化学腐蚀等多种因素的影响，这些因素与多物理场耦合相互交织，使得问题更加复杂，但目前对此类复杂情况的研究较少。实验研究方面，由于实验条件的限制，难以全面模拟实际运行中的各种工况，导致实验数据的代表性和可靠性存在一定局限性。而且，不同研究团队的实验条件和方法存在差异，使得实验结果之间的可比性较差，不利于研究成果的整合和推广。在冷却壁的材料选择和结构设计方面，虽然已经开展了一些研究，但如何综合考虑多物理场耦合效应，实现冷却壁的优化设计，以提高其性能和寿命，仍然是一个亟待解决的问题。1.3研究目标与内容1.3.1研究目标本研究旨在深入探究超临界二氧化碳燃煤发电系统冷却壁的多物理场耦合机理，全面揭示热传导、对流换热、流体流动以及化学反应等物理过程之间的相互作用机制，明确各物理场之间的耦合关系和影响规律。通过对多物理场耦合特性的深入分析，建立准确可靠的多物理场耦合数学模型和数值计算方法，为冷却壁的性能预测和优化设计提供坚实的理论基础和有效的技术手段。基于研究成果，提出一套科学合理、切实可行的冷却壁结构优化策略和运行参数优化方案，以显著提高冷却壁的换热效率，降低其温度分布的不均匀性，有效减少热应力集中现象，从而延长冷却壁的使用寿命，保障超临界二氧化碳燃煤发电系统的安全稳定运行，提高整个发电系统的能源利用效率。1.3.2研究内容对超临界二氧化碳的特殊物理性质进行深入研究，包括其在超临界状态下的密度、黏度、比热容、导热系数等随温度和压力的变化规律。通过实验研究和理论分析相结合的方法，获取准确的物性数据，为后续的多物理场耦合分析提供可靠的基础数据。例如，利用高精度的实验设备，测量不同温度和压力条件下超临界二氧化碳的物性参数，并与现有理论模型进行对比分析，验证和完善理论模型。深入研究冷却壁内的热传导过程，建立热传导数学模型，分析不同材料和结构的冷却壁在高温、高压环境下的热传导特性。考虑冷却壁材料的热物性随温度的变化，以及热辐射对热传导的影响，通过数值模拟和实验验证，揭示热传导过程中的关键影响因素和规律。例如，采用有限元方法对冷却壁的热传导进行数值模拟，分析不同材料和结构参数对热传导的影响，并通过实验测量冷却壁的温度分布，验证数值模拟结果的准确性。对冷却壁内超临界二氧化碳的对流换热过程进行全面研究，分析对流换热系数的影响因素，建立对流换热数学模型。考虑超临界二氧化碳在近临界区的特殊物性对对流换热的影响，以及冷却壁内流道结构对流体流动和换热的影响，通过实验研究和数值模拟，优化对流换热过程，提高换热效率。例如，通过实验研究不同工况下超临界二氧化碳在冷却壁内的对流换热特性，分析流速、温度、压力等因素对对流换热系数的影响，并利用数值模拟优化流道结构，提高对流换热效率。详细研究冷却壁内超临界二氧化碳的流体流动特性，包括流速分布、压力分布、流量分配等。建立流体流动数学模型，考虑流体的可压缩性、黏性以及冷却壁内流道的几何形状和粗糙度等因素，通过数值模拟和实验研究，分析流体流动过程中的压力损失和流动稳定性，优化流体流动状态，降低流动阻力。例如，采用计算流体力学（CFD）方法对冷却壁内的流体流动进行数值模拟，分析不同流道结构和运行参数对流体流动的影响，并通过实验测量流速和压力分布，验证数值模拟结果的可靠性。综合考虑热传导、对流换热和流体流动等物理过程，建立冷却壁多物理场耦合数学模型。采用合适的数值计算方法对耦合模型进行求解，分析多物理场之间的相互作用机制和耦合效应，研究不同工况下冷却壁的温度分布、热应力分布以及换热和流动性能。例如，将热传导、对流换热和流体流动模型进行耦合，利用数值计算方法求解耦合模型，分析多物理场之间的相互影响，并通过实验研究验证耦合模型的准确性。基于多物理场耦合研究成果，对冷却壁的结构进行优化设计。通过改变冷却壁的形状、尺寸、材料以及内部流道结构等参数，利用数值模拟和优化算法，以提高换热效率、降低热应力为目标，确定冷却壁的最优结构参数。例如，采用响应面法等优化算法，结合数值模拟，对冷却壁的结构参数进行优化，得到最优的结构设计方案。研究冷却壁在不同运行参数下的性能变化规律，包括温度、压力、流量等参数对冷却壁换热、流动和热应力的影响。通过数值模拟和实验研究，确定冷却壁的最佳运行参数范围，为超临界二氧化碳燃煤发电系统的安全稳定运行提供操作依据。例如，通过实验研究不同运行参数下冷却壁的性能，分析参数变化对冷却壁性能的影响，并利用数值模拟预测不同工况下冷却壁的性能，确定最佳运行参数范围。1.4研究方法与技术路线1.4.1研究方法数值模拟：运用计算流体力学（CFD）软件，如ANSYSFluent、COMSOLMultiphysics等，对超临界二氧化碳在冷却壁内的流动、传热以及多物理场耦合过程进行数值模拟。通过建立精确的几何模型和合理的物理模型，设置与实际工况相符的边界条件和初始条件，模拟不同运行参数下冷却壁内的温度分布、速度场、压力分布以及热应力分布等，深入分析多物理场之间的相互作用机制和耦合效应。例如，利用Fluent软件中的RNGk-ε湍流模型来模拟超临界二氧化碳的湍流流动，采用有限体积法对控制方程进行离散求解，得到冷却壁内详细的流场和温度场信息。实验研究：搭建超临界二氧化碳燃煤发电系统冷却壁实验平台，开展实验研究。实验平台主要包括超临界二氧化碳供应系统、加热系统、冷却壁试验段、测量系统等。通过实验测量不同工况下冷却壁内超临界二氧化碳的流速、温度、压力等参数，以及冷却壁的壁面温度、热流密度等，获取实际运行数据，用于验证数值模拟结果的准确性，并为理论分析提供实验依据。例如，使用高精度的温度传感器和压力传感器测量冷却壁不同位置的温度和压力，利用粒子图像测速（PIV）技术测量超临界二氧化碳的流速分布。理论分析：基于传热学、流体力学、热力学等基本理论，对冷却壁内的热传导、对流换热、流体流动以及化学反应等物理过程进行理论分析。建立相应的数学模型，推导相关的控制方程和边界条件，通过理论计算和分析，揭示各物理过程的内在规律和多物理场耦合的基本原理。例如，根据傅里叶定律建立冷却壁的热传导方程，依据牛顿冷却定律建立对流换热方程，运用纳维-斯托克斯方程描述超临界二氧化碳的流体流动，从理论上分析各物理量之间的关系和影响因素。1.4.2技术路线第一步：文献调研与理论基础研究：广泛查阅国内外关于超临界二氧化碳燃煤发电系统冷却壁多物理场耦合的相关文献资料，了解该领域的研究现状、发展趋势以及存在的问题。深入学习超临界二氧化碳的物理性质、传热学、流体力学、热力学等基础理论知识，为后续的研究工作奠定坚实的理论基础。第二步：物理模型建立：根据冷却壁的实际结构和运行工况，建立冷却壁的几何模型，并对其进行合理的简化和假设。综合考虑热传导、对流换热、流体流动以及化学反应等物理过程，建立冷却壁多物理场耦合的数学模型，确定模型的控制方程、边界条件和初始条件。第三步：数值模拟计算：将建立好的数学模型导入到CFD软件中进行数值模拟计算，选择合适的数值计算方法和求解器，对控制方程进行离散求解。通过模拟不同运行参数（如压力、温度、流量等）下冷却壁内的多物理场分布，分析各物理场之间的相互作用和耦合效应，得到冷却壁的温度分布、热应力分布、换热性能和流动性能等模拟结果。第四步：实验平台搭建与实验研究：按照实验方案搭建超临界二氧化碳燃煤发电系统冷却壁实验平台，对实验设备进行调试和校准，确保实验数据的准确性和可靠性。在不同工况下开展实验研究，测量冷却壁内超临界二氧化碳的流速、温度、压力等参数，以及冷却壁的壁面温度、热流密度等数据。对实验数据进行整理和分析，与数值模拟结果进行对比验证，评估数值模拟模型的准确性和可靠性。第五步：结果分析与讨论：综合数值模拟结果和实验数据，对冷却壁多物理场耦合特性进行深入分析和讨论。研究不同运行参数对冷却壁性能的影响规律，揭示多物理场耦合的内在机制和关键影响因素。根据分析结果，提出优化冷却壁性能的建议和措施。第六步：冷却壁结构优化与运行参数优化：基于多物理场耦合研究成果，以提高换热效率、降低热应力为目标，对冷却壁的结构进行优化设计。利用数值模拟和优化算法，如遗传算法、粒子群优化算法等，对冷却壁的形状、尺寸、材料以及内部流道结构等参数进行优化，确定冷却壁的最优结构参数。同时，研究冷却壁在不同运行参数下的性能变化规律，通过数值模拟和实验研究，确定冷却壁的最佳运行参数范围，为超临界二氧化碳燃煤发电系统的安全稳定运行提供操作依据。第七步：研究成果总结与展望：对整个研究工作进行总结，归纳研究成果，撰写研究报告和学术论文。总结研究过程中存在的问题和不足之处，提出未来进一步研究的方向和建议。二、超临界二氧化碳燃煤火力发电系统及冷却壁概述2.1超临界二氧化碳燃煤火力发电系统原理与特点2.1.1系统工作原理超临界二氧化碳燃煤发电系统以超临界二氧化碳作为工质，遵循布雷顿循环实现能量的高效转化。系统主要由燃煤锅炉、主压缩机、主加热器、高压透平、再热器、低压透平、回热器、冷却器等部件组成。在整个能量转换过程中，煤炭在锅炉中燃烧释放出大量的热能，产生高温高压的烟气。这部分热能首先通过中间传热回路传递给超临界二氧化碳工质，使工质的温度和压力升高。中间传热回路通常采用液态金属或熔盐作为中间介质，其运行于常压工况，有利于换热面的设计，降低了管道材料要求，同时增大了锅炉的热负荷调节范围，提高了机组的灵活性。被加热的超临界二氧化碳工质随后进入主压缩机进行绝热压缩，压力进一步提升。压缩后的工质进入主加热器，在其中吸收更多的热量，达到更高的温度和压力。高温高压的超临界二氧化碳工质接着进入高压透平，推动透平的叶片高速旋转，将工质的热能转化为机械能，驱动发电机发电。从高压透平排出的工质压力和温度有所降低，为了进一步提高能量利用效率，这些工质会进入再热器，再次吸收热量，温度升高后进入低压透平继续膨胀做功，进一步驱动发电机发电。完成做功后的超临界二氧化碳工质，压力和温度大幅下降，进入回热器。在回热器中，工质将自身携带的部分热量传递给即将进入主加热器的低温工质，实现热量的回收利用，提高系统的热效率。经过回热器后的工质进入冷却器，将剩余的热量传递给冷却水或其他冷却介质，使工质温度降低，回到初始状态。冷却后的超临界二氧化碳工质再次被主压缩机吸入，开始下一个循环。整个系统通过这样的闭式循环，不断地将煤炭燃烧产生的热能转化为电能。2.1.2系统特点与优势与传统的燃煤发电系统相比，超临界二氧化碳燃煤发电系统在效率方面具有显著优势。传统的蒸汽朗肯循环发电系统存在着工质相态变化带来的能量损失，而超临界二氧化碳发电系统以超临界二氧化碳为工质，避免了这种相态变化，从而减少了压缩功的消耗。超临界二氧化碳工质的密度接近液体，大于气体2个数量级，这使得其传热效率高，做功能力强。相关研究表明，在相同的热源温度下，超临界二氧化碳发电系统的效率可比传统的水蒸气朗肯循环发电系统提高3%-5%以上，能够更有效地利用煤炭燃烧释放的热能，减少能源的浪费。超临界二氧化碳发电系统具有良好的灵活性。该系统的启动和停止时间相对较短，能够更快地响应负载变化。这一特点使其在应对电网波动和应急情况时表现出色，能够更好地满足电力系统对灵活性的要求。由于超临界二氧化碳的密度较低，整个系统的体积更小，占地面积小，这使得它适用于各种空间有限的场所，如离岸风电、城市内部的热电联产等，拓宽了发电系统的应用场景。超临界二氧化碳无毒、不燃、稳定，对臭氧层无破坏，且廉价易得。在发电过程中，超临界二氧化碳工质不会产生二氧化硫、氮氧化物等污染物，与传统的燃煤发电相比，大大减少了对环境的污染。超临界二氧化碳发电系统还可以与碳捕获和封存技术相结合，进一步降低二氧化碳的排放，有助于实现全球的碳减排目标，对环境保护具有重要意义。在设备成本方面，超临界二氧化碳的密度大、压缩性好，使得系统中的设备结构紧凑、体积小。相同功率输出下，超临界二氧化碳发电系统所需的管道、阀门、透平等设备的尺寸相对较小，从而降低了设备的制造和安装成本。超临界二氧化碳的工作温度和压力相对较低，对设备材料的要求也相对较低，进一步降低了设备成本。而且，由于系统的高效率和灵活性，发电成本也相应降低，提高了系统的经济性。2.2冷却壁在系统中的作用与结构形式2.2.1冷却壁的作用冷却壁在超临界二氧化碳燃煤发电系统中起着至关重要的作用，对系统的安全稳定运行和高效能量转换有着不可或缺的影响。冷却壁是保护锅炉炉壁的关键屏障。在燃煤锅炉中，煤炭燃烧会释放出大量的热量，使炉膛内的温度极高，一般可达1000℃以上。如此高温环境对锅炉炉壁材料的耐高温性能提出了极高要求，如果炉壁直接暴露在高温烟气中，会导致炉壁材料的强度迅速下降，发生变形、损坏，甚至引发安全事故。冷却壁通过内部流动的超临界二氧化碳工质吸收炉膛内的高温热量，有效降低炉壁温度，使炉壁温度保持在材料允许的安全范围内。相关研究表明，安装冷却壁后，炉壁温度可降低数百度，大大提高了炉壁的安全性和可靠性，保障了锅炉的长期稳定运行。冷却壁在热量传递和发电效率提升方面发挥着核心作用。超临界二氧化碳工质具有良好的传热性能，冷却壁内的超临界二氧化碳通过对流换热的方式，将炉膛内的热量迅速传递出去。这些热量一部分被用于加热其他超临界二氧化碳工质，使其达到更高的温度和压力，为透平做功提供更多的能量；另一部分热量则被回收利用，通过回热器等设备预热即将进入主加热器的低温工质，提高了系统的热效率。研究数据显示，冷却壁的高效换热能够使系统的热效率提高5%-10%左右，有效减少了能源的浪费，提高了煤炭的利用效率。冷却壁还对维持系统的压力平衡和流体稳定流动起着重要作用。在超临界二氧化碳燃煤发电系统中，系统内的压力和流体流动状态对发电效率和设备寿命有着显著影响。冷却壁内的超临界二氧化碳工质在吸收热量的过程中，其压力和流速会发生变化，冷却壁的结构和布置需要保证工质在管内能够稳定流动，避免出现压力波动、流量分配不均等问题。合理设计的冷却壁能够使超临界二氧化碳工质在管内均匀流动，维持系统的压力平衡，确保整个发电系统的稳定运行。例如，通过优化冷却壁的管径、管长以及内部流道结构，可以有效减少压力损失，提高流体的流动稳定性，保障系统的高效运行。2.2.2常见结构形式在超临界二氧化碳燃煤发电系统中，冷却壁存在多种常见的结构形式，每种结构形式都有其独特的特点和适用场景。螺旋管圈冷却壁是一种较为常见的结构形式。这种冷却壁的管子从炉膛底部以一定的倾角沿着炉膛四周盘旋上升，一般盘旋1.5-2.0圈后，上部改为垂直上升管子。螺旋管圈冷却壁的主要优点在于其能够很好地适应变压运行的要求。由于管子根数相对较少，且采用了较高的质量流速，在所有的负荷范围内均能保证管壁温度不超过允许温度，即使在压力较低时，水动力也十分稳定。螺旋管圈绕炉膛周界盘旋上升，使得各根管子受热较为均匀，水冷壁出口温度偏差小。而且，它可以采用较大管径和壁厚的水冷壁管，增加了水冷壁的刚性。不过，螺旋管圈冷却壁也存在一些缺点。由于管内质量流速高，螺旋管圈展开长度几乎为垂直管圈长度的二倍，导致水冷壁阻力较大，增加了系统的能耗。其制造工艺比较复杂，制造工作量大，每绕一圈就有四个弯头，组装率低，管圈的支承结构和刚性梁结构也较为复杂，制造和安装周期长，后期的维护和维修也相对复杂。垂直管圈冷却壁也是一种广泛应用的结构形式，又可细分为一次上升式和多次上升-下降式。一次上升垂直管圈的所有管屏都是并联的，工质从省煤器引入炉底进口集箱，在管屏中一次向上流动至炉顶出口集箱。多次上升-下降管圈则是工质从炉底进入几片管屏，向上流动到炉顶后，经过下降管引到炉底，再在另外几片管屏中向上流动，根据不同情况可有几次上升下降。多次垂直上升-下降管圈工质具有较高的质量流速，但由于相邻管屏间工质温度不一样，会引起相邻管屏外侧两根相邻管子之间壁温差大，因此只适用于定压运行的锅炉。对于变压运行的超临界锅炉，通常采用一次上升的垂直管圈水冷壁。为了得到较高的质量流速，一般要求锅炉容量较大并采用较细的管径。同时，为了抑制亚临界压力下炉膛水冷壁传热恶化、强化管内侧换热和确保水冷壁管工作的安全性，在热负荷较高的部位采用内螺纹管。一次上升垂直管圈水冷壁与螺旋管圈相比，质量流速较低，管子总长也较短，水冷壁阻力较小，降低了给水泵的耗电量。其水冷壁本身结构、支撑和刚性梁结构简单，安装、维护和检修比较容易。然而，它也存在一些局限性，如受到机组容量的限制，目前的最小管径为28mm，由于管径的限制，对于容量较小的锅炉，无法保证必要的质量流速，一般认为，一次上升垂直管圈水冷壁适用的锅炉最小容量为500-600MW。而且，水冷壁管径较小并采用内螺纹管，管子的制造精度和价格较高，还需要在水冷壁入口装设节流圈，增加了水冷壁下集箱结构的复杂性，水冷壁出口的温度偏差也比螺旋管圈大，启动和低负荷时为了保持必要的质量流速，需装设再循环泵。还有一种是膜式冷却壁。膜式冷却壁由管子和扁钢焊接而成，形成一个整体的壁面，管子之间通过扁钢连接，使得冷却壁具有良好的密封性和整体性。这种结构形式能够有效减少炉膛内的漏风现象，提高锅炉的热效率。膜式冷却壁的刚性较好，能够承受一定的热应力和机械应力，在高温、高压的工作环境下具有较高的可靠性。由于其结构紧凑，占用空间小，有利于锅炉的小型化设计。但是，膜式冷却壁的制造工艺要求较高，焊接质量直接影响到冷却壁的性能和寿命。在运行过程中，如果出现局部过热或损坏，修复难度较大。2.3冷却壁工作条件与多物理场环境在超临界二氧化碳燃煤发电系统中，冷却壁运行时炉膛内的温度通常高达1000-1500℃，这使得冷却壁直接承受着极高的温度。例如，在某600MW超临界二氧化碳燃煤发电项目中，实际测量炉膛中心区域的温度可达1300℃左右，冷却壁表面直接受到高温烟气的辐射和冲刷，处于这样的高温环境中，冷却壁材料的性能会受到严重考验，其强度、硬度等力学性能会随着温度的升高而下降。在高温作用下，材料内部的晶体结构可能发生变化，导致材料的蠕变和疲劳性能变差，增加了冷却壁发生变形和损坏的风险。超临界二氧化碳发电系统的运行压力一般在20-30MPa之间。以某超临界二氧化碳燃煤发电示范工程为例，其运行压力稳定在25MPa左右。在如此高的压力下，冷却壁承受着巨大的压力载荷。高压会使冷却壁管受到径向和轴向的压力作用，可能导致管壁发生屈服、破裂等失效形式。高压还会影响冷却壁内超临界二氧化碳的物理性质，如密度、黏度等，进而影响其流动和传热特性。在燃煤锅炉中，煤炭燃烧释放出大量的热量，使得冷却壁受到非常强的热流作用。热流密度通常在100-500kW/m²之间。在一些大型燃煤锅炉中，燃烧器附近的冷却壁热流密度甚至可达800kW/m²以上。强热流会使冷却壁表面的温度迅速升高，导致冷却壁内部产生较大的温度梯度。例如，当热流密度为300kW/m²时，冷却壁表面温度与内部温度的差值可能达到数百度，这种温度梯度会引发热应力，对冷却壁的结构完整性造成威胁。强热流还会加剧冷却壁的热疲劳现象，使其在长期运行过程中更容易出现裂纹和损坏。在冷却壁内部，超临界二氧化碳的流动伴随着复杂的对流换热过程。超临界二氧化碳的特殊物性，如在近临界区比热容的急剧变化，使得对流换热过程呈现出高度的非线性和复杂性。当超临界二氧化碳的温度和压力接近临界状态时，其比热容会发生剧烈变化，导致对流换热系数大幅波动。在某超临界二氧化碳冷却壁实验中，当工质温度接近临界温度时，对流换热系数在短时间内变化了数倍，这使得冷却壁的换热性能难以稳定控制，增加了冷却壁局部过热的风险。冷却壁在高温环境下，其材料会与周围的气体、杂质等发生化学反应。例如，在燃煤锅炉中，冷却壁表面可能会与烟气中的二氧化硫、氧气等发生氧化、硫化等反应，形成腐蚀产物，导致冷却壁材料的腐蚀和损坏。在高温和强热流的作用下，冷却壁表面的氧化膜可能会被破坏，使得内部材料更容易受到腐蚀介质的侵蚀。而且，化学反应还可能改变冷却壁表面的物理性质，影响其传热和流动性能。三、冷却壁多物理场耦合基本理论3.1传热学理论基础3.1.1超临界二氧化碳传热特性超临界二氧化碳的传热特性与常规流体相比具有显著差异，这主要源于其在超临界状态下独特的物理性质。当二氧化碳处于超临界状态时，其密度、黏度、比热容和导热系数等物性参数随温度和压力的变化呈现出高度的非线性。在近临界区域，二氧化碳的比热容会发生急剧变化，可在短时间内增大数倍。这种比热容的剧烈变化会极大地影响超临界二氧化碳的传热性能。当超临界二氧化碳在管道中流动并与管壁进行换热时，在比热容急剧变化的区域，其吸收或释放热量的能力会显著改变，从而导致传热系数发生大幅波动。热流密度也是影响超临界二氧化碳传热的重要因素。在较高的热流密度下，超临界二氧化碳的传热过程会变得更加复杂。热流密度的增加会使管壁与流体之间的温度差增大，从而加剧了流体的湍流程度。在某超临界二氧化碳冷却实验中，当热流密度从100kW/m²增加到200kW/m²时，流体的湍流强度增加了30%，这使得传热边界层的厚度减小，传热系数增大。过高的热流密度可能会导致传热恶化现象的发生。当热流密度超过一定阈值时，管壁附近的流体可能会出现局部过热，形成气膜层，气膜层的导热系数远低于液态或气态二氧化碳，从而阻碍了热量的传递，使传热系数急剧下降。质量流量对超临界二氧化碳的传热也有着重要影响。随着质量流量的增加，超临界二氧化碳的流速增大，这会增强流体的对流换热能力。在某超临界二氧化碳流动传热实验中，当质量流量从500kg/(m²・s)增加到1000kg/(m²・s)时，对流换热系数提高了约40%。这是因为流速的增大使得流体与管壁之间的接触更加频繁，热量传递更加迅速。质量流量的增加还可以减小传热边界层的厚度，进一步提高传热效率。然而，质量流量的增加也会带来一些问题，如增加了系统的能耗和压力损失。压力对超临界二氧化碳的传热特性同样具有不可忽视的影响。压力的变化会改变超临界二氧化碳的物性参数，从而影响其传热性能。当压力升高时，超临界二氧化碳的密度增大，分子间的相互作用力增强，这使得其导热系数增大，传热性能得到提升。在某研究中，当压力从20MPa升高到25MPa时，超临界二氧化碳的导热系数提高了15%左右。压力的变化还会影响超临界二氧化碳的相变行为和临界特性，进而对传热产生影响。在近临界区域，压力的微小变化可能会导致物性参数的大幅波动，从而使传热过程变得不稳定。超临界二氧化碳在水平管和垂直管中的传热特性也存在差异。在水平管中，由于重力的影响，流体在管内的分布会出现不均匀现象，导致传热不均匀。靠近管子底部的流体流速相对较低，传热系数也相对较小；而靠近管子顶部的流体流速相对较高，传热系数也相对较大。在垂直管中，流体的流动方向与重力方向一致或相反，这会影响流体的湍流程度和传热性能。当流体向上流动时，重力会对流体产生一定的阻碍作用，使流体的湍流程度增加，传热系数增大；当流体向下流动时，重力会促进流体的流动，使流体的湍流程度减小，传热系数降低。3.1.2冷却壁传热过程分析冷却壁的传热过程是一个涉及多种传热方式和复杂物理现象的过程，对其进行深入分析对于理解冷却壁的工作原理和性能具有重要意义。热传导是冷却壁传热过程中的重要环节。在冷却壁内部，热量主要通过固体材料的晶格振动和自由电子的运动进行传递。冷却壁通常由金属材料制成，如碳钢、合金钢等，这些材料具有较高的导热系数，能够有效地传导热量。在高温环境下，冷却壁表面吸收炉膛内的热量后，热量会通过热传导的方式从壁面向内部传递。根据傅里叶定律，热传导的热流密度与温度梯度成正比，与材料的导热系数成反比。对于一维平板导热，热流密度q可以表示为：q=-λ(dT/dx)，其中λ为导热系数，dT/dx为温度梯度。在冷却壁中，由于温度分布不均匀，会产生温度梯度，从而导致热传导的发生。热传导的速率还受到材料的热物性参数（如比热容、密度等）的影响。比热容较大的材料能够吸收更多的热量，从而减缓温度的变化；密度较大的材料则会使热传导的速率降低。对流换热是冷却壁传热过程中的另一个关键环节。冷却壁内的超临界二氧化碳工质通过对流换热的方式与冷却壁壁面进行热量交换。对流换热是指流体中各部分之间由于相对运动而引起的热量传递现象，它是热对流和热传导共同作用的结果。根据牛顿冷却定律，对流换热的热流密度q可以表示为：q=h(Tw-Tf)，其中h为对流换热系数，Tw为壁面温度，Tf为流体温度。对流换热系数h是衡量对流换热强弱的重要参数，它受到多种因素的影响，如流体的流速、物性参数、壁面的粗糙度和几何形状等。在超临界二氧化碳冷却壁中，由于超临界二氧化碳的特殊物性，其对流换热过程呈现出高度的非线性和复杂性。在近临界区域，超临界二氧化碳的比热容急剧变化，会导致对流换热系数的大幅波动。流速的增加会增强流体的湍流程度，从而提高对流换热系数。热辐射在冷却壁传热过程中也占有一定的比例。在高温环境下，冷却壁表面会向周围环境发射热辐射，同时也会吸收周围环境的热辐射。热辐射是一种以电磁波的形式传递热量的方式，它不需要介质的存在，可以在真空中传播。根据斯蒂芬-玻尔兹曼定律，物体的辐射热流密度q可以表示为：q=εσT⁴，其中ε为物体的发射率，σ为斯蒂芬-玻尔兹曼常数，T为物体的绝对温度。在冷却壁传热过程中，热辐射的影响主要体现在高温区域，当冷却壁表面温度较高时，热辐射的作用会更加明显。在炉膛内高温区域，冷却壁表面的热辐射热流密度可占总热流密度的10%-20%左右。热辐射的大小还与周围环境的温度和发射率有关，周围环境温度越高，冷却壁表面吸收的热辐射就越多；周围环境发射率越大，冷却壁表面与周围环境之间的辐射换热就越强。在实际的冷却壁传热过程中，热传导、对流换热和热辐射这三种传热方式并不是孤立存在的，而是相互关联、相互影响的。冷却壁表面吸收炉膛内的高温热量后，热量首先通过热传导的方式从壁面向内部传递，然后通过对流换热的方式传递给冷却壁内的超临界二氧化碳工质，在这个过程中，热辐射也会参与热量的传递。而且，这三种传热方式之间的相互作用会随着工况的变化而发生改变。在高温、高压工况下，热辐射的作用会增强，对流换热和热传导的作用也会相应发生变化。因此，在研究冷却壁的传热过程时，需要综合考虑这三种传热方式的影响，建立准确的传热模型，以深入理解冷却壁的传热机理和性能。3.2流体力学理论基础3.2.1超临界二氧化碳流动特性超临界二氧化碳的流动特性与常规流体相比存在显著差异，这些差异主要源于其独特的物理性质以及在超临界状态下的特殊行为。超临界二氧化碳的黏度对其流动特性有着重要影响。在超临界状态下，二氧化碳的黏度随温度和压力的变化呈现出复杂的规律。当温度升高时，超临界二氧化碳的黏度通常会先降低后升高。在某一研究中，当压力保持在20MPa，温度从32℃升高到50℃时，黏度先从0.012mPa・s降低到0.008mPa・s，随后在温度继续升高时，黏度又逐渐上升。这是因为在较低温度范围内，分子热运动加剧，使得分子间的相互作用力减弱，从而导致黏度降低；而当温度进一步升高时，分子间的距离减小，相互作用力增强，黏度又开始增大。压力对超临界二氧化碳黏度的影响则表现为，随着压力的增加，黏度逐渐增大。当压力从10MPa增加到30MPa时，黏度从0.006mPa・s增大到0.018mPa・s。黏度的变化会直接影响超临界二氧化碳的流动阻力，黏度越大，流动阻力越大，流体的流动就越困难。密度是影响超临界二氧化碳流动特性的另一个重要因素。超临界二氧化碳的密度在临界点附近对温度和压力的变化非常敏感。当温度和压力接近临界值时，微小的变化就能导致密度的大幅波动。在某实验中，当压力在临界压力7.38MPa附近波动±0.5MPa时，密度变化可达20%-30%。密度的变化会影响超临界二氧化碳的惯性力和重力，从而改变其流动状态。在密度较大时，流体的惯性力较大，流动更倾向于保持稳定；而当密度较小时，重力的影响相对增大，可能会导致流体的流动出现不稳定现象。超临界二氧化碳的可压缩性也是其流动特性的一个关键方面。与常规流体相比，超临界二氧化碳具有较高的可压缩性。在压力变化时，其体积会发生明显的改变。当压力从15MPa降低到10MPa时，超临界二氧化碳的体积可能会增大50%-100%。这种可压缩性会导致在流动过程中出现压力波的传播和反射，影响流体的流动稳定性。在管道中，压力波的传播可能会引起局部压力的波动，进而影响流体的流速和流量分布。超临界二氧化碳在管道中的流动还会受到壁面粗糙度的影响。壁面粗糙度会增加流体与壁面之间的摩擦力，从而增大流动阻力。在某超临界二氧化碳管道流动实验中，当壁面粗糙度从0.01mm增加到0.05mm时，流动阻力增加了约30%。壁面粗糙度还会影响流体的边界层结构，使得边界层内的速度分布发生变化，进一步影响流体的流动特性。3.2.2冷却壁内流场分析冷却壁内流场的分析对于理解超临界二氧化碳在冷却壁内的流动行为以及优化冷却壁的设计和运行具有重要意义。在冷却壁内，超临界二氧化碳的流速分布呈现出一定的规律。靠近冷却壁壁面的流体流速较低，形成了一个速度边界层。这是因为壁面对流体存在摩擦力，使得流体的流速逐渐降低。在速度边界层内，流速从壁面处的零值逐渐增大到主流区的速度。而在冷却壁管道的中心区域，流体的流速较高，且分布相对均匀。根据相关理论和实验研究，对于圆形管道内的充分发展湍流流动，流速分布可以用幂律分布来近似描述，即u/umax=(1-r/R)^(1/n)，其中u为距管道中心距离为r处的流速，umax为管道中心处的最大流速，R为管道半径，n为与雷诺数有关的指数，通常在5-10之间。在超临界二氧化碳冷却壁中，由于其特殊的物性和流动状态，n的值可能会有所不同，需要通过实验或数值模拟来确定。压力分布也是冷却壁内流场分析的重要内容。随着超临界二氧化碳在冷却壁内的流动，由于流动阻力的存在，压力会逐渐降低。压力损失主要包括沿程压力损失和局部压力损失。沿程压力损失是由于流体与壁面之间的摩擦力以及流体内部的黏性力造成的，其大小与管道长度、流速、流体物性以及壁面粗糙度等因素有关。局部压力损失则是在管道的弯头、阀门、变径等部位，由于流体的流动方向和速度发生突然变化，产生了涡流和能量耗散，从而导致压力的降低。在某冷却壁实验中，通过测量不同位置的压力，发现沿程压力损失在整个压力损失中占比较大，约为70%-80%，而局部压力损失主要集中在管道的弯头处，占总压力损失的20%-30%。冷却壁内超临界二氧化碳的流量分配对于系统的稳定运行至关重要。在多通道冷却壁中，如果流量分配不均匀，会导致部分通道内的流体流速过低或过高，从而影响冷却效果和系统的安全性。流量分配不均匀的原因主要包括管道的几何形状差异、阻力系数不同以及入口条件的不均匀性等。为了实现均匀的流量分配，可以采取一些措施，如优化管道的布置和结构，使各通道的阻力尽量相等；在入口处设置均流器，改善入口条件的均匀性。通过数值模拟研究发现，在冷却壁入口处设置一个合适的均流器，可以使各通道内的流量偏差控制在5%以内，有效提高了流量分配的均匀性。冷却壁内的流场还会受到冷却壁结构形式的影响。不同结构形式的冷却壁，如螺旋管圈冷却壁、垂直管圈冷却壁和膜式冷却壁等，其内部的流道形状和几何参数不同，会导致流场特性存在差异。螺旋管圈冷却壁由于其螺旋上升的流道结构，会使流体在流动过程中产生离心力，从而影响流速分布和压力分布。而垂直管圈冷却壁的流道相对简单，流体的流动较为规则，但在管屏之间可能会存在流量分配不均匀的问题。膜式冷却壁由于其紧密的结构，流道内的流动阻力相对较大，但具有较好的密封性和整体性。3.3热应力与结构力学理论基础3.3.1热应力产生机制热应力是由于物体内部温度分布不均匀而产生的应力。在超临界二氧化碳燃煤发电系统冷却壁中，热应力的产生主要源于温度梯度的存在。当冷却壁受到高温烟气的加热时，其表面温度迅速升高，而内部温度由于热传导的延迟，升高速度较慢，从而在冷却壁内部形成了温度梯度。这种温度梯度会导致冷却壁不同部位的材料膨胀或收缩程度不同，由于材料之间的相互约束，无法自由膨胀或收缩，从而产生热应力。以一个简单的平板冷却壁为例，假设平板一侧受到高温烟气的加热，温度为T1，另一侧与超临界二氧化碳工质接触，温度为T2，且T1>T2。根据热膨胀原理，温度较高的一侧材料会有更大的膨胀趋势，而温度较低的一侧材料膨胀较小。在这种情况下，材料内部会产生应力，以平衡这种膨胀差异。对于各向同性的材料，热应力可以通过以下公式计算：\sigma=E\alpha\DeltaT其中，\sigma为热应力，E为材料的弹性模量，\alpha为材料的线性热膨胀系数，\DeltaT为温度变化。在冷却壁中，\DeltaT即为冷却壁内外表面的温度差。例如，对于某钢材制成的冷却壁，其弹性模量E=200GPa，线性热膨胀系数\alpha=1.2\times10^{-5}/^{\circ}C，当冷却壁内外表面温度差\DeltaT=200^{\circ}C时，根据上述公式可计算出热应力\sigma=200\times10^{9}\times1.2\times10^{-5}\times200=4.8\times10^{8}Pa。在实际的冷却壁中，温度分布往往是三维的，热应力的计算也更为复杂。需要考虑材料的各向异性、几何形状以及边界条件等因素。此时，通常采用有限元方法等数值计算方法来求解热应力。通过将冷却壁离散为多个小单元，在每个单元内应用热传导方程和力学平衡方程，求解出单元内的温度分布和应力分布，进而得到整个冷却壁的热应力分布。在某超临界二氧化碳燃煤发电系统冷却壁的有限元分析中，考虑了冷却壁的复杂几何形状和非均匀的温度分布，通过数值计算得到了冷却壁不同部位的热应力分布情况，发现热应力集中主要出现在冷却壁的拐角处和管子与管板的连接处，这些部位的热应力值明显高于其他部位。3.3.2结构力学分析方法在超临界二氧化碳燃煤发电系统冷却壁的研究中，结构力学分析方法用于评估冷却壁在各种载荷作用下的强度和稳定性。常用的结构力学分析方法包括材料力学方法、弹性力学方法和有限元方法。材料力学方法是一种基于简化假设的分析方法，适用于简单结构的初步分析。对于冷却壁中的管道结构，可以将其视为梁或圆柱壳，利用材料力学中的公式来计算应力和变形。对于承受内压的圆形管道，其周向应力\sigma_{\theta}和轴向应力\sigma_{z}可以分别用以下公式计算：\sigma_{\theta}=\frac{pD}{2t}\sigma_{z}=\frac{pD}{4t}其中，p为管道内压力，D为管道外径，t为管道壁厚。通过这些公式，可以初步估算管道在压力作用下的应力水平，判断其是否满足强度要求。例如，对于一个内径为100mm，壁厚为5mm，内压为25MPa的冷却壁管道，根据上述公式计算得到周向应力\sigma_{\theta}=\frac{25\times10^{6}\times110\times10^{-3}}{2\times5\times10^{-3}}=2.75\times10^{8}Pa，轴向应力\sigma_{z}=\frac{25\times10^{6}\times110\times10^{-3}}{4\times5\times10^{-3}}=1.375\times10^{8}Pa。将计算得到的应力值与材料的许用应力进行比较，可以判断管道是否安全。弹性力学方法是一种基于弹性理论的精确分析方法，能够考虑结构的复杂几何形状和边界条件。它通过求解弹性力学的基本方程，得到结构内的应力、应变和位移分布。对于冷却壁这种复杂结构，弹性力学的解析解往往很难得到，通常需要采用数值方法进行求解。在某冷却壁的弹性力学分析中，利用数值方法求解了弹性力学的基本方程，得到了冷却壁在热载荷和机械载荷共同作用下的应力分布，结果表明，弹性力学方法能够更准确地反映冷却壁内部的应力状态，但计算过程较为复杂，需要较高的数学和计算能力。有限元方法是一种广泛应用的数值分析方法，它将连续的结构离散为有限个单元，通过对每个单元的力学分析，得到整个结构的力学响应。有限元方法可以处理复杂的几何形状、材料特性和边界条件，能够准确地计算冷却壁的应力、应变和变形。在超临界二氧化碳燃煤发电系统冷却壁的分析中，利用有限元软件ANSYS对冷却壁进行建模和分析。首先，根据冷却壁的实际结构建立三维几何模型，并对其进行网格划分。然后，定义材料属性、边界条件和载荷工况，包括温度载荷、压力载荷等。最后，求解有限元方程，得到冷却壁的应力、应变和位移分布。通过有限元分析，可以直观地看到冷却壁在不同工况下的力学响应，为冷却壁的设计和优化提供重要依据。例如，通过有限元分析发现，在高温、高压工况下，冷却壁的某些部位出现了应力集中现象，通过优化冷却壁的结构设计，如改变管道的布置方式、增加加强筋等，可以有效地降低应力集中，提高冷却壁的强度和稳定性。四、冷却壁多物理场耦合特性研究4.1数值模拟研究方法与模型建立4.1.1数值模拟软件选择在对超临界二氧化碳燃煤火力发电系统冷却壁多物理场耦合特性的研究中，选择COMSOLMultiphysics软件作为数值模拟工具。该软件在多物理场耦合模拟领域展现出卓越的优势，能够有效应对冷却壁复杂的物理过程模拟需求。COMSOLMultiphysics具备强大的多物理场直接耦合分析能力，这使得它在处理冷却壁多物理场问题时具有显著优势。冷却壁工作过程涉及热传导、对流换热、流体流动以及热应力等多个物理场的相互作用，这些物理场之间存在复杂的耦合关系。而COMSOLMultiphysics允许用户在同一模型中，通过选择或自定义不同专业的偏微分方程进行任意组合，从而轻松实现这些多物理场的直接耦合分析。例如，在模拟冷却壁内超临界二氧化碳的流动与传热过程时，软件能够同时考虑流体流动的Navier-Stokes方程与传热的能量方程之间的耦合关系，精确模拟出超临界二氧化碳在冷却壁内的复杂流动和传热行为，以及它们之间的相互影响。该软件拥有丰富的物理场接口和预定义模型，涵盖了传热学、流体力学、结构力学等多个领域，这为冷却壁多物理场耦合模拟提供了便利。在研究冷却壁的热传导过程时，可以直接使用软件中的固体传热模块，该模块基于傅里叶定律，能够准确描述热量在冷却壁固体材料中的传导过程。在分析冷却壁内超临界二氧化碳的对流换热和流体流动时，软件的流体流动模块和传热模块能够提供全面的物理模型和计算方法。对于冷却壁的热应力分析，结构力学模块可以根据材料的力学性能和温度分布，准确计算出热应力的大小和分布情况。这些丰富的物理场接口和预定义模型，大大减少了用户自行建立复杂物理模型的工作量，提高了模拟的准确性和效率。COMSOLMultiphysics还具有出色的网格划分和自适应网格功能。在对冷却壁进行数值模拟时，网格划分的质量直接影响到模拟结果的准确性和计算效率。该软件提供了多种网格划分方法，包括结构化网格、非结构化网格和混合网格等，用户可以根据冷却壁的几何形状和物理场分布特点，选择合适的网格划分方法。软件的自适应网格功能能够根据模拟过程中物理场的变化，自动调整网格的精细度。在冷却壁内温度梯度较大的区域，自适应网格功能会自动加密网格，提高该区域的计算精度；而在物理场变化较为平缓的区域，则适当降低网格密度，减少计算量。通过这种方式，既保证了模拟结果的准确性，又提高了计算效率，优化了计算资源的使用。软件具备全面的第三方CAD导入功能，支持当前主流CAD软件格式文件的导入。在研究冷却壁多物理场耦合特性时，通常需要根据冷却壁的实际结构建立几何模型。利用COMSOLMultiphysics的CAD导入功能，可以直接将在其他CAD软件中设计好的冷却壁几何模型导入到软件中，避免了在软件中重新创建复杂几何模型的繁琐过程。而且，软件能够对导入的几何模型进行有效的处理和修复，确保模型的完整性和准确性，为后续的数值模拟提供可靠的几何基础。COMSOLMultiphysics拥有强大的后处理功能，可根据用户的需要进行各种数据、曲线、及动画的输出与分析。在完成冷却壁多物理场耦合模拟后，用户可以通过软件的后处理功能，直观地查看冷却壁内的温度分布、速度场、压力分布以及热应力分布等模拟结果。可以将这些结果以二维或三维图形、数据表格、曲线等形式输出，方便用户进行数据分析和比较。软件还支持动画输出功能，能够将模拟过程中的物理场变化以动画的形式展示出来，更加生动形象地呈现冷却壁多物理场耦合的动态过程，有助于用户深入理解多物理场之间的相互作用机制。4.1.2模型建立与网格划分在进行冷却壁多物理场耦合特性的数值模拟研究时，准确建立模型并进行合理的网格划分是确保模拟结果可靠性的关键步骤。首先，根据冷却壁的实际结构和尺寸，利用三维建模软件（如SolidWorks、AutoCAD等）建立冷却壁的几何模型。以常见的螺旋管圈冷却壁为例，在建模过程中，精确绘制螺旋上升的冷却水管，包括管子的直径、螺旋升角、盘旋圈数等参数，确保与实际结构一致。同时，详细构建冷却壁的壁体结构，考虑壁体的厚度、材质分布等因素。对于冷却壁与其他部件（如集箱、管道连接件等）的连接部分，也进行准确的建模，以反映实际的边界条件和热传递路径。完成几何模型的创建后，将其保存为COMSOLMultiphysics软件支持的格式（如.stl、.iges等），以便后续导入进行数值模拟。将几何模型导入COMSOLMultiphysics软件后，需要进行网格划分。在网格划分策略上，综合考虑模型的几何复杂度、物理场分布特点以及计算资源的限制，选择合适的网格类型和参数。对于冷却壁的壁体部分，由于其几何形状相对规则，且热传导过程较为平稳，采用结构化网格进行划分。结构化网格具有规则的网格布局，计算效率高，能够准确地捕捉壁体内的温度分布。在划分时，根据壁体的尺寸和温度梯度的大小，合理设置网格的尺寸。对于温度梯度较大的区域，如靠近高温烟气侧的壁面，适当减小网格尺寸，以提高计算精度；而在温度梯度较小的区域，增大网格尺寸，减少计算量。对于冷却壁内的冷却水管部分，由于其内部超临界二氧化碳的流动和传热过程较为复杂，采用非结构化网格进行划分。非结构化网格能够更好地适应复杂的几何形状和物理场变化，在冷却水管的弯曲部位、进出口等关键区域，可以灵活地加密网格，准确地模拟超临界二氧化碳的流动和传热特性。在划分冷却水管的网格时，根据超临界二氧化碳的流速、温度变化以及压力分布等因素，动态调整网格的密度。在流速变化较大、温度梯度明显的区域，如管道的弯头处和热负荷较高的部位，增加网格的数量，提高网格的分辨率，以确保模拟结果的准确性。为了进一步提高模拟精度，采用局部网格加密技术。在冷却壁与高温烟气接触的表面、冷却水管与壁体的连接处等关键部位，进行局部网格加密。这些部位通常存在较大的温度梯度和热应力集中现象，对冷却壁的性能和寿命影响较大。通过局部网格加密，可以更精确地计算这些区域的物理量分布，捕捉到细微的物理变化。在COMSOLMultiphysics软件中，可以通过定义局部网格控制区域，设置网格细化的级别和范围，实现对关键部位的网格加密。例如，在冷却壁与高温烟气接触的表面，将网格尺寸细化为其他区域的1/2或1/3，以提高该区域的计算精度。在完成网格划分后，对生成的网格进行质量检查。通过查看网格的形状质量、尺寸一致性、正交性等指标，评估网格的质量。确保网格单元的形状尽量规则，避免出现严重扭曲或畸形的单元。相邻网格单元的尺寸变化应保持在合理范围内，以减少数值计算中的误差。对于不符合质量要求的网格，进行调整和优化，直到满足模拟计算的要求。例如，如果发现某个区域的网格正交性较差，可以通过调整网格划分参数或使用网格优化工具，改善该区域的网格质量。4.1.3边界条件与参数设置在冷却壁多物理场耦合数值模拟中，准确设定边界条件和相关参数是保证模拟结果准确性的关键环节，它直接关系到模拟结果与实际工况的契合程度。在热边界条件方面，考虑到冷却壁一侧直接与高温烟气接触，另一侧与超临界二氧化碳工质换热，分别对这两侧设置不同的边界条件。对于与高温烟气接触的冷却壁表面，采用第三类边界条件，即对流换热边界条件。根据实际运行数据和相关研究，确定高温烟气与冷却壁表面之间的对流换热系数。某超临界二氧化碳燃煤发电系统中，高温烟气温度约为1200℃，与冷却壁表面的对流换热系数取值为500-800W/(m²・K)。通过该边界条件，能够准确模拟高温烟气向冷却壁传递热量的过程。对于冷却壁与超临界二氧化碳工质接触的表面，同样采用对流换热边界条件，根据超临界二氧化碳的流动特性和换热规律，确定其与冷却壁表面之间的对流换热系数。在不同的工况下，超临界二氧化碳的流速、温度等参数会发生变化，从而导致对流换热系数也相应改变。通过实验研究和理论分析，建立对流换热系数与这些参数之间的关联式，以便在模拟中根据实际工况准确设定对流换热系数。在流体边界条件方面，对于冷却壁内超临界二氧化碳的入口，设置质量流量入口边界条件。根据发电系统的设计参数和实际运行要求，确定超临界二氧化碳的入口质量流量。在某600MW超临界二氧化碳燃煤发电项目中，冷却壁内超临界二氧化碳的入口质量流量为100-150kg/s。通过设定质量流量入口边界条件，能够准确控制超临界二氧化碳进入冷却壁的流量，模拟其在冷却壁内的流动过程。对于超临界二氧化碳的出口，设置压力出口边界条件，根据系统的背压和运行压力要求，确定出口压力。一般情况下，出口压力略低于入口压力，以保证超临界二氧化碳能够顺利流出冷却壁。在结构力学边界条件方面，考虑到冷却壁在运行过程中受到自身重力、内部压力以及热应力等多种载荷的作用，对冷却壁的支撑部位设置固定约束边界条件。在冷却壁与集箱、管道连接件等部件的连接处，将这些部位的自由度进行约束，模拟实际的支撑情况。这样可以确保在模拟过程中，冷却壁能够按照实际的受力情况进行变形和应力分布计算。在参数设置方面，准确设定超临界二氧化碳的物性参数至关重要。超临界二氧化碳的密度、黏度、比热容、导热系数等物性参数随温度和压力的变化而变化，这些变化会显著影响其在冷却壁内的流动和传热特性。通过查阅相关文献资料和实验数据，获取超临界二氧化碳在不同温度和压力条件下的物性参数，并将其输入到数值模拟模型中。利用NISTREFPROP软件等物性计算工具，根据模拟的温度和压力范围，精确计算超临界二氧化碳的物性参数。对于冷却壁的材料参数，如导热系数、弹性模量、热膨胀系数等，根据所选用的材料类型，参考材料手册和相关标准，准确设定这些参数。如果冷却壁采用碳钢材料，其导热系数约为45W/(m・K)，弹性模量为200GPa，热膨胀系数为1.2×10⁻⁵/℃。在模拟过程中，还需要考虑材料参数随温度的变化，对于一些高温合金材料，其弹性模量和热膨胀系数会随着温度的升高而发生变化，通过建立材料参数与温度的函数关系，在模拟中准确反映这种变化。4.2模拟结果与分析4.2.1温度场分布特性通过数值模拟，得到了超临界二氧化碳燃煤发电系统冷却壁在不同工况下的温度场分布云图，图1展示了典型工况下冷却壁的温度场分布情况。从图中可以清晰地看到，冷却壁的温度分布呈现出明显的不均匀性。靠近高温烟气侧的冷却壁表面温度较高，达到了500-600℃，这是因为高温烟气直接对冷却壁表面进行加热，热量通过热传导和对流换热的方式传递到冷却壁内部。而靠近超临界二氧化碳工质侧的冷却壁表面温度较低，一般在100-200℃之间，这是由于超临界二氧化碳工质具有良好的吸热能力，能够有效地带走冷却壁表面的热量。在冷却壁的厚度方向上，温度呈现出逐渐降低的趋势。从高温烟气侧到超临界二氧化碳工质侧，温度梯度较大，这表明热量在冷却壁内的传递过程中存在较大的阻力。在某一位置处，冷却壁厚度为50mm，高温烟气侧温度为550℃，超临界二氧化碳工质侧温度为150℃，通过计算可得该位置处的温度梯度约为8000℃/m。这种较大的温度梯度会导致冷却壁内部产生热应力，对冷却壁的结构完整性产生影响。冷却壁的温度分布还受到超临界二氧化碳工质流速的影响。当工质流速增加时，对流换热系数增大，冷却壁与工质之间的换热效率提高，从而使冷却壁的温度降低。在某一模拟工况中，将超临界二氧化碳工质的流速从2m/s增加到4m/s，冷却壁的平均温度降低了约20℃。流速的增加也会使冷却壁内的温度分布更加均匀，减少温度梯度。这是因为流速的增加增强了工质的湍流程度，使得热量能够更均匀地传递。热流密度对冷却壁温度分布的影响也十分显著。当热流密度增大时，冷却壁表面吸收的热量增加，导致冷却壁的温度升高。在某实验中，将热流密度从200kW/m²增加到300kW/m²，冷却壁的最高温度升高了约50℃。热流密度的变化还会改变冷却壁内的温度梯度分布。在热流密度较大的区域，温度梯度也相应增大，这会加剧热应力的产生。冷却壁的结构形式也会对温度分布产生影响。不同结构形式的冷却壁，其内部的流道结构和换热面积不同，导致温度分布存在差异。螺旋管圈冷却壁由于其螺旋上升的流道结构，使得超临界二氧化碳工质在流动过程中能够更充分地与冷却壁壁面接触，换热效果较好，温度分布相对较为均匀。而垂直管圈冷却壁在管屏之间可能存在流量分配不均匀的问题，导致部分区域的换热效果较差，温度分布不均匀。[此处插入图1：典型工况下冷却壁的温度场分布云图]4.2.2速度场与压力场分布特性冷却壁内超临界二氧化碳的速度场分布呈现出一定的规律。在冷却壁管道的中心区域，超临界二氧化碳的流速较高，而靠近壁面的区域，由于壁面摩擦力的作用，流速逐渐降低，形成了速度边界层。图2展示了冷却壁内某一截面的速度矢量图，可以清晰地看到速度的分布情况。在管道中心处，超临界二氧化碳的流速可达5-8m/s，而在靠近壁面的边界层内，流速迅速降低，在壁面处流速趋近于零。速度场的分布还受到冷却壁结构形式的影响。对于螺旋管圈冷却壁，由于其螺旋上升的流道结构，超临界二氧化碳在流动过程中会受到离心力的作用，使得靠近外侧的流速相对较高，而靠近内侧的流速相对较低。在螺旋管圈冷却壁的某一截面处，通过数值模拟得到靠近外侧的流速比靠近内侧的流速高约1-2m/s。这种流速分布的不均匀性会影响冷却壁的换热效果，导致冷却壁温度分布的不均匀。冷却壁内的压力场分布也与速度场密切相关。随着超临界二氧化碳在冷却壁内的流动，由于流动阻力的存在，压力逐渐降低。图3展示了冷却壁内的压力分布云图。在冷却壁入口处，压力较高，一般在20-25MPa之间，而在出口处，压力降低到18-22MPa左右。压力损失主要包括沿程压力损失和局部压力损失。沿程压力损失是由于流体与壁面之间的摩擦力以及流体内部的黏性力造成的，其大小与管道长度、流速、流体物性以及壁面粗糙度等因素有关。局部压力损失则主要发生在管道的弯头、阀门、变径等部位，这些部位由于流体的流动方向和速度发生突然变化，产生了涡流和能量耗散，从而导致压力降低。在某冷却壁实验中，通过测量发现，沿程压力损失占总压力损失的70%-80%，而局部压力损失主要集中在管道的弯头处，占总压力损失的20%-30%。压力场的分布对冷却壁的结构和性能也有重要影响。过高的压力会增加冷却壁的负荷，对冷却壁的强度和密封性提出更高的要求。压力分布的不均匀可能会导致冷却壁局部受力不均，从而产生变形或损坏。在冷却壁的设计和运行过程中，需要合理控制压力场的分布，确保冷却壁的安全稳定运行。[此处插入图2：冷却壁内某一截面的速度矢量图][此处插入图3：冷却壁内的压力分布云图]4.2.3热应力场分布特性通过数值模拟得到了冷却壁在不同工况下的热应力场分布情况，图4展示了典型工况下冷却壁的热应力分布云图。从图中可以看出，冷却壁的热应力分布呈现出明显的不均匀性。在冷却壁的拐角处和管子与管板的连接处，热应力集中现象较为明显，这些部位的热应力值明显高于其他部位。在冷却壁的拐角处，热应力可达100-150MPa，而在管子与管板的连接处，热应力甚至可超过200MPa。热应力的大小和分布主要受到温度梯度的影响。由于冷却壁在高温烟气和超临界二氧化碳工质的作用下，温度分布不均匀，从而产生了温度梯度。温度梯度越大，热应力也就越大。在冷却壁的高温烟气侧，温度变化较为剧烈，温度梯度较大，因此热应力也较大。在某一位置处，温度梯度为5000℃/m，根据热应力计算公式，可计算出该位置处的热应力约为60MPa。冷却壁的材料属性也会对热应力分布产生影响。不同材料的弹性模量、热膨胀系数等参数不同，导致在相同的温度变化下，产生的热应力也不同。一般来说，弹性模量越大，热应力越大；热膨胀系数越大，热应力也越大。如果冷却壁采用弹性模量较高的材料，在相同的温度梯度下，热应力会比采用弹性模量较低的材料时更大。热应力对冷却壁的寿命有着重要影响。长期处于高应力状态下，冷却壁容易出现疲劳裂纹，随着裂纹的扩展，最终可能导致冷却壁的损坏。根据疲劳寿命理论，热应力的大小和循环次数会影响冷却壁的疲劳寿命。在某一工况下，通过疲劳寿命计算软件，结合冷却壁的热应力分布和工作循环次数，预测出冷却壁的疲劳寿命为5-8年。为了延长冷却壁的寿命，需要采取措施降低热应力，如优化冷却壁的结构设计，减少温度梯度；选择合适的材料，降低材料的热膨胀系数和弹性模量等。[此处插入图4：典型工况下冷却壁的热应力分布云图]4.3实验研究4.3.1实验装置与方案设计实验装置主要由超临界二氧化碳供应系统、加热系统、冷却壁试验段、测量系统以及数据采集与处理系统等部分组成。超临界二氧化碳供应系统用于提供稳定的超临界二氧化碳工质。该系统主要包括二氧化碳储罐、高压泵、预热器等设备。二氧化碳储罐储存液态二氧化碳，高压泵将液态二氧化碳加压至超临界状态，预热器则对超临界二氧化碳进行预热，使其达到实验所需的温度。在实验过程中，通过调节高压泵的流量和预热器的加热功率，精确控制超临界二氧化碳的压力和温度。加热系统用于模拟燃煤锅炉中的高温环境，为冷却壁试验段提供高温热源。加热系统采用电加热方式，由加热炉、温控器等组成。加热炉内部设置有加热丝，通过温控器控制加热丝的加热功率，使加热炉内部的温度达到实验设定的高温值。在加热炉内部，还安装有温度传感器，实时监测炉内温度，确保温度的稳定性和准确性。冷却壁试验段是实验的核心部分，用于模拟超临界二氧化碳燃煤发电系统中的冷却壁。试验段采用与实际冷却壁相同的结构和材料，以保证实验结果的真实性和可靠性。冷却壁试验段的尺寸根据实验需求进行设计，内部设置有超临界二氧化碳的流道，流道的形状和尺寸也与实际冷却壁一致。在冷却壁试验段的外壁面和内壁面，分别安装有多个温度传感器，用于测量壁面温度。在流道内，安装有压力传感器和流速传感器，用于测量超临界二氧化碳的压力和流速。测量系统由各种传感器组成，用于测量实验过程中的各种物理参数。温度传感器采用高精度的热电偶或热电阻，分别安装在冷却壁试验段的不同位置，包括外壁面、内壁面以及超临界二氧化碳流道内，用于测量不同位置的温度。压力传感器采用高精度的压力变送器，安装在超临界二氧化碳流道的入口和出口以及不同位置，用于测量超临界二氧化碳的压力。流速传感器采用电磁流量计或超声波流量计，安装在超临界二氧化碳流道内，用于测量超临界二氧化碳的流速。数据采集与处理系统用于采集和处理测量系统获取的实验数据。该系统由数据采集卡、计算机以及数据处理软件等组成。数据采集卡将传感器输出的电信号转换为数字信号，并传输到计算机中。计算机通过数据处理软件对采集到的数据进行实时监测、记录和分析。数据处理软件可以对数据进行滤波、平滑、插值等处理，以提高数据的准确性和可靠性。还可以根据实验数据绘制各种图表，如温度随时间的变化曲线、压力随流速的变化曲线等，直观地展示实验结果。实验方案设计如下：设置不同的实验工况，包括不同的超临界二氧化碳压力、温度、质量流量以及加热炉的热流密度等。在每个工况下，稳定运行实验装置，待各项参数稳定后，开始采集数据。每个工况下采集的数据持续时间不少于30分钟，以确保数据的稳定性和可靠性。在实验过程中，记录每个工况下冷却壁试验段的温度分布、超临界二氧化碳的流速和压力分布等数据。对采集到的数据进行整理和分析，研究不同工况下冷却壁的多物理场耦合特性。将实验结果与数值模拟结果进行对比，验证数值模拟模型的准确性和可靠性。4.3.2实验结果与数值模拟对比验证将实验得到的冷却壁温度分布数据与数值模拟结果进行对比，图5展示了在某一工况下，冷却壁沿厚度方向的温度分布对比情况。从图中可以看出，实验值与模拟值在整体趋势上基本一致，都呈现出从高温烟气侧到超临界二氧化碳工质侧温度逐渐降低的趋势。在高温烟气侧，实验测量的温度为480-520℃，模拟值为500-530℃，两者较为接近。在超临界二氧化碳工质侧，实验测量的温度为120-150℃，模拟值为130-160℃，也具有较好的一致性。对于冷却壁内超临界二氧化碳的流速分布，实验结果与数值模拟结果也具有较好的吻合度。图6展示了冷却壁内某一截面的流速分布对比情况。在管道中心区域，实验测量的流速为6-8m/s，模拟值为7-9m/s，两者相差不大。在靠近壁面的边界层区域，实验和模拟结果都显示流速逐渐降低，且变化趋势相似。在压力分布方面，实验测量的冷却壁入口压力为22-23MPa，出口压力为20-21MPa，数值模拟得到的入口压力为22.5-23.5MPa，出口压力为20.5-21.5MPa，实验值与模拟值的偏差在可接受范围内。实验结果和数值模拟结果都表明，压力沿冷却壁流动方向逐渐降