Comsol铌酸锂微盘模式求解。
几何建模这块儿，直接在COMSOL里撸个圆柱体就完事。直径设10微米，厚度0.5微米，注意Z轴方向要和晶体c轴对齐。材料库里的铌酸锂参数得手动调各向异性，别直接用默认值：
model.param.set('d', '10[um]'); // 直径参数化 model.geom.create('disk', 'Cylinder'); model.geom('disk').set('radius', 'd/2'); model.geom('disk').set('height', '500[nm]');这里偷懒用参数d代替具体数值，后面扫参改尺寸方便。网格划分建议用物理场控制，边界层得加到盘边缘，毕竟倏逝波区域对模式影响大。遇到过网格太粗导致模式泄露的情况，建议至少划分五层边界层网格。
边界条件这块儿最容易翻车。上下表面用完美电导体（PEC）近似金属电极，侧面用散射边界条件。注意铌酸锂的介电常数张量方向得和坐标系对齐，特别是当盘面不是主晶轴方向时：
model.physics('es').feature('d1').set('epsilonr', {'xx', 'ny^2', 'nz^2'}); //各向异性设置 model.pholution.create('sc', 'ScatteringBoundary', 'geom1');求解器设置建议用特征频率研究，模式数别设太多，3-5个足够。遇到过算出来Q值飘到1e6的情况，八成是边界条件没设对。记得打开场监视器看模式分布，正常TM模应该有明显的环形场分布。
结果后处理时重点关注有效折射率实部和Q值。用这个代码片段提取基模参数：
neff = mphglobal(model, 'real(ewfd.neff)'); //取实部 Q = 0.5*real(neff)./imag(neff); //品质因子估算 disp(['基模neff=',num2str(neff(1)),' Q=',num2str(Q(1))]);最后说个实战技巧：做参数扫描时别傻等，把研究步骤拆成特征频率计算和频域扫描两步走。遇到过直径扫到15微米时模式跳变的情况，这时候得手动选模式避免数据错乱。对了，温度变化对铌酸锂模式影响贼大，记得把热膨胀系数和热光系数打包进材料属性里。
