24小时热门版块排行榜    

查看: 987  |  回复: 6

UDorcer

新虫 (小有名气)

[求助] 求助UDF编写问题!!! 已有1人参与

小弟在编写UDF时,遇到以下问题:
在DEFINE_SOURCE宏里面,需要用一个与坐标有关的函数,例如return source=-5*x;
即一串代码里面需要用到坐标x,
请问怎样把坐标写进去?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tank_fire

铜虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
摘自14.5的help说明
/*******************************************************************/
    UDF for specifying an x-momentum source term in a spatially
   dependent porous media
*******************************************************************/

#include "udf.h"

#define C2 100.0

DEFINE_SOURCE(xmom_source,c,t,dS,eqn)
{
    real x[ND_ND];
    real con, source;
    C_CENTROID(x,c,t);
    con = C2*0.5*C_R(c,t)*x[1];
    source = -con*fabs(C_U(c, t))*C_U(c,t);
    dS[eqn] = -2.*con*fabs(C_U(c,t));
    return source;
}


    C_CENTROID(x,c,t);这句是把当前cell的中心位置赋值给x[ND_ND]   ND_ND如果是2D就是2
3D就是3
所以你要的坐标就是xy(z)   x[0] x[1] (x[2])
2楼2014-08-27 23:08:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

UDorcer

新虫 (小有名气)

引用回帖:
2楼: Originally posted by tank_fire at 2014-08-27 23:08:16
摘自14.5的help说明
/*******************************************************************/
    UDF for specifying an x-momentum source term in a spatially
   dependent porous media
************** ...

但是为什么你给的例子中说是沿着X轴源项的变化,却用了x[1],x轴不是x[0]吗
3楼2014-08-28 09:44:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

UDorcer

新虫 (小有名气)

引用回帖:
2楼: Originally posted by tank_fire at 2014-08-27 23:08:16
摘自14.5的help说明
/*******************************************************************/
    UDF for specifying an x-momentum source term in a spatially
   dependent porous media
************** ...

非常感谢你
4楼2014-08-28 09:58:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

UDorcer

新虫 (小有名气)

引用回帖:
2楼: Originally posted by tank_fire at 2014-08-27 23:08:16
摘自14.5的help说明
/*******************************************************************/
    UDF for specifying an x-momentum source term in a spatially
   dependent porous media
************** ...

我试了以下,还是不行。这个例子是x方向动量的损失,我需要的是质量变化,只不过用到了坐标轴
5楼2014-08-28 10:21:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tank_fire

铜虫 (小有名气)

【答案】应助回帖

要质量变化你可以把每步的网格里的 质量=密度×体积  存进用户自定义内存里C_UDMI(c,t,i)
算每步之间的差值
使用方法参见 3.2.10. User-Defined Memory (UDM) Macros  (ANSYS14.5 HELP-》FLUENT-》 UDF)
6楼2014-08-28 23:12:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

UDorcer

新虫 (小有名气)

内容已删除
7楼2014-08-29 09:11:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 UDorcer 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 17:36 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:14 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 14:05 by GkCVyODlx4mR
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 14:03 by GkCVyODlx4mR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 13:54 by GkCVyODlx4mR
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 06:00 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见