xmlns="http://www.w3.org/2000/svg"

style="display:

none;">

COMSOL，单个金纳米颗粒光热仿真，文章复现，波动光学，固体传热

/>

金纳米颗粒的光热效应挺有意思的——当特定波长的光照射时，电子集体振荡产生局域表面等离子体共振，这时候光能瞬间转化为热能。

咱们在COMSOL里复现这个现象时，得同时操作波动光学和固体传热两个模块。

比如设置金纳米球直径80nm，周围是水环境，用平面波垂直入射。

/>

材料属性设置是个技术活。

金的介电函数得用实验数据拟合的Drude模型：

%
eps_inf
1i*gamma*omega);

这里有个坑：当粒子尺寸小于电子平均自由程时，gamma需要根据尺寸修正。

比如80nm颗粒得把碰撞频率乘以1.5倍，否则吸收峰位置会偏移。

/>

COMSOL，单个金纳米颗粒光热仿真，文章复现，波动光学，固体传热

/>

波动光学模块里，边界条件要设置散射场公式。

有个骚操作是在球体表面添加端口激发，这样比平面波入射更容易收敛。

网格划分建议用四面体+边界层，特别是颗粒表面要做至少3层边界层网格，否则电场梯度算不准。

耦合传热时，热源项得用时间平均的焦耳损耗：

//
=
es.Ez));

这里有个细节：电磁场计算用的是频域，而传热是时域的话需要转换时间尺度。

通常直接做稳态传热计算更省事，毕竟激光脉宽远大于电子弛豫时间。

复现文献结果时常见翻车现场：

1. 温度场比原文小一个量级？检查是否漏了环境对流系数。

   水的自然对流系数大概在1e3-1e4

   W/(m²·K)之间

2. 吸收截面总对不上？试试把计算域扩大到5倍波长，PML层厚度设为半波长
3. 内存爆了？把对称性利用起来——对于球体，其实用二维轴对称模型就能搞定

最后秀个温度分布云图，再对比下不同粒径的共振峰位移，一篇像模像样的仿真复现就搞定了。

不过要当心，实验中的颗粒往往有尺寸分布，仿真结果太完美反而显得假，加点随机扰动更真实。