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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » Fluent中UDF自定义了一个热流密度
131/2返回列表12下一页
查看: 3475  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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)
}
请大家帮忙看看，是我的程序有问题吗？
回复此楼

» 猜你喜欢

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

5超7

至尊木虫 (著名写手)

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

既然这样的话，u应该放到loop里面
赞一下(1人) 回复此楼
我要奋斗到底
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 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 05:50 by 5BDX0d0WFp7t
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:49 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:01 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:25 by mibUvS8DDCwf
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:21 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 23:45 by vZfe6xYu34yj
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:33 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 22:45 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 by vZfe6xYu34yj
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
关闭
请填处理意见
关闭