 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2166 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
wangyuxuan05銀蟲 (初入文壇)
|
[求助]
液膜表面蒸發(fā)fluent求助！
|
||
|
我的課題是用FLUENT做分子蒸餾器的液膜流體傳質(zhì)傳熱過程模擬，模擬物系是DHA、油酸（oa）、花生四烯酸（ara），所選的幾何區(qū)域是高200mm，寬1mm的豎直液膜，操作壓力為1pa，上端進(jìn)口為mass-flow-inlet，左側(cè)為蒸發(fā)壁面，下端為壓力出口。目前遇到的問題是： 1、右側(cè)液膜表面的邊界條件不知道要如何選取 2、液膜表面涉及到多組分的蒸發(fā)，所以我想在這個(gè)邊界上加入一個(gè)udf 1）我只做縱向上的一維研究，把組分的蒸發(fā)速率編成udf，但是在導(dǎo)入fluent時(shí)總是有錯(cuò)誤，您能把我看一下程序嗎？謝謝！ #include "udf.h" DEFINE_SOURCE(evapor, cell, thread,j) { real nc=80000; real zl=3; real p; real m_e; real M[3]={328.49,282.46,302.45}; real T; real T0=573; real fm; real p0[3]={3200,2800,3000}; real x0[3]={0.2158,0.6268,0.1574}; real interval; real x; if (interval==1) { for (j=0;j<3;j++) { x=x0[j]; } T[cell]=T0; } else { for (i=1;i<=zl;i++) { fm=0; fm=fm+C_YI(cell,thread,i)/M[i-1]; } for (j=1;j<=zl;j++) { x[cell,j]=C_YI(cell,thread,j)/M[j-1]/fm; } T[cell]=C_T(cell,thread); } p[cell,j]=p0[j-1]; m_e(cell,thread,j)=0.9*15.8*x[cell,j]*P[cell,j]*M[j]/T[cell]; return m_e(cell,thread,i,j); } 2）我也嘗試著將論壇里水蒸發(fā)的udf修改成多組分的蒸發(fā)，但是成功導(dǎo)入U(xiǎn)DF后，fluent不能運(yùn)行，您能幫我把這個(gè)udf也看看嗎？謝謝！ #include "udf.h" #include "sg_mphase.h" #define T_SAT_dha 593 #define T_SAT_ara 533 #define T_SAT_oa 523 #define LAT_HT_dha 237425.8 #define LAT_HT_ara 285632.3 #define LAT_HT_oa 280781.3 DEFINE_SOURCE(liq_src_dha, cell, pri_th, dS, eqn) { Thread *mix_th, *sec_th; real m_dot_l_dha; mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th, 1); if(C_T(cell, mix_th)>=T_SAT_dha) { m_dot_l_dha = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT_dha)/T_SAT_dha; dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT_dha)/T_SAT_dha; } else {m_dot_l_dha= 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT_dha-C_T(cell,mix_th))/T_SAT_dha; dS[eqn] = 0.; } return m_dot_l_dha; } DEFINE_SOURCE(vap_src_dha, cell, sec_th, dS, eqn) { Thread * mix_th, *pri_th; real m_dot_v_dha; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th, 0); if(C_T(cell, mix_th)>=T_SAT_dha) {m_dot_v_dha = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT_dha)/T_SAT_dha; dS[eqn] = 0.; } else {m_dot_v_dha= -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT_dha-C_T(cell,mix_th))/T_SAT_dha; dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT_dha)/T_SAT_dha; } return m_dot_v_dha; } DEFINE_SOURCE(liq_src_ara, cell, pri_th, dS, eqn) { Thread *mix_th, *sec_th; real m_dot_l_ara; mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th, 1); if(C_T(cell, mix_th)>=T_SAT_ara) { m_dot_l_ara = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT_ara)/T_SAT_ara; dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT_ara)/T_SAT_ara; } else {m_dot_l_ara = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT_ara-C_T(cell,mix_th))/T_SAT_ara; dS[eqn] = 0.; } return m_dot_l_ara; } DEFINE_SOURCE(vap_src_ara, cell, sec_th, dS, eqn) { Thread * mix_th, *pri_th; real m_dot_v_ara; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th, 0); if(C_T(cell, mix_th)>=T_SAT_ara) {m_dot_v_ara = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT_ara)/T_SAT_ara; dS[eqn] = 0.; } else {m_dot_v_ara = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT_ara-C_T(cell,mix_th))/T_SAT_ara; dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT_ara)/T_SAT_ara;} return m_dot_v_ara; } DEFINE_SOURCE(liq_src_oa, cell, pri_th, dS, eqn) { Thread *mix_th, *sec_th; real m_dot_l_oa; mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th, 1); if(C_T(cell, mix_th)>=T_SAT_oa) { m_dot_l_oa= -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT_oa)/T_SAT_oa; dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT_oa)/T_SAT_oa; } else {m_dot_l_oa = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT_oa-C_T(cell,mix_th))/T_SAT_oa; dS[eqn] = 0.; } return m_dot_l_oa ; } DEFINE_SOURCE(vap_src_oa, cell, sec_th, dS, eqn) { Thread * mix_th, *pri_th; real m_dot_v_oa; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th, 0); if(C_T(cell, mix_th)>=T_SAT_oa) {m_dot_v_oa = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT_oa)/T_SAT_oa; dS[eqn] = 0.; } else {m_dot_v_oa = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT_oa-C_T(cell,mix_th))/T_SAT_oa; dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT_oa)/T_SAT_oa; } return m_dot_v_oa; |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 博士科研新路歷程 |
» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有69人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 新手求助聚砜鑄膜液表面遇水生成白色薄膜的問題
已經(jīng)有6人回復(fù)
- 【求助】氣液分離兩相流動(dòng)模擬fluent gambit出入口的設(shè)置
已經(jīng)有11人回復(fù)
- 【求助】用fluent模擬超重力氣液流體動(dòng)力模型和傳質(zhì)過程，中間有化學(xué)吸收反應(yīng)，fluent
已經(jīng)有23人回復(fù)
- 【求助】Fluent中的對(duì)稱問題
已經(jīng)有10人回復(fù)
- 【求助】Fluent 收斂問題
已經(jīng)有10人回復(fù)
|
本帖內(nèi)容被屏蔽 |
3樓2019-02-11 08:34:28
ggpass
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 148.8
- 帖子: 52
- 在線: 27.4小時(shí)
- 蟲號(hào): 4404894
- 注冊(cè): 2016-02-16
- 專業(yè): 化學(xué)工程及工業(yè)化學(xué)
2樓2019-01-10 11:17:57
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +7 | 麻辣魷魚 2026-03-27 | 8/400 |
|
|
|
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考博] 26申博 +3 | 加油沖�。� 2026-03-26 | 3/150 |
|
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 349求調(diào)劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|