变工况非定常流场测量的联合稀疏重构方法

《Energy and AI》：A Joint Sparse Reconstruction Method for Unsteady Flow-Field Measurement under Variable Operating Conditions

　　

摘要
涡轮机械中非定常流场的精确重构对于能源与推进系统的气动优化和状态监测至关重要，但时间分辨测量通常受限于硬件采样能力以及变工况下的基失配问题。本研究提出了一种多工况联合稀疏重构框架，融合了跨域数据同化、加权互相关准则和混合l_2,1-范数优化。基于不同攻角下流场间的相关性发展了一种通用联合本征正交分解（POD）基，以改善跨工况表示能力，同时行稀疏系数公式利用了各运行状态下共享的相干结构。该方法使用压气机叶栅在0°至10°五个攻角下的时间分辨粒子图像测速（PIV）数据进行验证。结果表明，所提出的基于POD的联合重构方法显著优于传统重构方法。在压缩比为0.6时，与传统基于POD的方法相比，均方误差降低了约89.6%。对于五工况包线，加权联合基能够捕获不同攻角状态下的公共结构，但在大攻角下由于更强的分离和工况特定的涡动力学特性，误差较大。一旦压缩比超过信息阈值，该方法对不同采样频率组合保持鲁棒性。所提出的框架为稀疏采样和变工况条件下高保真非定常流场重构提供了一种有效的物理信息数据同化策略。

**变工况非定常流场测量的联合稀疏重构方法：论文解读**

**1. 研究背景与问题**
现代能源系统，如先进推进系统、大型燃气轮机和超临界二氧化碳布雷顿循环，其核心部件（高负荷压气机与涡轮）内部的流动本质上是三维、非定常且强非线性的。叶尖泄漏涡、角区分离、尾迹脱落和大尺度流动分离等复杂二次流现象常常主导着气动场的相干结构。准确感知和重构这些非定常流场，对于气动优化设计、失速预警、智能状态监测、预测性维护以及数字孪生技术的开发至关重要。然而，在实际工业环境中获取高保真、高频的时空流场数据极其困难。传统非侵入式全局测量技术（如时间分辨粒子图像测速，PIV）虽能实现高空间分辨率，但受硬件限制、测量不确定性和高数据吞吐成本制约，难以同时提供足够高的时间分辨率。依据奈奎斯特-香农采样定理，为避免混叠并分辨高频湍流结构，采样频率至少需为感兴趣最高波动频率的两倍，这对传感器、成像系统、存储与后处理流程提出了严峻挑战。因此，实际测量通常在空间、时间或两者上均呈现稀疏性，数据稀疏问题已成为实时流场监测与数据同化的关键瓶颈。

**2. 研究目标与方法**
为克服硬件采样速率的物理限制，压缩感知（CS）、本征正交分解（POD）和动态模态分解等技术已被广泛用于流体测量数据的后处理与重构。然而，现有重构框架大多针对单一工况开发。在实际压气机中，流量、进气攻角及气动负荷的波动会导致流场变化；当攻角改变时，POD基可能无法充分表示另一工况下的流动结构，产生严重的基失配问题，本质上是领域偏移问题。针对该问题，研究人员提出了多工况联合稀疏重构框架。该框架构建跨工况的联合POD基以提供统一低维表示，并引入混合l_2,1-范数优化以利用模态系数的组稀疏性（行稀疏性），增强稀疏与噪声测量下的重构鲁棒性。研究采用了压气机叶栅水洞实验的时间分辨PIV数据库进行验证。数据库涵盖五个攻角工况（0°、2.5°、5°、7.5°、10°），每个工况包含1000个快照。实验在中低雷诺数（基于弦长的雷诺数Re_c=24,000）条件下进行，来流速度0.1 m/s，叶片弦长126.8 mm，栅距72 mm，安装角30°，叶栅为十二片C4型叶型，进口湍流度低于1%。PIV系统采用双脉冲Nd:YAG激光器（波长532 nm）、12位CCD相机（分辨率2072×2072像素），采样频率10 Hz，每个工况采集1000对图像。最终矢量间距约为1.15 mm×0.29 mm，图像采用多遍互相关算法处理，无效矢量经归一化中值检验后由局部插值替换。

**3. 关键技术方法**
该研究主要采用了以下关键技术方法：
- **加权联合相关POD基构建**：通过计算各工况空间相关矩阵的加权和构造联合相关矩阵，经特征值分解提取跨工况共享的联合POD基，实现对多工况流场结构的统一低维表示，权重系数可依据监测优先级调节。
- **混合l_2,1-范数联合稀疏优化**：采用混合l_2,1-范数替代独立l₁-范数最小化，在模态系数矩阵上施加组稀疏约束，使算法在不同工况间联合挑选共有相干结构模态，抑制任意删除重要模态，该优化策略与多测量向量（MMV）模型结合，提升多工况同步重构的鲁棒性。
- **多测量向量重构模型**：将各工况的稀疏观测方程装配为统一形式，利用联合感知矩阵对协同系数矩阵进行整体求解，充分发挥跨工况结构稀疏性的约束作用。
- **多工况实验数据库验证**：利用五攻角PIV数据构建联合基并评估重构性能，以压缩比（CR）为参数，系统考察方法在不同采样稀疏程度下的误差特性与鲁棒性。

**4. 研究结果**
**4.1 双工况基线重构性能（0°/5°）**
基于0°和5°双工况PIV数据，研究人员比较了所提出联合稀疏重构方法与常规重构方法的性能。结果显示，以联合POD基为稀疏表示基时，一旦CR超过约0.4，重构均方误差（MSE）保持在0.10~0.11的狭窄范围，表现出良好鲁棒性；而基于离散余弦变换（DCT）基的重构误差显著较高且不随CR呈收敛趋势。在CR=0.6时，传统POD基重构的MSE达0.96，而联合重构的MSE约为0.10，误差降低约89.6%，证明基于数据同化的联合稀疏框架在多工况流场恢复中具有显著优势。

**4.2 联合POD基的扩展**
研究者进一步构造了包含0°与5°双工况相关信息的联合POD基。与单独基于0°的POD基相比，联合POD基对5°工况的重构误差由约0.30降至约0.19，而0°工况的误差仅从0.103小幅升至0.115。这表明联合基在略微牺牲基线工况精度的前提下，显著改善了偏设计工况的表示能力。CR敏感性分析显示，低采样区（CR≤0.3）误差较高，存在临界信息阈值；当CR提高至0.3时，误差显著降低50%以上；CR达到0.6~0.95区间后误差趋于稳定（约0.18~0.20），说明误差下限主要受基函数表示能力限制而非采样不足。

**4.3 多工况联合重构与常规重构比较**
将框架推广至五工况包线（0°至10°）后，研究人员采用中心加权向量[0.5, 0.8, 1, 0.8, 0.5]构建加权联合相关矩阵并提取联合POD基。重构结果显示，0°、2.5°、5°、7.5°和10°五工况的MSE分别为0.083、0.145、0.158、0.384和0.297。小攻角和中攻角下重构精度较高，而大攻角误差增大，这与分离增强和工况特定的涡动力学特性相关。值得注意的是，7.5°误差略高于10°，研究人员将其归因于该过渡状态分离拓扑与相干涡组织的不规则切换。五工况情形下，低采样区（CR≤0.4）误差保持在1.0以上，表明需要更高采样信息才能提取可靠公共结构；当CR增至0.5后误差降至0.806，随后逐步下降至CR=1.0时的0.230，未出现明显收敛平台，说明多工况复杂性持续挑战联合基的表示能力。权重分配实验表明，增大中心工况权重可降低该工况误差（5°由0.158降至0.124），但以牺牲相邻工况精度为代价（0°增至0.125，10°增至0.401），体现有限表示空间内的精度权衡。

**4.4 采样策略鲁棒性**
为考察方法在实际测量条件下的稳定性，研究者评估了不同低采样频率组合和相位偏移对重构精度的影响。结果显示，当CR超过临界信息阈值（双工况约0.4、五工况约0.5）后，不同采样频率对和任意相位偏移下的MSE分布非常集中，方差极小，证明该方法在满足信息量要求的前提下对具体采样模式和通道相位配置不敏感，适应多通道硬件受限测量场景。

**5. 讨论与结论**
该研究提出的多工况联合稀疏重构框架通过联合相关POD基和l_2,1-范数优化，有效缓解了变工况下传统单工况重构的基失配问题。研究结论主要包括：(1) 基于POD的联合重构显著优于常规方法，CR=0.6时误差降低约89.6%；(2) 联合POD基通过融合跨工况相关性增强泛化能力，双工况下5°工况误差由0.30降至约0.19；(3) 五工况包线验证了加权联合基对公共结构的捕获能力，低攻角工况重构良好，高攻角误差较大；(4) 权重设置可实现目标工况精度与整体精度的权衡；(5) 重构精度强烈依赖CR，存在约0.4（双工况）和约0.5（五工况）的临界信息阈值，超过后对采样组合不敏感。整体而言，该框架为稀疏测量下高保真非定常涡轮机械流场重构提供了有效的物理信息数据同化策略。未来工作将聚焦自适应权重优化、三维重构扩展以及与物理信息学习方法的集成。