基于格子波兹曼方法与NND格式的高速流动模拟研究

基于格子波兹曼方法与NND格式的高速流动模拟研究一、引言1.1研究背景与意义在现代科学与工程领域，高速流动现象广泛存在，其模拟研究具有极其重要的价值。在航空航天领域，飞行器在大气层内高速飞行时，周围空气的高速流动特性直接影响飞行器的气动性能、飞行稳定性与安全性。例如，飞机在超音速飞行时，机翼和机身表面会产生复杂的激波结构，这些激波与边界层相互作用，可能导致飞行器的阻力急剧增加、升力系数下降，甚至引发结构的振动和疲劳破坏。准确模拟高速流动，能够为飞行器的外形设计、气动布局优化提供关键依据，从而降低飞行阻力、提高燃油效率，提升飞行器的整体性能和可靠性。在高超声速飞行器的研发中，了解高速气流在飞行器表面的流动特性，对于热防护系统的设计至关重要，可有效避免飞行器在高速飞行时因气动加热而导致的结构损坏。在能源领域，燃气轮机、航空发动机等动力设备内部的气体流动也处于高速状态。精确模拟这些高速流动过程，有助于优化发动机的内部结构，提高燃烧效率，降低污染物排放，增强能源利用效率。在火箭发动机的设计中，对推进剂燃烧室内高速气流的模拟，能够为发动机的性能提升和可靠性保障提供有力支持。在工业领域，高速风洞实验是研究高速流动特性的重要手段之一，但实验成本高昂、周期长，且受到实验条件的限制。相比之下，数值模拟方法具有成本低、灵活性高、可重复性强等优势，能够对复杂的高速流动现象进行深入研究，为实验提供理论指导和数据补充。格子波兹曼方法（LatticeBoltzmannMethod，LBM）作为一种新兴的计算流体力学方法，近年来受到了广泛关注。它基于介观尺度的粒子分布函数，通过模拟粒子在格子上的碰撞和迁移过程来求解流体动力学方程。LBM具有天然的并行性，能够充分利用现代计算机的多核计算资源，大大提高计算效率。其边界处理简单，易于处理复杂的几何边界条件，在模拟复杂流场时具有独特的优势。在模拟多孔介质中的流体流动时，LBM能够轻松处理多孔介质复杂的内部结构边界。然而，经典的LBM在模拟高速流动时存在一定的局限性，如在高马赫数下数值稳定性较差、无法准确捕捉激波等。NND（无波动、无自由参数耗散）格式是一种高精度的差分格式，具有良好的激波捕捉能力和数值稳定性。它通过对通量的分裂和重构，有效地抑制了数值解在激波附近的振荡，能够准确地模拟包含激波的复杂流场。将NND格式与格子波兹曼方法相结合，有望充分发挥两者的优势，为高速流动的模拟提供一种新的有效途径。这种结合可以弥补经典LBM在高速流动模拟中的不足，提高模拟的精度和可靠性，为航空航天、能源等领域的工程设计和科学研究提供更准确的数值模拟工具。对于研究高超声速飞行器在再入大气层时的复杂流场，结合后的方法能够更精确地模拟激波与边界层的相互作用，为飞行器的热防护设计和飞行控制提供更可靠的依据。1.2国内外研究现状格子波兹曼方法自20世纪80年代末被提出以来，在国内外得到了广泛的研究和应用。早期，研究主要集中在LBM的基本理论和模型的建立上，如离散速度模型的选择、碰撞算子的设计等。随着研究的深入，LBM在复杂流场模拟中的优势逐渐显现，其应用领域不断拓展。在国内，清华大学、北京大学、上海大学等高校的研究团队在LBM的理论和应用方面取得了一系列重要成果。清华大学的学者通过对LBM的理论分析，揭示了其与传统计算流体力学方法的内在联系，为LBM的进一步发展奠定了理论基础。上海大学的研究团队将LBM应用于多孔介质中的多相流模拟，成功地捕捉到了多相流的复杂界面现象，为相关领域的工程设计提供了有力支持。在国际上，美国、德国、法国等国家的科研机构也在积极开展LBM的研究工作。美国斯坦福大学的研究人员利用LBM模拟了微尺度下的流体流动，发现了一些传统理论无法解释的新现象，推动了微流体力学的发展。德国哥廷根大学的学者则将LBM与分子动力学方法相结合，研究了流体在纳米尺度下的行为，为纳米技术的发展提供了理论依据。然而，经典的LBM在模拟高速流动时存在一定的局限性，为了克服这些局限性，国内外学者提出了多种改进方法。一些学者通过改进碰撞算子，提高了LBM在高马赫数下的数值稳定性。文献提出了一种基于多松弛时间（MRT）的碰撞算子，该算子通过引入多个松弛时间，对不同的物理量进行独立的松弛，从而提高了LBM的数值稳定性和精度。还有一些学者通过引入高阶矩模型，增强了LBM对激波等复杂流动结构的捕捉能力。在国内，北京航空航天大学的研究团队提出了一种基于高阶矩的LBM模型，该模型通过考虑高阶矩的影响，能够更准确地捕捉激波和边界层等复杂流动结构。在国际上，法国巴黎高等师范学院的学者则通过引入非平衡态分布函数的高阶矩，改进了LBM的激波捕捉能力。NND格式由中国学者张涵信于20世纪80年代提出，该格式一经提出便受到了国内外计算流体力学领域的高度关注。NND格式基于无波动、无自由参数耗散的思想，通过对通量的分裂和重构，有效地抑制了数值解在激波附近的振荡，能够准确地模拟包含激波的复杂流场。在国内，中国空气动力研究与发展中心等科研机构在NND格式的理论研究和应用方面开展了大量工作。中国空气动力研究与发展中心的研究人员通过对NND格式的理论分析，提出了一系列改进措施，提高了格式的精度和计算效率。他们将NND格式应用于高超声速飞行器的流场模拟，准确地捕捉到了激波与边界层的相互作用，为飞行器的设计提供了重要参考。在国际上，美国、俄罗斯等国家的科研人员也对NND格式进行了深入研究。美国国家航空航天局（NASA）的研究团队将NND格式应用于航空发动机内部流场的模拟，取得了良好的效果。俄罗斯科学院的学者则将NND格式与其他数值方法相结合，研究了复杂外形飞行器的气动特性。近年来，将格子波兹曼方法与NND格式相结合的研究逐渐成为热点。这种结合旨在充分发挥LBM的并行性和边界处理优势以及NND格式的激波捕捉能力，为高速流动的模拟提供更有效的方法。上海大学的学者提出了一种将NND格式应用于格子玻尔兹曼方程的模型，通过使用Wang的NND模式，在二阶时间和空间精度上对经典算法进行了改进。该模型通过典型的一维和二维数值试验验证了其正确性，并对激波相互作用进行了数值模拟，结果与理论分析相符。国际上，一些研究团队也在探索LBM与NND格式的结合方式，尝试应用于不同的高速流动问题，如高超声速飞行器的绕流、激波管中的激波传播等。虽然目前该领域的研究取得了一定进展，但在模型的精度、稳定性以及适用范围等方面仍存在一些问题有待进一步解决。1.3研究目标与内容本研究旨在通过将格子波兹曼方法与NND格式相结合，改进高速流动的模拟效果，提高模拟的精度和稳定性，为相关工程领域提供更可靠的数值模拟工具。具体研究内容如下：格子波兹曼方法与NND格式的结合与改进：深入研究格子波兹曼方法和NND格式的基本原理，分析经典格子波兹曼方法在模拟高速流动时存在的局限性，以及NND格式在激波捕捉方面的优势。在此基础上，探索将NND格式引入格子波兹曼方法的有效途径，对结合后的算法进行改进和优化。通过对碰撞算子和通量计算方式的改进，提高算法在高马赫数下的数值稳定性和激波捕捉能力。研究不同参数对算法性能的影响，确定最优的参数设置，以提高模拟的精度和效率。数值试验与验证：利用改进后的算法，开展一系列典型的一维和二维高速流动数值试验，如激波管问题、超声速绕流问题等。将数值试验结果与理论解、实验数据或其他成熟数值方法的结果进行对比分析，验证改进算法的正确性和有效性。通过对比不同马赫数下的模拟结果，评估算法在不同高速流动条件下的性能表现。分析数值试验中出现的问题，进一步优化算法，提高其可靠性和通用性。复杂高速流动现象的模拟与分析：运用改进后的算法，对包含激波与边界层相互作用、激波干扰等复杂高速流动现象进行数值模拟。深入分析模拟结果，研究激波的形成、传播和相互作用规律，以及边界层对激波的影响。探讨复杂高速流动现象中的物理机制，为相关工程问题的解决提供理论支持。结合实际工程应用，如航空发动机内部流场、高超声速飞行器绕流等，进行数值模拟研究，为工程设计和优化提供有价值的参考。1.4研究方法与技术路线本研究将综合运用理论研究、数值试验和对比分析等多种方法，深入探究格子波兹曼方法与NND格式相结合在高速流动模拟中的应用，具体研究方法如下：理论研究：深入剖析格子波兹曼方法和NND格式的基本原理、数学模型以及适用条件。研究格子波兹曼方法在高速流动模拟中的局限性根源，如在高马赫数下数值稳定性问题的产生机制，以及NND格式在激波捕捉方面的优势原理，包括其通量分裂和重构的具体方式对激波捕捉的影响。探索将NND格式引入格子波兹曼方法的理论依据和可行途径，从数学层面分析结合后算法的稳定性、收敛性和精度等特性。数值试验：利用改进后的算法开展一系列典型的一维和二维高速流动数值试验。针对激波管问题，设置不同的初始条件和边界条件，模拟激波在管内的传播、反射和相互作用过程。对于超声速绕流问题，选择不同形状的物体，如圆柱、翼型等，模拟超声速气流绕过物体时的流场特性，包括激波的形成位置、强度以及边界层的发展情况。通过改变马赫数、雷诺数等参数，系统地研究不同高速流动条件下改进算法的性能表现。对比分析：将数值试验结果与理论解、实验数据或其他成熟数值方法的结果进行详细的对比分析。对于有理论解的问题，如简单激波管问题的解析解，直接比较数值解与理论解的差异，评估改进算法的准确性。对于有实验数据的情况，如超声速风洞实验中测量的流场参数，将数值模拟结果与实验数据进行对比，验证算法的可靠性。与其他成熟的数值方法，如有限体积法、有限差分法等在相同条件下的模拟结果进行比较，分析改进算法在精度、计算效率和稳定性等方面的优势与不足。本研究的技术路线如下：首先，深入研究格子波兹曼方法与NND格式的基本原理，分析经典格子波兹曼方法在模拟高速流动时的局限性以及NND格式在激波捕捉方面的优势。基于此，通过改进碰撞算子和通量计算方式，将NND格式引入格子波兹曼方法，对结合后的算法进行优化，确定最优的参数设置。接着，利用改进后的算法开展典型的一维和二维高速流动数值试验，将数值试验结果与理论解、实验数据或其他成熟数值方法的结果进行对比分析，验证改进算法的正确性和有效性。根据对比分析结果，进一步优化算法，提高其可靠性和通用性。最后，运用优化后的算法对包含激波与边界层相互作用、激波干扰等复杂高速流动现象进行数值模拟，深入分析模拟结果，研究复杂高速流动现象中的物理机制，为相关工程问题的解决提供理论支持和参考。二、格子波兹曼方法与NND格式基础2.1格子波兹曼方法（LBM）2.1.1LBM基本原理格子波兹曼方法（LBM）是一种基于介观尺度的计算流体力学方法，其基本原理源于对流体微观粒子动力学的模拟。与传统基于宏观连续介质假设的计算流体力学方法不同，LBM将流体视为由大量在离散格子节点上运动和相互作用的粒子组成。在LBM中，时间和空间都被离散化，流体的宏观行为通过粒子在格子上的碰撞和迁移过程来体现。具体而言，LBM定义了离散的速度集合，每个格子节点上的粒子具有特定的离散速度。在每个时间步，粒子首先按照各自的速度方向进行迁移，从一个格子节点移动到相邻的节点，这一过程模拟了流体的对流。粒子在节点处发生碰撞，碰撞过程改变了粒子的速度分布，从而模拟了流体的粘性和扩散等物理现象。通过对大量粒子的运动和相互作用进行统计平均，就可以得到流体的宏观物理量，如密度、速度和压力等。这种基于微观粒子动力学的模拟方式，使得LBM具有独特的优势。它能够直观地描述流体的微观行为，物理图像清晰，为理解复杂的流体现象提供了新的视角。其离散的计算方式天然适合并行计算，能够充分利用现代计算机的多核架构，大大提高计算效率。在处理复杂边界条件时，LBM也表现出较好的适应性，能够方便地处理不规则的几何边界。2.1.2LBM数学模型LBM的数学模型基于离散速度分布函数。在二维空间中，常用的离散速度模型为D2Q9模型，其中“D2”表示二维空间，“Q9”表示有9个离散速度方向。离散速度分布函数f_{i}(\mathbf{x},t)表示在时刻t、位置\mathbf{x}处具有速度\mathbf{c}_{i}的粒子分布密度，其中i=0,1,\cdots,8。LBM的演化方程可以表示为：f_{i}(\mathbf{x}+\mathbf{c}_{i}\Deltat,t+\Deltat)=f_{i}(\mathbf{x},t)+\Omega_{i}(\mathbf{x},t)其中，\Deltat为时间步长，\Omega_{i}(\mathbf{x},t)为碰撞算子，表示粒子在位置\mathbf{x}和时刻t的碰撞作用。碰撞算子通常采用BGK（Bhatnagar-Gross-Krook）近似，其形式为：\Omega_{i}(\mathbf{x},t)=-\frac{1}{\tau}[f_{i}(\mathbf{x},t)-f_{i}^{eq}(\mathbf{x},t)]其中，\tau为松弛时间，它与流体的运动粘度相关，\tau越大，流体的粘性越大；f_{i}^{eq}(\mathbf{x},t)为平衡态分布函数，其表达式为：f_{i}^{eq}(\rho,\mathbf{u})=\omega_{i}\rho\left[1+\frac{\mathbf{c}_{i}\cdot\mathbf{u}}{c_{s}^{2}}+\frac{(\mathbf{c}_{i}\cdot\mathbf{u})^{2}}{2c_{s}^{4}}-\frac{\mathbf{u}^{2}}{2c_{s}^{2}}\right]其中，\rho为流体密度，\mathbf{u}为流体速度，c_{s}为格子声速，在D2Q9模型中，c_{s}=\frac{1}{\sqrt{3}}；\omega_{i}为权重系数，对于D2Q9模型，\omega_{0}=\frac{4}{9}，\omega_{1-4}=\frac{1}{9}，\omega_{5-8}=\frac{1}{36}。通过对离散速度分布函数f_{i}(\mathbf{x},t)进行统计平均，可以得到流体的宏观物理量。流体密度\rho和速度\mathbf{u}的计算公式分别为：\rho=\sum_{i=0}^{8}f_{i}\rho\mathbf{u}=\sum_{i=0}^{8}\mathbf{c}_{i}f_{i}压力p可以通过状态方程p=c_{s}^{2}\rho计算得到。2.1.3LBM的特点与局限性LBM具有许多显著的特点，使其在计算流体力学领域得到广泛应用。首先，LBM具有高度的并行性。由于其计算过程主要是在每个格子节点上独立进行粒子的碰撞和迁移操作，不同节点之间的计算相互独立，因此非常适合并行计算。通过并行计算，可以大大缩短计算时间，提高计算效率，使其能够处理大规模的复杂流场问题。在模拟大型航空发动机内部复杂流场时，利用并行计算的LBM可以在较短时间内完成计算，为发动机的设计优化提供快速的数值模拟支持。其次，LBM的算法相对简单。它基于简单的粒子碰撞和迁移规则，不需要像传统计算流体力学方法那样对复杂的偏微分方程进行离散和求解，编程实现相对容易。这使得研究人员能够更快地开发和验证数值模型，降低了计算流体力学研究的门槛。对于一些非计算流体力学专业的研究人员，LBM简单的算法使其更容易上手，能够快速开展相关的研究工作。再者，LBM在处理复杂边界条件方面具有独特的优势。它可以通过简单的边界条件处理方法，如反弹边界条件、半反弹边界条件等，方便地处理各种不规则的几何边界。在模拟多孔介质中的流体流动时，LBM能够轻松处理多孔介质复杂的内部结构边界，准确地模拟流体在其中的流动特性。然而，LBM在模拟高速流动时也存在一些局限性。经典的LBM基于Boltzmann方程的低速近似，在高马赫数下，数值稳定性较差。随着马赫数的增加，LBM模拟结果中的数值振荡会逐渐加剧，导致计算结果不准确。经典LBM对激波等强间断的捕捉能力有限。激波是高速流动中的重要现象，其具有很强的非线性和间断性，经典LBM难以准确地捕捉激波的位置和强度，这限制了其在包含激波的高速流动模拟中的应用。在模拟高超声速飞行器绕流时，经典LBM无法准确模拟激波与边界层的相互作用，影响了对飞行器气动性能的准确评估。2.2NND格式2.2.1NND格式基本概念NND（Non-oscillatoryandNon-free-parameterDissipation）格式，即无波动无自由参数耗散格式，是一种在计算流体力学中具有重要地位的高精度差分格式。它的出现旨在有效解决传统差分格式在模拟包含激波等复杂流动现象时所面临的数值振荡和精度不足等问题。NND格式的核心思想是通过对通量的精细分裂和重构，实现对流动中各种物理量的准确捕捉和计算，特别是在激波附近，能够有效抑制数值解的非物理振荡，从而获得更加精确和稳定的数值结果。NND格式是由中心差分、二阶和一阶迎风格式混合组成的杂交型格式。中心差分格式在光滑流场区域能够提供较高的精度，它基于节点两侧的信息来计算导数，对于变化较为平缓的物理量分布能够准确地反映其变化趋势。在模拟无激波的均匀流场时，中心差分格式可以精确地计算速度、压力等物理量的梯度。然而，当遇到激波等强间断时，中心差分格式会产生数值振荡，因为它无法准确处理物理量在间断处的剧烈变化。二阶迎风格式则根据流动方向，利用上游节点的信息来计算导数。在超声速流动中，气流的信息是沿着流动方向传播的，二阶迎风格式能够更好地捕捉这种信息传播的特性，对于激波等具有方向性的流动结构有较好的适应性。它在一定程度上可以减少激波附近的数值振荡，但在某些情况下，仍可能存在一定的数值耗散，导致激波的分辨率不够高。一阶迎风格式是一种更为简单的迎风差分格式，它只考虑上游一个节点的信息。虽然一阶迎风格式的精度相对较低，但它具有较强的稳定性，能够有效地抑制数值振荡。在处理强间断问题时，一阶迎风格式可以提供一个稳定的基础，防止数值解出现严重的振荡和发散。NND格式巧妙地将这三种格式结合起来，充分发挥它们各自的优势。在光滑流场区域，主要采用中心差分格式来保证计算精度；在激波等间断附近，通过引入二阶和一阶迎风格式，根据流动的具体情况自动调整格式的权重，从而有效地抑制数值振荡，提高激波的捕捉能力。这种混合格式的设计，使得NND格式在模拟复杂流场时具有更高的精度和稳定性，能够准确地捕捉到激波的位置、强度以及其他复杂的流动结构。2.2.2NND格式的数值实现在数值实现方面，NND格式首先需要对控制方程进行离散化。以二维可压缩流体的Navier-Stokes方程为例，其守恒形式为：\frac{\partial\mathbf{U}}{\partialt}+\frac{\partial\mathbf{E}}{\partialx}+\frac{\partial\mathbf{F}}{\partialy}=\mathbf{S}其中，\mathbf{U}为守恒变量向量，\mathbf{E}和\mathbf{F}分别为x和y方向的通量向量，\mathbf{S}为源项向量。对于空间离散，NND格式采用有限差分法对通量进行计算。以x方向的通量\mathbf{E}为例，其离散形式可以表示为：\mathbf{E}_{i+\frac{1}{2},j}=\mathbf{E}_{i+\frac{1}{2},j}^{C}+\mathbf{E}_{i+\frac{1}{2},j}^{U}其中，\mathbf{E}_{i+\frac{1}{2},j}^{C}为中心差分部分的通量，\mathbf{E}_{i+\frac{1}{2},j}^{U}为迎风格式部分的通量。中心差分部分的通量\mathbf{E}_{i+\frac{1}{2},j}^{C}通常采用二阶中心差分计算，即：\mathbf{E}_{i+\frac{1}{2},j}^{C}=\frac{1}{2}(\mathbf{E}_{i+1,j}+\mathbf{E}_{i,j})迎风格式部分的通量\mathbf{E}_{i+\frac{1}{2},j}^{U}则根据流动方向和当地的马赫数等参数来确定。具体来说，首先需要对通量进行分裂，将其分为正向通量和反向通量。以Roe通量分裂为例，\mathbf{E}可以分裂为\mathbf{E}^{+}和\mathbf{E}^{-}，满足\mathbf{E}=\mathbf{E}^{+}+\mathbf{E}^{-}。然后，根据流动方向，选择合适的迎风模板来计算迎风格式部分的通量。在超声速气流从左向右流动的情况下，如果节点(i,j)处的气流速度为正，则\mathbf{E}_{i+\frac{1}{2},j}^{U}主要由\mathbf{E}_{i,j}^{+}和\mathbf{E}_{i+1,j}^{-}组成，通过适当的权重分配来计算。在时间推进方面，NND格式通常采用显式或隐式的时间积分方法。显式时间积分方法，如Runge-Kutta法，计算简单，易于实现，但时间步长受到稳定性条件的限制，对于复杂的大规模计算可能计算效率较低。隐式时间积分方法，如隐式向后差分法，虽然可以采用较大的时间步长，提高计算效率，但需要求解大型的线性方程组，计算复杂度较高。在实际应用中，需要根据具体问题的特点和计算资源的限制，选择合适的时间积分方法。2.2.3NND格式在高速流动模拟中的优势NND格式在高速流动模拟中具有多方面的显著优势。首先，在捕捉激波方面，NND格式表现出色。激波是高速流动中的重要现象，其具有很强的非线性和间断性，准确捕捉激波的位置和强度对于理解高速流动特性至关重要。传统的差分格式在激波附近容易产生数值振荡，导致激波的分辨率降低，无法准确反映激波的真实特性。而NND格式通过其独特的通量分裂和重构方式，能够有效地抑制激波附近的数值振荡，清晰地捕捉到激波的位置和强度。在模拟高超声速飞行器绕流时，NND格式可以准确地确定激波在飞行器表面的附着位置和激波的强度分布，为飞行器的气动设计提供关键数据。其次，NND格式能够提高计算精度。在高速流动中，物理量的变化往往较为剧烈，对计算精度要求很高。NND格式结合了中心差分、二阶和一阶迎风格式的优点，在光滑流场区域采用高精度的中心差分格式，在激波等间断区域通过合理的迎风差分来保证计算的稳定性和准确性，从而在整个流场中都能提供较高的计算精度。与传统的一阶迎风格式相比，NND格式在计算复杂流场的压力、速度等物理量时，能够更准确地反映其分布和变化，减少数值误差。NND格式还具有良好的稳定性。在高速流动模拟中，由于流动的复杂性和非线性，数值稳定性是一个关键问题。不稳定的数值计算可能导致计算结果发散，无法得到有效的物理信息。NND格式通过对通量的合理处理和格式的优化，增强了数值计算的稳定性，能够在较大的马赫数范围内保持稳定的计算。在模拟高马赫数的激波管流动时，NND格式能够稳定地计算激波的传播和相互作用过程，即使在长时间的计算中也不会出现数值不稳定的情况。三、基于NND格式的格子波兹曼方法改进3.1现有LBM在高速流动模拟中的问题分析经典的格子波兹曼方法（LBM）在模拟低速、低马赫数的流体流动时表现出良好的性能，但在处理高速流动问题时，存在一些明显的局限性，限制了其在高速流动模拟领域的广泛应用。在高马赫数下，经典LBM的数值稳定性较差。随着马赫数的增加，流体的可压缩性效应变得显著，流动中的非线性因素增强。经典LBM基于Boltzmann方程的低速近似，其模型中的一些假设在高速流动情况下不再成立。在高马赫数下，流体的密度、压力等物理量的变化更加剧烈，经典LBM的碰撞算子和演化方程难以准确地描述这些变化，导致数值解出现不稳定的振荡现象。当马赫数达到3以上时，经典LBM模拟结果中的速度和压力分布会出现明显的波动，无法收敛到稳定的解，严重影响了模拟结果的准确性和可靠性。经典LBM对激波的捕捉能力有限。激波是高速流动中的重要特征，它是一种强间断现象，伴随着物理量的剧烈变化。经典LBM采用的离散速度模型和碰撞算子在处理激波时存在缺陷，难以准确地捕捉激波的位置和强度。在激波附近，经典LBM的数值解会出现过度扩散或振荡的情况，导致激波的分辨率降低，无法真实地反映激波的特性。在模拟激波管中的激波传播时，经典LBM得到的激波位置与理论值存在较大偏差，激波的强度也被低估，无法准确地模拟激波与边界层的相互作用等复杂现象。经典LBM在模拟高速流动时的计算精度也有待提高。在高速流动中，物理量的变化梯度较大，对计算精度要求更高。经典LBM的离散格式和数值算法在处理这些大梯度变化时，容易引入数值误差，导致计算精度下降。在模拟高超声速飞行器绕流时，经典LBM对飞行器表面的压力分布和摩擦力系数的计算结果与实验数据存在较大误差，无法满足工程设计对精度的要求。3.2Qu模型及改进思路3.2.1Qu模型介绍Qu模型是一种为解决高速流动模拟问题而提出的改进型格子波兹曼模型，其核心在于将一阶时间和二阶空间准确的NND有限差分法巧妙地应用于经典格子波尔兹曼方程。在传统的格子波兹曼方法中，离散速度分布函数的演化方程主要由粒子的迁移和碰撞过程决定。而Qu模型在此基础上，引入NND有限差分法来处理方程中的对流项和扩散项，以提高对高速流动中复杂物理现象的模拟能力。具体而言，对于格子波兹曼方程中的对流项，传统方法通常采用简单的离散格式，在高速流动情况下，这种简单格式难以准确描述流体的对流运动，导致数值误差较大。Qu模型运用NND有限差分法，通过对通量的精确计算和分裂，能够更准确地捕捉流体的对流信息。在处理超声速气流的对流时，NND有限差分法可以根据气流的速度方向和大小，合理地分配通量，从而更精确地模拟气流的运动轨迹和速度变化。在扩散项的处理上，Qu模型同样利用NND有限差分法的高精度特性，对扩散过程进行更细致的模拟。在高速流动中，流体的扩散现象往往与对流相互耦合，传统的格子波兹曼模型难以准确刻画这种复杂的耦合关系。Qu模型通过NND有限差分法，能够更好地处理扩散项中的空间导数，从而提高对扩散过程的模拟精度。这种将NND有限差分法应用于经典格子波尔兹曼方程的方式，使得Qu模型在高速流动模拟中展现出独特的优势。它能够有效地抑制数值振荡，提高模拟结果的稳定性和准确性，特别是在处理激波等强间断现象时，表现出比传统格子波兹曼模型更强的捕捉能力。在模拟激波管中的激波传播时，Qu模型能够准确地捕捉到激波的位置和强度变化，其模拟结果与理论解和实验数据的吻合度更高。3.2.2基于Wang的NND模式的改进在Qu模型的基础上，进一步采用Wang的NND模式对经典算法进行改进，旨在进一步提升算法在二阶时间和空间精度上的性能。Wang的NND模式在通量计算和差分格式的选择上具有独特的优势，能够更精确地处理高速流动中的复杂物理过程。在通量计算方面，Wang的NND模式通过对通量的更细致分裂和重构，提高了计算的精度。传统的NND格式在通量分裂时，可能会忽略一些高阶项的影响，导致在处理高速流动中物理量的剧烈变化时精度不足。Wang的NND模式考虑了更多的高阶项信息，通过对这些高阶项的合理处理，能够更准确地反映物理量在空间和时间上的变化趋势。在计算激波附近的通量时，Wang的NND模式能够更精确地捕捉到激波前后物理量的突变，从而提高对激波的分辨率。在差分格式的选择上，Wang的NND模式针对二阶时间和空间精度进行了优化。在时间离散上，采用了更精确的二阶时间差分格式，相比于传统的一阶时间差分格式，能够更好地捕捉物理量随时间的变化细节。在空间离散上，通过合理地组合中心差分和迎风格式，根据流场的具体情况自动调整格式的权重，从而在保证计算稳定性的提高空间分辨率。在光滑流场区域，增加中心差分格式的权重，以提高计算精度；在激波等间断区域，加大迎风格式的权重，以抑制数值振荡。通过使用Wang的NND模式，改进后的算法在处理高速流动问题时，能够更准确地模拟物理量的分布和变化，提高了对激波、边界层等复杂流动结构的捕捉能力。在模拟高超声速飞行器绕流时，改进后的算法能够更清晰地分辨出激波与边界层的相互作用，准确地计算出飞行器表面的压力分布和摩擦力系数，为飞行器的气动设计提供更可靠的数据支持。3.3改进后模型的理论推导与分析改进后的基于NND格式的格子波兹曼模型在理论推导方面有着独特的过程和优势，其在时间和空间精度、稳定性等方面与原模型存在显著差异，这些差异决定了改进后模型在高速流动模拟中的性能提升。从理论推导来看，改进后的模型首先对格子波兹曼方程中的对流项和扩散项进行了重新处理。以二维流动为例，在原格子波兹曼方程中，对流项的离散通常采用简单的一阶迎风格式或中心差分格式，这在高速流动中难以准确描述流体的对流运动。改进后的模型运用Wang的NND模式，对对流项进行了更精确的离散。对于对流项\frac{\partial(\rhou_{i}u_{j})}{\partialx_{j}}（其中i,j=1,2表示空间方向），在原模型中，若采用简单的一阶迎风格式，其离散形式可能为：\left(\frac{\partial(\rhou_{i}u_{j})}{\partialx_{j}}\right)_{i,j}^{n}\approx\frac{(\rhou_{i}u_{j})_{i,j}^{n}-(\rhou_{i}u_{j})_{i-1,j}^{n}}{\Deltax}\quad(u_{j}\geq0)\left(\frac{\partial(\rhou_{i}u_{j})}{\partialx_{j}}\right)_{i,j}^{n}\approx\frac{(\rhou_{i}u_{j})_{i+1,j}^{n}-(\rhou_{i}u_{j})_{i,j}^{n}}{\Deltax}\quad(u_{j}\lt0)这种简单的迎风格式在激波等间断区域会产生较大的数值耗散，导致激波的分辨率降低。而在改进后的模型中，运用Wang的NND模式，首先对通量\rhou_{i}u_{j}进行分裂，将其分为正向通量和反向通量。以Roe通量分裂为例，\rhou_{i}u_{j}可以分裂为(\rhou_{i}u_{j})^{+}和(\rhou_{i}u_{j})^{-}，满足\rhou_{i}u_{j}=(\rhou_{i}u_{j})^{+}+(\rhou_{i}u_{j})^{-}。然后，根据流动方向和当地的马赫数等参数，采用不同的差分模板来计算通量。在超声速气流从左向右流动的情况下，对于节点(i,j)处的对流项离散，其表达式可能为：\left(\frac{\partial(\rhou_{i}u_{j})}{\partialx_{j}}\right)_{i,j}^{n}\approx\frac{(\rhou_{i}u_{j})^{+}_{i,j}^{n}-(\rhou_{i}u_{j})^{+}_{i-1,j}^{n}}{\Deltax}+\frac{(\rhou_{i}u_{j})^{-}_{i+1,j}^{n}-(\rhou_{i}u_{j})^{-}_{i,j}^{n}}{\Deltax}其中，(\rhou_{i}u_{j})^{+}_{i,j}^{n}和(\rhou_{i}u_{j})^{-}_{i,j}^{n}的计算会考虑到更多的高阶项信息，通过对这些高阶项的合理处理，能够更准确地反映物理量在空间和时间上的变化趋势。在计算(\rhou_{i}u_{j})^{+}_{i,j}^{n}时，会考虑到节点(i,j)及其相邻节点的速度、密度等信息，通过复杂的公式计算得到更精确的正向通量值。在扩散项的处理上，原模型的离散方式也存在一定的局限性。对于扩散项\frac{\partial(\mu\frac{\partialu_{i}}{\partialx_{j}})}{\partialx_{j}}（\mu为动力粘度），原模型的离散可能无法准确描述扩散过程中的非线性特性。改进后的模型利用Wang的NND模式，对扩散项中的空间导数进行了更精细的计算。通过引入高阶差分格式，能够更好地捕捉扩散过程中物理量的变化细节。在计算\frac{\partial(\mu\frac{\partialu_{i}}{\partialx_{j}})}{\partialx_{j}}时，采用四阶中心差分格式来计算\frac{\partialu_{i}}{\partialx_{j}}，再进行后续的计算，从而提高了对扩散过程的模拟精度。从时间和空间精度方面分析，原经典格子波兹曼模型在时间精度上通常为一阶，空间精度为二阶。这意味着在时间推进过程中，其对物理量随时间变化的描述相对粗糙，在空间离散上，对于复杂流场中物理量的梯度变化捕捉能力有限。而改进后的模型采用Wang的NND模式，在时间精度上提升到了二阶。这使得改进后的模型在时间推进过程中，能够更准确地捕捉物理量随时间的变化趋势。在模拟激波的传播过程时，二阶时间精度可以更精确地计算激波在不同时刻的位置和强度变化，减少时间离散带来的误差。在空间精度上，改进后的模型通过合理地组合中心差分和迎风格式，根据流场的具体情况自动调整格式的权重，进一步提高了空间分辨率。在光滑流场区域，增加中心差分格式的权重，以提高计算精度；在激波等间断区域，加大迎风格式的权重，以抑制数值振荡。这使得改进后的模型在整个流场中都能提供更高的计算精度，更准确地模拟物理量的分布和变化。在稳定性方面，原经典格子波兹曼模型在高马赫数下数值稳定性较差，随着马赫数的增加，数值振荡逐渐加剧，导致计算结果不准确。这是因为原模型基于Boltzmann方程的低速近似，在高马赫数下，其假设不再成立，模型难以准确描述流体的可压缩性效应和非线性因素。改进后的模型通过引入NND格式，增强了数值计算的稳定性。NND格式通过对通量的合理处理和格式的优化，有效地抑制了数值振荡。在高马赫数流动模拟中，改进后的模型能够稳定地计算激波的传播和相互作用过程，即使在长时间的计算中也不会出现数值不稳定的情况。这是因为改进后的模型在处理通量时，考虑了更多的物理因素，通过对通量的精确分裂和重构，减少了数值误差的积累，从而提高了稳定性。与原模型相比，改进后的模型在理论和性能上都有显著的差异和提升。在理论上，改进后的模型对对流项和扩散项的处理更加精细，考虑了更多的物理因素和高阶项信息，使得模型的理论基础更加完善。在性能上，改进后的模型在时间和空间精度、稳定性以及激波捕捉能力等方面都有明显的提高。这些优势使得改进后的模型更适合用于高速流动的模拟，能够为相关工程领域提供更准确、可靠的数值模拟结果。四、数值试验与结果验证4.1一维数值试验4.1.1试验设置为了验证改进后的基于NND格式的格子波兹曼模型在高速流动模拟中的有效性，进行了一维激波管数值试验。一维激波管是计算流体力学中的经典问题，能够很好地检验数值方法对激波等复杂流动现象的捕捉能力。试验中，激波管的长度设定为L=1.0，在x=0.5处设置一个薄膜，将激波管分为左右两个区域。初始时刻，左区域的气体状态设定为：密度\rho_{L}=1.0，速度u_{L}=0.0，压力p_{L}=1.0；右区域的气体状态设定为：密度\rho_{R}=0.125，速度u_{R}=0.0，压力p_{R}=0.1。气体的绝热指数\gamma=1.4。在边界条件的设置上，激波管的两端均采用非反射边界条件，以模拟气体的自由出流。这种边界条件能够有效地避免反射波对计算结果的影响，保证激波在管内的传播过程不受干扰。在数值计算过程中，时间步长\Deltat设定为0.001，空间步长\Deltax设定为0.002。这些参数的选择经过了多次调试和验证，能够在保证计算精度的前提下，提高计算效率。同时，为了确保计算结果的准确性，对计算区域进行了足够的网格划分，以满足数值计算的精度要求。4.1.2结果与分析利用改进后的模型对一维激波管问题进行数值模拟，并将计算结果与经典LBM以及精确解进行对比分析。图1展示了t=0.2时刻改进模型、经典LBM和精确解的密度分布对比。从图中可以清晰地看出，经典LBM的计算结果在激波附近出现了明显的数值振荡，激波的位置和强度与精确解存在较大偏差。经典LBM得到的激波位置比精确解偏右，激波强度也被低估，导致密度分布在激波附近出现了不合理的波动。这是由于经典LBM在处理激波等强间断时，其离散格式和碰撞算子无法准确地捕捉激波的特性，容易产生数值耗散和振荡。相比之下，改进后的模型能够准确地捕捉激波的位置和强度，计算结果与精确解吻合良好。在激波附近，改进后的模型没有出现明显的数值振荡，密度分布曲线光滑，能够真实地反映激波的物理特性。这得益于改进后的模型引入了NND格式，通过对通量的精确分裂和重构，有效地抑制了数值振荡，提高了激波的分辨率。在计算通量时，改进后的模型考虑了更多的高阶项信息，能够更准确地反映物理量在激波附近的变化趋势，从而提高了对激波的捕捉能力。图2展示了t=0.2时刻改进模型、经典LBM和精确解的速度分布对比。同样，经典LBM的计算结果在激波附近存在较大误差，速度分布出现了不合理的波动。而改进后的模型计算结果与精确解基本一致，能够准确地模拟激波传播过程中速度的变化。在激波后的区域，改进后的模型计算得到的速度值与精确解相符，准确地反映了激波对气体速度的影响。通过对一维激波管问题的数值模拟和结果分析，可以得出改进后的基于NND格式的格子波兹曼模型在高速流动模拟中具有更高的精度和更好的激波捕捉能力，能够有效地克服经典LBM在模拟高速流动时的局限性。4.2二维数值试验4.2.1试验设置为了进一步验证改进后的基于NND格式的格子波兹曼模型在复杂二维流场模拟中的性能，进行了二维圆柱绕流数值试验。二维圆柱绕流是流体力学中的经典问题，包含了多种复杂的流动现象，如边界层分离、尾流涡街的形成等，对数值方法的精度和稳定性要求较高。在几何模型方面，设置一个直径D=0.1的圆柱体，放置在一个长L=2.0、宽W=1.0的矩形计算域中心。圆柱体的中心坐标为(1.0,0.5)。计算域的入口位于x=0处，出口位于x=2.0处，上下边界分别位于y=0和y=1.0处。在网格划分上，采用非结构化三角形网格对计算域进行离散。为了准确捕捉圆柱体表面和近壁区域的流动细节，在圆柱体周围进行了网格加密。近壁区域的最小网格尺寸设定为\Deltax_{min}=\Deltay_{min}=0.001，随着远离圆柱体，网格尺寸逐渐增大，最大网格尺寸设定为\Deltax_{max}=\Deltay_{max}=0.02。通过这种网格划分方式，既能保证在关键区域的计算精度，又能控制计算量在合理范围内。在进行网格质量检查时，确保网格的最小内角大于30^{\circ}，最大纵横比小于10，以保证数值计算的稳定性。在参数设置方面，来流马赫数Ma设定为0.8，对应超声速流动。来流速度u_{\infty}根据马赫数和当地声速计算得到，假设来流温度T_{\infty}=300K，根据理想气体状态方程p=\rhoRT（其中R为气体常数，对于空气R=287J/(kg\cdotK)）和声速公式c=\sqrt{\gammaRT}（\gamma=1.4为空气的绝热指数），可计算得到当地声速c_{\infty}=\sqrt{1.4\times287\times300}\approx347m/s，则来流速度u_{\infty}=Ma\timesc_{\infty}=0.8\times347=277.6m/s。来流密度\rho_{\infty}=1.225kg/m^{3}，压力p_{\infty}=101325Pa。雷诺数Re根据来流速度、圆柱直径和运动粘度计算，假设运动粘度\nu=1.5\times10^{-5}m^{2}/s，则Re=\frac{u_{\infty}D}{\nu}=\frac{277.6\times0.1}{1.5\times10^{-5}}\approx1.85\times10^{6}。在边界条件的设置上，入口采用速度入口边界条件，给定来流速度、密度和温度；出口采用压力出口边界条件，给定出口压力；圆柱体表面采用无滑移壁面边界条件，即流体在壁面上的速度为零；上下边界采用对称边界条件，以模拟无限宽流场的情况。在时间步长的选择上，通过稳定性分析和多次调试，设定时间步长\Deltat=1\times10^{-6}s，以确保数值计算的稳定性。4.2.2结果与分析利用改进后的模型对二维圆柱绕流问题进行数值模拟，并将计算结果与经典LBM以及实验数据进行对比分析。图3展示了改进模型、经典LBM在t=0.1s时刻的流场马赫数云图对比。从图中可以明显看出，经典LBM的计算结果在圆柱体下游的尾流区域出现了较大的数值振荡，马赫数分布不均匀，无法准确地反映尾流的真实结构。经典LBM得到的尾流区域马赫数波动较大，存在不合理的高马赫数和低马赫数区域，这是由于经典LBM在处理边界层分离和尾流涡街等复杂流动现象时，其离散格式和碰撞算子无法准确地捕捉流动特性，导致数值解出现不稳定的振荡。相比之下，改进后的模型能够清晰地捕捉到圆柱绕流的流场结构。在圆柱体表面，改进后的模型准确地模拟了边界层的形成和发展，边界层厚度与理论值相符。在圆柱体下游的尾流区域，改进后的模型能够准确地捕捉到尾流涡街的形成和发展，涡街的形状、间距和脱落频率与实验数据吻合良好。在涡街的脱落过程中，改进后的模型计算得到的涡街间距为0.2D，与实验测量的0.21D非常接近，这得益于改进后的模型引入了NND格式，通过对通量的精确分裂和重构，有效地抑制了数值振荡，提高了对复杂流场结构的捕捉能力。图4展示了改进模型、经典LBM在t=0.1s时刻沿圆柱表面的压力系数分布对比。压力系数C_p定义为C_p=\frac{p-p_{\infty}}{\frac{1}{2}\rho_{\infty}u_{\infty}^{2}}，其中p为当地压力。从图中可以看出，经典LBM的计算结果在圆柱表面的压力系数分布与实验数据存在较大偏差，特别是在圆柱的背风面，压力系数的计算值与实验值相差较大。经典LBM在圆柱背风面得到的压力系数波动较大，无法准确地反映压力的真实分布，这是因为经典LBM在处理边界层分离后的复杂流动时，无法准确地计算压力的变化。而改进后的模型计算得到的压力系数分布与实验数据基本一致，能够准确地反映圆柱表面的压力变化。在圆柱的迎风面，改进后的模型计算得到的压力系数与实验值相符，准确地模拟了气流的冲击作用。在圆柱的背风面，改进后的模型能够准确地捕捉到压力系数的变化趋势，压力系数的最小值和位置与实验数据吻合，这表明改进后的模型能够更准确地模拟边界层分离后的压力分布。通过计算阻力系数C_D来进一步评估改进模型的准确性，阻力系数C_D定义为C_D=\frac{F_D}{\frac{1}{2}\rho_{\infty}u_{\infty}^{2}D}，其中F_D为作用在圆柱体上的阻力。改进后的模型计算得到的阻力系数为C_D=1.25，与实验测量值1.28较为接近，相对误差为2.3\%。而经典LBM计算得到的阻力系数为1.42，与实验值的相对误差为10.9\%。这表明改进后的模型在计算阻力系数方面具有更高的精度，能够更准确地预测作用在圆柱体上的阻力。通过对二维圆柱绕流问题的数值模拟和结果分析，可以得出改进后的基于NND格式的格子波兹曼模型在复杂二维流场模拟中具有更高的精度和更好的稳定性，能够有效地克服经典LBM在模拟复杂流场时的局限性，准确地捕捉流场结构、压力分布和阻力系数等关键物理量。4.3与其他算法的对比验证为了更全面地评估改进后的基于NND格式的格子波兹曼模型的性能，选择了其他常用于高速流动模拟的算法进行对比验证，主要从计算精度、效率和稳定性等方面展开分析。在计算精度方面，将改进后的模型与传统的有限体积法（FVM）和有限差分法（FDM）进行对比。以二维圆柱绕流问题为例，在相同的计算条件下，分别使用改进后的模型、FVM和FDM进行数值模拟。通过比较模拟得到的圆柱表面压力系数分布和尾流区域的速度分布等关键物理量与实验数据的吻合程度，来评估各算法的计算精度。从圆柱表面压力系数分布的对比结果来看，FVM在计算时，由于其对通量的计算方式，在圆柱的背风面等压力变化剧烈的区域，计算结果与实验数据存在一定偏差。FVM在计算背风面压力系数时，由于数值耗散的影响，压力系数的最小值与实验值相比偏低，且压力系数的变化趋势在某些位置也与实验数据不完全一致。FDM在处理复杂边界条件时，其离散格式的局限性导致在圆柱表面附近的计算精度受限，压力系数分布在边界层区域与实验数据的差异较大。相比之下，改进后的基于NND格式的格子波兹曼模型能够更准确地捕捉圆柱表面的压力变化，压力系数分布与实验数据的吻合度更高。在尾流区域的速度分布对比中，FVM和FDM在捕捉尾流涡街的特性时，都存在一定的误差。FVM得到的涡街间距与实验值相比偏大，涡街的强度也有所减弱；FDM则在涡街的位置和形状上与实验数据存在偏差。而改进后的模型能够准确地模拟尾流涡街的形成、发展和脱落过程，涡街的形状、间距和脱落频率与实验数据相符，充分体现了其在计算精度上的优势。在计算效率方面，对比改进后的模型与基于GPU加速的有限元法（FEM）。在模拟大型航空发动机内部复杂流场时，分别采用改进后的模型和基于GPU加速的FEM进行计算。记录两种算法在相同计算规模下的计算时间，以此来评估它们的计算效率。基于GPU加速的FEM虽然利用了GPU的并行计算能力，但由于其计算过程中需要求解大规模的线性方程组，计算复杂度较高，导致计算时间较长。在模拟一个包含数百万个计算单元的航空发动机流场时，基于GPU加速的FEM计算时间达到了数小时。而改进后的基于NND格式的格子波兹曼模型具有天然的并行性，其计算过程主要是在每个格子节点上独立进行粒子的碰撞和迁移操作，不同节点之间的计算相互独立，能够充分利用现代计算机的多核计算资源。在相同的计算规模下，改进后的模型计算时间仅为数十分钟，大大提高了计算效率。改进后的模型在内存占用方面也相对较低，对于大规模计算具有更好的适应性。在稳定性方面，将改进后的模型与高精度加权本质无振荡格式（WENO）进行对比。以高马赫数激波管问题为例，在马赫数为5的情况下，分别使用改进后的模型和WENO格式进行数值模拟。观察两种算法在长时间计算过程中的收敛情况和数值振荡情况，以评估它们的稳定性。WENO格式在高马赫数下，虽然具有较好的激波捕捉能力，但由于其对数值耗散的控制方式，在长时间计算时，容易出现数值振荡逐渐放大的情况，导致计算结果不稳定。在激波管问题模拟中，随着计算时间的增加，WENO格式计算得到的激波位置和强度出现了较大的波动，无法收敛到稳定的解。而改进后的基于NND格式的格子波兹曼模型通过对通量的合理处理和格式的优化，有效地抑制了数值振荡。在相同的高马赫数条件下，改进后的模型能够稳定地计算激波的传播和相互作用过程，即使在长时间的计算中，激波的位置和强度也能保持稳定，计算结果准确可靠。通过与其他算法在计算精度、效率和稳定性等方面的对比验证，可以得出改进后的基于NND格式的格子波兹曼模型在高速流动模拟中具有更优的综合性能，能够为相关工程领域提供更准确、高效和稳定的数值模拟结果。五、高速流动中激波相互作用的模拟5.1激波相互作用的物理现象与理论基础激波是高速流动中一种极为重要的物理现象，它的产生与流体的可压缩性以及速度变化密切相关。当流体的流速超过当地声速时，就会产生激波。在激波的传播过程中，其波面两侧的流体参数，如压力、密度、温度和速度等，会发生急剧的变化，这种变化是在极短的距离内完成的，呈现出一种强间断的特性。在超声速飞行器的飞行过程中，飞行器头部的气流受到强烈压缩，流速急剧降低，压力和温度急剧升高，从而形成激波。这种激波的产生会对飞行器的气动性能产生重大影响，例如增加飞行器的阻力，影响飞行器的稳定性和操纵性。激波的传播过程遵循一定的物理规律。在理想气体中，激波的传播速度与气体的性质、初始状态以及激波的强度等因素有关。根据激波的传播方向与气流方向的关系，激波可分为正激波和斜激波。正激波是指波面与气流方向垂直的激波，在正激波的作用下，气流通过激波后，速度会突然降低，压力、密度和温度会突然升高。斜激波则是波面与气流方向成一定角度的激波，斜激波的特性比正激波更为复杂，其强度和波角与气流的马赫数、偏转角等因素有关。在二维超声速气流流过楔形物体时，在楔形物体的头部会形成斜激波，斜激波的波角和强度会随着气流马赫数和楔形物体顶角的变化而变化。当多个激波在流场中相遇时，就会发生激波相互作用。激波相互作用的形式多种多样，主要包括正激波相互作用和斜激波相互作用。正激波相互作用通常发生在两个或多个正交于气流方向的激波之间。当两个正激波相互靠近并相遇时，它们会相互干扰，形成复杂的波系结构。在某些情况下，两个正激波可能会合并成一个更强的激波，这种合并过程会导致波后流体参数的进一步变化。在激波管实验中，当两个相向传播的正激波相遇时，它们会相互碰撞，合并成一个更强的激波，波后流体的压力和温度会急剧升高。斜激波相互作用则发生在斜向于气流方向的激波之间。斜激波相互作用的情况更为复杂，除了激波的合并，还可能出现激波的反射和折射现象。当斜激波遇到障碍物或其他激波时，会发生反射，反射激波的强度和方向与入射激波的特性以及障碍物的形状和位置有关。斜激波在不同介质的界面上传播时，还会发生折射，折射激波的传播方向和强度会根据介质的性质和界面的特性而变化。在高超声速气流流过具有复杂外形的飞行器时，飞行器表面的不同部位会产生多个斜激波，这些斜激波在传播过程中会相互作用，出现反射和折射现象，使得流场中的压力、温度和速度分布变得极为复杂。描述激波相互作用的理论基础主要包括流体力学的基本守恒方程，即连续性方程、动量方程和能量方程。连续性方程描述了流体在流动过程中质量守恒的特性，其表达式为：\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\mathbf{u})=0其中，\rho为流体密度，\mathbf{u}为流体速度矢量，t为时间。动量方程描述了流体动量的守恒，在无粘性流体中，动量方程简化为欧拉方程，其表达式为：\frac{\partial(\rho\mathbf{u})}{\partialt}+\nabla\cdot(\rho\mathbf{u}\mathbf{u})=-\nablap其中，p为流体压力。能量方程描述了流体能量的守恒，包括动能、内能等，其表达式为：\frac{\partial(\rhoe)}{\partialt}+\nabla\cdot(\rhoe\mathbf{u})=-\nabla\cdot(p\mathbf{u})其中，e为单位质量流体的总能量。在激波相互作用的研究中，还需要考虑激波的间断条件，即兰金-胡克尼奥（Rankine-Hugoniot）关系。该关系基于质量、动量和能量守恒定律，描述了激波前后流体参数的变化关系。对于正激波，兰金-胡克尼奥关系可以表示为一系列关于激波前后压力、密度、温度和速度等参数的等式。设激波前的参数为p_1、\rho_1、T_1、u_1，激波后的参数为p_2、\rho_2、T_2、u_2，则有：\rho_1u_1=\rho_2u_2p_1+\rho_1u_1^2=p_2+\rho_2u_2^2h_1+\frac{u_1^2}{2}=h_2+\frac{u_2^2}{2}其中，h为单位质量流体的焓。这些关系对于理解激波相互作用过程中流体参数的突变和守恒特性至关重要，为数值模拟和理论分析提供了重要的依据。5.2基于改进模型的激波模拟方法利用改进后的格子波兹曼-NND模型模拟激波相互作用时，需遵循特定的步骤和实现方法，以确保模拟的准确性和可靠性。首先，要对计算区域进行合理的离散化处理。对于二维激波相互作用问题，可采用结构化或非结构化网格对计算区域进行划分。若采用结构化网格，需根据激波的传播方向和相互作用区域的特点，合理确定网格的疏密程度。在激波可能发生相互作用的区域，如激波碰撞点附近，应加密网格，以提高对激波细节的捕捉能力。在一个包含两个斜激波相互作用的二维流场中，在激波相互作用区域，将网格尺寸设置为\Deltax=\Deltay=0.001，而在远离相互作用区域的流场其他部分，网格尺寸可适当增大至\Deltax=\Deltay=0.01。若采用非结构化网格，则可根据流场的复杂程度自动生成网格，在关键区域实现网格的自适应加密。在模拟高超声速飞行器表面的激波相互作用时，非结构化网格能够更好地贴合飞行器复杂的外形，在激波与飞行器表面相互作用的区域，自动加密网格，提高计算精度。接着，需对初始条件和边界条件进行精确设定。在初始条件方面，要根据激波的类型和产生机制，准确给定流场中各物理量的初始值。对于正激波相互作用的模拟，需给定激波前后的压力、密度、速度等物理量。假设在一个一维激波管中，初始时刻在x=0.5处存在一个正激波，激波前的压力p_1=1.0，密度\rho_1=1.0，速度u_1=0.0；激波后的压力p_2=10.0，密度\rho_2=6.16，速度u_2=2.69（根据激波关系计算得到）。在边界条件的设定上，要根据实际问题的物理背景，选择合适的边界条件。对于计算区域的入口和出口，可采用非反射边界条件，以避免反射波对激波相互作用的影响。在模拟超声速风洞中的激波相互作用时，风洞入口给定超声速来流的速度、压力和密度等参数，出口采用非反射边界条件，确保激波在传播过程中不受出口边界的干扰。对于固体壁面边界，可采用无滑移边界条件，即流体在壁面上的速度为零。在模拟激波与平板相互作用时，平板表面采用无滑移边界条件，准确模拟激波与壁面的相互作用过程。在数值计算过程中，要按照改进后的格子波兹曼-NND模型的演化方程进行计算。在每个时间步，首先进行粒子的迁移操作，根据离散速度分布函数f_{i}(\mathbf{x},t)和离散速度\mathbf{c}_{i}，计算粒子在迁移后的分布函数f_{i}(\mathbf{x}+\mathbf{c}_{i}\Deltat,t+\Deltat)。然后，进行碰撞操作，利用改进后的碰撞算子，如基于Wang的NND模式改进后的碰撞算子，计算碰撞后的分布函数。在碰撞算子的计算中，要考虑到通量的精确分裂和重构，以准确模拟激波相互作用过程中的物理现象。在计算激波相互作用区域的通量时，采用Roe通量分裂方法，将通量分为正向通量和反向通量，根据流动方向和当地的马赫数等参数，精确计算通量，从而提高对激波的捕捉能力。通过对分布函数的统计平均，得到流场的宏观物理量，如密度、速度和压力等。在模拟过程中，还需对计算结果进行实时监测和分析。可通过设置监测点，实时记录监测点处的物理量随时间的变化情况。在激波相互作用区域设置多个监测点，记录压力、密度等物理量的变化，以分析激波相互作用的过程和特征。利用可视化技术，如绘制马赫数云图、压力云图等，直观地展示流场中激波的传播和相互作用情况。通过可视化分析，能够更清晰地观察激波的形状、位置和相互作用的结果，为进一步的分析和研究提供依据。5.3模拟结果与分析利用改进后的格子波兹曼-NND模型对多种激波相互作用场景进行模拟，包括激波与激波、激波与边界的相互作用，通过对模拟结果的详细分析，揭示其中的流动特性和物理规律。在激波与激波相互作用的模拟中，以两个斜激波相互作用为例。图5展示了马赫数为3的两个斜激波相互作用过程中的压力云图。从图中可以清晰地观察到，当两个斜激波逐渐靠近并相互作用时，在相互作用区域形成了复杂的波系结构。在波系结构中，除了原有的两个斜激波外，还产生了反射激波和折射激波。反射激波和折射激波的强度和方向与原斜激波的特性以及相互作用的角度密切相关。在相互作用区域，压力分布呈现出复杂的变化，压力峰值明显高于周围区域。这是因为激波相互作用导致流体的压缩和能量的聚集，使得压力急剧升高。通过对压力云图的分析，可以准确地确定激波相互作用的位置和范围，以及波系结构的形成和演变过程。进一步分析激波相互作用区域的速度矢量分布。图6展示了同一时刻激波相互作用区域的速度矢量图。从图中可以看出，在激波相互作用区域，速度