首先，在物理上是守恒的。你的结果不守恒有两种情况。一是进出口差值很小，比方说1e-4，这样的话，可以认为是守恒的，这是由于数值计算的截断误差导致的。二是差值比较大，比方说1e-1，这样的话就是你的模拟过程有问题，要么没有完全收敛（这种可能性最大），只要多算些时间，就可能收敛了。要么是你的其他设置有问题，这就没法一一解释了。
