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

返回列表

lwzhao1990

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 2057.2
散金: 610
紅花: 3
帖子: 146
在線: 31.1小時(shí)
蟲號(hào): 660539
注冊(cè): 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小時(shí)
蟲號(hào): 663213
注冊(cè): 2008-11-27
專業(yè): 武術(shù)理論與方法

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

給你一份網(wǎng)上流傳比較多的計(jì)算冷凝蒸發(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;
}

希望能有點(diǎn)啟發(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小時(shí)
蟲號(hào): 779354
注冊(cè): 2009-05-25
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

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

呵呵斑竹先加分后減分還有獎(jiǎng)勵(lì)打折這么一說啊。。。。

贊一下(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小時(shí)
蟲號(hào): 660539
注冊(cè): 2008-11-23
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

這個(gè)我也有，是計(jì)算沸騰的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小時(shí)
蟲號(hào): 663213
注冊(cè): 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小時(shí)
蟲號(hào): 663213
注冊(cè): 2008-11-27
專業(yè): 武術(shù)理論與方法

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

都是相變

原理是相同的

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

贊一下(21人)

回復(fù)此樓

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

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

kuhailangyu

榮譽(yù)版主 (著名寫手)

小木蟲蛀蟲防疫站站長(zhǎng)

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

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

引用回帖:

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

原理是相同的

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

我還有個(gè)帶注解的版本，貌似一樣的！應(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ì)量源項(xiàng)UDF

{

Thread *mix_th, *sec_th;  ‘定義計(jì)算區(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;

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

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

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) ‘氣相質(zhì)量源項(xiàng)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ì)量源項(xiàng)對(duì)來自液相的質(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ì)量源項(xiàng)對(duì)自身的質(zhì)量轉(zhuǎn)移的偏導(dǎo)不為零

   }

  return m_dot_v;

}

DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) ‘混合模型能量源項(xiàng)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ù)號(hào)

{

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)移率相乘得能量源項(xiàng)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小時(shí)
蟲號(hào): 660539
注冊(cè): 2008-11-23
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

nono2009(金幣+0,VIP+0):如果不懂UDF的話，還是建議你好好學(xué)一學(xué)。不然運(yùn)行出了問題怎么辦？總不能讓人幫你算吧。從長(zhǎng)遠(yuǎn)講，對(duì)研究也是絕對(duì)有用。 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)助: (幼兒園)
在線:
蟲號(hào): 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小時(shí)
蟲號(hào): 928610
注冊(cè): 2009-12-15
專業(yè): 細(xì)胞生物學(xué)

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

看來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 的主題更新

返回列表

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

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

[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	8/400	2026-03-20 11:24 by kkcoco25
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[論文投稿] 申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 0703化學(xué)調(diào)劑 +10	妮妮ninicgb 2026-03-15	14/700	2026-03-19 22:59 by 學(xué)員8dgXkO
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 324分 085600材料化工求調(diào)劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 304求調(diào)劑 +6	司空. 2026-03-18	6/300	2026-03-18 23:03 by 星空星月
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠

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

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

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

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

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

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