[求助] 類似空話模型（cavitation）的udf，求大家?guī)兔纯磳Σ粚? 已有1人參與

在用fluent做管內(nèi)汽水兩相流。實(shí)際系統(tǒng)比較復(fù)雜，所以入口條件不是特別明確。猜測進(jìn)口是滿管飽和水，流動過程中發(fā)生汽化，然后又氣體又液化成水。但是因?yàn)檎麄€流域溫差很小，暫時不考慮汽化潛熱所以汽化模型應(yīng)該類似于空化。因?yàn)橛胿of模型來追蹤氣液相界，所以在使用vof的同時不能使用cavitation模型。就想著自己寫個udf。
這個udf是按照沸騰的udf改編的，并將最后的能量源項(xiàng)去除不考慮，求大家?guī)兔纯从袥]有什么問題。相變參考量用的壓力。主相是液相，次相是汽相。
編譯沒有問題，但是計(jì)算總是發(fā)散。求大家?guī)兔�。感激不盡。

#include "udf.h"
#include "sg_mphase.h"
#define P_SAT 12826

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); “混合區(qū)主相即液相”
  sec_th = THREAD_SUB_THREAD(mix_th, 1); “控制區(qū)次相即汽相”
  if(C_P(cell, mix_th)<P_SAT)          “混合區(qū)單元壓力小于汽化壓力，則液相發(fā)生汽化”
{
  m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
  fabs(C_P(cell, mix_th) -P_SAT)/P_SAT;
dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_P(cell, mix_th) - P_SAT)/P_SAT;
                              }
  else {                                           “混合區(qū)單元壓力大于汽化壓力，則汽相液化”
m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
            fabs(P_SAT-C_P(cell,mix_th))/P_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);    “混合區(qū)次相即汽相”
  pri_th = THREAD_SUB_THREAD(mix_th, 0);    “控制區(qū)主相即液相”
  if(C_P(cell, mix_th)<T_SAT)                   “混合區(qū)單元壓力小于汽化壓力，則液相汽化”
{
m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
            fabs(C_P(cell, mix_th) - P_SAT)/P_SAT;
dS[eqn] = 0.;
                           }
  else {                                                    “混合區(qū)單元壓力大于汽化壓力，則汽相液化”
m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
            fabs(P_SAT-C_P(cell,mix_th))/P_SAT;
dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_P(cell, mix_th) - P_SAT)/P_SAT;
   }
  return m_dot_v;
}

回復(fù)此樓

» 猜你喜歡

275求調(diào)劑已經(jīng)有10人回復(fù)
286求調(diào)劑已經(jīng)有6人回復(fù)
083000學(xué)碩274求調(diào)劑已經(jīng)有5人回復(fù)
【求調(diào)劑】085601材料工程專碩 | 總分272 | 已經(jīng)有5人回復(fù)
275求調(diào)劑已經(jīng)有13人回復(fù)
0856，材料與化工321分求調(diào)劑已經(jīng)有8人回復(fù)
085701環(huán)境工程，267求調(diào)劑已經(jīng)有12人回復(fù)
320分，材料與化工專業(yè)，求調(diào)劑已經(jīng)有7人回復(fù)
085701環(huán)境工程求調(diào)劑已經(jīng)有8人回復(fù)
295材料工程專碩求調(diào)劑已經(jīng)有5人回復(fù)

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

udf不同壁面熱流求助已經(jīng)有14人回復(fù)
UDF并行編譯求助-附程序段已經(jīng)有3人回復(fù)
[已完結(jié)]UDF計(jì)算species 模型出錯已經(jīng)有3人回復(fù)
fluent14編譯UDF時總是出錯，求高手�。。�！難道沒人遇到過嗎？？？已經(jīng)有11人回復(fù)
fluent 中二維模型中，材料屬性隨著溫度變化的udf源代碼？已經(jīng)有14人回復(fù)
急求：CSF表面張力模型的UDF代碼已經(jīng)有6人回復(fù)
fluent udf的應(yīng)用已經(jīng)有20人回復(fù)
求助！請問怎么用UDF控制mixture模型第二相進(jìn)口的volume fraction 已經(jīng)有6人回復(fù)
全部身家求大神幫忙解決fluent壓力脈沖怎么編寫UDF的問題已經(jīng)有8人回復(fù)
Fluent關(guān)于DPM模型UDF編程的求助已經(jīng)有4人回復(fù)
如何在udf里面判斷界面的位置�。╲of模型）已經(jīng)有9人回復(fù)
求助udf的輸入輸出已經(jīng)有16人回復(fù)
vof模型算氣液兩相逆流使用udf加動量源項(xiàng) 碰到問題請高手解答已經(jīng)有14人回復(fù)
求助！��！我想要在fluent中針對一個面進(jìn)行積分，那么UDF應(yīng)該用什么宏函數(shù)？怎么用？已經(jīng)有5人回復(fù)
求助：UDF模型的編譯問題已經(jīng)有15人回復(fù)
udf中的平方怎么寫已經(jīng)有5人回復(fù)
關(guān)于linux系統(tǒng)下使用UDF的問題已經(jīng)有7人回復(fù)
【求助】通氣空泡多相流（氣、水蒸汽、水）空化模型的UDF如何入手��？？？？已經(jīng)有4人回復(fù)
【求助】求UDF高手幫忙看看我的問題已經(jīng)有18人回復(fù)

1樓 2014-04-11 21:20:35

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

gyctju

金蟲 (正式寫手)

應(yīng)助: 314 (大學(xué)生)
金幣: 942.6
散金: 227
紅花: 48
沙發(fā): 3
帖子: 996
在線: 427.6小時
蟲號: 1569491
注冊: 2012-01-08
專業(yè): 分離過程

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
哇塞骰: 金幣+10, ★★★★★最佳答案, 先給十個，謝謝回答 2014-04-12 10:02:56

1，這個網(wǎng)上流傳很廣的udf例子，最多只能算個現(xiàn)象出來，這種方式算出來的相變速率（其實(shí)就是傳熱速率）與實(shí)際相差甚遠(yuǎn)，這跟他的一些常數(shù)取值有關(guān)（0.1，這個值的取法目前看不出出處）。
2，如果你將你研究對象的相變類比成空化的話，建議參考空化的理論模型編寫udf，幫助文檔里對采用的空化模型的數(shù)學(xué)表達(dá)式有詳細(xì)的描述。
3，空化模型有相應(yīng)的使用范圍，如空化體積分?jǐn)?shù)等，超出一定限度不僅計(jì)算結(jié)果可能不合理甚至收斂也是很大的問題。
4，新版本里，vof模型下已經(jīng)允許采用空化模型

贊一下

回復(fù)此樓

2樓2014-04-11 22:35:53

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

哇塞骰

銅蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 134.8
散金: 25
帖子: 152
在線: 43.2小時
蟲號: 2293475
注冊: 2013-02-22
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

2樓: Originally posted by gyctju at 2014-04-11 22:35:53
1，這個網(wǎng)上流傳很廣的udf例子，最多只能算個現(xiàn)象出來，這種方式算出來的相變速率（其實(shí)就是傳熱速率）與實(shí)際相差甚遠(yuǎn)，這跟他的一些常數(shù)取值有關(guān)（0.1，這個值的取法目前看不出出處）。
2，如果你將你研究對象的 ...

首先很感謝你這么詳細(xì)的答復(fù)，真的非常感謝。
對于0.1的傳熱系數(shù)，之前也見到過論壇的帖子，說這個0.1是不是靠譜不好說，據(jù)說取幾百設(shè)置上萬的。這個我準(zhǔn)備根據(jù)模擬效果改的。
之前看過幫助里關(guān)于空化是數(shù)學(xué)描述，但是看到不是很詳細(xì)，我想有必要再看看。
然后就是你說的新版本里vof可以使用空化了，因?yàn)橹翱吹臅①Y料什么的都是6.3的，都說用空化要在混合多相流模型中使用，我現(xiàn)在用的是14.5，不知道能否同時使用，還請指點(diǎn)一二。
再次感謝

贊一下(1人)

回復(fù)此樓

3樓2014-04-12 10:08:52

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

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

[考研] 265求調(diào)劑11408 +3	劉小鹿lu 2026-03-27	3/150	2026-03-27 20:53 by nihaoar
[考研] 291求調(diào)劑 +6	孅華 2026-03-22	6/300	2026-03-27 19:56 by WYUMater
[考研] 315分求調(diào)劑 +6	26考研上岸版26 2026-03-26	6/300	2026-03-27 19:54 by WYUMater
[考研] 295求調(diào)劑 +4	1428151015 2026-03-27	5/250	2026-03-27 19:53 by JourneyLucky
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 07化學(xué)280分求調(diào)劑 +10	722865 2026-03-23	10/500	2026-03-27 15:51 by Plutoqq
[考博] 26申博 +3	加油沖��！ 2026-03-26	3/150	2026-03-27 15:38 by cls512
[考研] 269專碩求調(diào)劑 +10	金恩貝 2026-03-21	10/500	2026-03-27 15:10 by caszguilin
[考研] 305求調(diào)劑 +5	哇盧卡庫 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 339求調(diào)劑 +4	烤麥芽 2026-03-27	5/250	2026-03-27 13:23 by 752105528
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚的貓
[考研] 一志愿河工大 081700 276求調(diào)劑 +4	地球繞著太陽轉(zhuǎn) 2026-03-23	4/200	2026-03-26 14:27 by zzll406
[考研] 材料調(diào)劑 +3	iwinso 2026-03-23	3/150	2026-03-25 11:29 by greychen00
[考研] 086003食品工程求調(diào)劑 +6	淼淼111 2026-03-24	6/300	2026-03-25 10:29 by 3Strings
[考研] 調(diào)劑 +4	13853210211 2026-03-24	4/200	2026-03-24 19:44 by ms629
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO