《Computer Physics Communications》：Implicit velocity correction schemes for scale-resolving simulations of incompressible flow: Stability, accuracy, and performance
编辑推荐：
尺度解析模拟（scale-resolving simulations）是准确捕捉非定常流动现象（如涡脱落、分离和层流—湍流转捩）所必需的，但其广泛应用仍受限于在物理相关时间区间内推进控制方程所需的高计算成本，尤其是对复杂几何外形周围的高雷诺数流动。计算工作量由
尺度解析模拟（scale-resolving simulations）是准确捕捉非定常流动现象（如涡脱落、分离和层流—湍流转捩）所必需的，但其广泛应用仍受限于在物理相关时间区间内推进控制方程所需的高计算成本，尤其是对复杂几何外形周围的高雷诺数流动。计算工作量由求解时间（time-to-solution）决定，即 Tsol = NΔt × TΔt，其中 NΔt 为时间步数，TΔt 为每步成本。本研究假定空间离散固定，以隔离时间推进格式对上述两个量的影响，仅对最大线性隐式时间步作一处例外说明。显式格式每步成本低，但受 Courant–Friedrichs–Lewy（CFL）条件约束，须采用小时间步长，因而 NΔt 很大。隐式格式放宽了这一限制并允许更大时间步，但会增加 TΔt；此外，增大时间步可能降低时间精度。因此，选择合适的时间积分策略需要在稳定性、精度与计算效率之间取得平衡。研究人员针对复杂高雷诺数曲面几何（其曲面施加严格的 CFL 限制）周围的不可压缩流动研究这一权衡。研究在速度校正格式（velocity correction schemes）类框架内开展，该类格式将不可压缩 Navier–Stokes 方程分裂为压力 Poisson 问题与速度分量的 Helmholtz 问题，是 Neko、ExaDG、Nek5000、NekRS 和 Nektar++ 等多个大规模不可压缩流动求解器的基础。多数实现采用半隐式（semi-implicit）形式，即扩散项隐式处理、对流项显式处理；该形式计算高效，但仍受显式对流项引起的 CFL 限制，对高阶离散和几何复杂区域该限制可能非常严苛。已有若干策略用于放宽该限制：全隐式（fully implicit）形式隐式处理对流项，导致非线性方程组，Baek 等采用 Aitken 松弛的定点迭代报告时间步稳定性提升可达两个数量级，但每步成本显著增加；算子积分因子分裂（operator-integration-factor-splitting，OIFS）方法通过在伪时间中对对流方程进行子步进（sub-stepping）来应对该限制，Sherwin 及其合作者在速度校正框架内应用该方法，可在引入与子步数成正比的额外成本的同时实现更大的外层时间步，半拉格朗日（semi-Lagrangian）形式通过追踪特征线允许时间步大幅增加，但主要针对层流研究；线性隐式（linear-implicit）方法通过单次 Picard 迭代线性化对流算子以避免非线性求解，以时变系统矩阵为代价扩展稳定性极限，能量稳定形式进一步增强稳定性，但每步需额外 Helmholtz 与 Poisson 求解，使线性求解器工作量实际翻倍。本研究聚焦子步进与线性隐式形式，旨在不付出非线性隐式求解器全部代价的前提下放宽 CFL 限制。与以往主要考察典型或层流构型的研究不同，研究人员采用统一的高阶谱/hp 单元离散，对复杂高雷诺数几何的隐式大涡模拟（implicit large-eddy simulation，iLES）系统比较这些格式。在固定空间近似（除线性隐式格式在 Δt ≥ 20 ΔtCFL 时的压力多项式阶数外）的条件下，研究人员隔离时间积分的影响，直接量化：（i）跨越初始瞬态的可达稳定性极限；（ii）作为时间步长函数的精度退化；（iii）准稳态区间内每对流时间单位的求解时间。据研究人员所知，这是首次在复杂几何上真实壁面解析大涡模拟（wall-resolved LES，WRLES）条件下对隐式速度校正策略开展的综合评估，为大规模不可压缩流动求解器中选择时间积分格式提供定量指导。同一三种格式在经典基准（圆柱绕流）上的补充比较见姊妹论文。对于湍流模拟，存在两个相关性能区间：统计量通过对准稳态湍流状态采样直至目标量（如阻力）收敛而获得，但模拟须先经历初始瞬态才能到达该区间，这引出另一性能目标，即缩短到达瞬态结束所需时间；在此阶段，较宽松的精度目标可能允许更大的稳定步长，前提是允许解在步长减小后、统计量采集前重新平衡，一旦开始采样，目标观测量对精度的要求即限制这一灵活性。本文成本测量在准稳态区间进行，对整个瞬态到采样流程计时仍属未来工作。作为代表性基准，研究人员考虑 Imperial Front Wing（IFW）几何，该几何为具有曲面与高雷诺数流动特征的多部件复杂翼型构型，完整几何已有实验研究并以 LES 数值研究，研究人员聚焦该几何的拉伸切片（eIFW），其保留基本流动物理同时降低计算复杂度；eIFW 的 WRLES 基准数据集见 Slaughter 等的工作，Ntoukas 等近期比较了该构型的显式与隐式 LES 方法。时间推进格式在开源谱/hp 单元框架 Nektar++ 中实现，该框架支持曲面几何上的高阶离散与可扩展并行执行；高阶格式的低数值耗散使本研究可采用隐式 LES 策略，其中亚格子尺度效应由离散化表示而不引入显式湍流模型。
研究背景方面，尺度解析模拟（scale-resolving simulations）是准确捕捉涡脱落、分离与层流—湍流转捩等非定常流动现象的必要手段，但其在物理相关时间区间上的高计算成本限制了广泛应用，尤其对复杂几何周围的高雷诺数流动。计算成本由求解时间 Tsol = NΔt × TΔt 决定。显式格式每步成本低但受 CFL 条件约束，时间步小、步数多；隐式格式放宽限制、允许大时间步，但提高每步成本且可能降低时间精度。因此时间积分策略需在稳定性、精度与效率之间权衡。现有多数速度校正（velocity correction）实现采用半隐式（semi-implicit）形式，扩散隐式、对流显式，仍受显式对流项引起的 CFL 限制，对高阶离散与复杂几何尤为严重。已有放宽策略包括：全隐式（fully implicit）形式导致非线性系统，Baek 等以 Aitken 松弛定点迭代报告时间步稳定性提升可达两个数量级，但每步成本大增；算子积分因子分裂（OIFS）方法以伪时间子步进（sub-stepping）对流方程，Sherwin 等在速度校正框架内实现更大外层时间步并引入与子步数成正比的额外成本，半拉格朗日（semi-Lagrangian）形式可大幅增大时间步但主要限于层流；线性隐式（linear-implicit）方法以单次 Picard 迭代线性化对流算子避免非线性求解，以时变系统矩阵为代价扩展稳定性极限，能量稳定形式进一步增强稳定性但每步额外 Helmholtz 与 Poisson 求解使线性求解工作量翻倍。以往研究多针对典型或层流构型，缺乏在真实壁面解析大涡模拟（WRLES）条件下对复杂几何的系统比较，这构成本研究的动因。
研究人员在统一的高阶谱/hp 单元离散框架下，对复杂高雷诺数几何的隐式大涡模拟（iLES）系统比较半隐式、子步进与线性隐式三种速度校正格式，固定空间近似（仅线性隐式在 Δt ≥ 20 ΔtCFL 时压力多项式阶数例外）以隔离时间积分影响，量化可达稳定性极限、精度随步长的退化以及准稳态区间每对流时间单位的求解时间。基准为 Imperial Front Wing（IFW）几何的拉伸切片（eIFW），格式在开源框架 Nektar++ 中实现。结论为：子步进格式在 20 倍显式时间步极限下保持稳定；线性隐式格式在等阶速度与压力空间下可达 100 倍；受稳定性约束的线性隐式加速比达 8.7 倍；力的一致性可能掩盖大的分离泡误差。其意义在于首次在复杂几何真实 WRLES 条件下综合评估隐式速度校正策略，为大规模不可压缩流动求解器的时间积分格式选择提供定量依据。
关键技术方法（不超过250字）：采用开源谱/hp 单元框架 Nektar++ 实现半隐式、子步进与线性隐式三种速度校正格式，在统一高阶离散下开展 eIFW 构型的隐式大涡模拟；以固定空间离散隔离时间积分影响，通过稳定性测试确定各格式可达时间步极限，以精度评估考察步长增大引起的偏差，并在准稳态区间测量每对流时间单位的求解时间；空间分辨率与统计收敛性另行验证，基准数据参照 Slaughter 等的 eIFW WRLES 数据集。
研究结果方面，稳定性极限部分表明，子步进格式在 20 倍显式时间步极限下仍保持稳定，线性隐式格式在采用等阶速度与压力空间时可达 100 倍显式极限，说明放宽 CFL 限制的潜力显著。精度评估部分表明，时间步增大导致精度退化，且力系数的一致性可能掩盖分离泡区域的大误差，提示仅以积分力验证不足。计算性能部分表明，受稳定性约束的线性隐式格式加速比可达 8.7 倍，但每步成本随格式不同而变化，需结合稳定性与精度综合权衡。空间分辨率与统计收敛性验证部分确认了所采用离散与采样足以支撑上述比较。
讨论与结论部分：研究人员指出湍流模拟存在两个性能区间，即从初始瞬态到准稳态的过渡阶段与统计采样阶段，前者可采用较宽松精度目标以允许更大稳定步长，后者受目标观测量精度约束；本文成本测量限于准稳态区间，完整瞬态到采样流程的计时为未来工作。研究强调在复杂几何真实 WRLES 条件下，格式选择应同时考虑稳定性极限、精度退化与每步成本，而非仅凭力的一致性判断；子步进与线性隐式格式均可在不同程度上放宽 CFL 限制，其中线性隐式在等阶空间下稳定性优势最大但需注意精度与成本代价。研究结论为：子步进格式在 20 倍显式时间步极限下保持稳定；线性隐式格式在使用等阶速度与压力空间时可达 100 倍；受稳定性约束的线性隐式加速比达 8.7 倍；力的一致性可能掩盖大的分离泡误差。
