水波弱隐身经由空间变换超材料实现

《Water》：Weak Cloaking in Water Waves via Spatial-Transformation Metamaterials

　　

当水波遭遇海洋结构物时会产生显著散射，进而影响海洋装备的稳定性与安全性。基于水波变换超材料（water-wave transformation metamaterials）的海洋结构物隐身为解决该问题提供了新思路。然而，在完美隐身所需的各向异性水深中，径向水深趋于无穷而周向水深在结构边缘趋于零，导致奇异性问题，并使结构附近的超材料参数趋于极端，难以制备。针对上述难题，研究人员基于空间变换超材料设计了一种用于水波弱散射的梳状超材料（comb-type metamaterial），实现对水波散射特性的优化控制。首先，研究人员为弱散射与完美隐身建立了统一的非线性空间变换关系，进而计算并分析了两种情形下的各向异性水深与等效重力加速度参数。随后采用有限元方法求解亥姆霍兹方程（Helmholtz equation），对比分析了弱散射与完美隐身情形下各向异性水深及超材料中的波场分布。结果表明：弱散射情形下径向与周向水深在结构边缘趋于常值，从而规避了完美隐身的奇异性；同时其水波控制性能与完美隐身相当。相较于完美隐身超材料，所设计的弱散射超材料在保持水波控制性能的同时显著降低了结构附近的水深参数，更易于工程制备。

研究背景方面，水波作为自由液面扰动现象具有折射、反射与共振等特征，海洋结构物引起的水波散射会改变流场并威胁装备安全。基于散射抵消法（SCM）与空间变换结合超材料（STM & MM）的隐身技术中，后者借助坐标变换下控制方程的形式不变性引导波前传播。然而典型完美隐身变换将物理点映射为有限区域，导致内边界处径向水深h_r→∞、周向水深h_θ→0的介质参数奇异性，无法物理实现，阻碍工程化。为此，研究人员提出弱隐身空间关系，以规避无限水深，并设计梳状超材料在保持波场控制性能的同时降低垂向尺寸。该文发表于《Water》。

关键技术方法上，研究人员采用浅水线性化推导得到以波面升高ζ(x₁,x₂;t)描述的亥姆霍兹方程；借助正交坐标变换?*=f(?)及其雅可比矩阵J，将空间变换吸收进介质参数，得到各向异性尺度因子τ_h=JJ^T/|J|与τ_g=|J|；为便于实现，令τ_g为常数并并入各向异性水深h_r、h_θ，用有效介质法将连续各向异性水深离散为K=25层交替各向同性子层a/b近似（层厚远小于波长）；以COMSOL Multiphysics 6.0有限元求解方程(10)，计算域10 m×10 m，入射波幅A₀=0.05 m、波长λ=3 m，外边界加完美匹配层（PML），监测半径r₀=8 m处的散射能衰减系数ε。

研究结果部分，按原文小节归纳如下。

在变换空间设计中，研究人员将原始空间环域α<r<γ映射到变换空间β<r^*<γ（α<β<γ），α=0对应完美隐身极限；该非线性变换使半径为β的结构散射场等价于半径为α的小结构，从而实现弱散射定义。

在各向异性介质分析中，由变换式(15)得τ_g弱散射取0.571、完美隐身取0.556；径向分量h_r/h₀随r单调增、周向分量h_θ/h₀单调减。内边界r=β处，弱散射h_r/h₀=0.0476、h_θ/h₀=21.0000，均有限；完美隐身则h_r/h₀=0、h_θ/h₀→∞，奇异性消除。有限元波场显示弱散射使半径2 m原结构总波场与小圆柱（α=0.5 m）相近，散射衰减系数ε_wi=ε_ss=0.16；完美隐身ε_pi≈0。

在弱散射超材料中，取K=25层，β=2 m、γ=4 m。最内层水深比τ_a(1)弱散射为17.98，完美隐身为38.61，降低约53.4%；次层降26.9%。简化各向异性参数下ε_ws=0.22，超材料ε_wm=0.25；完美隐身简化参数ε_ps=0.04，其超材料ε_pm=0.21，因内层奇异致离散近似偏差。远场r₀=8 m处弱散射超材料散射幅值全方向低于0.3，与完美隐身超材料（最大约0.35）同级，证明弱散射超材料以更缓水深梯度获可比抑制效果。

讨论与结论翻译部分，研究人员指出：基于空间变换超材料的水波隐身具保护离岸结构潜力，为避免完美隐身固有参数奇异性，提出弱隐身空间关系并依此设计梳状弱散射超材料。统一非线性变换框架中α<r<γ→β<r^*<γ，α=0退化为完美隐身；弱散射使半径β结构在散射特性上等价于半径α小结构。两情形τ_g近常值，h_r/h₀与h_θ/h₀径向趋势一致，但内边界处完美隐身h_r/h₀=0、h_θ/h₀→∞，弱散射皆有限，奇异性解除。有限元求解亥姆霍兹方程证实用弱散射各向异性介质可令原结构总波场等价于小结构。梳状超材料最内层水深比降逾50%，散射抑制性能与完美隐身超材料相近，垂向尺寸与水深梯度均缓和，提升空间隐身超材料可行性。未来可研隐身区与层厚等以降制备难度，并补实验验证。