跳到正文
半导体开放世界

有限元网格、研究类型与求解器收敛策略

Finite Element Mesh, Study Type and Solver Convergence Strategy（FEM / Mesh / Study / Solver）也称 有限元离散、网格收敛、稳态瞬态频域、直接迭代求解器、多物理场求解策略

有限元可信度来自正确弱式、可解析网格、匹配问题的研究类型和受控求解误差；求解器返回成功只说明离散代数系统达到某个准则。

专深出处 7 份资料
报告问题

速览

约 1 分钟

有限元用形函数在单元上近似连续场。网格应分辨薄层、梯度、曲率和波长，并用目标量随细化的变化证明趋稳。稳态、瞬态、频域和特征值回答不同问题；直接或迭代求解器的选择取决于规模、稀疏性、条件数、内存和预条件，而不是固定优劣。

系统理解

约 3 分钟

有限元先把偏微分方程乘测试函数并在域上积分，导数通过分部积分降低阶次，边界通量进入弱式；再以单元形函数展开未知场，形成稀疏矩阵。单元形状、阶次、积分规则和场之间的离散兼容性会改变误差与稳定性，网格数量本身不是质量指标。

网格应围绕目标量设计：薄层厚度、结区、尖角、接触边缘、热源、材料跳变和结构圆角常需局部加密；长而规则区域可用扫掠或映射网格。全局积分通常比局部峰值更快收敛，理想尖角的应力或电场可能随加密持续上升，必须以真实圆角、过程区或积分量重新定义可验证输出。

研究类型对应时间假设。稳态寻找时间导数为零的平衡，瞬态保留初始状态和时间尺度，频域围绕周期小信号响应，特征值求固有模式。把强非线性大信号直接套用线性频域，或从稳态结果推断启动过冲，都是问题定义错误。参数延续、载荷渐增和分步耦合可以改善求解路径，但必须证明最终解与路径无关或公开多稳态。

直接求解器通常稳健但内存随规模迅速增长；迭代求解器可扩展，却依赖谱性质、缩放和预条件。全耦合同时更新所有场，分离求解交替更新场组；前者不自动更准确，后者也不自动遗漏耦合。可信报告同时给出残差、守恒、容差敏感性、网格/时间步收敛和基准或实验验证。

承接几何、域和边界定义，向后连接TCAD离散校准、多物理耦合、优化和模型治理。

公式提示

离散系统常写K(u)u=f；线性化迭代求J(uᵏ)Δu=−R(uᵏ)。相对残差小只说明代数方程接近满足，仍需比较目标量Q在网格h、h/2、h/4下的变化。

关键参数

关键参数
参数单位含义
单元尺寸、质量、阶次与自由度m / — / n / count控制几何和场的离散误差
稳态、瞬态、频域与特征值研究— / s / Hz / Hz匹配目标时间假设
残差、相对/绝对容差与守恒误差— / — / %判断代数解和物理闭合
直接/迭代、预条件与矩阵内存— / — / GB决定稳健性和可扩展性
目标量网格/步长变化与基准误差%建立数值验证证据

常见误区

  • 单元越多结果必然越准确
  • 彩色云图光滑就说明网格独立
  • 直接求解器总比迭代求解器快且准确
  • 稳态结果可以代表启动瞬态
  • 残差低等于守恒和实验都正确

工程深入

12 条结论 · 1 道例题

有限元弱式

定义

有限元以单元形函数近似满足方程的连续场

弱式把控制方程与测试函数积分并将通量边界纳入，再把未知场展开成有限自由度的稀疏代数系统。12

成立条件 弱式和自然边界一致；形函数空间满足问题要求

网格设计

机制

网格必须分辨几何尺度和场变化尺度

薄层、结区、材料跳变、热源与接触边缘需要局部加密，远离梯度的规则区域可使用较粗或扫掠网格。13

成立条件 关键长度尺度已估计；单元质量和增长率受控

目标量收敛

边界

网格收敛必须追踪决策所用目标量

全局积分、端口量、平均温度和局部峰值的收敛速度不同，不能用一个易收敛量替所有输出签核。23

成立条件 至少三档系统细化；几何物理和容差保持一致

奇异点

边界

理想尖角局部峰值可能永不网格收敛

零半径角点、点载荷或材料不连续可产生数学奇异，应引入真实圆角、有限载荷区或改用积分/距离外目标量。23

成立条件 识别峰值位置和几何理想化；不把单节点极值作为材料签核

研究类型

定义

稳态瞬态频域和特征值回答不同问题

稳态求平衡，瞬态保留初值和时间演化，频域描述周期小扰动，特征值给固有模式；它们不能仅因计算方便而互换。14

成立条件 激励与线性化假设明确；时间尺度或频带覆盖目标

非线性路径

机制

载荷渐增和参数延续是求解路径而不是新物理

把强非线性逐步开启可改善初值和收敛，但应检查不同步长、方向或初值是否到达同一分支。15

成立条件 记录延续参数和路径；检查多稳态和回滞

线性求解

边界

直接与迭代求解器没有脱离问题规模的绝对优劣

直接法通常稳健但耗内存，迭代法可扩展却依赖条件数、缩放和预条件；两者都应在相同误差准则下比较。12

成立条件 矩阵规模和稀疏性可记录；残差与目标量容差一致

耦合求解

边界

全耦合和分离求解改变算法而不应改变合格物理解

全耦合同时更新所有未知量，分离法交替更新场组；充分收敛且解唯一时应得到一致目标量。14

成立条件 相同离散和最终容差；检查耦合残差与守恒

残差边界

边界

代数残差低不能证明离散误差小

求解器可精确求得粗网格上的错误近似，必须分别检查容差敏感性、网格/时间步收敛和物理守恒。23

成立条件 残差按变量缩放；目标量对容差和离散分别扫描

后处理

边界

插值平滑和颜色范围会改变视觉但不会增加数值信息

云图应同时报告网格、单位、色标、极值定义和取样方式，平滑显示不能隐藏不连续或奇异。12

成立条件 保留原始场和派生表达式；比较时固定色标和参考面

收敛算例

公式

相邻网格变化减小可证明趋稳但不是严格误差上界

三档网格的目标量可用于判断趋势，若阶次和渐近区未知，不能把最后一次变化直接宣称为真实误差。23

δ12=3.03%，δ23=1.00%

δ
相邻网格目标量相对变化，单位 %

成立条件 Q=9.6、9.9、10.0；同一模型和容差

例题：动手算一次教学假设，可复算

粗中细网格的热流为9.6、9.9、10.0 W，相邻变化多少？

已知

  1. Q1=9.6 W
  2. Q2=9.9 W
  3. Q3=10.0 W

步骤

  1. δ12=|9.9−9.6|/9.9=3.03%，粗到中网格。
  2. δ23=|10.0−9.9|/10.0=1.00%，中到细网格。
结果 相邻变化约3.03%与1.00%。

该积分量正在趋稳。

不能外推到 仍需检查更细网格、局部峰值、求解容差和实验偏差。

验证层级

关联

数值验证和实验确认必须分开

网格、时间步和基准解回答方程是否算对，实验回答方程材料和边界是否代表真实对象，二者缺一不可。26

成立条件 基准与实验数据独立；误差来源分项报告

出处

7 份资料 · 7 处引用
  1. 1
  2. 2
  3. 3
  4. 4

    The MEMS Module

    COMSOL · 2025

    查看这个来源 →
  5. 5

    The Semiconductor Interface

    COMSOL · 2025

    查看这个来源 →
  6. 6

    Multiscale Modeling and Validation of Semiconductor Materials and Devices

    NIST Semiconductor Metrology Program · 2026

    查看这个来源 →
  7. 7

    Charon User Manual: Semiconductor Device Simulation

    Sandia National Laboratories · 2021

    • Charon User Manual §§13.3–13.4，physics, material models and discretization 用于本页概述
    查看这个来源 →