 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1517 | 回復(fù): 4 | ||
guoli8832700新蟲 (初入文壇)
|
[求助]
求各位fluent大神，UDF曳力問題，幫幫忙，感激不盡。。。 已有1人參與
|
|
書上看到的一個(gè)流化床模擬，代碼我能看懂，但是為什么要用兩個(gè)氣固交換系數(shù)，書上是氣相選擇的custom_drag_syam，固相選擇custom_drag_ihme。還有一個(gè)問題，第一個(gè)用的是syamlal-o’brien模型，但是第二個(gè)我查了下，不是常用的幾個(gè)模型啊，沒搞懂，書上也沒說。大神幫忙看下，非常感謝。。。。 其中UDF代碼如下： /* This is a routine for customizing default Syamlal drag law in Fluent 6. The default drag law uses 0.8 (for void<=0.85) and 2.65 (void>0.85) for bfac. This is for a min fluid vel of 25 cm/s. The current drag law has been tuned for a min fluid vel of 8 cm/s and uses 0.28 and 9.07 for these parameters. */ #include "udf.h" #include "sg_mphase.h" # define pi 4.*atan(1.) #define diam2 3.e-4 DEFINE_EXCHANGE_PROPERTY(custom_drag_syam, cell, mix_thread, s_col, f_col) { Thread *thread_g, *thread_s; real x_vel_g, x_vel_s, y_vel_g, y_vel_s, abs_v, slip_x, slip_y, rho_g, rho_s, mu_g, reyp, afac, bfac, void_g, vfac, fdrgs, taup, k_g_s; /* find the threads for the gas (primary) and solids (secondary phases). These phases appear in columns 2 and 1 in the Interphase panel respectively*/ thread_g = THREAD_SUB_THREAD(mix_thread, s_col);/*gas phase*/ thread_s = THREAD_SUB_THREAD(mix_thread, f_col);/* solid phase*/ /* find phase velocities and properties*/ x_vel_g = C_U(cell, thread_g); y_vel_g = C_V(cell, thread_g); x_vel_s = C_U(cell, thread_s); y_vel_s = C_V(cell, thread_s); slip_x = x_vel_g - x_vel_s; slip_y = y_vel_g - y_vel_s; rho_g = C_R(cell, thread_g); rho_s = C_R(cell, thread_s); mu_g = C_MU_L(cell, thread_g); /*compute slip*/ abs_v = sqrt(slip_x*slip_x + slip_y*slip_y); /*compute reynolds number*/ reyp = rho_g*abs_v*diam2/mu_g; /* compute particle relaxation time */ taup = rho_s*diam2*diam2/18./mu_g; void_g = C_VOF(cell, thread_g);/* gas vol frac*/ /*compute drag and return drag coeff, k_g_s*/ afac = pow(void_g,4.14); if(void_g<=0.85) bfac = 0.281632*pow(void_g, 1.28); else bfac = pow(void_g, 9.076960); vfac = 0.5*(afac-0.06*reyp+sqrt(0.0036*reyp*reyp+0.12*reyp*(2.*bfac- afac)+afac*afac)); fdrgs = void_g*(pow((0.63*sqrt(reyp)/vfac+4.8*sqrt(vfac)/vfac),2))/24.0; k_g_s = (1.-void_g)*rho_s*fdrgs/taup; return k_g_s; } DEFINE_EXCHANGE_PROPERTY(custom_drag_ihme, cell, mix_thread, s_col, f_col) { Thread *thread_g, *thread_s; real x_vel_g, x_vel_s, y_vel_g, y_vel_s, abs_v, slip_x, slip_y, rho_g, rho_s, mu_g, reyp, cd, eg, void_g, k_g_s; /* find the threads for the gas (primary) and solids (secondary phases). These phases appear in columns 2 and 1 in the Interphase panel respectively*/ thread_g = THREAD_SUB_THREAD(mix_thread, s_col);/*gas phase*/ thread_s = THREAD_SUB_THREAD(mix_thread, f_col);/* solid phase*/ /* find phase velocities and properties*/ x_vel_g = C_U(cell, thread_g); y_vel_g = C_V(cell, thread_g); x_vel_s = C_U(cell, thread_s); y_vel_s = C_V(cell, thread_s); slip_x = x_vel_g - x_vel_s; slip_y = y_vel_g - y_vel_s; rho_g = C_R(cell, thread_g); rho_s = C_R(cell, thread_s); mu_g = C_MU_L(cell, thread_g); /*compute slip*/ abs_v = sqrt(slip_x*slip_x + slip_y*slip_y); /*compute reynolds number*/ reyp = rho_g*abs_v*diam2/mu_g; cd = (24./(reyp+SMALL)) + 5.48*pow((reyp+SMALL),-0.573) + 0.36; void_g = C_VOF(cell, thread_g);/* gas vol frac*/ eg = pow(void_g,-2.65); k_g_s = (3./4.)*(cd*(1.-void_g)*abs_v*rho_g*eg)/diam2; return k_g_s; } |
回復(fù)此樓» 猜你喜歡
2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有142人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有30人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有64人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- fluent UDF代碼調(diào)用預(yù)定義程序的問題（如fortran程序）
已經(jīng)有15人回復(fù)
- 64位ANSYS FLUENT不能compile UDF的問題
已經(jīng)有20人回復(fù)
- 【求助】FLUENT UDF 氣液兩相流 汽蝕 汽化潛熱(問題詳見內(nèi)容)
已經(jīng)有21人回復(fù)
- 好心人幫幫忙吧，為什么這個(gè)UDF導(dǎo)入fluent以后，在邊界條件處沒有顯示~~~
已經(jīng)有17人回復(fù)
- UDF載入fluent之后出錯(cuò)，求指點(diǎn)！
已經(jīng)有7人回復(fù)
- 曳力函數(shù)的的UDF
已經(jīng)有15人回復(fù)
- 【討論】Fluent邊界條件及虛擬質(zhì)量力
已經(jīng)有6人回復(fù)
- 【求助】linux版fluent使用udf并行問題請(qǐng)教
已經(jīng)有4人回復(fù)
- 【求助】如何改變Fluent中的曳力模型（UDF）
已經(jīng)有12人回復(fù)
liubao2012
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2126.3
- 散金: 10
- 紅花: 1
- 帖子: 208
- 在線: 123.9小時(shí)
- 蟲號(hào): 1886137
- 注冊(cè): 2012-07-10
- 性別: GG
- 專業(yè): 化學(xué)反應(yīng)工程
2樓2016-01-01 01:45:11
ZxrzxR0508
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 0.5
- 帖子: 10
- 在線: 13小時(shí)
- 蟲號(hào): 5568095
- 注冊(cè): 2017-02-13
- 性別: GG
- 專業(yè): 化學(xué)工程及工業(yè)化學(xué)
3樓2018-03-06 15:33:12
2017kaixin
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 68.6
- 帖子: 11
- 在線: 53.2小時(shí)
- 蟲號(hào): 6312437
- 注冊(cè): 2017-04-13
- 專業(yè): 燃燒學(xué)
4樓2018-04-17 10:52:23
jiajinghua
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 664
- 帖子: 137
- 在線: 7小時(shí)
- 蟲號(hào): 698525
- 注冊(cè): 2009-02-08
- 專業(yè): 化學(xué)反應(yīng)工程
5樓2018-04-18 09:07:08
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 guoli8832700 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +3 | 西南交通專材355 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +8 | 手機(jī)用戶 2026-03-17 | 8/400 |
|
|
|
[考研]
 271材料工程求調(diào)劑
+7
|
.6lL 2026-03-18 | 7/350 |
|
|
|
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +3 | Charlieyq 2026-03-19 | 3/150 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
|
[考研] 321求調(diào)劑 +8 | 何潤采123 2026-03-18 | 10/500 |
|
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研]
11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|