 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1063 | 回復: 3 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
求組大俠給我看看鄙人編寫的UDF，感激不盡……小弟沒有金幣啊嗚嗚
|
||
|
#include "udf.h" #include "sg_mphase.h" #define A 16.0114 #define B 3456.80 #define C -8.67 #define L 0.00000125.0 #define POP 101325.0 #define K 0.1 #define C10H22_NUM 0 #define LIQ_PHASE_ID 0 #define GAS_PHASE_ID 1 DEFINE_SOURCE(liq_zhengkuiwan,c,liq_thread,dS,eqn) { int i; cell_t c; Thread*gas_thread,*mixture_thread; Domain*mixture_domain,*liq_domain,*gas_domain; Material *gas_material,*species; real c10h22_mole_fract,P_c10h22,total_mole,Mw[MAX_SPE_EQNS]; mixture_domain=Get_Domain(1); liq_domain=DOMAIN_SUB_DOMAIN(mixture_domain,LIQ_PHASE_ID); gas_domain=DOMAIN_SUB_DOMAIN(mixture_domain,GAS_PHASE_ID); //感覺沒什么用，第一次編寫udf就都寫上了 mixture_thread=THREAD_SUPER_THREAD(liq_thread); gas_thread=THREAD_SUB_THREAD(mixture_thread,GAS_PHASE_ID); gas_material=THREAD_MATERIAL(gas_thread); real m_dot_liq,c_s,c_g,m_kuosan; real x[ND_ND],y=x[1]; real ps,v; //ps為液滴表面飽和壓力 v為單個網格的體積 thread_loop_c(mixture_thread,mixture_domain) //在整個二維控制區(qū)查詢混合物的單元線 { begin_c_loop(c,mixture_thread) //在單元線上查詢單元 {if(0.01<=C_VOF(c,liq_thread)&&C_VOF(c,liq_thread)<1.0) //查詢氣液界面處的網格，我選了0.01到1之間，因為當液滴的體積分數小于0.01時會產生沸騰蒸發(fā) { ps=133*exp(A-B/(C_T(c,liq_th)+C); //Antoine蒸氣壓方程求得液滴表面飽和壓力ps c_s=ps/(C_RGAS(c,gas_thread)*C_T(c,liq_thread)); //求得液滴表面濃度 C_RGAS(c,gas_thread)氣體常數這樣可以嗎？？也可以直接輸入8.3134 total_mole=0; mixture_species_loop(gas_material,species,i) { Mw=MATERIAL_PROP(species,PROP_mwi); total_mole+=C_YI(c,gas_thread,i)/Mw; //求得氣液界面網格內氣體的總摩爾數 } c10h22_mole_fract=(C_YI(c,gas_thread,C10H22_NUM)/Mw[C10H22_NUM])/total_mole; //求得氣液界面內c10h22的摩爾分數 P_c10h22=POP*c10h22_mole_fract; //求得氣液界面內c10h22的分壓 c_g=P_c10h22/(C_RGAS(c,gas_thread)*C_T(c,liq_thread)); //求得氣液界面網格內氣體中c10h22的濃度 m_kuosan=c_s-c_g; //氣液界面網格內液滴表面c10h22濃度和氣體中c10h22濃度的差值，從而產生擴散 v=L*L*L; m_dot_liq=142*K*y*m_kuosan*C_DIFF_EFF(c,gas_thread,i)/v*C_VOF(c,liq_th)*pow(1-C_VOF(c,liq_thread),2.0)*(1-(1-C_VOF(c,liq_thread))*c10h22_mole_fract); dS[eqn]= //最終源項 dS還沒有求，比較麻煩 } } end_c_loop(c,mixture_thread) return m_dot_liq; } [ Last edited by jimoxingzhe on 2013-4-9 at 16:14 ] |
回復此樓» 猜你喜歡
289求調劑
已經有5人回復
-
一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩
已經有5人回復
-
一志愿西安交通大學 學碩 354求調劑
已經有3人回復
-
317求調劑
已經有11人回復
-
287求調劑
已經有6人回復
-
一志愿南昌大學，327分，材料與化工085600
已經有7人回復
-
材料080500調劑求收留
已經有8人回復
-
353求調劑
已經有3人回復
-
北科281學碩材料求調劑
已經有3人回復
-
A區(qū)線材料學調劑
已經有3人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 鄙人手上有三個化合物，求過路大俠鑒定有沒有做過晶體的研究。。。感激不盡
已經有13人回復
3樓2013-04-09 17:20:57
2樓2013-04-09 17:20:17
匿名
用戶注銷 (著名寫手)
- 應助: 84 (初中生)
- 金幣: 5015.3
- 散金: 1030
- 紅花: 35
- 帖子: 2148
- 在線: 235.4小時
- 蟲號: 0
- 注冊: 2011-10-29
- 性別: GG
- 專業(yè): 伊斯蘭教
4樓2013-04-13 22:51:05
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學 學碩 354求調劑 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 287求調劑 +5 | 晨昏線與星海 2026-03-19 | 6/300 |
|
|
[考研] 一志愿南昌大學，327分，材料與化工085600 +7 | Ncdx123456 2026-03-19 | 7/350 |
|
|
|
[考研] 材料學求調劑 +3 | Stella_Yao 2026-03-20 | 3/150 |
|
|
|
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調劑） +3 | 123阿標 2026-03-17 | 3/150 |
|
|
|
[考研] 求調劑 +4 | Mqqqqqq 2026-03-19 | 4/200 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
|
[考研] 304求調劑 +5 | 曼殊2266 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調劑 +3 | 葵梓衛(wèi)隊 2026-03-18 | 5/250 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
|
[考研] 0703化學調劑，求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
|
[考研] 一志愿華中科技大學，080502，354分求調劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
|
[考研] 08工科 320總分 求調劑 +5 | 梨花珞晚風 2026-03-17 | 5/250 |
|
|
|
[考研] 311求調劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研]
 0854，計算機類招收調劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 301求調劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 機械專碩325，尋找調劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 326求調劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|