COMSOL中的求解器（1）—— 方程形式

1. 流程

COMSOL中将PDE转成ODE(瞬态仿真)，再通过对时间项离散，最后获得稀疏矩阵方程，通过求解器求解。而稳态仿真则跳过上述时间离散的过程，其余与瞬态仿真求解一致。

流程如下：

瞬态：

稳态:

2. 隐式ODE，及其离散形式

其中，D：质量矩阵；K：刚度矩阵；L{Ω}：载荷；L{partial_Ω}：边界载荷。

将隐式方程L(U对时间的导数, U, t) = 0 进行离散，获得代数系统(代数方程)。

时间离散默认使用后向差分法BDF(隐式欧拉法)，其稳定性好。COMSOL会根据精度，自动在1阶到5阶中选择，下图中的G中的q即代表阶数。

 3. 非线性求解(牛顿法)

如果物性参数具有非线性，比如随温度变化。那么需要求解的就是一个非线性问题，要用牛顿法求解。

F’是标量形式，矢量形式为Jacobian矩阵。

然而牛顿法的求解是否成功，与初始值的选取有关。若初始值没选正确，会导致迭代时发散。因此，COMSOL中可以选择带阻尼的牛顿法求解。

 在COMSOL中可以选择constant 阻尼系数，也可以自动选择阻尼系数。

另外，Jacobian矩阵的更新会很耗时，因此默认使用最少次数更新Jacobian矩阵。如果遇到不收敛，可以手动调整为每次迭代都更新Jacobian矩阵。

下面是非线性求解器设置的技巧：

 4. 瞬态及稳态的方程小结