版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

songsny

銅蟲 (初入文壇)

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

[求助] 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;

}

}

}

回復此樓

» 猜你喜歡

07化學280分求調(diào)劑已經(jīng)有3人回復
考研化學308分求調(diào)劑已經(jīng)有8人回復
070300化學求調(diào)劑已經(jīng)有7人回復
生物學一志愿985，分數(shù)349求調(diào)劑已經(jīng)有9人回復
341求調(diào)劑(一志愿湖南大學070300) 已經(jīng)有4人回復
336化工調(diào)劑已經(jīng)有5人回復
279分求調(diào)劑一志愿211 已經(jīng)有17人回復
291求調(diào)劑已經(jīng)有4人回復
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有12人回復
一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑已經(jīng)有5人回復

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

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

3_1415926

金蟲 (小有名氣)

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

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(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
加上去就沒事了

贊一下

回復此樓

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

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

淺淺笑王子

鐵蟲 (正式寫手)

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

引用回帖:

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客戶端

贊一下

回復此樓

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

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

3_1415926

金蟲 (小有名氣)

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

引用回帖:

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)就行了

贊一下

回復此樓

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

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

3_1415926

金蟲 (小有名氣)

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

引用回帖:

4樓: Originally posted by 3_1415926 at 2016-03-08 10:44:09
這行不用，是用來說明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)就行了...

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

贊一下

回復此樓

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

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

淺淺笑王子

鐵蟲 (正式寫手)

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

好的，謝謝

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

回復此樓

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

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

songsny

銅蟲 (初入文壇)

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

引用回帖:

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

謝謝~按照你的提示改過了，編譯可以通過，可是一計算就提示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 */這兩行指針有問題，可是也不知道怎么修改，正在糾結中

贊一下

回復此樓

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

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

3_1415926

金蟲 (小有名氣)

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

引用回帖:

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應該是來自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中計算（這個應該是可行的）

贊一下

回復此樓

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

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

songsny

銅蟲 (初入文壇)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

wo_不想念

銅蟲 (初入文壇)

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

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

贊一下

回復此樓

我在幸福里等你。。。

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

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

相關版塊跳轉我要訂閱樓主 songsny 的主題更新

返回列表

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

[考研] 279分求調(diào)劑一志愿211 +15	chaojifeixia 2026-03-19	17/850	2026-03-23 18:30 by allen-yin
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +5	jiajunX 2026-03-22	5/250	2026-03-23 18:07 by YMU施老師
[考研] 336求調(diào)劑 +3	收到VS 2026-03-20	3/150	2026-03-23 13:23 by Iveryant
[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 生物學調(diào)劑 +5	Surekei 2026-03-21	5/250	2026-03-22 14:39 by tcx007
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學員8dgXkO
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學員8dgXkO
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	14/700	2026-03-21 13:38 by lature00
[考研] 材料學學碩080502 337求調(diào)劑-一志愿華中科技大學 +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 二本跨考鄭大材料306英一數(shù)二 +3	z1z2z3879 2026-03-17	3/150	2026-03-21 02:29 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 一志愿南昌大學，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 086500 325 求調(diào)劑 +3	領帶小熊 2026-03-19	3/150	2026-03-20 18:38 by 盡舜堯1
[考研] 收復試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷