为何我的continuity方程残差总是下不去？

aleisia

我的问题是，一个大空间（计算域）里面有个小房间，小房间的外面某处有个源项释放species， 200g/m3的速率，小房间的另一边外面有一个车，吸风，处理species，再把处理好的干净空气提供给小房间的室内各个小隔间。
车和小房间是有输风管干管和支管连接的。
小房间外顶上还有三个风扇，旁边一个水平面上还有一个风机，吹风。
我的边界条件是：风机和风扇的进口是fan，出口是interior。species源项是指定这个体积为fluid，source term。六个表面是interior。
车的吸风口我也设为fan，各个输风管的支管设为interior。
所有体都是互相split成的，所以截出来的面都是默认的interior。
我还在输风管的干管上加了一个sink项，以-199g/m3的速率移除那些species。剩下的这部分1g/m3，就让它随着管道被风带到各个小隔间里。

如此开始计算，continuity残差无论如何都下不到10（-1），可是看速度残差都在10（-4），species 残差在10（-6）档次。
我试过采用不同的初始值，比如表压改为60Pa再重新开始initialization后再计算等等，都无法把continuity给收敛好。
我网上搜过，有的人认为continuity不到10（-3）也不一定意味着不收敛，要监控各个关心的面，如果那些面的量守恒了，也可认为收敛。
而我的一个很有经验的朋友和我说，continuity还是要达到10（-3）以下的好。
请问大家，我该如何改进？
是边界条件我设置的有误？还是什么别的地方？
我的网格应该是没问题的，该通的都通，该是interior的地方都是interior。
谢谢！