comsol水声超表面反射系数与反射相位计算。
打开模型树先给几何结构来点硬核配置。假设咱们设计的是锯齿状超表面单元，用AppendAxisymmetric搞个二维轴对称模型省点计算量。材料属性直接上内置的液态水，密度和声速参数别照搬默认值，实测海域数据记得用Parameter节点喂进去。
model.param.set('rho_water', '1025[kg/m^3]', '声学介质密度'); model.param.set('c_water', '1500[m/s]', '水中声速');边界条件设置是重头戏。入射边界挂上Background Pressure Field，散射场计算记得勾选"Compute scattered field"。重点来了——在超表面边界启用阻抗边界条件，这里藏着相位调控的秘密：
% 设置周期性边界阻抗 surf_imp = model.physics('acpr').feature.create('imp1', 'ImpedanceBoundary', 2); surf_imp.set('Z', '1.5*rho_water*c_water*(1 + 0.3i)'); //复数阻抗控制相位网格划分建议用自由四面体+边界层，声学计算对网格质量敏感得很。频率扫描范围设置在20kHz-50kHz，步长别太抠门，毕竟水下超表面工作带宽要紧。
跑完仿真别急着关窗口，后处理藏着宝贝。反射系数得用场计算器玩点花样：
% 提取反射声压幅值 refl_coeff = mphinterp(model,{'pr'},'coord',[0;0.1],'dataset','dset1')/p0;相位计算更讲究，得找准参考点。建议在入射波路径上取两个监测点做相位差：
phase1 = angle(mphinterp(model,{'p'},'coord',[0,0.5])); phase2 = angle(mphinterp(model,{'p'},'coord',[0,0.6])); delta_phase = unwrap(phase2 - phase1);遇到过反射相位跳变的坑吗？那是计算器自动取主值惹的祸，记得用unwrap函数把相位展开。有个骚操作是直接导出数据到MATLAB做复数运算：
[~,p_data] = mphmean(model,{'p'},'surface','surface1'); reflection_phase = angle(p_data);最后整个参数化扫描，把超表面齿高、周期这些参数设为变量，批量跑完数据导出来做曲面拟合。别光盯着绝对相位值，相对相位梯度才是调控波前的关键。记住，COMSOL的活做完只是开始，拿Python或MATLAB做后处理才能解锁超表面的真正玩法。
