| 查看: 1172 | 回复: 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。 |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有264人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
2楼2017-09-26 19:16:32
20