瞬态模拟连续性方程维持在10-2到10-3之间，达不到10-3，该怎么解决

做液氧的泄露模拟，初始化风场后开始泄露，2.5kg/s，时间步为0.001s，刚开始可以收敛到10-3，泄露0.1s后，连续性方程残差逐渐升高超过10-3。计算域为40m✖️200m✖️2m，泄露口较小，附近加密网格，网格总数70万。要求的残差最大10-3，收敛因子也已经调小，时间步长已经是0.001s，采用piso算法，接下来是需要继续加密网格吗？还有在泄露计算过程中会出现，在多个网格中湍流粘度大于10的5次方，问题出在了哪里...