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步,分布函数一直保持不变。