应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » UDF中DEFINE_ADJUST宏中使用C_UDSI_G
371/4返回列表1234下一页
查看: 10155  |  回复: 36
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

duandiancqu

金虫 (初入文坛)

[求助] UDF中DEFINE_ADJUST宏中使用C_UDSI_G 已有5人参与

我在用DEFINE_ADJUST宏时，只要加上C_UDSI_G(c,t,0），就提示错误，
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.

将C_UDSI_G()替换成常数时，程序能够执行，是不是在这种情况下有什么要求啊，那位大神给我解释下，很纠结。
程序代码如下：
DEFINE_ADJUST(mhd_adjust,d)
{
        Thread *t;
        cell_t c;

        if (! Data_Valid_P())
        return;

        thread_loop_c(t,d)
        {
                begin_c_loop(c,t)
                {
                                               
                        C_UDMI(c,t,0)=-1*C_UDSI_DIFF(c,t,0)*C_UDSI_G(c,t,0)[0];
                        C_UDMI(c,t,1)=-1*C_UDSI_DIFF(c,t,0)*C_UDSI_G(c,t,0)[1];
                        C_UDMI(c,t,2)=-1*C_UDSI_DIFF(c,t,0)*C_UDSI_G(c,t,0)[2];
                }
                end_c_loop(c,t)
        }
}
回复此楼

» 收录本帖的淘帖专辑推荐

等离子体相关

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-07-02 22:09:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

759599290

银虫 (小有名气)

【答案】应助回帖

最近我也测试了使用DEFINE_ADJUST宏来获得某一变量的二阶导数。心得如下：
1、梯度变量在FLUENT计算过程中是暂存变量，计算完成后将自动释放。如果要获取，需使用如下方法强制求解器保存梯度变量已供后用：
在FLUENT主窗中输入：solve--set--expert，出现的是否保存梯度变量询问语句回答为yes
2、梯度变量在首次计算中是不存在的。即使在计算之前，UDS已经使用一定方法进行过初始化，UDS_G仍将是空值，在第一步计算时内存中不存在这一变量，因此无论如何设置，第一步计算都不能调用C_UDSI_G宏。计算进行了一步或更多步之后才可以。
3、C_UDSI_G调出的变量并不一定是内部计算使用的梯度值。在某些复杂边界处，Fluent内部有其他未公开的处理函数，C_UDSI_G只在计算区域中的绝大部分里非常接近内部计算使用的梯度值。
仅个人看法，未必正确，盼有所助。
赞一下(9人) 回复此楼
中国热管理协会特约专家，中国热设计网高级顾问，东莞鸿艺电子技术总监
23楼2014-01-13 08:50:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

759599290

银虫 (小有名气)

【答案】应助回帖

★ ★
xiegangmai: 金币+2, 谢谢参与！ 2013-03-30 22:27:23
建议使用NNULL函数判断线程中是否存有C_UDSI_G(c,t,i)数据。data_valid_p和NNULL函数判断的侧重点可能并不一致。
赞一下(2人) 回复此楼
中国热管理协会特约专家，中国热设计网高级顾问，东莞鸿艺电子技术总监
18楼2013-03-30 16:57:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tommyan

新虫 (初入文坛)
这个代码貌似是非法访问。
C_UDSI_G(c,t,0) 返回的是一个矢量，不能直接用在乘法表达式。
可以尝试：
求矢量模：NV_MAG(C_UDSI_G(c,t,0))再乘，
或者
real g[ND_ND];
g=C_UDSI_G(c,t,0);
使用g[0],g[1] g[2]
赞一下(2人) 回复此楼
25楼2014-01-30 18:25:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

gyctju

金虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
臭水沟: 金币+1, 谢谢交流~~ 2012-07-04 20:39:50
可以尝试在不链接udf的情况下先计算几步，使得uds0的值合理。
赞一下 回复此楼
2楼2012-07-04 12:31:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duandiancqu

金虫 (初入文坛)
引用回帖:
2楼: Originally posted by gyctju at 2012-07-04 12:31:48
可以尝试在不链接udf的情况下先计算几步，使得uds0的值合理。

我试过了，不行。。。
回复此楼
3楼2012-07-04 14:45:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lucklhz

木虫 (正式写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
臭水沟: 金币+1, 谢谢交流~~ 2012-07-04 20:40:06
楼主，看你程序，你是想取你的UDS自定义标量uds0在空间个三个方向的梯度吗？还有，这种表达方式C_UDSI_G(c,t,0)[1]，是不是有问题？C_UDSI_G(c,t,i)默认的在三个方向上的梯度是一样的
赞一下 回复此楼
4楼2012-07-04 19:43:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duandiancqu

金虫 (初入文坛)
引用回帖:
4楼: Originally posted by lucklhz at 2012-07-04 19:43:32
楼主，看你程序，你是想取你的UDS自定义标量uds0在空间个三个方向的梯度吗？还有，这种表达方式C_UDSI_G(c,t,0)，是不是有问题？C_UDSI_G(c,t,i)默认的在三个方向上的梯度是一样的

我是想取载个方向上的梯度，而且C_UDSI_G（c,t,0）[1]没有问题吧,帮助上就有这个啊，C_UDSI_G(c,t,0)，i=1,2,3的
赞一下 回复此楼
5楼2012-07-04 22:42:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yulinjun1986

新虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
臭水沟: 金币+1, 谢谢交流~~ 2012-07-05 17:20:16
我想问的是，你之前编写uds了吗？是不是编写的不对呀？如果没有编写的话，那可能就是需要编写uds了……
赞一下 回复此楼
6楼2012-07-05 16:26:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duandiancqu

金虫 (初入文坛)
引用回帖:
6楼: Originally posted by yulinjun1986 at 2012-07-05 16:26:29
我想问的是，你之前编写uds了吗？是不是编写的不对呀？如果没有编写的话，那可能就是需要编写uds了……

写了的，就是一个电势方程，求出电势，再利用欧姆定律求出电流密度
赞一下 回复此楼
7楼2012-07-05 21:21:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duandiancqu

金虫 (初入文坛)
引用回帖:
6楼: Originally posted by yulinjun1986 at 2012-07-05 16:26:29
我想问的是，你之前编写uds了吗？是不是编写的不对呀？如果没有编写的话，那可能就是需要编写uds了……

我现在把那个去掉了，可以计算，能量方程就是不收敛。。。。
赞一下 回复此楼
8楼2012-07-05 21:25:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yulinjun1986

新虫 (小有名气)

【答案】应助回帖

★ ★
xiegangmai: 金币+2, 谢谢参与！ 2012-07-06 22:25:11
引用回帖:
8楼: Originally posted by duandiancqu at 2012-07-05 21:25:21
我现在把那个去掉了，可以计算，能量方程就是不收敛。。。。...

你用的是非稳还是稳态？要是非稳的话，试试把时间步长调小些。另外就是要调整松弛因子了。你试试吧，如果不加uds的话，应该没问题才对。我的udf也遇到了你的问题，一加入和质量相关的源相就提示这个错误，我也相当无奈，正在求解中…希望大家相互帮助吧。
赞一下 回复此楼
9楼2012-07-06 09:12:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duandiancqu

金虫 (初入文坛)
引用回帖:
9楼: Originally posted by yulinjun1986 at 2012-07-06 09:12:43
你用的是非稳还是稳态？要是非稳的话，试试把时间步长调小些。另外就是要调整松弛因子了。你试试吧，如果不加uds的话，应该没问题才对。我的udf也遇到了你的问题，一加入和质量相关的源相就提示这个错误，我也相当 ...

谢谢楼上了，我算的是稳态的
赞一下 回复此楼
10楼2012-07-06 09:17:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 duandiancqu 的主题更新
371/4返回列表1234下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 h4CP7TrQR8Lg 2026-08-22 4/200 2026-08-23 20:46 by 0p9y1isJp2jM
[基金申请] 让我中一个面上吧！ +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 19:43 by lXgU0RvNTDC9
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 19:34 by lXgU0RvNTDC9
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 18:58 by lXgU0RvNTDC9
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 18:46 by lXgU0RvNTDC9
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 h4CP7TrQR8Lg 2026-08-22 3/150 2026-08-23 16:30 by lXgU0RvNTDC9
[基金申请] 什么时候开奖？ +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
信息提示
关闭
请填处理意见
关闭