版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

pancard

鐵蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 4
紅花: 2
帖子: 61
在線: 33.8小時
蟲號: 1271842
注冊: 2011-04-20
性別: GG
專業(yè): 零件成形制造

[求助] VOF體積分數(shù)的梯度為什么為0呢？？

分別使用如下兩個UDF計算VOF中體積分數(shù)的梯度，兩只情況下梯度值都是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)
   }

}
}
水-空氣相分布

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : FLUENTWORK.rar

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

歡迎電磁仿真、焊接、快速成型的朋友們來交流

1樓 2012-06-06 21:37:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gyctju

金蟲 (正式寫手)

應(yīng)助: 314 (大學(xué)生)
金幣: 942.6
散金: 227
紅花: 48
沙發(fā): 3
帖子: 996
在線: 427.6小時
蟲號: 1569491
注冊: 2012-01-08
專業(yè): 分離過程

★
臭水溝: 金幣+1, 謝謝交流~~ 2012-06-07 19:22:06

看得出樓主的意圖是用C_UDSI_G來求vof的梯度，用define_on_demand是用來做后處理吧，如果是，建議用tec來做

贊一下

回復(fù)此樓

2樓2012-06-07 10:05:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

pancard

鐵蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 4
紅花: 2
帖子: 61
在線: 33.8小時
蟲號: 1271842
注冊: 2011-04-20
性別: GG
專業(yè): 零件成形制造

自己頂

回復(fù)此樓

歡迎電磁仿真、焊接、快速成型的朋友們來交流

3樓2012-06-07 14:55:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

pancard

鐵蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 4
紅花: 2
帖子: 61
在線: 33.8小時
蟲號: 1271842
注冊: 2011-04-20
性別: GG
專業(yè): 零件成形制造

引用回帖:

2樓: Originally posted by gyctju at 2012-06-07 10:05:06
看得出樓主的意圖是用C_UDSI_G來求vof的梯度，用define_on_demand是用來做后處理吧，如果是，建議用tec來做

并不是用來作后處理，主要是define_on_demand調(diào)用方便，用來調(diào)試啊。

贊一下

回復(fù)此樓

歡迎電磁仿真、焊接、快速成型的朋友們來交流

4樓2012-06-07 14:56:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

pancard

鐵蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 4
紅花: 2
帖子: 61
在線: 33.8小時
蟲號: 1271842
注冊: 2011-04-20
性別: GG
專業(yè): 零件成形制造

再頂起：
讀了兩年博，發(fā)現(xiàn)中國學(xué)術(shù)不行的最重要一點就是：導(dǎo)師不搞學(xué)術(shù)。
理由有：
1 學(xué)術(shù)不能全靠博士，學(xué)生是流水的兵，不會形成有效的積累；
2 沒有積累，新來的博士都是從0開始，沒有“巨人的肩膀”可站；
3 導(dǎo)師不搞學(xué)術(shù)，導(dǎo)師學(xué)術(shù)水平下降，從而導(dǎo)致外行（導(dǎo)師）領(lǐng)導(dǎo)內(nèi)行（博士），瞎指揮。

贊一下

回復(fù)此樓

歡迎電磁仿真、焊接、快速成型的朋友們來交流

5樓2012-06-07 16:19:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

pancard

鐵蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 4
紅花: 2
帖子: 61
在線: 33.8小時
蟲號: 1271842
注冊: 2011-04-20
性別: GG
專業(yè): 零件成形制造

發(fā)錯地了，

回復(fù)此樓

歡迎電磁仿真、焊接、快速成型的朋友們來交流

6樓2012-06-07 16:20:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zlx198200

新蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 878.3
散金: 759
紅花: 5
帖子: 103
在線: 106.1小時
蟲號: 759444
注冊: 2009-04-28
專業(yè): 多相流熱物理學(xué)

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝參與！ 2012-06-07 23:39:56

你網(wǎng)格也太粗了吧，界面這么不清晰，頭文件加了么？
類似sg_mphase.h

贊一下

回復(fù)此樓

7樓2012-06-07 21:55:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zlx198200

新蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 878.3
散金: 759
紅花: 5
帖子: 103
在線: 106.1小時
蟲號: 759444
注冊: 2009-04-28
專業(yè): 多相流熱物理學(xué)

【答案】應(yīng)助回帖

★ ★
xiegangmai: 金幣+2, 謝謝參與！ 2012-06-07 23:40:12

#define BOUNDARY_ID_fluid 2 /*全局變量，對整個流體區(qū)域*/

/*************/
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)

贊一下

回復(fù)此樓

8樓2012-06-07 22:29:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zlx198200

新蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 878.3
散金: 759
紅花: 5
帖子: 103
在線: 106.1小時
蟲號: 759444
注冊: 2009-04-28
專業(yè): 多相流熱物理學(xué)

★
xiegangmai: 金幣+1, 謝謝參與！ 2012-06-07 23:40:21

頭文件還有好幾個，你自己看看幫助吧，你這個UDF我覺得自己都沒寫清楚，SUB,SUPER這些不要瞎用，

贊一下

回復(fù)此樓

9樓2012-06-07 22:31:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

pancard

鐵蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 4
紅花: 2
帖子: 61
在線: 33.8小時
蟲號: 1271842
注冊: 2011-04-20
性別: GG
專業(yè): 零件成形制造

自己解決了，過程很復(fù)雜，詳見
http://wenku.baidu.com/view/1018444f581b6bd97e19ea3c.html

贊一下

回復(fù)此樓

歡迎電磁仿真、焊接、快速成型的朋友們來交流

10樓2014-04-24 21:01:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 pancard 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 材料專碩英一數(shù)二306 +4	z1z2z3879 2026-03-18	4/200	2026-03-18 20:00 by 楤哥
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +9	秋有木北 2026-03-14	9/450	2026-03-18 19:49 by macy2011
[考研] 0703化學(xué)調(diào)劑 +7	妮妮ninicgb 2026-03-15	11/550	2026-03-18 19:48 by macy2011
[考研] 0703化學(xué) 305求調(diào)劑 +3	FY_yy 2026-03-14	3/150	2026-03-18 19:40 by macy2011
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	9/450	2026-03-18 15:05 by stone_128
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 331求調(diào)劑（0703有機化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考研] 268求調(diào)劑 +7	好運連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655

24小時熱門版塊排行榜

[求助] VOF體積分數(shù)的梯度為什么為0呢？？

» 本帖附件資源列表

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖