| 1 | 1/1 | 返回列表 |
| 查看: 675 | 回復: 0 | ||
xxfreedomxx新蟲 (初入文壇)
|
[求助]
求助各位大神,我的wen-yu曳力UDF不能收斂
|
|
求助各位大神,如下是我的UDF,我發(fā)現(xiàn)自己寫的別的曳力公式和這個wen-yu的UDF是一個情況,我做的是固液攪拌體系,500個時間步之后第二項速度在10-3以上就收斂不了了,慢慢增大,調節(jié)收斂因子,時間步長,湍流模型,都不行,單純用自帶模型就沒有問題,選擇的時候加入了brucato的模型進行修正。在這里卡了兩個星期了。 #include "udf.h" #define pi 4.*atan(1.) DEFINE_EXCHANGE_PROPERTY(s_l_drag_IL_sl, cell, mix_thread, s_col, f_col) { Thread *thread_l, *thread_s; real abs_v, slip_x, slip_y, slip_z; real reyp, void_s, rho_l, rho_s, mu_l, fdrag, k_l_s, void_l, diam2, Res; /* find the threads for the liquid(primary) */ /* and solids (secondary phases) */ thread_l = THREAD_SUB_THREAD(mix_thread, f_col); /* liquid phase */ thread_s = THREAD_SUB_THREAD(mix_thread, s_col); /* solid phase */ /* find phase velocities and properties*/ rho_l = C_R(cell, thread_l); mu_l = C_MU_L(cell, thread_l); rho_s = C_R(cell, thread_s); diam2 = C_PHASE_DIAMETER(cell, thread_s); slip_x = C_U(cell, thread_l) - C_U(cell, thread_s); slip_y = C_V(cell, thread_l) - C_V(cell, thread_s); slip_z = C_W(cell, thread_l) - C_W(cell, thread_s); /*compute slip*/ abs_v = sqrt(slip_x*slip_x + slip_y*slip_y + slip_z*slip_z); /*compute Reynold's number*/ reyp = rho_l*abs_v*diam2/mu_l; /* compute particle relaxation time */ void_s = C_VOF(cell, thread_s);/* solid vol frac*/ void_l = C_VOF(cell, thread_l);/* liquid vol frac*/ Res=void_l*reyp; if(Res<=1000){ fdrag = 24.*(1. + 0.15*pow(Res, 0.687)/Res; } else(Res>1000){ fdrag = 0.44; } k_l_s = 0.75*void_s*void_l*mu_l*fdrag*abs_v*pow(void_l, -2.65)/diam2; return k_l_s; } ![]() |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學/290求調劑/本科經(jīng)歷豐富/工科也可 +6 | 丹青奶蓋 2026-03-26 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 尋找調劑 +5 | 倔強芒? 2026-03-21 | 8/400 |
|
|
[考研] 一志愿北京化工大學材料與化工 264分各科過A區(qū)國家線 +6 | 哈哈157349 2026-03-21 | 6/300 |
|
|
[考研] 309求調劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 求調劑 +3 | 蘆lty 2026-03-25 | 4/200 |
|
|
[考研] 334分 一志愿武理-080500 材料求調劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 求調劑 +3 | QiMing7 2026-03-25 | 3/150 |
|
|
[考研] 招08考數(shù)學 +8 | laoshidan 2026-03-20 | 17/850 |
|
|
[考研] 285求調劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 0854電子信息求調劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
[考研] 食品專碩 一志愿雙一流 328 +3 | xiaom99 2026-03-21 | 4/200 |
|
|
[考研] 291求調劑 +3 | HanBeiNingZC 2026-03-24 | 3/150 |
|
|
[考研] 材料專碩331求調劑 +4 | 鮮當牛 2026-03-24 | 4/200 |
|
|
[考研] 一志愿吉大化學322求調劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
[考研] 材料/農(nóng)業(yè)專業(yè),07/08開頭均可,過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調劑 +4 | 幸運的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 求調劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 260求調劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 0703化學297求調劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 材料 271求調劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|