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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

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;

}

}

}

回復此樓

» 猜你喜歡

南京大學化學調(diào)劑已經(jīng)有14人回復
一志愿北京工業(yè)大學，324分求調(diào)劑已經(jīng)有6人回復
一志愿河北工業(yè)大學0817化工278分求調(diào)劑已經(jīng)有14人回復
求調(diào)劑已經(jīng)有5人回復
281求調(diào)劑已經(jīng)有6人回復
26申博自薦已經(jīng)有6人回復
一志愿哈爾濱工業(yè)大學材料與化工方向336分已經(jīng)有10人回復
求調(diào)劑有機化學考研356分已經(jīng)有5人回復
300求調(diào)劑，材料科學英一數(shù)二已經(jīng)有10人回復
材料與化工328分調(diào)劑已經(jīng)有7人回復

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 */這兩行指針有問題，可是也不知道怎么修改，正在糾結(jié)中

贊一下

回復此樓

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的回帖

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

返回列表

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

[考研] 一志愿華東師范大學有機化學專業(yè)，初試351分，復試被刷求調(diào)劑! +5	真名有冰 2026-03-29	6/300	2026-03-29 20:53 by 唐沐兒
[考研] 337求調(diào)劑 +4	《樹》 2026-03-29	4/200	2026-03-29 20:52 by 唐沐兒
[考研] 291求調(diào)劑 +7	Y-cap 2026-03-29	8/400	2026-03-29 19:53 by klasasda
[考研] 一志愿北京化工大學材料與化工（085600）296求調(diào)劑 +17	稻妻小編 2026-03-26	17/850	2026-03-29 19:07 by 學員8dgXkO
[考研] 一志愿：西北大學，英一數(shù)一408-284分求調(diào)劑 +4	12.27 2026-03-27	4/200	2026-03-29 14:40 by zhshch
[考研] 332求調(diào)劑 +8	蕉蕉123 2026-03-28	8/400	2026-03-29 10:46 by 周梓丹
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6	芝士雪baoo 2026-03-28	8/400	2026-03-29 08:16 by 松花缸1201
[考研] 081200-11408-276學碩求調(diào)劑 +6	崔wj 2026-03-26	6/300	2026-03-29 01:11 by hanserlol
[考研] 本科新能源科學與工程，一志愿華理能動285求調(diào)劑 +7	AZMK 2026-03-28	11/550	2026-03-28 21:01 by xxxsssccc
[考研] 壓國家一區(qū)線，求導師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 085600，材料與化工321分求調(diào)劑 +9	大饞小子 2026-03-28	9/450	2026-03-28 14:56 by 神馬都不懂
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 315分求調(diào)劑 +7	26考研上岸版26 2026-03-26	7/350	2026-03-28 04:05 by fmesaito
[考研] 0703化學求調(diào)劑，各位老師看看我�。�！ +5	祁祺祺 2026-03-25	5/250	2026-03-27 21:44 by 東方豬豬
[考研] 求調(diào)劑 +4	零八# 2026-03-27	4/200	2026-03-27 18:07 by yu221
[考研] 考研調(diào)劑 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 085600材料與化工調(diào)劑 +10	A-哆啦Z夢 2026-03-23	16/800	2026-03-27 15:13 by caszguilin
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當當
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito