24小时热门版块排行榜     意得辑论文翻译润色,首单减免200元

查看: 150  |  回复: 3
【悬赏金币】回答本帖问题,作者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的回帖

佩佩LIU

铜虫 (初入文坛)

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

????神仙啊  怎么就自己解决了  能教教我吗
4楼2018-09-29 20:49:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 291571607 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 读博士困惑 +16 我就是我1 2019-01-16 17/850 2019-01-21 13:33 by 飞飞=飞飞
[硕博家园] 坚持 +15 培培奋斗ing 2019-01-21 16/800 2019-01-21 11:42 by 199303
[论文投稿] 研究生二年级了,我的第一篇论文,老师看过之后建议我投nature nanotechnology。 +58 小猫真可爱 2019-01-18 71/3550 2019-01-21 10:56 by mmYYxx
[硕博家园] 爱江山更爱美人 +13 余生请多 2019-01-19 13/650 2019-01-21 10:49 by 十九小妹
[硕博家园] 潍坊妹纸,诚征男友 +19 冷月茕钟 2019-01-20 36/1800 2019-01-21 08:54 by 赵帅。
[基金申请] 国基金面上项目要求的5篇代表作是近几年的? +4 jbz1124 2019-01-20 5/250 2019-01-21 06:38 by Aaronyin
[考博] 大方广圓覺了義經 +4 青木桴蓠 2019-01-20 8/400 2019-01-20 22:31 by 青木桴蓠
[硕博家园] 关于国内Top 和国外二流,多日的纠结束。。。 +37 江湖学术 2019-01-17 67/3350 2019-01-20 17:02 by 格林达姆少校
[论文投稿] sci自重 10+3 srxmimi 2019-01-20 8/400 2019-01-20 16:22 by nono2009
[文学芳草园] 我们的纪念 +8 131820 2019-01-16 14/700 2019-01-20 14:11 by Pardon168
[硕博家园] 关于硕士选导师 +11 清灵7066 2019-01-16 15/750 2019-01-20 00:31 by kl5top
[海外博后] 收到捷克的offer 求教捷克生活咋样 +5 bhl321023 2019-01-19 10/500 2019-01-19 23:05 by 海之蓝
[有机交流] 求指导一个反应机理 20+6 南昌吴彦祖 2019-01-17 24/1200 2019-01-19 08:07 by 南昌吴彦祖
[育儿交流] 关于备孕的经验分享。。。后续会分享有关孕期,以及生产经验,以及带宝宝经验 +7 兜兜里有糖12 2019-01-15 26/1300 2019-01-18 08:19 by 兜兜里有糖12
[考博] 材料丢了 +5 猫喵 2019-01-16 9/450 2019-01-17 07:29 by 猫喵
[硕博家园] 怎么看xps数据 +3 cj516294429 2019-01-15 5/250 2019-01-16 20:19 by cj516294429
[基金申请] 国家自然科学基金资助项目变更管理规程(试行) +4 天天向上3035 2019-01-16 6/300 2019-01-16 12:17 by fallenrain
[硕博家园] 对于我上一天帖子补充 +3 余生请多 2019-01-16 3/150 2019-01-16 11:02 by 宝藏箱
[硕博家园] 各位好 +3 Aliceyonggan 2019-01-16 3/150 2019-01-16 10:12 by hdddage
[有机交流] 氢氧化钠加少量水能产生燃烧 +4 Huisekuai 2019-01-14 6/300 2019-01-15 16:07 by ghj-hebust
信息提示
请填处理意见