多矩有限体积浸入边界法：原理、应用与展望.docxVIP

多矩有限体积浸入边界法：原理、应用与展望

一、引言

1.1研究背景与意义

在计算流体力学（CFD）领域，准确高效地求解复杂几何区域内的流动问题始终是核心关注点。从航空航天中飞行器的复杂外形绕流，到生物医学中人体血管内的血液流动，再到能源领域中发动机内部的燃烧流动等，这些实际工程和科学问题都涉及复杂几何区域内的流动现象。准确模拟这类流动，对于优化设计、提高性能、降低成本以及深入理解相关物理机制具有重要意义。例如，在航空航天领域，精确掌握飞行器复杂外形周围的流场特性，有助于优化飞行器的气动外形，减少阻力、提高升力，进而提升飞行效率和燃油经济性；在生物医学领域，深入了解血液在复杂血管结构中的流动情况，对于心血管疾病的诊断、治疗方案的制定以及医疗器械的研发有着关键作用。

传统上，解决复杂几何区域内流动问题常采用适体的曲线坐标系或非结构网格方法。适体曲线坐标系通过构建与几何边界贴合的坐标系统，使计算网格能够较好地适应复杂边界形状。非结构网格则具有更强的灵活性，能够根据几何形状自由地生成网格单元，无需像结构网格那样遵循规则的排列方式。然而，这两种方法都存在明显的局限性。适体曲线坐标系在处理复杂几何时，坐标变换的构建过程繁琐复杂，需要耗费大量的时间和人力，而且在计算过程中，由于坐标变换引入的额外计算，使得离散求解过程比笛卡尔网格复杂，增加了计算的难度和时间成本。非结构网格虽然在网格生成上具有更大的灵活性，但其质量往往受到几何复杂程度的显著影响。当几何形状极为复杂时，非结构网格容易出现网格扭曲、质量下降等问题，这不仅会降低计算精度，甚至可能导致计算结果的不收敛，无法得到有效的数值解。此外，在处理运动或变形边界条件时，这两种方法都面临着巨大的挑战。对于适体曲线坐标系和非结构网格，当边界发生运动或变形时，网格通常需要随边界变化进行动态调整。这涉及到复杂的网格重构算法，不仅计算量巨大，而且容易引入数值误差，进一步增加了计算的复杂性和不确定性。

浸入边界法（ImmersedBoundaryMethod，IBM）为解决上述问题提供了一种新的有效途径。与传统方法不同，浸入边界法在笛卡尔坐标上离散求解流体控制方程，通过在控制方程中添加相应源项来体现浸入边界的作用。这种方法无需生成复杂的贴体网格，避免了坐标变换和网格重构的难题，大大简化了计算过程，提高了计算效率。自1972年CharlesPeskin为模拟血液在可收缩的心脏瓣膜中流动而首次提出浸入边界法以来，该方法凭借其简单、高效的显著特点，在计算流体力学领域得到了广泛的应用和深入的发展。目前，浸入边界法已成功应用于模拟血液流动、湍流的直接数值模拟、多相流动以及流固耦合等诸多领域，展现出强大的生命力。

多矩有限体积法是依据多矩思想发展起来的新型流体力学数值格式，近年来在多个研究领域的数值模式构建中得到应用。其中，VSIAM3（Volume/SurfaceIntegratedAverageMulti-MomentMethod）格式通过使用压力投影算法，开发了可压、不可压流体的统一数值解法。数值实验结果表明，该算法具有准确、高效及鲁棒性强等特点。多矩有限体积法的独特优势在于，通过定义更多自由度，能够在单网格内实现高阶空间重构，从而提高数值模拟的精度。同时，VSIAM3方法使用不同类型的矩实现速度、压力在交错网格上的耦合，不同变量的同类矩均定义在同一网格点上，这种独特的耦合方式为浸入边界法的实施提供了便利条件。将多矩有限体积法与浸入边界法相结合，形成多矩有限体积浸入边界法，有望充分发挥两种方法的优势，为复杂几何区域内流动问题的求解提供更高效、更准确的数值模拟方法。这对于推动计算流体力学在各个领域的应用，解决实际工程和科学研究中的复杂流动问题具有重要的理论意义和实际应用价值。

1.2国内外研究现状

多矩有限体积法的研究方面，国内外学者取得了一系列重要成果。国外一些研究团队致力于拓展多矩有限体积法在不同物理场景下的应用。如[具体团队1]通过对多矩有限体积法的深入研究，将其应用于高超声速流动的模拟，成功捕捉到了复杂的激波结构和流动细节，验证了该方法在处理极端流动条件下的有效性。[具体团队2]则在多矩有限体积法的基础上，开发了适用于多相流模拟的数值算法，通过引入新的界面捕捉和相平衡模型，实现了对多相流中复杂界面行为和相间相互作用的准确模拟。在国内，也有众多学者对多矩有限体积法展开了深入研究。例如，[具体团队3]针对多矩有限体积法的数值稳定性和精度进行了系统分析，通过改进数值格式和优化算法参数，提高了该方法在复杂流动模拟中的可靠性和准确性。[具体团队4]则将多矩有限体积法与并行计算技术相结合，实现了大规模复杂流动问题的高效求解，显著提升了计算效率，为解决实际工程中的大规模计