20170920
Case 1
測試均勻流算例怀估，初始值為激波結(jié)構(gòu)下游平衡態(tài)损姜，更新213步后条获，在xc = 28.25忠荞，yc = 0.25處，離散速度分布函數(shù)uIndex = 200, vIndex = 198處fT帅掘，gT正負(fù)交替振蕩委煤。監(jiān)測文件（44.25,0.25,200,1），（44.25,0.75,200,1）顯示fT和gT也出現(xiàn)振蕩修档。
Case 2
測試均勻流算例碧绞，初始值為激波結(jié)構(gòu)上游平衡態(tài)，更新240步后吱窝，在xc = 33.75, yc = 0.25處讥邻，DDF uIndex = 200, vIndex = 195處fT迫靖，gT正負(fù)交替震蕩。監(jiān)測文件（44.25,0.25,200,1）兴使，（44.25,0.75,200,1）顯示fT和gT也出現(xiàn)振蕩系宜。
由于pointwise導(dǎo)出網(wǎng)格過程中存在誤差，導(dǎo)致網(wǎng)格點(diǎn)的坐標(biāo)會有1E-15量級的誤差发魄，統(tǒng)一修正界面面積為常值0.5盹牧，單元體積為0.25。
Case 3
測試均勻流算例励幼，初始值為激波結(jié)構(gòu)上游平衡態(tài)汰寓，監(jiān)測文件（44.25,0.25,200,1），（44.25,0.75,200,1）顯示fT和gT仍然存在振蕩赏淌，但只在兩個確定的值之間交替變化踩寇。對于監(jiān)測點(diǎn)監(jiān)測文件（44.25,0.25,200,1）fT對應(yīng)的兩個值為
3.176594418654839e-78
3.176594307434596e-78
界面fh對應(yīng)的兩個值為
3.176594309398118e-78
3.176594416691326e-78
gT對應(yīng)的兩個值為
1.588297209327433e-78
1.588297153717284e-78
界面gh對應(yīng)的兩個值為
1.588297154699046e-78
1.588297208345676e-78
計(jì)算1000步之后仍然保持穩(wěn)定，不會出現(xiàn)負(fù)值六水。
Case 4
測試均勻流算例俺孙，初始值為激波結(jié)構(gòu)下游平衡態(tài)，500步之后掷贾，監(jiān)測文件（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
解除對界面面積和網(wǎng)格體積的限定港准，同時在計(jì)算界面法向矢量時不進(jìn)行0值判定(注釋SetZero函數(shù))旨剥，測試均勻流算例，初始值為激波結(jié)構(gòu)上游平衡態(tài)浅缸，監(jiān)測文件顯示轨帜，fT,gT發(fā)散。
Case 6
限制界面面積和網(wǎng)格體積為恒定衩椒，測試激波結(jié)構(gòu)算例蚌父，監(jiān)測文件顯示，fT,gT發(fā)散毛萌。
Case 7
限制界面面積和網(wǎng)格體積為恒定苟弛，初始值為激波結(jié)構(gòu)下游平衡態(tài)，界面處所有分布函數(shù)亦為平衡態(tài)阁将，跳過宏觀量（格心和界面）更新步驟膏秫，結(jié)果顯示，殘差一直為0做盅，監(jiān)測文件顯示荔睹，fT,gT,fh,gh在70步之后出現(xiàn)負(fù)值狸演，發(fā)散。
進(jìn)一步限制(aTP = 1,bTP = 0),(ah = 1,bh = 0),(aBP = 1,bBP = 0)僻他，使得在演化過程中與平衡態(tài)有關(guān)的系數(shù)都為0宵距，排除平衡態(tài)分布函數(shù)的影響，測試結(jié)果表明吨拗，70步之后發(fā)散满哪。
不對通量和宏觀量進(jìn)行更新，解除系數(shù)限制（a*,b*）劝篷，更新500步哨鸭，分布函數(shù)一直保持不變。
