24小时热门版块排行榜

返回列表

zhufeiyu_1

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 304.9
帖子: 88
在线: 7.8小时
虫号: 13775313
注册: 2019-01-08
专业: 传热传质学

[交流] 添加连续性方程质量源项UDF编译无错误，进行计算时出现错误已有3人参与

本人目前所做内容为研究地埋管周围土壤热湿耦合传递特性，现添加连续性方程源项UDF编译无错误，进行计算时出现错误：
Error:  received a fatal signal (Segmentation fault).

Error:  received a fatal signal (Segmentation fault).
Error Object: #f

哪位大神能帮忙看看是什么地方得错误，万分感谢！！！
连续性方程源项UDF如下：
#include"udf.h"

DEFINE_SOURCE(source_mass,cell,t,dS,eqn)

{
      real A[ND_ND];
      real source;
      real temp0,temp1,temp2;
      real vol;
      real vof;
      real a0[ND_ND],a1[ND_ND],a2[ND_ND];
      real diff;
      real xv[ND_ND],yv[ND_ND],zv[ND_ND];
      real x;

      Domain *domain;
      Thread *c_thread;
      Thread *subthread;
      cell_t c;
      face_t f;
      Domain *d=Get_Domain(2);    //第一相对应的子域 domain_id 是 2
      Thread *mixture_thread = THREAD_SUPER_THREAD(subthread);

      temp0=C_T_G(cell,c_thread)[0];
      temp1=C_T_G(cell,c_thread)[1];
      temp2=C_T_G(cell,c_thread)[2];
      vol=C_VOLUME(cell,c_thread); //控制体体积；
      vof=C_UDSI(cell,c_thread,0);    //控制体内含湿量；
      x=vof*12.53;
      diff=(7*(-2.41-0.002*pow(vof,-1.75))/0.0718)*(-1.425*pow(10,-4)-4.76*pow(10,-7)*C_T(c,t))*9.48*pow(10,x)/pow(10,10);    //等温扩散系数
      xv[0]=1;xv[1]=0;xv[2]=0;
      yv[0]=0;yv[1]=1;yv[2]=0;
      zv[0]=0;zv[1]=0;zv[2]=1;

      begin_f_loop(f,c_thread)  //历遍混合相的所有面
      {
            F_AREA(A,f,c_thread);    //F_AREA可以被用于返回一个实数的面积向量
            a0[ND_ND]+=temp0*NV_DOT(A,xv);             //x方向的温度梯度与控制体x方向投影面积的乘积
            a1[ND_ND]+=temp1*NV_DOT(A,yv);
            a2[ND_ND]+=temp2*NV_DOT(A,zv);
      }
      end_f_loop(f,c_thread)

      source=diff*(a0[0]+a1[1]+a2[2])/vol;
      dS[eqn]=0;
      return source;
}

回复此楼