中国报告大厅网讯，化学机械抛光(CMP)技术作为半导体制造中实现晶圆全局平坦化的核心工艺，其输控系统的可靠性直接决定芯片制造良品率。隔膜阀凭借介质与驱动机构完全隔离的结构特性，成为超洁净研磨液输送系统的关键控制元件。2026年全球隔膜阀市场规模预计达到6.5亿美元，其中半导体行业EP级隔膜阀市场规模同比增长12.8%，达到6.9亿元。在这一产业背景下，研磨液在隔膜阀内的稳定性问题日益凸显——当研磨液颗粒发生团聚形成大于0.5μm的聚团时，将导致抛光晶圆表面产生缺陷，引发电气故障并显著降低制造良品率。聚团不仅破坏硅片研磨抛光的曲度、平坦度和平行度，还会加速系统过滤器失效及隔膜阀膜片磨损。深入探究隔膜阀内研磨液颗粒的团聚-分散机理，对于优化隔膜阀结构设计、提升研磨液配方稳定性具有重要的工程价值与科学意义。

  一、隔膜阀内研磨液流动特性的CFD-DEM耦合建模方法

  《2025-2030年中国隔膜阀行业市场深度研究及发展前景投资可行性分析报告》指出，研磨液作为固液多相流体系，其流动特性研究需兼顾连续相与离散相的相互作用。采用计算流体动力学(CFD)与离散元方法(DEM)双向耦合框架，液相通过有限体积法求解纳维-斯托克斯方程，颗粒相则通过软球碰撞模型追踪运动轨迹。耦合算法采用CPU-GPU异构并行计算策略，CFD求解器运行于CPU设备，DEM求解器利用GPU设备加速，通过数据交换实现流固耦合。

  在颗粒受力分析层面，建立包含流场作用力、接触作用力及DLVO作用力的综合模型。流场作用力涵盖曳力、压力梯度力、升力、虚拟质量力及湍流扩散力，其中曳力系数根据颗粒形状分别选用Schiller-Naumann模型(球形颗粒)、Marheineke-Wegener模型(椭球颗粒)及Haider-Levenspiel模型(多面体颗粒)。接触作用力采用Hertzian弹簧-阻尼模型计算法向力，Mindlin-Deresiewicz模型计算切向力，并引入JKR粘附力模型表征颗粒表面能效应。DLVO作用力包括范德华吸引力与双电层静电斥力，通过自定义Rocky SDK模块实现代码嵌入，其中范德华力作用距离设定为1-10nm，双电层力通过Debye-Huckel参数与zeta电位计算。

  聚团识别采用配位数过滤算法：当颗粒配位数CRN≥1时判定为接触状态，通过连通性搜索构建聚团拓扑结构。聚团表征参数包括分形维数Df(描述结构紧密程度)、回转半径Rg(表征空间尺度)、球形度φ(反映形状规整性)及团聚率λ(定量团聚强度)。团聚率定义为发生接触团聚的颗粒数目占计算域内颗粒总数的比率，计算公式为λ=NA/ND，其中NA为团聚颗粒数，ND为颗粒总数。

  二、隔膜阀结构特征与研磨液流动死区分析

  以美国Parker公司MV-20-04系列气动隔膜阀为研究对象，建立三维几何模型。该隔膜阀采用机械加工PTFE阀体，最大阀口开度2.38mm，内径3mm，碗式阀座直径15mm，倾角45°。经SpaceClaim抽取流体域后，采用Fluent Meshing生成计算网格，最小网格尺寸300μm以满足网格粒径比≥3的要求。

  研磨液物性参数基于硅片粗磨工艺实际配比确定：去离子水质量分数81.95%、三氧化二铝磨粒10%、螯合剂2%、分散剂2%、pH调节剂4%、消泡剂0.05%。计算得混合相密度1111kg/m³，粘度1.28mPa·s，固相体积分数2.785%。颗粒相选用α-Al₂O₃材料，密度4100kg/m³，杨氏模量434GPa，泊松比0.22，表面能3mJ/m²。颗粒粒径分布60-100μm，符合对数正态分布特征，其中80μm颗粒占比最高达40%，60μm与100μm颗粒各占10%，70μm与90μm颗粒各占20%。

  数值模拟设置质量流量入口2L/min，对应雷诺数11789，湍流模型选用RNG k-ε模型。颗粒相设置质量流量0.222kg/min，入口速度4.587m/s，DEM时间步长1.33×10⁻⁷s，数值软化因子0.0001以加速计算，总模拟时长1s。

  流场分析揭示隔膜阀内存在显著的速度分布不均与流动死区。入口直管段流速均匀约4.6m/s，流经45°弯管时因惯性效应在远离曲率中心侧形成高速区(7m/s)，靠近曲率中心侧出现低速回流。流体进入碗式阀腔后速度骤降至1m/s以下，在阀杆底部冲击后分为上下两股：上股流体在阀腔左侧形成涡流区1，与入口射流共同构成死区A;下股流体沿壁面上升形成涡流区2和3，在圆柱阀杆右侧形成死区B。出口90°弯管处流道收缩，流速回升至7m/s，但靠近曲率中心侧因逆压梯度产生局部空化。

  颗粒运动轨迹与流场高度耦合。在45°弯管远离曲率中心侧，颗粒因惯性力和压力梯度力聚集，局部体积分数达17%，碰撞频率9979Hz。阀腔死区A内颗粒滞留时间最长0.56s，死区B内数密度更高但滞留时间较短。碗式阀座底部因流速低、颗粒补充快，形成体积分数超过40%的高浓度区。出口直管段受二次流影响，局部体积分数仍达15%。

  三、隔膜阀内颗粒形状对团聚行为的影响机制

  在暂不考虑双电层力的条件下，研究四种不同形状颗粒的团聚特性：多面体颗粒A(球形度0.900，10个顶点)、多面体颗粒B(球形度0.958，20个顶点)、椭球颗粒C(球形度0.993，纵横比1.2)及球形颗粒D(球形度1.0)。

  流动稳定后，四种颗粒在隔膜阀内的总数分别维持在32万、30万、27.5万、30万个，团聚颗粒数目分别为10000、6000、5000、4000个，对应团聚率3.1%、2.1%、2.0%、1.4%。数据表明团聚率与颗粒球形度呈显著负相关：球形度越低，接触面积越大，聚团越难在流场剪切作用下破碎分散。

  隔膜阀行业技术特点分析指出，多面体颗粒A、B因面接触特性，接触粘附力与范德华力合力更大，形成的聚团分形维数1.3-1.8，呈线性或网状结构，配位数最高达5，聚团尺寸大且稳定。椭球颗粒C的曳力最大，流动性最差，局部体积分数峰值达50%，但碰撞恢复力大，聚团易分散。球形颗粒D的接触面积最小，碰撞恢复力最大，团聚率最低，聚团几乎均为二元结构。

  粒径分布分析显示，60μm最小颗粒在阀内滞留数占比从初始17.73%升至40%，团聚颗粒中60μm颗粒占比超过32%。小颗粒因惯性小、跟随性差，更易在流场死区滞留并发生频繁碰撞，从而促进团聚。80μm颗粒占比从36.79%降至18%，大颗粒因惯性大、跟随性好，快速通过阀腔不易滞留。

  聚团空间分布呈现显著的区域差异性。直管段剪切强度高，聚团基本为二元结构且数量稀少;45°弯管远离曲率中心侧聚团数目增加;阀腔死区因湍流剪切弱、颗粒聚集时间长，形成大量多元聚团，其中多面体颗粒A、B的聚团数目明显多于椭球颗粒C和球形颗粒D。阀口与90°弯管收缩处应变速率超过3000s⁻¹，聚团流经此处经历强烈剪切破碎，导致出口管道内基本为二元聚团。

  四、隔膜阀内双电层力对团聚的抑制效应

  通过自定义模块实现双电层力模型，研究不同zeta电位(30mV、40mV)对团聚的抑制作用。双电层力计算公式基于Gouy-Chapman-Stern理论，考虑紧密层与扩散层的电荷分布，通过泊松-玻尔兹曼方程求解电势衰减特性。德拜长度设定20nm，相对介电常数80，屏蔽距离200μm。

  受力分析表明，双电层力的量级(10⁻⁶N)显著大于范德华力(10⁻⁹N)与接触粘附力(10⁻⁷N)，在颗粒聚集区域形成有效势垒。