基于VOF的数值波浪水槽中采用EnKF方法的序贯水波重构

《Journal of Ocean Engineering and Science》：Sequential water wave reconstruction in VOF-based numerical wave tanks with the EnKF approach

　　

现有的相分辨波浪重构方法主要依赖势流理论，这严重限制了它们捕捉强非线性现象（如波浪破碎动力学）的能力。为突破这一瓶颈，本研究开发了一个鲁棒的序贯多相数据同化框架，该框架将集合卡尔曼滤波（EnKF）与基于VOF的CFD求解器相结合。一个关键特征是在同化前进行两相位表示转换，将不连续的原始VOF场映射为光滑的水平集符号距离场，从而增强模态基的紧凑性并抑制非物理的数值振荡。随后，对光滑界面和速度场的联合集合应用多变量本征正交分解（多元POD），将高维状态空间压缩为具有统一模态系数的低维子空间，这些系数自然保留了跨变量的物理耦合和运动学一致性。为缓解重复稀疏观测下的集合塌缩，在降阶系数空间中嵌入了松弛到先验扰动膨胀（RTPP）方法，利用CFD验证的先验运动学作为物理锚点，沿物理一致轨迹维持样本多样性。该框架通过三个代表性波浪条件下的孪生实验进行验证：五阶Stokes波、JONSWAP不规则波以及浅水斜坡上的极端卷破波。参数分析确认了其对观测噪声和传感器稀疏性的鲁棒性。减小同化间隔有效抑制了CFD求解器中快速非线性模型偏差的累积，确保了直至严重波浪破碎起始的高结构保真度。这一集成策略为数值波浪水槽与实际传播状态的同步提供了一条精确、物理一致的路径。

**论文解读文章**

**研究背景与问题**
当前，相分辨波浪重构技术主要依赖势流理论，其局限性在于无法准确捕捉波浪破碎等强非线性现象。然而，在船舶与海洋工程中，基于纳维-斯托克斯（Navier-Stokes）方程的高保真粘性流求解器（如计算流体动力学CFD）是模拟这些复杂动力学的关键，但其与真实物理状态之间的差异导致数值波浪水槽难以实现确定性复现。现有数据同化方法多基于势流模型，难以处理多相流中VOF（流体体积法）场的不连续性。因此，研究人员旨在开发一种集成EnKF（集合卡尔曼滤波）与VOF-based CFD求解器的序贯数据同化框架，解决波浪重构中的非线性瓶颈，实现高保真相分辨波浪场重建。该研究发表于《Journal of Ocean Engineering and Science》。

**主要关键技术方法**
研究人员采用以下关键技术：1) 构建初始波浪参数集合，通过随机扰动目标波形参数（如波高）生成N=50个成员；2) 两相位表示转换，将VOF场转换为连续的水平集（level set）符号距离场，再在重构后逆转换回VOF场；3) 多变量POD（本征正交分解），对水平集和速度场的联合集合进行奇异值分解，提取共享模态系数的低维子空间（r=9或更多）；4) 在降阶系数空间内执行EnKF分析，更新模态系数，并结合RTPP（松弛到先验扰动膨胀）方法维持集合多样性；5) 将更新后的系数映射回物理空间，更新CFD求解器的内部场。所有验证基于数值波浪水槽的孪生实验，无物理样本队列。

**研究结果**

**3.1 规则波重构**
针对五阶Stokes波规则波，研究人员通过孪生实验设置真实波高H_t=0.13 m、基线波高H_b=0.1 m，构造初始集合并执行单步同化（t₀=20 s）。采用99%累积方差阈值截断模态（r=9），利用10个均匀分布波高计（观测噪声σ_o=0.002 m）进行EnKF更新。结果表明，同化后DA均值与真实波面几乎完全重合，且速度场（u_x、u_z）重构误差极低，证实了多变量POD的跨变量耦合机制。对比直接使用VOF场的方案，水平集辅助方法在相同能量阈值下所需模态更少且重构误差更低（L₁范数）。参数敏感性分析显示，观测噪声增大或传感器间距增大（至6 m）时误差仍可控，波陡变化时重构误差呈U形趋势，在中间陡度（0.0337）达到最优，验证了框架的鲁棒性。

**3.2 不规则波重构**
针对JONSWAP不规则波（峰值周期T_p=2.5 s，有效波高H_s=0.075 m），研究人员设置真实与基线使用完全随机相位，集合成员也随机化相位。序贯同化间隔Δt_DA=1 s，采用20个波高计，RTPP松弛因子κ=0.5。在6个代表性时刻（(t-t₀)/T_p=0,2,4,6,8,10）的空间剖面显示，DA均值逐步跟踪真实波面，尤其在远场区域误差显著减小；RTPP膨胀后的集合保持光滑物理结构，避免传统膨胀引入的数值噪声。速度场重构在近场较弱，但随序贯更新逐渐改善。L₁误差时间演化表明，在初始快速收敛后，误差保持有界，无发散趋势，证明框架在随机边界下实现稳定跟踪。

**3.3 卷破波重构**
在浅水斜坡（初始水深0.4 m，坡长16 m）上，研究人员采用余弦波（周期T=5 s，H_t=0.15 m，H_b=0.12 m）生成卷破波。动态波高计集群（间距Δx=0.5 m，Δz=0.1 m）自适应跟踪波峰，同化间隔Δt_DA=0.25 s。从t=10 s至12 s（9个步骤）的序列显示，DA均值在卷破前（t=12 s）精确重构了翻卷唇和封闭空气囊的复杂拓扑，速度场（u_x、u_z）与真实几乎一致。对比无RTPP的控制实验，无RTPP时因集合塌缩导致最终时刻严重失真，证实RTPP是防止滤波发散的关键。另外，比较Δt_DA=0.5 s与0.25 s，较小间隔有效抑制了非线性模型偏差，在破碎前保持高保真度。

**总结与讨论**
研究人员在结论部分指出，该研究成功开发并验证了一个鲁棒的序贯多相数据同化框架，用于VOF-based数值波浪水槽的相分辨波浪重构。通过集成EnKF与有限体积纳维-斯托克斯求解器，框架克服了传统势流方法在捕捉强非线性波浪破碎动力学方面的长期局限。关键创新在于两相位表示转换机制，将高度不连续的VOF场映射为连续光滑的水平集场，显著增强模态基紧凑性并消除数值截断噪声；多变量POD方法使速度与界面场共享统一模态系数，自然保持跨变量物理耦合和运动学一致性；RTPP膨胀方案在降阶子空间中利用CFD验证的先验运动学维持集合方差。框架在三种波浪条件（五阶Stokes规则波、JONSWAP不规则波、浅水斜坡卷破波）下通过孪生实验验证，参数分析确认其对观测噪声和传感器稀疏性的鲁棒性。减小同化间隔有效抑制了CFD求解器中非线性模型偏差的快速累积，确保了直至严重波浪破碎起始的高结构保真度。该工作为构建与真实传播状态同步的高保真“数字孪生”波浪水槽提供了物理一致途径。未来研究将集中在：1) 集成降阶建模或高性能并行计算加速前向CFD模拟；2) 在物理实验室波浪水槽中进行实验验证；3) 将框架推广至三维波浪场，实现多向海况的高保真重构。