双曲守恒律方程中自适应人工粘性高阶迎风格式的构建与应用.docxVIP

双曲守恒律方程中自适应人工粘性高阶迎风格式的构建与应用

一、引言

1.1研究背景与目的

双曲守恒律方程作为偏微分方程领域的重要分支，在现代科学与工程领域中扮演着举足轻重的角色。这类方程广泛应用于描述各种守恒物理量在时空中的演变，如流体力学中的质量、动量和能量守恒，交通流理论中的车辆密度守恒，以及电磁学中的电荷守恒等。以流体力学为例，双曲守恒律方程的典型代表——欧拉方程和纳维-斯托克斯方程，是研究流体流动现象的基础。从航空航天领域的飞行器空气动力学设计，到水利工程中的水流模拟，这些方程为理解和预测流体行为提供了关键的数学框架。在交通流研究中，双曲守恒律方程用于描述车辆在道路上的流动，帮助优化交通规划和管理，缓解交通拥堵。

尽管双曲守恒律方程在理论和应用上具有重要意义，但其数值求解面临着诸多挑战。由于方程的非线性特性，解在演化过程中可能产生间断，如激波和接触间断等，这使得传统的数值方法难以准确捕捉和处理这些复杂的流动现象。为了克服这些困难，高阶迎风格式应运而生。高阶迎风格式通过利用更多的邻域信息进行数值离散，能够在保证计算稳定性的同时，提高数值解的精度和分辨率，更好地捕捉流场中的复杂结构和间断。然而，高阶迎风格式在实际应用中也存在一些问题，如数值振荡和色散误差等，这些问题在处理强间断时尤为突出，可能导致数值解的非物理行为，影响计算结果的可靠性。

为了进一步提升高阶迎风格式的性能，本研究引入自适应人工粘性方法。自适应人工粘性能够根据流场的局部特征，动态地调整人工粘性系数，在流场光滑区域保持较小的粘性，以减少数值耗散对精度的影响；而在流场存在强间断或大梯度区域，自动增加粘性，有效地抑制数值振荡，保证计算的稳定性。这种自适应调节机制使得数值格式能够更好地适应流场的变化，提高对复杂流动现象的模拟能力。通过深入研究双曲守恒律方程的带有自适应人工粘性的高阶迎风格式，本研究旨在为相关领域的数值模拟提供更高效、准确和稳定的计算方法，推动科学研究和工程应用的发展。

1.2国内外研究现状

在双曲守恒律方程数值解法的发展历程中，众多学者做出了卓越贡献，推动了该领域的不断进步。自Courant、Friedrichs和Lewy等人在二十世纪三十年代对双曲守恒律方程解的存在性和唯一性进行开创性研究，并给出著名的CFL条件以来，数值解法的研究如雨后春笋般展开。早期的研究主要集中在低阶格式，如一阶迎风格式和中心差分格式。一阶迎风格式虽然具有简单稳定的优点，但其数值耗散较大，导致计算精度受限，尤其在处理复杂流场时，难以准确捕捉流场的精细结构。中心差分格式虽具有二阶精度，但在高对流占优的情况下，容易出现数值不稳定现象，解的振荡问题严重影响了计算结果的可靠性。

随着对计算精度要求的不断提高，高阶迎风格式逐渐成为研究热点。二阶迎风格式通过对一阶迎风格式的改进，在一定程度上提高了计算精度，减少了数值耗散。例如，基于Taylor展开的待定系数方法构建的二阶迎风格式，在处理一些简单流场时表现出较好的性能。然而，对于复杂的双曲守恒律方程，二阶迎风格式仍显不足。三阶迎风型格式进一步提高了离散精度，但格式的稳定性条件变得更加苛刻，且在实际应用中，其迎风性相对较差，需要更加精细的处理来确保计算的可靠性。QUICK格式，即对流运动的二次迎风插值格式，采用了迎风思想对界面值进行修正，在一定程度上改善了计算精度和稳定性，但在处理边界条件时，面临着复杂的问题，需要特殊的技巧来保证计算的准确性。

在自适应人工粘性的研究方面，近年来取得了显著进展。传统的人工粘性方法通常采用固定的粘性系数，这种方法在整个流场中施加相同的粘性，无法根据流场的局部特征进行动态调整。在流场光滑区域，过大的粘性会导致不必要的数值耗散，降低计算精度；而在强间断区域，固定的粘性系数又可能不足以抑制数值振荡，影响计算的稳定性。为了解决这些问题，自适应人工粘性方法应运而生。学者们提出了多种自适应人工粘性的判据和计算方法，例如根据波数识别方法，由流场密度的高阶导数求出局部等效无量纲波数，再将其与格式的修正无量纲波数进行对比，从而确定人工粘性系数的大小。当等效无量纲波数在格式的完全解析范围内时，人工粘性系数取零或很小的值，以减少数值耗散；当等效无量纲波数超出解析范围时，人工粘性系数根据一定规则逐渐增大，直至最大值，有效地抑制数值振荡。这种自适应调节机制使得人工粘性的施加更加合理，提高了数值格式对复杂流场的适应性。

国内的研究团队在双曲守恒律方程的数值解法方面也取得了丰硕成果。长安大学封建湖教授团队长期致力于双曲守恒律方程高分辨率熵稳定格式的研究，通过加入限制器，构造出了高分辨率熵稳定/熵相容的数值格式，并通过对一维Euler方程Lax激波管问题及爆震波问题等的数值试验，验证了所提方法的有效性