纠正SLD结冰参数：融合结冰风洞数据和数值先验的岭回归方法

《Aerospace》：Correcting SLD Icing Parameters: A Ridge Regression Method Fusing Icing Wind Tunnel Data and Numerical Priors

　　

为了探究FENSAP-ICE数值模拟在超冷大水滴（SLD）条件下的显著系统性偏差，研究人员提出了一种融合数值先验和领域知识辅助特征的岭回归校正方法，用于冰形几何参数。将FENSAP-ICE对八个冰形几何参数的预测作为数值先验纳入机器学习模型，将学习目标从“从零预测”转变为“校正系统性偏差”。基于SLD结冰物理机制构建了六个工程辅助特征，为小样本建模提供具有物理意义的可调维度。针对下表面结冰极限处的零冰情况，设计了一种结合逻辑回归（Logistic Regression）分类和岭回归（Ridge Regression）的两阶段模型。在29组NACA0012翼型SLD结冰风洞实验数据上通过留一交叉验证（Leave-One-Out Cross-Validation, LOOCV）进行评估，与FENSAP-ICE预测相比，改进系统将总冰面积（Atotal）的对称平均绝对百分比误差（sMAPE）从79.47%降低到34.65%，下表面冰角（θ_l）从105.19%降低到28.08%，上表面结冰极限（x_u）从61.75%降低到35.48%，平均冰厚（t_avg）从42.12%降低到20.89%。消融实验进一步揭示，仅引入数值先验即可将预测误差降低约10个百分点，成为主要性能驱动因素。所提方法具有低计算成本和强物理一致性，为小样本条件下的SLD冰形预测提供了实用的偏差校正框架。此外，为解决小样本交叉验证可能带来的乐观偏差，还额外采用了嵌套交叉验证（nested cross-validation）和多个线性基线模型来验证所提校正方法的鲁棒性和相对竞争力。

**论文解读：基于岭回归和数值先验的SLD结冰冰形偏差校正方法**

**研究背景与问题**

飞机结冰是航空安全的持续威胁，尤其在超冷大水滴（SLD，即中值体积直径MVD超过50μm的过冷水滴）条件下，水滴在撞击前会发生显著变形、破碎和飞溅，导致冰形预测更为复杂。Broeren等人指出，一旦冰层扩展到防冰系统保护区域之外，几分钟内升力损失可超过30%，凸显了准确预测SLD冰形的紧迫性。当前工程实践中，主流数值结冰代码如FENSAP-ICE在SLD条件下存在系统性低估，尤其是上下表面结冰极限，直接采用会导致防冰区域设计不足，产生严重残留结冰。现有物理信息机器学习方法主要通过嵌入控制方程约束或领域迁移来缓解实验数据稀缺问题，但尚未系统应用于校正SLD数值模拟偏差。为此，研究人员开展本研究，将FENSAP-ICE预测作为数值先验，利用岭回归学习其与结冰风洞测量之间的系统性偏差，以低计算成本校正关键冰形参数。该论文发表在《Aerospace》上。

**研究内容与结论**

研究人员提出了一种基于岭回归的改进SLD冰形预测系统，融合FENSAP-ICE数值先验和领域知识辅助特征，在29组NACA0012翼型SLD结冰风洞实验数据上通过留一交叉验证进行评估。结果表明，该系统显著降低了关键冰形参数的预测误差：总冰面积（Atotal）的sMAPE从79.47%降至34.65%，下表面冰角（θ_l）从105.19%降至28.08%，上表面结冰极限（x_u）从61.75%降至35.48%，下表面结冰极限（x_l）从59.00%降至37.21%，平均冰厚（t_avg）从42.12%降至20.89%。消融实验表明，仅引入数值先验即可降低约10个百分点误差，是主要性能驱动；工程辅助特征进一步改善局部参数。嵌套交叉验证和线性基线模型对比证实了校正的鲁棒性和相对竞争力。该方法具有低计算成本和强物理一致性，适用于小样本条件下的SLD冰形预测，可直接嵌入现有结冰模拟流程，用于飞机防除冰系统初步设计。

**主要技术方法**

研究人员采用以下关键技术方法：（1）数值先验：将FENSAP-ICE对八个冰形几何参数（上表面最大厚度t_u、下表面最大厚度t_l、总冰面积A_total、上表面结冰极限x_u、下表面结冰极限x_l、上表面冰角θ_u、下表面冰角θ_l、平均冰厚t_avg）的预测作为先验输入机器学习模型。（2）工程辅助特征：基于SLD结冰物理机制构建六个无量纲工程特征，涵盖水滴惯性、热力学冻结状态、破碎趋势及上下表面不对称性。（3）两阶段模型：针对下表面结冰极限x_l，先使用逻辑回归分类判断是否存在可测量冰，再对存在冰的情况用岭回归预测。（4）岭回归核心算法：采用L2正则化抑制过拟合，通过留一交叉验证网格搜索最优正则化参数λ∈[10^-2, 10²]。（5）验证策略：采用留一交叉验证（LOOCV）评估，并额外使用嵌套交叉验证和多个线性基线模型（Lasso、普通最小二乘、弹性网络）检验鲁棒性。样本来源：中国空气动力研究与发展中心FL-16结冰风洞的NACA0012翼型（1 m弦长）29组实验条件，攻角0°–4.9°，静温253.15–268.55 K，液态水含量0.21–0.50 g/m³，MVD 19–526 μm，结冰时间120–1350 s。

**研究结果**

**3. Results and Analysis**

**3.1 改进系统与FENSAP-ICE的sMAPE对比**
通过留一交叉验证，改进系统在总冰面积（A_total）、下表面冰角（θ_l）、上表面结冰极限（x_u）、下表面结冰极限（x_l）和平均冰厚（t_avg）上的sMAPE均显著低于FENSAP-ICE原始预测，分别从79.47%降至34.65%、105.19%降至28.08%、61.75%降至35.48%、59.00%降至37.21%、42.12%降至20.89%。对于上表面最大厚度（t_u）、下表面最大厚度（t_l）和上表面冰角（θ_u），改进系统与FENSAP-ICE精度相当或略有提升。

**3.2 散点图分析**
各参数预测值与实验值的散点图显示，大多数点紧密聚集在理想对角线附近，无结构性偏差。总冰面积和平均冰厚的预测鲁棒性最为突出，仅个别案例偏离超过一个标准单位。局部厚度参数（t_u、t_l）在高LWC或长结冰时间条件下稍有偏差，但整体可接受。结冰极限参数（x_u、x_l）中，x_l的散点呈现均匀分布，两阶段策略有效避免了负值或近零值。冰角参数（θ_u、θ_l）中，θ_u存在轻微低估倾向，θ_l偏差较小，符合下表面冰角受翼型约束更严格的物理现实。

**3.3 消融实验**
以t_u、x_u和t_avg为代表，比较三种特征配置：M1（仅原始条件参数）、M2（M1+FENSAP先验）、M3（M2+工程辅助特征）。M2相比M1，t_u的sMAPE降低5.71个百分点，x_u降低4.18个百分点，t_avg降低10个百分点，证明数值先验是主要性能驱动。M3相比M2，t_u进一步降低，x_u降低至49.18%，t_avg改善较小，表明工程特征对局部参数影响更显著。

**3.4 扩展验证与鲁棒性分析**
嵌套交叉验证显示，标准LOOCV的乐观偏差可忽略，嵌套与标准sMAPE基本一致。与Lasso、普通最小二乘（OLS）和弹性网络（ElasticNet）对比，岭回归在多数参数上保持竞争力或更优。岭回归系数稳定性分析表明，大多数参数的系数在交叉验证折间保持符号一致，但部分参数（如A_total、x_u、x_l的某些FENSAP先验系数）变化较大，需谨慎解释。残差分析显示残差围绕零散点分布，无明显结构。

**讨论与结论**

研究人员在讨论部分指出，改进系统通过校正结冰极限和总冰面积，缓解了原始FENSAP-ICE预测中防冰区域设计不足的风险。由于校正是一种计算成本极低的后处理步骤，可直接嵌入现有结冰模拟流程，适用于飞机防除冰系统的初步设计。研究受限于实验样本数量（29组）和速度条件（仅61 m/s和95 m/s两个水平），未来工作将扩展数据规模、引入不确定性量化，并探索更直接捕捉SLD微物理过程的特征。

**结论翻译**：针对超冷大水滴（SLD）条件下数值结冰模拟的系统性偏差，本文提出了一种基于岭回归和物理先验的改进冰形预测系统。在29组NACA0012翼型结冰风洞实验数据上通过留一交叉验证评估，主要结论如下：（1）数值先验的引入是主要性能驱动因素；消融实验表明，仅此即可将平均冰厚预测误差降低约10个百分点。在此基础上构建的六个工程辅助特征可进一步降低局部冰形参数（如上表面最大厚度）的误差约2个百分点。两者分别实现基线校正和局部精化。（2）与原始FENSAP-ICE预测相比，改进系统降低了所有关键参数的sMAPE：总冰面积A_total从79.47%降至34.65%，下表面冰角θ_l从105.19%降至28.08%，上表面结冰极限x_u从61.75%降至35.48%，下表面结冰极限x_l从59.00%降至37.21%，平均冰厚t_avg从42.12%降至20.89%。（3）本研究受限于实验样本数量和覆盖的速度条件范围，评估可能包含乐观偏差。未来工作将扩展数据规模、引入不确定性量化，并探索更直接捕捉SLD微物理过程的特征。（4）通过嵌套交叉验证、与其他线性基线模型的对比以及系数稳定性分析，证实了所提校正的鲁棒性。嵌套评估显示仅存在微小乐观偏差；校正方法在多数目标参数上优于或持平基线模型；大部分岭回归系数在交叉验证折间保持符号一致，少数系数幅度变化较大，需谨慎解释。（5）通过校正结冰极限和总冰面积，所提方法缓解了原始FENSAP-ICE预测中防冰区域设计不足的风险。由于校正是一种计算成本极低的后处理步骤，可直接嵌入现有结冰模拟流程，因此非常适用于飞机防除冰系统的初步设计。