航空发动机叶片冷却设计方法论文

航空发动机叶片冷却设计方法论文一.摘要

航空发动机作为飞机的核心动力装置，其性能直接关系到飞行器的整体效能与安全性。叶片冷却技术是提升航空发动机性能的关键手段，特别是在高温、高负荷的运行环境下，有效的冷却设计能够显著延长叶片使用寿命，提高发动机推重比。本研究以某型先进军用航空发动机为背景，针对其高压涡轮叶片冷却系统进行了深入分析。研究方法主要包括数值模拟与实验验证相结合的技术路线，通过建立高精度的CFD模型，对叶片内部冷却通道的流动与传热特性进行了详细模拟，同时结合风洞实验数据，对模型的准确性进行了验证与修正。研究发现，叶片表面微结构的优化设计能够显著提升冷却效率，特别是在叶片压力面与吸力面不同区域的温度分布上，微结构能够有效降低热应力集中现象。此外，通过优化冷却气的流量分配，可以在保证冷却效果的同时降低冷却损失，从而提高发动机的总体效率。研究结果表明，采用多级冷却策略与智能变循环设计相结合的方法，能够显著提升叶片的耐久性与可靠性。基于上述发现，本论文提出了针对性的叶片冷却优化方案，为实际工程设计提供了重要的理论依据和技术参考，对于推动航空发动机技术的进一步发展具有重要的实践意义。

二.关键词

航空发动机；叶片冷却；数值模拟；微结构设计；热应力；冷却效率

三.引言

航空发动机作为现代航空工业的基石，其性能水平是衡量一个国家综合国力的重要标志之一。随着航空技术的飞速发展，对发动机性能的要求日益苛刻，特别是在军用飞机领域，推重比和可靠性成为核心指标。在极高的工作温度和巨大的热负荷作用下，发动机内部叶片的冷却问题成为限制性能提升的关键瓶颈。叶片是航空发动机中最为关键且最为恶劣工作的部件之一，它直接承受着数千转每分钟的高速旋转和数千摄氏度的高温燃气冲刷，其工作环境之严苛，在航空航天领域堪称典范。在这样的条件下，叶片内部的温度梯度极大，容易引发热应力、蠕变和热疲劳等问题，进而导致叶片过早失效，严重威胁飞行安全。因此，如何有效降低叶片表面温度，提高其耐久性，成为航空发动机设计领域必须解决的核心难题。

叶片冷却技术自航空发动机诞生之初便伴随着其发展历程，历经了从简单到复杂、从被动到主动的演变过程。早期的冷却方式主要依赖于简单的气膜冷却，即通过在叶片表面开设冷却气孔，将冷却气直接喷射到热端区域进行降温。然而，随着发动机推力需求的不断增长，燃气温度和热负荷持续攀升，这种简单的冷却方式已难以满足要求，叶片表面温度过高的问题日益突出。为了应对这一挑战，研究人员逐渐探索出更为先进的冷却技术，如内部通道冷却、冲击冷却、气膜冷却与内部冷却相结合的多级冷却策略等。现代先进航空发动机普遍采用复杂的内部冷却系统，通过精心设计的冷却通道网络，将冷却气输送到叶片内部最需要降温的位置，并通过叶片表面的微结构（如扰流柱、锯齿槽等）进一步强化传热效果。

近年来，随着计算流体力学（CFD）和计算热力学（CHT）技术的飞速发展，对叶片冷却系统的设计和优化手段也得到了极大的提升。CFD技术能够在计算机上模拟冷却气在叶片内部通道和表面微结构中的流动与传热过程，为研究人员提供了直观、高效的分析工具。通过CFD模拟，可以精确预测叶片表面的温度分布、冷却气体的流动损失以及微结构对传热效率的影响，从而指导冷却系统的优化设计。同时，实验验证仍然是叶片冷却研究不可或缺的一环，通过风洞实验可以获取真实的流动和热环境数据，对CFD模型进行验证和修正，确保模拟结果的准确性。数值模拟与实验验证相结合的方法，已经成为现代叶片冷却设计中不可或缺的环节。

尽管叶片冷却技术已经取得了长足的进步，但在面对未来更高推力、更高效率、更长时间运行的发动机需求时，仍然面临诸多挑战。例如，如何在保证冷却效果的同时，最大限度地降低冷却气体的流动损失，以提升发动机的总体效率；如何在复杂的几何结构和高速旋转环境下，精确预测热应力分布，避免叶片发生失效；如何开发新型微结构材料和设计方法，进一步提升冷却性能和耐久性等。这些问题不仅涉及流体力学、传热学、材料科学等多个学科的交叉融合，也对数值模拟技术和实验研究方法提出了更高的要求。

本研究聚焦于某型先进军用航空发动机的高压涡轮叶片冷却系统，旨在通过数值模拟和实验验证相结合的方法，深入分析叶片内部冷却通道的流动与传热特性，探索优化叶片冷却性能的有效途径。具体而言，本研究将重点关注以下几个方面：首先，建立高精度的CFD模型，详细模拟冷却气在叶片内部通道和表面微结构中的流动与传热过程，分析不同设计参数对冷却效果的影响；其次，结合风洞实验数据，对CFD模型进行验证和修正，确保模拟结果的准确性；再次，通过优化叶片表面的微结构设计，探索提升冷却效率、降低热应力的有效方法；最后，基于上述研究结果，提出针对性的叶片冷却优化方案，为实际工程设计提供重要的理论依据和技术参考。本研究的意义在于，通过对叶片冷却系统的深入分析和优化，可以有效提升叶片的耐久性和可靠性，进而提高航空发动机的整体性能，对于推动我国航空发动机技术的进一步发展具有重要的理论意义和实践价值。

本研究的主要假设是：通过优化叶片内部冷却通道的设计和表面微结构的形状与布局，可以显著提升冷却效率，降低叶片表面的最高温度和温度梯度，从而有效缓解热应力集中现象，延长叶片的使用寿命。为了验证这一假设，本研究将采用数值模拟和实验验证相结合的方法，对叶片冷却系统进行系统性的研究和分析。通过这一研究，期望能够为航空发动机叶片冷却设计提供新的思路和方法，为我国航空发动机技术的自主可控做出贡献。

四.文献综述

航空发动机叶片冷却技术的研究历史悠久，伴随着发动机性能的提升而不断发展。早期的研究主要集中在简单冷却方式的效果分析上。20世纪中叶，随着航空发动机推力需求的增加，气膜冷却作为一种有效的冷却技术被广泛应用。Bakker等人对简单气膜孔的冷却效果进行了深入研究，分析了孔径、孔距等因素对冷却效率的影响，为气膜冷却的初步设计提供了理论依据。随后，研究人员开始探索更复杂的冷却结构，如多排气膜孔、内部通道冷却与气膜冷却相结合的多级冷却策略等。这些研究主要集中在优化冷却气的流量分配和改善冷却气与热气的掺混效果上，以提升冷却效率并降低流动损失。

随着计算流体力学（CFD）技术的兴起，叶片冷却的研究进入了新的阶段。CFD技术能够精确模拟冷却气在叶片内部通道和表面微结构中的流动与传热过程，为冷却系统的设计和优化提供了强大的工具。Shih等人利用CFD技术研究了不同微结构形状（如扰流柱、锯齿槽等）对叶片表面传热的影响，发现这些微结构能够显著提升冷却效率。随后，nhiềunghiêncứutiếptục深化了微结构设计的研究，如Li等人通过优化扰流柱的排布和形状，进一步提升了冷却效果。这些研究为叶片冷却的微结构设计提供了重要的理论指导。

然而，尽管CFD技术在叶片冷却研究中取得了显著的进展，但在某些方面仍然存在局限性。首先，CFD模拟需要大量的计算资源和时间，对于复杂几何结构的叶片，模拟计算量巨大，难以满足实际工程设计的快速需求。其次，CFD模型的准确性依赖于网格质量和边界条件的设定，对于一些复杂流动现象（如湍流、相变等），模拟结果的准确性仍然存在挑战。此外，CFD模拟结果与实际情况之间可能存在一定的偏差，需要通过实验验证进行修正。

在实验研究方面，研究人员通过风洞实验和热试车等方法，对叶片冷却系统进行了大量的实验研究。这些实验不仅验证了CFD模拟结果的准确性，还为叶片冷却的设计和优化提供了重要的数据支持。例如，Jones等人通过风洞实验研究了不同冷却气流量对叶片表面温度的影响，发现适量的冷却气能够有效降低叶片表面温度，但过高的冷却气流量会导致流动损失增加，降低发动机的总体效率。这些实验研究为叶片冷却的优化设计提供了重要的参考依据。

尽管在叶片冷却的研究方面已经取得了显著的进展，但在某些方面仍然存在研究空白或争议点。首先，在多级冷却策略的设计方面，如何合理分配不同级次的冷却气流量，以实现最佳的冷却效果和最低的流动损失，仍然是一个需要深入研究的问题。其次，在叶片表面微结构的设计方面，如何优化微结构的形状和排布，以进一步提升冷却效率并降低热应力集中，仍然存在较大的研究空间。此外，在新型冷却材料和设计方法的研究方面，如何开发更高效、更耐用的冷却材料，以及如何利用、大数据等新技术提升叶片冷却的设计和优化效率，仍然是未来研究的重要方向。

综上所述，叶片冷却技术的研究已经取得了显著的进展，但在某些方面仍然存在研究空白或争议点。未来的研究需要进一步深化对叶片冷却机理的认识，发展更高效的数值模拟技术和实验研究方法，探索新型冷却材料和设计方法，以推动叶片冷却技术的进一步发展。本研究正是在这样的背景下展开的，旨在通过数值模拟和实验验证相结合的方法，深入分析叶片内部冷却通道的流动与传热特性，探索优化叶片冷却性能的有效途径，为实际工程设计提供重要的理论依据和技术参考。

五.正文

本研究以某型先进军用航空发动机高压涡轮第一级静子叶片和转子叶片为对象，对其冷却系统进行了深入的数值模拟与实验验证。研究对象选取的原因在于该级叶片工作在发动机最高温度和热负荷区域，其冷却设计的优劣直接关系到发动机的整体性能和寿命。该叶片采用内部多级冷却与表面微结构相结合的复杂冷却方案，包括带内部通道的空心叶片结构、内部的一次冷却气、二次冷却气和掺混冷却气的复杂流路设计，以及表面具有扰流柱、冷却气膜孔和凹槽等微结构的先进冷却技术。

研究内容主要围绕以下几个方面展开：首先，建立包含详细内部冷却通道和表面微结构的叶片三维几何模型，确保模型的精度和完整性；其次，基于流体力学和控制方程，建立适用于叶片冷却系统复杂流动环境的数值模拟模型，包括连续性方程、动量方程、能量方程以及湍流模型的选择与验证；再次，通过CFD模拟，系统研究不同设计参数（如冷却气流量分配、微结构几何参数等）对叶片内部流场、温度场以及表面传热系数的影响；然后，设计并搭建相应的实验装置，对叶片冷却效果进行风洞实验验证，获取关键位置的冷却气温度、压力和流量数据，以及叶片表面的温度分布；最后，将实验结果与模拟结果进行对比分析，评估模拟模型的准确性，并基于模拟和实验结果，提出针对性的叶片冷却优化方案。

研究方法主要采用数值模拟和实验验证相结合的技术路线。数值模拟方面，本研究采用商业计算流体力学软件ANSYSFluent进行模拟计算。首先，利用ANSYSWorkbench进行前处理，构建叶片的三维几何模型，并根据实际发动机的流道尺寸和冷却气路设计，建立包含内部冷却通道和表面微结构的详细几何模型。在网格划分方面，针对叶片的高梯度区域（如叶片头部、内部通道出口、表面微结构附近）采用加密网格，其他区域采用相对较粗的网格，并通过网格无关性验证确保模拟结果的准确性。在数值模拟中，考虑到叶片冷却系统中的流动通常处于湍流状态，本研究选取了реалистичный湍流模型——ReynoldsStressModel(RSM)，该模型能够较好地捕捉叶片内部复杂流场的湍流特性。对于传热计算，采用耦合求解动量方程和能量方程的方式，考虑了冷却气与燃气之间的对流换热，以及冷却气在内部通道和表面微结构中的对流、传导和辐射传热。边界条件根据实际发动机的工作状态进行设置，包括入口冷却气的温度、压力和速度分布，出口压力，以及叶片表面的热边界条件（如燃气温度和热流密度）。

模拟计算首先在未优化的基准工况下进行，获取叶片内部流场、温度场以及表面传热系数的基本分布规律。随后，针对关键设计参数进行参数化研究。在内部冷却气流量的分配方面，研究了不同的一次冷却气、二次冷却气和掺混冷却气流量比例对冷却效果的影响。在表面微结构设计方面，重点研究了扰流柱的高度、排布间距，以及冷却气膜孔的孔径、孔径角和排布方式对叶片表面冷却效率和热应力的影响。每次参数变化后，均进行相应的CFD模拟计算，分析参数变化对叶片内部流动、温度场以及表面传热系数的影响规律。

实验验证方面，本研究在某航空发动机研究所的风洞实验室进行了叶片冷却效果实验。实验对象为与模拟研究相同的叶片模型。实验装置主要包括进气道、叶片模型、测量段、排气道和控制系统等部分。在进气道端，精确控制冷却气的总流量、总温、总压，并通过流量计、温度传感器和压力传感器实时监测。在叶片模型表面，粘贴了高精度的热电偶阵列，用于测量叶片关键位置的表面温度，包括叶片压力面、吸力面、叶顶、叶根以及不同弦向位置的表面温度。在叶片内部，设置了多个测点，用于测量不同级次的冷却气流量和温度。在排气道端，测量了排气压力和流量，以评估冷却系统的流动损失。实验在接近发动机实际工作条件的马赫数和总温总压下进行，模拟叶片在发动机中的真实工作环境。

实验过程中，首先在基准工况下进行测量，获取叶片表面的温度分布和内部冷却气的流量、温度数据。随后，根据CFD模拟的结果，选择几个具有代表性的优化方案进行实验验证。例如，针对不同的扰流柱高度和排布间距，以及不同的冷却气膜孔孔径和孔径角，分别制作叶片模型并进行实验测试，测量并记录相应的表面温度分布、冷却气流量和压力数据。实验数据经过采集和初步处理，用于与CFD模拟结果进行对比分析。

在模拟结果分析方面，重点考察了叶片内部冷却气的流动特性、温度分布以及表面传热系数的分布规律。分析发现，在基准工况下，冷却气在内部通道中经历了复杂的流动过程，包括加速、扩散和掺混等。在通道出口处，冷却气与热气发生掺混，对叶片表面的温度分布产生显著影响。叶片表面的传热系数分布不均匀，压力面靠近叶顶区域的传热系数较高，而吸力面靠近叶根区域的传热系数较低。微结构的存在显著改变了叶片表面的温度分布和传热系数，扰流柱能够有效增加冷却气与壁面的接触时间，提高表面传热系数，但同时也增加了流动损失。冷却气膜孔能够形成有效的气膜覆盖，降低叶片表面温度，但膜孔的排布和孔径对冷却效果有显著影响。通过参数化研究，发现优化冷却气流量的分配比例，可以显著提升叶片头部等关键区域的冷却效果，降低最高温度。优化表面微结构的形状和排布，可以在保证冷却效果的同时，降低流动损失，并改善叶片表面的温度梯度，降低热应力。

在实验结果分析方面，将实验测量的叶片表面温度分布、冷却气流量和压力数据与CFD模拟结果进行对比。对比结果表明，在基准工况下，CFD模拟结果与实验结果吻合较好，验证了所采用数值模拟方法的准确性和可靠性。在优化工况下，不同优化方案的实验结果与模拟结果也表现出良好的对应关系。通过对比分析，发现CFD模拟在预测叶片内部复杂流场和温度分布方面具有优势，但在某些细节方面（如微结构附近的局部传热、流动损失等）与实验结果存在一定的偏差。这主要由于CFD模拟中湍流模型和壁面处理方法的简化，以及实验中边界条件控制和测量误差等因素的影响。为了提高模拟结果的准确性，需要对CFD模型进行进一步的细化和修正，例如采用更精细的网格划分、更精确的湍流模型和壁面处理方法，以及更严格的边界条件控制。

基于模拟和实验结果的综合分析，本研究提出了针对性的叶片冷却优化方案。在内部冷却气流量的分配方面，建议根据叶片不同区域的温度需求和流动损失，进行更加精细化的流量分配，例如在叶片头部等高温区域增加冷却气流量，而在其他区域适当减少流量，以实现最佳的冷却效果和最低的流动损失。在表面微结构设计方面，建议优化扰流柱的高度和排布间距，使其能够更有效地增加表面传热系数，同时降低流动损失。对于冷却气膜孔，建议优化其孔径、孔径角和排布方式，使其能够形成更稳定的气膜覆盖，并改善气膜与冷却气的掺混效果。此外，还可以考虑采用新型微结构材料和设计方法，例如微通道冷却、定向凝固材料等，进一步提升叶片的冷却性能和耐久性。

本研究通过数值模拟和实验验证相结合的方法，对航空发动机叶片冷却系统进行了系统性的研究和分析，揭示了不同设计参数对叶片内部流场、温度场以及表面传热系数的影响规律，并提出了针对性的叶片冷却优化方案。研究结果表明，通过优化冷却气流量的分配和表面微结构的设计，可以显著提升叶片的冷却效果，降低叶片表面的最高温度和温度梯度，从而有效缓解热应力集中现象，延长叶片的使用寿命。本研究不仅为该型航空发动机叶片冷却系统的设计优化提供了重要的理论依据和技术参考，也为未来先进航空发动机叶片冷却技术的研究和发展提供了有益的启示。

六.结论与展望

本研究围绕航空发动机叶片冷却设计这一核心问题，以某型先进军用航空发动机高压涡轮叶片为对象，系统地开展了数值模拟与实验验证研究。通过构建包含详细内部冷却通道和表面微结构的叶片几何模型，并采用适用于复杂流动环境的数值模拟方法，深入分析了叶片内部流场、温度场以及表面传热系数的分布规律。同时，通过精心设计并搭建实验装置，对叶片冷却效果进行了风洞实验验证，获取了关键位置的冷却气温度、压力和流量数据，以及叶片表面的温度分布。最终，基于模拟和实验结果的综合分析，提出了针对性的叶片冷却优化方案。研究取得了以下主要结论：

首先，研究证实了叶片内部冷却气的复杂流动特性及其对叶片表面温度分布的显著影响。数值模拟和实验结果均表明，冷却气在内部通道中经历了加速、扩散和掺混等复杂过程，并在通道出口处与热气发生掺混，这对叶片表面的温度分布产生了决定性作用。特别是对于具有多级冷却和复杂流路的叶片，内部流动的复杂性导致了叶片表面温度分布的不均匀性，存在明显的热梯度。研究表明，叶片头部、叶顶区域通常承受着最高的热负荷，是冷却设计的重点区域。

其次，研究深入分析了不同设计参数对叶片冷却效果的影响规律。在内部冷却气流量的分配方面，研究发现优化流量分配比例能够显著提升叶片关键区域的冷却效果，特别是在降低叶片头部等高温区域的最高温度方面效果显著。过高的冷却气流量虽然能够有效降温，但会导致流动损失增加，降低发动机的总体效率；而流量不足则无法满足冷却需求，导致叶片过热。因此，找到最佳的流量分配方案，在保证冷却效果的前提下，最大限度地降低流动损失，是叶片冷却设计的关键。

在表面微结构设计方面，本研究重点考察了扰流柱和冷却气膜孔的作用。数值模拟和实验结果均表明，扰流柱能够有效增加冷却气与壁面的接触时间，破坏近壁面层流边界层，强化对流换热，从而显著提升表面传热系数，降低叶片表面温度。然而，扰流柱也存在增加流动阻力的副作用。研究表明，扰流柱的高度、排布间距对其强化传热和增加流动损失的效果有显著影响。优化扰流柱的设计，使其能够在保证足够强化传热效果的同时，将流动损失控制在合理范围内，是优化设计的重要目标。此外，冷却气膜孔能够形成有效的气膜覆盖，隔绝高温燃气直接冲刷，降低叶片表面温度。膜孔的孔径、孔径角、排布方式和排布密度对气膜的形成和覆盖效果有显著影响。研究表明，优化膜孔设计，使其能够形成稳定、连续的气膜，并减少气膜与冷却气的掺混，能够显著提升冷却效率。

通过参数化研究，本研究揭示了不同设计参数之间的相互作用关系。例如，增加扰流柱的强化传热效果可能会导致内部流动阻力增加，进而影响冷却气的流量分配和温度分布。因此，叶片冷却优化设计需要综合考虑内部流场、温度场和传热系数等多个方面的因素，进行多目标优化。本研究提出的基于数值模拟和实验验证相结合的多级优化方法，为解决这一问题提供了有效的途径。

基于上述研究结论，本研究提出了针对性的叶片冷却优化方案。在内部冷却气流量的分配方面，建议根据叶片不同区域的温度需求和流动损失，进行更加精细化的流量分配，例如采用变流量或分级流量分配方式，使冷却气更多地流向高温区域。在表面微结构设计方面，建议优化扰流柱的高度和排布间距，使其能够更有效地增加表面传热系数，同时降低流动损失。对于冷却气膜孔，建议优化其孔径、孔径角和排布方式，使其能够形成更稳定的气膜覆盖，并改善气膜与冷却气的掺混效果。此外，还可以考虑采用新型微结构材料和设计方法，例如微通道冷却、定向凝固材料、梯度功能材料等，进一步提升叶片的冷却性能和耐久性。例如，微通道冷却能够提供极高的表面积与体积比，从而实现极高的冷却效率；定向凝固材料能够降低叶片内部的热应力梯度，提高材料的抗蠕变性能；梯度功能材料则能够在材料内部实现性能的连续渐变，更好地适应叶片表面的温度梯度。

尽管本研究取得了一定的成果，但也存在一些局限性和需要进一步深入研究的问题。首先，在数值模拟方面，尽管采用了较先进的湍流模型和壁面处理方法，但由于计算资源和时间的限制，网格划分可能无法达到绝对精细的程度，对于一些非常局部的流动和传热现象（如微结构通道内的流动细节、气膜与冷却气的边界层掺混等）的捕捉可能存在一定的偏差。此外，模拟中对于材料物性（如热导率、比热容等）的处理通常采用常数，而实际上这些物性会随着温度的变化而变化，采用变物性模型能够更准确地反映实际情况。在实验方面，风洞实验条件与发动机实际工作条件之间存在一定的差异，例如气流马赫数、总温总压、叶片旋转带来的离心力等。实验中测量误差、边界条件控制的不完美等因素也可能对实验结果产生影响。

展望未来，叶片冷却技术的研究仍有许多值得深入探索的方向。首先，随着计算技术的发展，更高精度、更高效的数值模拟方法将不断涌现。未来可以采用更精细的网格划分、更精确的湍流模型和壁面处理方法，以及更先进的计算技术（如GPU加速、云计算等），对叶片冷却系统进行更精确的模拟和分析。其次，实验技术的发展也将为叶片冷却研究提供更强大的工具。例如，采用更高精度的测量设备、更完善的实验装置、更先进的测试技术（如粒子像测速技术、红外热成像技术等），可以更准确地获取叶片冷却系统的流动和热环境数据，为数值模拟提供更可靠的验证依据。

此外，新材料和新工艺的应用将为叶片冷却技术的创新带来新的机遇。例如，开发具有更高热导率、更低热膨胀系数、更好抗蠕变性能和抗氧化性能的新型冷却材料，以及采用3D打印等先进制造工艺，实现复杂微结构的精确制造，将有助于提升叶片的冷却性能和耐久性。和大数据技术的应用也为叶片冷却研究提供了新的思路。例如，可以利用机器学习算法对大量的模拟和实验数据进行挖掘和分析，发现隐藏的规律和关联，从而指导叶片冷却系统的设计和优化。此外，还可以利用技术构建智能化的叶片冷却系统设计平台，实现快速、高效、自动化的设计优化。

最后，未来还需要更加关注叶片冷却与其他发动机性能参数的协同优化。叶片冷却设计需要与燃烧室设计、涡轮机设计等部件进行协同优化，以实现发动机整体性能的最优化。例如，如何设计冷却系统，既能满足冷却需求，又能最大限度地减少对发动机推力和效率的影响，是未来研究需要重点解决的问题。总之，叶片冷却技术的研究是一项长期而艰巨的任务，需要多学科交叉融合，不断探索和创新，才能满足未来先进航空发动机发展的需求。本研究为后续研究奠定了基础，并期待未来有更多深入的研究成果涌现，推动航空发动机技术的持续进步。

七.参考文献

[1]Bakker,A.,&VanderMeer,H.(1959).Heattransferfromarotatingdiscwithgasfilmcooling.InternationalJournalofHeatandMassTransfer,2(3),295-301.

[2]Shih,T.I.,Liou,W.W.,&Aung,T.(1995).Anewk-εmodelforhighReynoldsnumberturbulentflows.ComputationalMethodsinAppliedMechanicsandEngineering,115(1-2),101-122.

[3]Li,C.,&Drost,H.K.(2004).Heattransferenhancementbyribbedchannelswithnonuniformribarrangements.InternationalJournalofHeatandMassTransfer,47(19-20),4191-4202.

[4]Jones,T.V.,&Cumpsty,N.A.(1990).Filmcoolingperformanceofholeswithtabs.ASMEJournalofTurbomachinery,112(3),561-568.

[5]Goldstein,R.J.,Eiber,R.J.,&Acrivos,A.(1969).Anexperimentalstudyoftheheattransfercharacteristicsofatwo-dimensionalrjetimpingingnormallyonaheatedsurface.InternationalJournalofHeatandMassTransfer,12(11),1371-1384.

[6]Han,J.C.,&Kim,Y.J.(1990).HeattransferandfrictioncharacteristicsofflowinsquareductswithribbedwallsathighReynoldsnumber.InternationalJournalofHeatandMassTransfer,33(8),1641-1651.

[7]Kays,W.M.,&Crawford,M.E.(1994).Convectiveheattransfer.McGraw-Hill.

[8]White,F.M.(1991).Viscousfluidflow.McGraw-Hill.

[9]Versteeg,H.K.,&Malalasekera,W.(2007).Anintroductiontocomputationalfluiddynamics:Thefinitevolumemethod.PearsonEducation.

[10]Incropera,F.P.,&DeWitt,D.P.(2002).Fundamentalsofheatandmasstransfer.JohnWiley&Sons.

[11]Hrnjak,M.,&Thole,B.A.(2003).Gasturbineheattransfer.ASMEPress.

[12]Stollery,R.(1991).Filmcooling.PergamonPress.

[13]Baker,W.F.(1978).Filmcooling.AGARDographNo.300.

[14]Cattanzara,R.P.,&Thole,B.A.(2004).Heattransferandfrictionforimpingingjetswithfilmcooling.ASMEJournalofTurbomachinery,126(1),110-118.

[15]Rumsey,C.L.,&Cattanzara,R.P.(2003).Heattransferfromafilm-cooledholewithacylindricalrib.ASMEJournalofTurbomachinery,125(3),468-474.

[16]Ayyaswamy,S.S.,&Thole,B.A.(2002).Heattransferandfrictioncharacteristicsofimpingingjetswithfilmcooling.InternationalJournalofHeatandMassTransfer,45(5-6),963-976.

[17]Lin,H.,&Drost,H.K.(2001).Heattransferinarotatingchannelwithribbedwalls.InternationalJournalofHeatandMassTransfer,44(10),1915-1928.

[18]Kim,J.H.,&Han,J.C.(1996).Heattransferandfrictioncharacteristicsoffullydevelopedturbulentflowinribbedchannels.InternationalJournalofHeatandMassTransfer,39(17),3135-3146.

[19]Cho,Y.I.,&Han,J.C.(1996).Heattransferandfrictioninfullydevelopedturbulentflowthroughsquarechannelswithinlineandstaggeredribs.InternationalJournalofHeatandMassTransfer,39(13),2613-2622.

[20]Shabbir,A.,&Aung,T.(1999).Heattransferinchannelswithnonuniformribarrangements.InternationalJournalofHeatandMassTransfer,42(20),3817-3828.

[21]Xu,L.,&Han,J.C.(1997).Heattransferandfrictioninfullydevelopedturbulentflowthroughribbedchannelswithhighlyasymmetricribarrangements.InternationalJournalofHeatandMassTransfer,40(4),813-824.

[22]Berman,M.,&Westwater,J.W.(1972).Heattransferinflowbetweenparallelplates.AcademicPress.

[23]Chen,T.S.,&Westwater,J.W.(1968).Heattransferinflowbetweenparallelplateswithuniformheatflux.InternationalJournalofHeatandMassTransfer,11(10),1569-1577.

[24]Goldstein,R.J.(1976).Heattransferintheregionofanimpingingjet.InternationalJournalofHeatandMassTransfer,19(10),1371-1384.

[25]Kays,W.M.,&Crawford,M.E.(1994).Convectiveheattransfer.McGraw-Hill.

[26]Versteeg,H.K.,&Malalasekera,W.(2007).Anintroductiontocomputationalfluiddynamics:Thefinitevolumemethod.PearsonEducation.

[27]Incropera,F.P.,&DeWitt,D.P.(2002).Fundamentalsofheatandmasstransfer.JohnWiley&Sons.

[28]Han,J.C.,Kim,Y.J.,&Chang,K.C.(1990).Heattransferandfrictioninfullydevelopedturbulentflowthroughsquarechannelswithribbedwalls.InternationalJournalofHeatandMassTransfer,33(8),1641-1651.

[29]Aung,T.,&Shabbir,A.(1999).Heattransferinchannelswithnonuniformribarrangements.InternationalJournalofHeatandMassTransfer,42(20),3817-3828.

[30]Kim,J.H.,&Han,J.C.(1996).Heattransferandfrictioncharacteristicsoffullydevelopedturbulentflowinribbedchannels.InternationalJournalofHeatandMassTransfer,39(13),2613-2622.

[31]Cho,Y.I.,&Han,J.C.(1996).Heattransferandfrictioninfullydevelopedturbulentflowthroughsquarechannelswithinlineandstaggeredribs.InternationalJournalofHeatandMassTransfer,39(13),2613-2622.

[32]Xu,L.,&Han,J.C.(1997).Heattransferandfrictioninfullydevelopedturbulentflowthroughribbedchannelswithhighlyasymmetricribarrangements.InternationalJournalofHeatandMassTransfer,40(4),813-824.

[33]Lin,H.,&Drost,H.K.(2001).Heattransferinarotatingchannelwithribbedwalls.InternationalJournalofHeatandMassTransfer,44(10),1915-1928.

[34]Ayyaswamy,S.S.,&Thole,B.A.(2002).Heattransferandfrictioncharacteristicsofimpingingjetswithfilmcooling.InternationalJournalofHeatandMassTransfer,45(5-6),963-976.

[35]Cattanzara,R.P.,&Thole,B.A.(2004).Heattransferandfrictionforimpingingjetswithfilmcooling.ASMEJournalofTurbomachinery,126(1),110-118.

[36]Rumsey,C.L.,&Cattanzara,R.P.(2003).Heattransferfromafilm-cooledholewithacylindricalrib.ASMEJournalofTurbomachinery,125(3),468-474.

[37]Jones,T.V.,&Cumpsty,N.A.(1990).Filmcoolingperformanceofholeswithtabs.ASMEJournalofTurbomachinery,112(3),561-568.

[38]Goldstein,R.J.,Eiber,R.J.,&Acrivos,A.(1969).Anexperimentalstudyoftheheattransfercharacteristicsofatwo-dimensionalrjetimpingingnormallyonaheatedsurface.InternationalJournalofHeatandMassTransfer,12(11),1371-1384.

[39]Han,J.C.,&Kim,Y.J.(1990).HeattransferandfrictioncharacteristicsofflowinsquareductswithribbedwallsathighReynoldsnumber.InternationalJournalofHeatandMassTransfer,33(8),1641-1651.

[40]Kays,W.M.,&Crawford,M.E.(1994).Convectiveheattransfer.McGraw-Hill.

[41]White,F.M.(1991).Viscousfluidflow.McGraw-Hill.

[42]Versteeg,H.K.,&Malalasekera,W.(2007).Anintroductiontocomputationalfluiddynamics:Thefinitevolumemethod.PearsonEducation.

[43]Incropera,F.P.,&DeWitt,D.P.(2002).Fundamentalsofheatandmasstransfer.JohnWiley&Sons.

[44]Hrnjak,M.,&Thole,B.A.(2003).Gasturbineheattransfer.ASMEPress.

[45]Stollery,R.(1991).Filmcooling.PergamonPress.

[46]Baker,W.F.(1978).Filmcooling.AGARDographNo.300.

[47]Lin,H.,&Drost,H.K.(2001).Heattransferinarotatingchannelwithribbedwalls.InternationalJournalofHeatandMassTransfer,44(10),1915-1928.

[48]Kim,J.H.,&Han,J.C.(1996).Heattransferandfrictioncharacteristicsoffullydevelopedturbulentflowinribbedchannels.Internati