简介：本资源是一份面向微电子、MEMS及微纳制造领域科研人员与高年级本科生的二维湿法化学蚀刻仿真教学资料，聚焦层流条件下铜基底腔体蚀刻过程的多物理场耦合建模。内容系统解析蚀刻液对流-扩散传输、铜表面氧化溶解反应动力学及腔体几何形变的动态反馈机制，可支撑工艺参数优化（如流速、浓度、蚀刻时间）与形状保真度分析。资源为单文件PDF文档（914KB），完整涵盖模型定义、控制方程推导（含质量输运方程、层流Navier-Stokes简化形式、移动边界通量条件）、边界设置说明及图示化流程解析，附有典型仿真结果示意图与关键参数取值依据。目前已有95人学习下载，适合开展COMSOL等平台仿真实践、理解湿法蚀刻中流体动力学与电化学反应协同作用的研究者深入学习与复现。
1. 为什么二维层流蚀刻仿真不能只看反应速率？——铜腔体形变才是真实瓶颈
湿法化学蚀刻在微电子产线里不是“倒进去搅一搅”就完事的工艺。哪怕你把CuCl₂浓度、温度、pH全调到教科书最优值，实际蚀刻出的腔体仍可能一边深一边浅、底部钝化、侧壁圆角失控——问题不出在化学反应本身，而出在流体与界面的动态耦合上。本例用二维几何模型直击这个常被忽略的核心：当蚀刻液以层流状态（Re < 2000）掠过暴露铜面时，边界层发展、反应物输运受限、副产物堆积这三者会随腔体深度实时反馈，导致蚀刻前沿不再是均匀推进，而是形成典型的“入口加速-中部滞留-深腔衰减”分布。图2中t=10000 s的CuCl₂浓度云图显示：入口区浓度梯度陡峭，蚀刻速率高；而腔体底部浓度趋近于零，反应几近停滞。这不是模型缺陷，而是真实物理——层流下扩散主导传质，流速越低、腔越深，反应物越难抵达底部。因此，本仿真不只输出“蚀刻了多深”，更揭示几何形变如何反作用于流场，进而重构整个反应动力学边界。适合正在调试MEMS压力传感器蚀刻均匀性、光电子器件沟道陡直度，或需要预判IC掩模版下蚀刻偏差的工艺工程师；也适合COMSOL初学者理解“移动边界+对流-扩散+流固耦合”三物理场协同建模的落地逻辑。
2. 三层物理场耦合建模：从方程到COMSOL节点映射
湿法蚀刻的二维仿真本质是三个强耦合物理过程的数值求解：化学组分输运（对流-扩散）、层流流动（Navier-Stokes）、几何形变（移动边界）。三者并非独立求解再拼接，而是通过共享变量（如表面反应速率r_surface）和共用网格（deformed mesh）实现瞬态耦合。下面拆解每个物理场的数学内核与COMSOL实现关键点，避免常见建模陷阱。
2.1 输运稀释物种（tds）：浓度驱动蚀刻速率的底层逻辑
蚀刻反应速率由局部CuCl₂浓度cCuCl₂决定，遵循线性动力学：r_surface = k·cCuCl₂（式2）。该速率直接作为Deformed Geometry模块的输入，驱动边界移动。但浓度本身由对流-扩散方程控制（式1）：
$$\frac{\partial c}{\partial t} + \mathbf{u}\cdot\nabla c = \nabla\cdot(D\nabla c)$$
其中D=1×10⁻⁹ m²/s为CuCl₂在水中的扩散系数，远小于典型流速下的Péclet数（Pe≈uL/D > 10⁴），说明对流主导输运，但边界层内扩散不可忽略。
注意：在COMSOL中，Transport of Diluted Species接口必须绑定Laminar Flow的velocity field（spf），否则u为零，退化为纯扩散模型，完全无法复现入口加速效应。操作路径：
Transport Properties 1 → Convection → u list → Velocity field (spf)。
2.1.1 边界条件设置：区分“固定浓度”与“反应通量”
- 入口/顶部边界（B1, B3）：设为Concentration，c₀=cCuCl₂_bulk（本例取1 mol/m³）。这是蚀刻液初始浓度，也是反应物供给上限。
- 移动蚀刻面（B4, B6, B13, B15, B16）：必须用Flux边界，而非Concentration。因为此处浓度未知，且通量由反应动力学决定：J₀ = r_surface = k·cCuCl₂。若误设为Concentration=0，模型将强制该面浓度为零，导致反应提前终止。
- 其余边界（B2, B5, B7等）：No Flux（默认），符合物理——无反应、无渗透的绝缘壁面。
# COMSOL Model Builder 中 Flux 边界的关键参数设置（对应 Modeling Instructions 第11页） # 在 Flux 1 节点下： # Species: cCuCl₂ ✓ # Inward Flux (J0,cCuCl₂): r_surface # 注意：r_surface 是用户定义变量，非内置函数 # r_surface 定义见 Variables 节点：r_surface = -k_f * cCuCl₂ # 负号表示消耗该代码行后必须紧接说明：r_surface的负号源于COMSOL约定——正向通量指向域内，而蚀刻是物质消耗，故通量方向向外，需加负号。若遗漏，模型将错误地向铜面注入CuCl₂，导致腔体“愈合”而非蚀刻。
2.2 层流流动（spf）：速度场如何被蚀刻形变重塑
层流方程（式3,4）求解稳态/瞬态速度场u和压力p，但其边界条件直接受蚀刻形变影响：
- 移动壁面（蚀刻面）：应用No Slip条件，即u=0。这意味着流体紧贴铜面静止，所有剪切力在此耗散。
- 入口（B1）：设为Inflow，指定流速（本例未显式给出，由参数文件chemical_etching_parameters.txt定义，典型值0.001 m/s）。
- 出口（B14）：Outflow，允许压力自然释放，避免回流。
- 固定壁面（B3）：Wall，设Translational velocity为Manual，但utr=0，即静止壁。
关键洞察在于：腔体深度增加→流道截面积减小→局部流速升高→边界层变薄→传质增强，但此正反馈在深腔中被扩散限制扼杀。图3箭头图证实：t=10000 s时，腔底流速趋近于零，因流体被“挤”向两侧，中心形成低速死区。
提示：若忽略Deformed Geometry耦合，直接在静态几何上跑spf，所得流场恒定不变，无法捕捉“蚀刻越深、流场越畸变”的真实反馈。必须启用Moving Mesh并绑定spf的Wall Movement。
2.3 变形几何（dg）：用网格位移量化铜的“溶解”
Deformed Geometry模块不求解物理，而是根据表面法向速度υ计算网格节点位移，使计算域随蚀刻实时变形。其核心方程（式5）为：
$$\mathbf{v} = \frac{D M}{\rho} \mathbf{n}$$
其中D为扩散系数，M为摩尔质量（Cu=63.5 g/mol），ρ为铜密度（8960 kg/m³）。但本例采用更直接的工程简化：υ = v_surface = -r_surface·M_Cu/ρ_Cu（见Variables定义表）。
该简化跳过扩散项，将蚀刻视为表面反应控制，符合“反应速率限制型蚀刻”前提（原文明确“fast, usually reaction rate limited”）。
# Variables 节点中 v_surface 的完整定义（Modeling Instructions 第9页） # Name: v_surface # Expression: -r_surface * M_Cu / rho_Cu # Unit: m/s # Description: Surface normal velocity (negative for etching) # # Name: r_surface # Expression: -k_f * cCuCl2 # Unit: mol/(m²·s) # Description: Etching flux (negative for consumption)参数k_f（蚀刻速率常数）是标定关键。若k_f过大，腔体瞬间蚀穿；过小则形变不明显。本例k_f=1×10⁻⁶ m/s（隐含于parameters.txt），对应典型Cu/CuCl₂体系在室温下的实测值。务必检查k_f单位是否为m/s——若误用cm/s，结果将放大100倍。
3. 几何构建与参数驱动：从草图到可复现实验
二维模型的几何精度直接决定仿真可信度。本例采用参数化建模，所有尺寸由全局参数控制，便于批量研究不同掩模宽度、初始腔深的影响。几何构建非简单画矩形，而是通过边界选择策略精准定义蚀刻面，这是后续移动边界条件生效的前提。
3.1 参数化几何：四步构建带掩模的铜基底
模型包含三个矩形（r1掩模、r2初始腔、r3流体域）和一个圆角（fil1），全部尺寸由参数文件驱动：
| 参数名 | 数值 | 单位 | 物理意义 |
|---|---|---|---|
| d_mask | 1e-3 | m | 掩模宽度（覆盖铜区） |
| h_mask | 0.5e-3 | m | 掩模高度（y方向厚度） |
| h_boundary_layer | 0.2e-3 | m | 入口边界层厚度（影响初始流场） |
| h_seed_cavity | 0.5e-3 | m | 初始腔半径（图1中圆形凹陷） |
构建步骤严格按顺序执行，尤其注意坐标系原点与选择逻辑：
- 掩模矩形（r1）：宽=d_mask，高=h_mask，位置默认(0,0)，覆盖铜基底上方。
- 初始腔矩形（r2）：宽=高=h_seed_cavity，但x=-h_seed_cavity, y=-h_seed_cavity——将矩形左下角置于(-0.5mm,-0.5mm)，使其右上角恰在原点，形成“裸露铜面”的起始区域。
- 圆角（fil1）：仅选r2的Point 1（左下角）和Point 2（右下角），半径=h_seed_cavity/2。此举将矩形底边圆化，模拟蚀刻起始时的平滑过渡，避免尖角奇点导致网格失真。
- 流体域矩形（r3）：宽=2×d_mask，高=3×h_mask，x=-d_mask/2, y=h_mask——确保完全覆盖掩模及上方流体区，为层流提供充足发展长度。
警告：若在Step 3中误选r2所有4个顶点，圆角将作用于整个矩形，破坏“仅底部暴露”的物理设定。必须严格按文档选Points 1&2（即底边两端）。
3.2 移动边界定义：五条边的精确捕获
蚀刻仅发生在未被掩模覆盖的铜表面，即r2矩形的底边及两侧边（因圆角存在，实际为曲线）。COMSOL中需手动选择这些边界：
- Boundary 4：r2底边（圆角连接处）
- Boundary 6：r2右侧边（从圆角终点垂直向上）
- Boundary 13, 15, 16：r2与r3交界处的流体域边界，共同构成蚀刻面外延。
此选择逻辑源于几何布尔运算后的边界编号规则。若直接用“Select by Expression”输入y<0，会误选r1下方的无效边界。唯一可靠方法是：
- 点击
Zoom Extents全览几何； - 按住Ctrl逐个点击上述5条边；
- 在Explicit节点的Selection列表中确认编号无误。
-- 验证移动边界选择的SQL式检查（COMSOL不支持SQL，但逻辑等效） SELECT boundary_id, geometry_entity, description FROM model_boundaries WHERE boundary_id IN (4,6,13,15,16) AND geometry_entity = 'r2' OR geometry_entity = 'r3' AND description LIKE '%etch%';该伪代码强调：边界必须属于r2或r3，且描述含“etch”关键词。实际操作中，右键Boundary查看属性，确认其“Entity”列显示r2或r3。
3.3 材料与参数加载：避免硬编码的工程实践
材料直接选用COMSOL内置Water, liquid，其密度ρ=998 kg/m³、粘度μ=0.001 Pa·s已预置，无需手动输入。但关键参数k_f,D,cCuCl₂_bulk等必须从外部文件加载：
Parameters 1 → Load from File → chemical_etching_parameters.txt- 该文件内容为纯文本键值对：
k_f 1e-6 [m/s] D 1e-9 [m^2/s] cCuCl2_bulk 1 [mol/m^3] M_Cu 0.0635 [kg/mol] rho_Cu 8960 [kg/m^3]
提示：若手动在Parameters界面输入数值，易遗漏单位或小数点（如k_f=1e-6误为1e-5），导致蚀刻速率偏差10倍。参数文件确保可追溯、可版本管理。
4. 求解器配置与结果验证：识别收敛失败的三大信号
本模型为强非线性瞬态问题（移动边界+浓度依赖速率），求解器配置不当将导致不收敛或虚假结果。必须放弃默认设置，针对性调整。
4.1 时间步长与求解序列：先稳态后瞬态的必要性
Modeling Instructions要求先运行Stationary研究（Step 1），再添加Time Dependent。这是关键技巧：
- Stationary求解：固定几何，求解初始流场u和浓度场c，为瞬态提供高质量初值。若跳过此步，Time Dependent从零场开始，前几个时间步因剧烈变化而发散。
- Time Dependent设置：
Output times: range(0,500,10000)表示从0到10000 s，每500 s输出一次。但实际求解步长由自动时序控制，用户需监控：Relative tolerance: 1e-3（默认1e-2太松，浓度梯度处易失真）Maximum step size: 10（秒）——防止单步跨越关键形变阶段
# Study 1 → Time Dependent → Settings → Solver Configurations # Time Stepping: Intermediate # Relative Tolerance: 0.001 # Maximum Step: 10 [s] # Initial Step: 0.1 [s] # 强制首步精细解析起始反应4.2 收敛诊断：三类失败信号及修复方案
| 信号 | 原因 | 修复动作 |
|---|---|---|
| Error: Failed to find a solution | 移动边界导致网格畸变，雅可比矩阵奇异 | 在Deformed Geometry → Free Deformation → Mesh smoothing type选Hyperelastic（已配置），并增加Mesh quality improvement迭代次数至5 |
| Warning: Maximum number of Newton iterations reached | 浓度场剧烈变化，非线性求解器不收敛 | 在Stationary研究中，Fully Coupled求解器 →Nonlinear controller→Damping factor设为0.7（默认1.0易振荡） |
| Result: Cavity depth oscillates wildly | 时间步长过大，未能捕捉蚀刻前沿突变 | 缩小Maximum step size至1 s，并启用Events检测蚀刻面速度突变，触发自适应步长 |
4.3 结果验证：用三张图交叉印证物理一致性
单一结果图易误导，必须三图联动验证：
- 图2（浓度）：入口区高浓度→高蚀刻速率；腔底浓度趋零→速率衰减。若腔底仍有显著浓度，说明扩散系数D过大或流速过低。
- 图3（速度）：箭头长度对数缩放，清晰显示腔底流速≈0。若箭头均匀分布，说明未启用Moving Mesh，流场未响应形变。
- 图4（网格）：观察r2区域网格是否随蚀刻向下拉伸。若网格保持刚性，证明Prescribed Normal Mesh Velocity未正确绑定v_surface。
验证操作：在Mesh Plot中，右键→Show Element Quality，最小雅可比值应>0.3。若<0.1，需在Free Deformation节点增加Mesh smoothing iterations。
5. 形变-流场-浓度耦合的进阶分析：提取蚀刻均匀性量化指标
单纯观察t=10000 s的最终形貌不够，工艺优化需量化“均匀性”。本节提供三个可直接在COMSOL中提取的指标，用于评估不同流速、浓度下的蚀刻质量。
5.1 腔体轮廓陡直度（Sidewall Angle）
定义为蚀刻面法向与垂直方向夹角。在Arrow Surface图中，右键→Add Evaluation→Line Integration，沿左侧蚀刻面（Boundary 6）积分：
Expression: acos(abs(ny))（ny为y方向法向分量）Unit: deg- 物理意义：理想各向同性蚀刻应为0°（垂直），>5°表明侧壁倾斜，源于入口流速过高导致侧向冲刷。
5.2 底部浓度均匀性指数（Bottom Uniformity Index）
在腔底（y=-h_seed_cavity处）创建Cut Plane，计算浓度标准差σ与均值μ之比：
sigma(cCuCl2)/avg(cCuCl2)- 阈值：若>0.3，表明底部传质严重不均，需增大流速或减小腔深。
5.3 蚀刻前沿速度分布（Front Velocity Profile）
利用Derived Values→Surface Maximum，在移动边界上求v_surface最大值，再用Surface Average求均值：
max(v_surface)vsavg(v_surface)- 解读：比值>3表示蚀刻前沿高度不均，入口区速率远超底部，是深腔蚀刻失效的早期预警。
技巧：将上述三个指标导出为CSV，用Python绘制
流速 vs 均匀性指数曲线，可快速定位工艺窗口。例如，当入口流速从0.0005 m/s增至0.002 m/s，陡直度从2°升至8°，而底部均匀性指数从0.25恶化至0.45——证明过高的流速牺牲了形貌保真度。
本文还有配套的精品资源，点击获取
数据的实战指南)