 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1531 | 回復: 4 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
guoli8832700新蟲 (初入文壇)
|
[求助]
求各位fluent大神，UDF曳力問題，幫幫忙，感激不盡。。。 已有1人參與
|
|
|
書上看到的一個流化床模擬，代碼我能看懂，但是為什么要用兩個氣固交換系數(shù)，書上是氣相選擇的custom_drag_syam，固相選擇custom_drag_ihme。還有一個問題，第一個用的是syamlal-o’brien模型，但是第二個我查了下，不是常用的幾個模型啊，沒搞懂，書上也沒說。大神幫忙看下，非常感謝。。。。 其中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; } |
回復此樓» 猜你喜歡
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復
-
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！
已經(jīng)有1人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有144人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復
-
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復
-
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院
已經(jīng)有26人回復
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復
-
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院
已經(jīng)有8人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- fluent UDF代碼調(diào)用預定義程序的問題（如fortran程序）
已經(jīng)有15人回復
- 64位ANSYS FLUENT不能compile UDF的問題
已經(jīng)有20人回復
- 【求助】FLUENT UDF 氣液兩相流 汽蝕 汽化潛熱(問題詳見內(nèi)容)
已經(jīng)有21人回復
- 好心人幫幫忙吧，為什么這個UDF導入fluent以后，在邊界條件處沒有顯示~~~
已經(jīng)有17人回復
- UDF載入fluent之后出錯，求指點！
已經(jīng)有7人回復
- 曳力函數(shù)的的UDF
已經(jīng)有15人回復
- 【討論】Fluent邊界條件及虛擬質(zhì)量力
已經(jīng)有6人回復
- 【求助】linux版fluent使用udf并行問題請教
已經(jīng)有4人回復
- 【求助】如何改變Fluent中的曳力模型（UDF）
已經(jīng)有12人回復
ZxrzxR0508
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 0.5
- 帖子: 10
- 在線: 13小時
- 蟲號: 5568095
- 注冊: 2017-02-13
- 性別: GG
- 專業(yè): 化學工程及工業(yè)化學
3樓2018-03-06 15:33:12
liubao2012
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2126.3
- 散金: 10
- 紅花: 1
- 帖子: 208
- 在線: 123.9小時
- 蟲號: 1886137
- 注冊: 2012-07-10
- 性別: GG
- 專業(yè): 化學反應(yīng)工程
2樓2016-01-01 01:45:11
2017kaixin
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 68.6
- 帖子: 11
- 在線: 53.2小時
- 蟲號: 6312437
- 注冊: 2017-04-13
- 專業(yè): 燃燒學
4樓2018-04-17 10:52:23
jiajinghua
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 664
- 帖子: 137
- 在線: 7小時
- 蟲號: 698525
- 注冊: 2009-02-08
- 專業(yè): 化學反應(yīng)工程
5樓2018-04-18 09:07:08
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 生物學 377分 +4 | zzll03 2026-03-31 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
|
[考研] 調(diào)劑申請 +5 | 張張張張zy 2026-03-31 | 5/250 |
|
|
[有機交流]
 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機理
 10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
|
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調(diào)劑！ +7 | Zuhui0306 2026-03-25 | 7/350 |
|
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
|
[考研] 292求調(diào)劑 +13 | 是妍子也是研子 2026-03-30 | 13/650 |
|
|
|
[考研] 293求調(diào)劑 +3 | 末未mm 2026-03-30 | 5/250 |
|
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
|
[考研] 310求調(diào)劑 +10 | 爭取九點睡 2026-03-30 | 10/500 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 調(diào)劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
|
[考研] 復試調(diào)劑，一志愿南農(nóng)083200食品科學與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 中國科學院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|