基于PHG平台的不可压流体并行自适应有限元算法：理论、实现与应用.docxVIP

基于PHG平台的不可压流体并行自适应有限元算法：理论、实现与应用

一、引言

1.1研究背景与意义

不可压流体计算在众多科学与工程领域中扮演着举足轻重的角色，对推动现代科技的进步发挥着关键作用。在航空航天领域，飞机、火箭等飞行器在飞行过程中，周围的空气可视为不可压流体。精确计算不可压流体的流动特性，如升力、阻力等参数，对于飞行器的设计至关重要。通过优化飞行器的外形设计，使其在不可压流体环境中具有更优异的空气动力学性能，能够提高飞行效率、降低能耗并增强飞行安全性。在汽车工程领域，汽车在行驶时与周围空气的相互作用涉及不可压流体计算。研究汽车周围的流场分布，有助于减少空气阻力，降低燃油消耗，同时提升汽车的行驶稳定性和操控性能，为汽车的轻量化和节能设计提供重要依据。在能源领域，风力发电作为一种清洁能源，风电机组的叶片设计需要深入了解不可压流体（风）的流动规律。通过对叶片周围流场的精确计算，优化叶片形状和结构，能够提高风能捕获效率，增加发电功率，降低发电成本。在水利工程中，大坝、桥梁等水工建筑物的设计需要考虑水流的作用。利用不可压流体计算，分析水流对建筑物的压力、冲击力等，确保建筑物在复杂水流条件下的结构安全和稳定性。

随着科学技术的不断发展，各领域对不可压流体计算的精度和效率提出了更高的要求。传统的计算方法在面对复杂的几何形状和物理过程时，往往难以满足实际需求。自适应有限元算法作为一种高效的数值计算方法，能够根据计算区域的物理特征和误差分布，自动调整网格的疏密程度，从而在保证计算精度的同时，有效提高计算效率。然而，随着计算规模的不断扩大和问题复杂度的增加，并行计算成为实现高效不可压流体计算的必然趋势。

PHG平台作为我国自主研发的并行自适应有限元程序开发平台，具有重要的战略意义和实用价值。它为不可压流体并行自适应有限元算法的研究和应用提供了强有力的支撑，能够有效解决大规模科学计算中的诸多难题。PHG平台具备支撑面向万亿至十亿亿次级超级计算机的并行有限元应用程序开发的能力，并且支持各类国产高性能计算机。这使得基于PHG平台开发的并行有限元程序能够充分利用国产超级计算机的强大计算能力，实现大规模不可压流体问题的高效求解。在实际应用中，PHG平台已在多个领域展现出显著优势。在集成电路设计领域，基于PHG平台开发的并行有限元程序成功实现了大规模三维集成电路热应力的大规模并行计算，为三维集成电路的设计提供了高精度的应力计算模型，有力推动了我国集成电路技术的发展。在电磁模拟领域，基于PHG平台研制的电磁模拟软件FEMAG，能够实现近2万处理器核、愈8亿自由度的数值模拟，为大型变压器铁损模拟研究提供了重要工具，受到应用合作单位的高度评价。

1.2国内外研究现状

在不可压流体有限元算法方面，国内外学者取得了丰硕的研究成果。早期，有限元方法主要应用于简单的不可压流体问题，随着理论的不断完善和计算技术的发展，逐渐能够处理复杂的几何形状和物理现象。在稳定性分析方面，学者们提出了多种方法来确保有限元格式的稳定性，如混合有限元方法、稳定化有限元方法等。混合有限元方法通过引入额外的变量，改善了压力和速度的耦合关系，有效提高了格式的稳定性；稳定化有限元方法则通过添加适当的稳定项，抑制了数值振荡，保证了计算结果的可靠性。在误差估计方面，基于残差的误差估计器得到了广泛研究和应用。这种误差估计器通过分析计算结果的残差，能够准确地评估计算误差，并为自适应网格加密提供依据。自适应策略的研究也取得了重要进展，包括基于误差指标的自适应策略、基于后验误差估计的自适应策略等。这些自适应策略能够根据计算区域的误差分布，自动调整网格的疏密程度，从而在保证计算精度的同时，提高计算效率。

在PHG平台的应用方面，国内研究取得了一系列重要成果。在中国科学院的相关研究中，基于PHG平台开发了大规模三维集成电路热应力的大规模并行计算工具，成功应用于三维集成电路载流子沟道迁移率计算、TSV及互连线寿命分析、水冷结构稳定性分析等关键领域。该工具在神威太湖之光超级计算机1024CPU核组上进行了并行强可扩展性和弱可扩展性测试，强可扩展性效率达到72.42%，弱可扩展性效率为81.4%，展示了PHG平台在大规模科学计算中的强大能力。国内还基于PHG平台研制了国内第一个可在数万处理器核上运行的具有“线性”计算复杂度的并行寄生参数提取工具（ParAFEMImp和ParAFEMCap），以及我国唯一一个自主研发的支持超大规模并行计算和h-p自适应的电磁模拟软件FEMAG。这些成果表明，PHG平台在推动我国高性能计算应用和算法研究方面发挥了重要作用。在国际上，虽然类似的并行自适应有限元平台也有一定的发展，但PHG平台凭借其对国产计算机体