24小时热门版块排行榜     先丰纳米 买纳米材料就点这里

查看: 100  |  回复: 2
【悬赏金币】回答本帖问题,作者291571607将赠送您 6 个金币

291571607

新虫 (初入文坛)

[求助] 修改边界剪切力为速度的二阶导数出现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))算的话,又会使剪切力直接为零,球球各位大佬帮帮忙,感激不尽。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

291571607

新虫 (初入文坛)

ok,解决了,谢谢大家
2楼2018-05-28 09:02:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星星在我家门

木虫 (正式写手)

热爱互联网

引用回帖:
2楼: Originally posted by 291571607 at 2018-05-28 09:02:45
ok,解决了,谢谢大家

求问大佬,怎么解决的,我改了一个星期了,改不对
于无声听惊雷
3楼2018-06-14 18:46:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 291571607 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 还有半年毕业,各种愁 +7 Me诺言liu 2018-06-23 17/850 2018-06-23 14:07 by 啄木鸟先森
[基金申请] 时间点 +13 xiaomuch163 2018-06-19 17/850 2018-06-23 13:46 by 金蛋2012
[硕博家园] 读博 +10 zhangbarton 2018-06-23 10/500 2018-06-23 11:12 by jean枫叶
[公派出国] 海外个人申请高水平 +4 lucky0630 2018-06-22 4/200 2018-06-23 10:50 by zhaoyang14
[硕博家园] 有准备19年出国读博的虫虫吗?坐标西安 +14 张小可^ 2018-06-22 30/1500 2018-06-23 09:21 by 少根筋的大姐姐
[找工作] 中船重工702所到底怎么样 +6 aluba123 2018-06-22 8/400 2018-06-23 09:07 by aluba123
[找工作] 大连,杭州怎么选? +29 yourkathy 2018-06-17 31/1550 2018-06-23 08:04 by surry
[基金申请] 国家社科基金什么时候出项目号呀? +9 xiaokoala 2018-06-20 15/750 2018-06-22 18:12 by luckyken66
[博后之家] 港大博后求合租 +5 sdkj199017 2018-06-18 11/550 2018-06-22 16:01 by sdkj199017
[博后之家] 北京理工大学“复杂系统智能控制与决策”国家重点实验室招聘博士后若干 +5 yuch0006 2018-06-21 7/350 2018-06-22 11:53 by yuch0006
[药学] 做天药真的就没有发展前途吗? +10 阮思行 2018-06-21 18/900 2018-06-22 08:26 by along9227
[硕博家园] 国外读博士可以带老公和宝宝去吗? +9 风雪桃花 2018-06-21 11/550 2018-06-21 22:46 by fireme33
[硕博家园] 春节一起去旅游啊 +21 yangbing78 2018-06-20 33/1650 2018-06-21 19:47 by yangbing78
[硕博家园] 有人做环糊精吗,求大神带飞 +3 渣渣鱼 2018-06-20 12/600 2018-06-20 21:53 by 乐观的小虫
[找工作] 盐城师范学院新进博士待遇 +4 wenxiaoju 2018-06-19 8/400 2018-06-20 20:28 by xmy1983
[高分子] 请问这个图是用什么软件画出来的? +3 允允允儿~~~ 2018-06-19 3/150 2018-06-20 13:36 by 月只蓝
[论文投稿] 论文返修回去了,求接受 +16 pku周轩 2018-06-19 19/950 2018-06-20 10:28 by ycylovepsf
[论文投稿] 大修有多大可能性接收? +14 yeshuaichen 2018-06-17 15/750 2018-06-19 16:50 by 我乘着风
[公派出国] 求问有没有人知道没有拿到硕士学位可以通过CSC在欧洲申请读博吗? +6 xhtong 2018-06-17 15/750 2018-06-19 09:39 by xhtong
[留学DIY] 澳洲 格里菲斯大学 洁净环境与能源中心 招收博士 +4 robert_lau 2018-06-18 8/400 2018-06-19 09:20 by robert_lau
信息提示
请填处理意见