24小时热门版块排行榜    

查看: 1181  |  回复: 4
【悬赏金币】回答本帖问题,作者打工预备役将赠送您 5 个金币

打工预备役

新虫 (小有名气)

[求助] udf颗粒沉积中求第一层网格高度函数 已有1人参与

想用ds来求距离壁面第一层网格的高度
BOUNDARY_FACE_GEOMETRY(f,t,A,ds,es,A_by_es,dr0);
tauwall1=C_MU_EFF(c,tcell)(du/ds);
wallfricv1=sqrt(tauwall1/C_R(c,tcell));
if (yplus<11.25)
{
tauwall2=C_MU_EFF(c,tcell)C_STRAIN_RATE_MAG(c,tcell);
wallfricv2=sqrt(tauwall2/C_R(c,tcell));
}
else
{
wallfricv2 = (1./0.41)log(yplus9.);
}
ds=C_WALL_DIST(c,tcell);
wallfricv4=(VISCyplus)/(dsC_R(c,tcell));

```
但这里面的ds好像是以前版本的不太正确,现在想用网格中心到壁面距离乘2来表示,但我在帮助文档里没找到求距离的函数,请问有人知道吗?求指导
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何弃疗啊哥

铜虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
1)定义一个数组x[ND_ND],然后x = x[0],y = y[1]。分别表示网格中心点在x和y方向的位置

2)或者就是C_WALL_DIST(c,t),求的是网格中心点到壁面的距离,我用的是第二种。需要注意的是,如果是上下两侧都有壁面,C_WALL_DIST求得的是网格中心与距离最近的那个壁面的距离

» 本帖已获得的红花(最新10朵)

2楼2023-01-03 10:36:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

打工预备役

新虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by 何弃疗啊哥 at 2023-01-03 10:36:06
1)定义一个数组x,然后x = x,y = y。分别表示网格中心点在x和y方向的位置

2)或者就是C_WALL_DIST(c,t),求的是网格中心点到壁面的距离,我用的是第二种。需要注意的是,如果是上下两侧都有壁面,C_WALL_DI ...

那请问这一句ds=C_WALL_DIST(c,tcell);是不是就可以直接用啊?假设距离壁面就是最近,乘以2就是网格高度

发自小木虫Android客户端
3楼2023-01-03 12:07:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何弃疗啊哥

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by 打工预备役 at 2023-01-03 12:07:55
那请问这一句ds=C_WALL_DIST(c,tcell);是不是就可以直接用啊?假设距离壁面就是最近,乘以2就是网格高度
...

C_WALL_DIST本身是可以直接用的,c和t跟你的udf指针相对应就行。其他的应该也没啥问题
4楼2023-01-03 14:06:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

打工预备役

新虫 (小有名气)

引用回帖:
4楼: Originally posted by 何弃疗啊哥 at 2023-01-03 14:06:44
C_WALL_DIST本身是可以直接用的,c和t跟你的udf指针相对应就行。其他的应该也没啥问题...

好嘞,谢谢您解答了,感谢感谢

发自小木虫Android客户端
5楼2023-01-03 14:08:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 打工预备役 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 max7zaruFVS8 2026-10-06 3/150 2026-10-07 16:51 by kUQbiRsIeOpF
[考博] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 z1aZaKkVi71s 2026-10-06 3/150 2026-10-07 16:31 by kUQbiRsIeOpF
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 z1aZaKkVi71s 2026-10-06 3/150 2026-10-07 15:41 by kUQbiRsIeOpF
[博后之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 z1aZaKkVi71s 2026-10-06 3/150 2026-10-07 15:41 by kUQbiRsIeOpF
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 z1aZaKkVi71s 2026-10-05 3/150 2026-10-07 15:09 by kUQbiRsIeOpF
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 FlhJF7cq3uT3 2026-10-04 8/400 2026-10-07 11:11 by EXuX9l7TjjIY
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 FlhJF7cq3uT3 2026-10-04 9/450 2026-10-07 11:07 by EXuX9l7TjjIY
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 FlhJF7cq3uT3 2026-10-04 8/400 2026-10-07 11:01 by EXuX9l7TjjIY
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 FlhJF7cq3uT3 2026-10-04 7/350 2026-10-07 10:51 by EXuX9l7TjjIY
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 FlhJF7cq3uT3 2026-10-04 11/550 2026-10-07 10:29 by EXuX9l7TjjIY
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 FlhJF7cq3uT3 2026-10-04 8/400 2026-10-07 10:21 by EXuX9l7TjjIY
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 FlhJF7cq3uT3 2026-10-04 9/450 2026-10-07 10:11 by EXuX9l7TjjIY
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 FlhJF7cq3uT3 2026-10-04 13/650 2026-10-07 10:09 by EXuX9l7TjjIY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 FlhJF7cq3uT3 2026-10-04 7/350 2026-10-07 10:07 by EXuX9l7TjjIY
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 FUSxPqiIHRjH 2026-10-04 10/500 2026-10-07 10:01 by EXuX9l7TjjIY
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 FUSxPqiIHRjH 2026-10-04 11/550 2026-10-07 09:51 by EXuX9l7TjjIY
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FUSxPqiIHRjH 2026-10-04 12/600 2026-10-07 09:49 by EXuX9l7TjjIY
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 FUSxPqiIHRjH 2026-10-04 9/450 2026-10-06 14:28 by x7yLzftaCUEL
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 FlhJF7cq3uT3 2026-10-05 4/200 2026-10-06 05:08 by z1aZaKkVi71s
[教师之家] 烦躁 +5 水冰月月野兔 2026-10-03 5/250 2026-10-05 21:59 by 火星超人xi
信息提示
请填处理意见