多指示子驱动的杂交WENO格式：理论、算法与应用的深度剖析

多指示子驱动的杂交WENO格式：理论、算法与应用的深度剖析一、引言1.1研究背景与意义在现代科学与工程计算领域，数值模拟技术已成为不可或缺的研究手段，尤其在计算流体力学（CFD）中，对于各类复杂流动现象的精确模拟至关重要。随着航空航天、能源动力、水利工程等行业的快速发展，对CFD数值模拟的精度和效率提出了更高要求。在这些应用场景中，涉及到的流体流动往往包含激波、边界层、漩涡等复杂的流动结构，这些结构相互作用，使得流场呈现出高度的非线性和非定常特性。准确捕捉这些复杂流动结构，对于理解流体的物理机制、优化工程设计具有重要意义。例如，在航空发动机的设计中，精确模拟燃烧室中的燃烧过程和气流流动，可以提高燃烧效率，降低污染物排放；在水利工程中，准确预测河道中的水流速度和水位变化，有助于防洪减灾和水资源合理利用。加权本质无振荡（WENO）格式作为一类高精度数值方法，自提出以来在CFD及其他相关领域得到了广泛应用。WENO格式能够在解的光滑区域达到高阶精度，有效减少数值耗散和色散误差，从而更准确地分辨流场中的细微结构；在解的间断附近，它又能保持本质无振荡的特性，避免出现非物理的数值振荡，这对于准确捕捉激波等强间断具有重要作用。然而，在实际应用中，单一的WENO格式往往难以满足所有复杂流动问题的需求。不同的流动问题具有不同的特点，例如，在一些流动中，激波的强度和位置变化剧烈，而在另一些流动中，边界层的特性对计算结果影响较大。因此，为了进一步提高数值模拟的精度和效率，发展基于多种指示子的杂交WENO格式具有重要的现实意义。通过引入多种指示子，可以根据流场的局部特征，自适应地选择合适的WENO格式或其组合，从而在不同的流动区域发挥各种格式的优势，实现对复杂流动的更精确模拟。1.2WENO格式发展历程回顾WENO格式的发展是数值计算领域的一个重要历程，它的诞生源于对高精度、无振荡数值格式的追求。20世纪80年代，随着计算流体力学的兴起，人们对数值格式的精度和稳定性提出了更高要求。传统的低阶格式，如一阶迎风格式和二阶中心差分格式，在处理复杂流动时存在明显的局限性，如数值耗散大、无法准确捕捉间断等问题。为了解决这些问题，Harten等人在1987年提出了本质无振荡（ENO）格式，ENO格式通过在重构过程中选择最光滑的模板，能够在间断附近保持无振荡特性，并且在光滑区域达到高阶精度。然而，ENO格式在实际应用中也存在一些问题，例如在光滑区域，它对模板的选择较为保守，导致计算效率较低，而且逻辑判断较为复杂，不利于并行计算。为了改进ENO格式的不足，Liu、Osher和Chan在1994年首次提出了加权本质无振荡（WENO）格式的概念，并由Jiang和Shu在1996年进一步完善。WENO格式的基本思想是利用所有子模板的加权凸组合来构造最终的数值通量。在解的光滑区域，各个子模板的权值逼近某个理想权值，从而使格式获得最高阶的精度；在含间断的区域，那些含间断的子模板的权值接近于0，使得格式具有基本无振荡（ENO）的性质。这种设计思想使得WENO格式在光滑区域和间断附近都能表现出良好的性能，既提高了计算精度，又保证了数值稳定性。此后，WENO格式得到了广泛的研究和发展。众多学者从不同角度对WENO格式进行了改进和拓展。在精度提升方面，研究人员通过优化光滑因子、改进权值计算方法等手段，不断提高WENO格式的收敛阶数和数值精度。例如，Borges等人在2008年提出了WENO-Z格式，通过引入新的光滑度指标，使得格式在光滑解的情况下具有更好的性能。在格式的适应性和鲁棒性方面，学者们提出了各种改进措施，以增强WENO格式对不同类型流动问题的处理能力。例如，针对复杂几何形状的网格，发展了适用于曲线坐标系的WENO格式；针对多相流、湍流等复杂流动，对WENO格式进行了相应的改进和应用研究。随着计算机技术的发展，WENO格式也在不断适应并行计算和大规模计算的需求，通过优化算法结构和数据存储方式，提高了计算效率和可扩展性。如今，WENO格式已经成为计算流体力学、天体物理学、燃烧理论等众多领域中求解双曲守恒律方程和对流占优问题的重要数值方法之一，其应用范围不断拓展，研究也在持续深入。1.3杂交WENO格式研究现状随着对复杂流动现象研究的不断深入，单一的WENO格式在某些情况下难以满足高精度和高效率的要求，因此杂交WENO格式应运而生。杂交WENO格式通过结合多种不同的WENO格式或其他数值格式，充分发挥各格式的优势，以提高对复杂流动的模拟能力。目前，杂交WENO格式的研究已经取得了一定的成果。在格式的构建方面，研究者们提出了多种不同的杂交策略。例如，将不同阶数的WENO格式进行杂交，在光滑区域使用高阶WENO格式以获得高精度，在间断附近切换到低阶WENO格式以保证稳定性和无振荡性。还有将WENO格式与其他类型的格式，如中心差分格式、TVD（TotalVariationDiminishing）格式等进行杂交，利用不同格式在不同流动特性下的优势，实现对复杂流场的更精确模拟。在指示子的应用方面，为了实现格式的自适应切换，多种指示子被引入到杂交WENO格式中。常见的指示子包括基于流场物理量梯度的指示子、基于局部光滑度的指示子以及基于激波检测的指示子等。基于流场物理量梯度的指示子通过计算流场中物理量（如速度、压力等）的梯度来判断流场的变化剧烈程度，当梯度较大时，认为可能存在间断或复杂流动结构，从而切换到更适合处理这种情况的格式；基于局部光滑度的指示子则通过分析局部区域的光滑性来确定格式的选择，在光滑区域采用高精度格式，在非光滑区域采用更稳健的格式；基于激波检测的指示子专门用于检测激波的位置和强度，当检测到激波时，采用能够有效捕捉激波的格式。然而，现有杂交WENO格式的研究仍存在一些不足之处。一方面，不同指示子的性能和适用范围存在差异，如何合理选择和组合指示子，以实现对各种复杂流动的有效判断和格式切换，仍然是一个有待深入研究的问题。另一方面，在杂交格式的稳定性和收敛性分析方面，目前的研究还不够完善，缺乏统一的理论框架来保证杂交格式在各种情况下的可靠性。此外，随着计算需求的不断增加，如何进一步提高杂交WENO格式的计算效率，降低计算成本，也是需要解决的重要问题。因此，对基于多种指示子的杂交WENO格式的研究具有重要的必要性和紧迫性，有望为复杂流动的数值模拟提供更有效的方法。1.4本文研究内容与创新点本文主要致力于基于多种指示子的杂交WENO格式的研究，旨在构建一种能够更高效、精确地模拟复杂流动的数值方法。具体研究内容包括：首先，系统地研究和分析现有的各种指示子，包括它们的计算方法、性能特点以及适用范围，在此基础上，提出一种新的多指示子融合策略。通过将不同类型的指示子进行有机结合，充分利用它们各自的优势，实现对流场更全面、准确的判断。例如，将基于物理量梯度的指示子与基于局部光滑度的指示子相结合，既能捕捉流场中物理量的剧烈变化，又能准确判断局部区域的光滑程度，从而更合理地选择数值格式。其次，根据提出的多指示子融合策略，构建全新的杂交WENO格式。在格式构建过程中，详细分析不同WENO格式在不同流动条件下的性能表现，选择合适的WENO格式进行杂交。通过精心设计格式的切换机制，确保在不同流场区域能够无缝地切换到最适合的格式，从而提高整个格式的精度和稳定性。同时，对构建的杂交WENO格式进行严格的理论分析，包括精度分析、稳定性分析和收敛性分析，从理论上保证格式的可靠性。在数值实验方面，选取一系列具有代表性的一维和多维算例，包括Euler方程、Burgers方程、浅水波方程以及带有源项的浅水中污染物输运方程等，对所提出的杂交WENO格式进行全面的验证和测试。通过与其他经典的数值格式进行对比，详细分析杂交WENO格式在精度、稳定性、计算效率等方面的优势和不足。针对数值实验中发现的问题，进一步优化和改进杂交WENO格式，以提高其实际应用能力。本文的创新点主要体现在以下两个方面：一是在指示子的选取和融合策略上，提出了一种全新的多指示子融合方法，该方法能够更全面、准确地反映流场的局部特征，为格式的自适应切换提供更可靠的依据。二是在杂交WENO格式的构建上，通过精心设计格式的切换机制和优化格式组合，使得构建的杂交WENO格式在复杂流动模拟中具有更高的精度、更好的稳定性和更强的适应性，能够有效解决现有杂交WENO格式存在的一些问题，为复杂流动的数值模拟提供了一种新的有效手段。二、WENO格式理论基础2.1双曲守恒律方程与WENO格式双曲守恒律方程是一类在科学与工程领域中广泛应用的偏微分方程，其基本形式在一维情况下可表示为：\frac{\partial\boldsymbol{u}}{\partialt}+\frac{\partial\boldsymbol{f}(\boldsymbol{u})}{\partialx}=0其中，\boldsymbol{u}为守恒变量向量，例如在流体力学中，对于理想气体的欧拉方程，\boldsymbol{u}=\begin{pmatrix}\rho\\\rhou\\E\end{pmatrix}，分别表示密度、动量和总能量；\boldsymbol{f}(\boldsymbol{u})是通量函数向量，与守恒变量相关，对应上述欧拉方程的通量函数为\boldsymbol{f}(\boldsymbol{u})=\begin{pmatrix}\rhou\\\rhou^2+p\\u(E+p)\end{pmatrix}，这里p为压力，u为速度。双曲守恒律方程描述了物理量在空间和时间上的守恒特性，广泛应用于流体力学、电磁学、交通流等领域，如在流体力学中用于描述流体的流动，在电磁学中用于描述电磁波的传播。WENO格式作为求解双曲守恒律方程的一种高精度数值方法，其基本原理是基于非线性加权的思想。在重构数值通量时，WENO格式考虑多个不同的低阶子模板。以五阶WENO格式为例，通常会选取三个三阶子模板。对于每个子模板，通过线性重构方法得到低阶数值通量。然后，引入光滑度指标来衡量每个子模板的光滑程度。光滑度指标的计算基于子模板上函数值的导数信息，例如对于五点模板的三阶子模板，其光滑度指标\beta_{i}（i=0,1,2分别对应不同子模板）的计算式为：\beta_{0}=(\overline{u}_{j-2}-2\overline{u}_{j-1}+\overline{u}_{j})^2+\frac{13}{12}(\overline{u}_{j-2}-4\overline{u}_{j-1}+3\overline{u}_{j})^2\beta_{1}=(\overline{u}_{j-1}-2\overline{u}_{j}+\overline{u}_{j+1})^2+\frac{13}{12}(\overline{u}_{j-1}-\overline{u}_{j+1})^2\beta_{2}=(\overline{u}_{j}-2\overline{u}_{j+1}+\overline{u}_{j+2})^2+\frac{13}{12}(3\overline{u}_{j}-4\overline{u}_{j+1}+\overline{u}_{j+2})^2其中\overline{u}_{k}表示网格点k上的函数值。根据光滑度指标，计算出每个子模板的非线性权重\omega_{i}：\omega_{i}=\frac{\alpha_{i}}{\sum_{k=0}^{2}\alpha_{k}}\alpha_{i}=\frac{d_{i}}{(\epsilon+\beta_{i})^2}这里d_{i}是理想权重，\epsilon是一个小常数，用于避免分母为零的情况，通常取10^{-6}到10^{-10}之间的值。最后，将这些非线性权重与低阶数值通量相结合，得到最终的高阶数值通量。通过这种方式，在解的光滑区域，光滑度指标较小的子模板对应的权重会趋近于理想权重，使得格式能够达到高阶精度；在含间断的区域，含有间断的子模板的光滑度指标较大，其权重会趋近于0，从而保证格式在间断附近的本质无振荡特性。与传统的数值格式相比，如一阶迎风格式和二阶中心差分格式，WENO格式在精度和间断捕捉能力上具有明显优势。一阶迎风格式虽然能够稳定地计算，但数值耗散较大，导致解的分辨率较低，无法准确捕捉流场中的细微结构；二阶中心差分格式在光滑区域具有较好的精度，但在间断附近会产生数值振荡，影响解的准确性。而WENO格式通过其独特的非线性加权机制，既能够在光滑区域保持高阶精度，又能在间断附近有效抑制数值振荡，准确捕捉激波等强间断，因此在复杂流动的数值模拟中得到了广泛应用。2.2一维WENO格式详解在一维WENO格式中，模板的选择是构建格式的基础。以五阶WENO格式为例，考虑在均匀网格x_{j}=j\Deltax（j=\cdots,-1,0,1,\cdots）上求解双曲守恒律方程。通常选取三个三阶子模板，分别为：S_{0}=\{x_{j-2},x_{j-1},x_{j}\}S_{1}=\{x_{j-1},x_{j},x_{j+1}\}S_{2}=\{x_{j},x_{j+1},x_{j+2}\}这些子模板覆盖了不同的网格点组合，为后续的权重计算和数值通量重构提供了基础。权重计算是一维WENO格式的关键步骤。首先，对于每个子模板，需要计算其低阶数值通量。以子模板S_{0}为例，使用线性重构方法得到低阶数值通量f_{0,j+\frac{1}{2}}，假设守恒变量在子模板上的分布可以用线性多项式来近似，通过对该线性多项式进行积分等运算，可以得到f_{0,j+\frac{1}{2}}的表达式。类似地，可以得到子模板S_{1}和S_{2}的低阶数值通量f_{1,j+\frac{1}{2}}和f_{2,j+\frac{1}{2}}。然后，计算每个子模板的光滑度指标\beta_{i}（i=0,1,2），如前文所述，光滑度指标是基于子模板上函数值的导数信息来衡量子模板的光滑程度。以\beta_{0}为例，其计算式为：\beta_{0}=(\overline{u}_{j-2}-2\overline{u}_{j-1}+\overline{u}_{j})^2+\frac{13}{12}(\overline{u}_{j-2}-4\overline{u}_{j-1}+3\overline{u}_{j})^2其中\overline{u}_{k}表示网格点k上的守恒变量值。根据光滑度指标，计算非线性权重\omega_{i}。首先定义\alpha_{i}：\alpha_{i}=\frac{d_{i}}{(\epsilon+\beta_{i})^2}其中d_{i}是理想权重，对于五阶WENO格式，通常取d_{0}=\frac{1}{10}，d_{1}=\frac{6}{10}，d_{2}=\frac{3}{10}；\epsilon是一个小常数，用于避免分母为零的情况，一般取值在10^{-6}到10^{-10}之间。然后，非线性权重\omega_{i}为：\omega_{i}=\frac{\alpha_{i}}{\sum_{k=0}^{2}\alpha_{k}}最终的数值通量f_{j+\frac{1}{2}}通过将非线性权重与低阶数值通量进行加权组合得到：f_{j+\frac{1}{2}}=\omega_{0}f_{0,j+\frac{1}{2}}+\omega_{1}f_{1,j+\frac{1}{2}}+\omega_{2}f_{2,j+\frac{1}{2}}在计算过程中，需要注意边界条件的处理。对于周期边界条件，在边界处可以将边界点视为与另一侧的对应点相连，从而保证计算的连续性；对于非周期边界条件，如Dirichlet边界条件（给定边界上的函数值）或Neumann边界条件（给定边界上函数的导数），需要根据具体的边界条件对边界点的数值通量计算进行特殊处理。例如，在Dirichlet边界条件下，直接使用给定的边界函数值参与数值通量的计算；在Neumann边界条件下，通过对边界点的导数信息进行近似处理，以保证格式在边界处的精度和稳定性。2.3多维WENO格式拓展从一维WENO格式拓展到多维情况时，通常采用张量积的方法。以二维为例，考虑双曲守恒律方程：\frac{\partial\boldsymbol{u}}{\partialt}+\frac{\partial\boldsymbol{f}(\boldsymbol{u})}{\partialx}+\frac{\partial\boldsymbol{g}(\boldsymbol{u})}{\partialy}=0其中\boldsymbol{u}是守恒变量向量，\boldsymbol{f}(\boldsymbol{u})和\boldsymbol{g}(\boldsymbol{u})分别是x和y方向的通量函数向量。在构建多维WENO格式时，可以分别在x和y方向上应用一维WENO格式的思想。对于x方向，类似于一维情况，选择合适的模板，计算光滑度指标和权重，得到x方向的数值通量f_{i+\frac{1}{2},j}；对于y方向，同样选择模板，计算相关参数，得到y方向的数值通量g_{i,j+\frac{1}{2}}。然后，通过某种组合方式将两个方向的数值通量结合起来，得到二维情况下的数值通量。一种常见的组合方式是采用交替方向法，即在每个时间步中，先在x方向上更新守恒变量，再在y方向上进行更新。然而，这种拓展面临着一些挑战。首先是计算量的大幅增加，在多维情况下，模板的数量和计算复杂度都显著提高。以二维五阶WENO格式为例，每个方向上选择三个子模板，那么在二维空间中，组合后的模板数量会大幅增加，相应的光滑度指标和权重计算量也会急剧上升。其次，在处理复杂边界条件和非结构网格时，难度明显加大。对于复杂边界条件，如不规则边界的几何形状，难以像一维情况那样简单地处理边界点；在非结构网格中，由于网格的不规则性，模板的定义和权重计算变得更加复杂，需要特殊的处理方法来保证格式的精度和稳定性。为了解决这些问题，研究者们提出了一系列策略。针对计算量增加的问题，采用并行计算技术，将计算任务分配到多个处理器上，提高计算效率。例如，使用MPI（MessagePassingInterface）等并行编程模型，实现计算任务的并行化。在处理复杂边界条件和非结构网格时，发展了适用于非结构网格的WENO格式，如基于三角形或四面体网格的WENO格式。对于复杂边界条件，可以采用边界拟合网格技术，使网格与边界形状相适应，或者使用浸入边界法等方法，将边界条件的处理与内部网格的计算分开，从而降低处理难度，保证格式在多维复杂情况下的有效应用。2.4带有源项方程的WENO格式应用在实际应用中，许多物理问题的控制方程不仅包含对流项，还带有源项，如浅水波方程和浅水中污染物输运方程。以浅水波方程为例，其在一维情况下可表示为：\frac{\partialh}{\partialt}+\frac{\partial(hu)}{\partialx}=0\frac{\partial(hu)}{\partialt}+\frac{\partial(hu^{2}+\frac{1}{2}gh^{2})}{\partialx}=-gh\frac{\partialz_{b}}{\partialx}其中h是水深，u是流速，g是重力加速度，z_{b}是河床地形。这里-gh\frac{\partialz_{b}}{\partialx}就是源项，它反映了河床地形对水流的影响。对于带有源项的方程，在应用WENO格式时，需要对源项进行适当的处理。一种常见的方法是将源项视为一个独立的部分，与对流项分别进行计算。在时间离散上，可以采用分步算法，先使用WENO格式计算对流项对守恒变量的影响，得到一个中间值；然后，根据源项的具体形式，采用合适的数值方法计算源项对中间值的修正，得到最终的守恒变量值。例如，对于浅水波方程中的源项，可以采用有限差分法对\frac{\partialz_{b}}{\partialx}进行离散，然后将离散后的源项与对流项的计算结果相结合。再以浅水中污染物输运方程为例，其一般形式为：\frac{\partialc}{\partialt}+\frac{\partial(uc)}{\partialx}=D\frac{\partial^{2}c}{\partialx^{2}}+S其中c是污染物浓度，u是水流速度，D是扩散系数，S是源项，如污染物的排放或降解。在这种情况下，除了对流项\frac{\partial(uc)}{\partialx}和源项S外，还包含扩散项D\frac{\partial^{2}c}{\partialx^{2}}。在应用WENO格式时，对于对流项可以采用WENO格式进行离散，对于扩散项可以采用中心差分等方法进行离散，而源项则根据其具体物理意义进行相应的数值处理。通过这种方式，能够有效地将WENO格式应用于带有源项的复杂方程，实现对相关物理过程的准确模拟，为实际工程和科学研究提供有力的数值工具。2.5时间离散方法适配在数值求解双曲守恒律方程时，时间离散方法的选择与WENO格式的适配性对计算效率和精度有着重要影响。常见的时间离散方法包括显式方法和隐式方法。显式方法如Runge-Kutta方法，以三阶TVD（TotalVariationDiminishing）Runge-Kutta方法为例，其基本步骤如下：\boldsymbol{u}^{(1)}=\boldsymbol{u}^{n}+\DeltatL(\boldsymbol{u}^{n})\boldsymbol{u}^{(2)}=\frac{3}{4}\boldsymbol{u}^{n}+\frac{1}{4}\boldsymbol{u}^{(1)}+\frac{1}{4}\DeltatL(\boldsymbol{u}^{(1)})\boldsymbol{u}^{n+1}=\frac{1}{3}\boldsymbol{u}^{n}+\frac{2}{3}\boldsymbol{u}^{(2)}+\frac{2}{3}\DeltatL(\boldsymbol{u}^{(2)})其中\boldsymbol{u}^{n}表示n时刻的守恒变量，\Deltat是时间步长，L(\boldsymbol{u})是空间离散后的算子，包含了WENO格式计算得到的数值通量。显式方法的优点是计算简单，易于编程实现，并且在解的光滑区域能够保持较高的精度。然而，其稳定性条件较为苛刻，时间步长受到Courant-Friedrichs-Lewy（CFL）条件的限制，即\Deltat\leqC\frac{\Deltax}{v_{max}}，其中C是CFL数，一般取小于1的值，\Deltax是空间步长，v_{max}是流场中的最大波速。这意味着在一些波速较大或空间步长较小的情况下，显式方法需要采用非常小的时间步长，从而导致计算量大幅增加。隐式方法如向后欧拉法，其时间离散公式为：\frac{\boldsymbol{u}^{n+1}-\boldsymbol{u}^{n}}{\Deltat}+L(\boldsymbol{u}^{n+1})=0隐式方法的优势在于稳定性好，时间步长不受CFL条件的严格限制，可以采用较大的时间步长进行计算，从而提高计算效率。但是，隐式方法需要求解一个关于\boldsymbol{u}^{n+1}的非线性方程组，计算复杂度较高，通常需要使用迭代方法如牛顿迭代法等来求解，这增加了计算的复杂性和计算时间。在选择时间离散方法与WENO格式适配时，需要综合考虑计算效率和精度。对于一些对流占优且流场变化较为平缓的问题，可以选择显式的Runge-Kutta方法，通过合理调整CFL数和空间步长，在保证精度的前提下提高计算效率。例如，在模拟一些低速流体流动问题时，显式方法能够较好地满足计算需求。而对于一些含有强间断或波速变化较大的复杂流场问题，为了保证稳定性和计算效率，隐式方法可能更为合适。例如，在模拟激波等强间断问题时，隐式方法能够在较大时间步长下保持数值解的稳定性，虽然求解非线性方程组会增加一定的计算量，但相比显式方法在小时间步长下的大量计算，总体计算效率可能更高。此外，还可以考虑采用混合时间离散方法，如在解的光滑区域使用显式方法，在间断附近或波速变化三、多种指示子分析3.1常见指示子类型及原理在基于多种指示子的杂交WENO格式研究中，常见的指示子类型主要包括坏单元指示子和光滑度指示子，它们在数值模拟复杂流动时发挥着关键作用。坏单元指示子旨在标记间断解附近的单元，这些单元在数值计算中容易产生振荡，影响计算结果的准确性和稳定性。以KXRCF（Koren,vanLeer,Roe,andChakravarthyFluxLimiter）坏单元指示子为例，其工作原理基于对数值解的梯度分析。通过计算单元内数值解的梯度，判断该单元内解的变化剧烈程度。在间断附近，数值解的梯度会急剧增大，KXRCF指示子能够捕捉到这种变化，将对应的单元标记为坏单元。具体来说，KXRCF指示子通过比较相邻单元间数值解的变化率，结合一些经验参数来确定单元是否为坏单元。在激波等强间断问题中，KXRCF指示子可以有效地标记出激波附近的坏单元，为后续的格式调整提供依据。然而，KXRCF指示子存在依赖于问题本身参数的局限性，在不同的流动问题中，需要根据数值试验对这些参数进行调整，增加了使用的复杂性。另一种常见的坏单元指示子是TVB（TotalVariationBounded）指示子。TVB指示子的原理基于总变差有界的概念，通过控制数值解的总变差来判断单元是否为坏单元。在数值计算中，总变差反映了数值解在空间上的变化总量。如果某个单元内数值解的总变差超过了一定的阈值，TVB指示子就会将其标记为坏单元。在模拟含有复杂间断结构的流动时，TVB指示子能够根据总变差的变化，准确地识别出坏单元，从而保证数值解在间断附近的稳定性。但是，TVB指示子也面临着参数调整的问题，不同的流动场景需要不同的阈值参数，这在一定程度上限制了其通用性。光滑度指示子则侧重于衡量流场局部区域的光滑程度，为在光滑区域选择高精度格式、在非光滑区域选择稳健格式提供依据。以基于多项式重构的光滑度指示子为例，其计算原理是通过对局部区域内的数值解进行多项式重构，然后根据重构多项式的导数信息来计算光滑度指标。假设在一个局部模板上，对数值解进行三次多项式重构，通过计算该多项式的一阶导数和二阶导数的平方和等方式来得到光滑度指标。在光滑区域，数值解的变化较为平缓，重构多项式的导数较小，光滑度指标也较小；而在非光滑区域，数值解变化剧烈，导数较大，光滑度指标相应增大。在模拟边界层流动时，边界层内的速度、温度等物理量变化剧烈，基于多项式重构的光滑度指示子能够准确地检测到边界层区域的非光滑性，从而为选择合适的数值格式提供准确的判断。还有一种基于小波变换的光滑度指示子。其工作原理是利用小波变换能够将信号分解为不同频率成分的特性，对数值解进行小波变换，分析不同频率成分的能量分布来判断光滑度。在光滑区域，低频成分的能量占主导；而在非光滑区域，高频成分的能量会显著增加。通过计算高频成分与低频成分的能量比值等方式，可以得到基于小波变换的光滑度指示子的值。在处理含有复杂漩涡结构的流动时，基于小波变换的光滑度指示子能够有效地检测到漩涡区域的非光滑特性，因为漩涡内部的流动具有强烈的非线性和非均匀性，会产生丰富的高频成分，该指示子能够敏锐地捕捉到这些变化，为格式的选择提供重要参考。3.2指示子性能对比与分析为了深入了解不同指示子的性能，通过一系列数值实验对其在捕捉间断、判断光滑区域等方面的能力进行对比分析。在捕捉间断方面，选取经典的激波管算例，如Lax激波管问题。对于坏单元指示子，KXRCF指示子在该算例中能够较好地标记出激波所在的单元，使得在这些单元上可以采用更适合捕捉间断的格式，从而有效抑制激波附近的数值振荡。然而，由于其参数依赖于具体问题，在不同强度激波的情况下，需要对参数进行调整，否则可能会出现过度标记或标记不足的情况。例如，当激波强度较弱时，若参数设置不当，KXRCF指示子可能会将一些正常单元误标记为坏单元，增加不必要的计算成本；而当激波强度较强时，可能会因为参数不合适而未能准确标记出所有含间断的单元，导致激波附近的数值振荡无法得到有效控制。TVB指示子在Lax激波管算例中，能够根据总变差有界的原则，较为准确地识别出激波附近的坏单元。与KXRCF指示子相比，TVB指示子在参数调整相对简单的情况下，也能在一定程度上保证对间断的有效捕捉。但是，由于其对总变差的控制较为严格，在一些复杂流场中，可能会将一些虽然总变差稍大但并非真正间断的单元标记为坏单元，从而影响计算效率。在存在多个弱间断相互作用的流场中，TVB指示子可能会因为对总变差的敏感，将一些弱间断附近的正常单元误判为坏单元，使得格式切换过于频繁，增加计算量。在判断光滑区域方面，采用一个简单的一维光滑流场算例进行测试。基于多项式重构的光滑度指示子在该算例中能够准确地判断出光滑区域，其计算得到的光滑度指标在光滑区域数值很小，表明该区域的解是光滑的，适合采用高阶WENO格式以获得高精度的数值解。然而，当流场中存在局部微小的扰动时，由于多项式重构对这些小扰动的敏感性，可能会导致光滑度指标出现一定的波动，从而影响对光滑区域的准确判断。在流场中存在小尺度的噪声时，基于多项式重构的光滑度指示子可能会将这些噪声区域误判为非光滑区域，导致在这些区域采用不必要的低阶稳健格式，降低计算精度。基于小波变换的光滑度指示子在光滑流场算例中，能够通过分析小波变换后的频率成分，清晰地分辨出光滑区域和非光滑区域。在光滑区域，低频成分占据主导，其光滑度指示子的值较小；而在非光滑区域，高频成分的能量增加，指示子的值增大。该指示子对光滑区域的判断相对稳定，受局部小扰动的影响较小。但是，小波变换的计算过程相对复杂，需要较高的计算成本，这在一定程度上限制了其在大规模计算中的应用。在三维复杂流场的数值模拟中，由于需要对大量的数值解数据进行小波变换，基于小波变换的光滑度指示子的计算时间会显著增加，可能会影响整个计算的效率。3.3新型指示子的提出与验证基于现有指示子在实际应用中存在的不足，提出一种新型指示子，旨在更全面、准确地反映流场的局部特征，提高杂交WENO格式的性能。新型指示子融合了物理量梯度信息和局部光滑度信息，通过构建一个综合的判断函数来实现对流场的有效判断。具体来说，新型指示子首先计算流场中物理量（如速度、压力等）的梯度，以捕捉流场中物理量的剧烈变化，判断是否存在间断或强梯度区域。然后，结合基于多项式重构的局部光滑度信息，分析局部区域的光滑程度。通过将这两种信息进行合理的加权组合，得到一个综合的指示子值。在判断激波时，不仅考虑激波处物理量的突变（通过梯度反映），还考虑激波附近区域的光滑性（通过局部光滑度反映），从而更准确地识别激波的位置和范围。从理论分析角度来看，新型指示子具有以下优势。在光滑区域，由于物理量梯度较小且局部光滑度高，综合指示子值较小，能够准确地判断该区域为光滑区域，从而可以选择高阶WENO格式进行计算，保证高精度。在间断附近，物理量梯度急剧增大，同时局部光滑度降低，综合指示子值会显著增大，能够有效地识别出间断区域，进而切换到更适合处理间断的格式，保证数值解的稳定性。与传统指示子相比，新型指示子避免了单一指示子的局限性，能够更全面地反映流场的复杂特性。为了验证新型指示子的有效性和优越性，进行了一系列数值实验。选取了多个具有代表性的算例，包括一维和二维的激波问题、边界层问题以及含有复杂漩涡结构的流场问题。在一维激波管算例中，新型指示子能够准确地标记出激波所在的单元，并且与KXRCF和TVB等传统坏单元指示子相比，在参数调整方面更加简单和鲁棒。新型指示子能够在不同强度激波的情况下，无需频繁调整参数，就可以准确地识别出激波附近的坏单元，有效抑制激波附近的数值振荡，提高了计算结果的准确性。在二维边界层算例中，新型指示子能够精确地判断边界层区域的非光滑性，与基于多项式重构的光滑度指示子和基于小波变换的光滑度指示子相比，新型指示子在判断边界层区域时更加准确和稳定。在边界层内，物理量的梯度变化较大，同时局部光滑度较低，新型指示子能够综合考虑这两个因素，准确地识别出边界层区域，从而为在边界层区域选择合适的数值格式提供可靠依据。在边界层厚度变化的情况下，新型指示子依然能够准确地捕捉到边界层的位置和范围，使得在边界层区域采用的格式能够更好地适应边界层的特性，提高了边界层模拟的精度。在含有复杂漩涡结构的流场算例中，新型指示子同样表现出良好的性能。它能够清晰地分辨出漩涡区域的非光滑特性，并且能够准确地判断漩涡与周围光滑区域的边界。与其他指示子相比，新型指示子在处理复杂漩涡结构时，能够更全面地反映漩涡内部和周围流场的特征，从而使得在该区域的数值模拟更加准确。在漩涡强度和位置变化的情况下，新型指示子能够及时捕捉到这些变化，为格式的自适应切换提供准确的信号，保证了对复杂漩涡结构的有效模拟。通过这些数值实验，充分验证了新型指示子在提高杂交WENO格式对复杂流场模拟能力方面的有效性和优越性。四、基于多种指示子的杂交WENO格式构建4.1杂交WENO格式基本思想杂交WENO格式的核心在于巧妙地融合线性方法和WENO格式，以实现对复杂流场的高效模拟。线性方法，如中心差分格式和迎风格式，具有计算简单、效率高的显著优势。在流场的光滑区域，物理量变化较为平缓，线性方法能够快速且准确地给出近似解，有效减少计算量。例如，在模拟均匀来流的简单流场时，中心差分格式可以高效地计算出物理量的空间导数，得到较为精确的数值解。然而，当流场中出现间断，如激波时，线性方法由于其固有的数值耗散和色散特性，容易在间断附近产生非物理的数值振荡，导致解的失真。WENO格式则在间断捕捉方面表现出色。它通过非线性加权的方式，结合多个低阶子模板来重构数值通量。在解的光滑区域，WENO格式能够自动调整权重，使得各个子模板的贡献趋近于理想情况，从而达到高阶精度，有效减少数值耗散和色散误差，精确分辨流场中的细微结构。在模拟边界层流动时，WENO格式可以准确捕捉边界层内速度、温度等物理量的急剧变化，清晰地分辨出边界层的厚度和特性。在间断附近，WENO格式能够通过对光滑度指标的计算和权重的调整，使含有间断的子模板的权值趋近于0，从而保持本质无振荡的特性，准确捕捉激波等强间断。为了充分发挥两者的优势，杂交WENO格式利用指示子对流场区域进行分类。指示子作为判断流场特性的关键工具，能够根据流场中物理量的变化特征，如梯度、光滑度等，准确地识别出光滑区域和间断区域。基于物理量梯度的指示子通过计算流场中速度、压力等物理量的梯度，当梯度超过一定阈值时，判断该区域可能存在间断；基于局部光滑度的指示子则通过分析局部区域内物理量的变化是否平缓，来判断该区域的光滑程度。在模拟激波与边界层相互作用的复杂流场时，指示子能够准确地标记出激波所在的区域以及边界层的位置。在光滑区域，选择计算简单、效率高的线性方法进行数值计算，以提高计算效率；在间断区域，切换到能够有效捕捉间断的WENO格式，确保数值解的准确性和稳定性。通过这种方式，杂交WENO格式实现了对复杂流场的精确模拟，既保证了计算效率，又提高了数值解的质量。4.2基于不同指示子的格式构建过程基于多种指示子构建杂交WENO格式是一个系统而严谨的过程，主要包括网格单元分类、格式选择和数值通量计算等关键步骤。在网格单元分类阶段，以坏单元指示子KXRCF为例，对于每个网格单元，通过计算单元内数值解的梯度来判断其是否为坏单元。假设在一个一维网格中，网格单元I_j=[x_{j-\frac{1}{2}},x_{j+\frac{1}{2}}]，计算该单元内数值解u在x_j处的梯度\nablau_j，如果|\nablau_j|超过了预先设定的阈值\delta，则将该单元标记为坏单元，即认为该单元内可能存在间断。对于基于局部光滑度的指示子，如基于多项式重构的光滑度指示子，在一个局部模板\{x_{j-2},x_{j-1},x_{j},x_{j+1},x_{j+2}\}上，对数值解u进行三次多项式重构P(x)，然后计算P(x)的一阶导数P^\prime(x)和二阶导数P^{\prime\prime}(x)，通过公式\beta=\int_{x_{j-\frac{1}{2}}}^{x_{j+\frac{1}{2}}}(P^\prime(x))^2+(P^{\prime\prime}(x))^2dx计算光滑度指标\beta，若\beta超过一定阈值\epsilon，则认为该单元为非光滑单元。通过这些指示子的计算，将计算区域内的网格单元准确地分为光滑单元和间断单元，为后续的格式选择提供依据。根据网格单元的分类结果进行格式选择。对于被标记为光滑单元的网格，选择线性方法进行数值计算。在光滑区域，采用二阶中心差分格式计算数值通量。对于x_{j+\frac{1}{2}}处的数值通量f_{j+\frac{1}{2}}，其计算公式为f_{j+\frac{1}{2}}=\frac{1}{2}(f(u_{j+1})+f(u_{j}))，其中f(u)是通量函数，u_j是网格点j上的守恒变量值。这种格式计算简单，在光滑区域能够快速得到较为准确的数值解，提高计算效率。对于间断单元，选择合适的WENO格式。以五阶WENO格式为例，在间断单元上，首先确定三个三阶子模板S_0=\{x_{j-2},x_{j-1},x_{j}\}，S_1=\{x_{j-1},x_{j},x_{j+1}\}，S_2=\{x_{j},x_{j+1},x_{j+2}\}，然后分别计算每个子模板的低阶数值通量f_{i,j+\frac{1}{2}}（i=0,1,2），再根据前文所述的方法计算光滑度指标\beta_i和非线性权重\omega_i，最终得到该间断单元在x_{j+\frac{1}{2}}处的数值通量f_{j+\frac{1}{2}}=\omega_{0}f_{0,j+\frac{1}{2}}+\omega_{1}f_{1,j+\frac{1}{2}}+\omega_{2}f_{2,j+\frac{1}{2}}，通过这种方式，在间断区域能够有效地捕捉间断，保证数值解的稳定性和准确性。在数值通量计算方面，无论是线性方法还是WENO格式计算得到的数值通量，都需要与时间离散方法相结合，以推进数值解在时间上的演化。在时间离散上采用三阶TVDRunge-Kutta方法，对于守恒变量\boldsymbol{u}，其时间推进公式为：\boldsymbol{u}^{(1)}=\boldsymbol{u}^{n}+\DeltatL(\boldsymbol{u}^{n})\boldsymbol{u}^{(2)}=\frac{3}{4}\boldsymbol{u}^{n}+\frac{1}{4}\boldsymbol{u}^{(1)}+\frac{1}{4}\DeltatL(\boldsymbol{u}^{(1)})\boldsymbol{u}^{n+1}=\frac{1}{3}\boldsymbol{u}^{n}+\frac{2}{3}\boldsymbol{u}^{(2)}+\frac{2}{3}\DeltatL(\boldsymbol{u}^{(2)})其中\boldsymbol{u}^{n}表示n时刻的守恒变量，\Deltat是时间步长，L(\boldsymbol{u})是空间离散后的算子，包含了根据网格单元类型选择的线性方法或WENO格式计算得到的数值通量。通过这种时间离散方法，将不同格式计算得到的数值通量融入到守恒变量的时间推进过程中，实现对复杂流场的动态模拟。4.3格式的稳定性与收敛性分析从理论分析角度来看，对于杂交WENO格式的稳定性，基于能量方法进行分析。考虑一个离散的守恒律方程的半离散形式\frac{d\boldsymbol{u}_j}{dt}=-\frac{1}{\Deltax}(F_{j+\frac{1}{2}}-F_{j-\frac{1}{2}})，其中\boldsymbol{u}_j是网格点j上的守恒变量，F_{j+\frac{1}{2}}是j+\frac{1}{2}处的数值通量，\Deltax是空间步长。定义离散能量E^n=\frac{1}{2}\sum_{j}\Deltax|\boldsymbol{u}_j^n|^2，对其求时间导数\frac{dE^n}{dt}=\sum_{j}\Deltax\boldsymbol{u}_j^n\frac{d\boldsymbol{u}_j^n}{dt}。将半离散形式代入可得\frac{dE^n}{dt}=-\sum_{j}\boldsymbol{u}_j^n(F_{j+\frac{1}{2}}-F_{j-\frac{1}{2}})。对于杂交WENO格式，在光滑区域使用线性方法，在间断区域使用WENO格式，通过分析数值通量F_{j+\frac{1}{2}}的性质，证明在一定的CFL（Courant-Friedrichs-Lewy）条件下，\frac{dE^n}{dt}\leq0，从而保证格式的稳定性。在满足CFL数C\leq1时，即\Deltat\leqC\frac{\Deltax}{v_{max}}，其中v_{max}是流场中的最大波速，杂交WENO格式的离散能量不会随时间增加，保证了数值解的稳定性。在收敛性分析方面，基于Lax-Richtmyer等价定理，该定理表明在一定条件下，一个适定的线性偏微分方程的数值格式，如果它是相容的且稳定的，那么它就是收敛的。对于杂交WENO格式，首先证明其相容性。相容性要求数值格式在网格尺寸趋于0时，能够逼近原偏微分方程。通过对光滑区域和间断区域分别进行Taylor展开分析，验证杂交WENO格式的数值通量在网格尺寸\Deltax\rightarrow0和时间步长\Deltat\rightarrow0时，能够准确地逼近原双曲守恒律方程的通量函数，从而满足相容性条件。在光滑区域，线性方法的数值通量与原方程通量的误差在O(\Deltax^p)（p为线性方法的精度阶数）量级；在间断区域，WENO格式的数值通量与原方程通量的误差在O(\Deltax^q)（q为WENO格式在间断区域保持稳定时的精度阶数）量级。结合前面证明的稳定性，根据Lax-Richtmyer等价定理，可以得出杂交WENO格式是收敛的。为了进一步验证格式的稳定性和收敛性，进行数值实验。在一维Euler方程的激波管算例中，采用不同的CFL数进行计算。当CFL数取0.5时，计算结果稳定，激波的位置和强度能够准确捕捉，数值解没有出现明显的振荡；当CFL数逐渐增大接近1时，仍然能够保持稳定，但当CFL数超过1时，数值解开始出现振荡，表明格式的稳定性受到破坏，验证了理论分析中CFL条件对稳定性的限制。在收敛性验证方面，通过加密网格，计算不同网格尺寸下的数值解与精确解的误差。在光滑区域，随着网格尺寸的减小，线性方法计算得到的数值解与精确解的误差按照理论分析的收敛阶数逐渐减小；在间断区域，WENO格式计算得到的数值解与精确解的误差也按照其在间断区域的收敛特性逐渐减小，从而验证了杂交WENO格式的收敛性。通过理论分析和数值实验，全面地验证了基于多种指示子的杂交WENO格式的稳定性和收敛性，为其在实际工程和科学研究中的应用提供了坚实的理论基础和实践依据。五、数值实验与结果分析5.1一维问题模拟5.1.1Euler方程算例为了验证基于多种指示子的杂交WENO格式在捕捉激波等方面的优势，选取经典的一维Euler方程算例——Lax激波管问题。该算例初始条件设定为：在x<0.5区域，\rho=0.445，u=0.698，p=3.528；在x\geq0.5区域，\rho=0.5，u=0，p=0.571。计算区域为[0,1]，采用均匀网格离散，网格数量分别取N=100、N=200和N=400，以研究不同网格分辨率下格式的性能。时间离散采用三阶TVDRunge-Kutta方法，CFL数取为0.4。将杂交WENO格式与传统WENO-JS格式的计算结果进行对比。在密度分布结果中，当网格数量N=100时，传统WENO-JS格式在激波附近出现了一定程度的数值振荡，虽然能够捕捉到激波的大致位置，但激波的过渡区域较宽，导致激波的分辨率较低。而杂交WENO格式由于采用了多种指示子来判断流场特性，在激波附近能够准确地切换到合适的格式，有效抑制了数值振荡，激波的过渡区域更窄，更清晰地捕捉到了激波的位置和强度。随着网格数量增加到N=200和N=400，两种格式的计算精度都有所提高，但杂交WENO格式始终保持着更好的激波捕捉能力，其计算结果更接近精确解，在激波和接触间断处的分辨率明显高于传统WENO-JS格式。在压力分布方面，当N=100时，传统WENO-JS格式在接触间断处的压力分布存在一定的偏差，不能准确地反映出压力的突变。杂交WENO格式通过指示子对间断区域的准确判断，在接触间断处能够更准确地捕捉到压力的变化，压力分布与精确解更为吻合。当网格加密到N=200和N=400时，杂交WENO格式在压力分布的计算上依然表现出色，能够更精确地模拟压力在间断处的变化，进一步验证了其在捕捉复杂流动间断方面的优势。通过该算例可以看出，基于多种指示子的杂交WENO格式在处理一维Euler方程的激波和接触间断等复杂流动结构时，相较于传统WENO格式具有更高的精度和更好的间断捕捉能力，能够更准确地模拟流场特性。5.1.2浅水波方程算例为分析新格式在处理水波传播和反射等问题时的精度和稳定性，选取一维溃坝问题作为浅水波方程的算例。初始条件为：在x<0.5区域，水深h=1，流速u=0；在x\geq0.5区域，水深h=0.1，流速u=0。计算区域为[0,1]，分别采用N=100、N=200和N=400的均匀网格进行离散，时间离散采用三阶TVDRunge-Kutta方法，CFL数取0.4。当N=100时，传统WENO格式在模拟水波传播过程中，波前出现了一定程度的数值振荡，并且在水波反射区域，水位的计算值与理论值存在较大偏差。而杂交WENO格式利用指示子准确判断水波传播过程中的流场变化，在波前能够有效抑制数值振荡，使波前更加光滑；在水波反射区域，能够更准确地捕捉到水位的变化，计算结果更接近理论值。随着网格加密到N=200和N=400，杂交WENO格式的优势更加明显，其计算结果在水波传播的各个阶段都能更精确地跟踪水位和流速的变化，数值振荡进一步减小，稳定性更高。在计算效率方面，由于杂交WENO格式在光滑区域采用了计算简单的线性方法，虽然在间断区域切换格式会增加一定的计算量，但总体计算时间相较于传统WENO格式并没有显著增加。在N=100时，杂交WENO格式的计算时间略高于传统WENO格式，但随着网格加密，两者的计算时间差距逐渐缩小。当N=400时，杂交WENO格式的计算时间仅比传统WENO格式多了约5%，考虑到其在精度和稳定性上的显著提升，这种计算效率的损失是可以接受的。通过该算例可知，基于多种指示子的杂交WENO格式在处理一维浅水波方程的水波传播和反射问题时，具有更高的精度和更好的稳定性，同时在计算效率上也具有一定的竞争力。5.1.3浅水中污染物输运方程算例以浅水中污染物输运方程为例，检验新格式在模拟污染物扩散和输运过程中的性能。算例初始条件设定为：在x<0.5区域，污染物浓度c=1；在x\geq0.5区域，污染物浓度c=0。扩散系数D=0.01，源项S=0。计算区域为[0,1]，采用N=100、N=200和N=400的均匀网格离散，时间离散采用三阶TVDRunge-Kutta方法，CFL数取0.4。当N=100时，传统WENO格式在模拟污染物扩散时，浓度分布出现了明显的数值振荡，尤其是在浓度梯度较大的区域，振荡更为严重，导致对污染物扩散的模拟不准确。杂交WENO格式通过多种指示子对浓度梯度和光滑度的判断，在浓度梯度较大区域能够自动切换到合适的格式，有效抑制了数值振荡，浓度分布更加平滑，更准确地反映了污染物的扩散情况。随着网格数量增加到N=200和N=400，杂交WENO格式的计算精度进一步提高，在不同时刻的污染物浓度分布与理论解的吻合度更高。在长时间模拟中，传统WENO格式由于数值振荡的积累，会导致污染物浓度的计算结果出现较大偏差，甚至可能出现不符合物理规律的情况。杂交WENO格式凭借其稳定的计算性能，能够在长时间模拟中保持对污染物输运过程的准确模拟，浓度分布始终符合污染物扩散的物理特性。通过该算例可以得出，基于多种指示子的杂交WENO格式在模拟浅水中污染物输运方程时，能够有效抑制数值振荡，提高计算精度，更准确地模拟污染物在浅水中的扩散和输运过程，为相关的环境科学研究和工程应用提供了更可靠的数值方法。5.2多维问题模拟5.2.1Burgers方程算例为展示新格式在多维复杂流场模拟中的能力和效果，选取二维Burgers方程算例。方程形式为：\frac{\partialu}{\partialt}+u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy}=0\frac{\partialv}{\partialt}+u\frac{\partialv}{\partialx}+v\frac{\partialv}{\partialy}=0初始条件设定为：u(x,y,0)=\sin(\pix)\cos(\piy)，v(x,y,0)=-\cos(\pix)\sin(\piy)。计算区域为[0,1]\times[0,1]，采用均匀网格离散，分别考虑50\times50、100\times100和200\times200三种网格分辨率。时间离散采用三阶TVDRunge-Kutta方法，CFL数取0.3。在50\times50的网格分辨率下，传统WENO格式在模拟流场时，速度场的等值线出现了明显的锯齿状，这表明存在较大的数值误差，尤其是在流场变化剧烈的区域，如涡旋中心附近，速度的计算值与理论值偏差较大。杂交WENO格式通过多种指示子对流场的判断，在流场变化剧烈区域能够自动调整格式，速度场的等值线更加光滑，更准确地描绘出了涡旋的形状和位置，在涡旋中心附近的速度计算值也更接近理论值。随着网格分辨率提高到100\times100和200\times200，杂交WENO格式的优势更加显著，能够更清晰地分辨出流场中的复杂结构，如小尺度的涡旋和速度梯度变化较大的区域，其计算结果与理论解的吻合度更高。在计算复杂流场的能量时，传统WENO格式由于数值误差的存在，计算得到的能量随时间的变化曲线出现了较大的波动，不能准确反映流场能量的真实变化。杂交WENO格式计算得到的能量曲线更加平滑，能够准确地跟踪流场能量的衰减过程，与理论分析的能量变化趋势一致。通过该算例可知，基于多种指示子的杂交WENO格式在处理二维Burgers方程这类多维复杂流场问题时，能够有效提高计算精度，更准确地模拟流场中的复杂结构和能量变化，展现出了在多维复杂流场模拟中的强大能力和良好效果。5.2.2Euler方程算例模拟多维Euler方程，分析新格式在处理复杂流动现象时的表现，如激波与边界层的相互作用。以二维激波与边界层相互作用算例为例，计算区域为[0,2]\times[0,1]，在x=0处设置为激波入口，激波强度为M=1.2，在y=0处设置为固体壁面，模拟边界层的形成。采用非均匀网格离散，在边界层和激波附近进行网格加密，分别考虑不同的网格加密程度。时间离散采用三阶TVDRunge-Kutta方法，CFL数取0.35。在较低的网格加密程度下，传统WENO格式在模拟激波与边界层相互作用时，激波的位置和形状出现了一定的偏差，边界层内的速度和压力分布也存在较大误差。杂交WENO格式利用指示子能够准确判断激波和边界层区域的流场特性，在激波附近有效捕捉激波的位置和强度，在边界层内能够更准确地模拟速度和压力的变化。在边界层内，传统WENO格式计算得到的速度剖面与理论的对数律分布存在明显差异，而杂交WENO格式计算得到的速度剖面更接近对数律分布，能够准确地反映边界层内的速度变化规律。在较高的网格加密程度下，杂交WENO格式依然能够保持良好的计算性能，对激波与边界层相互作用的复杂流动现象进行更精确的模拟。通过纹影图可以清晰地看到，杂交WENO格式计算得到的激波形状更加清晰、准确，激波与边界层的相互作用细节也能得到更好的展现。在压力分布方面，杂交WENO格式能够准确地捕捉到激波与边界层相互作用区域的压力突变和压力分布的变化，与实验结果和理论分析更为吻合。通过该算例表明，基于多种指示子的杂交WENO格式在处理多维Euler方程中激波与边界层相互作用这类复杂流动现象时，具有更高的精度和更好的模拟能力，能够为相关的空气动力学研究和工程应用提供更可靠的数值模拟结果。5.2.3浅水波方程算例研究多维浅水波方程的模拟结果，验证新格式在处理地形复杂的水波问题时的适应性和准确性。以二维溃坝问题为例，初始条件为在[0.3,0.7]\times[0.3,0.7]区域内水深h=1，流速u=v=0，其余区域水深h=0.1。计算区域为[0,1]\times[0,1]，采用非均匀网格离散，在水波传播和反射的关键区域进行网格加密，分别考虑不同的地形条件，如平底、有障碍物的地形。时间离散采用三阶TVDRunge-Kutta方法，CFL数取0.4。在平底地形条件下，当采用较低的网格加密程度时，传统WENO格式在模拟水波传播过程中，波前出现了数值振荡，并且在水波反射区域，水位的计算值与理论值存在一定偏差。杂交WENO格式通过指示子对水波传播和反射区域流场的准确判断，在波前有效抑制了数值振荡，使波前更加光滑，在水波反射区域能够更准确地捕捉到水位的变化，计算结果更接近理论值。随着网格加密程度的提高，杂交WENO格式的优势更加明显，能够更精确地模拟水波在平底地形上的传播和反射过程。在有障碍物的地形条件下，传统WENO格式在模拟水波与障碍物相互作用时，障碍物附近的水位和流速计算出现了较大误差，无法准确捕捉到水波的绕射和反射现象。杂交WENO格式利用指示子能够敏锐地感知到障碍物附近流场的复杂变化，在该区域自动切换到合适的格式，准确地模拟了水波与障碍物的相互作用，清晰地捕捉到了水波的绕射和反射波的传播，障碍物附近的水位和流速分布与理论分析和实验结果更为一致。通过该算例可以得出，基于多种指示子的杂交WENO格式在处理多维浅水波方程中地形复杂的水波问题时，具有良好的适应性和准确性，能够有效地模拟水波在不同地形条件下的传播、反射和绕射等现象，为水利工程和海洋科学等领域的相关研究提供了更有效的数值模拟工具。5.2.4浅水中污染物输运方程算例考察新格式在多维空间中模拟污染物输运的能力，分析其对复杂