| 5 | 1/1 | 返回列表 |
| 查看: 3191 | 回復: 8 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
linyi7445962新蟲 (初入文壇)
|
[求助]
有關(guān)水蒸汽冷凝的UDF 已有1人參與
|
||
|
一個球體模型,往里通入混合工質(zhì) 其中包含水蒸汽 采用了多組分傳輸來確定混合工質(zhì)的組成 需要考慮混合工質(zhì)中水蒸氣的冷凝情況 在網(wǎng)上查了下冷凝的UDF 都是多相流的 所以加入多相流模型 首相為混合工質(zhì) 次相為水(液體) 在邊界條件的fluid設(shè)置中的mixture選項中框選source terms 在engery中導入能量項 在首相中也框選source terms 并且在水蒸氣項中導入vap_src 同理在次相中導入liq_src 設(shè)置好進行迭代 但是出現(xiàn)幾種錯誤 1 能迭代但沒產(chǎn)生質(zhì)量轉(zhuǎn)移 2不能進行迭代 提示錯誤 好像發(fā)散之類的 因為剛開始學習UDF 對于這方面不太理解 對于THREAD_SUPER_THREAD(pri_th)以及THREAD_SUB_THREAD(mix_th,1) 括號中的參數(shù)如何確定 以及指針到底指向哪里不是很清楚 在FLUENT HELP中看了下 也不是太明白 希望能幫忙解答下 或者有這方面經(jīng)驗的能交流下 以下是采用的UDF 由于剛進論壇 金幣數(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ā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學大學 ,080500材料科學與工程學碩,總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +7 | 手機用戶 2026-03-17 | 7/350 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 復試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|