24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1245  |  回复: 3
【悬赏金币】回答本帖问题,作者木风10086将赠送您 5 个金币

木风10086

新虫 (小有名气)

[求助] fluent边界附近源项的UDF已有1人参与

在固体wall边界附近的一层单元内添加源项,本意是当单元压力大于某一值、单元上wall的那条面温度大于某一值时就释放源项,代码如下:
可是运行的时候总是整个wall的温度大于设定值时才释放,而不是单元内的与wall对应的那个面的温度大于值就在该单元释放,问题在哪?
C_T到底是整个wall的温度还是单元内那个面的温度?
#include "udf.h"
DEFINE_SOURCE(c_e, cell, thread, dS, eqn)
{
    real source;
  int n;
   double st;
   face_t f;
   Thread *t;
  Domain *domain=Get_Domain(1);
   t=Lookup_Thread(domain,15);  

begin_f_loop(f,t)
{
       st=F_T(f,t);
    if ( C_P(cell,thread) >= 120000. && st>=310.)
    {
      source =1.10921e11;
      dS[eqn] = 0;
     }
    else
    {
      source = 0;
      dS[eqn] = 0;
     }
  return source;
}   
end_f_loop(f,thread)
}
回复此楼
心情不爽!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木风10086

新虫 (小有名气)

上面写错了,是F_T到底是整个wall的温度还是单元内那个面的温度?
心情不爽!
2楼2020-10-31 23:47:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

佳廷主夫

铁杆木虫 (著名写手)

虫虫

【答案】应助回帖

感谢参与,应助指数 +1
应该是你加载在哪个面上就在哪个面上(三维在体上)
源项中好像没有判断边界的语句  自然默认是所有面或体了
才情似水,假期如梦
3楼2020-11-01 11:13:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

木风10086

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 佳廷主夫 at 2020-11-01 11:13:58
应该是你加载在哪个面上就在哪个面上(三维在体上)
源项中好像没有判断边界的语句  自然默认是所有面或体了

t=Lookup_Thread(domain,15);  指向的就是wall边界,现在是想在单元内判断,但它总是判断整个wall。重新改了下循环,貌似是可以达到要求了。仍在确认中。
心情不爽!
4楼2020-11-01 22:57:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 木风10086 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见