 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1065 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
求組大俠給我看看鄙人編寫的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為單個(gè)網(wǎng)格的體積 thread_loop_c(mixture_thread,mixture_domain) //在整個(gè)二維控制區(qū)查詢混合物的單元線 { begin_c_loop(c,mixture_thread) //在單元線上查詢單元 {if(0.01<=C_VOF(c,liq_thread)&&C_VOF(c,liq_thread)<1.0) //查詢氣液界面處的網(wǎng)格，我選了0.01到1之間，因?yàn)楫?dāng)液滴的體積分?jǐn)?shù)小于0.01時(shí)會產(chǎn)生沸騰蒸發(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)氣體常數(shù)這樣可以嗎？？也可以直接輸入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; //求得氣液界面網(wǎng)格內(nèi)氣體的總摩爾數(shù) } c10h22_mole_fract=(C_YI(c,gas_thread,C10H22_NUM)/Mw[C10H22_NUM])/total_mole; //求得氣液界面內(nèi)c10h22的摩爾分?jǐn)?shù) P_c10h22=POP*c10h22_mole_fract; //求得氣液界面內(nèi)c10h22的分壓 c_g=P_c10h22/(C_RGAS(c,gas_thread)*C_T(c,liq_thread)); //求得氣液界面網(wǎng)格內(nèi)氣體中c10h22的濃度 m_kuosan=c_s-c_g; //氣液界面網(wǎng)格內(nèi)液滴表面c10h22濃度和氣體中c10h22濃度的差值，從而產(chǎn)生擴(kuò)散 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]= //最終源項(xiàng) dS還沒有求，比較麻煩 } } end_c_loop(c,mixture_thread) return m_dot_liq; } [ Last edited by jimoxingzhe on 2013-4-9 at 16:14 ] |
回復(fù)此樓» 猜你喜歡
308求調(diào)劑
已經(jīng)有3人回復(fù)
-
321求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有8人回復(fù)
-
330求調(diào)劑
已經(jīng)有4人回復(fù)
-
304求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑
已經(jīng)有6人回復(fù)
-
南京大學(xué)化學(xué)376求調(diào)劑
已經(jīng)有3人回復(fù)
-
考研調(diào)劑求學(xué)校推薦
已經(jīng)有5人回復(fù)
-
323求調(diào)劑
已經(jīng)有3人回復(fù)
-
317求調(diào)劑
已經(jīng)有12人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 鄙人手上有三個(gè)化合物，求過路大俠鑒定有沒有做過晶體的研究。。。感激不盡
已經(jīng)有13人回復(fù)
匿名
用戶注銷 (著名寫手)
- 應(yīng)助: 84 (初中生)
- 金幣: 5015.3
- 散金: 1030
- 紅花: 35
- 帖子: 2148
- 在線: 235.4小時(shí)
- 蟲號: 0
- 注冊: 2011-10-29
- 性別: GG
- 專業(yè): 伊斯蘭教
4樓2013-04-13 22:51:05
2樓2013-04-09 17:20:17
3樓2013-04-09 17:20:57
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|