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

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 各位大侠帮我看看这个UDF到底哪错了？编译没有问题，一初始化计算就致命错误！！
【调剂】北京石油化工学院2024年16个专业接受调剂
31/1返回列表
查看: 966  |  回复: 2
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

459051993

版主 (职业作家)

小木虫打酱油者

[交流] 各位大侠帮我看看这个UDF到底哪错了？编译没有问题，一初始化计算就致命错误！！

#include "udf.h"
#include "sg.h"
#include "sg_mphase.h"
#include "flow.h"
#include "mem.h"
/**************************************************************/
/* UDF for specifying an interfacail area density */
/**************************************************************/
DEFINE_ADJUST(area_density, domain)
{
Thread *t;
Thread **pt;
cell_t c;
Domain *pDomain = DOMAIN_SUB_DOMAIN(domain,P_PHASE);
real voidx, voidy, voidz=0;
real evap_rate = 0.0;
real int_area = 0.0;
real h_t_rate = 0.0;
{
Alloc_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
Scalar_Reconstruction(pDomain, SV_VOF,-1,SV_VOF_RG,NULL);
Scalar_Derivatives(pDomain,SV_VOF,-1,SV_VOF_G,SV_VOF_RG,
Vof_Deriv_Accumulate);
}
{
Alloc_Storage_Vars(domain, SV_T_RG, SV_T_G, SV_NULL);
T_derivatives(domain);
Free_Storage_Vars(domain, SV_T_RG, SV_NULL);
}
mp_thread_loop_c (t,domain,pt)
if (FLUID_THREAD_P(t))
{
Thread *tp = pt[P_PHASE];
begin_c_loop (c,t)
{
C_UDMI(c,t,0) = (C_VOF_G(c,tp)[0]*C_T_G(c,t)[0]+
C_VOF_G(c,tp)[1]*C_T_G(c,t)[1]+C_VOF_G(c,tp)[2]*C_T_G(c,t)[2]);
C_UDMI(c,t,0) = (C_VOF_G(c,tp)[0]*C_T_G(c,t)[0]+
C_VOF_G(c,tp)[1]*C_T_G(c,t)[1]);
C_UDMI(c,t,3) = NV_MAG(C_VOF_G(c,tp))*C_VOLUME(c,t);
if(C_VOF(c,tp)>0 && C_VOF(c,tp)<1)
{evap_rate = evap_rate + C_UDMI(c, t, 4);
int_area = int_area + C_UDMI(c,t,3);
h_t_rate = h_t_rate + C_UDMI(c,t,6); }
C_UDMI(c,t,7) = C_UDMI(c,t,6)/C_UDMI(c,t,3);
}
end_c_loop (c,t)
}
/* real avg_ht_rate = h_t_rate/int_area;*/
/* if(N_ITER%30==0)
printf ("%e\n",evap_rate);*/
Free_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL);
Free_Storage_Vars(domain, SV_T_G, SV_NULL);
}
DEFINE_SOURCE(gas, cell, thread, dS, eqn)
{
real x[ND_ND];
real source;
Thread *tm = THREAD_SUPER_THREAD(thread);
Thread **pt = THREAD_SUB_THREADS(tm);
real Kl = C_K_L(cell, pt[1])*C_VOF(cell, pt[1]),
Kg = C_K_L(cell, pt[0])*C_VOF(cell, pt[0]);
real L = 2270e3;
source = -(Kl+Kg)*C_UDMI(cell,tm,0) / L;
C_UDMI(cell, tm, 1) = source;
C_UDMI(cell,tm,4) = source*C_VOLUME(cell,tm);
C_UDMI(cell,tm,5) = C_VOLUME(cell,tm);
C_UDMI(cell, tm, 2) = -source*L;
C_UDMI(cell,tm,6) = source*L*C_VOLUME(cell,tm);
dS[eqn] =0;
return source;
}
DEFINE_SOURCE(liquid, cell, thread, dS, eqn)
{
real x[ND_ND];
real source;
Thread *tm = THREAD_SUPER_THREAD(thread);
Thread **pt = THREAD_SUB_THREADS(tm);
source = -C_UDMI(cell, tm, 1);
dS[eqn] = 0;
return source;
}
DEFINE_SOURCE(energy, cell, thread, dS, eqn)
{
real x[ND_ND];
real source;
Thread *tm = thread;
source = C_UDMI(cell, tm, 2);
dS[eqn] = 0;
return source;
}
DEFINE_SOURCE(energy_tsat, cell, thread, dS, eqn)
{
real x[ND_ND];
real source;
Thread *tm = thread;
Thread **pt = THREAD_SUB_THREADS(tm);
if(C_VOF(cell, pt[1]) == 0)
{
source = pow(10,30)*(373.15 - C_T(cell,tm));
dS[eqn] = -pow(10,30);
}
else source = 0.0;
return source;
}
回复此楼

» 猜你喜欢
进化论进行中
1楼2019-01-12 20:06:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lovess0928

银虫 (小有名气)
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
459051993: 金币+5, 谢谢 2019-03-21 20:38:49
楼主大大你好，我用了你的UDF试算了一下，发现UDF是没有问题的，不知是否是你的算例设置是否出错了。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

459051993
2楼2019-02-20 16:24:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

459051993

版主 (职业作家)

小木虫打酱油者
送红花一朵
引用回帖:
2楼: Originally posted by lovess0928 at 2019-02-20 16:24:00
楼主大大你好，我用了你的UDF试算了一下，发现UDF是没有问题的，不知是否是你的算例设置是否出错了。

你好，确实是我弄错了，应该是并行运算，谢谢啦
赞一下 回复此楼
进化论进行中
3楼2019-03-21 20:38:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 459051993 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件) 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 研二光催化6月底四篇二区什么水平 5+3 wjtab 2024-04-22 8/400 2024-04-24 00:48 by wjtab
[考博] 24年 申博 化学/材料 一作6篇sci +3 wangyp123 2024-04-23 4/200 2024-04-23 23:38 by 莹亮的星空
[论文投稿] 期刊推荐 20+3 木颜尘ip 2024-04-22 5/250 2024-04-23 23:03 by topedit
[公派出国] 東京大学 应用化学专业柳田研究室 招收修士/博士生 +5 蝶梦123 2024-04-19 5/250 2024-04-23 22:40 by 90472jvvh
[论文投稿] 最近遇到这样一个问题 +4 asd123gfa689 2024-04-18 11/550 2024-04-23 20:14 by asd123gfa689
[有机交流] 紧急求助，有谁用过三甲基硅醇钾TMSOK吗？ 1000+3 genius2008 2024-04-19 8/400 2024-04-23 16:57 by czyzsu
[基金申请] 前同事拿我以前未中的国自然标书申请今年的国自然，现在基金委发来调查函 +28 bjdxyxy 2024-04-19 44/2200 2024-04-23 16:29 by 淀粉搬运工
[考研] 本科211 277分求调剂 +3 0854通信工程 2024-04-18 6/300 2024-04-23 15:39 by H.W.Y
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +4 我属驴核动力驴 2024-04-22 4/200 2024-04-23 12:18 by lyfbangong
[电化学] 有偿收购废弃锂硫电池 +3 汤圆喵 2024-04-19 5/250 2024-04-23 10:41 by Happy C
[教师之家] 国内科研环境差导致国内科研水平差的根源 +8 zju2000 2024-04-18 8/400 2024-04-23 10:40 by Quakerbird
[论文投稿] Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？ 10+3 yifeng11 2024-04-20 7/350 2024-04-23 09:29 by bobvan
[论文投稿] 国家重点研发计划论文可以发MDPI里的期刊吗 +6 lin_2018 2024-04-17 6/300 2024-04-22 17:32 by 20094958
[找工作] 化工专业 10+3 xiao太阳94 2024-04-18 7/350 2024-04-22 13:48 by mlnjut
[考研] 博士生招生！！ +3 2458827884 2024-04-19 3/150 2024-04-22 12:37 by 钰波love
[有机交流] 兄弟们帮我看看这两个结构怎么合成 +6 xl2088131 2024-04-17 8/400 2024-04-21 22:39 by zyqchem
[硕博家园] 材料工程，二本硕士，值得读吗 +3 简风23 2024-04-20 6/300 2024-04-21 12:19 by 简风23
[考研] 298求调剂 +10 孙大大@ 2024-04-17 11/550 2024-04-20 18:33 by 371122gyd
[基金申请] 估计今年青基又没戏 +11 忆念7 2024-04-18 11/550 2024-04-20 17:29 by mygod_nju
[论文投稿] SCI 评审意见回复 模板 15+3 Kilig0317 2024-04-18 4/200 2024-04-19 08:46 by bobvan
信息提示
关闭
请填处理意见
关闭