脉冲激光打靶看着就是"啪"的一下，可真正发生过的事情集中在纳秒到微秒这个窗口里：光子被自由电子吸收、晶格升温、熔融、气化、原子电离、等离子体长大、激光被屏蔽、羽流膨胀、离子复合。实验能测的往往只有终态——坑多深、谱线多强、羽流跑多快，中间的电子密度和电子温度曲线探针很难伸进去。这就是COMSOL 脉冲激光诱导等离子体仿真模型值得花时间搭的原因：它把看不见的中间过程摊开给你看。这篇东西写给两类人：一类是刚接手激光加工或激光光谱课题、需要交出一个能跑的模型的人；另一类是已经把模型跑起来、但结果"看着像那么回事又不敢信"的人。下面我按自己搭这类模型的顺序讲，从物理时间轴、接口选型、参数换算、网格时间步，一直到踩过的坑。

1. 动手之前先把物理时间轴捋清楚

1.1 从光子吸收到等离子体屏蔽的五个阶段

搭模型最容易犯的错，是一上来就打开软件堆物理场。我一般会先在纸上画一条时间轴，把脉冲激光打到固体靶上的过程切成五段，每一段对应不同的主导物理，这样后面选接口才有依据。

第一段是光吸收与自由电子加热。金属靶里光子主要被自由电子吸收，穿透深度通常只有十几到几十纳米，脉宽 10 ns 的话，表层电子温度可以在很短时间内冲到远高于晶格温度，这就是所谓双温效应的来源。第二段是晶格升温与熔融，能量从电子传给晶格，热扩散长度只有亚微米到微米量级，熔池在几十纳秒内形成。第三段是气化与电离，表面温度越过沸点后开始剧烈蒸发，蒸发的原子继续吸收后续激光，被多光子或碰撞电离，电子密度爬升。第四段是等离子体屏蔽，电子密度一旦接近临界密度，激光就进不去了，能量沉积位置从靶面转移到等离子体内部，形成吸收波。第五段是膨胀、冷却与复合，脉冲结束后等离子体绝热膨胀，电子温度掉下来，谱线由强转弱。

这五段在模型里不一定都要算全。真正决定你搭什么样模型的，是第二小节那个问题：你到底想回答什么。

1.2 建模目标的三种定位，决定你该省掉哪些物理

我见过太多人一上手就想"全都要"，结果模型跑到一半发散，最后连温度场都没拿到。实际项目里我会先把目标分成三档，档位不同，砍掉的物理完全不同。

定性看趋势这一档，目标是比较两种波长哪个烧蚀深、或者光斑大小对熔池形状的影响。这种需求下，我通常只保留传热加一个等效热源，等离子体用一层等效吸收区代替，网格可以放到微米级，一个算例几分钟到几十分钟。半定量对比这一档，需要给出熔池深度、羽流速度、电子密度的量级，这时候必须加分相变、流体和简化等离子体输运，网格到亚微米，时间步到皮秒量级。定量标定这一档最费劲，要求电子温度和密度的绝对值能和发射谱线反演结果对得上，需要引入重粒子输运、详细化学动力学、甚至双温模型，算例按天计。

提示：目标档位最好在项目开始时就写进文档，后面每加一个物理场都回头看一眼——它属于哪一档？不属于就删掉。这条规矩救过我不止一次。

1.3 一个常被忽略的前置动作：量纲与量级估算

软件打开之前，我习惯先做十分钟量级估算，这不是走形式，它能直接告诉你网格要划多细、时间步能放多大，也能提前判断某些物理到底该不该算。

举个真实用过的例子：波长 1064 nm、脉宽 10 ns、单脉冲能量 10 mJ、光斑 1/e² 直径 200 μm。先算峰值能量密度，面积是 π×(1e-4)² ≈ 3.14e-8 m²，除以面积得到 3.2e5 J/m²，也就是 32 J/cm² 左右。铁的纳秒脉冲烧蚀阈值大致在 1 J/cm² 这个量级，也就是说这个条件下远超阈值，烧蚀一定发生，模型里必须包含相变和气化。再算热扩散长度，铁的热扩散系数约 2.3e-5 m²/s，用 √(2Dτ) 估算得到 0.67 μm，意味着温度场在亚微米尺度上变化剧烈，网格在靶面法向必须做到百纳米以下。

再看等离子体侧，如果假设电子密度到 1e20 m⁻³、电子温度 1 eV，德拜长度只有 0.74 μm；如果电子密度被推到 1e21 cm⁻³ 量级，德拜长度直接掉到 0.2 nm。这个估算的意义很直接：等离子体鞘层在常规网格下根本解析不了，所以模型里必须走准中性假设，或者对鞘层做特殊处理。不做这一步估算，后面网格怎么划都是瞎猜。

2. 模型骨架：维度、接口与耦合关系的取舍

2.1 维度选择：1D、2D 轴对称还是 3D

维度不是越简单越好，也不是越复杂越准，它取决于你的物理场有没有天然的对称性。激光光斑是圆形的，材料是各向同性的平板，这种情形下 2D 轴对称几乎总是首选——它能保留径向的径向膨胀、径向温度梯度和羽流的侧向扩展，计算量又只有三维的几十分之一。

1D 什么时候用？只关心靶面法向温度历史、或者只想快速试几个脉宽参数的时候。1D 模型的迭代速度极快，我经常用它先扫一遍能量密度范围，找出感兴趣的工作点，再拿结果去定 2D 模型的边界。缺点是它完全丢失了光斑边缘效应和羽流的侧向流动，光斑小到接近热扩散长度时误差会很大。

3D 我只在两种情况下用：一是光斑是长条形或者有整形，二是靶材有各向异性或结构图案。三维模型的内存和时间开销非常现实，一个 200 μm 见方、网格到 100 nm 的三维域，单元数轻松上百万，再叠加流体和等离子体，普通工作站很难在合理时间内跑完。真要上三维，我一般先把等离子体退化成等效吸收层，等温度场和流场稳定了再局部加密。

2.2 物理场接口的组合方案

接口选择的核心逻辑是：谁主导能量守恒，谁就必须被显式求解。下面这张表是我常用的三套组合，从轻到重，按需取用。

轻量方案的搭法是：传热接口里加"沉积光束功率"域条件，直接给高斯光束的功率、脉宽波形、光斑半径和吸收系数，软件内部按 Beer-Lambert 沿光束方向衰减。这套方案的好处是不用划分光波网格，缺点是无法描述等离子体屏蔽——能量沉积位置永远钉在靶面。

中等方案要在传热基础上加层流或者两相流，靶面用动网格跟踪烧蚀后退，相变用表观热容法处理。这里有个关键取舍：熔池里的马兰戈尼对流要不要加？加了更真实，但对流项会让求解难度陡增，我一般先不加，跑通之后再作为敏感性项打开。

完整方案的等离子体部分，我倾向用漂移-扩散加局部能量近似，而不是直接上全动理学。全动理学（比如双项近似玻尔兹曼）能算电子能量分布函数，但和流体耦合起来代价极高；在电子密度 1e20~1e22 m⁻³ 这个区间，麦克斯韦分布近似通常够用，除非你专门研究非平衡效应。

2.3 激光热源的两种处理方式

这是整个模型里最容易出分歧的地方，我把它拆成两条路线讲清楚。

路线一：Beer-Lambert 体热源。把激光当成一束沿深度指数衰减的能量流，热源项写成 Q = (1-R)·I(t)·α·exp(-αz)，其中 R 是反射率、α 是吸收系数、I(t) 是时间波形。这条路线的物理前提是介质在光传播方向上足够厚、侧向散射可忽略。优点是计算便宜、稳定、容易收敛。缺点也很明显：等离子体形成后吸收系数会随电子密度剧烈变化，需要用等效吸收模型去近似，否则会低估屏蔽。

路线二：电磁波频域加吸收介质。直接求解麦克斯韦方程，用复折射率描述材料吸收。这条路线的物理更完整，能处理多层膜、周期性结构、甚至等离激元效应。但代价是网格必须解析波长：1064 nm 在真空中的周期是 3.5 fs，按每波长 10 到 20 个单元算，网格尺度要到几十纳米，再乘上微米级的计算域，单元数会失控。

我的实际做法是混合：传热和流体用 Beer-Lambert 热源，只有在需要精确算吸收分布或者研究表面微结构的时候，才在局部子模型里用全波。两者不矛盾，子模型算出的吸收分布可以作为父模型热源的输入。

2.4 高温物性的补全

常温物性参数在激光加工这个场景下几乎全部失效。铁在 1 eV 下的电导率、热导率、吸收率跟室温差好几倍，等离子体的输运系数更是随电子温度变化几个数量级。

我的处理顺序是：先用常数物性跑通，确认能量守恒和量级正确；然后把热导率、比热、密度换成随温度变化的分段线性插值，插值区间要覆盖到气化温度以上；最后给等离子体区的电导率和热导率加温度依赖。这里有个细节，分段线性插值的外推行为很难控制，温度一旦超出表格上限，软件可能按最后一段斜率一路外推，出现物理上荒谬的值。稳妥做法是给表格两端加平台段，并且在求解器里设置温度上限监控。

另外，气化温度以上的物性数据往往缺失严重。我的经验是：宁可让参数在超高温区保持常数并记录下来，也不要随便找一组数据硬插——错误的数据比没有数据更危险，因为它不会报错，只会给你一个看起来很合理的错误结果。

3. 参数计算：把实验室条件翻译成边界条件

3.1 能量密度、功率密度与峰值功率的换算

这一节看着基础，实际出错率极高。实验同门口头说的是"能量密度 30 J/cm²"，这个数到底是峰值还是平均值、光斑半径按 1/e² 还是 1/e 定义，不同人的习惯完全不同，直接抄进模型里偏差可能达到两倍。

高斯光束的峰值能量密度和总能量的关系是 F_peak = 2E/(πw²)，其中 w 是 1/e² 半径。如果对方给的是 1/e 半径，需要先换算。下面这套参数定义我一般直接写进模型参数表里，后面所有表达式都引用它，避免混用：

E0 = 10[mJ] // 单脉冲能量 tau = 10[ns] // 脉宽（FWHM 需额外换算） w0 = 100[um] // 1/e^2 光斑半径 Fp = 2*E0/(pi*w0^2) // 峰值能量密度 Ip = Fp/tau // 峰值功率密度

按这组数算下来，Fp ≈ 6.4e5 J/m²（64 J/cm²），Ip ≈ 6.4e13 W/m²。注意这里脉宽用的是方波等效，如果实验给的是高斯时间波形的 FWHM，峰值功率密度还要再乘一个约 1.06 的系数。这个差别在阈值附近会直接影响"烧不烧得动"的判断，别忽略。

还有一点：脉冲时间波形在模型里不能用理想方波，阶跃边沿会让求解器在开关瞬间疯狂缩步长。我用的是平滑阶跃或者分段三次多项式，上升沿给到脉宽的 5% 到 10%，既保留物理又不折腾求解器。

3.2 等离子体特征尺度：德拜长度、临界密度、趋肤深度

这三个量决定了你的模型能不能"解析"等离子体，必须提前算。

德拜长度按 λ_D ≈ 0.74 × √(T_e[eV] × 1e14 / n_e[cm⁻³]) μm 估算。电子密度 1e16 cm⁻³、电子温度 1 eV 时约 74 nm；密度到 1e18 cm⁻³ 时只剩 7.4 nm。这意味着除非你把网格划到纳米级，否则德拜鞘层永远在网格之下，只能靠准中性假设或者鞘层边界条件来处理。

临界密度按 n_c ≈ 1.1e27 × (λ/1μm)⁻² m⁻³ 算，1064 nm 对应约 1e27 m⁻³，也就是 1e21 cm⁻³。当电子密度接近这个值时，等离子体频率接近激光频率，激光无法继续传播，反射和吸收同时剧增。模型里如果电子密度算出来超过临界密度好几倍，而热源还在往深处灌能量，那基本可以判定屏蔽项没做对。

趋肤深度在过密等离子体里约为 δ = c/ω_pe，电子密度 1e27 m⁻³ 时对应几十纳米量级。这个尺度同样小于常规网格，所以实际建模时不追求解析，而是用等效吸收系数把能量沉积"抹"在几个单元的厚度里，保证总吸收能量对得上就行。

3.3 网格与时间步的硬约束

网格划分我遵循"三个尺度取最小"的原则：光学吸收尺度、热扩散尺度、流场边界层尺度。以 1064 nm 打铁为例，光学穿透深度十几纳米，热扩散长度 0.67 μm，流场边界层在微米量级——真正卡脖子的是光学吸收，但全域按十几纳米划网格不现实。

解决办法是边界层网格加梯度自适应：靶面法向第一层厚度取 20 到 50 nm，增长率 1.2，铺 15 到 20 层，把吸收和温度梯度都罩住；径向按光斑半径的 1/50 划，保证高斯分布的形状不失真；远场区域用自由四面体或映射网格快速放大。这样 2D 轴对称域（200 μm × 400 μm）单元数能控制在 3 万到 8 万，是能接受的范围。

时间步的约束来自两处。传热侧，显式格式受限于 Δt < Δx²/(2D)，取 Δx = 50 nm、D = 2.3e-5 m²/s，得到 Δt < 5.4e-11 s，也就是 54 ps；隐式格式没有这个限制，但为了捕捉脉宽内的瞬态，我一般取脉宽的 1/500 到 1/1000，10 ns 脉宽对应 10 到 20 ps。流体侧受 CFL 条件约束，Δt < Δx/u，声速 1000 m/s、网格 0.2 μm 时 Δt < 0.2 ns，比传热宽松。真正紧的是等离子体侧，介电弛豫时间和电子碰撞频率决定了漂移-扩散方程的刚性，求解器会自动缩步长，这时候与其硬扛，不如在物理上做准稳态近似。

注意：见过不少人为了"精度"把时间步设成 1 fs，结果一个 10 ns 的算例要跑一千万步，实际上完全没必要。脉宽内的温度演化特征时间是纳秒，皮秒步长已经足够。

3.4 烧蚀阈值与相变潜热的处理

烧蚀在模型里有三种处理层次，选哪个取决于你要不要看坑的形状。

最简的是焓-孔隙率法：把固液共存区当成一个糊状区，用表观热容把潜热摊进去，液相分数从 0 到 1 平滑过渡，超过气化温度的部分用能量"扣除"来等效材料移除。这种方法不需要动网格，稳定、便宜，缺点是坑的形状出不来，表面始终是平的。

进一层是动网格：靶面设一个随烧蚀速率后退的边界，后退速度由表面能量平衡决定，v = (I_abs - 传导与辐射损失)/(ρ·L_v)。这里 L_v 是气化潜热，量级在 1e7 J/kg，算出来的速度通常是米每秒到几十米每秒，10 ns 内后退几十到几百纳米。动网格的问题是网格畸变，后退量大时必须配合重划分或者映射。

再进一层是水平集或相场追踪界面，能处理熔池飞溅、深孔形成这类复杂形貌。这一层我只在专门研究形貌演化时才用，耦合量太大，容易失控。

阈值判断有个实用技巧：不要用硬阶跃函数，用平滑过渡的软化函数把"烧蚀/不烧蚀"过渡带拉开，宽度取阈值的 5% 到 10%。物理上这是在近似材料表面的不均匀性和脉冲的时间分布，数值上是给求解器留活路。

4. 跑通第一个算例：分阶段实操流程

4.1 第一步：只算传热，拿到温度场

我几乎从不一次性建完整个模型。第一步只做瞬态传热，几何就是一块平板，热源用 Beer-Lambert，边界给对流和辐射。

建完之后先做三件事检查：一是能量守恒，把热源积分出来的总能量和 E0 乘吸收率对一下，差太多说明表达式写错了；二是看靶面中心温度曲线，形状应该是快速上升、脉宽结束附近到峰、然后按热扩散时间常数衰减，如果出现锯齿或者负温度，说明时间步或者热源符号有问题；三是看熔池深度随时间的变化，和解析估算的 √(Dτ) 对比，量级对不上就回头查吸收率和热导率。

这一步跑通大概需要半天到一天，但它给你的是后面所有耦合的基础。温度场错误，后面加什么场都是错的。

4.2 第二步：加流体，看熔池与羽流

在传热基础上加层流，把熔融区当成低黏度液体，表面加马兰戈尼应力项（表面张力随温度的梯度乘以温度切向梯度）。这一步最容易崩的地方是物性突变：固相黏度极大、液相黏度很小，直接给会产生刚性问题，我一般把固相区域的黏度设成一个很大的有限值（比如 1e6 Pa·s），靠糊状区平滑过渡。

另一个坑是自由表面。羽流膨胀需要一个开放边界，但如果把整个上边界都设成开放，环境气体会倒灌，流场变得很脏。我的做法是给一个远场压力边界，并且把计算域的上边界放到光斑尺寸的 5 到 10 倍以外，让膨胀波有地方走。

马兰戈尼对流要不要开，我的经验是先关掉跑一遍，记录熔池形状作为基准；再打开跑一遍，对比深度和宽度差异。差异小于 10% 就可以关掉省算力，差异大就必须开——这说明表面张力驱动是这个工况的主导机制，关掉等于把物理丢了。

4.3 第三步：加等离子体与化学

等离子体部分我建议先在 1D 或者简化 2D 上单独跑通，再嵌回主模型。

单独跑的时候，域内初始电子密度不能给 0，给个种子值 1e15 到 1e17 m⁻³，来源可以是环境电离或者残留电荷。边界条件里，靶面给二次电子发射或者热电子发射，远场给接地或者零通量。电子能量方程里必须包含弹性碰撞、激发、电离三类损失项，少一项电子温度就会偏高。

化学部分先用简化机理，比如只保留三到五个主要反应，跑通之后再逐步加。机理越全，刚性越强，收敛越难。我见过有人一上来挂三十个反应，结果电子密度在 1e-10 和 1e10 之间来回跳，那不是物理，那是数值噪声。

嵌回主模型时，等离子体区和中性气体区的界面要特别处理。我的做法是设一个密度阈值，超过阈值的区域启用等离子体输运，以下的区域只用中性流体，交界处用平滑过渡函数混合输运系数。这样避免在全域求解刚性的漂移-扩散方程。

4.4 求解器设置与收敛控制

全耦合求解这类问题是自找麻烦，我用的是分离式求解，按物理场的响应时间排序：先解传热（快），再解流场（慢），最后解等离子体和化学（最慢）。每个子步的松弛因子从 0.5 开始，收敛困难时降到 0.2。

时间和空间的离散，我一般这样配：时间用 BDF 自由阶次，最大阶数 2，初始步长给 1 ps，允许求解器自己放大；空间用二次拉格朗日单元，等离子体区域的电子密度和电子温度用对数变量，数量级跨度太大时对数化几乎是必须的。

推荐配置（2D 轴对称，10 ns 脉宽）： 时间步进 : BDF 自由，最大阶 2，初始 1 ps 线性求解器 : 直接求解器（PARDISO），规模过大改迭代 分离式耦合 : 传热 -> 层流 -> 等离子体，松弛 0.3~0.5 变量缩放 : 电子密度、电子温度取对数

非线性发散时，第一反应不该是降时间步，而是看残差曲线是哪个变量先炸。如果电子温度先炸，多半是能量方程里少项或者单位错了；如果压力先炸，多半是边界条件过约束。

4.5 后处理：电子密度、温度、发射谱线

后处理阶段我最常看的四个量：靶面中心温度随时间曲线、熔池深度与宽度、电子密度和电子温度的时空分布、以及轴向的等离子体膨胀速度。

发射谱线是验证模型最有力的工具。给定电子温度和密度，用碰撞辐射模型算出某条谱线的强度，再和实验的光谱仪数据对比。这里要注意谱线展宽机制：多普勒展宽反映离子温度，斯塔克展宽反映电子密度，仪器展宽必须从测量值里扣掉。对比的时候不要只比峰值强度，要比半高宽，因为强度受绝对标定影响太大。

动画方面，我习惯输出电子密度的对数等值面随时间演化，能直观看到等离子体从靶面剥离、向外膨胀、最后分裂成团块的过程。这个过程和实验的高速摄影对得上，模型的可信度就上了一个台阶。

5. 踩过的坑与排查速查表

5.1 不收敛与报错的典型场景

第一个高频坑是"找不到一致的初始值"。表现是求解开始就报错，残差从第一步就是 NaN。原因往往是边界条件互相矛盾，比如靶面给了固定温度，同时热源又在同一个边界上注入能量；或者流场的入口给了速度、出口又给了压力，中间还有个壁面把流道堵死。排查办法是把边界条件逐个禁用，看哪个一关就正常。

第二个坑是负密度、负温度。漂移-扩散方程在强场下会把数值推出物理区间。解决办法是给变量加下限约束，或者直接对数化。对数化之后要小心初始值不能给 0，给个 1e10 这种小量兜底。

第三个坑是动网格畸变导致网格质量为负。烧蚀后退量大、或者熔池表面变形剧烈时必现。我的处理是限制单步最大位移（不超过相邻单元尺寸的 20%），配合自动重划分，重划分频率不要太高，否则引入的插值误差会累积。

第四个坑是求解到某个时刻突然变慢。通常是某个物理场出现了刚性，比如化学反应的速率常数在某个温度下突然指数增长。看求解日志里的步长变化，如果步长被压缩到 1e-15 量级，基本可以确定是刚性问题，需要检查反应速率表达式或者加稳态近似。

5.2 收敛了但物理上不对的坑

比不收敛更可怕的是"跑完了、图很漂亮、结论是错的"。我踩过几个典型的。

坑一：能量不守恒但看不出来。热源的能量密度表达式里用了峰值功率密度，同时又乘了脉宽做时间积分，结果总能量算成了两倍。检查方法是把热源对时间和体积做双重积分，和设定能量对比，误差应小于 5%。

坑二：等离子体屏蔽被算成了"透明"。电子密度明明超过了临界密度，激光还是照常往靶面灌能量。原因是热源项没有和电子密度耦合。修正方法是在吸收系数里加入自由电子贡献项，让吸收系数随电子密度上升而增大，同时在电子密度接近临界值时给一个陡降的透过率。

坑三：羽流膨胀速度快了一个量级。常见原因是把电子温度直接当成了重粒子温度。电子温度可以到几个 eV，重粒子温度只有几千 K，两者差一个量级。如果流体方程里用的温度场是电子温度，声速就错了。正确做法是流体用重粒子温度，等离子体能量方程单独解电子温度。

坑四：熔深和实验差三倍。排查顺序是：吸收率是否正确（金属在熔融前后吸收率会跳变）、相变潜热是否包含、热导率是否用了高温值、等离子体屏蔽是否漏算。我遇到过的实际原因是吸收率用了室温值 0.3，而实际熔融后接近 0.7，能量直接少算一半。

5.3 常见问题速查表

提示：这张表我贴在显示器边上，每次模型出问题先扫一遍，能省掉大量盲目试错的时间。

6. 模型验证与扩展方向

6.1 拿什么验证：三个可对比的量

模型跑通不等于模型可信。我一般找三个量做交叉验证，成本从低到高。

烧蚀坑形貌最容易拿到，实验做完打白光干涉或者共聚焦显微镜就有。对比时看深度、直径、深径比三个指标，误差在 20% 以内我认为可以接受，因为烧蚀本身对材料不均匀性极其敏感。等离子体发射谱线的时间积分次之，需要光谱仪，对比的是谱线相对强度和展宽，能间接约束电子密度。羽流前沿速度最难，需要高速摄影或者阴影法，但它的约束力最强，因为它直接反映膨胀动力学。

三个量里如果有两个对不上，不要急着调参数，先回去检查物理假设，尤其是吸收率和屏蔽项。调参数凑结果是最容易走上歧路的一步。

6.2 参数化扫描与实验设计

模型稳定之后，最有价值的事情是做参数扫描。我一般扫三类：能量密度（从阈值附近到明显烧蚀）、脉宽（从皮秒到百纳秒）、光斑半径。扫描的时候用批量求解，把结果存成表格，画等值线图看趋势。

扫描设计要和实验对齐，否则模型给出一堆漂亮的曲线，实验没法验证。我的习惯是先和做实验的同事确认他们能改哪几个参数、测量精度多少，然后只在那个范围内扫，别浪费算力去算一个实验根本做不到的工况。

6.3 后续可扩展的方向

这个模型搭好之后，往下可以走三条路。一是加双温模型，把电子和晶格温度分开，研究飞秒到皮秒尺度的非平衡烧蚀；二是把等离子体部分换成详细的重粒子输运，做谱线定量反演；三是把靶材换成多层膜或者复合材料，研究界面处的热应力与分层。

第三条路我最近用得比较多，因为多层结构里的热阻和反射率变化会让烧蚀阈值出现明显偏移，纯实验很难拆解每一层的贡献，模型反而更有优势。扩展的时候注意一点：每加一层物理，都要重新做一次网格和时间步的敏感性分析，别指望原来的设置还能用。

我个人在实际操作中的体会是，这类多物理场模型最大的成本不在算力，而在"判断哪些物理可以省"。省对了，模型既快又准；省错了，跑得再久也是错的。所以每次动手前那十分钟的量级估算，千万别跳过。