comsol三次谐波仿真光学仿真最近在折腾非线性光学仿真的时候第三次被三次谐波生成的问题卡脖子了。COMSOL这玩意儿就像个傲娇的猫主子参数调不对分分钟给你摆烂。今天就跟大伙唠唠怎么用波动方程模块驯服这个磨人的小妖精。先打开电磁波频域模块别手贱选成射频模块别问我怎么知道的。材料库里的石英玻璃参数得手动改非线性极化率这里有个骚操作model.param.set(chi3, 1e-22*(m/V)^2); // 三阶非线性系数 model.component(comp1).material(mat1).propertyGroup(NonlinearOptical).set(chi3, {chi3});这串代码相当于给材料打鸡血注意单位换算能把人逼疯。有次手滑写成m²/V²结果谐波强度直接突破大气层老板还以为我发现了新物理现象...网格划分是重头戏基波波长532nm的情况下三次谐波177nm得用极致变态网格。分享个魔改版自适应剖分设置with model.mesh as mes: mes.create(mesh1, geom1) mes.feature(size).set(custom, on) mes.feature(size).set(hcurve, 5) # 曲率部位加密 mes.feature(size).set(hmin, lambda0/20)别被lambda0骗了这里得用三次谐波的波长做基准。有回偷懒用基波波长结果场分布图糊得像马赛克被组里戏称量子打码技术。comsol三次谐波仿真光学仿真求解器设置才是真正的玄学现场。碰到不收敛别急着摔键盘试试这个魔改版迭代配置model.study(std1).feature(time).set(tlist, range(0,2e-15,100e-15)); model.solver(sol1).feature(tDef).set(maxiter, 50); model.solver(sol1).feature(tDef).set(damp, 0.7);阻尼系数调到0.7能治各种不服就像给数值振荡吃了镇静剂。记得监控场能量曲线要是出现心电图式的波动八成是时间步长该剁手了。后处理阶段整点骚操作用LiveLink联动Matlab提取场分布Ez_THG mphinterp(model,Ez,coord,x_coord,dataset,dset3); imagesc(real(Ez_THG)); colormap(jet(256)); colorbar(Ticks,linspace(-1e7,1e7,5),Limits,[-1e7 1e7]);注意谐波场的实部虚部要分开处理有次把两者直接相加生成的花纹像极了毕加索抽象画还被隔壁组要去做PPT封面...最后说个血泪教训千万别在参数化扫描时同时开六个并行线程否则你会看到CPU温度成功突破沸点主机箱变身暖脚器的奇幻场景。建议用批处理模式分批次跑保命要紧。