24小时热门版块排行榜

>论坛更新日志 (4363)
>休闲灌水 (440)
>虫友互识 (379)
>导师招生 (292)
>文献求助 (259)
>博后之家 (168)
>论文投稿 (145)
>招聘信息布告栏 (128)
>硕博家园 (61)
>学术会议 (50)
>考博 (44)
>基金申请 (36)
>育儿交流 (32)
>教师之家 (32)
>有机交流 (27)
>找工作 (26)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

董啊

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 112.9
帖子: 54
在线: 7.8小时
虫号: 26497402

[交流] fluent中梯度UDF源项

我想编写源项𝑆=∇(𝜌𝑣(ℎ_𝑙−&#119867

)，fluent里好像没有对一堆式子求梯度的宏，那我想用C_UDSI=𝜌𝑣(ℎ_𝑙−&#119867

，再用UDM=C_UDSI_G，最后用源项S=UDM，不知行得通吗？或者直接把这个式子开出来，但这样又麻烦又难写，且不知正确。如果可行，UDS如何设置。请各位前辈畅言，下面是我的UDF：
DEFINE_ADJUST(adjust_gradient, domain)
{
Thread *t;
real a,b,d;
cell_t c;
face_t f;
domain = Get_Domain(domain_ID);
/* Fill UDS with the variable. */
thread_loop_c (t,domain)
   {
   begin_c_loop (c,t)
      {
         a=C_U(c,t);
         b=C_V(c,t);
         d=sqrt(a*a+b*b);
      C_UDSI(c,t,0) = C_R(c,t)*d*((C_H(c,t)-(1-C_LIQF(c,t))*g*C_T(c,t))/C_LIQF(c,t)-C_H(c,t));
      }
   end_c_loop (c,t)
   }
  }
#include"udf.h"

DEFINE_SOURCE(guye_source,cell,thread,dS,eqn)
{
real source;
cell_t c;
Thread *t;
C_UDMI(c,t,0)=NV_MAG(C_UDSI_G(c,t,0));
source=C_UDMI(c,t,0);
dS[eqn]=0.0;
return source;
}