找回密码
 注册
查看: 2224|回复: 6

fluent迭代突然停止

[复制链接]
发表于 2014-11-6 18:17:53 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
本帖最后由 学习小能手_二流 于 2014-11-6 18:22 编辑


S

S

 楼主| 发表于 2014-11-6 18:19:08 | 显示全部楼层
怎么贴图啊,坑爹

点评

检查一下你浏览器的“兼容性视图”是否开启,开启以后直接在编辑页面点击上方工具栏里的图片按钮就可以贴图了。  详情 回复 发表于 2014-11-6 19:18
 楼主| 发表于 2014-11-6 18:24:32 | 显示全部楼层
/*迭代时间不同,入口速度也不同*/

#include "udf.h"
#include "mem.h"
#include "metric.h"
int kount=0;
real m;
real r=350.;
real kp=0.01,ki=0.01;
real temp_last;
/* real temp;*/
DEFINE_ADJUST(demo_calc,domain)
{
kount++;
printf("kount=%d\n",kount);
}
DEFINE_ADJUST(best_point_temp,domain)
{
cell_t c;
Thread *t;
real xc[ND_ND];
real temp;
/* loop over all cell threads in the domain*/
thread_loop_c(t,domain)
{
/* loop over all cells*/
begin_c_loop(c,t)
  {
     C_CENTROID(xc,c,t);
    if (kount<50.)
        {
         if(sqrt(ND_SUM(pow(xc[0]-0.07,2.),pow(xc[1]-0.07,2.)))<0.001)
      {
         m=1;
        temp=r-C_T(c,t);
        temp_last=r-C_T_M1(c,t);
       m+=kp*(temp-temp_last)+ki*temp;
       printf("temp_last=%d\n",temp_last);
       }
     else
      {
       return;
       }
        }
          else
      {
       temp_last=r-C_T(c,t);
           printf("temp_last=%n \n",temp_last);
           m=2;
       }
    }
     end_c_loop(c,t)
  }
}

DEFINE_PROFILE(x_velocity,thread,nv)
{
float x[2];
  float y;
  face_t f;
begin_f_loop(f,thread)
{
y=x[1];            /* x[0]=x; x[1]=y; x[2]=z */


        F_PROFILE(f,thread,nv)=m;
}
end_f_loop(f,thread)
}


 楼主| 发表于 2014-11-6 18:26:14 | 显示全部楼层
迭代完50步之后再迭代就出现这个问题,去掉
           temp_last=r-C_T(c,t);
           printf("temp_last=%n \n",temp_last);
就可以迭代,我想知道只是为什么呀。
哭瞎,跪求大神回复!!!!
发表于 2014-11-6 19:18:20 | 显示全部楼层

检查一下你浏览器的“兼容性视图”是否开启,开启以后直接在编辑页面点击上方工具栏里的图片按钮就可以贴图了。
 楼主| 发表于 2014-11-6 19:29:07 | 显示全部楼层
周华 发表于 2014-11-6 19:18
检查一下你浏览器的“兼容性视图”是否开启,开启以后直接在编辑页面点击上方工具栏里的图片按钮就可以贴 ...

站长回复我,小小鸡冻。这个贴图解决了。站长能回答为什么迭代软件会停止工作吗,郁闷,就加了两行语句。
发表于 2014-11-6 22:31:19 | 显示全部楼层
检查一下语法,看看变量类型是否符合。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表