20170920

Case 1

测试均匀流算例，初始值为激波结构下游平衡态，更新213步后，在xc = 28.25，yc = 0.25处，离散速度分布函数uIndex = 200, vIndex = 198处fT，gT正负交替振荡。监测文件（44.25,0.25,200,1），（44.25,0.75,200,1）显示fT和gT也出现振荡。

Case 2

测试均匀流算例，初始值为激波结构上游平衡态，更新240步后，在xc = 33.75, yc = 0.25处，DDF uIndex = 200, vIndex = 195处fT，gT正负交替震荡。监测文件（44.25,0.25,200,1），（44.25,0.75,200,1）显示fT和gT也出现振荡。

由于pointwise导出网格过程中存在误差，导致网格点的坐标会有1E-15量级的误差，统一修正界面面积为常值0.5，单元体积为0.25。

Case 3

测试均匀流算例，初始值为激波结构上游平衡态，监测文件（44.25,0.25,200,1），（44.25,0.75,200,1）显示fT和gT仍然存在振荡，但只在两个确定的值之间交替变化。对于监测点监测文件（44.25,0.25,200,1）fT对应的两个值为

3.176594418654839e-78
3.176594307434596e-78

界面fh对应的两个值为

3.176594309398118e-78
3.176594416691326e-78

gT对应的两个值为

1.588297209327433e-78
1.588297153717284e-78

界面gh对应的两个值为

1.588297154699046e-78
1.588297208345676e-78

计算1000步之后仍然保持稳定，不会出现负值。

Case 4

测试均匀流算例，初始值为激波结构下游平衡态，500步之后，监测文件（44.25,0.25,200,1），（44.25,0.75,200,1）显示fT为恒定值无振荡，为9.184982020200940e-68。
界面fh的值与fT相同，gT在两个值之间跳变，分别为

5.487069987901249e-68
5.487069988318061e-68

界面处gh变化值为

5.487069988307720e-68
5.487069987911585e-68

Case 5

解除对界面面积和网格体积的限定，同时在计算界面法向矢量时不进行0值判定(注释SetZero函数)，测试均匀流算例，初始值为激波结构上游平衡态，监测文件显示，fT,gT发散。

Case 6

限制界面面积和网格体积为恒定，测试激波结构算例，监测文件显示，fT,gT发散。

Case 7

限制界面面积和网格体积为恒定，初始值为激波结构下游平衡态，界面处所有分布函数亦为平衡态，跳过宏观量（格心和界面）更新步骤，结果显示，残差一直为0，监测文件显示，fT,gT,fh,gh在70步之后出现负值，发散。

进一步限制(aTP = 1,bTP = 0),(ah = 1,bh = 0),(aBP = 1,bBP = 0)，使得在演化过程中与平衡态有关的系数都为0，排除平衡态分布函数的影响，测试结果表明，70步之后发散。

不对通量和宏观量进行更新，解除系数限制（a*,b*），更新500步，分布函数一直保持不变。