24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

DJ-Fish

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 61
帖子: 29
在线: 53.6小时
虫号: 1728315
注册: 2012-03-31
专业: 内流流体力学

[求助] 500金币，帮忙找错，UDF并行

单机下能正常使用，单机多核并行报错：
单机UDF如下
/*UDF功能：从出口取出颗粒质量流量给到入口*/
#include "udf.h"
#include "mem.h"
#include "sg.h"
#include "sg_mphase.h"
#include "flow.h"
#include "metric.h"

#define UDF_FILENAME "hxout_s_flux.dat"/*每次开始新的算例要修改名字否则会接着原来的文件写入数据*/

real last_ts1=-1;
DEFINE_PROFILE(inlet_1,t,index)
{

Domain *domain;
Thread *t_s,*t_mix;
face_t f;
double flux,flux_tot,veloc;
real ro_s;
real A[ND_ND];
real area;
real current_time;

current_time=RP_Get_Real("flow-time"

;/*取得当前的计算时间*/
if(last_ts1!=current_time)
{
last_ts1=current_time;
domain=Get_Domain(1);/*1：混合相，2：主相，3：第二相*/
t_mix=Lookup_Thread(domain,5);/*5：出口边界ID*/
t_s=THREAD_SUB_THREAD(t_mix,1);/*0：主相，1：第二相，以此类推*/

begin_f_loop(f,t_mix)
{
flux=F_FLUX(f,t_s);
flux_tot+=flux;
}
end_f_loop(f,t_mix)

begin_f_loop(f,t)
{
ro_s=F_R(f,t);
F_AREA(A,f,t);
area+=NV_MAG(A);
}
end_f_loop(f,t)
veloc=flux_tot/ro_s/area/0.3/2.0;
begin_f_loop(f,t)
{
F_PROFILE(f,t,index)=veloc;
}
end_f_loop(f,t)

}
}

回复此楼

1楼 2013-01-19 19:23:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

DJ-Fish

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 61
帖子: 29
在线: 53.6小时
虫号: 1728315
注册: 2012-03-31
专业: 内流流体力学

小妹愚钝，望高手帮忙指点，不胜感激！！！

赞一下

回复此楼

4楼2013-01-19 19:26:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 16 个回答

DJ-Fish

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 61
帖子: 29
在线: 53.6小时
虫号: 1728315
注册: 2012-03-31
专业: 内流流体力学

以下是小妹修改的其中一个版本

#include "udf.h"
#include "mem.h"
#include "sg.h"
#include "sg_mphase.h"
#include "flow.h"
#include "metric.h"

#define UDF_FILENAME "hxout_s_flux4.0.dat"

real last_ts1=-1;
real veloc;

DEFINE_PROFILE(outlet_pressure,t,index)
{
real flux;

#if !RP_HOST
face_t f;
Thread *t_s;
#endif

#if !RP_HOST

t_s=THREAD_SUB_THREAD(t,1);/*0：主相，1：第二相，以此类推*/

begin_f_loop(f,t)
{
F_PROFILE(f,t,index)=-4000.0;
flux+=F_FLUX(f,t_s);
}
end_f_loop(f,t)
#endif

#if RP_NODE
flux=PRF_GRSUM1(flux);
#endif
#if RP_NODE
flux=PRF_GRSUM1(flux);
#endif

node_to_host_real_1(flux);
veloc=flux/2678.0/0.004145/0.3/2.0;

}

DEFINE_PROFILE(inlet,thread_in,index)
{

face_t inlet_face;

#if !RP_HOST
real time=RP_Get_Real("flow-time"

;
begin_f_loop(inlet_face,thread_in)
{

      F_PROFILE(inlet_face,thread_in,index)=veloc;

}
      end_f_loop(inlet_face,thread_in)
#endif

/*printf("%f\n",inlet_yi); */
}

赞一下

回复此楼

2楼2013-01-19 19:24:47

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖