| 查看: 773 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者291571607将赠送您 6 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
修改边界剪切力为速度的二阶导数出现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))算的话，又会使剪切力直接为零，球球各位大佬帮帮忙，感激不尽。 |
回复此楼» 猜你喜欢
SCAPS-1D学习交流
已经有2人回复
-
深圳大学跟哈工大联合博士
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有114人回复
-
Journal of alloys and compounds期刊
已经有2人回复
-
请教一个概率问题
已经有0人回复
-
铁氧体材料在太赫兹波段的响应时间
已经有0人回复
-
optisworks 的软件和教程
已经有0人回复
-
求助FAPbI3 cif文件
已经有0人回复
-
计算开关的响应时间
已经有0人回复
-
今年面上预计又打酱油了
已经有14人回复
4楼2018-09-29 20:49:50
20