版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

songsny

銅蟲 (初入文壇)

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

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

想用UDF中的宏 DEFINE_DPM_HEAT_MASS來模擬液滴對SO2的吸收，就照著例子自己改了改，可是編譯出錯(cuò)，又用原來的例子直接編譯，發(fā)現(xiàn)同樣的錯(cuò)誤，代碼如下：錯(cuò)誤提示為..\..\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
求大神指點(diǎn)怎么修改

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ù)此樓

» 猜你喜歡

284求調(diào)劑已經(jīng)有10人回復(fù)
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
07化學(xué)280分求調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有12人回復(fù)
求材料，環(huán)境專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
335求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
一志愿吉大化學(xué)322求調(diào)劑已經(jīng)有4人回復(fù)
環(huán)境學(xué)碩288求調(diào)劑已經(jīng)有8人回復(fù)
341求調(diào)劑(一志愿湖南大學(xué)070300) 已經(jīng)有6人回復(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小時(shí)
蟲號: 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)
你少了個(gè)參量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小時(shí)
蟲號: 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)
你少了個(gè)參量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小時(shí)
蟲號: 1853296
注冊: 2012-06-09
性別: GG
專業(yè): 工程熱力學(xué)

引用回帖:

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

這行不用，是用來說明Z的含義。主要是這個(gè)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的回帖

3_1415926

金蟲 (小有名氣)

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

引用回帖:

4樓: Originally posted by 3_1415926 at 2016-03-08 10:44:09
這行不用，是用來說明Z的含義。主要是這個(gè)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)就行了...

按照你的UDF，應(yīng)該是改成DEFINE_DPM_HEAT_MASS(multivap,p,Cp,hgas,hvap,cvap_surf,Z,dydt,dzdt)

贊一下

回復(fù)此樓

5樓2016-03-08 10:45:42

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

淺淺笑王子

鐵蟲 (正式寫手)

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

好的，謝謝

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

回復(fù)此樓

6樓2016-03-08 11:09:29

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

songsny

銅蟲 (初入文壇)

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

引用回帖:

5樓: Originally posted by 3_1415926 at 2016-03-08 10:45:42
按照你的UDF，應(yīng)該是改成DEFINE_DPM_HEAT_MASS(multivap,p,Cp,hgas,hvap,cvap_surf,Z,dydt,dzdt)...

謝謝~按照你的提示改過了，編譯可以通過，可是一計(jì)算就提示Error: FLUENT received fatal signal (ACCESS_VIOLATION) ，我懷疑是cell_t c0 = RP_CELL(&p->cCell); /* cell and thread */ Thread *t0 = RP_THREAD(&p->cCell); /* where the particle is in */這兩行指針有問題，可是也不知道怎么修改，正在糾結(jié)中

贊一下

回復(fù)此樓

7樓2016-03-09 10:34:29

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

3_1415926

金蟲 (小有名氣)

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

引用回帖:

7樓: Originally posted by songsny at 2016-03-09 10:34:29
謝謝~按照你的提示改過了，編譯可以通過，可是一計(jì)算就提示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這個(gè)函數(shù)在高版本ansys fluent中有點(diǎn)變化，多了個(gè)變量Z，不太確定其他是否也有改變
另外是，高版本（我這邊是14.5）中
#include "surf.h" /* RP_CELL and RP_THREAD are defined in surf.h */
你加這個(gè)頭文件進(jìn)去試試
或者用修改前的udf在fluent 6.3中計(jì)算（這個(gè)應(yīng)該是可行的）

贊一下

回復(fù)此樓

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

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

songsny

銅蟲 (初入文壇)

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

引用回帖:

8樓: Originally posted by 3_1415926 at 2016-03-09 22:16:46
DPM的UDF我也沒怎么用過
不過我看了一下，你的UDF應(yīng)該是來自fluent 6.3的幫助文檔里面的
之前的問題是DEFINE_DPM_HEAT_MASS這個(gè)函數(shù)在高版本ansys fluent中有點(diǎn)變化，多了個(gè)變量Z，不太確定其他是否也有改變
另 ...

在6.3中可以計(jì)算了，就是算一段時(shí)間后會(huì)溢出，應(yīng)該是我的程序內(nèi)容有問題，我再改一改，謝謝！

贊一下

回復(fù)此樓

9樓2016-03-11 09:24:48

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

wo_不想念

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 248.6
帖子: 22
在線: 19小時(shí)
蟲號: 2132169
注冊: 2012-11-17
性別: GG
專業(yè): 制造系統(tǒng)與自動(dòng)化

請問樓主讀懂這段程序了嗎。。。我也正在研究這段，還沒搞懂

贊一下

回復(fù)此樓

我在幸福里等你。。。

10樓2016-03-11 11:35:18

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 songsny 的主題更新

返回列表

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

[考研] 291 求調(diào)劑 +4	化工2026屆畢業(yè)?/a> 2026-03-21	5/250	2026-03-23 16:46 by 化工2026屆畢業(yè)?/a>
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt
[考研] 北科281學(xué)碩材料求調(diào)劑 +8	tcxiaoxx 2026-03-20	9/450	2026-03-23 12:16 by tcxiaoxx
[考研] 070300，一志愿北航320求調(diào)劑 +3	Jerry0216 2026-03-22	5/250	2026-03-23 09:16 by 。。堂堂
[考研] 318求調(diào)劑 +4	plum李子 2026-03-21	7/350	2026-03-22 14:17 by ColorlessPI
[考研] 354求調(diào)劑 +7	Tyoumou 2026-03-18	10/500	2026-03-22 11:11 by 人來盛
[考研] 0856材料專碩353求調(diào)劑 +4	NIFFFfff 2026-03-20	4/200	2026-03-22 09:49 by 2026paper
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無際的草原
[考研] 求調(diào)劑 +4	要好好無聊 2026-03-21	4/200	2026-03-21 18:57 by 學(xué)員8dgXkO
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 299求調(diào)劑 +6	△小透明* 2026-03-17	6/300	2026-03-21 02:42 by JourneyLucky
[考研] 332求調(diào)劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e