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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » VOF体积分数的梯度为什么为0呢？？
71/1返回列表
查看: 2354  |  回复: 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的回帖

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的回帖

pancard

铁虫 (小有名气)
自己解决了，过程很复杂，详见
http://wenku.baidu.com/view/1018444f581b6bd97e19ea3c.html
赞一下 回复此楼
欢迎电磁仿真、焊接、快速成型的朋友们来交流
10楼2014-04-24 21:01:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pancard 的主题更新
71/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 初始合伙人来啦！（生物试剂耗材标准品） +13 欢快的小科研人 2024-06-15 24/1200 2024-06-18 13:08 by 欢快的小科研人
[找工作] 药学硕士找不到工作，打算去做科研助理了 +6 pom戴墨镜 2024-06-14 19/950 2024-06-18 13:06 by lialee22
[考博] 这个博士要读吗 +9 Sea Breeze 2024-06-16 18/900 2024-06-18 12:57 by Sea Breeze
[教师之家] 请问事业编制和年薪制冲突吗？ +11 ZHONGWU_U 2024-06-14 11/550 2024-06-18 12:26 by zqx1990
[高分子] 烧瓶内合成聚酯 +3 大帝国乐 2024-06-17 7/350 2024-06-18 11:10 by 大帝国乐
[有机交流] 求教最后一步问题出在哪 20+3 苏州浩恒 2024-06-15 4/200 2024-06-18 08:33 by like_2005
[基金申请] F口信息学部拿面上，大概需要什么样的成果 +7 _奋黎_ 2024-06-16 12/600 2024-06-18 07:41 by ld_cauchy
[公派出国] 去英国的小伙伴儿都在哪儿租的房子呀？ +5 65syn 2024-06-14 8/400 2024-06-17 19:57 by 爱吃菜的虫
[教师之家] 每次骚扰女学生的都是院系领导，而不是普通教师，小编们要注意措辞正确 +9 zju2000 2024-06-15 11/550 2024-06-16 14:49 by appleapple2
[基金申请] 博后基金，博管会会提前知道消息吗？ +4 yuyiang 2024-06-13 4/200 2024-06-16 11:40 by yangyuzhong4
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[论文投稿] 投了一篇4区的SCI，审稿人一个拒稿，一个小修，编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[考研] 物理化学一对一辅导 +3 林大diao 2024-06-12 5/250 2024-06-14 20:57 by 林大diao
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 75批博后基金 +10 kyukitu 2024-06-13 13/650 2024-06-14 10:31 by kyukitu
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[硕博家园] 机械研究生如何拿到年薪40+w +13 阿巴阿巴哦哦 2024-06-11 15/750 2024-06-13 15:40 by 113745685
[论文投稿] with editor日期变更 +3 慎独的小花卷 2024-06-12 8/400 2024-06-13 11:00 by 慎独的小花卷
[基金申请] 博后特助这周出结果吗？往年都是啥时候啊？ +13 jsqy 2024-06-12 17/850 2024-06-12 19:55 by Lynn212
信息提示
关闭
请填处理意见
关闭