| 3 | 1/1 | 返回列表 |
| 查看: 1367 | 回復(fù): 2 | ||
哇塞骰銅蟲 (小有名氣)
|
[求助]
類似空話模型(cavitation)的udf,求大家?guī)兔纯磳Σ粚? 已有1人參與
|
|
在用fluent做管內(nèi)汽水兩相流。實際系統(tǒng)比較復(fù)雜,所以入口條件不是特別明確。猜測進口是滿管飽和水,流動過程中發(fā)生汽化,然后又氣體又液化成水。但是因為整個流域溫差很小,暫時不考慮汽化潛熱所以汽化模型應(yīng)該類似于空化。因為用vof模型來追蹤氣液相界,所以在使用vof的同時不能使用cavitation模型。就想著自己寫個udf。 這個udf是按照沸騰的udf改編的,并將最后的能量源項去除不考慮,求大家?guī)兔纯从袥]有什么問題。相變參考量用的壓力。主相是液相,次相是汽相。 編譯沒有問題,但是計算總是發(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; } |
|
1,這個網(wǎng)上流傳很廣的udf例子,最多只能算個現(xiàn)象出來,這種方式算出來的相變速率(其實就是傳熱速率)與實際相差甚遠,這跟他的一些常數(shù)取值有關(guān)(0.1, 這個值的取法目前看不出出處)。 2,如果你將你研究對象的相變類比成空化的話,建議參考空化的理論模型編寫udf,幫助文檔里對采用的空化模型的數(shù)學(xué)表達式有詳細的描述。 3,空化模型有相應(yīng)的使用范圍,如空化體積分數(shù)等,超出一定限度不僅計算結(jié)果可能不合理甚至收斂也是很大的問題。 4,新版本里,vof模型下已經(jīng)允許采用空化模型 |
銅蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
平樂樂樂 2026-03-26 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 機械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +16 | RichLi_ 2026-03-25 | 16/800 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
[考研] 274求調(diào)劑 +13 | 顧九笙要謙虛 2026-03-24 | 19/950 |
|
|
[考研] 321求調(diào)劑 +3 | 璞玉~~ 2026-03-25 | 3/150 |
|
|
[考研] 招08考數(shù)學(xué) +8 | laoshidan 2026-03-20 | 17/850 |
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
[考研] 材料學(xué)求調(diào)劑 +6 | Stella_Yao 2026-03-20 | 6/300 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +9 | 苑豆豆 2026-03-20 | 9/450 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[考研] 305分求調(diào)劑(食品工程) +5 | Sxy112 2026-03-21 | 7/350 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +3 | 你好明天你好 2026-03-23 | 3/150 |
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|