双向耦合粒子追踪稳态求解器找到未定义的值？

问题描述

模型同时包含静电场和带电粒子追踪模块，当迭代终止方法中的迭代次数＞1时，求解器中稳态求解器报错：

找到未定义的值。

-详细信息：在 中找到未定义的值方程残余矢量。有xxxxxx个自由度，在变量comp1. V 的矢量中为NaN /Inf。

位于坐标：（xx，xx），... 

解决措施

迭代次数低于2次，意味着此时粒子模型并未考虑自身电荷密度变化对原静电场的扰动，所以计算不会报错，一旦迭代次数上升，带电粒子本身的电荷密度扰动被带入静电场泊松方程中来。入口位置的微小数值差异可能会导致后续迭代发生剧烈波动甚至不收敛的情况。

措施一：缩短仿真步长

• 可试着尝试把仿真结束时间步长（0, t_term, t_end中的t_end）缩短后，观察模型是否可正常运行进行多次迭代求解并观察仿真日志中的时间步长变化，对比相较于之前报错模型的步长与该正常运行模型步长之间数量级差异。

措施二：调整粒子入射电流

• 对于等离子而言，应该满足一个最基本属性即：电中性。调整入射电流的目的在于尽可能缩小在入口处离子和电子之间数密度差异，若二者数值差异过大将会产生剧烈的扰动电场，导致变量矩阵内部出现NaN/If。