 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1062 | 回復(fù): 3 | ||
jimoxingzhe新蟲(chóng) (初入文壇)
|
[求助]
求組大俠給我看看鄙人編寫(xiě)的UDF，感激不盡……小弟沒(méi)有金幣啊嗚嗚
|
|
#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); //感覺(jué)沒(méi)什么用，第一次編寫(xiě)udf就都寫(xiě)上了 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ū)查詢(xún)混合物的單元線(xiàn) { begin_c_loop(c,mixture_thread) //在單元線(xiàn)上查詢(xún)單元 {if(0.01<=C_VOF(c,liq_thread)&&C_VOF(c,liq_thread)<1.0) //查詢(xún)氣液界面處的網(wǎng)格，我選了0.01到1之間，因?yàn)楫?dāng)液滴的體積分?jǐn)?shù)小于0.01時(shí)會(huì)產(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還沒(méi)有求，比較麻煩 } } end_c_loop(c,mixture_thread) return m_dot_liq; } [ Last edited by jimoxingzhe on 2013-4-9 at 16:14 ] |
回復(fù)此樓» 猜你喜歡
材料與化工專(zhuān)碩調(diào)劑
已經(jīng)有6人回復(fù)
-
086500 325 求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料080500調(diào)劑求收留
已經(jīng)有7人回復(fù)
-
一志愿南昌大學(xué)，327分，材料與化工085600
已經(jīng)有6人回復(fù)
-
環(huán)境工程調(diào)劑
已經(jīng)有9人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有9人回復(fù)
-
287求調(diào)劑
已經(jīng)有5人回復(fù)
-
22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù)
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有3人回復(fù)
-
317求調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 鄙人手上有三個(gè)化合物，求過(guò)路大俠鑒定有沒(méi)有做過(guò)晶體的研究。。。感激不盡
已經(jīng)有13人回復(fù)
jimoxingzhe
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 17.8
- 帖子: 4
- 在線(xiàn): 6.7小時(shí)
- 蟲(chóng)號(hào): 2109486
- 注冊(cè): 2012-11-06
2樓2013-04-09 17:20:17
jimoxingzhe
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 17.8
- 帖子: 4
- 在線(xiàn): 6.7小時(shí)
- 蟲(chóng)號(hào): 2109486
- 注冊(cè): 2012-11-06
3樓2013-04-09 17:20:57
匿名
用戶(hù)注銷(xiāo) (著名寫(xiě)手)
- 應(yīng)助: 84 (初中生)
- 金幣: 5015.3
- 散金: 1030
- 紅花: 35
- 帖子: 2148
- 在線(xiàn): 235.4小時(shí)
- 蟲(chóng)號(hào): 0
- 注冊(cè): 2011-10-29
- 性別: GG
- 專(zhuān)業(yè): 伊斯蘭教
4樓2013-04-13 22:51:05
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 jimoxingzhe 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶(hù)573181 2026-03-20 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-19 | 9/450 |
|
|
|
[考研] 274求調(diào)劑 +8 | S.H1 2026-03-18 | 8/400 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|