 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1743 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
sunzhijian新蟲(chóng) (初入文壇)
|
[求助]
關(guān)于網(wǎng)上空氣蒸汽冷凝UDF celltest的問(wèn)題 已有2人參與
|
||
|
在網(wǎng)上找了一個(gè)空氣蒸汽冷凝UDF 用于模擬濕空氣冷凝的問(wèn)題，在模擬過(guò)程中出現(xiàn)Cpp1.obj : error LNK2019: 無(wú)法解析的外部符號(hào) celltest，該符號(hào)在函數(shù) temp_info 中被引用，Cpp1.obj : error LNK2019: 無(wú)法解析的外部符號(hào) cell_T，該符號(hào)在函數(shù) vapor_m_source 中被引。不知道這個(gè)celltest代表什么意思，請(qǐng)高手予以幫助，或給個(gè)可用的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; } } |
回復(fù)此樓» 猜你喜歡
290求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
327求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600 286分 材料求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑
已經(jīng)有15人回復(fù)
-
085601材料工程找調(diào)劑
已經(jīng)有15人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有15人回復(fù)
-
南京大學(xué)化學(xué)調(diào)劑
已經(jīng)有14人回復(fù)
-
一志愿北京工業(yè)大學(xué)，324分求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
月只藍(lán): 金幣-100, 應(yīng)助指數(shù)-2, 屏蔽內(nèi)容, 違規(guī)存檔, 多次用完全一樣、無(wú)實(shí)質(zhì)性幫助的內(nèi)容，回復(fù)不同應(yīng)助帖！ 2017-03-24 11:15:28
感謝參與，應(yīng)助指數(shù) +1
月只藍(lán): 金幣-100, 應(yīng)助指數(shù)-2, 屏蔽內(nèi)容, 違規(guī)存檔, 多次用完全一樣、無(wú)實(shí)質(zhì)性幫助的內(nèi)容，回復(fù)不同應(yīng)助帖！ 2017-03-24 11:15:28
|
本帖內(nèi)容被屏蔽 |
3樓2016-06-24 08:57:13
mycc
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +181 - 應(yīng)助: 230 (大學(xué)生)
- 金幣: 2589.6
- 紅花: 35
- 帖子: 531
- 在線(xiàn): 58.8小時(shí)
- 蟲(chóng)號(hào): 881523
- 注冊(cè): 2009-10-23
- 性別: GG
- 專(zhuān)業(yè): 流體力學(xué)
- 管轄: 仿真模擬
2樓2016-06-22 21:50:30
1046418028
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 22.5
- 帖子: 9
- 在線(xiàn): 12.2小時(shí)
- 蟲(chóng)號(hào): 3712240
- 注冊(cè): 2015-03-05
- 專(zhuān)業(yè): 傳熱傳質(zhì)學(xué)
4樓2017-10-24 17:38:51
18283607156
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 182
- 紅花: 1
- 帖子: 24
- 在線(xiàn): 9.4小時(shí)
- 蟲(chóng)號(hào): 7035247
- 注冊(cè): 2017-08-07
- 性別: GG
- 專(zhuān)業(yè): 高等教育學(xué)
5樓2017-12-27 11:25:02
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚(yú) 2026-03-26 | 9/450 |
|
|
|
[考研] 343求調(diào)劑 +5 | 愛(ài)羈絆 2026-03-28 | 5/250 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書(shū)。 2026-03-28 | 3/150 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲(chóng)085600 2026-03-27 | 8/400 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺(jué)_16 2026-03-24 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
|
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5 | SYA！ 2026-03-23 | 5/250 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣(mài)魚(yú) 2026-03-26 | 3/150 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研] 化工專(zhuān)碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|