| 查看: 1107 | 回复: 1 | ||
[求助]
udf曳力相关问题求教。
|
|
写了一个曳力的UDF，但是开始计算就发现问题了。 同样的模型，加入了UDF（wen yu曳力系数）之后计算出错。利用自带的计算就没有问题 # Divergence detected in AMG solver: mp-x-momentum -> Decreasing coarsening group size! # Divergence detected in AMG solver: mp-x-momentum -> Increasing relaxation sweeps! # Divergence detected in AMG solver: pressure correction -> Turning off correction scaling! # Divergence detected in AMG solver: pressure correction -> Increasing relaxation sweeps! udf 如下： #include "udf.h" #define diam2 9.e-2 DEFINE_EXCHANGE_PROPERTY(LS_drag, cell, mix_thread, s_col, l_col) { Thread *thread_l, *thread_s; thread_l = THREAD_SUB_THREAD(mix_thread, l_col); thread_s = THREAD_SUB_THREAD(mix_thread, s_col); real abs_v, slip_x, slip_y, slip_z; real rho_l, mu_l, reyp, void_s, void_l, f_drag, k_l_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); rho_l = C_R(cell, thread_l); mu_l = C_MU_L(cell, thread_l); void_s = C_VOF(cell, thread_s); void_l = C_VOF(cell, thread_l); abs_v = sqrt(slip_x*slip_x + slip_y*slip_y + slip_z*slip_z); reyp = rho_l*abs_v*diam2/mu_l; f_drag = (reyp <= 1000) ? (24.0*(1.0 + 0.15*pow(reyp, 0.687))/reyp) :0.44; k_l_s = 0.75*rho_l*void_s*f_drag*abs_v/diam2; return k_l_s; } 其中换算了单位，所以颗粒的直径单位是mm。 |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有78人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有28人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有48人回复
2楼2017-09-26 19:16:32
