24小时热门版块排行榜    

查看: 2399  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

拾贝壳的少年

新虫 (初入文坛)


[交流] fluent UDF DEFINE_SOURCE 问题

各位虫友,大家好。
我现在的研究课题是一模拟电阻加热炉的升温过程,加热元件是电加热棒,现在已经打开辐射模型,在加热棒的加热功率随炉子里的热电偶温度的变化而变化。现在在加热棒的固体区域设置能量源项。但现在遇到的问题是DEFINE_SOURCE(name,c,t,ds,eqn)定义中的 SΦ=A+BΦ,dS[eqn]=dS/dΦ,这个自变量Φ可不可以是另一计算域的变量,即能量源项中的体积热流密度是另一个计算域中某点温度的函数(反馈调节),如果体积热流密度是该计算域中温度的函数的话,udf 比较好写,但是如果是另一个计算域某一点温度的函数的话不知可不可以实现?
现在课题卡在这了,求各位大神不吝赐教!
#include"udf.h"
DEFINE_SOURCE(energy_source,c,t,dS,eqn)
{
  Domain *domain=Get_Domain(37);//获得炉子中监控位置dormain的ID值。
  real temp_ave;               
  real temp_sum=0;              
  int  count=0;
  cell_t c1;
  Thread *thread_thermocouple;  
  real source;

//下面这个循环是为了得到监控区域的cell的平均温度
thread_loop_c(thread_thermocouple,domain)
{
  begin_c_loop(c1,thread_thermocouple)
    {
      ++count;
      temp_sum+=C_T_M1(c1,thread_thermocouple);
     }
  end_c_loop(c1,thread_thermocouple)
}
  temp_ave=temp_sum/count;//监控位置的平均温度

//如果温度监控位置的平均温度小于973K,则体积热流密度为3.468*10^7
if(temp_ave<973)
{
   source=3.468*pow(10,7);
   dS[eqn]=0;
  }
/*如果温度监控位置的平均温度大于973K小于1023K,则体积热流密度为-884.56*temp_ave+904904.88 */
  else if(temp_ave<1023)
{
   source=-884.56*temp_ave+904904.88;
   dS[eqn]=-884.56;
  }
//如果温度监控位置的平均温度大于1023K,则体积热流密度为
  else
  {
    source=0;
    dS[eqn]=0;
   }
  return source;
}
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

北落北海

新虫 (初入文坛)



拾贝壳的少年(金币+1): 谢谢参与
楼主请问你解决了么   求指教~
5楼2017-10-28 13:49:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

拾贝壳的少年

新虫 (初入文坛)


金币不多,请各位大神帮帮忙。
3楼2017-07-07 09:28:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

4dff

新虫 (小有名气)



拾贝壳的少年(金币+1): 谢谢参与
http://muchong.com/t-11125735-1可以参考一下这个
6楼2018-12-10 09:54:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赤道以北晴

新虫 (初入文坛)



拾贝壳的少年(金币+1): 谢谢参与
引用回帖:
3楼: Originally posted by 拾贝壳的少年 at 2017-07-07 09:28:25
金币不多,请各位大神帮帮忙。

楼主,您可以指点一下吗,用了和您类似的udf 计算,编译成功,但是一开始就报错了
7楼2019-04-03 20:24:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2017-07-13 02:45   回复  
拾贝壳的少年(金币+1): 谢谢参与
发自小木虫IOS客户端
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +7 ziyangfang 2026-08-25 10/500 2026-09-01 09:32 by qvbaoxiao
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 面上函评意见出来了,像什么等级? 20+3 Tsingking1 2026-08-27 14/700 2026-09-01 08:29 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +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
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +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
[基金申请] 怎么查啊 +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
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国合里面能看到了 +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
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见