comsol声子晶体复能带模型 PDE模块 声学黑洞 复能带模型 实能带与虚能带的绘制 参考论文 前两个是论文图后四个是模型及结果图。 可根据模型设置进行其他二维结构的分析复能带这玩意儿搞声子晶体的肯定不陌生但用COMSOL PDE模块手搓模型的时候总有几个坑得自己趟。声学黑洞结构本身自带波导衰减特性把它和二维晶格结合做复能带分析相当于给传统能带曲线加了虚部参数——这时候实部对应频率虚部直接反映能量损耗强度比单纯实能带多了层物理信息。COMSOL里新建PDE模块时记得勾选频域瞬态混合求解别用默认的频域稳态重点在于波动方程要改写成复数形式。比如声压场的控制方程得手动添加虚数项% 弱形式PDE设置COMSOL弱形式模板 epsilon 0.05; % 损耗因子 test(p)*j*omega*(1 - 1i*epsilon)*p test(gradx_p)*c^2*gradx_p test(grady_p)*c^2*grady_p 0;这里的虚数项1i*epsilon就是引入损耗的关键相当于在亥姆霍兹方程里加了个复波数。网格划分建议用三角形边界层特别是声学黑洞区域需要加密到λ/6以下否则高虚部频段容易发散。comsol声子晶体复能带模型 PDE模块 声学黑洞 复能带模型 实能带与虚能带的绘制 参考论文 前两个是论文图后四个是模型及结果图。 可根据模型设置进行其他二维结构的分析数据提取阶段在能带图扫描后导出k点和对应的复频率。用Matlab做后处理时注意频率虚部的符号处理import numpy as np data np.loadtxt(band_data.dat) k_points data[:,0] # 标准化波矢 real_freq data[:,1] # 实部GHz单位 imag_freq data[:,2] # 虚部衰减系数 # 绘制复能带 plt.figure(figsize(10,6)) plt.scatter(k_points, real_freq, cimag_freq, cmapviridis, s15) plt.colorbar(labelImaginary Part (1e8 Np/m)) plt.xlabel(Wave vector Γ-X-M-Γ) plt.ylabel(Frequency (GHz))颜色映射的虚部数值范围建议取对数坐标特别是当存在数量级差异时。有个骚操作是把高损耗区域用半透明红色层叠加在实能带上这样既能看清带隙又突出衰减特性。换其他二维结构比如蜂窝状或Kagome晶格时重点调整晶格基矢和原胞参数。比如改成正方形晶格只需修改周期性边界条件的波矢分量% COMSOL周期性边界设置 k_x k0 * [1, 0]; % 原波矢分量 k_y k0 * [0, 1]; % 改为正方晶格但要注意某些对称性高的结构可能导致虚部能带出现简并这时候需要微调单元内的散射体尺寸打破对称性。实测发现当散射体偏心超过0.12a时虚部能带分裂现象明显。最后提醒复能带计算时间比常规能带多出30%-50%建议先用粗网格试跑观察趋势锁定关键频段后再上精细网格。遇到过虚部结果震荡的情况八成是网格不够密或者PML层吸收边界没设好——把PML的拉伸系数从2调到3.5通常能稳下来。