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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

lwzhao1990

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 2057.2
散金: 610
紅花: 3
帖子: 146
在線: 31.1小時
蟲號: 660539
注冊: 2008-11-23
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

[交流] 【求助】冷凝的UDF程序已有7人參與

哪位仁兄有冷凝的UDF程序？上傳一下，謝謝了！

回復(fù)此樓

1樓 2009-10-29 09:13:02

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

回帖支持 ( 顯示支持度最高的前 50 名 )

tacaro

至尊木蟲 (知名作家)

小木蟲大能貓騎士勛章

仿真EPI: 1
應(yīng)助: 21 (小學(xué)生)
貴賓: 0.06
金幣: 33632.5
散金: 1128
紅花: 17
沙發(fā): 2
帖子: 6301
在線: 1292.2小時
蟲號: 663213
注冊: 2008-11-27
專業(yè): 武術(shù)理論與方法

★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
nono2009(金幣+3,VIP+0):今天表現(xiàn)不錯！ 10-29 12:01

給你一份網(wǎng)上流傳比較多的計算冷凝蒸發(fā)的udf

#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 373
#define LAT_HT 1.e3

DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn)
{

  Thread *mix_th, *sec_th;
  real m_dot_l;

  mix_th = THREAD_SUPER_THREAD(pri_th);
  sec_th = THREAD_SUB_THREAD(mix_th, 1);

  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;

dS[eqn] = -0.1*C_R(cell, pri_th)*
   fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;
                           }
  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.;

   }

  return m_dot_l;

}

DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn)

{

  Thread * mix_th, *pri_th;
  real m_dot_v;

  mix_th = THREAD_SUPER_THREAD(sec_th);
  pri_th = THREAD_SUB_THREAD(mix_th, 0);

  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;
dS[eqn] = 0.;
                           }
  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;

   }

  return m_dot_v;

}

DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)
{
  Thread *pri_th, *sec_th;
  real m_dot;
  pri_th = THREAD_SUB_THREAD(mix_th, 0);
  sec_th = THREAD_SUB_THREAD(mix_th, 1);

  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;

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;}

  return LAT_HT*m_dot;
}

希望能有點啟發(fā)

贊一下(13人)

回復(fù)此樓

2樓2009-10-29 11:53:19

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

nirvanatt

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 72.5
紅花: 1
帖子: 175
在線: 40小時
蟲號: 779354
注冊: 2009-05-25
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
nono2009(金幣+0,VIP+0):都是版上的老蟲子了，活躍一下氣氛。常來交流，獎勵多多。呵呵 10-29 14:55

呵呵斑竹先加分后減分還有獎勵打折這么一說啊。。。。

贊一下(13人)

回復(fù)此樓

6樓2009-10-29 14:51:22

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

lwzhao1990

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 2057.2
散金: 610
紅花: 3
帖子: 146
在線: 31.1小時
蟲號: 660539
注冊: 2008-11-23
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

這個我也有，是計算沸騰的UDF，不知道冷凝的有沒有？

贊一下(1人)

回復(fù)此樓

3樓2009-10-29 13:13:39

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

tacaro

至尊木蟲 (知名作家)

小木蟲大能貓騎士勛章

仿真EPI: 1
應(yīng)助: 21 (小學(xué)生)
貴賓: 0.06
金幣: 33632.5
散金: 1128
紅花: 17
沙發(fā): 2
帖子: 6301
在線: 1292.2小時
蟲號: 663213
注冊: 2008-11-27
專業(yè): 武術(shù)理論與方法

nono2009(金幣+0,VIP+0):你說是，偶也不信。呵呵 10-30 20:07

6樓不是我mj……

贊一下(11人)

回復(fù)此樓

8樓2009-10-30 19:02:53

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

普通回帖

tacaro

至尊木蟲 (知名作家)

小木蟲大能貓騎士勛章

仿真EPI: 1
應(yīng)助: 21 (小學(xué)生)
貴賓: 0.06
金幣: 33632.5
散金: 1128
紅花: 17
沙發(fā): 2
帖子: 6301
在線: 1292.2小時
蟲號: 663213
注冊: 2008-11-27
專業(yè): 武術(shù)理論與方法

小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
nono2009(金幣-1,VIP+0):濫竽充數(shù)，獎勵打折！ 10-29 13:35
kuhailangyu(金幣+0,VIP+0):多多參與，知道你是周伯通類型的了！哈哈 10-29 13:46

都是相變

原理是相同的

改下udf里幾個方程應(yīng)該就差不多了把

贊一下(21人)

回復(fù)此樓

4樓2009-10-29 13:16:03

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

kuhailangyu

榮譽版主 (著名寫手)

小木蟲蛀蟲防疫站站長

仿真EPI: 1
應(yīng)助: 0 (幼兒園)
貴賓: 5.006
金幣: 11615.2
紅花: 11
帖子: 1716
在線: 14.7小時
蟲號: 214037
注冊: 2006-03-09
性別: GG
專業(yè): 能源化工
管轄: 仿真模擬

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
nono2009(金幣+2,VIP+0):謝謝K版提供加強版！ 10-29 13:47

引用回帖:

Originally posted by tacaro at 2009-10-29 13:16:
都是相變

原理是相同的

改下udf里幾個方程應(yīng)該就差不多了把

我還有個帶注解的版本，貌似一樣的！應(yīng)該更能看得懂吧，哈哈

#include "udf.h"  ‘包括常規(guī)宏

#include "sg_mphase.h"  ‘包括體積分?jǐn)?shù)宏CVOF(C,T)

#define T_SAT 373 ‘定義蒸發(fā)溫度100℃

#define LAT_HT 1.e3 ‘定義蒸發(fā)潛熱J/Kg

DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) ‘液相質(zhì)量源項UDF

{

Thread *mix_th, *sec_th;  ‘定義計算區(qū)線指針

  real m_dot_l; ‘定義液相質(zhì)量轉(zhuǎn)移 kg/(m2.s)

  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)  ‘如果液相單元的溫度高于蒸發(fā)溫度，液相向氣相的質(zhì)量質(zhì)量轉(zhuǎn)移

{

  m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*

  fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;

dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;

            ‘ 定義源項對質(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;

‘如果指向混合區(qū)液相的單元溫度小于蒸發(fā)溫度，氣相向液相的質(zhì)量轉(zhuǎn)移，液相得

dS[eqn] = 0.;  ‘由于是氣相向液相轉(zhuǎn)移，所以液相的質(zhì)量源項對質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零

   }

  return m_dot_l;

}

DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) ‘氣相質(zhì)量源項UDF

{

  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) ‘如果混合區(qū)單元的溫度高于蒸發(fā)溫度，液相向氣相的質(zhì)量質(zhì)量轉(zhuǎn)移

{

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;

dS[eqn] = 0.; ‘由于是液相向氣相轉(zhuǎn)移，所以氣相的質(zhì)量源項對來自液相的質(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;

‘如果指向混合區(qū)的單元溫度小于蒸發(fā)溫度，氣相向液相的質(zhì)量轉(zhuǎn)移，氣相失

dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;

由于是氣相向液相轉(zhuǎn)移，所以氣相的質(zhì)量源項對自身的質(zhì)量轉(zhuǎn)移的偏導(dǎo)不為零

   }

  return m_dot_v;

}

DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) ‘混合模型能量源項UDF

{

  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) ‘如果混合區(qū)的單元溫度高于蒸發(fā)溫度。質(zhì)量轉(zhuǎn)移由液相向氣相轉(zhuǎn)移，吸熱，質(zhì)量轉(zhuǎn)移量前有負(fù)號

{

m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*

            fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;

dS[eqn] =  -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;

                           }

  else {                                              相反，氣相向液相轉(zhuǎn)移則放熱

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;}

  return LAT_HT*m_dot;  氣化潛熱與質(zhì)量轉(zhuǎn)移率相乘得能量源項W/m3

}

贊一下(12人)

回復(fù)此樓

有能力有脾氣的人--懷才不遇；有能力沒脾氣的人--春風(fēng)得意；沒能力有脾氣的人--一事無成；沒能力沒脾氣的人--貴人相助。

5樓2009-10-29 13:45:07

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

lwzhao1990

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 2057.2
散金: 610
紅花: 3
帖子: 146
在線: 31.1小時
蟲號: 660539
注冊: 2008-11-23
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

nono2009(金幣+0,VIP+0):如果不懂UDF的話，還是建議你好好學(xué)一學(xué)。不然運行出了問題怎么辦？總不能讓人幫你算吧。從長遠(yuǎn)講，對研究也是絕對有用。 10-30 20:11

高手能不能幫忙改一下，UDF方面真的是不懂啊。。。。。。。。。。。。。

贊一下(10人)

回復(fù)此樓

7樓2009-10-30 18:51:14

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

wangzhuang04

應(yīng)助: (幼兒園)
在線:
蟲號: 992615

非常感謝大家

回復(fù)此樓

9樓2010-04-15 10:11:27

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

千古龍飛地

木蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 4481.2
散金: 200
紅花: 2
帖子: 343
在線: 148.5小時
蟲號: 928610
注冊: 2009-12-15
專業(yè): 細(xì)胞生物學(xué)

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流

看來udf要好好學(xué)習(xí)啊

贊一下(1人)

回復(fù)此樓

10樓2010-05-10 20:04:51

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

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

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

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

[考研] 265求調(diào)劑 +8	梁梁校校 2026-03-19	8/400	2026-03-20 09:08 by 每天只擺一小會
[考研] 085600材料與化工求調(diào)劑 +17	enenenhui 2026-03-13	18/900	2026-03-20 08:55 by xingguangj
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +20	rare12345 2026-03-18	20/1000	2026-03-20 08:42 by 無際的草原
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 296求調(diào)劑 +3	www_q 2026-03-18	6/300	2026-03-19 22:28 by zhq0425
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 324分 085600材料化工求調(diào)劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 346求調(diào)劑[0856] +3	WayneLim327 2026-03-16	6/300	2026-03-19 11:21 by WayneLim327
[考研] 274求調(diào)劑 +6	S.H1 2026-03-18	6/300	2026-03-19 09:34 by 花店相見
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky

24小時熱門版塊排行榜

[交流] 【求助】冷凝的UDF程序 已有7人參與

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

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

[交流] 【求助】冷凝的UDF程序已有7人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: