20170921
8:20 11:40 2:50 + 0:30 = 3:20
Case 1
重新測試均勻流算例垄惧，發(fā)散，不知道哪里改錯了怪得。
查看在另一臺服務器上的程序發(fā)現(xiàn)通砍，昨天測試算例有誤。
cell.fT[i][j] = aTP*cell.fBP[i][j] - bTP*cell.fT[i][j] + cell.DtSlashVolume*cell.fFlux[i][j];
cell.gT[i][j] = aTP*cell.gBP[i][j] - bTP*cell.gT[i][j] + cell.DtSlashVolume*cell.gFlux[i][j];
cell.fT[i][j] = aTP*cell.fBP[i][j] - bTP*cell.fT[i][j] + 4.0*dt*cell.fFlux[i][j];
cell.gT[i][j] = aTP*cell.gBP[i][j] - bTP*cell.gT[i][j] + 4.0*dt*cell.gFlux[i][j];
昨天的算例固定cell volume為0.25薯定，所以dt/volume = 4.0*dt始绍，但是1、2行的內(nèi)容應該注釋掉话侄，昨天測試不發(fā)散的算例都沒有注釋1亏推、2行，這就意味著年堆，fT被更新了兩次吞杭，結果錯誤。
Case 2
注釋上述1变丧、2行代碼芽狗，更新網(wǎng)格文件，保證格點坐標無誤差痒蓬，重新測試均勻流童擎，初始化為激波結構上游平衡態(tài)滴劲，164步之后監(jiān)測點發(fā)散，與之前結果不同的是顾复，fT, fh在46步到145值恒定為9.184982020200943e-68 班挖，gT, gh在52步到145值恒定為5.487069988109647e-68，之后出現(xiàn)振蕩芯砸，進而發(fā)散萧芙。
14:20 18:00 3:10
Case 3
Case 31 ubuntu001
- 出現(xiàn)振蕩的分布函數(shù)基本都位于速度空間的邊緣，絕對值較大假丧，在更新界面通量時双揪，并沒有乘以時間，而是在更新fT的時候將通量乘以dt/V包帚，推測是fh乘以大數(shù)引入了機器誤差渔期，所以將dt移到xi_n_dS中，限制與fh相乘的值渴邦。
- 跳過宏觀量更新擎场，將xi_n_dS的值乘以小數(shù)dt，結果發(fā)散几莽。
- 跳過宏觀量更新，修改界面法矢量求解方法宅静，不歸一化章蚣，220步發(fā)散。
- 修改監(jiān)測點姨夹。
Case 32 ubuntu002
- 跳過宏觀量更新纤垂，界面分布函數(shù)和格心分布函數(shù)設為定值，系數(shù)保持恒定磷账，fT在160步出現(xiàn)振蕩峭沦，gT在219步未發(fā)生振蕩。
- 添加通量輸出函數(shù)逃糟。
- 添加監(jiān)測點周圍單元fT,gT的輸出吼鱼。
- 重新測試，監(jiān)測點位置不佳绰咽，更換監(jiān)測點（12.75,0.25,199,3）菇肃。
- 將xi_n_dS的絕對值設為1，符號保持不變取募，
if(xi·n >=0) 1 else -1
Case 33 ubuntu003
- 跳過通量更新琐谤，保持宏觀量更新，系數(shù)保持恒定玩敏，2500步分布函數(shù)保持恒定斗忌。
20:20 22:30
