添加阻尼系数并修改时步

在计算过程中增加阻尼系数以阻止流量计算结果的震荡。添加阻尼系数后时步最大可取至1800，比较奇怪的地方是当时步取900时计算结果错误。当时步去3600时，如果流量变化不剧烈，计算结果还可以，如果流量变化很剧烈，计算结果就会发生错误。这证明了虽然隐式算法是无条件稳定的，但是在选取时间步长时仍然需要考虑边界条件，不能任意选取。