24小时热门版块排行榜     意得辑英文润色,师兄用过都说好

查看: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 沈阳化X研究院,有黑恶势力! +3 不太冷190 2018-10-21 3/150 2018-10-21 20:31 by 海上的小船
[考研] 中科院化学所 +3 隔壁小明 2018-10-21 4/200 2018-10-21 18:39 by 奇妙能力L
[考博] 考博名额 20+3 刘 帅 2018-10-16 12/600 2018-10-21 16:21 by ?快乐人生
[博后之家] 博士后跟别的临时工作有什么区别? +6 如猫似狗 2018-10-20 7/350 2018-10-21 15:29 by sun881014
[硕博家园] 搞科研怎么这么难呢 +5 娟娟 2018-10-18 6/300 2018-10-21 15:27 by 柳絮虫
[硕博家园] 博士的感情 +14 我就是我1 2018-10-15 38/1900 2018-10-21 14:42 by Pardon168
[论文投稿] 记录投稿时间线~ +3 Bubu_Xie 2018-10-21 4/200 2018-10-21 12:53 by disoly
[硕博家园] 做实验的意义在哪里 +8 Alfalfl 2018-10-19 8/400 2018-10-21 11:31 by dlzhao516
[论文投稿] 如何成为一位高产的SCI学者 (金币+5) +35 NES润色 2018-10-19 38/1900 2018-10-20 18:12 by 紫夜星痕
[仿真模拟] 请问有人会aspen plus么?求助大神。 +4 志在爱颖 2018-10-17 6/300 2018-10-20 15:19 by 还是回家好啊
[硕博家园] 职业瓶颈期,要不要继续回去读书学习 +14 qisi8395 2018-10-18 15/750 2018-10-20 14:20 by yongerli
[精细化工] 布艺沙发清洗剂 +4 合肥威尔 2018-10-18 7/350 2018-10-19 12:51 by sljchem
[硕博家园] 导师不管,自己本身基础就差,论文写不出来,伤心啊 +12 向日葵CO 2018-10-18 19/950 2018-10-18 20:45 by 飞行鸟
[考博] 老师什么么意思啊 +14 爱7号的me 2018-10-16 22/1100 2018-10-18 19:05 by 比得兔嘿嘿
[考博] 审核制读博 +5 tell_more 2018-10-18 5/250 2018-10-18 18:49 by xzs528
[硕博家园] 版主 +4 yangzi1023 2018-10-16 8/400 2018-10-17 09:25 by yangzi1023
[论文投稿] OPTIK大修提交了,祈祷 +3 飞飞的晨曦 2018-10-15 3/150 2018-10-17 03:06 by 皇家小虫
[论文投稿] 投JPS,审稿五个多月拒搞 +4 只能说我忘了 2018-10-16 7/350 2018-10-16 17:20 by 只能说我忘了
[海外博后] 8个月的欧洲博士后值得去吗 +14 星垂野沃 2018-10-15 15/750 2018-10-16 12:11 by imchentang
[硕博家园] 23岁的青春别的妹子纠结用那支口红我纠结做实验还是写paper +7 liwei_chen 2018-10-15 10/500 2018-10-15 19:45 by sosososo1114
信息提示
请填处理意见