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

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

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

返回列表

下一頁(yè)

fengyezhi

新蟲(chóng) (小有名氣)

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

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

水沸騰UDF導(dǎo)入后顯示：line 66:parse error,好不容易數(shù)到第66行，還是沒(méi)發(fā)現(xiàn)錯(cuò)誤....還請(qǐng)大家?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)對(duì)質(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)對(duì)質(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)對(duì)來(lái)自液相的質(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)對(duì)自身的質(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ù)號(hào)*/
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ù)此樓

» 收錄本帖的淘帖專(zhuān)輯推薦

多相流UDF

我的收藏

» 猜你喜歡

考研調(diào)劑已經(jīng)有4人回復(fù)
281求調(diào)劑已經(jīng)有4人回復(fù)
0805 316求調(diào)劑已經(jīng)有6人回復(fù)
085601求調(diào)劑總分293英一數(shù)二已經(jīng)有3人回復(fù)
08工學(xué)調(diào)劑已經(jīng)有17人回復(fù)
340求調(diào)劑已經(jīng)有4人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)
食品專(zhuān)碩一志愿雙一流 328 已經(jīng)有4人回復(fù)
材料調(diào)劑已經(jīng)有6人回復(fù)
300分，材料，求調(diào)劑，英一數(shù)二已經(jīng)有5人回復(fù)

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

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

47425804

木蟲(chóng) (正式寫(xiě)手)

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

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

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

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

贊一下(1人)

回復(fù)此樓

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

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

賈曉林

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 2.8
散金: 2
帖子: 24
在線: 4.4小時(shí)
蟲(chóng)號(hào): 929402
注冊(cè): 2009-12-16
性別: GG
專(zhuān)業(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ò)誤或者你的相名稱沒(méi)改過(guò)來(lái)？我沒(méi)發(fā)現(xiàn)錯(cuò)誤呵呵

贊一下(1人)

回復(fù)此樓

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

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

fengyezhi

新蟲(chóng) (小有名氣)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

fengyezhi

新蟲(chóng) (小有名氣)

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

引用回帖:

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

這個(gè)語(yǔ)法錯(cuò)誤我已經(jīng)改過(guò)來(lái)了，謝謝你呵呵。但是還是算不下去，我用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

木蟲(chóng) (正式寫(xiě)手)

仿真EPI: 1
應(yīng)助: 24 (小學(xué)生)
金幣: 631.1
散金: 505
紅花: 27
帖子: 807
在線: 181.2小時(shí)
蟲(chóng)號(hào): 905802
注冊(cè): 2009-11-18
專(zhuān)業(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è)語(yǔ)法錯(cuò)誤我已經(jīng)改過(guò)來(lái)了，謝謝你呵呵。但是還是算不下去，我用mixture模型，每次引入liq質(zhì)量源項(xiàng)時(shí)就會(huì)報(bào)錯(cuò)，計(jì)算不下去。只引入vap源項(xiàng)和能量源項(xiàng)就可以算，不知道為什么。我的有傳熱

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

贊一下(1人)

回復(fù)此樓

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

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

倚天屠龍love

金蟲(chóng) (小有名氣)

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

建議樓主看下韓占忠09年出的那本關(guān)于Fluent的書(shū)的139頁(yè)，里面那個(gè)UDF程序和你的一樣，你對(duì)比一下吧~
另外，我想請(qǐng)教一下樓主這個(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

新蟲(chóng) (初入文壇)

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

樓主的問(wèn)題解決了嗎？我也是做氣液相變換熱的，剛開(kāi)始沒(méi)多久，麻煩問(wèn)一下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)出來(lái)的呀？謝謝啦

贊一下

回復(fù)此樓

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

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

小小鳥(niǎo)110

銅蟲(chóng) (正式寫(xiě)手)

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

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

贊一下

回復(fù)此樓

追求心智的成熟

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

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

gearboy78

禁蟲(chóng) (小有名氣)

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

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

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

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

返回列表

下一頁(yè)

普通表情龍兔虎貓

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

[考研] 食品專(zhuān)碩一志愿雙一流 328 +3	xiaom99 2026-03-21	4/200	2026-03-24 21:20 by lailaisimei
[考研] 材料與化工考研調(diào)劑 +7	孅華 2026-03-22	7/350	2026-03-24 21:04 by greychen00
[考研] 07化學(xué)280分求調(diào)劑 +6	722865 2026-03-23	6/300	2026-03-24 20:58 by allen-yin
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +5	niko- 2026-03-22	5/250	2026-03-23 22:01 by fuyu_
[考研] 化學(xué)308分求調(diào)劑 +3	你好明天你好 2026-03-23	3/150	2026-03-23 20:11 by macy2011
[考研] 一志愿武理材料工程348求調(diào)劑 +6	￣^￣゜汗 2026-03-19	9/450	2026-03-23 19:53 by pswait
[考研] 336求調(diào)劑 +4	收到VS 2026-03-20	4/200	2026-03-23 19:02 by macy2011
[考研] 考研化學(xué)308分求調(diào)劑 +7	你好明天你好 2026-03-23	8/400	2026-03-23 18:39 by macy2011
[論文投稿] 急發(fā)核心期刊論文 +3	賢達(dá)問(wèn)津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹
[考研] 291 求調(diào)劑 +4	化工2026屆畢業(yè)?/a> 2026-03-21	5/250	2026-03-23 16:46 by 化工2026屆畢業(yè)?/a>
[考研] 311求調(diào)劑 +6	冬十三 2026-03-18	6/300	2026-03-22 20:18 by edmund7
[考研] 298求調(diào)劑一志愿211 +3	上岸6666@ 2026-03-20	3/150	2026-03-22 15:50 by ColorlessPI
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 22408 344分求調(diào)劑一志愿華電計(jì)算機(jī)技術(shù) +4	solanXXX 2026-03-20	4/200	2026-03-20 23:49 by alg094825
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 086500 325 求調(diào)劑 +3	領(lǐng)帶小熊 2026-03-19	3/150	2026-03-20 18:38 by 盡舜堯1