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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 还有谁现在还没收到评审意见的吗? +3 cathy1025 2018-08-18 3/150 2018-08-18 09:26 by glock117
[教师之家] 高校入职体检血压高会不会有影响? +14 wjswsl 2018-08-14 16/800 2018-08-17 10:57 by 司马缸APP
[基金申请] 工才个人可查E01 +5 z剑指苍天 2018-08-17 5/250 2018-08-17 06:48 by zzzgyx
[论文投稿] 计算机应用被增刊的来聊一下 +4 1846106498 2018-08-15 5/250 2018-08-16 17:56 by 1846106498
[基金申请] 基金凉凉 +13 banzhu最牛 2018-08-15 13/650 2018-08-16 16:44 by jayfelling
[基金申请] 地学口,个人账号有变化吗 +6 雷峰塔镇河妖 2018-08-16 8/400 2018-08-16 16:17 by leon29
[金属] 合金XRD分析 +4 sheqianwei 2018-08-16 10/500 2018-08-16 13:29 by sheqianwei
[基金申请] 去年合作项目的漏洞还能用吗?请中了的验证一下 +5 gaoyu_8848 2018-08-16 6/300 2018-08-16 12:58 by gaoyu_8848
[基金申请] 诸位,我昨天发的帖子说的没错吧。有些人就是喜欢作妖,恐天下不乱。 +8 chem.metal 2018-08-16 9/450 2018-08-16 12:14 by leiqhd
[基金申请] 基金问题询问 +5 益达核桃露 2018-08-16 7/350 2018-08-16 12:07 by wxyx2004
[硕博家园] 生活总在一天天变好 +8 怡慕 2018-08-11 8/400 2018-08-16 11:55 by simple__plan
[文学芳草园] 最后的日子 +5 myrtle 2018-08-15 6/300 2018-08-16 10:22 by 转身又一年
[硕博家园] 对巴拿赫空间很熟悉的?? +4 CapHope 2018-08-15 4/200 2018-08-16 00:07 by 王鹏飞3826
[基金申请] 未收到消息 +8 myqzw123 2018-08-15 9/450 2018-08-15 23:19 by 琵琶123
[基金申请] 定了!结果公布 +30 深缘浅笑 2018-08-14 31/1550 2018-08-15 08:07 by rumsfld
[硕博家园] 有谁了解索氏提取法吗? +5 海盗11111 2018-08-14 6/300 2018-08-15 07:07 by trywpl
[论文投稿] 记录突破 +20 cwwuchao 2018-08-14 25/1250 2018-08-14 15:09 by JasmineDeng-
[硕博家园] 共勉 +4 loveDYT 2018-08-13 4/200 2018-08-13 18:01 by 坑爹理查德
[基金申请] 有没有人上会过了 然后最后公布却发现没过的? +31 candy1117 2018-08-11 68/3400 2018-08-13 09:19 by mqstca
[论文投稿] 有投过 ACS表面的么? +8 南有先生 2018-08-12 8/400 2018-08-12 18:56 by yunxiang123
信息提示
请填处理意见