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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » Fluent中UDF自定义了一个热流密度
131/2返回列表12下一页
查看: 3460  |  回复: 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一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:52 by o3EMxHUAQDtS
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:23 by o3EMxHUAQDtS
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗？ +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (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
信息提示
关闭
请填处理意见
关闭