版塊導(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í)熱門版塊排行榜

返回列表

fengyezhi

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 282.7
散金: 29
紅花: 5
帖子: 230
在線: 49.2小時(shí)
蟲號: 1053472
注冊: 2010-07-07
專業(yè): 多相流熱物理學(xué)

[求助] 水沸騰UDF總是報(bào)錯(cuò)

水沸騰UDF導(dǎo)入后顯示：line 66:parse error,好不容易數(shù)到第66行，還是沒發(fā)現(xiàn)錯(cuò)誤....還請大家?guī)蛶兔�，看看到底哪里出錯(cuò)了，謝啦
#include "udf.h"       /*包括常規(guī)宏*/
#include "sg_mphase.h" /*包括體積分?jǐn)?shù)宏*/
#define T_SAT 373       /*定義飽和溫度*/
#define LAT_HT 1.e3    /*定義水蒸汽潛熱*/

DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn)  /*定義液相源項(xiàng)*/
{
  Thread *mix_th, *sec_th; /*混合相、第二相定義計(jì)算區(qū)域指針*/
  real m_dot_l;          /*蒸發(fā)冷凝速率  定義液相質(zhì)量轉(zhuǎn)移*/
  mix_th = THREAD_SUPER_THREAD(pri_th); /*指向混合區(qū)的主相即液相的指針*/
  sec_th = THREAD_SUB_THREAD(mix_th, 1); /* 指向單相控制區(qū)的氣相的指針*/
  if(C_T(cell, mix_th)>=T_SAT){
m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
            fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;  /*如果液相單元的溫度高于蒸發(fā)溫度，液相向氣相的質(zhì)量轉(zhuǎn)移*/
dS[eqn] = -0.1*C_R(cell, pri_th)*
   fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*定義源項(xiàng)對質(zhì)量轉(zhuǎn)移偏導(dǎo)*/
                           }
  else {
m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
            fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
dS[eqn] = 0.; /*于是氣相向液相轉(zhuǎn)移，所以液相的質(zhì)量源項(xiàng)對質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零*/
   }

  return m_dot_l;
}

DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /*定義氣相源項(xiàng)*/
{
  Thread * mix_th, *pri_th; /*混合相、第一相*/
  real m_dot_v;
  mix_th = THREAD_SUPER_THREAD(sec_th); /*指向混合區(qū)的第二相即氣相的指針*/
  pri_th = THREAD_SUB_THREAD(mix_th, 0);  /*指向單相控制區(qū)的液相的指針，液相為主相*/
  if(C_T(cell, mix_th)>=T_SAT){
m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
            fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; /*如果混合區(qū)單元的溫度高于蒸發(fā)溫度，液相向氣相的質(zhì)量質(zhì)量轉(zhuǎn)移*/
dS[eqn] = 0.; /*由于是液相向氣相轉(zhuǎn)移，所以氣相的質(zhì)量源項(xiàng)對來自液相的質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零*/
                           }
  else {
m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
            fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
dS[eqn] = -0.1*C_R(cell, sec_th)*
   fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;  /*由于是氣相向液相轉(zhuǎn)移，所以氣相的質(zhì)量源項(xiàng)對自身的質(zhì)量轉(zhuǎn)移的偏導(dǎo)不為零*/
   }

  return m_dot_v;
}

DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) /*混合模型能量源項(xiàng)*/
{
  Thread *pri_th, *sec_th;
  real m_dot;
  pri_th = THREAD_SUB_THREAD(mix_th, 0); /*指向混合區(qū)的液相的指針*/
  sec_th = THREAD_SUB_THREAD(mix_th, 1); /*指向混合區(qū)的氣相的指針*/
  if(C_T(cell, mix_th)>=T_SAT){
m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
            fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;  /*如果混合區(qū)的單元溫度高于蒸發(fā)溫度。質(zhì)量轉(zhuǎn)移由液相向氣相轉(zhuǎn)移，吸熱，質(zhì)量轉(zhuǎn)移量前有負(fù)號*/
dS[eqn] = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
                           }
  else {
m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
            fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
dS[eqn] = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;}  /*相反，氣相向液相轉(zhuǎn)移則放熱*/
  return LAT_HT*m_dot; /*能量大小氣化潛熱與質(zhì)量轉(zhuǎn)移率相乘得能量源項(xiàng)*/

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

多相流UDF

我的收藏

» 猜你喜歡

274求調(diào)劑已經(jīng)有7人回復(fù)
北科281學(xué)碩材料求調(diào)劑已經(jīng)有11人回復(fù)
085600 材料與化工 329分求調(diào)劑已經(jīng)有3人回復(fù)
302求調(diào)劑已經(jīng)有3人回復(fù)
312求調(diào)劑已經(jīng)有5人回復(fù)
網(wǎng)絡(luò)空間安全0839招調(diào)劑已經(jīng)有4人回復(fù)
329求調(diào)劑已經(jīng)有3人回復(fù)
【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額已經(jīng)有8人回復(fù)
085602 289分求調(diào)劑已經(jīng)有7人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有4人回復(fù)

1樓 2011-05-16 16:01:39

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

47425804

木蟲 (正式寫手)

仿真EPI: 1
應(yīng)助: 24 (小學(xué)生)
金幣: 631.1
散金: 505
紅花: 27
帖子: 807
在線: 181.2小時(shí)
蟲號: 905802
注冊: 2009-11-18
專業(yè): 化學(xué)反應(yīng)工程

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

★ ★
xiegangmai(金幣+2): 謝謝應(yīng)助！ 2011-05-16 22:30:34
fengyezhi(金幣+5): 2011-05-18 08:31:37

如果你的程序完全是這么寫的話，那么就是最后面少了一個(gè)大括號~~~

贊一下(1人)

回復(fù)此樓

2樓2011-05-16 17:34:17

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

賈曉林

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 2.8
散金: 2
帖子: 24
在線: 4.4小時(shí)
蟲號: 929402
注冊: 2009-12-16
性別: GG
專業(yè): 流體力學(xué)

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

★
fengyezhi(金幣+5): 2011-05-18 08:29:38
xiegangmai(金幣+1): 謝謝應(yīng)助！ 2011-05-19 21:29:42

是不是你的首相和第二項(xiàng)設(shè)置錯(cuò)誤或者你的相名稱沒改過來？我沒發(fā)現(xiàn)錯(cuò)誤呵呵

贊一下(1人)

回復(fù)此樓

3樓2011-05-17 09:17:59

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

fengyezhi

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 282.7
散金: 29
紅花: 5
帖子: 230
在線: 49.2小時(shí)
蟲號: 1053472
注冊: 2010-07-07
專業(yè): 多相流熱物理學(xué)

引用回帖:

Originally posted by 賈曉林 at 2011-05-17 09:17:59:
是不是你的首相和第二項(xiàng)設(shè)置錯(cuò)誤或者你的相名稱沒改過來？我沒發(fā)現(xiàn)錯(cuò)誤呵呵

我用多相流模型算的，liquid和vapor兩相，但是每次算的時(shí)候liquid殘差都很大，后來就算不下去了....我的模型里有傳熱

贊一下

回復(fù)此樓

4樓2011-05-18 08:31:25

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

fengyezhi

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 282.7
散金: 29
紅花: 5
帖子: 230
在線: 49.2小時(shí)
蟲號: 1053472
注冊: 2010-07-07
專業(yè): 多相流熱物理學(xué)

引用回帖:

Originally posted by 47425804 at 2011-05-16 17:34:17:
如果你的程序完全是這么寫的話，那么就是最后面少了一個(gè)大括號~~~

這個(gè)語法錯(cuò)誤我已經(jīng)改過來了，謝謝你呵呵。但是還是算不下去，我用mixture模型，每次引入liq質(zhì)量源項(xiàng)時(shí)就會(huì)報(bào)錯(cuò)，計(jì)算不下去。只引入vap源項(xiàng)和能量源項(xiàng)就可以算，不知道為什么。我的有傳熱

贊一下

回復(fù)此樓

5樓2011-05-18 08:33:22

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

47425804

木蟲 (正式寫手)

仿真EPI: 1
應(yīng)助: 24 (小學(xué)生)
金幣: 631.1
散金: 505
紅花: 27
帖子: 807
在線: 181.2小時(shí)
蟲號: 905802
注冊: 2009-11-18
專業(yè): 化學(xué)反應(yīng)工程

★ ★
xiegangmai(金幣+2): 鼓勵(lì)討論交流！ 2011-05-19 21:29:24

引用回帖:

Originally posted by fengyezhi at 2011-05-18 08:33:22:
這個(gè)語法錯(cuò)誤我已經(jīng)改過來了，謝謝你呵呵。但是還是算不下去，我用mixture模型，每次引入liq質(zhì)量源項(xiàng)時(shí)就會(huì)報(bào)錯(cuò)，計(jì)算不下去。只引入vap源項(xiàng)和能量源項(xiàng)就可以算，不知道為什么。我的有傳熱

那就是編譯順利通過啦~~
可否詳述下報(bào)錯(cuò)信息。
另外，有相變的情況下殘差就是容易很大，建議適當(dāng)減小時(shí)間步長。

贊一下(1人)

回復(fù)此樓

6樓2011-05-19 09:54:01

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

倚天屠龍love

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1140.9
帖子: 101
在線: 38小時(shí)
蟲號: 1194893
注冊: 2011-01-23
性別: GG
專業(yè): 流體力學(xué)

建議樓主看下韓占忠09年出的那本關(guān)于Fluent的書的139頁，里面那個(gè)UDF程序和你的一樣，你對比一下吧~
另外，我想請教一下樓主這個(gè)UDF能用于VOF模型么？樓主現(xiàn)在的沸騰換熱的計(jì)算進(jìn)行的怎么樣了，我現(xiàn)在也在搞這個(gè)，希望能和樓主交流一下！

贊一下

回復(fù)此樓

7樓2011-10-28 19:28:18

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

yanxiujuan

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 20.1
散金: 10
帖子: 22
在線: 46.5小時(shí)
蟲號: 1502550
注冊: 2011-11-21
專業(yè): 技術(shù)經(jīng)濟(jì)學(xué)

樓主的問題解決了嗎？我也是做氣液相變換熱的，剛開始沒多久，麻煩問一下m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT是怎么導(dǎo)出來的呀？謝謝啦

贊一下

回復(fù)此樓

8樓2012-10-17 18:52:41

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

小小鳥110

銅蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 84.2
紅花: 4
帖子: 416
在線: 67.8小時(shí)
蟲號: 2209408
注冊: 2012-12-27
性別: GG
專業(yè): 多相流熱物理學(xué)

#define LAT_HT 1.e3 /*定義水蒸汽潛熱*/ 請問這一句中的潛熱，為什么是1000呢，怎么查啊

贊一下

回復(fù)此樓

追求心智的成熟

9樓2014-03-25 20:42:32

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

gearboy78

禁蟲 (小有名氣)

本帖內(nèi)容被屏蔽

10樓2014-03-26 08:25:27

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

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

返回列表

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

[考研] 材料與化工考研調(diào)劑 +9	孅華 2026-03-22	9/450	2026-03-25 13:09 by cmz0325
[考研] 284求調(diào)劑 +15	Zhao anqi 2026-03-22	15/750	2026-03-25 12:51 by wht0531
[考研] 081700 調(diào)劑 267分 +10	迷人的哈哈 2026-03-23	10/500	2026-03-25 12:12 by userper
[考研] 材料調(diào)劑 +3	iwinso 2026-03-23	3/150	2026-03-25 11:29 by greychen00
[考研] 299求調(diào)劑 +7	shxchem 2026-03-20	9/450	2026-03-25 10:41 by lbsjt
[考研] 287求調(diào)劑 +10	晨昏線與星海 2026-03-19	11/550	2026-03-25 10:35 by userper
[考研] 311求調(diào)劑 +3	冬十三 2026-03-24	3/150	2026-03-24 21:31 by peike
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +15	吃吃吃才有意義 2026-03-19	16/800	2026-03-24 19:51 by 了了了了。。
[考研] 0854 考研調(diào)劑招生了！AI 方向 +5	pk3725069 2026-03-19	17/850	2026-03-24 17:30 by zhouxuan..
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +5	leaflight 2026-03-24	5/250	2026-03-24 16:25 by laoshidan
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta
[考研] 工科0856求調(diào)劑 +5	沐析汀汀 2026-03-21	5/250	2026-03-23 17:56 by 海瑟薇-
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚.. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +9	niko- 2026-03-19	10/500	2026-03-22 16:08 by ColorlessPI
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 317求調(diào)劑 +5	申子申申 2026-03-19	9/450	2026-03-20 22:26 by JourneyLucky
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky