1. 晶圆表面颗粒物仿真难点到底在哪晶圆表面颗粒物生成与沉积是半导体前道和后道都会遇到的顽固问题。一片12英寸晶圆表面只要有几十颗亚微米级颗粒落在关键图形区后续光刻的对准、薄膜的附着、键合的界面就可能出问题。产线上通常靠清洗、吹扫、洁净室等级控制去压这个指标但真要回答“这些颗粒是从哪来的、为什么会落在这个位置、多大的颗粒最难去除”光靠显微镜数颗粒是不够的。这正是EDEM这类离散元仿真工具能插进来的地方把颗粒当成一颗一颗独立的个体去追踪把晶圆当成几何体把颗粒与晶圆、颗粒与颗粒之间的碰撞、粘附、滚动全部按力学模型算出来。我用EDEM做这个入门案例目的不是复现某条真实产线而是搭一套最小可用的算例框架晶圆摆在那儿上方颗粒工厂源源不断生成颗粒颗粒在重力或简化气流作用下往下走落到晶圆表面后靠粘附力留下来最后统计表面颗粒的数量、面密度和粒径分布。这套流程走通一遍你后面想加静电场、加流场、加晶圆旋转或翻转都是在同一个骨架上加模块。适合谁看刚接触离散元、手里有软件但不知道第一个案例该建什么的同学也适合做洁净控制、设备搬运、晶圆传片机构设计的工程师想粗略估一下“这个动作会不会扬尘、颗粒落下来多少”。先说一句可能让人不太舒服但必须提前说明的话这个案例的物理尺度极其不友好。晶圆直径300 mm污染颗粒粒径常在0.1 μm到5 μm之间差了三到四个数量级而DEM求解的时间步长又必须小于颗粒自身弹性波传播的特征时间。这两件事叠加意味着“按真实尺度老老实实算一秒”在普通工作站上是不可行的。所以入门案例的核心不是把参数填得多准而是先搞明白哪些量可以缩放、哪些量绝对不能乱动。下面我按自己实际搭这套算例的顺序把设计思路、参数来源、操作细节和踩过的坑一层层铺开。2. 案例搭建前的整体设计动手点鼠标之前得先把物理问题翻译成计算模型。这一步做不好后面所有的参数调整都是在错的框架里打补丁。我通常会把一个DEM算例拆成四个问题计算域多大、颗粒代表什么、时间算多长、输出什么指标。这四个问题回答完几何、工厂、求解器设置基本就自动确定了。2.1 建模范围与简化取舍真实场景里晶圆躺在一个片盒、传片手臂或者工艺腔体里周围有气流、有机械振动、有静电。全都要算不动全都不要结果没意义。我的取舍原则是只保留对这个指标有直接贡献的因素。入门案例里我只保留晶圆本体、一个简化包围盒当作计算域边界同时模拟腔体侧壁、颗粒工厂、重力场。气流先不加因为加了就得耦合耦合就得引入流体网格和时间尺度匹配入门阶段容易把自己绕进去。晶圆本体我用直径100 mm的缩比圆盘代替真实的300 mm厚度按同比例缩到0.3 mm左右。为什么缩比因为颗粒数量、接触检测的开销跟几何尺寸和颗粒数量强相关。把直径砍到三分之一同样面密度下颗粒数少了将近九成跑起来轻松很多而颗粒在表面的沉积行为、粘附判据、粒径分布规律不会因为晶圆变小而改变。只要在后处理里把“颗/mm²”这个归一化指标算出来缩比模型和全尺寸模型是可比的。这里必须明确这是基于常见仿真实践做的等效处理不是真实尺寸复现写报告时要把缩放关系写清楚别让人误读成实物数据。包围盒我做成一个开口长方体底面就是晶圆所在平面稍微往下一点四周壁面设成可以吸收颗粒或者直接让颗粒离开计算域避免颗粒堆在角落里反复弹跳浪费时间。颗粒工厂放在晶圆正上方约80 mm处做成一个薄片状的虚拟体积让颗粒从近似一个面上往下释放这样落点分布比较均匀不会全都堆在圆心。2.2 输入参数清单与来源离散元算例最容易被质疑的就是“你这参数哪来的”。我习惯在动手前先列一张表把每个参数、取值、来源和敏感度写清楚跑完再回头看哪些参数其实不敏感可以粗放哪些必须标定。参数入门取值来源与说明敏感度颗粒密度2650 kg/m³二氧化硅类颗粒的常见量级中颗粒半径20~50 μm缩比折算后的等效粒径非真实污染颗粒高颗粒泊松比0.22脆性氧化物类材料的常见范围低颗粒剪切模量3×10¹⁰ Pa为放大时间步有意识降低非真实值高恢复系数颗粒-晶圆0.5常规金属/氧化物碰撞的中间值中静摩擦系数0.3干燥表面常见范围中滚动摩擦系数0.02微颗粒滚动阻力需标定高JKR表面能0.05 J/m²粘附强度的核心可调参数极高重力加速度9.81 m/s²常数低时间步长Rayleigh步长的20%稳定性要求高这张表里有两个值我要特别强调因为它们直接决定了算例能不能跑、结果像不像。第一个是剪切模量。真实二氧化硅的剪切模量在30 GPa量级如果老老实实填进去Rayleigh时间步会小到0.5 ns左右算1秒需要二十亿步根本不现实。常规做法是把剪切模量人为降低一到两个数量级时间步可以放大约3倍而颗粒的宏观堆积、沉降、碰撞行为变化不大——这个技巧在离散元圈子里叫“模量缩放”是被广泛接受的工程近似手段但前提是颗粒的变形量不能大到穿透且碰撞速度不能太高。第二个是JKR表面能它是控制“颗粒落上去粘不粘”的总开关后面我会专门讲怎么标定。2.3 三招让算例跑得动尺度矛盾摆在面前必须做取舍。我的做法是三招叠加而不是单靠某一招硬扛。第一招是粒径放大。把代表污染颗粒的粒径从亚微米放大到20~50 μm。这样一来时间步能放大几十倍颗粒数量也降下来了。代价是失去了真实小颗粒的布朗运动效应和范德华力主导的物理图像所以这个算例只能用来研究“宏观沉降与粘附趋势”不能用来预测真实亚微米颗粒的沉积率。这一点必须在结论里写明白否则就是误导。第二招是时间窗截断。不去算“颗粒从生成到稳态沉积几个小时”只算“释放后0.5~1秒内颗粒的落点分布”。在这个时间窗里颗粒基本完成了自由下落和初次碰撞粘附与否已经能看出来。稳态覆盖率这类指标可以用多个短时间窗的结果做统计外推或者用“一次性释放N颗”的方式代替“持续释放”。第三招是接触检测优化。EDEM的邻居搜索靠网格分箱网格尺寸设得合适能省大量时间。我一般让网格边长略大于最大颗粒直径太小会导致一个颗粒跨多个网格反复检索太大会让每格里的候选颗粒过多。另外把“只在需要时才开启的接触模型”关掉比如入门阶段不需要传热、不需要电荷、不需要破碎这些模块挂着不用也吃内存。三招叠完之后一台普通多核工作站跑几千颗颗粒、一秒物理时间大概几个小时能出结果属于可接受范围。如果你连这个都嫌慢那就把颗粒数再砍一半先看趋势别一上来就追求统计意义上的漂亮曲线。3. 几何建模与颗粒工厂实操几何和颗粒工厂是这个案例里最“手工”的部分也是最容易出问题的地方。软件里点几下就能建个圆柱但要让颗粒乖乖从工厂里出来、按预期落下去、落在晶圆上不弹飞细节全在这儿。3.1 晶圆与容器的几何处理晶圆我用软件自带的圆柱体基元建直径100 mm、高度0.3 mm位置放在原点上表面z0。这里有个细节晶圆的侧面和上下表面要分别命名。EDEM里几何体的每个面section可以单独指定接触参数和运动状态如果整个圆柱只当成一个面后面想单独统计“落在上表面的颗粒”就很麻烦。我在导入或创建几何之后会在几何编辑器里把上表面单独拎出来命名成“Wafer_Top”下表面和侧壁另起名字。这个习惯在后期做接触统计时省事很多。包围盒我用长方体基元拼四面墙加一个底板底板放在晶圆下方2 mm处避免和晶圆共面导致接触判定混乱。四面墙的高度大概150 mm顶部敞开。墙体的接触参数可以设成“颗粒一碰就消失”用颗粒删除或者让颗粒飞出边界后自动移除这样颗粒不会在角落里反复弹跳浪费计算资源。注意几何体之间千万不要有重叠或者共面。共面会导致接触法向计算出现奇异表现为颗粒在某个平面上疯狂抖动或者直接穿过去。我最早搭这个算例的时候底板和晶圆下表面放在同一个z坐标结果颗粒穿模穿得一塌糊涂排查了大半天才发现是共面问题。错开0.5 mm以上的间隙就能解决。几何单位也要确认。软件默认可能是米也可能是毫米导入的CAD模型和手动建的基元如果单位不一致会出现“晶圆大得像足球场”或者“颗粒小得看不见”的情况。我一般统一用米制SI所有尺寸在建模时就换算好避免中途改单位导致参数混乱。3.2 颗粒工厂的生成逻辑颗粒工厂是这个案例的名字来源——“颗粒物生成”。EDEM里的工厂分动态和静态两种静态工厂在仿真开始瞬间一次性把颗粒全部生成完动态工厂则在指定时间内按速率持续生成。研究沉降落点分布我更倾向用动态工厂因为一次性生成几千颗颗粒会在初始时刻造成严重的重叠和爆炸式排斥颗粒初速度乱飞落点完全失真。工厂的核心设置有三项生成体积、生成速率、粒径分布。生成体积我做成一个位于晶圆正上方80 mm处、尺寸120 mm×120 mm×2 mm的薄盒。薄一点的好处是颗粒初始位置集中在一个近似平面内落点分布更接近“面源沉降”而不是“体源扩散”。生成速率设成每秒钟500颗总生成2000颗物理时间上大概4秒生成完——但因为我只算1秒实际是分多次短算例跑每次生成几百颗累加统计。粒径分布用正态分布均值35 μm标准差8 μm并给上下截断20~50 μm。截断很重要正态分布理论上会出现负半径软件虽然会拦但生成失败会在日志里刷一堆警告看着心烦还影响排查真正的错误。初速度我给了向下的0.1 m/s纯粹是为了让颗粒快点离开工厂区域避免在工厂体积内互相碰撞干扰。这个速度远小于重力在1秒内产生的速度增量约9.8 m/s所以对最终落点影响可以忽略属于“给个初速度省时间”的工程习惯。提示工厂生成颗粒时如果位置和已存在的颗粒重叠EDEM会延迟生成或者直接报错。薄盒工厂在小粒径、高生成速率下特别容易出现这个问题。解决办法要么把生成速率降下来要么把工厂体积在厚度方向加大一点给颗粒留出互不干涉的空间。3.3 接触模型与JKR参数设定接触模型选错整个算例就白搭。颗粒与晶圆、颗粒与颗粒之间我用的是Hertz-Mindlin无滑移基础模型叠加JKR粘附模型。基础模型负责法向弹性力、切向摩擦力和滚动阻力JKR负责在颗粒与表面分离时施加一个额外的拉脱力模拟范德华力主导的粘附。JKR模型的核心输入就是表面能。软件里填的这个值物理含义是单位面积上的粘附能实际计算拉脱力时按“拉脱力 1.5π·γ·R*”估算其中R*是等效半径。我可以给你算一笔账让你直观感受这个力有多大取颗粒半径20 μm、表面能0.05 J/m²拉脱力大约是1.5×3.14×0.05×2×10⁻⁵ ≈ 4.7×10⁻⁶ N。而同样一颗半径20 μm、密度2650 kg/m³的颗粒重力大约只有 4.4×10⁻¹⁰ N。也就是说粘附力比重力大了四个数量级。这就是为什么颗粒一旦贴上晶圆靠重力、靠轻轻吹气根本弄不下来也是为什么后道清洗工序要用超声、兆声、刷洗这些强制手段。这个数量级关系直接决定了参数设定的思路只要JKR表面能给到“粘得住”的量级颗粒落地后基本就是粘住不会弹跳。所以恢复系数、摩擦系数在这个算例里的敏感度其实不高真正要调的是表面能。我在初次调试时把表面能设得偏小0.005 J/m²结果颗粒落地后弹几下就滚走了表面一颗都留不住看着像“仿真失败”其实是粘附力不够属于参数问题不是模型问题。颗粒与颗粒之间我也开了JKR因为落下来的颗粒如果聚成团团块的沉降速度和粘附行为跟单颗粒不一样。但如果你的算例颗粒数很多颗粒间粘附会让它们迅速抱团成一个大球反而不利于观察单颗粒沉积。这种情况可以把颗粒间表面能设小一点或者干脆关掉颗粒间粘附只保留颗粒-晶圆的粘附。这是典型的“按研究目的取舍”没有标准答案。4. 求解器设置、运行与监控参数填完点开始之前求解器这几项必须过一遍。很多新手跑出来的结果诡异问题就出在这几个默认值上。4.1 时间步长与计算域网格时间步长我不用软件自动推荐的默认值而是自己算一遍。Rayleigh时间步的估算公式是 t_R πR·√(ρ/G) / (0.1631ν 0.8766)取它的20%作为实际时间步。以半径20 μm、密度2650 kg/m³、剪切模量3×10¹⁰ Pa、泊松比0.22代入Rayleigh步长约在10⁻⁷ s量级20%就是2×10⁻⁸ s。这个数意味着什么算1秒物理时间需要五千万步。即便颗粒数只有两千也是极其漫长的。所以我前面说的“时间窗截断”不是偷懒是必然选择。我的实际做法是把总时长设成0.2秒几百万步跑一晚上能出结果。想看更长时间的沉积就分段跑把上一段结束时的颗粒状态作为下一段的初始状态EDEM支持从上次结果继续或者干脆用不同的随机种子多跑几个短算例做统计平均。计算域网格邻居搜索网格的边长我设成最大颗粒直径的1.2倍左右也就是60 μm。设太小一个颗粒会横跨多个格子接触检索时反复遍历设太大每个格子里塞的候选颗粒太多检索量爆炸。这个值在颗粒粒径分布很宽的时候要特别小心如果最大颗粒是最小颗粒的十倍网格尺寸就很难兼顾通常按最大粒径设接受一点效率损失。保存间隔也要设合理。设得太密结果文件几个G后处理打开都卡设得太疏颗粒落地瞬间的动态过程就丢了。我一般按总时长的百分之一设保存间隔0.2秒的算例每2毫秒存一帧一百帧足够看清沉降过程。4.2 计算过程中的观察点点了开始不等于就可以去干别的。前几分钟必须盯着看因为绝大多数建模错误会在最初几帧就暴露出来。我通常关注三件事。第一颗粒从工厂出来的时候是不是“炸开”的。如果一生成就四散飞射说明生成位置重叠严重颗粒间产生了巨大的排斥力得回去调工厂体积或生成速率。第二颗粒下落过程中有没有穿模。正常的轨迹是一条平滑的抛物线如果某颗粒突然消失或者出现在晶圆下方就是穿模多半是时间步太大或者几何法向有问题。第三颗粒落到晶圆表面后的行为。正常的现象是接触瞬间有轻微反弹然后在一两个时间步内被粘附力“拽住”之后基本静止或者缓慢滚动。如果颗粒一直弹、一直滚说明粘附不够如果颗粒接触瞬间就“定住”且位置诡异可能是接触刚度太大导致数值刚性问题。软件的实时监控窗口可以看颗粒总数、接触数、平均速度这些曲线。我一般会盯“与晶圆接触的颗粒数”这条曲线它应该是阶梯式上升的每一级台阶对应一批颗粒落地。如果曲线上下剧烈跳动说明有颗粒在反复脱离和重新接触粘附参数或者时间步需要复查。5. 后处理把“落在晶圆上”变成可量化数据仿真跑完画面好看没用得拿出能写进报告的数和图。这一步我一般分三层直接读数、空间统计、力学估算。5.1 接触统计与区域划分软件自带的分析模块里几何体的接触数据可以直接导出。我关心的第一个指标是与晶圆上表面接触的颗粒数以及这个数随时间的变化。这条曲线能告诉你沉积速率和饱和趋势。第二个指标是空间分布。我用网格分箱工具把晶圆上表面划分成若干个小格统计每个格子里接触颗粒的数量换算成面密度颗/mm²。网格划分要跟晶圆尺寸匹配100 mm直径的晶圆划10 mm见方的格子大概七八十个有效格太少看不出分布规律太多每格样本太少噪声大。第三个指标是粒径分布。把落在晶圆上的颗粒按粒径分档统计你会发现一个很有意思的现象大颗粒落得快、动能大、容易反弹或滚走小颗粒落得慢、粘附力相对重量更占优势、更容易留下来。最终表面上的粒径分布和释放时的粒径分布并不一致这个“选择性沉积”恰恰是仿真的价值所在靠实验中数颗粒很难数清楚。5.2 用脚本算面密度和径向分布软件导出的颗粒位置数据是CSV字段包括ID、坐标、速度、半径。我用Python做二次统计比在软件里点来点去灵活得多。import pandas as pd import numpy as np # 读取导出的颗粒结果数据 df pd.read_csv(wafer_particles.csv) r_wafer 0.050 # 晶圆半径 50 mm单位与导出数据保持一致 z_tol 0.0005 # 表面判定带0.5 mm z_top 0.0 # 晶圆上表面高度 # 第一步筛出高度上贴近上表面、且水平投影在圆内的颗粒 on_top df[(df[Z] z_top) (df[Z] z_top z_tol)].copy() on_top[r] np.hypot(on_top[X], on_top[Y]) on_top on_top[on_top[r] r_wafer] # 面密度 area np.pi * r_wafer ** 2 density len(on_top) / area print(f表面颗粒数: {len(on_top)}) print(f面密度: {density:.6f} 颗/mm^2) # 第二步按径向分环统计看是不是边缘多、中心少 bins np.linspace(0, r_wafer, 11) on_top[ring] pd.cut(on_top[r], bins, labelsFalse) ring_stat on_top.groupby(ring).agg( count(ID, count), mean_d(Radius, lambda x: 2 * x.mean()) ) print(ring_stat)这段脚本里z_tol的选取是个经验活。取得太小颗粒因为接触变形微微陷入晶圆坐标落在表面下方就被漏掉了取得太大把正在下落的颗粒也算进去沉积数虚高。我的经验是取最大颗粒半径的十分之一左右配合接触数据交叉验证一下数量对不对得上。径向分布那张表往往能看出问题。如果边缘环的颗粒数明显高于中心可能是工厂的初始分布不均也可能是包围盒侧壁反弹把颗粒“推”向了边缘。这种系统性偏差不是物理现象是建模瑕疵得回去改。5.3 粘附力与去除难度的快速估算后处理做到这里其实还可以再往前推一步算一下每颗颗粒受到的粘附力与自身重力的比值。这个比值直接对应“这颗颗粒有多难去掉”。比值越大清除越困难。做法是把拉脱力量级估算和颗粒重量估算对每一档粒径都算一遍做成表。你会发现粒径越小这个比值越夸张粒径增大到几十微米以上重力才开始和粘附力同量级颗粒才可能出现“抖一抖就掉”的情况。这个结论对工艺的意义很直接小颗粒一旦落上几乎不可能靠物理振动去除必须靠化学或者流体力学手段削弱它的附着力。这就是为什么清洗工艺里表面活性剂、兆声波这些手段重要——它们不是“吹掉”颗粒是降低界面能、削弱粘附力。6. 常见问题与排查备忘这套算例我从零搭到能出图中间踩过的坑不少。下面按症状归类给一份可以直接对照的排查清单。6.1 颗粒穿模、飞散、莫名消失穿模是所有DEM算例的头号问题症状是颗粒穿过晶圆或者墙壁出现在计算域外。原因通常有三个时间步太大、接触刚度太高、几何法向错误。排查顺序我建议从时间步开始先把步长砍到原来的四分之一如果穿模消失就是步长问题如果还在穿检查几何体之间有没有重叠或共面都没问题再去看接触刚度的设置是不是超出了材料合理范围。颗粒生成时“炸开”是另一类常见现象本质是初始位置重叠。工厂体积相对于颗粒粒径太小、生成速率太快都会导致新颗粒生在老颗粒体内软件强行把它们推开产生巨大的初始速度。把生成速率砍半或者把工厂在生成方向上的厚度加大通常能解决。6.2 颗粒一颗都不粘或者全粘成一坨这两个极端都是粘附参数的问题但方向相反。一颗不粘要么是JKR表面能设小了要么是恢复系数太大导致颗粒一碰就弹飞。可以先试试把表面能提高一个数量级看趋势如果立刻粘住了就说明量级不对再回来细调。我最初就吃过这个亏按直觉设了个“看起来很小”的表面能结果什么也粘不住后来一算拉脱力和重力的比值发现小了两个数量级。全粘成一坨是颗粒间粘附太强颗粒释放后迅速抱团落下来是一个团块而不是分散颗粒。解决办法是把颗粒间的表面能单独调小甚至置零只保留颗粒-晶圆的粘附。很多新手以为接触模型参数必须两两一致其实完全可以根据研究目的区别对待。6.3 算到一半卡死或者慢得离谱先看是不是时间步突然被拉小了。有些情况软件会根据最大速度自动减小步长比如有颗粒获得了异常大的速度穿模前的征兆步长就会被压得极小计算直接卡住。遇到这种情况去看日志里的最大速度找出异常颗粒的位置基本能定位到几何问题。再看颗粒数是不是超了预期。工厂按速率持续生成如果没设总量上限颗粒会一直累积跑着跑着就上万颗速度自然崩。我习惯在工厂里明确设置总生成量并在求解器里开启“超出计算域自动移除”把飞出去的颗粒清掉。6.4 常见问题速查表现象最可能原因快速验证处理办法颗粒穿过晶圆时间步过大步长减半复跑按Rayleigh步长20%重设颗粒生成即四散生成位置重叠降低生成速率观察加厚工厂体积或降速表面一颗都留不住表面能过小提高一个数量级试按拉脱力/重力比值反推颗粒抱团成大块颗粒间粘附过强关掉颗粒间JKR分离颗粒间与颗粒-壁面参数计算后期急剧变慢颗粒数失控看实时颗粒总数设总量上限与移除边界结果文件打不开保存帧太密看文件体积按总时长1%设保存间隔落点集中在边缘侧壁反弹或分布不均看径向分布表侧壁改吸收边界工厂居中提示表格里每一条都是我在实际调试中真真切切遇到的不是抄来的清单。真跑算例的时候把这几个现象和参数对应的关系记熟能省下大量瞎试的时间。最后说点个人的使用体会。这个算例最大的价值不在结果本身——缩比模型算出来的面密度数字没人会拿去做工艺判定。它的价值在于把“颗粒怎么来、怎么落、怎么粘”这条因果链条在计算机里跑通一遍让你对每个环节的关键参数建立起手感。等以后遇到真的问题比如某台设备搬片时扬尘、某个腔体里颗粒回吸到晶圆上你能迅速判断该建什么模型、该重点调哪个参数、结论的可信边界在哪。入门案例就该这么用不求准求通。
1. 晶圆表面颗粒物仿真难点到底在哪晶圆表面颗粒物生成与沉积是半导体前道和后道都会遇到的顽固问题。一片12英寸晶圆表面只要有几十颗亚微米级颗粒落在关键图形区后续光刻的对准、薄膜的附着、键合的界面就可能出问题。产线上通常靠清洗、吹扫、洁净室等级控制去压这个指标但真要回答“这些颗粒是从哪来的、为什么会落在这个位置、多大的颗粒最难去除”光靠显微镜数颗粒是不够的。这正是EDEM这类离散元仿真工具能插进来的地方把颗粒当成一颗一颗独立的个体去追踪把晶圆当成几何体把颗粒与晶圆、颗粒与颗粒之间的碰撞、粘附、滚动全部按力学模型算出来。我用EDEM做这个入门案例目的不是复现某条真实产线而是搭一套最小可用的算例框架晶圆摆在那儿上方颗粒工厂源源不断生成颗粒颗粒在重力或简化气流作用下往下走落到晶圆表面后靠粘附力留下来最后统计表面颗粒的数量、面密度和粒径分布。这套流程走通一遍你后面想加静电场、加流场、加晶圆旋转或翻转都是在同一个骨架上加模块。适合谁看刚接触离散元、手里有软件但不知道第一个案例该建什么的同学也适合做洁净控制、设备搬运、晶圆传片机构设计的工程师想粗略估一下“这个动作会不会扬尘、颗粒落下来多少”。先说一句可能让人不太舒服但必须提前说明的话这个案例的物理尺度极其不友好。晶圆直径300 mm污染颗粒粒径常在0.1 μm到5 μm之间差了三到四个数量级而DEM求解的时间步长又必须小于颗粒自身弹性波传播的特征时间。这两件事叠加意味着“按真实尺度老老实实算一秒”在普通工作站上是不可行的。所以入门案例的核心不是把参数填得多准而是先搞明白哪些量可以缩放、哪些量绝对不能乱动。下面我按自己实际搭这套算例的顺序把设计思路、参数来源、操作细节和踩过的坑一层层铺开。2. 案例搭建前的整体设计动手点鼠标之前得先把物理问题翻译成计算模型。这一步做不好后面所有的参数调整都是在错的框架里打补丁。我通常会把一个DEM算例拆成四个问题计算域多大、颗粒代表什么、时间算多长、输出什么指标。这四个问题回答完几何、工厂、求解器设置基本就自动确定了。2.1 建模范围与简化取舍真实场景里晶圆躺在一个片盒、传片手臂或者工艺腔体里周围有气流、有机械振动、有静电。全都要算不动全都不要结果没意义。我的取舍原则是只保留对这个指标有直接贡献的因素。入门案例里我只保留晶圆本体、一个简化包围盒当作计算域边界同时模拟腔体侧壁、颗粒工厂、重力场。气流先不加因为加了就得耦合耦合就得引入流体网格和时间尺度匹配入门阶段容易把自己绕进去。晶圆本体我用直径100 mm的缩比圆盘代替真实的300 mm厚度按同比例缩到0.3 mm左右。为什么缩比因为颗粒数量、接触检测的开销跟几何尺寸和颗粒数量强相关。把直径砍到三分之一同样面密度下颗粒数少了将近九成跑起来轻松很多而颗粒在表面的沉积行为、粘附判据、粒径分布规律不会因为晶圆变小而改变。只要在后处理里把“颗/mm²”这个归一化指标算出来缩比模型和全尺寸模型是可比的。这里必须明确这是基于常见仿真实践做的等效处理不是真实尺寸复现写报告时要把缩放关系写清楚别让人误读成实物数据。包围盒我做成一个开口长方体底面就是晶圆所在平面稍微往下一点四周壁面设成可以吸收颗粒或者直接让颗粒离开计算域避免颗粒堆在角落里反复弹跳浪费时间。颗粒工厂放在晶圆正上方约80 mm处做成一个薄片状的虚拟体积让颗粒从近似一个面上往下释放这样落点分布比较均匀不会全都堆在圆心。2.2 输入参数清单与来源离散元算例最容易被质疑的就是“你这参数哪来的”。我习惯在动手前先列一张表把每个参数、取值、来源和敏感度写清楚跑完再回头看哪些参数其实不敏感可以粗放哪些必须标定。参数入门取值来源与说明敏感度颗粒密度2650 kg/m³二氧化硅类颗粒的常见量级中颗粒半径20~50 μm缩比折算后的等效粒径非真实污染颗粒高颗粒泊松比0.22脆性氧化物类材料的常见范围低颗粒剪切模量3×10¹⁰ Pa为放大时间步有意识降低非真实值高恢复系数颗粒-晶圆0.5常规金属/氧化物碰撞的中间值中静摩擦系数0.3干燥表面常见范围中滚动摩擦系数0.02微颗粒滚动阻力需标定高JKR表面能0.05 J/m²粘附强度的核心可调参数极高重力加速度9.81 m/s²常数低时间步长Rayleigh步长的20%稳定性要求高这张表里有两个值我要特别强调因为它们直接决定了算例能不能跑、结果像不像。第一个是剪切模量。真实二氧化硅的剪切模量在30 GPa量级如果老老实实填进去Rayleigh时间步会小到0.5 ns左右算1秒需要二十亿步根本不现实。常规做法是把剪切模量人为降低一到两个数量级时间步可以放大约3倍而颗粒的宏观堆积、沉降、碰撞行为变化不大——这个技巧在离散元圈子里叫“模量缩放”是被广泛接受的工程近似手段但前提是颗粒的变形量不能大到穿透且碰撞速度不能太高。第二个是JKR表面能它是控制“颗粒落上去粘不粘”的总开关后面我会专门讲怎么标定。2.3 三招让算例跑得动尺度矛盾摆在面前必须做取舍。我的做法是三招叠加而不是单靠某一招硬扛。第一招是粒径放大。把代表污染颗粒的粒径从亚微米放大到20~50 μm。这样一来时间步能放大几十倍颗粒数量也降下来了。代价是失去了真实小颗粒的布朗运动效应和范德华力主导的物理图像所以这个算例只能用来研究“宏观沉降与粘附趋势”不能用来预测真实亚微米颗粒的沉积率。这一点必须在结论里写明白否则就是误导。第二招是时间窗截断。不去算“颗粒从生成到稳态沉积几个小时”只算“释放后0.5~1秒内颗粒的落点分布”。在这个时间窗里颗粒基本完成了自由下落和初次碰撞粘附与否已经能看出来。稳态覆盖率这类指标可以用多个短时间窗的结果做统计外推或者用“一次性释放N颗”的方式代替“持续释放”。第三招是接触检测优化。EDEM的邻居搜索靠网格分箱网格尺寸设得合适能省大量时间。我一般让网格边长略大于最大颗粒直径太小会导致一个颗粒跨多个网格反复检索太大会让每格里的候选颗粒过多。另外把“只在需要时才开启的接触模型”关掉比如入门阶段不需要传热、不需要电荷、不需要破碎这些模块挂着不用也吃内存。三招叠完之后一台普通多核工作站跑几千颗颗粒、一秒物理时间大概几个小时能出结果属于可接受范围。如果你连这个都嫌慢那就把颗粒数再砍一半先看趋势别一上来就追求统计意义上的漂亮曲线。3. 几何建模与颗粒工厂实操几何和颗粒工厂是这个案例里最“手工”的部分也是最容易出问题的地方。软件里点几下就能建个圆柱但要让颗粒乖乖从工厂里出来、按预期落下去、落在晶圆上不弹飞细节全在这儿。3.1 晶圆与容器的几何处理晶圆我用软件自带的圆柱体基元建直径100 mm、高度0.3 mm位置放在原点上表面z0。这里有个细节晶圆的侧面和上下表面要分别命名。EDEM里几何体的每个面section可以单独指定接触参数和运动状态如果整个圆柱只当成一个面后面想单独统计“落在上表面的颗粒”就很麻烦。我在导入或创建几何之后会在几何编辑器里把上表面单独拎出来命名成“Wafer_Top”下表面和侧壁另起名字。这个习惯在后期做接触统计时省事很多。包围盒我用长方体基元拼四面墙加一个底板底板放在晶圆下方2 mm处避免和晶圆共面导致接触判定混乱。四面墙的高度大概150 mm顶部敞开。墙体的接触参数可以设成“颗粒一碰就消失”用颗粒删除或者让颗粒飞出边界后自动移除这样颗粒不会在角落里反复弹跳浪费计算资源。注意几何体之间千万不要有重叠或者共面。共面会导致接触法向计算出现奇异表现为颗粒在某个平面上疯狂抖动或者直接穿过去。我最早搭这个算例的时候底板和晶圆下表面放在同一个z坐标结果颗粒穿模穿得一塌糊涂排查了大半天才发现是共面问题。错开0.5 mm以上的间隙就能解决。几何单位也要确认。软件默认可能是米也可能是毫米导入的CAD模型和手动建的基元如果单位不一致会出现“晶圆大得像足球场”或者“颗粒小得看不见”的情况。我一般统一用米制SI所有尺寸在建模时就换算好避免中途改单位导致参数混乱。3.2 颗粒工厂的生成逻辑颗粒工厂是这个案例的名字来源——“颗粒物生成”。EDEM里的工厂分动态和静态两种静态工厂在仿真开始瞬间一次性把颗粒全部生成完动态工厂则在指定时间内按速率持续生成。研究沉降落点分布我更倾向用动态工厂因为一次性生成几千颗颗粒会在初始时刻造成严重的重叠和爆炸式排斥颗粒初速度乱飞落点完全失真。工厂的核心设置有三项生成体积、生成速率、粒径分布。生成体积我做成一个位于晶圆正上方80 mm处、尺寸120 mm×120 mm×2 mm的薄盒。薄一点的好处是颗粒初始位置集中在一个近似平面内落点分布更接近“面源沉降”而不是“体源扩散”。生成速率设成每秒钟500颗总生成2000颗物理时间上大概4秒生成完——但因为我只算1秒实际是分多次短算例跑每次生成几百颗累加统计。粒径分布用正态分布均值35 μm标准差8 μm并给上下截断20~50 μm。截断很重要正态分布理论上会出现负半径软件虽然会拦但生成失败会在日志里刷一堆警告看着心烦还影响排查真正的错误。初速度我给了向下的0.1 m/s纯粹是为了让颗粒快点离开工厂区域避免在工厂体积内互相碰撞干扰。这个速度远小于重力在1秒内产生的速度增量约9.8 m/s所以对最终落点影响可以忽略属于“给个初速度省时间”的工程习惯。提示工厂生成颗粒时如果位置和已存在的颗粒重叠EDEM会延迟生成或者直接报错。薄盒工厂在小粒径、高生成速率下特别容易出现这个问题。解决办法要么把生成速率降下来要么把工厂体积在厚度方向加大一点给颗粒留出互不干涉的空间。3.3 接触模型与JKR参数设定接触模型选错整个算例就白搭。颗粒与晶圆、颗粒与颗粒之间我用的是Hertz-Mindlin无滑移基础模型叠加JKR粘附模型。基础模型负责法向弹性力、切向摩擦力和滚动阻力JKR负责在颗粒与表面分离时施加一个额外的拉脱力模拟范德华力主导的粘附。JKR模型的核心输入就是表面能。软件里填的这个值物理含义是单位面积上的粘附能实际计算拉脱力时按“拉脱力 1.5π·γ·R*”估算其中R*是等效半径。我可以给你算一笔账让你直观感受这个力有多大取颗粒半径20 μm、表面能0.05 J/m²拉脱力大约是1.5×3.14×0.05×2×10⁻⁵ ≈ 4.7×10⁻⁶ N。而同样一颗半径20 μm、密度2650 kg/m³的颗粒重力大约只有 4.4×10⁻¹⁰ N。也就是说粘附力比重力大了四个数量级。这就是为什么颗粒一旦贴上晶圆靠重力、靠轻轻吹气根本弄不下来也是为什么后道清洗工序要用超声、兆声、刷洗这些强制手段。这个数量级关系直接决定了参数设定的思路只要JKR表面能给到“粘得住”的量级颗粒落地后基本就是粘住不会弹跳。所以恢复系数、摩擦系数在这个算例里的敏感度其实不高真正要调的是表面能。我在初次调试时把表面能设得偏小0.005 J/m²结果颗粒落地后弹几下就滚走了表面一颗都留不住看着像“仿真失败”其实是粘附力不够属于参数问题不是模型问题。颗粒与颗粒之间我也开了JKR因为落下来的颗粒如果聚成团团块的沉降速度和粘附行为跟单颗粒不一样。但如果你的算例颗粒数很多颗粒间粘附会让它们迅速抱团成一个大球反而不利于观察单颗粒沉积。这种情况可以把颗粒间表面能设小一点或者干脆关掉颗粒间粘附只保留颗粒-晶圆的粘附。这是典型的“按研究目的取舍”没有标准答案。4. 求解器设置、运行与监控参数填完点开始之前求解器这几项必须过一遍。很多新手跑出来的结果诡异问题就出在这几个默认值上。4.1 时间步长与计算域网格时间步长我不用软件自动推荐的默认值而是自己算一遍。Rayleigh时间步的估算公式是 t_R πR·√(ρ/G) / (0.1631ν 0.8766)取它的20%作为实际时间步。以半径20 μm、密度2650 kg/m³、剪切模量3×10¹⁰ Pa、泊松比0.22代入Rayleigh步长约在10⁻⁷ s量级20%就是2×10⁻⁸ s。这个数意味着什么算1秒物理时间需要五千万步。即便颗粒数只有两千也是极其漫长的。所以我前面说的“时间窗截断”不是偷懒是必然选择。我的实际做法是把总时长设成0.2秒几百万步跑一晚上能出结果。想看更长时间的沉积就分段跑把上一段结束时的颗粒状态作为下一段的初始状态EDEM支持从上次结果继续或者干脆用不同的随机种子多跑几个短算例做统计平均。计算域网格邻居搜索网格的边长我设成最大颗粒直径的1.2倍左右也就是60 μm。设太小一个颗粒会横跨多个格子接触检索时反复遍历设太大每个格子里塞的候选颗粒太多检索量爆炸。这个值在颗粒粒径分布很宽的时候要特别小心如果最大颗粒是最小颗粒的十倍网格尺寸就很难兼顾通常按最大粒径设接受一点效率损失。保存间隔也要设合理。设得太密结果文件几个G后处理打开都卡设得太疏颗粒落地瞬间的动态过程就丢了。我一般按总时长的百分之一设保存间隔0.2秒的算例每2毫秒存一帧一百帧足够看清沉降过程。4.2 计算过程中的观察点点了开始不等于就可以去干别的。前几分钟必须盯着看因为绝大多数建模错误会在最初几帧就暴露出来。我通常关注三件事。第一颗粒从工厂出来的时候是不是“炸开”的。如果一生成就四散飞射说明生成位置重叠严重颗粒间产生了巨大的排斥力得回去调工厂体积或生成速率。第二颗粒下落过程中有没有穿模。正常的轨迹是一条平滑的抛物线如果某颗粒突然消失或者出现在晶圆下方就是穿模多半是时间步太大或者几何法向有问题。第三颗粒落到晶圆表面后的行为。正常的现象是接触瞬间有轻微反弹然后在一两个时间步内被粘附力“拽住”之后基本静止或者缓慢滚动。如果颗粒一直弹、一直滚说明粘附不够如果颗粒接触瞬间就“定住”且位置诡异可能是接触刚度太大导致数值刚性问题。软件的实时监控窗口可以看颗粒总数、接触数、平均速度这些曲线。我一般会盯“与晶圆接触的颗粒数”这条曲线它应该是阶梯式上升的每一级台阶对应一批颗粒落地。如果曲线上下剧烈跳动说明有颗粒在反复脱离和重新接触粘附参数或者时间步需要复查。5. 后处理把“落在晶圆上”变成可量化数据仿真跑完画面好看没用得拿出能写进报告的数和图。这一步我一般分三层直接读数、空间统计、力学估算。5.1 接触统计与区域划分软件自带的分析模块里几何体的接触数据可以直接导出。我关心的第一个指标是与晶圆上表面接触的颗粒数以及这个数随时间的变化。这条曲线能告诉你沉积速率和饱和趋势。第二个指标是空间分布。我用网格分箱工具把晶圆上表面划分成若干个小格统计每个格子里接触颗粒的数量换算成面密度颗/mm²。网格划分要跟晶圆尺寸匹配100 mm直径的晶圆划10 mm见方的格子大概七八十个有效格太少看不出分布规律太多每格样本太少噪声大。第三个指标是粒径分布。把落在晶圆上的颗粒按粒径分档统计你会发现一个很有意思的现象大颗粒落得快、动能大、容易反弹或滚走小颗粒落得慢、粘附力相对重量更占优势、更容易留下来。最终表面上的粒径分布和释放时的粒径分布并不一致这个“选择性沉积”恰恰是仿真的价值所在靠实验中数颗粒很难数清楚。5.2 用脚本算面密度和径向分布软件导出的颗粒位置数据是CSV字段包括ID、坐标、速度、半径。我用Python做二次统计比在软件里点来点去灵活得多。import pandas as pd import numpy as np # 读取导出的颗粒结果数据 df pd.read_csv(wafer_particles.csv) r_wafer 0.050 # 晶圆半径 50 mm单位与导出数据保持一致 z_tol 0.0005 # 表面判定带0.5 mm z_top 0.0 # 晶圆上表面高度 # 第一步筛出高度上贴近上表面、且水平投影在圆内的颗粒 on_top df[(df[Z] z_top) (df[Z] z_top z_tol)].copy() on_top[r] np.hypot(on_top[X], on_top[Y]) on_top on_top[on_top[r] r_wafer] # 面密度 area np.pi * r_wafer ** 2 density len(on_top) / area print(f表面颗粒数: {len(on_top)}) print(f面密度: {density:.6f} 颗/mm^2) # 第二步按径向分环统计看是不是边缘多、中心少 bins np.linspace(0, r_wafer, 11) on_top[ring] pd.cut(on_top[r], bins, labelsFalse) ring_stat on_top.groupby(ring).agg( count(ID, count), mean_d(Radius, lambda x: 2 * x.mean()) ) print(ring_stat)这段脚本里z_tol的选取是个经验活。取得太小颗粒因为接触变形微微陷入晶圆坐标落在表面下方就被漏掉了取得太大把正在下落的颗粒也算进去沉积数虚高。我的经验是取最大颗粒半径的十分之一左右配合接触数据交叉验证一下数量对不对得上。径向分布那张表往往能看出问题。如果边缘环的颗粒数明显高于中心可能是工厂的初始分布不均也可能是包围盒侧壁反弹把颗粒“推”向了边缘。这种系统性偏差不是物理现象是建模瑕疵得回去改。5.3 粘附力与去除难度的快速估算后处理做到这里其实还可以再往前推一步算一下每颗颗粒受到的粘附力与自身重力的比值。这个比值直接对应“这颗颗粒有多难去掉”。比值越大清除越困难。做法是把拉脱力量级估算和颗粒重量估算对每一档粒径都算一遍做成表。你会发现粒径越小这个比值越夸张粒径增大到几十微米以上重力才开始和粘附力同量级颗粒才可能出现“抖一抖就掉”的情况。这个结论对工艺的意义很直接小颗粒一旦落上几乎不可能靠物理振动去除必须靠化学或者流体力学手段削弱它的附着力。这就是为什么清洗工艺里表面活性剂、兆声波这些手段重要——它们不是“吹掉”颗粒是降低界面能、削弱粘附力。6. 常见问题与排查备忘这套算例我从零搭到能出图中间踩过的坑不少。下面按症状归类给一份可以直接对照的排查清单。6.1 颗粒穿模、飞散、莫名消失穿模是所有DEM算例的头号问题症状是颗粒穿过晶圆或者墙壁出现在计算域外。原因通常有三个时间步太大、接触刚度太高、几何法向错误。排查顺序我建议从时间步开始先把步长砍到原来的四分之一如果穿模消失就是步长问题如果还在穿检查几何体之间有没有重叠或共面都没问题再去看接触刚度的设置是不是超出了材料合理范围。颗粒生成时“炸开”是另一类常见现象本质是初始位置重叠。工厂体积相对于颗粒粒径太小、生成速率太快都会导致新颗粒生在老颗粒体内软件强行把它们推开产生巨大的初始速度。把生成速率砍半或者把工厂在生成方向上的厚度加大通常能解决。6.2 颗粒一颗都不粘或者全粘成一坨这两个极端都是粘附参数的问题但方向相反。一颗不粘要么是JKR表面能设小了要么是恢复系数太大导致颗粒一碰就弹飞。可以先试试把表面能提高一个数量级看趋势如果立刻粘住了就说明量级不对再回来细调。我最初就吃过这个亏按直觉设了个“看起来很小”的表面能结果什么也粘不住后来一算拉脱力和重力的比值发现小了两个数量级。全粘成一坨是颗粒间粘附太强颗粒释放后迅速抱团落下来是一个团块而不是分散颗粒。解决办法是把颗粒间的表面能单独调小甚至置零只保留颗粒-晶圆的粘附。很多新手以为接触模型参数必须两两一致其实完全可以根据研究目的区别对待。6.3 算到一半卡死或者慢得离谱先看是不是时间步突然被拉小了。有些情况软件会根据最大速度自动减小步长比如有颗粒获得了异常大的速度穿模前的征兆步长就会被压得极小计算直接卡住。遇到这种情况去看日志里的最大速度找出异常颗粒的位置基本能定位到几何问题。再看颗粒数是不是超了预期。工厂按速率持续生成如果没设总量上限颗粒会一直累积跑着跑着就上万颗速度自然崩。我习惯在工厂里明确设置总生成量并在求解器里开启“超出计算域自动移除”把飞出去的颗粒清掉。6.4 常见问题速查表现象最可能原因快速验证处理办法颗粒穿过晶圆时间步过大步长减半复跑按Rayleigh步长20%重设颗粒生成即四散生成位置重叠降低生成速率观察加厚工厂体积或降速表面一颗都留不住表面能过小提高一个数量级试按拉脱力/重力比值反推颗粒抱团成大块颗粒间粘附过强关掉颗粒间JKR分离颗粒间与颗粒-壁面参数计算后期急剧变慢颗粒数失控看实时颗粒总数设总量上限与移除边界结果文件打不开保存帧太密看文件体积按总时长1%设保存间隔落点集中在边缘侧壁反弹或分布不均看径向分布表侧壁改吸收边界工厂居中提示表格里每一条都是我在实际调试中真真切切遇到的不是抄来的清单。真跑算例的时候把这几个现象和参数对应的关系记熟能省下大量瞎试的时间。最后说点个人的使用体会。这个算例最大的价值不在结果本身——缩比模型算出来的面密度数字没人会拿去做工艺判定。它的价值在于把“颗粒怎么来、怎么落、怎么粘”这条因果链条在计算机里跑通一遍让你对每个环节的关键参数建立起手感。等以后遇到真的问题比如某台设备搬片时扬尘、某个腔体里颗粒回吸到晶圆上你能迅速判断该建什么模型、该重点调哪个参数、结论的可信边界在哪。入门案例就该这么用不求准求通。
