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

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

24小時(shí)熱門版塊排行榜

返回列表

下一頁(yè)

pancard

鐵蟲 (小有名氣)

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

[求助] VOF體積分?jǐn)?shù)的梯度為什么為0呢？？

分別使用如下兩個(gè)UDF計(jì)算VOF中體積分?jǐn)?shù)的梯度，兩只情況下梯度值都是0？
請(qǐng)高手指教，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)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(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小時(shí)
蟲號(hào): 1569491
注冊(cè): 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小時(shí)
蟲號(hào): 1271842
注冊(cè): 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小時(shí)
蟲號(hào): 1271842
注冊(cè): 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小時(shí)
蟲號(hào): 1271842
注冊(cè): 2011-04-20
性別: GG
專業(yè): 零件成形制造

再頂起：
讀了兩年博，發(fā)現(xiàn)中國(guó)學(xué)術(shù)不行的最重要一點(diǎn)就是：導(dǎo)師不搞學(xué)術(shù)。
理由有：
1 學(xué)術(shù)不能全靠博士，學(xué)生是流水的兵，不會(huì)形成有效的積累；
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小時(shí)
蟲號(hào): 1271842
注冊(cè): 2011-04-20
性別: GG
專業(yè): 零件成形制造

發(fā)錯(cuò)地了，

回復(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小時(shí)
蟲號(hào): 759444
注冊(cè): 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小時(shí)
蟲號(hào): 759444
注冊(cè): 2009-04-28
專業(yè): 多相流熱物理學(xué)

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

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

#define BOUNDARY_ID_fluid 2 /*全局變量，對(duì)整個(gè)流體區(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); /*基本相是不是默認(rèn)1*/
   t=Lookup_Thread(domain,2)  /*混合相指針*/
   t1=THREAD_SUB_THREAD(t,0); /*一相*/
   t2=THREAD_SUB_THREAD(t,1);/*另一項(xiàng)*/

   /* 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小時(shí)
蟲號(hào): 759444
注冊(cè): 2009-04-28
專業(yè): 多相流熱物理學(xué)

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

頭文件還有好幾個(gè)，你自己看看幫助吧，你這個(gè)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小時(shí)
蟲號(hào): 1271842
注冊(cè): 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 的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓

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

[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +9	SUICHILD 2026-03-12	9/450	2026-03-19 09:35 by 花店相見
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4	1孫悟空 2026-03-17	4/200	2026-03-18 17:59 by fivewind
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 學(xué)碩285求調(diào)劑 +13	Wisjxn 2026-03-12	46/2300	2026-03-14 10:33 by JourneyLucky
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321

24小時(shí)熱門版塊排行榜

[求助] VOF體積分?jǐn)?shù)的梯度為什么為0呢？？

» 本帖附件資源列表

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

» 猜你喜歡

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

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

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

[求助] VOF體積分?jǐn)?shù)的梯度為什么為0呢？？