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

24小时热门版块排行榜    

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

铁虫 (小有名气)
自己解决了，过程很复杂，详见
http://wenku.baidu.com/view/1018444f581b6bd97e19ea3c.html
赞一下 回复此楼
欢迎电磁仿真、焊接、快速成型的朋友们来交流
10楼2014-04-24 21:01:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

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

铁虫 (小有名气)
引用回帖:
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的回帖
查看全部 11 个回答
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 2025考博 +7 自强不息a?a 2024-06-15 12/600 2024-06-17 22:51 by 自强不息a?a
[基金申请] 太卷了 +12 shuigubio 2024-06-17 23/1150 2024-06-17 22:38 by master_110
[硕博家园] 关于硕博连读的一些疑问？ +4 Lwenter 2024-06-14 5/250 2024-06-17 20:42 by Lwenter
[基金申请] 化学会年会改成一年开一次了？ +7 babu2015 2024-06-12 15/750 2024-06-17 20:28 by 落笔画枫
[考博] 这个博士要读吗 +9 Sea Breeze 2024-06-16 17/850 2024-06-17 19:50 by zeolitess
[基金申请] 最新消息：2024国自然下载文件名变了 +15 dream200606 2024-06-16 28/1400 2024-06-17 19:04 by 6543yes
[论文投稿] 审稿问题：为什么荧光激发波长和紫外吸收波长差的大？ 10+5 sdawege 2024-06-14 10/500 2024-06-17 18:54 by HH-探针
[有机交流] 车间生产，真空度很高，温度很高，但减压蒸馏速度很慢。 10+12 召唤鬼泣lL 2024-06-13 38/1900 2024-06-17 16:46 by 科研虫子小助手
[基金申请] Nature 11日发文，中国著名学者们称造假迫不得已 +6 babu2015 2024-06-14 6/300 2024-06-16 18:37 by 喜欢和一氧化二�
[教师之家] 每次骚扰女学生的都是院系领导，而不是普通教师，小编们要注意措辞正确 +9 zju2000 2024-06-15 11/550 2024-06-16 14:49 by appleapple2
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
[考博] 上海交大招收材料化学方向科研助理/“申请考核”博士（已招满） +3 灵梦and紫 2024-06-12 4/200 2024-06-15 20:58 by 1822836277
[论文投稿] 投稿时忘记修改一作 +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
[基金申请] 工材E口JQ有消息了吗 +4 babyduck 2024-06-11 4/200 2024-06-14 17:23 by firepick
[基金申请] E12面上申请 +4 汉风之遗 2024-06-13 4/200 2024-06-14 15:28 by 天外飞去来
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[基金申请] 75批博后基金 +10 kyukitu 2024-06-13 13/650 2024-06-14 10:31 by kyukitu
[论文投稿] 最近写了一篇控制优化领域的文章，可以投哪里啊？有没有水一些的期刊推荐 +7 香瓜木香 2024-06-12 13/650 2024-06-14 07:05 by 香瓜木香
[基金申请] 博后特助这周出结果吗？往年都是啥时候啊？ +13 jsqy 2024-06-12 17/850 2024-06-12 19:55 by Lynn212
信息提示
关闭
请填处理意见
关闭