| 5 | 1/1 | 返回列表 |
| 查看: 3211 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
linyi7445962新蟲 (初入文壇)
|
[求助]
有關(guān)水蒸汽冷凝的UDF 已有1人參與
|
||
|
一個球體模型,往里通入混合工質(zhì) 其中包含水蒸汽 采用了多組分傳輸來確定混合工質(zhì)的組成 需要考慮混合工質(zhì)中水蒸氣的冷凝情況 在網(wǎng)上查了下冷凝的UDF 都是多相流的 所以加入多相流模型 首相為混合工質(zhì) 次相為水(液體) 在邊界條件的fluid設(shè)置中的mixture選項(xiàng)中框選source terms 在engery中導(dǎo)入能量項(xiàng) 在首相中也框選source terms 并且在水蒸氣項(xiàng)中導(dǎo)入vap_src 同理在次相中導(dǎo)入liq_src 設(shè)置好進(jìn)行迭代 但是出現(xiàn)幾種錯誤 1 能迭代但沒產(chǎn)生質(zhì)量轉(zhuǎn)移 2不能進(jìn)行迭代 提示錯誤 好像發(fā)散之類的 因?yàn)閯傞_始學(xué)習(xí)UDF 對于這方面不太理解 對于THREAD_SUPER_THREAD(pri_th)以及THREAD_SUB_THREAD(mix_th,1) 括號中的參數(shù)如何確定 以及指針到底指向哪里不是很清楚 在FLUENT HELP中看了下 也不是太明白 希望能幫忙解答下 或者有這方面經(jīng)驗(yàn)的能交流下 以下是采用的UDF 由于剛進(jìn)論壇 金幣數(shù)太少 懸賞金少請見諒 #include "udf.h" #include "sg_mphase.h" #define T_SAT 373 #define LAT_HT 1.e3 DEFINE_SOURCE(vap_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,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=0; } else { m_dot_l=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn]=-0.1*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; } return m_dot_l; } DEFINE_SOURCE(liq_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,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=-0.1*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; } else { m_dot_v=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn]=0; } 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,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)/T_SAT; } else { m_dot=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT; dS[eqn]=0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)/T_SAT;} return LAT_HT*m_dot; } |
Fluent中UDF |
新蟲 (小有名氣)
新蟲 (初入文壇)
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-30 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 +14 | kikiki7 2026-03-30 | 14/700 |
|
|
[考研] 266分,求材料相關(guān)專業(yè)調(diào)劑 +6 | 哇呼哼呼哼 2026-03-30 | 8/400 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
[有機(jī)交流] 考研調(diào)劑 +8 | watb 2026-03-26 | 8/400 |
|
|
[考研] 342求調(diào)劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
[考研] 一志愿上海理工能源動力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑。。。 +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|