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)會(huì)有1E-15量級(jí)的誤差，統(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仍然存在振蕩，但只在兩個(gè)確定的值之間交替變化。對(duì)于監(jiān)測點(diǎn)監(jiān)測文件（44.25,0.25,200,1）fT對(duì)應(yīng)的兩個(gè)值為

3.176594418654839e-78
3.176594307434596e-78

界面fh對(duì)應(yīng)的兩個(gè)值為

3.176594309398118e-78
3.176594416691326e-78

gT對(duì)應(yīng)的兩個(gè)值為

1.588297209327433e-78
1.588297153717284e-78

界面gh對(duì)應(yīng)的兩個(gè)值為

1.588297154699046e-78
1.588297208345676e-78

計(jì)算1000步之后仍然保持穩(wěn)定，不會(huì)出現(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在兩個(gè)值之間跳變，分別為

5.487069987901249e-68
5.487069988318061e-68

界面處gh變化值為

5.487069988307720e-68
5.487069987911585e-68

Case 5

解除對(duì)界面面積和網(wǎng)格體積的限定，同時(shí)在計(jì)算界面法向矢量時(shí)不進(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ā)散。

不對(duì)通量和宏觀量進(jìn)行更新，解除系數(shù)限制（a*,b*），更新500步，分布函數(shù)一直保持不變。