飞秒激光加工中的双温方程与热力耦合仿真

1. 项目背景与核心问题

飞秒激光加工技术因其超短脉冲特性（10^-15秒量级）在精密制造领域展现出独特优势。当激光作用于金属表面时，电子与晶格子系统会呈现显著的非平衡态，这种瞬态过程用传统热传导模型已无法准确描述。我们的研究聚焦于构建双温方程（Two-Temperature Model, TTM）与热力耦合模型的多物理场仿真框架，揭示飞秒激光烧蚀过程中能量传递与力学响应的耦合机制。

关键挑战：电子-声子弛豫时间尺度（皮秒级）与激光脉冲宽度（飞秒级）的差异导致传统连续介质假设失效，需要建立跨尺度耦合模型。

2. 模型构建与理论基础

2.1 双温方程控制方程组

电子子系统与晶格子系统的能量交换通过以下耦合微分方程描述：

code复制∂Te/∂t = ∇·(ke∇Te)/Ce - G(Te-Tl)/Ce + Q/Ce
∂Tl/∂t = G(Te-Tl)/Cl

其中：

• Te, Tl：电子与晶格温度（K）
• ke：电子热导率（W/m·K），具有Te依赖性
• Ce, Cl：电子与晶格热容（J/m³·K）
• G：电子-声子耦合系数（W/m³·K）
• Q：激光热源项（W/m³），采用高斯时空分布：

code复制Q(r,z,t) = 0.94J/(tpπr0²) × (1-R) × exp[-2(r/r0)² - (z/δ) - 4ln2(t/tp)²]

2.2 热力耦合实现方案

在COMSOL中通过以下模块联动实现耦合：

1. 热传导模块：求解双温方程
2. 固体力学模块：计算热膨胀应力
3. 移动网格模块：跟踪烧蚀界面演化
4. 自定义PDE接口：处理电子热导率非线性项

关键耦合参数设置：

matlab复制% 铜材料参数示例
G = 2.6e16;  % 电子-声子耦合系数
ke_0 = 400;  % 稳态热导率
Ce = γTe;    % 电子热容（γ=96 J/m³K²）
Cl = 3.5e6;  % 晶格热容（恒定值）

3. COMSOL仿真实现细节

3.1 几何建模与网格策略

采用二维轴对称模型简化计算：

• 计算域：20×10 μm（径向×轴向）
• 边界条件：
  • 激光入射面：热流密度边界
  • 其他表面：热对流+辐射复