24小时热门版块排行榜

返回列表

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

pancard

铁虫 (小有名气)

应助: 4 (幼儿园)
金币: 4
红花: 2
帖子: 61
在线: 33.8小时
虫号: 1271842
注册: 2011-04-20
性别: GG
专业: 零件成形制造

[求助] VOF体积分数的梯度为什么为0呢？？

分别使用如下两个UDF计算VOF中体积分数的梯度，两只情况下梯度值都是0？
请高手指教，UDF是否有什么问题？
cas文件及UDF在附件中。
UDF2:

#include "udf.h"

DEFINE_ON_DEMAND(on_demand_calc)
{
Domain *subdomain;
/* declare domain pointer since it is not passed as an argument to the DEFINE macro  */
Thread *cell_thread;
cell_t cell;
float a1,a2,x[2];
subdomain=Get_Domain(2); /* returns phase with ID=2 domain pointer£¬meanning primary phase*/
/* Loop over all cell threads in the domain */
thread_loop_c(cell_thread,subdomain)
{
   /* Loop over all cells  */
   begin_c_loop(cell,cell_thread)
   {
      C_UDSI(cell,cell_thread,0)=C_VOF(cell,cell_thread);
}
   end_c_loop(cell,cell_thread)
   begin_c_loop(cell,cell_thread)
   {
      C_CENTROID(x,cell,cell_thread);
C_UDMI(cell,cell_thread,0)=NV_MAG(C_UDSI_G(cell,cell_thread,0));
a1=C_VOF(cell,cell_thread);
a2=C_UDMI(cell,cell_thread,0);
printf("location:%f,%f,vof:%f,gradient:%f\n",x[0],x[1],a1,a2);
   }
   end_c_loop(cell,cell_thread)
}
}


UDF3：
#include "udf.h"
DEFINE_ON_DEMAND(grad)
{
  Domain *mixture_domain, *sub_domain;
  Thread *cell_thread;
  cell_t cell;
  int phase_domain_index;
  float a1,a2,x[2];
  mixture_domain=Get_Domain(1);

sub_domain_loop(sub_domain,mixture_domain,phase_domain_index)
  {
   /* loop if primary phase  */
/* if (DOMAIN_ID(sub_domain) == 2)    */

   /* loop over all cell threads in the primary phase domain  */
   thread_loop_c(cell_thread,sub_domain)
   {
      /* loop over all cells in primary phase cell threads  */
      begin_c_loop_all(cell,cell_thread)
      {
      C_CENTROID(x,cell,cell_thread);
   C_UDSI(cell,cell_thread,0)=C_VOF(cell,cell_thread);
   C_UDMI(cell,cell_thread,0)=NV_MAG(C_UDSI_G(cell,cell_thread,0));
   a1=C_VOF(cell,cell_thread);
   a2=C_UDMI(cell,cell_thread,0);
   printf("location:%f,%f,vof:%f,gradient:%f\n",x[0],x[1],a1,a2);
      }
      end_c_loop_all(cell,cell_thread)
   }

}
}
水-空气相分布

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : FLUENTWORK.rar

2012-06-06 21:35:30, 35.79 K

欢迎电磁仿真、焊接、快速成型的朋友们来交流

1楼 2012-06-06 21:37:37

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

zlx198200

新虫 (小有名气)

应助: 4 (幼儿园)
金币: 878.3
散金: 759
红花: 5
帖子: 103
在线: 106.1小时
虫号: 759444
注册: 2009-04-28
专业: 多相流热物理学

【答案】应助回帖

★ ★
xiegangmai: 金币+2, 谢谢参与！ 2012-06-07 23:40:12

#define BOUNDARY_ID_fluid 2 /*全局变量，对整个流体区域*/

/*************/
DEFINE_ON_DEMAND(on_demand_calc)
{
Domain *subdomain;
Thread *t;
cell_t c;
float a1,a2,x[ND_ND],xx,yy;
FILE *fpc;

  domain=Get_Domain(1); /*基本相是不是默认1*/
   t=Lookup_Thread(domain,2)  /*混合相指针*/
   t1=THREAD_SUB_THREAD(t,0); /*一相*/
   t2=THREAD_SUB_THREAD(t,1);/*另一项*/

   /* Loop over all cell threads in the domain */
thread_loop_c(c,t)
{
   C_CENTROID(x,c,t);
   xx=x[0];
   yy=x[1];

   C_UDSI(c,t,0)=C_VOF(c,t);
   C_UDMI(c,t,0)=NV_MAG(C_UDSI_G(c,t,0));
   a1=C_UDSI(c,t,0);

   fpc=fopen("gradient.txt","a"