It remains only the matrix operations at element level, by a strategy of feedback correction, the incompressibility condition is maintained in the filled portion at each time step.

• 对于各类注射填充的模拟,此算法避免了压力场的全局求解,矩阵操作仅在单元一级进行,通过反馈修正,在各时间增量步使已填充的区域满足不可压缩条件。