高超声速气动热数值模拟：方法、挑战与应用进展

高超声速气动热数值模拟：方法、挑战与应用进展一、绪论1.1研究背景与意义随着航空航天技术的飞速发展，高超声速飞行器已成为国际航空领域的研究热点和前沿。高超声速飞行器，通常指飞行速度大于5马赫（即5倍音速）的飞行器，因其具备高速、高效和高机动性等显著优势，在军事和民用领域均展现出了巨大的应用潜力。在军事领域，高超声速飞行器能够突破现有防空系统的拦截，实现对敌方目标的快速精确打击，显著提升作战效能和战略威慑力；在民用领域，高超声速飞行器可大幅缩短长途旅行时间，实现全球快速运输，推动航空运输业的变革性发展，同时在太空探索、气象监测等领域也具有重要应用价值。然而，高超声速飞行器在高速飞行过程中，会与周围空气发生剧烈的相互作用，产生极为严重的气动热现象。当飞行器以高超声速飞行时，其头部、机翼前缘、机身表面等部位的空气会被急剧压缩，形成强烈的激波，导致气体温度急剧升高。同时，飞行器表面与高速气流之间的摩擦也会进一步加剧热量的产生，使得飞行器表面承受极高的热流密度和温度。这种极端的气动热环境会对飞行器的结构完整性、材料性能和飞行性能产生严重威胁。例如，过高的温度可能导致飞行器结构材料软化、变形甚至熔化，降低结构强度，危及飞行安全；还可能影响飞行器上电子设备、推进系统等关键部件的正常工作，进而影响飞行器的整体性能和任务完成能力。因此，准确预测和有效应对气动热问题，成为高超声速飞行器设计和发展面临的关键挑战之一。为了解决高超声速飞行器的气动热问题，气动热数值模拟技术应运而生，并发挥着至关重要的作用。气动热数值模拟是基于计算流体力学（CFD）、传热学、热力学等多学科理论，通过数值计算方法对高超声速飞行器在飞行过程中产生的气动热现象进行模拟和分析。它能够提供飞行器表面的热流密度分布、温度场分布等详细信息，帮助研究人员深入理解气动热的产生机制和演化规律。与传统的实验研究方法相比，气动热数值模拟具有诸多优势。一方面，数值模拟不受实验条件的限制，可以方便地对不同飞行工况、不同飞行器外形和不同热防护方案进行研究，大大提高了研究效率和灵活性；另一方面，数值模拟可以获取实验难以测量的流场细节信息，为理论分析和工程设计提供有力支持。同时，数值模拟还可以在飞行器设计的早期阶段进行，通过对不同设计方案的气动热性能进行预测和评估，为优化设计提供依据，从而减少实验次数，降低研发成本。在高超声速飞行器的设计过程中，气动热数值模拟结果对于热防护系统的设计和优化具有关键指导作用。热防护系统是高超声速飞行器抵御气动热的关键部件，其性能直接关系到飞行器的安全和可靠性。通过气动热数值模拟，能够准确预测飞行器表面的热载荷，进而根据热载荷分布情况合理选择热防护材料、设计热防护结构，确保热防护系统能够有效地保护飞行器结构免受高温的破坏。例如，根据数值模拟得到的热流密度和温度分布，确定在高温区域采用耐高温性能更好的材料，在热流密度较大的部位优化热防护结构形式，以提高热防护效率。此外，数值模拟还可以对热防护系统的性能进行评估和验证，通过对比模拟结果和实际需求，及时发现热防护系统设计中存在的问题并进行改进，确保热防护系统在各种飞行条件下都能可靠工作。综上所述，开展高超声速气动热数值模拟研究，对于深入理解高超声速飞行器的气动热物理机制，提高飞行器的设计水平和热防护能力，推动高超声速技术的发展和应用具有重要的理论意义和实际工程价值。它不仅有助于解决高超声速飞行器研发过程中的关键技术难题，提升我国在航空航天领域的核心竞争力，还将为未来高超声速飞行器在军事和民用领域的广泛应用奠定坚实的技术基础。1.2国内外研究现状在高超声速气动热数值模拟领域，国外起步较早，取得了一系列具有开创性的研究成果。美国作为该领域的先驱者，凭借其强大的科研实力和雄厚的资金支持，开展了众多具有深远影响力的研究项目。在20世纪60年代的“阿波罗”计划中，美国便运用数值模拟方法对返回舱再入过程中的气动热进行了研究，通过不断改进数值算法和模型，成功预测了返回舱表面的热流分布和温度变化，为返回舱的热防护系统设计提供了关键依据。此后，美国国家航空航天局（NASA）持续加大对高超声速气动热研究的投入，先后实施了一系列高超声速飞行器研发项目，如X-43、X-51等。在这些项目中，数值模拟技术得到了广泛而深入的应用。科研人员采用高精度的数值方法对飞行器复杂外形下的三维流场进行模拟，精确捕捉激波与边界层的相互作用，深入研究气动热的产生和传递机制。同时，结合先进的实验技术对模拟结果进行验证和校准，不断完善数值模拟方法和模型，显著提高了气动热数值模拟的准确性和可靠性。欧洲各国在高超声速气动热数值模拟方面也展现出了卓越的研究实力。法国国家航空航天研究院（ONERA）长期致力于高超声速领域的研究，在数值模拟方法和应用研究方面取得了丰硕成果。ONERA研发了一系列先进的计算流体力学（CFD）软件，如elsA等，这些软件具备强大的数值计算能力和良好的工程适用性，能够对高超声速飞行器的复杂流场进行高效模拟。通过对不同类型高超声速飞行器的数值模拟研究，ONERA深入分析了气动热的影响因素和变化规律，为欧洲高超声速飞行器的设计和发展提供了重要的技术支持。德国在高超声速气动热研究方面也有着深厚的技术积累，其在数值算法、网格生成技术和多物理场耦合模拟等方面开展了大量前沿研究工作。德国航空航天中心（DLR）研发的TAU软件在高超声速流场模拟中表现出色，能够精确处理复杂的几何外形和物理现象，为德国及欧洲其他国家的高超声速项目提供了重要的数值模拟工具。俄罗斯在高超声速领域同样有着深厚的技术底蕴和丰富的研究经验。俄罗斯的科研机构和高校在高超声速气动热数值模拟方面开展了系统而深入的研究。他们注重理论研究与工程实践相结合，在数值方法的理论基础和应用技术方面取得了多项突破。例如，在求解高超声速流动的控制方程时，俄罗斯学者提出了一系列高效的数值算法，能够有效提高计算效率和精度。同时，俄罗斯在高超声速飞行器的热防护材料和结构设计方面也取得了显著成果，通过数值模拟与实验研究相结合的方式，优化热防护系统的设计，提高飞行器在极端气动热环境下的生存能力。近年来，国内在高超声速气动热数值模拟领域也取得了长足的进步。随着我国航空航天事业的快速发展，对高超声速技术的研究需求日益迫切，国内众多科研机构和高校纷纷加大在该领域的研究投入，取得了一系列具有国际影响力的研究成果。中国空气动力研究与发展中心作为我国空气动力学研究的核心机构，在高超声速气动热数值模拟方面发挥了重要的引领作用。该中心自主研发了多款具有自主知识产权的CFD软件，如FDS等。这些软件针对高超声速流动的特点进行了优化设计，能够准确模拟高超声速飞行器复杂外形下的流场特性和气动热现象。通过对多种高超声速飞行器模型的数值模拟研究，该中心深入分析了不同飞行条件下的气动热分布规律，为我国高超声速飞行器的设计和热防护系统的优化提供了关键的技术支持。北京航空航天大学、南京航空航天大学等高校在高超声速气动热数值模拟领域也开展了大量富有成效的研究工作。这些高校依托其在航空航天领域的学科优势和人才资源，组建了专业的研究团队，深入开展高超声速气动热的基础理论和应用技术研究。在数值模拟方法方面，高校研究团队不断探索创新，提出了一系列针对高超声速流动的高效数值算法，如高精度的有限体积法、有限差分法和谱方法等。这些算法在提高计算精度和效率方面取得了显著成效，能够更准确地模拟高超声速飞行器的气动热现象。同时，高校研究团队还注重多学科交叉融合，将数值模拟与实验研究、材料科学、结构力学等学科相结合，开展高超声速飞行器热防护系统的一体化设计研究。通过数值模拟与实验验证的相互迭代，不断优化热防护系统的设计方案，提高其性能和可靠性。在数值模拟技术的应用方面，国内科研人员已将其广泛应用于我国多种高超声速飞行器的研制过程中。例如，在我国新型高超声速导弹和航天飞行器的设计中，数值模拟技术发挥了至关重要的作用。通过数值模拟，科研人员能够在设计阶段对飞行器的气动热性能进行准确预测和评估，为飞行器的外形优化、热防护系统设计提供科学依据。同时，数值模拟还可以用于分析不同热防护材料和结构的性能，指导热防护系统的选材和结构优化，提高热防护系统的防护效果和可靠性。通过数值模拟与实验研究的紧密结合，我国在高超声速飞行器的研制过程中，有效缩短了研发周期，降低了研发成本，提高了研发效率和产品质量。1.3研究目的与内容本研究旨在通过对高超声速气动热数值模拟方法的深入研究，改进和完善现有的数值模拟技术，提高模拟的准确性和效率，为高超声速飞行器的热防护系统设计和优化提供更加可靠的依据。具体而言，本研究拟攻克高超声速气动热数值模拟中的关键难题，如高精度数值算法的开发、复杂流场的准确模拟、多物理场耦合效应的处理等，以提升数值模拟的精度和可靠性。此外，本研究还致力于拓展高超声速气动热数值模拟的应用范围，将其应用于不同类型高超声速飞行器的设计和分析中，为实际工程应用提供有力支持。为实现上述研究目的，本研究将围绕以下几个方面展开具体内容的研究：高超声速气动热物理机制研究：深入剖析高超声速飞行器在飞行过程中产生气动热的物理机制，包括激波与边界层的相互作用、粘性耗散、化学反应等因素对气动热的影响。通过理论分析和数值模拟相结合的方法，建立准确的物理模型，揭示气动热的产生和演化规律。数值模拟方法研究：对现有的数值模拟方法进行系统研究和比较，如有限体积法、有限差分法、谱方法等。针对高超声速气动热数值模拟的特点和需求，优化和改进数值算法，提高计算精度和效率。研究数值模拟中的关键技术问题，如网格生成、边界条件处理、湍流模型选择等，提出有效的解决方案。多物理场耦合效应研究：考虑高超声速飞行器在飞行过程中涉及的多种物理场耦合效应，如流场与温度场、结构场的耦合等。建立多物理场耦合的数学模型和数值求解方法，研究耦合效应对气动热和飞行器性能的影响。数值模拟软件研发与验证：基于上述研究成果，研发具有自主知识产权的高超声速气动热数值模拟软件。对软件进行功能测试和验证，通过与实验数据和已有研究成果的对比分析，评估软件的准确性和可靠性。工程应用研究：将研发的数值模拟软件应用于实际高超声速飞行器的设计和分析中，如对飞行器的不同外形设计方案进行气动热性能评估，为外形优化提供依据。研究不同热防护材料和结构在高超声速气动热环境下的性能表现，为热防护系统的设计和优化提供指导。二、高超声速气动热基础理论2.1物理机制2.1.1激波与边界层相互作用当高超声速飞行器在大气层中飞行时，其飞行速度远高于当地声速，飞行器前方的空气无法及时避让，被急剧压缩，从而形成激波。激波是一种强间断面，在激波前后，气体的压力、密度、温度等参数会发生剧烈的突变。以正激波为例，气体经过正激波后，压力会急剧升高数倍甚至数十倍，密度也会显著增大，而速度则会大幅降低。这种参数的突变导致气体的内能急剧增加，进而使温度迅速升高。例如，当飞行器以马赫数8飞行时，正激波后的气体温度可升高至数千开尔文。边界层是指在固体表面附近，由于流体粘性的作用，流速从壁面处的零逐渐变化到主流速度的薄层。在高超声速流动中，边界层内的气体与飞行器表面发生强烈的摩擦，产生大量的热量。同时，激波与边界层之间存在着复杂的相互作用，这种相互作用会进一步加剧气动热现象。当激波与边界层相互作用时，会导致边界层分离、再附等复杂流动现象的出现。边界层分离会使气流在分离区域形成漩涡，增加了流动的复杂性和能量损失，从而产生更多的热量。而边界层再附则会使气流在再附点处产生强烈的冲击，导致局部热流密度急剧升高。例如，在高超声速飞行器的机翼前缘，激波与边界层的相互作用常常导致该部位的热流密度远高于其他部位，对机翼的热防护提出了极高的要求。2.1.2粘性耗散与热传导粘性耗散是指在流体流动过程中，由于流体粘性的存在，流体内部各层之间存在相对运动，从而产生内摩擦力，这种内摩擦力做功将机械能转化为热能的现象。在高超声速流动中，气体的粘性作用不可忽略，粘性耗散产生的热能对气动热有重要贡献。当高超声速气流流过飞行器表面时，边界层内的气体由于粘性作用，速度逐渐降低，其动能不断转化为热能，使得边界层内气体的温度升高。粘性耗散产生的热能会使边界层内气体的温度分布发生变化，靠近壁面处的气体温度升高更为显著，形成较大的温度梯度。热传导是指由于温度差的存在，热量从高温区域向低温区域传递的现象。在高超声速气动热问题中，热传导在热量传递过程中起着重要的作用。一方面，边界层内由于粘性耗散产生的热能会通过热传导的方式向飞行器表面传递，使飞行器表面温度升高。热传导的速率与材料的热导率、温度梯度等因素有关。对于热导率较高的材料，热量能够更快速地传递到表面，导致表面温度升高更快。另一方面，飞行器表面吸收的热量也会通过热传导在结构内部进行传递，影响结构的温度分布。如果结构内部的热传导性能不佳，热量会在局部积聚，导致局部温度过高，可能对结构造成损害。例如，在高超声速飞行器的热防护系统中，热传导性能良好的隔热材料能够有效地阻止热量向结构内部传递，降低结构的温度，从而保护结构的安全。2.1.3真实气体效应在高超声速飞行条件下，气体所处的环境温度极高，气体分子的热运动非常剧烈，此时气体的性质会偏离理想气体的假设，呈现出真实气体效应。真实气体效应主要包括气体分子的振动激发、解离、电离以及化学反应等现象。当气体温度升高到一定程度时，气体分子的振动自由度被激发，分子内部的振动能量增加。以空气为例，当温度达到2500K左右时，氧气和氮气分子的振动自由度开始被激发。随着温度的进一步升高，气体分子会发生解离，即分子分解为原子。在2000K-4000K范围内，氧分子会逐渐解离，并与氮气发生化学反应。当温度超过4000K时，氮分子也会发生解离。当温度达到9000K以上时，氧、氮原子会发生电离，产生自由电子和离子，使气体变成导电的等离子体。这些真实气体效应会对气动热产生显著的影响。一方面，气体分子的振动激发、解离、电离等过程都需要吸收能量，这会导致气体内能的变化，进而影响气体的热力学性质和传热特性。例如，气体分子的解离会吸收大量的热量，使得激波后的气体温度升高幅度减小，从而在一定程度上缓解了气动热问题。另一方面，真实气体效应还会改变激波的位置与形状、边界层的结构和特性等流场参数。例如，气体电离产生的等离子体鞘套会对电磁波产生吸收和散射作用，影响飞行器的通信和雷达探测性能，同时也会改变飞行器表面的热流分布。此外，化学反应的发生会产生新的物质，改变气体的组成成分，进一步影响气体的热力学和输运性质，从而对气动热产生复杂的影响。2.2数学模型2.2.1流动控制方程在高超声速气动热数值模拟中，流动控制方程是描述流体运动的基本方程，其核心是纳维-斯托克斯（N-S）方程。N-S方程基于质量守恒、动量守恒和能量守恒定律，全面地描述了粘性流体的运动规律。在直角坐标系下，三维非定常N-S方程的守恒形式可表示为：\frac{\partial\boldsymbol{U}}{\partialt}+\frac{\partial\boldsymbol{E}_{i}}{\partialx_{i}}+\frac{\partial\boldsymbol{F}_{i}}{\partialy_{i}}+\frac{\partial\boldsymbol{G}_{i}}{\partialz_{i}}=\boldsymbol{H}其中，\boldsymbol{U}为守恒变量向量，\boldsymbol{U}=[\rho,\rhou,\rhov,\rhow,\rhoe]^{T}，\rho为流体密度，u、v、w分别为x、y、z方向的速度分量，e为单位质量流体的总能量。\boldsymbol{E}_{i}、\boldsymbol{F}_{i}、\boldsymbol{G}_{i}分别为无粘通量向量和粘性通量向量，\boldsymbol{H}为源项向量。无粘通量向量\boldsymbol{E}_{i}在x方向的表达式为\boldsymbol{E}_{i}=[\rhou,\rhou^{2}+p,\rhouv,\rhouw,(\rhoe+p)u]^{T}，其中p为压力。粘性通量向量\boldsymbol{F}_{i}在x方向的表达式较为复杂，涉及到流体的粘性系数和速度梯度等。源项向量\boldsymbol{H}则包含了各种体积力和其他物理过程的影响。在高超声速流动中，根据具体问题的特点和研究需求，N-S方程可以进行适当的简化。例如，在某些情况下，当流动满足定常条件时，时间导数项\frac{\partial\boldsymbol{U}}{\partialt}为零，方程可简化为定常N-S方程，从而大大减少计算量。在研究高超声速飞行器外部绕流问题时，若忽略飞行器表面的质量注入和化学反应等因素，源项向量\boldsymbol{H}中的相应项可设为零。当流动特征尺度较大且流体粘性相对较小时，可采用无粘的欧拉方程进行近似计算。欧拉方程是N-S方程忽略粘性项后的简化形式，其形式相对简单，计算效率较高，在一些对粘性效应要求不高的工程应用中具有广泛的应用。然而，在高超声速流动中，边界层内的粘性效应显著，此时必须使用完整的N-S方程才能准确描述流动特性。例如，在计算高超声速飞行器机翼表面的热流密度时，边界层内的粘性耗散和热传导对热流密度的分布有着重要影响，若使用欧拉方程则无法准确捕捉这些物理现象。N-S方程适用于各种粘性流体的流动问题，无论是低速流动还是高超声速流动，都能从理论上对其进行描述。在高超声速领域，N-S方程能够准确捕捉激波与边界层的相互作用、粘性耗散等复杂物理现象。例如，在模拟高超声速飞行器再入大气层的过程中，N-S方程可以精确计算激波的位置和强度，以及边界层内的速度、温度和压力分布，为研究气动热现象提供了重要的理论基础。然而，N-S方程是一组高度非线性的偏微分方程，其求解过程非常复杂，通常需要借助数值方法进行求解。随着计算机技术的不断发展，各种高效的数值算法应运而生，如有限体积法、有限差分法等，使得N-S方程在高超声速气动热数值模拟中的应用成为可能。2.2.2热传导方程热传导方程是描述热量在物体内部传递规律的重要方程，它基于傅里叶定律，揭示了温度场与热流密度之间的关系。在各向同性的介质中，三维非稳态热传导方程的一般形式为：\rhoc_{p}\frac{\partialT}{\partialt}=\nabla\cdot(k\nablaT)+q_{v}其中，\rho为物体的密度，c_{p}为定压比热容，T为温度，t为时间，k为热导率，q_{v}为内热源强度。方程左边\rhoc_{p}\frac{\partialT}{\partialt}表示单位体积物体内由于温度随时间变化而储存的热量；右边第一项\nabla\cdot(k\nablaT)表示通过热传导方式传入单位体积物体内的热量，其中\nabla为梯度算子，k\nablaT表示热流密度向量，其方向与温度梯度方向相反，大小与温度梯度成正比；右边第二项q_{v}表示单位体积内热源产生的热量。热传导方程在高超声速气动热问题中有着广泛的应用。在高超声速飞行器的热防护系统设计中，需要准确预测热量在热防护材料内部的传递过程，以评估热防护系统的性能。此时，热传导方程可以用于计算热防护材料内部的温度分布，为热防护材料的选择和结构设计提供依据。例如，对于采用陶瓷基复合材料作为热防护材料的高超声速飞行器，通过求解热传导方程，可以得到在不同飞行工况下，陶瓷基复合材料内部的温度随时间和空间的变化规律。根据这些温度分布信息，可以确定热防护材料的厚度是否满足要求，以及是否需要采取额外的隔热措施。在研究高超声速飞行器结构的热响应时，热传导方程也起着关键作用。飞行器结构在气动热的作用下，其内部的温度会发生变化，进而导致结构材料的力学性能改变。通过求解热传导方程，可以获得结构内部的温度场，再结合结构力学的相关理论，就可以分析结构在温度载荷作用下的应力、应变分布，评估结构的强度和稳定性。热传导方程的求解方法有多种，包括解析法、数值法等。解析法适用于一些简单的几何形状和边界条件的问题，通过数学推导可以得到温度分布的精确解析解。然而，在实际的高超声速气动热问题中，物体的几何形状往往非常复杂，边界条件也多种多样，解析法的应用受到很大限制。因此，数值法成为求解热传导方程的主要手段。常用的数值方法有有限差分法、有限元法和有限体积法等。有限差分法是将求解区域离散成网格，通过差商近似代替导数，将热传导方程转化为代数方程组进行求解；有限元法则是将求解区域划分为有限个单元，通过插值函数逼近温度场，再利用变分原理建立求解方程；有限体积法是基于控制体积的思想，将热传导方程在每个控制体积上进行积分，得到离散的代数方程。这些数值方法各有优缺点，在实际应用中需要根据具体问题的特点选择合适的方法。2.2.3热辐射模型在高超声速气动热环境中，当气体温度极高时，热辐射成为不可忽视的热量传递方式。热辐射是物体由于自身温度而向外发射电磁波的现象，它不需要任何介质即可在真空中传播。热辐射模型在高超声速气动热数值模拟中起着关键作用，用于考虑辐射换热对温度场和热流分布的影响。在高超声速流动中，气体分子的热运动剧烈，激发态分子通过辐射跃迁回到基态时会释放出光子，从而产生热辐射。同时，高温气体对外部辐射也具有吸收和散射作用，使得辐射换热过程变得复杂。常见的热辐射模型包括灰体辐射模型和谱带模型。灰体辐射模型假设物体的辐射特性与波长无关，发射率和吸收率为常数。在该模型中，辐射热流密度可根据斯蒂芬-玻尔兹曼定律计算：q_{r}=\varepsilon\sigma(T^{4}-T_{0}^{4})其中，q_{r}为辐射热流密度，\varepsilon为物体的发射率，\sigma为斯蒂芬-玻尔兹曼常数，T为物体表面温度，T_{0}为周围环境温度。灰体辐射模型适用于一些对辐射精度要求不高的情况，计算相对简单，能够快速估算辐射换热的大致量级。然而，在实际的高超声速气动热问题中，气体的辐射特性与波长密切相关，灰体辐射模型的准确性难以满足要求。此时，需要采用更精确的谱带模型。谱带模型考虑了气体辐射的波长选择性，将整个辐射光谱划分为多个谱带，分别计算每个谱带内的辐射特性。常用的谱带模型有窄谱带模型和宽谱带模型。窄谱带模型将光谱划分为非常窄的谱带，能够精确计算辐射特性，但计算量较大；宽谱带模型则对光谱进行更粗的划分，在一定程度上平衡了计算精度和效率。在模拟高超声速飞行器再入大气层时，飞行器周围的高温气体包含多种成分，如氧气、氮气等，这些气体在不同波长下的辐射特性差异显著。采用谱带模型可以更准确地考虑这些气体的辐射特性，从而得到更精确的辐射换热计算结果。热辐射模型在高超声速气动热数值模拟中的应用，能够显著提高模拟结果的准确性。通过考虑辐射换热，数值模拟可以更真实地反映高超声速飞行器在高温环境下的热状态，为热防护系统的设计和优化提供更可靠的依据。例如，在设计高超声速飞行器的热防护材料时，准确的热辐射模型可以帮助确定材料所需承受的辐射热流，从而选择合适的耐高温、抗辐射材料。在分析飞行器表面的温度分布时，考虑辐射换热可以更准确地预测高温区域的位置和温度值，为热防护结构的布局提供指导。三、数值模拟方法3.1常见数值方法在高超声速气动热数值模拟中，为了求解复杂的控制方程，需要采用合适的数值方法将连续的物理问题离散化，转化为可在计算机上求解的代数方程组。目前，常用的数值方法主要有有限体积法、有限差分法和谱方法，它们在离散控制方程、求解原理以及适用范围等方面各有特点。3.1.1有限体积法有限体积法（FiniteVolumeMethod，FVM）是一种基于控制体积积分的数值方法，其基本思想是将计算区域划分为一系列互不重叠的控制体积，把待解的偏微分方程对每个控制体积进行积分，从而得到一组离散方程。在有限体积法中，控制体积是应用控制方程和守恒定律的最小单位，通过对控制体积的积分，使得离散方程在每个控制体积内都满足守恒原理，这是有限体积法的一个重要优势。以二维可压缩N-S方程为例，在一般曲线坐标系下，其守恒形式为：\frac{1}{J}\frac{\partial\boldsymbol{Q}}{\partialt}+\frac{\partial(\widetilde{\boldsymbol{F}}-\widetilde{\boldsymbol{F}}_{v})}{\partial\xi}+\frac{\partial(\widetilde{\boldsymbol{G}}-\widetilde{\boldsymbol{G}}_{v})}{\partial\eta}=0其中，J为坐标变换雅克比行列式，\boldsymbol{Q}为守恒变量，\widetilde{\boldsymbol{F}}、\widetilde{\boldsymbol{G}}为无粘通量，\widetilde{\boldsymbol{F}}_{v}、\widetilde{\boldsymbol{G}}_{v}为粘性通量。对上述方程运用基于单元中心型的有限体积法进行空间离散，将计算区域划分为网格，每个网格单元的中心定义为一个节点。对于无粘通量，可采用合适的通量分裂格式，如M-AUSMPW+格式，通过MUSCL方法选用Vanleer平均限制器进行空间重构，以达到二阶精度；粘性通量则采用二阶中心格式离散。在有限体积法中，对控制方程在控制体积上积分时，需要将积分项转化为控制体积界面上的通量。以无粘通量为例，根据高斯散度定理，控制体积内的通量积分可转化为控制体积界面上的通量积分。假设控制体积的界面为S，则有：\int_{V}\frac{\partial\widetilde{\boldsymbol{F}}}{\partial\xi}dV=\oint_{S}\widetilde{\boldsymbol{F}}\cdot\boldsymbol{n}dS其中，\boldsymbol{n}为界面S的单位法向量。通过这种方式，将偏微分方程转化为关于控制体积界面通量的代数方程。在实际计算中，需要对界面通量进行近似计算，常用的方法有中心差分、迎风差分等。例如，在中心差分格式中，界面通量取相邻节点通量的平均值；迎风差分格式则根据流动方向，选择上游节点的通量来计算界面通量。迎风差分格式能够更好地捕捉激波等间断现象，在高超声速流动模拟中得到广泛应用。有限体积法在高超声速模拟中具有诸多优点。一方面，由于其基于守恒原理，能够准确地满足质量、动量和能量守恒定律，这对于高超声速流动中复杂的物理过程模拟至关重要。例如，在模拟高超声速飞行器再入大气层的过程中，准确的守恒性能够保证对激波与边界层相互作用、粘性耗散等现象的模拟精度。另一方面，有限体积法对不规则边界具有较好的适应性，可以方便地处理复杂的飞行器外形。通过生成贴体网格，能够更好地贴合飞行器表面的几何形状，提高计算精度。同时，有限体积法可以使用与流场相同的计算网格，方便地与流场求解的程序相结合，在高超声速气动热数值模拟中得到了广泛的应用。3.1.2有限差分法有限差分法（FiniteDifferenceMethod，FDM）是一种经典的数值方法，其基本原理是将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程中的导数用差商代替，推导出含有离散点上有限个未知数的差分方程组，通过求解差分方程的解来得到微分方程的近似数值解。有限差分法的核心在于通过泰勒级数展开式将变量的导数写成变量在不同时间或空间点值的差分形式。对于高超声速流动控制方程中的导数项，如\frac{\partialu}{\partialx}，可以采用不同的差分格式进行近似。常见的差分格式有前向差分、后向差分和中心差分。以前向差分为例，根据泰勒级数展开，u(x+\Deltax)可表示为u(x)+\frac{\partialu}{\partialx}\Deltax+\frac{\partial^{2}u}{\partialx^{2}}\frac{\Deltax^{2}}{2}+\cdots，则\frac{\partialu}{\partialx}的前向差分近似为\frac{u(x+\Deltax)-u(x)}{\Deltax}，其截断误差为O(\Deltax)。后向差分的近似为\frac{u(x)-u(x-\Deltax)}{\Deltax}，同样具有一阶截断误差。中心差分格式通过将u(x+\Deltax)和u(x-\Deltax)的泰勒展开式相减得到，\frac{\partialu}{\partialx}的中心差分近似为\frac{u(x+\Deltax)-u(x-\Deltax)}{2\Deltax}，其截断误差为O(\Deltax^{2})，精度比前向差分和后向差分更高。对于二阶导数\frac{\partial^{2}u}{\partialx^{2}}，常用的差分近似为\frac{u(x+\Deltax)-2u(x)+u(x-\Deltax)}{\Deltax^{2}}，截断误差为O(\Deltax^{2})。在高超声速气动热数值模拟中，利用有限差分法求解控制方程时，首先需要根据问题的特点和计算精度要求选择合适的差分格式，对控制方程进行离散。以二维非定常N-S方程为例，在笛卡尔坐标系下，将时间和空间进行离散，时间步长为\Deltat，空间步长在x和y方向分别为\Deltax和\Deltay。对控制方程中的各项导数采用相应的差分格式进行近似，将偏微分方程转化为关于网格节点上变量的差分方程组。然后，根据初始条件和边界条件，通过迭代求解差分方程组，得到网格节点上的数值解。有限差分法具有计算简单、程序易实现的优点，尤其适用于规则区域的数值计算。在高超声速流动模拟中，对于一些简单几何形状的飞行器，如平板、圆柱等，有限差分法能够快速有效地求解控制方程，得到流场参数的分布。然而，有限差分法也存在一些局限性。一方面，它对不规则边界的处理能力较弱，当计算区域的边界形状复杂时，生成合适的差分网格较为困难，可能会导致计算精度下降。另一方面，有限差分法的精度在很大程度上依赖于网格的精细程度，为了提高计算精度，需要加密网格，这会导致计算量急剧增加，对计算机的存储和计算能力提出较高要求。3.1.3谱方法谱方法（SpectralMethod）是一种基于函数逼近的数值方法，其基本原理是将求解函数展开成一系列已知的正交函数的线性组合，通过求解展开系数来得到原问题的近似解。常用的正交函数有傅里叶级数、切比雪夫多项式、勒让德多项式等。以傅里叶谱方法为例，假设待求解的函数u(x)在区间[-L,L]上，可将其展开为傅里叶级数：u(x)=\sum_{k=-\infty}^{\infty}a_{k}e^{i\frac{k\pi}{L}x}其中，a_{k}为傅里叶系数，可通过对u(x)与e^{-i\frac{k\pi}{L}x}在区间[-L,L]上进行积分得到：a_{k}=\frac{1}{2L}\int_{-L}^{L}u(x)e^{-i\frac{k\pi}{L}x}dx在实际计算中，通常只取有限项进行截断，即u(x)\approx\sum_{k=-N}^{N}a_{k}e^{i\frac{k\pi}{L}x}，其中N为截断阶数。将待求解的偏微分方程中的函数用正交函数展开后，代入方程中，利用正交函数的性质，将偏微分方程转化为关于展开系数的代数方程组。例如，对于一个含有\frac{\partialu}{\partialx}的偏微分方程，对u(x)的傅里叶展开式求导可得\frac{\partialu}{\partialx}=\sum_{k=-\infty}^{\infty}a_{k}i\frac{k\pi}{L}e^{i\frac{k\pi}{L}x}，将其代入原方程，通过对各项进行积分运算，得到关于a_{k}的代数方程。然后，求解这些代数方程，得到展开系数a_{k}，进而得到原函数u(x)的近似解。谱方法的优点是具有极高的精度，对于具有光滑解的问题，随着截断阶数的增加，谱方法的解能够以指数速度收敛到精确解。这使得谱方法在一些对精度要求极高的高超声速气动热问题中具有独特的优势，如研究高超声速边界层的稳定性和转捩问题时，谱方法能够准确捕捉边界层内的小尺度流动结构和物理量的细微变化。然而，谱方法的适用范围相对有限，主要适用于具有光滑解的问题。当问题中存在奇点、间断或复杂的边界条件时，谱方法的精度会受到严重影响。此外，谱方法的计算量通常较大，尤其是在处理高维问题时，计算成本会显著增加，这在一定程度上限制了其在实际工程中的广泛应用。3.2数值方法对比与选择有限体积法、有限差分法和谱方法在高超声速气动热数值模拟中各有优劣，其适用性与具体的高超声速流动场景密切相关。有限体积法由于其基于守恒原理，能够精确满足质量、动量和能量守恒定律，在模拟高超声速流动中复杂的物理过程时具有显著优势。在模拟高超声速飞行器再入大气层时，飞行器周围的流场存在激波与边界层的强烈相互作用、粘性耗散等复杂现象，有限体积法能够准确捕捉这些物理过程，保证模拟结果的准确性。该方法对不规则边界具有良好的适应性，能够方便地处理复杂的飞行器外形。通过生成贴体网格，有限体积法可以更好地贴合飞行器表面的几何形状，提高计算精度。在模拟具有复杂外形的高超声速飞行器时，如乘波体飞行器，有限体积法能够根据飞行器的外形特点生成合适的网格，准确计算流场参数和气动热分布。然而，有限体积法在处理一些高精度要求的问题时，可能需要更细的网格，从而导致计算量增加。有限差分法计算简单，程序易于实现，尤其适用于规则区域的数值计算。在高超声速流动模拟中，对于一些简单几何形状的飞行器，如平板、圆柱等，有限差分法能够快速有效地求解控制方程，得到流场参数的分布。在研究高超声速平板边界层流动时，有限差分法可以方便地对平板表面的流场进行离散和计算，快速得到边界层内的速度、温度等参数的分布。但是，有限差分法对不规则边界的处理能力较弱，当计算区域的边界形状复杂时，生成合适的差分网格较为困难，可能会导致计算精度下降。此外，有限差分法的精度在很大程度上依赖于网格的精细程度，为了提高计算精度，需要加密网格，这会显著增加计算量，对计算机的存储和计算能力提出较高要求。谱方法具有极高的精度，对于具有光滑解的问题，随着截断阶数的增加，谱方法的解能够以指数速度收敛到精确解。这使得谱方法在一些对精度要求极高的高超声速气动热问题中具有独特的优势，如研究高超声速边界层的稳定性和转捩问题时，谱方法能够准确捕捉边界层内的小尺度流动结构和物理量的细微变化。然而，谱方法的适用范围相对有限，主要适用于具有光滑解的问题。当问题中存在奇点、间断或复杂的边界条件时，谱方法的精度会受到严重影响。在模拟高超声速飞行器表面存在激波等间断现象时，谱方法可能无法准确捕捉激波的位置和强度，导致模拟结果偏差较大。此外，谱方法的计算量通常较大，尤其是在处理高维问题时，计算成本会显著增加，这在一定程度上限制了其在实际工程中的广泛应用。在选择数值方法时，需要综合考虑多种因素。对于高超声速飞行器的工程设计应用，通常需要在计算精度、计算效率和计算成本之间进行权衡。如果飞行器外形复杂，且对计算精度要求较高，有限体积法是较为合适的选择。它能够在保证计算精度的前提下，较好地处理复杂外形，满足工程设计的需求。若研究问题主要关注简单几何形状飞行器的基本流动特性，且对计算效率要求较高，有限差分法可以快速提供初步的计算结果，为后续研究提供参考。而对于一些对精度要求极高，且问题具有光滑解的基础研究，如高超声速边界层稳定性的理论研究，谱方法则能够发挥其高精度的优势，提供准确的模拟结果。在实际应用中，还可以根据具体情况对不同的数值方法进行组合使用，充分发挥各自的优点，以提高数值模拟的效果。四、数值模拟关键问题与处理策略4.1网格生成技术在高超声速气动热数值模拟中，网格生成是一项至关重要的基础工作，其质量和特性直接影响到数值模拟的精度、效率以及计算结果的可靠性。合适的网格能够准确地捕捉流场的细节特征，如激波的位置和强度、边界层的厚度和结构等，从而为气动热的精确计算提供保障。同时，高效的网格生成方法能够减少计算资源的消耗，提高模拟的效率，使得大规模的数值模拟成为可能。根据网格的拓扑结构和生成方式，可将其分为结构网格、非结构网格以及结合两者优点的混合网格，它们各自具有独特的特点和适用场景。4.1.1结构网格生成结构网格是一种具有规则拓扑结构的网格形式，其网格节点在空间上呈有序排列，如同矩形或六面体网格在规则区域中的分布。在二维情况下，结构网格可看作是由一系列平行的直线或曲线相互交织而成，每个网格单元都具有明确的邻接关系和几何形状；在三维情况下，结构网格则类似于由六面体单元紧密堆积而成，每个单元的面、边和顶点都与相邻单元有清晰的连接。这种规则的结构使得结构网格在数据存储和计算过程中具有较高的效率，能够方便地进行网格节点的编号和寻址，减少内存的占用。结构网格的生成方法主要包括代数方法和微分方程方法。代数方法是基于简单的数学变换和插值原理来生成网格。以二维区域为例，可通过对边界进行参数化处理，然后利用双线性插值或其他代数插值方法，在边界内部生成规则的网格。这种方法的优点是计算简单、速度快，能够快速生成初始网格。然而，代数方法对于复杂几何形状的适应性较差，当计算区域的边界存在弯曲、扭曲或不规则形状时，生成的网格容易出现变形、重叠或疏密分布不均匀的问题。微分方程方法则是通过求解椭圆型、双曲型或抛物型微分方程来生成网格。以椭圆型方程为例，将计算区域的边界条件作为方程的约束，通过迭代求解方程，使得网格节点在满足边界条件的同时，在区域内部达到某种优化的分布。这种方法生成的网格具有良好的正交性和光滑性，能够较好地适应复杂的几何形状，在边界层区域和复杂外形的模拟中表现出较高的精度。但是，微分方程方法的计算过程相对复杂，需要较大的计算量和较长的计算时间，对计算机的性能要求较高。在高超声速飞行器的数值模拟中，对于一些外形相对规则的部件，如平板机翼、圆柱机身等，结构网格具有显著的优势。以平板机翼的模拟为例，可采用结构网格对机翼表面和周围流场进行划分。在机翼表面，通过合理设置网格间距，能够精确捕捉边界层内的速度和温度梯度变化；在机翼周围的流场中，结构网格能够清晰地描述激波的传播和反射，准确计算激波与边界层的相互作用。在圆柱机身的模拟中，结构网格可以方便地沿着圆柱的轴线和圆周方向进行布置，有效地模拟圆柱表面的压力分布和气动热传递过程。由于结构网格的规则性，在进行数值计算时，可以采用高效的数值算法，如快速傅里叶变换（FFT）等，进一步提高计算效率。4.1.2非结构网格生成非结构网格是一种没有规则拓扑结构的网格形式，其网格节点和单元的分布较为灵活，能够适应各种复杂的几何形状。非结构网格的单元形状多样，包括三角形、四面体、多边形等，这些单元可以根据几何模型的特点进行自由组合和排列。在模拟复杂外形的高超声速飞行器时，非结构网格能够精确地贴合飞行器表面的曲线和曲面，在机翼的前缘、后缘、机身的拐角等复杂部位，非结构网格可以通过局部加密的方式，准确捕捉流场的细节特征。非结构网格的生成算法主要有Delaunay三角剖分算法、推进阵面法等。Delaunay三角剖分算法是基于点集的一种网格生成方法，其基本思想是将给定的离散点集进行三角剖分，使得生成的三角形网格满足Delaunay条件，即任意一个三角形的外接圆内不包含其他的离散点。通过这种方式生成的三角形网格具有较好的质量，三角形的形状较为规则，能够有效地减少数值计算中的误差。推进阵面法是从物体的边界开始，逐步向内部推进生成网格。首先在物体表面生成一层初始的网格，然后将这层网格作为推进阵面，沿着与边界垂直的方向向内部推进，在推进过程中不断生成新的网格单元。这种方法能够较好地控制网格的生长方向和密度分布，在边界层区域可以生成高质量的网格。在高超声速飞行器的数值模拟中，非结构网格在处理复杂外形时具有明显的优势。在模拟具有复杂外形的高超声速飞行器时，如带有复杂进气道、尾喷管和操纵面的飞行器，非结构网格能够根据飞行器的外形特点，灵活地生成网格。在进气道的模拟中，非结构网格可以准确地描述进气道内部的复杂流道形状，捕捉气流在进气道内的压缩、膨胀和分离等现象；在尾喷管的模拟中，非结构网格能够精细地刻画尾喷管出口的流场结构，计算喷流与周围空气的相互作用。此外，非结构网格还可以方便地进行局部网格加密和自适应调整，当流场中出现激波、边界层分离等复杂现象时，能够根据流场的变化自动调整网格的疏密程度，提高模拟的精度。4.1.3混合网格技术混合网格技术是结合了结构网格和非结构网格的优点而发展起来的一种网格生成方法。在高超声速气动热数值模拟中，不同区域的流场特性和对网格的要求存在差异，混合网格技术能够根据这些差异，在不同区域采用不同类型的网格，从而实现计算精度和效率的优化。在边界层区域，由于边界层内的速度和温度梯度变化剧烈，需要采用高质量的结构化网格来准确捕捉这些变化。结构化网格在边界层内可以沿着壁面法线方向进行加密，提高对边界层内物理量变化的分辨率。在远离边界层的区域，流场相对较为均匀，对网格的分辨率要求相对较低，此时可以采用非结构网格来填充，非结构网格的灵活性能够适应复杂的流场拓扑结构，同时减少网格数量，降低计算成本。在模拟高超声速飞行器时，可在飞行器表面的边界层区域生成结构化的棱柱网格，这种网格能够很好地贴合飞行器表面，并且在壁面法线方向具有较高的分辨率，能够准确计算边界层内的粘性耗散和热传导。在边界层外的区域，采用非结构的四面体网格进行填充，四面体网格可以根据流场的复杂程度进行灵活布置，在激波等关键区域进行局部加密，准确捕捉激波的位置和强度。通过这种混合网格的方式，既保证了边界层区域的计算精度，又提高了整体的计算效率。在实际应用中，混合网格技术还需要解决不同类型网格之间的连接和过渡问题，确保网格的连续性和光滑性，以避免在数值计算中产生误差。4.2边界条件处理在高超声速气动热数值模拟中，边界条件的准确设定对于模拟结果的准确性和可靠性起着至关重要的作用。边界条件的设置不仅影响流场的计算结果，还直接关系到气动热的预测精度。合理的边界条件能够准确地反映高超声速飞行器在实际飞行过程中的物理环境，为数值模拟提供真实可靠的边界信息。本节将详细介绍壁面边界条件和远场边界条件的设置与应用。4.2.1壁面边界条件壁面边界条件主要包括无滑移条件、等温条件和绝热条件，这些条件的选择取决于具体的物理问题和模拟需求。无滑移条件是壁面边界条件中最常用的一种，它基于流体的粘性特性，假设流体与壁面之间没有相对滑移，即壁面处流体的速度等于壁面的速度。在数学表达式上，对于静止的壁面，有u=v=w=0，其中u、v、w分别为流体在x、y、z方向的速度分量。无滑移条件在模拟高超声速飞行器的气动热问题中具有重要意义，它能够准确描述边界层内流体与壁面之间的相互作用。在边界层内，由于粘性的作用，流体速度从壁面处的零逐渐增加到主流速度，无滑移条件确保了这种速度变化的连续性。在模拟高超声速飞行器机翼表面的流场时，无滑移条件能够准确捕捉边界层内的速度梯度和粘性耗散，从而为计算机翼表面的热流密度提供准确的流场信息。等温条件是指在壁面处设定一个固定的温度值，即T=T_{w}，其中T为壁面处的温度，T_{w}为给定的壁面温度。这种边界条件适用于一些对壁面温度有明确要求的情况，在研究高超声速飞行器热防护系统的性能时，常常需要给定热防护层表面的温度，以评估热防护系统对飞行器结构的保护效果。通过设置等温条件，可以准确计算热量在热防护层内部的传递过程，为热防护系统的设计和优化提供依据。在模拟高超声速飞行器再入大气层的过程中，如果已知热防护层表面的最高允许温度，就可以将该温度作为等温条件输入到数值模拟中，分析热防护层内部的温度分布和热应力情况，从而判断热防护层是否能够满足飞行器的热防护需求。绝热条件则假设壁面与流体之间没有热量交换，即壁面处的热流密度为零，数学表达式为q_{w}=0，其中q_{w}为壁面处的热流密度。绝热条件在一些情况下能够简化计算，并且在某些实际问题中具有一定的合理性。在研究高超声速飞行器表面的局部热现象时，如果关注的是边界层内的流动特性和气动热的产生机制，而对壁面与外界的热量交换不太关心，就可以采用绝热条件。在模拟高超声速飞行器头部激波与边界层相互作用时，绝热条件可以帮助我们更专注于激波和边界层内的物理过程，而不必考虑壁面热交换对这些过程的影响。在某些热防护材料的初步研究中，也可以先采用绝热条件来分析材料在气动热环境下的热响应，为后续更精确的研究提供基础。4.2.2远场边界条件远场边界条件用于描述计算域边界远离飞行器的区域的流动状态，它对模拟精度有着重要的影响。在高超声速气动热数值模拟中，常用的远场边界条件包括压力远场条件和特征边界条件。压力远场条件是指在远场边界上给定压力、密度和速度等参数，以确定远场的流动状态。在高超声速流动中，通常假设远场为均匀来流，即给定远场的马赫数M_{\infty}、压力p_{\infty}和温度T_{\infty}等参数。根据理想气体状态方程和等熵流动关系，可以由这些给定参数计算出远场的密度\rho_{\infty}和速度u_{\infty}等其他参数。压力远场条件的设置相对简单，计算效率较高，在一些对计算精度要求不是特别高的工程应用中得到了广泛的应用。在初步设计阶段对高超声速飞行器的气动热性能进行估算时，采用压力远场条件可以快速得到流场和气动热的大致分布，为后续的详细设计提供参考。然而，压力远场条件在处理一些复杂的流动现象时存在一定的局限性，它无法准确捕捉激波在远场的传播和反射等现象，可能会导致模拟结果的误差。特征边界条件是基于特征线理论来设置远场边界条件的。在高超声速流动中，通过分析控制方程的特征线，可以确定在远场边界上哪些物理量是已知的，哪些是需要通过内场计算结果外推得到的。特征边界条件能够更准确地考虑波的传播特性，对于捕捉激波在远场的行为具有更好的效果。在模拟高超声速飞行器的复杂流场时，激波在远场的传播和反射会对飞行器周围的流场产生重要影响，采用特征边界条件可以更精确地模拟这些现象，提高模拟结果的准确性。在研究高超声速飞行器的进气道流场时，激波与进气道壁面的相互作用以及激波在远场的传播会影响进气道的性能，特征边界条件能够准确捕捉这些现象，为进气道的设计和优化提供更可靠的依据。然而，特征边界条件的计算相对复杂，需要对控制方程进行深入的分析和处理，计算成本较高。4.3湍流模型选择在高超声速气动热数值模拟中，湍流模型的选择对模拟结果的准确性和可靠性起着关键作用。高超声速流动中的湍流现象极为复杂，涉及到多尺度的涡旋运动、强烈的粘性耗散以及激波与湍流的相互作用等。不同的湍流模型基于不同的假设和理论，对这些复杂物理现象的描述能力各不相同。因此，深入了解常见湍流模型的原理、特点以及适用范围，对于准确模拟高超声速气动热具有重要意义。4.3.1常见湍流模型Spalart-Allmaras模型是一种一方程湍流模型，由Spalart和Allmaras于1992年提出。该模型通过求解一个关于涡粘性的输运方程来模拟湍流。其基本方程为：\frac{\partial\tilde{\nu}}{\partialt}+u_{j}\frac{\partial\tilde{\nu}}{\partialx_{j}}=c_{b1}\left(1-f_{t2}\right)\tilde{S}\tilde{\nu}+\frac{1}{\sigma}\left[\frac{\partial}{\partialx_{j}}\left((\nu+\tilde{\nu})\frac{\partial\tilde{\nu}}{\partialx_{j}}\right)+c_{b2}\left(\frac{\partial\tilde{\nu}}{\partialx_{j}}\right)^{2}\right]-c_{w1}f_{w}\left(\frac{\tilde{\nu}}{d}\right)^{2}其中，\tilde{\nu}为修正的涡粘性，u_{j}为速度分量，\tilde{S}为与应变率相关的参数，\sigma、c_{b1}、c_{b2}、c_{w1}等为模型常数，f_{t2}、f_{w}为经验函数，d为到壁面的距离。该模型的特点是计算量相对较小，收敛性好。由于只需求解一个方程，在计算资源有限的情况下，能够快速得到计算结果。它适用于航空领域的外流场模拟，特别是绕流过程，对于粗网格也能给出较为合理的结果。在模拟高超声速飞行器的外部绕流时，Spalart-Allmaras模型能够较好地捕捉边界层的发展和分离，计算出飞行器表面的压力和摩擦力分布。然而，该模型对复杂流动的预测能力有限，在处理涉及强分离、复杂激波与湍流相互作用的流动时，模拟结果可能存在较大误差。Standardk-ε模型是一种经典的二方程湍流模型，通过求解湍流动能k和湍流耗散率\varepsilon的输运方程来描述湍流。其湍流动能k的方程为：\frac{\partial(\rhok)}{\partialt}+\frac{\partial(\rhoku_{i})}{\partialx_{i}}=\frac{\partial}{\partialx_{j}}\left[\left(\mu+\frac{\mu_{t}}{\sigma_{k}}\right)\frac{\partialk}{\partialx_{j}}\right]+G_{k}-\rho\varepsilon湍流耗散率\varepsilon的方程为：\frac{\partial(\rho\varepsilon)}{\partialt}+\frac{\partial(\rho\varepsilonu_{i})}{\partialx_{i}}=\frac{\partial}{\partialx_{j}}\left[\left(\mu+\frac{\mu_{t}}{\sigma_{\varepsilon}}\right)\frac{\partial\varepsilon}{\partialx_{j}}\right]+C_{1\varepsilon}\frac{\varepsilon}{k}G_{k}-C_{2\varepsilon}\rho\frac{\varepsilon^{2}}{k}其中，\rho为密度，\mu为分子粘性系数，\mu_{t}为湍流粘性系数，\sigma_{k}、\sigma_{\varepsilon}、C_{1\varepsilon}、C_{2\varepsilon}为模型常数，G_{k}为湍动能的产生项。该模型适用于高雷诺数的完全湍流，包含粘性热、浮力、压缩性等选项，适合完全发展的湍流流动，如管道流动和外部绕流。在高超声速管道流动模拟中，Standardk-ε模型能够准确计算流场中的湍动能和耗散率分布，进而得到管道壁面的热流密度。它基于Boussinesq假设，将雷诺应力与平均速度梯度联系起来，使得计算相对简单。但该模型在处理低雷诺数流动、强旋流和各向异性湍流时，精度会有所下降。RNGk-ε模型在Standardk-ε模型的基础上，对\varepsilon方程进行了改进。通过重整化群理论，考虑了湍流的小尺度效应，在\varepsilon方程中引入了一个附加项，以改善对高应变流动的模拟能力。其\varepsilon方程为：\frac{\partial(\rho\varepsilon)}{\partialt}+\frac{\partial(\rho\varepsilonu_{i})}{\partialx_{i}}=\frac{\partial}{\partialx_{j}}\left[\alpha_{\varepsilon}\mu_{eff}\frac{\partial\varepsilon}{\partialx_{j}}\right]+C_{1\varepsilon}\frac{\varepsilon}{k}\left(1-\frac{\eta}{\eta_{0}}\right)G_{k}-C_{2\varepsilon}\rho\frac{\varepsilon^{2}}{k}其中，\alpha_{\varepsilon}为与湍流普朗特数相关的参数，\eta为与应变率和旋转率相关的参数，\eta_{0}为经验常数。RNGk-ε模型适用于预测中等强度的旋流和低雷诺数流动。在模拟高超声速飞行器进气道内的流动时，由于进气道内存在复杂的旋流和边界层分离现象，RNGk-ε模型能够更好地捕捉这些流动特征，比Standardk-ε模型具有更高的精度。该模型考虑了湍流的各向异性和小尺度效应，对复杂流动的适应性更强。然而，由于控制方程中额外的功能和非线性，RNGk-ε模型比Standardk-ε模型多消耗10-15%的CPU时间。4.3.2模型评估与选择在高超声速气动热数值模拟中，不同的流动情况对湍流模型的要求各异。对于简单的高超声速外部绕流，如平板绕流，流动相对规则，没有明显的分离和复杂的激波与湍流相互作用。此时，Spalart-Allmaras模型因其计算量小、收敛性好的特点，能够快速准确地模拟边界层的发展，计算出平板表面的热流密度，是较为合适的选择。而在模拟高超声速飞行器的进气道流动时，进气道内存在强烈的压缩、旋流和边界层分离现象，流动呈现出高度的复杂性和各向异性。这种情况下，RNGk-ε模型由于考虑了小尺度效应和各向异性，能够更好地捕捉进气道内的复杂流动特征，提供更准确的模拟结果。在选择湍流模型时，还需要考虑计算资源的限制。如果计算资源有限，无法进行大规模的计算，那么计算量较小的Spalart-Allmaras模型或Standardk-ε模型可能更为合适。它们可以在有限的资源下快速得到计算结果，虽然精度可能相对较低，但能够为工程设计提供初步的参考。相反，如果计算资源充足，且对模拟精度要求较高，对于复杂的高超声速流动问题，可以选择计算量较大但精度更高的模型，如RNGk-ε模型或雷诺应力模型（RSM）。RSM模型通过求解雷诺应力的输运方程，能够更准确地描述湍流的各向异性，在处理复杂流动时具有明显的优势。但RSM模型的计算成本较高，需要消耗大量的计算资源和时间。为了选择合适的湍流模型，还可以通过与实验数据对比来评估模型的准确性。将不同湍流模型的模拟结果与实验测量数据进行比较，分析模型对流动参数和气动热分布的预测能力。如果模拟结果与实验数据吻合较好，则说明该模型能够较好地描述所研究的流动现象，具有较高的可靠性。在模拟高超声速飞行器机翼的气动热时，可以将不同湍流模型计算得到的机翼表面热流密度与风洞实验测量结果进行对比。如果某个模型的计算结果与实验数据在热流密度的大小和分布趋势上都较为一致，那么该模型在模拟该机翼的气动热问题时就是一个较好的选择。通过不断地对比和评估，可以逐渐积累经验，提高对湍流模型选择的准确性和合理性。五、数值模拟案例分析5.1高超声速飞行器外形模拟5.1.1钝头体气动热模拟为了验证数值模拟模型与方法的有效性，以钝头体为研究对象进行气动热模拟。钝头体是高超声速飞行器中常见的外形部件，其头部的气动热问题尤为突出，对其进行模拟分析具有重要的代表性和参考价值。在本次模拟中，设定钝头体的头部半径为0.1m，长度为1m。来流条件为马赫数M=8，来流温度T_{\infty}=270K，来流压力p_{\infty}=5000Pa。采用有限体积法对控制方程进行离散求解，使用结构网格对钝头体周围的流场进行划分，在钝头体表面和边界层区域进行网格加密，以提高计算精度。通过数值模拟，得到了钝头体表面的热流密度分布和温度分布。从热流密度分布结果来看，钝头体头部驻点处的热流密度最高，这是由于在驻点处，气流速度急剧降为零，动能全部转化为内能，导致温度急剧升高，从而产生极高的热流密度。随着远离驻点，热流密度逐渐降低。在钝头体的肩部，由于激波与边界层的相互作用，热流密度出现了一个局部峰值。将模拟得到的热流密度分布与经典的费-里德尔（Fay-Riddell）平衡边界层驻点热流密度公式计算结果进行对比。费-里德尔公式是计算钝头体驻点热流密度的经典公式，其表达式为：q_{s}=0.76\rho_{\infty}^{0.5}u_{\infty}^{3}r_{n}^{0.5}\frac{h_{s}-h_{w}}{h_{s}-h_{r}}其中，q_{s}为驻点热流密度，\rho_{\infty}为来流密度，u_{\infty}为来流速度，r_{n}为头部半径，h_{s}为驻点焓，h_{w}为壁面焓，h_{r}为恢复焓。经过计算对比，数值模拟得到的驻点热流密度与费-里德尔公式计算结果在误差允许范围内基本一致，验证了数值模拟方法在计算驻点热流密度方面的准确性。从温度分布结果来看，钝头体头部温度最高，向尾部逐渐降低。在头部区域，由于热流密度大，热量传递迅速，温度梯度较大；在尾部区域，热流密度较小，温度变化相对平缓。将模拟得到的温度分布与相关实验数据进行对比。在相同的来流条件和钝头体几何参数下，实验测量得到了钝头体表面的温度分布。通过对比发现，数值模拟得到的温度分布与实验数据在整体趋势上吻合良好，在头部和肩部等关键部位的温度值也较为接近。虽然存在一定的误差，但误差范围在合理区间内，这可能是由于实验测量误差、数值模拟中的模型简化以及计算精度等因素导致的。总体而言，通过与经典理论公式和实验数据的对比，验证了所采用的数值模拟模型与方法在钝头体气动热模拟中的有效性和可靠性，能够准确地预测钝头体表面的热流密度和温度分布，为高超声速飞行器钝头体部件的热防护设计提供了有力的支持。5.1.2复杂外形飞行器模拟为了进一步研究高超声速飞行器在实际飞行中的气动热特性，对具有复杂外形的飞行器进行数值模拟。该复杂外形飞行器包含机身、机翼、尾翼以及进气道等多个部件，其外形设计旨在满足高超声速飞行的气动性能要求。在模拟过程中，来流马赫数设定为M=6，来流温度T_{\infty}=250K，来流压力p_{\infty}=8000Pa。考虑到飞行器外形的复杂性，采用混合网格技术对计算区域进行网格划分，在飞行器表面和边界层区域使用结构化的棱柱网格，以精确捕捉边界层内的流动和传热特性；在远离飞行器表面的区域使用非结构化的四面体网格，以提高网格生成的效率和对复杂流场的适应性。通过数值模拟，得到了飞行器表面不同部位的气动热分布特征。在飞行器头部，由于气流的强烈压缩和驻点效应，热流密度和温度均达到很高的值。头部的热流密度峰值远高于其他部位，这是因为在头部，气流速度急剧降为零，动能大量转化为热能，导致温度急剧升高，进而产生极高的热流密度。在机翼前缘，由于其处于气流的迎风面，且曲率变化较大，激波与边界层的相互作用强烈，使得该部位的热流密度也较高。机翼前缘的热流密度分布呈现出明显的不均匀性，在靠近头部的区域热流密度较大，随着向后缘延伸逐渐减小。这是由于靠近头部的区域，激波强度较大，与边界层的相互作用更剧烈，产生的热量更多；而向后缘延伸时，激波强度逐渐减弱，热流密度也随之降低。机身表面的热流密度和温度分布相对较为均匀，但在一些局部区域，如机身与机翼的连接处、进气道附近等，由于气流的干扰和分离，热流密度会出现局部升高的现象。在机身与机翼的连接处，气流的流动受到机翼的影响，形成了复杂的流场结构，导致局部热流密度增加。进气道附近的热流密度升高则是因为进气道内的气流压缩和膨胀过程会对周围流场产生干扰，使得该区域的气动热现象加剧。尾翼部分的热流密度和温度相对较低，这是由于尾翼处于飞行器的后方，受到的气流冲击相对较弱，且尾翼表面的气流速度相对较低，粘性耗散和激波作用产生的热量较少。然而，在尾翼的前缘和后缘，由于气流的绕流和分离，热流密度仍会出现一定程度的升高。通过对复杂外形飞行器不同部位气动热分布特征的分析，可以为飞行器的热防护系统设计提供有针对性的指导。对于热流密度和温度较高的部位，如头部、机翼前缘等，应采用耐高温性能更好的热防护材料，并优化热防护结构设计，以提高其热防护能力。在机身与机翼连接处、进气道附近等局部热流密度升高的区域，可以通过改进结构设计，减少气流的干扰和分离，降低热流密度；或者采用局部强化的热防护措施，如增加隔热层厚度、使用热导率更低的材料等。对于尾翼部分，虽然整体热环境相对较好，但仍需关注前缘和后缘等关键部位的热防护，确保尾翼在高超声速飞行条件下的结构完整性和可靠性。5.2不同飞行条件下的模拟5.2.1不同马赫数为了深入研究不同马赫数下高超声速飞行器的气动热特性，对马赫数分别为5、7、9的飞行条件进行数值模拟。保持来流温度T_{\infty}=260K，来流压力p_{\infty}=6000Pa不变，采用与前文相同的数值模拟方法和网格划分策略。模拟结果表明，随着马赫数的增加，飞行器表面的热流密度和温度显著升高。当马赫数从5增加到7时，飞行器头部驻点处的热流密度增加了约1.5倍，温度升高了约300K；当马赫数进一步增加到9时，驻点热流密度又增加了约1.8倍，温度升高了约500K。这是因为马赫数的增加意味着飞行器飞行速度的提高，气流与飞行器表面的相互作用更加剧烈。在高马赫数下，激波强度增强，气体被压缩得更厉害，导致激波后的气体温度和压力大幅升高。同时，边界层内的粘性耗散也随着马赫数的增加而加剧，产生更多的热量。在高马赫数下，边界层内的速度梯度更大，粘性力做功更多，将更多的机械能转化为热能，从而使热流密度和温度升高。不同马赫数下，热流密度和温度的分布规律也有所变化。在低马赫数时，热流密度和温度分布相对较为均匀，从头部到尾部的变化较为平缓。随着马赫数的增加，头部和机翼前缘等部位的热流密度和温度峰值更加突出，与其他部位的差异增大。在马赫数为9时，头部驻点处的热流密度是机身中部热流密度的5倍以上，温度也高出约800K。这是由于在高马赫数下，激波与边界层的相互作用在头部和机翼前缘等部位更加剧烈，导致这些部位的热流密度和温度急剧升高。同时，高马赫数下边界层的厚度相对变薄，使得热量更集中在表面附近，加剧了热流密度和温度分布的不均匀性。马赫数的变化还会影响飞行器表面的压力分布。随着马赫数的增加，飞行器表面的压力系数增大，压力分布更加不均匀。在头部和机翼前缘等部位，压力系数显著增大，这是由于激波的压缩作用导致这些部位的压力升高。压力分布的变化又会反过来影响热流密度和温度分布，因为压力梯度会影响边界层内的流动和传热特性。在压力梯度较大的区域，边界层内的流动更加复杂，粘性耗散增加，从而导致热流密度和温度升高。5.2.2不同攻角为了探究不同攻角对高超声速飞行器气动热的影响，对攻角分别为0°、5°、10°的飞行条件进行数值模拟。保持来流马赫数M=6，来流温度T_{\infty}=250K，来流压力p_{\infty}=8000Pa不变，采用相同的数值模拟方法和网格划分策略。模拟结果显示，攻角对飞行器表面的气动热分布和热流密度有显著影响。当攻角为0°时，飞行器表面的气动热分布相对较为对称，热流密度在头部和机翼前缘较高，向尾部逐渐降低。随着攻角的增加，飞行器迎风面和背风面的气动热差异逐渐增大。在攻角为5°时，迎风面的热流密度明