24小时热门版块排行榜    

查看: 2951  |  回复: 12

Iwinok

新虫 (初入文坛)

[交流] Fluent中UDF自定义了一个热流密度 已有4人参与

想在一个空壁上加载正弦变化的热流密度,用UDF编了加载之后,一初始化就报错。错误提示:
Fluent received fatal signal (Access_VIOLATION)
我写的程序:
#include "udf.h"
#define pi 3.1415926
#define F1 20000.0  
DEFINE_PROFILE(heat_flux,t,i)
{
  face_t f;
  real time=CURRENT_TIME;
  real temp=F_T(f,t);
  real u=0.2497+0.001471*temp-2.097e-6*pow(temp,2.0);
  begin_f_loop(f,t)
{
  F_PROFILE(f,t,i)=(F1+F1*sin(2*pi*F1*time))*u;
}
end_f_loop(f,t)
}
请大家帮忙看看,是我的程序有问题吗?
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

5超7

至尊木虫 (著名写手)

菜鸟


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by Iwinok at 2013-01-16 13:51:20
u是温度相关的一个系数

既然这样的话,u应该放到loop里面
我要奋斗到底
6楼2013-01-16 14:09:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

5超7

至尊木虫 (著名写手)

菜鸟


小木虫: 金币+0.5, 给个红包,谢谢回帖
u是什么东西?

[ Last edited by 5超7 on 2013-1-16 at 11:07 ]
我要奋斗到底
2楼2013-01-16 11:04:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meercy

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我看到两个问题,不知道有没有用:
1."2.097e"中e是什么?
2.我觉得real temp=F_T(f,t);应该放在loop里面
3楼2013-01-16 12:50:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Iwinok

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by meercy at 2013-01-16 12:50:44
我看到两个问题,不知道有没有用:
1."2.097e"中e是什么?
2.我觉得real temp=F_T(f,t);应该放在loop里面

谢谢啊,这个确实是要方在loop里。那time也是要放在里面吧
4楼2013-01-16 13:50:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Iwinok

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 5超7 at 2013-01-16 11:04:40
u是什么东西?

u是温度相关的一个系数
5楼2013-01-16 13:51:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Iwinok

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by 5超7 at 2013-01-16 14:09:52
既然这样的话,u应该放到loop里面...

是不是热流密度只能用解释型的,能用编译的吗?
7楼2013-01-17 16:37:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

5超7

至尊木虫 (著名写手)

菜鸟


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by Iwinok at 2013-01-17 16:37:13
是不是热流密度只能用解释型的,能用编译的吗?...

你试试不就知道了吗?
我要奋斗到底
8楼2013-01-17 17:00:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Iwinok

新虫 (初入文坛)

修改了下,现在能运行了,但是怎么只加载了一个点。
UDF如下:
#include "udf.h"
#define e 2.718281828459
#define pi 3.1415926
#define F1 15000.
#define F2 20000.
#define r1 0.0015
#define r2 0.001
DEFINE_PROFILE(heat_flux,t,i)
{
  face_t f;
  real time=CURRENT_TIME;
  real u;
  real F3;
  real temp=F_T(f,t);
if(temp<400.)
{
F_PROFILE(f,t,i)=1.0e7;
}
else
begin_f_loop(f,t)
{
  u=0.2497+0.001471*temp-2.097e-6*pow(temp,2.0);
  F3=F1+F2*sin(2*pi*F2*time);
  if(F3<0.)
  {
    F3=0.;
  }
  F_PROFILE(f,t,i)=4*pi*F2*F3*u*(pow(r1,3)-pow(r2,3))/(3*pow(r1,2)*(pow(r1,2)-pow(r2,2)));
}
end_f_loop(f,t)
}

运行后就一处有热流密度

9楼2013-01-18 14:56:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Iwinok

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by 5超7 at 2013-01-17 17:00:37
你试试不就知道了吗?...

能帮我看看吗
10楼2013-01-18 15:11:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Iwinok 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 311求调剂 +5 冬十三 2026-03-15 5/250 2026-03-15 18:38 by 无际的草原
[考研] 304求调剂 +5 小熊joy 2026-03-14 5/250 2026-03-14 21:07 by peike
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 0703化学求调剂 +5 很老实人 2026-03-09 5/250 2026-03-14 02:57 by JourneyLucky
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 泣血叩求调剂恩,愿以丹心报师恩 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
请填处理意见