gambit里三维图可为边设置边界条件吗
浏览:63046 回答:10
做三维圆柱绕流,想检测圆柱体某一截面(就是一个圆环)上的升、阻力系数,因在fluent里只能直接检测圆柱体的升、阻力系数,而不能设定面或边的,所以想在gambit里把要检测的圆环(线)设出来,但在边界条件设置面板中选边给圆环(线)设定边界条件,在输出网格时,会提示:warn: Boundary entity xxx does not contain any valid entity and is not written.
要能设该怎样设啊?
或者有没有办法直接在fluent里检测圆柱体某一截面(就是一个圆环边)上的升、阻力系数?
谢谢。
再问个弱问题。文献中得到的圆柱体截面上的压力系数(图4和9)是取的圆环(线)上的值吧?
但我得到的圆环上的平均和脉动压力系数值比文献上小很多,均值在-1.5~1(文献是-2.5~1),脉动值在0~0.3(文献是0~0.6),但我的压力系数瞬时值是在-3~1,不知怎么求均值和脉动值后差这么多,是不是我的求解思路错了,请帮我看下:
1)每隔一定时间步以ASCII格式从Fluent中输出圆环上的压力系数;
2)用Matlab编程处理这些数据:
a. 将不同时刻的压力系数作为列构成一个矩阵Y(i,j),这个矩阵上的每一行对应圆环上相同位置处在不同时刻的压力系数,i对应的是圆环上的点数,j对应的是保存的数据个数。
b. 均值X(i,1)是求上面得到的矩阵的每一行的和再除以列数得到的,即算术平均;
c. 脉动值x'(i,1)是矩阵的每一行的元素减去对应的均值(减去同一个数),对差值求平方(Y(i,j)- X(i,1))^2,再对其(每一行)平方和的均值(和除以列数)开方。即,标准差
d. 用‘cart2pol’命令将圆环上的笛卡尔坐标转换为极坐标(弧度),再转化为度,就可画出相应的平均和脉动压力系数图。这步应有问题,极坐标是以右边为0开始的,但流动是左边为0. 这部分的程序如下:
[the,r]=cart2pol(z10ms(:,1),z10ms(:,2));
% z10ms矩阵的前两列是坐标位置,第三列是均值
theta=the*180/pi;
dx=[theta,(z10ms(:,3)];
dx=sortrows(dx);
dx(:,1)=dx(:,1)+180; %平移
plot(dx(:,1),dx(:,2),'k-');