应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 500金币，帮忙找错，UDF并行
51/1返回列表
查看: 2259  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

DJ-Fish

铁虫 (初入文坛)

[求助] 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

铁虫 (初入文坛)
修改后的UDF 在加上这一句：node_to_host_real_1(flux);之后 报错 ：浮点溢出
赞一下 回复此楼
3楼2013-01-19 19:25:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

DJ-Fish

铁虫 (初入文坛)
以下是小妹修改的其中一个版本

#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的回帖

DJ-Fish

铁虫 (初入文坛)
小妹愚钝，望高手帮忙指点，不胜感激！！！
赞一下 回复此楼
4楼2013-01-19 19:26:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

【答案】应助回帖

★
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢参与！ 2013-01-24 19:44:44
对于你的这个代码，并行与串行的代码应该是一样的。你不妨直接把在单核中可用的代码不作任何改动拿到多核中试试
赞一下 回复此楼
5楼2013-01-20 19:43:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
关闭
请填处理意见
关闭