| 查看: 1412 | 回复: 3 | ||
[求助]
修改边界剪切力为速度的二阶导数出现segmemtation fault
|
|
我想将我的边界层剪切力修改为t=d2u/dx2,u为边界层上一点的速度，x为切向方向。就是将剪切力设定为速度在切向方向的二阶导但是在直接使用C_UDSI_G这个宏的时候会出现segmentation fault。下面是我的代码： # include "udf.h" # define domain_ID 2 DEFINE_ADJUST(adjust_gradient, domain) { Thread *t; cell_t c; face_t f; domain = Get_Domain(domain_ID); /* Fill UDS with the variable. */ thread_loop_c(t,domain) begin_c_loop (c,t) { C_UDSI(c,t,0) = C_DUDX(c,t); } end_c_loop (c,t) } DEFINE_ON_DEMAND(store_gradient) { Domain *domain; face_t f; cell_t c; Thread *t; domain=Get_Domain(1); /* Fill the UDM with magnitude of gradient. */ thread_loop_c(t,domain) begin_c_loop (c,t) { C_UDMI(c,t,0) = C_UDSI_G(c,t,0)[0]; } end_c_loop (c,t) } DEFINE_PROFILE(shear_stressx, t,position) { face_t f; cell_t c0; real xw[ND_ND],xc[ND_ND],taux; real a = 0.00012; Thread *t0= t->t0; begin_f_loop(f,t) { c0= F_C0(f,t); F_CENTROID(xw,f,t); C_CENTROID(xc,c0,t0); taux=-a*C_UDMI(c0,t0,0); F_PROFILE(f,t,position) = taux; } end_f_loop(f,t) } 在execute on demand 的时候就会出现SEGMENTATION 错误，但是如果将C_UDMI(c,t,0) = C_UDSI_G(c,t,0)[0];改成C_UDMI(c,t,0) =NV_MAG(C_UDSI_G(c,t,0))算的话，又会使剪切力直接为零，球球各位大佬帮帮忙，感激不尽。 |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有137人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
2楼2018-05-28 09:02:45
星星在我家门
木虫 (著名写手)
热爱互联网
- 应助: 21 (小学生)
- 金币: 2024.1
- 散金: 1232
- 红花: 8
- 沙发: 1
- 帖子: 1070
- 在线: 173.8小时
- 虫号: 3229413
- 注册: 2014-05-23
- 专业: 传递过程
3楼2018-06-14 18:46:42
4楼2018-09-29 20:49:50
