修正质量流量密度计算方法

修正质量流量密度计算方法，能够计算得到正确的结果。还存在两点问题：1.流量计算结果趋势正确，但是存在震荡。2.时步设置为600时计算结果发散。可以尝试采用设置阻尼系数的方法来解决这两个问题。