应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » VOF体积分数的梯度为什么为0呢？？
111/2返回列表12下一页
查看: 2360  |  回复: 10
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者pancard将赠送您 1 个金币

pancard

铁虫 (小有名气)

[求助] 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)
     }

   }
}
水-空气相分布
回复此楼

» 本帖附件资源列表

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

» 收录本帖的淘贴专辑推荐

udf相关 传热传质

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
欢迎电磁仿真、焊接、快速成型的朋友们来交流
1楼2012-06-06 21:37:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gyctju

金虫 (正式写手)
★
臭水沟: 金币+1, 谢谢交流~~ 2012-06-07 19:22:06
看得出楼主的意图是用C_UDSI_G来求vof的梯度，用define_on_demand是用来做后处理吧，如果是，建议用tec来做
赞一下 回复此楼
2楼2012-06-07 10:05:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pancard

铁虫 (小有名气)
自己顶
回复此楼
欢迎电磁仿真、焊接、快速成型的朋友们来交流
3楼2012-06-07 14:55:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pancard

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by gyctju at 2012-06-07 10:05:06
看得出楼主的意图是用C_UDSI_G来求vof的梯度，用define_on_demand是用来做后处理吧，如果是，建议用tec来做

并不是用来作后处理，主要是define_on_demand调用方便，用来调试啊。
赞一下 回复此楼
欢迎电磁仿真、焊接、快速成型的朋友们来交流
4楼2012-06-07 14:56:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pancard

铁虫 (小有名气)
再顶起：
读了两年博，发现中国学术不行的最重要一点就是：导师不搞学术。
理由有：
1 学术不能全靠博士，学生是流水的兵，不会形成有效的积累；
2 没有积累，新来的博士都是从0开始，没有“巨人的肩膀”可站；
3 导师不搞学术，导师学术水平下降，从而导致外行（导师）领导内行（博士），瞎指挥。
赞一下 回复此楼
欢迎电磁仿真、焊接、快速成型的朋友们来交流
5楼2012-06-07 16:19:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pancard

铁虫 (小有名气)
发错地了，
回复此楼
欢迎电磁仿真、焊接、快速成型的朋友们来交流
6楼2012-06-07 16:20:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zlx198200

新虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢参与！ 2012-06-07 23:39:56
你网格也太粗了吧，界面这么不清晰，头文件加了么？
类似sg_mphase.h
赞一下 回复此楼
7楼2012-06-07 21:55:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zlx198200

新虫 (小有名气)

【答案】应助回帖

★ ★
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";
         fprinth(fpc,"%f          %f          %f\n", xx,yy,a1);
         fclose(fpc);
   }
   end_c_loop(c,t)
赞一下 回复此楼
8楼2012-06-07 22:29:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zlx198200

新虫 (小有名气)
★
xiegangmai: 金币+1, 谢谢参与！ 2012-06-07 23:40:21
头文件还有好几个，你自己看看帮助吧，你这个UDF我觉得自己都没写清楚，SUB,SUPER这些不要瞎用，
赞一下 回复此楼
9楼2012-06-07 22:31:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pancard

铁虫 (小有名气)
自己解决了，过程很复杂，详见
http://wenku.baidu.com/view/1018444f581b6bd97e19ea3c.html
赞一下 回复此楼
欢迎电磁仿真、焊接、快速成型的朋友们来交流
10楼2014-04-24 21:01:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pancard 的主题更新
111/2返回列表12下一页
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 夏至，要不要硕博联谊 +6 我是王小帅 2024-06-21 8/400 2024-06-26 21:05 by zxw3578
[有机交流] 做什么表征可以检测塑料中碳氢氧的含量 5+3 pzr的sci之路 2024-06-25 7/350 2024-06-26 19:32 by wrgeng
[硕博家园] 回家两天，不想呆了 +7 368ghnf 2024-06-22 13/650 2024-06-26 19:07 by ZZZemmm
[有机交流] 有机物的酸酐如何准确表征 10+3 方酱 2024-06-24 7/350 2024-06-26 17:53 by 宁静远行
[考研] 刚当完兵回来想考研 +7 五官在线 2024-06-23 18/900 2024-06-26 17:05 by 1158057902
[硕博家园] 联培博士文章第一单位署名问题交流~ +9 橙成成c 2024-06-23 21/1050 2024-06-26 16:45 by 橙成成c
[找工作] 老博士选高校 5+5 woodhead8386 2024-06-24 8/400 2024-06-26 16:35 by 大大熊呀
[基金申请] 博后面上今天下午会公布吗？大家有无消息？ +9 地理学1995 2024-06-24 13/650 2024-06-26 16:02 by kyukitu
[考博] 申博好难 +6 自强不息a?a 2024-06-24 6/300 2024-06-25 23:02 by 考研吒儿
[考博] 没读上博，好焦虑！ +6 wangzhe_bs 2024-06-24 8/400 2024-06-25 21:15 by wangzhe_bs
[基金申请] 焦虑没有毛线用，默默前行是王道！ +4 漠上藜梭 2024-06-24 8/400 2024-06-25 14:32 by 漠上藜梭
[有机交流] 高温酯化反应喷料 20+4 hl24678 2024-06-21 7/350 2024-06-25 09:45 by ghj-hebust
[金属] EBSD的解析率只有10% +3 wallace6666 2024-06-20 7/350 2024-06-24 16:52 by wallace6666
[基金申请] 说博后基金7月出的真打电话了吗？ +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[基金申请] 青年和面上，哪个上会难度更大 +12 今晚推荐22 2024-06-21 18/900 2024-06-24 11:08 by 半简体
[有机交流] 三颈瓶加热 5+4 wzjwx888 2024-06-20 4/200 2024-06-24 10:11 by w87437985
[论文投稿] 论文提交二审还有三天就三个月了，连续问了编辑部几次 10+3 大王叫我来寻山� 2024-06-22 9/450 2024-06-24 08:50 by 大王叫我来寻山�
[公派出国] 博士csc联培会看重第一学历学校层次吗 +4 也就这样 2024-06-23 4/200 2024-06-24 08:18 by 晓目崇
[博后之家] 在国内某高校做全职博士后2年，现在找到新的单位，出站或退站对新工作有什么影响？ +10 nxplfcc 2024-06-20 10/500 2024-06-22 07:52 by 徐长安
[基金申请] 听大佬说今年信息口本子数量大幅增加？ +8 wutzxt 2024-06-21 9/450 2024-06-21 19:58 by wutzxt
信息提示
关闭
请填处理意见
关闭