24小时热门版块排行榜

返回列表

嘟嘟张

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8
帖子: 13
在线: 7.1小时
虫号: 3324490
注册: 2014-07-15
专业: 多相流热物理学

[求助] 两相流UDF编写求助已有1人参与

二维圆柱水滴撞击问题，编写UDF读取壁面第二层网格的水滴速度、密度和体积分数等，现在已经找到第二层网格了，但是读取值的时候还是有问题，求各位大神指导！

#include "udf.h"

DEFINE_ON_DEMAND(test)
{
Domain *d;
d = Get_Domain(3);/*water droplet phase*/

face_t f;
cell_t c;
cell_t  c0, c1= -1;  /* 定义相邻面index,是个整数,内部面矢量从c0指向c1 */
Thread *t0, *t1 = NULL; /*pointer to gouping of cell*/

real NV_VEC(A);
real NV_VEC(A0);
real NV_VEC(A1);
real NV_VEC(A2);
real NV_VEC(A3);
real NV_VEC(A4);

int ID = 18;
Thread *t = Lookup_Thread(d, ID);

Thread *tm = THREAD_SUPER_THREAD(t);
/* Thread tm will be a superthread, i.e., mixture */
Thread **pt = THREAD_SUB_THREADS(tm);

begin_f_loop(f, t)
{
F_AREA(A, f, t);
c0 = F_C0(f, t);
t0 = F_C0_THREAD(f, t);
c1 = F_C1(f, t);    /* Get cell on other side of face */
t1 = F_C1_THREAD(f, t);

//循环C里面的每一个face 先确定 boundary face, fabf ；读取法向量nbf；对每一个face，读法向量；找到对面的face；确定对面的C，c2 计算通量，UC2，//

int i = 3;
int *n = 0;

c_face_loop(c0, t0, i) /* loops over all faces of a cell */
{
face_t f0, f1, f2, f3;
Thread *tf0, *tf1, *tf2, *tf3;
f0 = C_FACE(c0, t0, 0);
tf0 = C_FACE_THREAD(c0, t0, 0);
f1 = C_FACE(c0, t0, 1);
tf1 = C_FACE_THREAD(c0, t0, 1);
f2 = C_FACE(c0, t0, 2);
tf2 = C_FACE_THREAD(c0, t0, 2);
f3 = C_FACE(c0, t0, 3);
tf3 = C_FACE_THREAD(c0, t0, 3);

F_AREA(A0, f0, tf0); /*get vector A value*/
F_AREA(A1, f1, tf1); /*get vector A value*/
F_AREA(A2, f2, tf2); /*get vector A value*/
F_AREA(A3, f3, tf3); /*get vector A value*/
/*确定f1是第一层网格上平行于壁面的face*/

if (BOUNDARY_FACE_THREAD_P(tf1)==0) /*Get interior face values*/
{
cell_t  c2, c3 ;
Thread *t2, *t3;

c2 = F_C0(f1, tf1);
t2 = F_C0_THREAD(f1, tf1);
c3 = F_C1(f1, tf1);    /* Get cell on other side of face */
t3 = F_C1_THREAD(f1, tf1);

F_AREA(A4, f1, tf1); /*get vector A value*/

/*确定c3是壁面第2层网格*/

C_UDMI(c0, t0, 0) = C_VOF(c0, t0);
C_UDMI(c0, t0, 1) = C_U(c0, t0);
C_UDMI(c0, t0, 2) = C_R(c0, t0);
/*能正常读取c0的值，但是c3的读不出来*/

C_UDMI(c0, t0, 3) = C_U(c3, t3);
C_UDMI(c0, t0, 4) = C_R(c3, t3)；
C_UDMI(c0, t0, 5) = C_VOF(c3, t3)；

}
else
{

}

}

}
end_f_loop(f, t)

}

回复此楼