24小时热门版块排行榜    

查看: 1741  |  回复: 4

8厘米阳光

新虫 (初入文坛)

[求助] 新手求助!!!编写一个UDF,能够导入fluent,但开始计算就出现错误 已有1人参与

新手求助!!!编写一个UDF,式子中含有根号,能够导入fluent,但开始计算就出现错误
UDF如下:
/****************************************************************************************
  UDF that computes specific heat and sets the sensible enthalpy
  to the referenced value
***************************************************************************************/
#include"udf.h"
DEFINE_SPECIFIC_HEAT(cp,T,Tref,h,yi)
{
        float cp;
        if (297<=T<=302)
           cp=814*sqrt(-T*T+599*T-89694)+4065;
        else
           cp=4065;
           *h=cp*(T-Tref);
        return cp;
}
出现的错误是:floating point exception
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sanguteduo

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
你把302代入根号下试试,结果是0,必然溢出,可以把89694改成89693.99么,哈哈哈,应该不影响,只是个比导热,且不在大比热区。
信仁立
2楼2015-08-09 16:05:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

8厘米阳光

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by sanguteduo at 2015-08-09 16:05:13
你把302代入根号下试试,结果是0,必然溢出,可以把89694改成89693.99么,哈哈哈,应该不影响,只是个比导热,且不在大比热区。

谢谢~试了一下 还是出现相同的错误
3楼2015-08-09 17:14:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CFD_light

铁虫 (正式写手)

呵呵,这个要很小心才不出错的
模拟的最高境界:同形不同理
4楼2015-08-10 07:43:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

8厘米阳光

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by CFD_light at 2015-08-10 07:43:31
呵呵,这个要很小心才不出错的

我现在换成正弦的了 计算结果跟以前别人做的有出入 我不知道在编正弦的UDF的时候 是不是有什么要注意的
5楼2015-09-14 11:33:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 8厘米阳光 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-31 02:13 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-31 01:16 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-31 01:15 by hZiFeudZyoGR
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-31 00:32 by hZiFeudZyoGR
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 12/600 2026-08-31 00:31 by hZiFeudZyoGR
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:21 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 4FFAWE8HcgUD 2026-08-29 3/150 2026-08-30 20:10 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:57 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:48 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:46 by hZiFeudZyoGR
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jCd0dEvKHShX 2026-08-29 6/300 2026-08-30 19:13 by hZiFeudZyoGR
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 18:20 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 国际合作可查了,中了面上 (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 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见