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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 关于网上空气蒸汽冷凝UDF celltest的问题
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1752  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sunzhijian

新虫 (初入文坛)

[求助] 关于网上空气蒸汽冷凝UDF celltest的问题 已有2人参与

在网上找了一个空气蒸汽冷凝UDF 用于模拟湿空气冷凝的问题，在模拟过程中出现Cpp1.obj : error LNK2019: 无法解析的外部符号 celltest，该符号在函数 temp_info 中被引用，Cpp1.obj : error LNK2019: 无法解析的外部符号 cell_T，该符号在函数 vapor_m_source 中被引。不知道这个celltest代表什么意思，请高手予以帮助，或给个可用的UDF。
#include "udf.h" /*head files*/
#include "mem.h"
#include "flow.h"
#pragma comment(lib,"ws2_32.lib"
#define PI 3.141592653 /*constant*/
#define Mv 18.01534 /*molecular weight of vapor(kg/kmol)*/
#define Ma 28.966 /*molecular weight of air(kg/kmol) */
int *a;
real *b;
real *d;
real *q;
real *g;
real *l;
real s1;
real s2;
real s3;
real s0;
int counter=0; /*counter*/
int j; /*index of interface wall cell*/
int k;
real latent_heat(real T)
{
real t;
real Hlg;
t=T-273.15;
if(t<10.0)
Hlg=2477208.8;
if((t>=10.0)&&(t<20.0))
Hlg=2477208.8+(t-10.0)*(2453549.6-2477208.8)/10.0;
if((t>=20.0)&&(t<30.0))
Hlg=2453549.6+(t-20.0)*(2429882.8-2453549.6)/10.0;
if((t>=30.0)&&(t<40.0))
Hlg=2429882.8+(t-30.0)*(2406058.7-2429882.8)/10.0;
if((t>=40.0)&&(t<50.0))
Hlg=2406058.7+(t-40.0)*(2382038.2-2406058.7)/10.0;
if((t>=50.0)&&(t<60.0))
Hlg=2382038.2+(t-50.0)*(2357691.0-2382038.2)/10.0;
if((t>=60.0)&&(t<70.0))
Hlg=2382038.2+(t-60.0)*(2333080.9-2382038.2)/10.0;
if((t>=70.0)&&(t<80.0))
Hlg=2333080.9+(t-70.0)*(2308065.7-2333080.9)/10.0;
if((t>=80.0)&&(t<90.0))
Hlg=2308065.7+(t-80.0)*(2282560.3-2308065.7)/10.0;
if((t>=90.0)&&(t<100.0))
Hlg=2282560.3+(t-90.0)*(2256472.9-2282560.3)/10.0;
if((t>=100.0)&&(t<110.0))
Hlg=2256472.9+(t-100.0)*(2229704.3-2256472.9)/10.0;
if((t>=110.0)&&(t<120.0))
Hlg=2229704.3+(t-110.0)*(2202149.7-2229704.3)/10.0;
if((t>=120.0)&&(t<130.0))
Hlg=2202149.7+(t-120.0)*(2173700.0-2202149.7)/10.0;
if((t>=130.0)&&(t<140.0))
Hlg=2173700.0+(t-130.0)*(2144243.7-2173700.0)/10.0;
if((t>=140.0)&&(t<150.0))
Hlg=2144243.7+(t-140.0)*(2113667.6-2144243.7)/10.0;
if((t>=150.0)&&(t<160.0))
Hlg=2113667.6+(t-150.0)*(2081855.9-2113667.6)/10.0;
if((t>=160.0)&&(t<170.0))
Hlg=2081855.9+(t-160.0)*(2048687.3-2081855.9)/10.0;
if((t>=170.0)&&(t<180.0))
Hlg=2048687.3+(t-170.0)*(2014031.4-2048687.3)/10.0;
if((t>=180.0)&&(t<190.0))
Hlg=2014031.4+(t-180.0)*(1977745.0-2014031.4)/10.0;
if((t>=190.0)&&(t<200.0))
Hlg=1977745.0+(t-190.0)*(1939668.5-1977745.0)/10.0;
if((t>=200.0)&&(t<210.0))
Hlg=1939668.0+(t-200.0)*(1899623.0-1939668.0)/10.0;
if((t>=210.0)&&(t<220.0))
Hlg=1899623.0+(t-210.0)*(1857409.0-1899623.0)/10.0;
if((t>=220.0)&&(t<230.0))
Hlg=1857409.0+(t-220.0)*(1812799.75-1857409.0)/10.0;
if((t>=230.0)&&(t<240.0))
Hlg=1812799.75+(t-230.0)*(1765537.21-1812799.75)/10.0;
if((t>=240.0)&&(t<250.0))
Hlg=1765537.21+(t-230.0)*(1715325.26-1765537.21)/10.0;
if(t>=250.0)
Hlg=1715325.26;
return Hlg;
}
real s_enthalpy(real T)
{
real t;
real H;
t=T-273.15;
if(t<10.0)
H=2519000.0;
if((t>=10.0)&&(t<20.0))
H=2519000.0+(t-10.0)*(2538000.0-2519000.0)/10.0;
if((t>=20.0)&&(t<30.0))
H=2538000.0+(t-20.0)*(2556000.0-2538000.0)/10.0;
if((t>=30.0)&&(t<40.0))
H=2556000.0+(t-30.0)*(2574000.0-2556000.0)/10.0;
if((t>=40.0)&&(t<50.0))
H=2574000.0+(t-40.0)*(2596000.0-2574000.0)/10.0;
if((t>=50.0)&&(t<60.0))
H=2596000.0+(t-50.0)*(2609000.0-2596000.0)/10.0;
if((t>=60.0)&&(t<70.0))
H=2609000.0+(t-60.0)*(2609000.0-2609000.0)/10.0;
if((t>=70.0)&&(t<80.0))
H=2626000.0+(t-70.0)*(2643000.0-2626000.0)/10.0;
if((t>=80.0)&&(t<90.0))
H=2643000.0+(t-80.0)*(2660000.0-2643000.0)/10.0;
if((t>=90.0)&&(t<100.0))
H=2660000.0+(t-90.0)*(2676000.0-2660000.0)/10.0;
if((t>=100.0)&&(t<110.0))
H=2676000.0+(t-100.0)*(2691000.0-2676000.0)/10.0;
if((t>=110.0)&&(t<120.0))
H=2691000.0+(t-110.0)*(2706000.0-2691000.0)/10.0;
if((t>=120.0)&&(t<130.0))
H=2706000.0+(t-120.0)*(2720000.0-2706000.0)/10.0;
if((t>=130.0)&&(t<140.0))
H=2720000.0+(t-130.0)*(2734000.0-2720000.0)/10.0;
if((t>=140.0)&&(t<150.0))
H=2734000.0+(t-140.0)*(2747000.0-2734000.0)/10.0;
if((t>=150.0)&&(t<160.0))
H=2747000.0+(t-150.0)*(2758000.0-2747000.0)/10.0;
if((t>=160.0)&&(t<170.0))
H=2758000.0+(t-160.0)*(2769000.0-2758000.0)/10.0;
if((t>=170.0)&&(t<180.0))
H=2769000.0+(t-170.0)*(2778000.0-2769000.0)/10.0;
if((t>=180.0)&&(t<190.0))
H=2778000.0+(t-180.0)*(2786000.0-2778000.0)/10.0;
if((t>=190.0)&&(t<200.0))
H=2786000.0+(t-190.0)*(2793000.0-2786000.0)/10.0;
if((t>=200.0)&&(t<=210.0))
H=2793000.0+(t-200.0)*(2798000.0-2793000.0)/10.0;
if((t>=210.0)&&(t<=220.0))
H=2798000.0+(t-210.0)*(2802000.0-2798000.0)/10.0;
if(t>=220.0)
H=2802000.0;
return H;
}

DEFINE_DELTAT(my_timestep,d)
{
real time_step;
real flow_time=RP_Get_Real("flow-time";
if(flow_time<4.0)
time_step=0.01;
else
time_step=0.05;
return time_step;
}

DEFINE_PROFILE(cwalltemp,t,i)
{
real x[ND_ND];
real s,ma,mb,mc,md,me;
face_t f;
ma=8.284;
mb=-25.585;
mc=22.699;
md=-5.3462;
me=350.87;
begin_f_loop(f,t)
{
F_CENTROID(x,f,t);
s=x[1];
F_PROFILE(f,t,i) = ma*s*s*s*s+mb*s*s*s+mc*s*s+md*s+me;
}
end_f_loop(f,t)
}



DEFINE_PROFILE(mass_flow,t,i)
{
face_t f;
begin_f_loop(f,t)
{
F_PROFILE(f,t,i) =-s0;
}
end_f_loop(f,t)
}



DEFINE_ON_DEMAND(INIT)
{
int thread_ID;
Domain *domain;
Thread *tt;
int i;
face_t f; /*face index on the wall*/
cell_t c0;
real A[ND_ND];
i=0;
thread_ID=6;
domain=Get_Domain(1);
tt=Lookup_Thread(domain,thread_ID);
/*calculate the total number of faces on the wall*/
begin_f_loop(f,tt)
{
c0=F_C0(f,tt);
counter++;
}
end_f_loop(f,tt)
a=(int*)malloc(counter*sizeof(int));
b=(real*)malloc(counter*sizeof(real));
d=(real*)malloc(counter*sizeof(real)); /* temp info*/
q=(real*)malloc(counter*sizeof(real)); /* temp info at previous time*/
g=(real*)malloc(counter*sizeof(real)); /* mass fraction info at previous time*/
l=(real*)malloc(counter*sizeof(real)); /* pressure info at previous time*/
begin_f_loop(f,tt)
{
c0=F_C0(f,tt);
a=c0;
F_AREA(A,f,tt);
b=NV_MAG(A);
i++;
}
end_f_loop(f,tt)
}








DEFINE_ADJUST(k_init,domain)
{
k=0;
}
DEFINE_ADJUST(temp_at_interface,domain)
{
Thread *f_thread;
face_t face;
int thread_ID;
int i=0;
thread_ID=6; /*ID of g_s_wall is 6*/
f_thread=Lookup_Thread(domain,thread_ID);
begin_f_loop(face,f_thread)
{
d=F_T(face,f_thread);
i++;
}
end_f_loop(face,f_thread);
}


DEFINE_ADJUST(avg,domain)
{
int i=0;
real tavg = 0;
real wavg = 0;
real Pavg = 0;
real temp,volume,vol_tot=0,w_h2o,Pa,condmass;
real sum=0;
Domain *d;
Thread *t;
cell_t c;
real x[ND_ND];
d = Get_Domain(1);
thread_loop_c(t,d)
{
begin_c_loop(c,t)
{
C_CENTROID(x,c,t);
volume = C_VOLUME(c,t); /* get cell volume */
temp = C_T(c,t); /* get cell temperature */
w_h2o= C_YI(c,t,i); /* get cell mass fraction of h2o */
Pa=C_P(c,t)+402145.725; /* get cell Absolute Pressure */
condmass=C_UDMI(c,t,0);
if ((-0.001<x[1])&(x[1]<2.001))
{
volume = C_VOLUME(c,t);
temp = C_T(c,t);
w_h2o= C_YI(c,t,i);
Pa=C_P(c,t)+402145.725;
condmass=C_UDMI(c,t,0);
}
else
{
temp =0;
volume = 0;
w_h2o=0;
Pa=0;
condmass=0;
}
sum+=condmass;
vol_tot+=volume;
tavg+=temp*volume;
wavg+=w_h2o*volume;
Pavg+=Pa*volume;
}
end_c_loop(c,t)
tavg /= vol_tot;
wavg /= vol_tot;
Pavg /= vol_tot;
}
s0=sum;
s1=tavg;
s2=wavg;
s3=Pavg;
}


DEFINE_EXECUTE_AT_END(temp_info)
{
Domain *domain;
Thread *f_thread;
Thread *c_thread;
face_t face;
cell_t cell;
int thread_ID;
int i=0;
int mf=0;
thread_ID=6; /*ID OF g_s_wall is 6*/
domain=Get_Domain(1);
f_thread=Lookup_Thread(domain,thread_ID);
begin_f_loop(face,f_thread)
{
q=F_T(face,f_thread);
i++;
}
end_f_loop(face,f_thread)
i=0;
thread_ID=2;/*ID of gas is 2*/
c_thread=Lookup_Thread(domain,thread_ID);
begin_c_loop(cell,c_thread)
{
if(cellTest(cell))
{
g=C_YI(cell,c_thread,mf);
l=C_P(cell,c_thread);
i++;
}
}
end_c_loop(cell,c_thread);
}


DEFINE_SOURCE(vapor_m_source, c, t, dS, eqn)
{
        real source = 0;
        int i; /*index of vapor*/
        real w; /*vapor mass fraction*/
        real w_delt; /*vapor mass fraction at previous time step*/
        real area; /*area of wall boundary*/
        real V; /*volume of neighbouring cell*/
        real Tw; /*temperature of wall, K*/
        real Tw_delt; /*temperature of fluid at previous time step, K*/
        real T; /*temperature of fluid, K*/
        real Hlg; /*latent heat of vapor(J/Kg)*/
        real Cu; /*adjustable coefficient in Uchida correlation*/
        real Tref; /*reference temperature*/
        real Cp_mix = 1500; /*spefic heat of vapor at constant pressure, J/(Kg.K)*/
        real Cp_air; /*spefic heat of air at constant pressure and referenct temperature, J/(Kg.K)*/
        real P; /*pressure of neighbouring cells*/
        real P_delt; /*pressure of neighbouring cells at previous time step, Pa*/
        real Pv; /*pressure of vapor pressure in the neighbouring cells*/
        real Ps; /*saturation pressure of vapor*/
        real hv; /*vapor specific enthalpy, J/Kg*/
        real Vx; /*axial velocity*/
        real Vy; /*radial velocity*/
        Cu = 0.98; /*adjustable coefficient*/
        Tref = 273.15;
        Cp_air = 1005.0;
        i = 0;
        T = s1;
        w = s2;
        P = s3;
        Vx = C_U(c, t);
        Vy = C_V(c, t);
        V = C_VOLUME(c, t);
        if (cellTest(c))
        {
                area = b[j];
                Tw = d[j];
                Tw_delt = q[j];
                w_delt = g[j];
                w = C_YI(c, t, i);
                Cp_mix = C_CP(c, t);
                P_delt = l[j] + 402145.725;
                /*Pv=Ma*w*P/(Ma*w+Mv*(1-w));*/
                Pv = Ma*w_delt*P_delt / (Ma*w_delt + Mv*(1 - w_delt));
                Ps = 2000 * exp(18.5916 - 3991.11 / (Tw_delt - 39.31)) / 15;
                hv = s_enthalpy(Tw_delt);
                if (Pv >= Ps)
                {
                        Hlg = latent_heat(Tw);
                        source = -Cu*(10189.3 + 0.0904164*P - (4314.4 + 0.046537*P)*log10(100 * (1 - w)))*area*pow((T - Tw), 0.4) / (Hlg*V);/*source=-Cu*pow(x/(1-x),0.8)*area*(T-Tw)/(Hlg*V)*/ /*Uchida correlation*/
                        dS[eqn] = 0.;
                }
                else
                {
                        source = 0.;
                }
                C_UDMI(c, t, 0) = source*V;
                C_UDMI(c, t, 1) = source*(T*Cp_mix - Tref*Cp_air)*V; /*C_UDMI(c,t,1)=source*hv;*/
                C_UDMI(c, t, 2) = source*Vx*V; /*moment source in axis direction*/
                C_UDMI(c, t, 3) = source*Vy*V; /*moment source in xradius direction*/
                return source;
        }



        DEFINE_SOURCE(vapor_species_source, c, t, dS, eqn);
        {
                real source = 0;
                real V;
                V = C_VOLUME(c, t);
                source = C_UDMI(c, t, 0) / V;
                dS[eqn] = 0.0;
                return source;
        }


        DEFINE_SOURCE(vapor_h_source, c, t, dS, eqn);
        {
                real source = 0;
                real V;
                V = C_VOLUME(c, t);
                source = C_UDMI(c, t, 1) / V;
                dS[eqn] = 0.0;
                return source;
        }


        DEFINE_SOURCE(Yaxi_mom_source, c, t, dS, eqn);
        {
                real source = 0;
                real V;
                V = C_VOLUME(c, t);
                source = C_UDMI(c, t, 2);
                dS[eqn] = 0.0;
                return source;
        }

        DEFINE_SOURCE(Xrad_mom_source, c, t, dS, eqn);
        {
                real source = 0;
                real V;
                V = C_VOLUME(c, t);
                source = C_UDMI(c, t, 3);
                dS[eqn] = 0.0;
                return source;
        }


        DEFINE_DIFFUSIVITY(vapor_diffusivity, c, t, i); /*vapor diffusivity in the air*/
        {
                real diff = 0;
                real T;
                real P;
                T = C_T(c, t);
                P = C_P(c, t) + 402145.725;
                diff = 4.7931*pow(10.0, -5.0)*pow(T, 1.81) / P;
                return diff;
        }
        int dd;
                int cellText(dd); /*judge if the given cell belongs to the neighbouring cells*/
        {
                int aa = 0;
                int i;
                for (i = 0; i < counter; i++)
                {
                        if (dd == a)
                        {
                                aa = 1;
                                j = i;
                                k++;
                                break;
                        }
                }
                return aa;
        }
}
回复此楼

» 猜你喜欢
1楼 2016-06-22 17:12:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guo718822580

新虫 (初入文坛)
楼主，你好
最近在做蒸发冷凝，
我也是看到这段代码，其中有一些不明：
if((t>=10.0)&&(t<20.0))
Hlg=2477208.8+(t-10.0)*(2453549.6-2477208.8)/10.0;
Hlg是一个什么值，他应用的是什么方程写的呢；
期待您的回答，谢谢~
赞一下 回复此楼
6楼2018-08-14 16:15:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

mycc

专家顾问 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
celltest估计是自己写的宏，来判断这个cell处于什么状态，因此，你这个UDF还不完整
赞一下 回复此楼
项目汇总：http://www.cfluid.com/forum.php?mod=viewthread&amp;tid=114340&amp;extra=
2楼2016-06-22 21:50:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小七工作室1

新虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
月只蓝: 金币-100, 应助指数-2, 屏蔽内容, 违规存档, 多次用完全一样、无实质性帮助的内容，回复不同应助帖！ 2017-03-24 11:15:28
本帖内容被屏蔽
3楼2016-06-24 08:57:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1046418028

新虫 (初入文坛)
你好，你的问题解决了吗，怎样解决的，我最近也在做相变冷凝的课题
赞一下 回复此楼
4楼2017-10-24 17:38:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296材料专硕求调剂 +4 202451007219 2026-04-02 5/250 2026-04-02 09:47 by 晴空210210
[考研] 各位老师好，我的一志愿为北京科技大学085601材料专硕 +13 Koxui 2026-03-28 13/650 2026-04-02 09:35 by 哒哒哒呱呱呱
[考研] 261求B区调剂 +4 明仔· 2026-04-01 6/300 2026-04-02 09:13 by 6781022
[考研] 一志愿北京科技，085601总分305求调剂 +9 半生瓜！ 2026-04-01 11/550 2026-04-02 08:28 by Wang200018
[考研] 一志愿北交大材料工程，总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 362求调剂 +13 西南交材料专硕3 2026-03-31 13/650 2026-04-01 17:38 by JYD2011
[考研] 285求调剂 +5 FZAC123 2026-03-30 5/250 2026-04-01 15:50 by 韩雨涵
[考研] 求调剂 +4 DADA怪 2026-03-31 4/200 2026-04-01 14:30 by ZXlzxl0425
[考研] 0710生物学求调剂！ +6 叙述文 2026-03-31 6/300 2026-04-01 09:39 by JourneyLucky
[考研] 0856材料化工调剂 总分330 +18 zhubinhao 2026-03-27 18/900 2026-04-01 09:37 by oooqiao
[考研] 299求调剂 +8 嗯嗯嗯嗯2 2026-03-27 8/400 2026-03-31 18:23 by lizhi8172
[考研] 材料与化工272求调剂 +25 阿斯蒂芬2004 2026-03-28 25/1250 2026-03-31 16:27 by hypershenger
[考研] 085404 22408 315分 +5 zhuangyan123 2026-03-31 6/300 2026-03-31 13:48 by limeifeng
[考研] 南京大学化学调剂 +11 景随风 2026-03-29 16/800 2026-03-31 10:14 by herarysara
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
[考研] 085701求调剂初试286分 +5 secret0328 2026-03-28 5/250 2026-03-30 12:54 by fangnagu
[考研] 081200-314 +3 LILIQQ 2026-03-27 4/200 2026-03-28 09:41 by 保护地球你我做�
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
信息提示
关闭
请填处理意见
关闭