《Journal of Fluid Mechanics》：Searching for invariant solutions to wall-bounded flows using resolvent-based optimisation
编辑推荐：
研究人员提出了一种基于Galerkin投影的新型变分优化框架，用于在满足不可压缩性和无滑移边界条件（Boundary Conditions）的约束下，寻找Navier-Stokes方程的精确非线性解。针对传统压力泊松方程（Pressure-Poisson Eq
研究人员提出了一种基于Galerkin投影的新型变分优化框架，用于在满足不可压缩性和无滑移边界条件（Boundary Conditions）的约束下，寻找Navier-Stokes方程的精确非线性解。针对传统压力泊松方程（Pressure-Poisson Equation）求解中难以同时满足Neumann边界条件和Dirichlet边界条件，导致壁面残差非零及速度场偏离切空间的问题，该方法通过将速度场和局部残差场在散度自由且正交的基函数上进行展开，从根本上消除了压力梯度项的计算。研究人员利用基于层流剖面的分辨析（Resolvent Analysis）生成的响应模态（Response Modes）作为投影基，结合伪谱法（Pseudo-spectral Method）计算非线性项，成功重构了旋转平面库埃特流（Rotating Plane Couette Flow, RPCF）在Re=50时的平衡态解，并在Re=400时发现了周期解。该方法的全局残差（Global Residual）可降至10-12量级，为计算流体力学中难以捕获的非线性相干结构提供了新的数值途径。
研究背景与意义
寻找Navier-Stokes方程的精确解对于理解湍流中的相干结构和动力学行为至关重要。然而，传统的牛顿类方法（如Newton-GMRES-hookstep）虽然收敛快，但在高维参数空间中极易因初始猜测偏离而发散，且涉及大型线性算子的求解，计算成本极高。此外，在数值迭代过程中，如何严格保证速度场的不可压缩性（Divergence-free）和无滑移边界条件是长期存在的难点。现有方法通常分别求解速度和压力泊松方程，这往往导致边界条件的冲突，使得残差无法真正归零。因此，研究人员致力于开发一种更为稳健且能严格满足物理约束的优化算法，以期在复杂的流动状态（如旋转湍流）中发现新的精确解。该研究发表于《Journal of Fluid Mechanics》。
关键技术方法
研究人员采用了一种基于Galerkin投影的变分优化策略。首先，构建一个满足不可压缩性和无滑移边界条件的正交基函数集合，该基函数由分辨析（Resolvent Analysis）中的左奇异向量（即响应模态）构成。其次，将速度场和残差场在该基函数上进行傅里叶展开，通过投影消除压力项。在计算非线性对流项时，采用伪谱法（Pseudo-spectral Method）结合3/2去混叠规则以提高效率。优化过程使用L-BFGS（Limited-memory Broyden-Fletcher-Goldfarb-Shanno）拟牛顿算法结合Hager-Zhang线搜索来最小化全局残差。验证案例选用了旋转平面库埃特流（RPCF），这是一种在反气旋旋转下具有有限雷诺数线性不稳定的经典模型。
研究结果
一般特征与约束处理
研究人员详细分析了传统优化方法的缺陷：标准压力泊松方程无法同时兼容Neumann和Dirichlet边界条件，导致壁面残差不为零，破坏了梯度的下降方向。为解决此问题，研究人员引入了散度自由的基函数 ψkm，该基函数天然满足无滑移边界条件。由于基函数的散度为零，在进行内积运算时，压力梯度项的贡献自动为零（??p,ψ?=0），从而在不求解压力场的情况下实现了速度与压力的完全耦合更新。
Galerkin投影与变分优化
研究人员建立了速度场 u和残差场 r的谱展开式。通过引入Leray投影算子，将Navier-Stokes方程转化为仅关于速度的演化方程。优化目标函数定义为残差的范数 R。通过对系数 akm求导，得到了残差关于系数的梯度表达式。由于基函数的正交性，推导出的梯度表达式中不包含任何压力相关的项，这使得优化过程在数学上更加简洁且物理上自洽。
分辨析与模态构造
为了获得符合要求的基函数，研究人员采用了分辨析方法。通过对线性化的Navier-Stokes算子进行奇异值分解（SVD），提取出响应模态 ψkm和强制模态 ?km。这些模态构成了完备的正交基。研究发现，尽管分辨析通常在非零频率下定义，但在零频（k=0）下对算子进行SVD同样能得到有效的平均流基函数。奇异值的快速衰减特性允许研究人员使用低秩近似来大幅降低计算维度。
旋转平面库埃特流（RPCF）
RPCF的控制方程在Navier-Stokes基础上增加了科里奥利力项，其无量纲参数为雷诺数 Re和旋转数 Ro。研究发现，当 Ro=0.5时流动最具线性不稳定性（Recrit≈20.7）。随着雷诺数的增加，流动经历了一系列分岔：从层流转捩为平衡态（如S1-S3），再分岔为周期解，最终进入混沌状态。DNS（直接数值模拟）结果显示，在 Re=50时存在多个具有不同展向卷数的稳定平衡态分支。
精确解的寻找
在 Re=50,Ro=0.5的条件下，研究人员成功重构了已知的平衡态解S1，优化后的残差低于 10?12，与DNS结果高度吻合。随后，通过改变初始扰动的展向波数，研究人员发现了新的平衡态解S3、S4(o)和S5(o)，证明了该优化框架在探索状态空间方面的鲁棒性。在 Re=400时，研究人员从含噪声的初始猜测出发，成功捕捉到了一个周期为 T≈25.05的周期解，其动力学特征表现为流向卷的周期性收缩与膨胀。
优化性能评估
相比于传统的牛顿法，变分优化虽然在接近最小值时呈现线性收敛（而牛顿法为二次收敛），但其单次迭代成本极低（仅需一次矩阵-向量乘法）。在 Re=400的高分辨率计算中，该方法表现出了极高的内存效率。尽管在处理更高雷诺数时面临存储响应模态的内存挑战，但该方法的鲁棒性使其成为混合策略（前期变分优化+后期牛顿法细化）的理想前端。
优化的条件数分析
研究人员推导了残差函数的Hessian算子，发现其谱特性与分辨析算子的奇异值直接相关：Hessian的特征值等于分辨析奇异值的平方的倒数（μωm=σωm?2）。这意味着，截断分辨析展开中最小奇异值对应的模态（即最活跃的响应模态），实际上是在去除Hessian中最大特征值带来的病态方向。这一理论发现不仅解释了为何基于分辨析的基函数能够加速收敛，也为后续的预处理技术提供了理论依据。
结论
研究人员成功开发了一种基于Galerkin投影的变分优化方法，该方法通过利用分辨析生成的响应模态作为基函数，完美解决了Navier-Stokes方程求解中物理约束难以满足的难题。该方法在旋转平面库埃特流中成功找到了高精度（10?12）的平衡态和周期解，验证了其在寻找复杂流动精确解方面的有效性和稳健性。研究进一步从理论上揭示了分辨析模态与优化问题条件数之间的内在联系，指出截断大奇异值模态是改善Hessian矩阵条件数的关键。这项研究不仅为计算流体力学提供了一种强大的新型数值工具，也深化了对流动稳定性与优化动力学关系的理解。
