COMSOL，单个金纳米颗粒光热仿真，文章复现，波动光学，固体传热
金纳米颗粒的光热效应挺有意思的——当特定波长的光照射时，电子集体振荡产生局域表面等离子体共振，这时候光能瞬间转化为热能。咱们在COMSOL里复现这个现象时，得同时操作波动光学和固体传热两个模块。比如设置金纳米球直径80nm，周围是水环境，用平面波垂直入射。
材料属性设置是个技术活。金的介电函数得用实验数据拟合的Drude模型：
% COMSOL内置材料库的近似表达式 eps_inf = 1.54; omega_p = 1.37e16; % 等离子体频率 gamma = 1.2e14; % 碰撞频率 epsilon = eps_inf - omega_p^2/(omega^2 + 1i*gamma*omega);这里有个坑：当粒子尺寸小于电子平均自由程时，gamma需要根据尺寸修正。比如80nm颗粒得把碰撞频率乘以1.5倍，否则吸收峰位置会偏移。
波动光学模块里，边界条件要设置散射场公式。有个骚操作是在球体表面添加端口激发，这样比平面波入射更容易收敛。网格划分建议用四面体+边界层，特别是颗粒表面要做至少3层边界层网格，否则电场梯度算不准。
耦合传热时，热源项得用时间平均的焦耳损耗：
// 在COMSOL的弱形式方程中插入热源 Q = 0.5*real(dot(es.Jx, es.Ex) + dot(es.Jy, es.Ey) + dot(es.Jz, es.Ez));这里有个细节：电磁场计算用的是频域，而传热是时域的话需要转换时间尺度。通常直接做稳态传热计算更省事，毕竟激光脉宽远大于电子弛豫时间。
复现文献结果时常见翻车现场：
- 温度场比原文小一个量级？检查是否漏了环境对流系数。水的自然对流系数大概在1e3-1e4 W/(m²·K)之间
- 吸收截面总对不上？试试把计算域扩大到5倍波长，PML层厚度设为半波长
- 内存爆了？把对称性利用起来——对于球体，其实用二维轴对称模型就能搞定
最后秀个温度分布云图，再对比下不同粒径的共振峰位移，一篇像模像样的仿真复现就搞定了。不过要当心，实验中的颗粒往往有尺寸分布，仿真结果太完美反而显得假，加点随机扰动更真实。
