亚波长超声聚焦仿真全解析：从COMSOL建模到HIFU与MEMS应用

做超声聚焦仿真这些年我见过最典型的“翻车现场”就是新建模型、画一个凹球面、给辐射面加一个声压边界然后点求解一气呵成。结果出来一看焦点完全不在设计位置声轴上甚至出现了一串假焦点。这类问题几乎都出在同一个地方——没有认真对待“亚波长”这三个字。所谓亚波长超声聚焦就是让焦点尺寸或者器件特征尺寸小于工作波长或者工作在波长量级附近。换到生物超声里这直接关系到高强度聚焦超声HIFU的疗效精度也关系到MEMS超声换能器这类微型传感器能不能真正把声能送进指定位置。COMSOL Multiphysics是这个领域最常见的仿真工具但它能给出什么结果取决于你如何选物理场、如何建几何、如何划网格、如何设置边界条件。这篇文章我会完整拆解一遍亚波长超声聚焦仿真的思路从模型选型、换能器建模、参数设置到网格划分、求解器选择、后处理再聊聊实际过程中容易踩的坑。内容面向正在做超声换能器、超声传感器、MEMS器件、生物超声应用的朋友也包括那些刚接触COMSOL声学仿真、想快速上手又不想走弯路的新手。1. 亚波长超声聚焦为什么仿真比“莽上去”靠谱1.1 先搞懂亚波长卡在哪条线上很多人对“亚波长”有误解以为只要尺寸小就叫亚波长。声波在介质里传播时波长等于声速除以频率λ c / f。以1 MHz的超声波为例在水中声速约1500 m/s算下来波长就是1.5 mm在人体软组织中声速约1540 m/s波长也差不多在这个量级。如果你用一个直径只有零点几毫米的换能器去辐射1 MHz的声波这就是典型的亚波长器件。亚波长聚焦的含义更严格一些。它不只是说器件尺寸小而是指焦点区域的空间尺度突破了常规衍射极限。传统聚焦透镜的焦点半径会被工作波长限制想做小就得提高频率。但在亚波长聚焦场景里可以在不提高频率的前提下通过特殊的结构设计、相位调制或者近场效应把声能压缩到比波长更小的区域里。仿真在这个问题上几乎是不可替代的。因为亚波长声场在空间上高度局域化实验测量极其麻烦水听器探头本身有尺寸会扰动声场光学观测又需要透明的传播介质和复杂的测振系统。而在COMSOL里你可以得到声场内任意一个点的声压、质点速度、声强甚至能直接切一条线看焦点处-6 dB宽度是多少这在实验里很难一次测准。1.2 生物超声、HIFU、MEMS都绕不开这件事高强度聚焦超声是亚波长仿真最典型的应用场景。HIFU的治疗原理是让换能器把超声波聚焦到体内的病灶区域通过声能转化为热能让组织温度在很短时间内上升到凝固性坏死的程度。治疗成败的评判标准一是焦点处是否达到足够的声压和温升二是焦点之外的健康组织是否受到热损伤。所以焦点尺寸、旁瓣水平、声场的空间分布这些物理量在设计阶段就必须清清楚楚。把仿真目标分解一下就有三条线第一是算准声场拿到焦点位置和声压分布第二是算准吸收结合组织对超声的吸收系数摸清产热分布第三是算准温度场通常会耦合经典的Pennes生物热传导方程。声场一维不够准后面的温度场全都会偏这就是亚波长仿真的压力所在。MEMS超声换能器更是把“亚波长”刻进基因。CMUT和pMUT这类器件的工作频率通常在几MHz到几十MHz空气中波长只有零点几毫米而换能器本身尺寸往往也只有几百微米甚至几十微米。器件几何尺寸和波长处于同一个量级甚至更小换能器表面振动与声场之间的耦合就不能用平面波近似必须做完整的声固耦合仿真才能预测它的远场方向性、灵敏度和声压输出。2. 仿真前必须想清楚的三件事2.1 频域还是瞬态先按你的工作方式选COMSOL声学模块里最常用的是“压力声学频域”接口它假设声场是稳态单频谐波求解结果是某个频率下的空间分布。对多数换能器设计和HIFU方案评估来说频域就足够了。你关心1 MHz连续波激励下的焦点声压频域求解给你一个确定的复声压幅值计算量小收敛也比瞬态容易控制。瞬态求解适合另一类问题脉冲超声、谐波成像、声波在多层介质里的传播过程。脉冲信号包含多个频率成分频域求解需要一个个频率去扫再叠加起来效率反而低。如果你要观察声波遇到组织边界后的反射波到达时间或者要模拟空化气泡在负压相中的行为就必须用瞬态。选择还有个原则只要条件允许先用频域算通一遍模型确认几何、材料、边界条件都没问题再考虑瞬态。我在初学阶段跳过频域直接做瞬态结果一个二维轴对称模型算了整整一夜最后还因为边界反射问题不得不返工得不偿失。2.2 换能器建模等效电路、压电-声耦合到底怎么取舍这是超声仿真里最容易纠结的环节。换能器的驱动信号是电压通过压电效应转化为机械振动再由机械振动向介质辐射声波。你要在COMSOL里复现这个过程有三种层次的选择。最简单、也最常用的做法不给换能器建几何只把辐射面设置成“法向加速度”或“法向速度”边界条件。比如你有一块圆片想在1 MHz下获得0.1 m/s的表面振速直接给那个面一个边界加速度就行。这种做法的好处是省事、鲁棒适合先研究声场传播和聚焦规律的场景。缺点是拿不到换能器本身的电学特性也不知道驱动电压要加多少才能达到这个振速。中等复杂度的做法在换能器位置加一个“压力声学-结构边界”耦合把压电层建模成结构域周围是流体域再给压电层加终端电压或电荷边界条件。COMSOL有专门的“压电效应”多物理场耦合你需要输入压电材料的密度、弹性矩阵、压电常数矩阵、相对介电常数矩阵甚至还要设置极化方向和电极边界。这一层能给你非常接近实际的电-机-声转化效率也能看阻抗曲线。还有一种常被搜索的“换能器等效电路”方法。等效电路是把超声换能器在谐振频率附近的机电行为简化成RLC电路再和声学辐射阻抗联系起来。COMSOL可以通过全局常微分方程或者电路接口把等效电路接入声学模型。“和加载组织后的等效电路”是这方向的升级版换能器压到人体组织上时组织相当于给它增加了一个声负载谐振频率和阻抗曲线会发生明显偏移通过仿真可以提前看到这个偏移量指导匹配层设计。2.3 材料参数差一点焦点就会骗你我见过最多的“仿真看起来没问题但实验差别很大”的案例最后都查到了材料参数不准确上。水的声速取1500 m/s还是1480 m/s在宏观声学里无关痛痒但在亚波长聚焦问题里声速误差直接改变相位分布焦点位置和强度都会变。基础参数表中至少要有这几项传播介质的密度、声速、动态粘度如果考虑损耗、声吸收系数。生物软组织中声速大约1540 m/s密度约1050 kg/m³吸收系数通常随频率近似线性增长比如0.5~1 dB/(cm·MHz)。这意味着10 MHz时组织对声波的吸收是1 MHz的十倍以上焦点区的产热分布差异会非常大。如果做HIFU还要设置吸收系数的方式。COMSOL里可以在声学模块的“吸收”选项中设置幂律衰减系数常见写为α α₀ · f^bb通常约1~2。拿这个参数试算一下焦点幅值比理想无损情况低20%~30%都很常见。别慌这是物理真实反而是高频组织穿透深度受限的直接体现。3. 动手实操从几何搭建到后处理全流程3.1 几何和边界条件二维轴对称省一半算力大多数超声聚焦换能器是旋转对称结构一个同心圆片或者球冠凹面绕中心轴旋转后形状不变。这种问题在COMSOL里优先选择“二维轴对称”空间维度把三维问题压缩成二维的r-z平面计算量直接降一个量级以上。几何上建议建三个域换能器辐射面前的液体传播域、后面接的组织域或流体-组织双层、最外层包裹完美匹配层PML。以典型的1 MHz HIFU为例换能器采用凹球面设计开口直径60 mm曲率半径40 mm那就意味着焦点大约在距换能器表面40 mm的地方。f-number约为0.67属于强聚焦。边界条件的设置按物理场来换能器凹面设为“法向加速度”边界幅值根据输出声功率反推液体域与组织域交界面激活连续边界默认就是声压和法向速度自动连续最外层采用完美匹配层吸收向外的声波模拟无限大传播域。PML层不是随意画一个框就行建议厚度至少一个波长内部不需要设置其他边界条件COMSOL会自动处理。3.2 网格划分的亚波长法则再快也要守住底线网格是亚波长声学仿真里最不能将就的环节。声学波动问题要求网格能分辨波长经验法则是线性单元每波长至少10个单元二阶单元每波长至少5个单元。换句话说1 MHz在水中波长为1.5 mm线性网格最大尺寸不要超过0.15 mm。如果你的焦点区域要求亚波长分辨率那么焦点附近的网格必须更密。我会在焦点周围加一个局部加密区网格尺寸取到λ/20甚至更小比如0.05 mm量级。道理很简单要分析焦点-6 dB宽度是否达到亚波长水平后处理时要在这片区域取多个点做曲线拟合网格太粗会直接抹掉焦点峰值。网格过渡也是经验活。均匀细网格在二维轴对称模型里带来的自由度增长还不算恐怖但三维模型就非常可观了。正确做法是在换能器表面到焦点这条传播路径上用规则网格或边界层网格保证声传播方向上的离散误差一致离开主声束的区域逐步放大网格尺寸。网格尺寸突变会引起数值反射所以放大的比例控制在1.3~1.5倍以内。3.3 求解器设置与后处理别只盯着声压场频域求解在COMSOL里可以直接用默认的PARDISO或MUMPS。如果做频率扫描建议使用“辅助扫描”从0.8 MHz扫到1.2 MHz步长20 kHz。每步求解速度都很快几十个频率点几分钟就能算完。扫描结束后可以用一维绘图组画“焦点处声压-频率”曲线直接读出谐振峰和-3 dB带宽。后处理阶段至少要做三件事第一二维云图看声压场的空间分布确认焦点位置落在几何设计的位置附近声轴上没有第二焦点的异常高值第二过焦点沿径向切一条线读取-6 dB声束宽度这是评价亚波长聚焦是否实现的核心指标第三如果关注热效应把声压场作为热源加载到传热模块。COMSOL里可以在瞬态传热方程中直接插入声源项表达式是有效声压的平方乘以吸收系数再乘以频率相关因子公式可以从“生物热学”相关文档里直接找到。有一类细节容易被忽略后处理里的“声压级”默认参考值通常是20 μPa这是空气声学的标准在水声和生物超声里需要换成1 μPa或者直接看Pa数值。不然你导出的焦点声压级单位都拿不准后续文章和报告里会很麻烦。4. 翻车现场仿真发散、反射污染、焦点漂移4.1 仿真发散第一步不是加密网格COMSOL报“不收敛”或者残差不下降大多数人第一反应是网格不够细立刻全局加密结果计算量翻了几倍问题依旧。以我自己的经验发散首先要检查三件事。第一是材料参数和几何尺度单位是否一致。COMSOL的几何建模里默认单位是米但不少人习惯用毫米画图换能器半径40画成了40 m这就不是仿真问题是数学灾难。第二是边界条件定义是否合理。辐射面的法向加速度幅值如果给得过大局部声压直接到MPa量级线性声学假设已经不成立求解当然容易震荡。第三是PML是否接触了强散射体。PML离换能器太近倏逝波还没衰减就被PML截断会产生虚假反射。如果这些都查过还是发散再考虑网格。网格加密要局部做先在焦点附近加密再看残差有没有明显下降。一把梭全模型加密只是自欺欺人。4.2 边界反射严重干扰焦点声压声波到达计算域边界后如果没有被吸收会反射回焦点区域和入射声场叠加形成驻波图案。这个现象在声压云图上看就是焦点附近出现一圈圈干涉波纹焦点幅值忽高忽低。解决反射问题首选是完美匹配层。但PML用不好反而添乱PML内部不能有几何突变不能有散射体厚度至少一个波长外部边界默认就当作“无反射”。在二维轴对称模型里PML画成圆环包围计算域内外半径差至少一个波长。实测下来PML厚度取1.5个波长焦点声压的误差可以忽略。还有一个补充手段给计算域外层加一个衰减区域也就是人为设置较强的声吸收系数让声波在到达外边界前就被“消耗”掉。这种做法相当于海绵消声水池里的吸声尖劈好处是稳健坏处是它会吸收掉一部分本来会反射回焦点的真实声能所以在生物组织里不太合适更适合纯水介质验证模型时使用。4.3 亚波长焦点“看起来很小”却算不对仿真云图里焦点画得很漂亮但一到定量提取就出问题这是后处理设置引起的。常见原因是数据集的采样分辨率不够。你虽然在计算网格上加密了但后处理云图默认用的还是计算网格插值如果焦点区域网格不够密峰值声压会被平滑掉。这就要回到网格划分阶段确定焦点大概在什么位置可以通过几何声学快速估计在焦点附近画一个很小的矩形或圆形区域单独把网格加密到λ/20。这个局部加密区还能让你在提取-6 dB宽度时拥有足够密度的数据点画出来的曲线才光滑。另一个坑是“相位参考面”不一致。频域结果的声压是复数后处理里直接画“abs(p)”没问题但如果你要看声束聚焦后的相位分布一定要统一参考面。否则相位云图会呈现奇怪的黑白条纹让人误以为是声场干涉。5. 从超声换能器到MEMS微型化之后的仿真变化5.1 CMUT和pMUT两条路线怎么选MEMS超声换能器把超声领域往微型化方向推了一大步。CMUT电容式微机械超声换能器靠静电力驱动本质是一层薄膜和一个背电极中间是真空或空气间隙。它的最大优势是带宽宽同一颗器件可以同时覆盖多个工作频率工艺和CMOS比较兼容。坏处是驱动力与电压平方有关非线性效应明显且耦合系数偏低发射灵敏度需要靠阵列规模补。pMUT压电式微机械超声换能器用的是氮化铝、PZT这类压电薄膜靠逆压电效应驱动。它的驱动效率高、功耗低适合电池供电的便携设备在指纹识别、接近检测、内窥超声里用得很多。缺点是带宽窄工作频率主要由结构谐振决定想要多频工作就得堆叠多个不同尺寸的单元。放在仿真里CMUT通常需要做“静电-结构-声学”三场耦合要把薄膜的预应力和空气间隙的电场都建模出来pMUT更偏向“压电-结构-声学”的经典多物理场问题。两者在COMSOL里都有成熟的多物理场耦合接口但要注意薄膜厚度远小于波长的特性结构振动模态会和声场强烈耦合背衬、封装结构都会影响谐振频率和方向性。5.2 小尺寸换能器的声-固-压电耦合仿真要点MEMS换能器尺寸到了亚毫米甚至几十微米二维轴对称近似经常失效。一个圆形pMUT单元可以用轴对称模型快速估算谐振频率但阵列的方向性计算就必须转三维。三维模型的自由度随尺寸成倍上涨尤其要扫频时跑一次要等很久。经验是先做特征频率分析选出目标模态再在目标频率附近做窄带扫频避免大范围盲扫。仿真里还要特别关注“基底效应”。MEMS换能器不是浮在真空中它贴在一个硅基底或者PCB板上。基底和结构的耦合会产生寄生振动模式进而在辐射声场里叠加杂散分量。建模时在换能器周围加一小圈基底层和封装结构能比理想固定边界更准确地预测实际声场。还有声学阻抗匹配的问题。微型换能器辐射面积小辐射阻抗很低与介质严重失配这决定了器件的带宽和灵敏度上限。仿真中可以直接算出声辐射阻抗随频率的变化曲线再拿它指导匹配层厚度和材料选择。这个设计流程和大型HIFU换能器其实是相通的差别只是尺度。6. 最后分享一点做亚波长仿真的大实话从实际项目经验来看亚波长超声聚焦仿真最忌讳一上来就追求“全耦合、高精度”。正确的推进路线是先做一维或二维轴对称的声场分析确认聚焦行为和设计参数匹配再加换能器模型和生物组织评估实际工作条件下焦点有没有偏移、吸收会不会让声压下降过多最后才考虑做温度场耦合、脉冲瞬态分析、优化扫描这类重活。我每次接到新项目都会强制自己走完“纯流体声场 → 入射边界 → 换能器耦合 → 组织耦合 → 多物理场”这条路。看起来麻烦实际上是省时间的。因为每加一层物理结果发生明显变化时你能立刻定位是哪一个环节造成的而不是等整套模型算完才发现问题藏在某一个边界条件里。另外一个小技巧所有模型参数都用一个全局参数表统一管理声速、密度、工作频率、焦距、换能器半径这些尽量写清楚来源。后续调参数、跑参数化扫描、写报告都会从容很多。超声换能器从设计到验证是一个反复迭代的过程仿真能帮你压缩几轮样品试错但前提是每一个环节都要足够扎实。