| 查看: 3080 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
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
» 猜你喜欢
085601材料工程找调剂
已经有11人回复
-
总分293求调剂
已经有9人回复
-
334分 一志愿武理 材料求调剂
已经有5人回复
-
327求调剂
已经有6人回复
-
329求调剂
已经有10人回复
-
275求调剂
已经有15人回复
-
291求调剂
已经有6人回复
-
309求调剂
已经有7人回复
-
297求调剂
已经有11人回复
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历
已经有6人回复
classic1573
新虫 (初入文坛)
- 应助: 4 (幼儿园)
- 金币: 5.4
- 散金: 100
- 帖子: 32
- 在线: 70.2小时
- 虫号: 2419446
- 注册: 2013-04-15
- 专业: 内流流体力学
11楼2014-08-11 22:31:07
2楼2012-06-07 10:05:06
4楼2012-06-07 14:56:17
5楼2012-06-07 16:19:49
