24小时热门版块排行榜    

查看: 2423  |  回复: 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的回帖

拾贝壳的少年

新虫 (初入文坛)


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

北落北海

新虫 (初入文坛)



拾贝壳的少年(金币+1): 谢谢参与
楼主请问你解决了么   求指教~
5楼2017-10-28 13:49:27
已阅   回复此楼   关注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客户端
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-17 02:02 by Tql5LQhh5rLK
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 i3KBRD19ERkG 2026-09-16 3/150 2026-09-17 00:50 by Tql5LQhh5rLK
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-14 7/350 2026-09-16 19:27 by rRvxEj3zorCy
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-16 18:00 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:15 by i3KBRD19ERkG
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 17:12 by i3KBRD19ERkG
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-16 16:24 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 14:39 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 07:48 by i0IGWqQuchto
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 04:51 by LH5NkK5Ud4bh
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 02:23 by DgNGHc3h5tPl
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:26 by zNcEhTcH7ihq
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见