songsny

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 112
帖子: 39
在線: 28.2小時
蟲號: 2381059
注冊: 2013-03-27
性別: MM
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

[求助] UDF編譯出錯 DEFINE_DPM_HEAT_MASS 已有2人參與

想用UDF中的宏 DEFINE_DPM_HEAT_MASS來模擬液滴對SO2的吸收，就照著例子自己改了改，可是編譯出錯，又用原來的例子直接編譯，發(fā)現(xiàn)同樣的錯誤，代碼如下：錯誤提示為..\..\src\multivap.c(6) : warning C4003: not enough actual parameters for macro 'DEFINE_DPM_HEAT_MASS'
..\..\src\multivap.c(7) : error C2055: expected formal parameter list, not a type list
..\..\src\multivap.c(25) : error C2109: subscript requires array or pointer type
..\..\src\multivap.c(25) : error C2106: '+=' : left operand must be l-value
..\..\src\multivap.c(26) : error C2223: left of '->energy' must point to struct/union
..\..\src\multivap.c(65) : error C2109: subscript requires array or pointer type
..\..\src\multivap.c(65) : error C2106: '-=' : left operand must be l-value
..\..\src\multivap.c(66) : error C2223: left of '->species' must point to struct/union
..\..\src\multivap.c(68) : error C2109: subscript requires array or pointer type
..\..\src\multivap.c(68) : error C2106: '-=' : left operand must be l-value
..\..\src\multivap.c(70) : error C2223: left of '->energy' must point to struct/union
求大神指點怎么修改

CODE:

/***********************************************************************

UDF for defining the heat and mass transport for

multicomponent particle vaporization

***********************************************************************/

#include "udf.h"

DEFINE_DPM_HEAT_MASS(multivap,p,Cp,hgas,hvap,cvap_surf,dydt,dzdt)

{

int ns;

int nc = TP_N_COMPONENTS( p ); /* number of particle components */

cell_t c0 = RP_CELL(&p->cCell); /* cell and thread */

Thread *t0 = RP_THREAD(&p->cCell); /* where the particle is in */

Material *gas_mix = THREAD_MATERIAL( t0 ); /* gas mixture material */

Material *cond_mix = p->injection->material;/* particle mixture material */

cphase_state_t *c = &(p->cphase); /* cell info of particle location */

real molwt[MAX_SPE_EQNS]; /* molecular weight of gas species */

real Tp = P_T(p); /* particle temperature */

real mp = P_MASS(p); /* particle mass */

real molwt_bulk = 0.; /* average molecular weight in bulk gas */

real Dp = DPM_DIAM_FROM_VOL( mp / P_RHO(p) ); /* particle diameter */

real Ap = DPM_AREA(Dp); /* particle surface */

real Pr = c->sHeat * c->mu / c->tCond; /* Prandtl number */

real Nu = 2.0 + 0.6 * sqrt( p->Re ) * pow( Pr, 1./3. ); /* Nusselt number */

real h = Nu * c->tCond / Dp; /* Heat transfer coefficient */

real dh_dt = h * ( c->temp - Tp ) * Ap; /* heat source term */

dydt[0] += dh_dt / ( mp * Cp );

dzdt->energy -= dh_dt;

{

Material *sp;

mixture_species_loop(gas_mix,sp,ns)

{

molwt[ns] = MATERIAL_PROP(sp,PROP_mwi); /* molecular weight of gas species */

molwt_bulk += C_YI(c0,t0,ns) / molwt[ns]; /* average molecular weight */

}

}

/* prevent division by zero */

molwt_bulk = MAX(molwt_bulk,DPM_SMALL);

for( ns = 0; ns < nc; ns++ )

{

/* gas species index of vaporization */

int gas_index = TP_COMPONENT_INDEX_I(p,ns);

if( gas_index >= 0 )

{

/* condensed material */

Material * cond_c = MIXTURE_COMPONENT(cond_mix, ns );

/* vaporization temperature */

real vap_temp = MATERIAL_PROP(cond_c,PROP_vap_temp);

/* diffusion coefficient */

real D = MATERIAL_PROP_POLYNOMIAL( cond_c, PROP_binary_diffusivity, c->temp);

/* Schmidt number */

real Sc = c->mu / ( c->rho * D );

/* mass transfer coefficient */

real k = ( 2. + 0.6 * sqrt(p->Re) * pow( Sc, 1./3. ) ) * D / Dp;

/* bulk gas concentration */

real cvap_bulk = c->pressure / UNIVERSAL_GAS_CONSTANT / c->temp

* c->yi[gas_index] / molwt_bulk / solver_par.molWeight[gas_index];

/* vaporization rate */

real vap_rate = k * molwt[gas_index] * Ap * ( cvap_surf[ns] - cvap_bulk );

/* only condensation below vaporization temperature */

if( 0. < vap_rate && Tp < vap_temp )

vap_rate = 0.;

dydt[1+ns] -= vap_rate;

dzdt->species[gas_index] += vap_rate;

/* dT/dt = dh/dt / (m Cp)*/

dydt[0] -= hvap[gas_index] * vap_rate / ( mp * Cp );

/* gas enthalpy source term */

dzdt->energy += hgas[gas_index] * vap_rate;

}

}

}

回復(fù)此樓

» 猜你喜歡

材料專業(yè)求調(diào)劑已經(jīng)有7人回復(fù)
材料專碩英一數(shù)二306 已經(jīng)有8人回復(fù)
一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑已經(jīng)有6人回復(fù)
請教下大家 2026年國家基金申請是雙盲審嗎？已經(jīng)有4人回復(fù)
一志愿南航材料專317分求調(diào)劑已經(jīng)有3人回復(fù)
環(huán)境學(xué)碩288求調(diào)劑已經(jīng)有7人回復(fù)
一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 已經(jīng)有4人回復(fù)
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑已經(jīng)有8人回復(fù)
化學(xué)308分求調(diào)劑已經(jīng)有3人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有11人回復(fù)

1樓 2016-03-04 16:35:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

3_1415926

金蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: 1994.3
紅花: 6
帖子: 159
在線: 57.3小時
蟲號: 1853296
注冊: 2012-06-09
性別: GG
專業(yè): 工程熱力學(xué)

引用回帖:

7樓: Originally posted by songsny at 2016-03-09 10:34:29
謝謝~按照你的提示改過了，編譯可以通過，可是一計算就提示Error: FLUENT received fatal signal (ACCESS_VIOLATION) ，我懷疑是cell_t c0 = RP_CELL(&p->cCell); /* cell and thread */ Thread *t0 = ...

DPM的UDF我也沒怎么用過
不過我看了一下，你的UDF應(yīng)該是來自fluent 6.3的幫助文檔里面的
之前的問題是DEFINE_DPM_HEAT_MASS這個函數(shù)在高版本ansys fluent中有點變化，多了個變量Z，不太確定其他是否也有改變
另外是，高版本（我這邊是14.5）中
#include "surf.h" /* RP_CELL and RP_THREAD are defined in surf.h */
你加這個頭文件進去試試
或者用修改前的udf在fluent 6.3中計算（這個應(yīng)該是可行的）

贊一下

回復(fù)此樓

8樓2016-03-09 22:16:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 13 個回答

3_1415926

金蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: 1994.3
紅花: 6
帖子: 159
在線: 57.3小時
蟲號: 1853296
注冊: 2012-06-09
性別: GG
專業(yè): 工程熱力學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
songsny: 金幣+10, ★★★★★最佳答案 2016-03-09 10:37:33

DEFINE_DPM_HEAT_MASS (name, p, C_p, hgas, hvap, cvap_surf, Z, dydt, dzdt)
你少了個參量Z
real Z， Compressibility
加上去就沒事了

贊一下

回復(fù)此樓

2樓2016-03-04 20:53:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

淺淺笑王子

鐵蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 59.5
散金: 527
紅花: 4
帖子: 370
在線: 215.2小時
蟲號: 3468509
注冊: 2014-10-11
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

2樓: Originally posted by 3_1415926 at 2016-03-04 20:53:34
DEFINE_DPM_HEAT_MASS (name, p, C_p, hgas, hvap, cvap_surf, Z, dydt, dzdt)
你少了個參量Z
real Z， Compressibility
加上去就沒事了

real Z，??Compressibility
這一行也寫在下面么

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

3樓2016-03-07 19:17:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

3_1415926

金蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: 1994.3
紅花: 6
帖子: 159
在線: 57.3小時
蟲號: 1853296
注冊: 2012-06-09
性別: GG
專業(yè): 工程熱力學(xué)

引用回帖:

3樓: Originally posted by 淺淺笑王子 at 2016-03-07 19:17:10
real Z，??Compressibility
這一行也寫在下面么
...

這行不用，是用來說明Z的含義。主要是這個DEFINE_DPM_HEAT_MASS(multivap,p,Cp,hgas,hvap,cvap_surf,dydt,dzdt)少了參量Z，改成DEFINE_DPM_HEAT_MASS (name, p, C_p, hgas, hvap, cvap_surf, Z, dydt, dzdt)就行了

贊一下

回復(fù)此樓

4樓2016-03-08 10:44:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 材料專業(yè)求調(diào)劑 +7	hanamiko 2026-03-18	7/350	2026-03-23 20:50 by baobaoye
[考研] 一志愿南航材料專317分求調(diào)劑 +3	炸呀炸呀炸薯條 2026-03-23	3/150	2026-03-23 20:47 by pswait
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte
[考研] 333求調(diào)劑 +6	87639 2026-03-21	10/500	2026-03-23 10:41 by Iveryant
[考研] 352求調(diào)劑 +3	大米飯！ 2026-03-22	3/150	2026-03-22 23:28 by king123！
[考研] 317求調(diào)劑 +12	申子申申 2026-03-19	18/900	2026-03-22 22:23 by luoyongfeng
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋懷
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	4/200	2026-03-22 11:04 by 搏擊518
[基金申請] 山東省面上項目限額評審 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3	李軼男003 2026-03-20	3/150	2026-03-21 13:28 by 搏擊518
[考研] 346求調(diào)劑[0856] +4	WayneLim327 2026-03-16	7/350	2026-03-21 04:02 by JourneyLucky
[考研] 機械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11	adaie 2026-03-19	11/550	2026-03-20 21:30 by laoshidan
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh