版塊導(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í)熱門版塊排行榜

返回列表

youngxhhhh

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1.8
帖子: 36
在線: 40.1小時(shí)
蟲號(hào): 3077194
注冊(cè): 2014-03-21
專業(yè): 傳遞過程

[求助] 氣相在液相中的擴(kuò)散

想用fluent模擬氣相在液相中的擴(kuò)散，需要編寫關(guān)于質(zhì)量源項(xiàng)的udf，同時(shí)還需設(shè)定uds，但是沒有得到想要的結(jié)果，求大神們幫忙~不勝感激~
以下是自己編寫的程序，另外就是迭代過程中uds殘差一直為零，這樣肯定是有問題的，但不知道什地方出了問題，還請(qǐng)大神指點(diǎn)！
#include "udf.h"
real ci=1.0;
real D=1.76e-09;
real rhog=1.8183;
real cl,ug,ul,ur,kl,drb,rb,rb2,a,Vl,sl;

int ts;
DEFINE_SOURCE(liq, cell, thread_l, dS, eqn)
{
Thread *mix, *thread_g;
mix=THREAD_SUPER_THREAD(thread_l);
thread_g=THREAD_SUB_THREAD(mix, 1);
ts=RP_Get_Integer("time-step"

;
if(ts==1)
{
      sl=0.0;
dS[eqn]=0.0;
}
else
{
rb=0.002;
ug=C_V(cell, thread_g);
ul=C_V(cell, thread_l);
ur=ug-ul;
if(ur!=ug)
{
if(ur!=ul)
{
      cl=C_UDSI(cell, thread_l,0);
      kl=2*sqrt(D*ur/M_PI/rb/2.0);
      drb=kl*(ci-cl)/rhog;
      rb2=rb-drb;
      a=4.0*M_PI*pow(rb2,2.0);
      Vl=C_VOLUME(cell,mix);
      sl=kl*a*(ci-cl)/Vl;
      dS[eqn]=0.0;
      rb=rb2;
}
else
{
      sl=0.0;
      dS[eqn]=0.0;
}
}
else
{
sl=0.0;
dS[eqn]=0.0;
}
}
return sl;
}

DEFINE_SOURCE(gas, cell, thread_g, dS, eqn)
{
Thread *mix, *thread_l;
mix=THREAD_SUPER_THREAD(thread_g);
thread_l=THREAD_SUB_THREAD(mix, 0);
ts=RP_Get_Integer("time-step"

;
if(ts==1)
{
      sl=0.0;
dS[eqn]=0.0;
}
else
{
rb=0.002;
ug=C_V(cell, thread_g);
ul=C_V(cell, thread_l);
ur=ug-ul;
if(ur!=ug)
{
if(ur!=ul)
{
      cl=C_UDSI(cell, thread_l,0);
        kl=2*sqrt(D*ur/M_PI/rb/2.0);
      drb=kl*(ci-cl)/rhog;
      rb2=rb-drb;
      a=4.0*M_PI*pow(rb2,2.0);
        Vl=C_VOLUME(cell,mix);
      sl=-kl*a*(ci-cl)/Vl;
      dS[eqn]=-8.0*M_PI*kl*rb2*(ci-cl)/Vl;
      rb=rb2;
}
else
{
      sl=0.0;
      dS[eqn]=0.0;
}
}
else
{
sl=0.0;
dS[eqn]=0.0;
}
}
return sl;
}

DEFINE_RW_FILE(write, fp)
{
fprintf(fp, "%d", rb);
}

DEFINE_RW_FILE(reader, fp)
{
fscanf(fp, "%d", &rb);
}

回復(fù)此樓

» 猜你喜歡

284求調(diào)劑已經(jīng)有10人回復(fù)
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
07化學(xué)280分求調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有12人回復(fù)
求材料，環(huán)境專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
335求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
一志愿吉大化學(xué)322求調(diào)劑已經(jīng)有4人回復(fù)
環(huán)境學(xué)碩288求調(diào)劑已經(jīng)有8人回復(fù)
341求調(diào)劑(一志愿湖南大學(xué)070300) 已經(jīng)有6人回復(fù)

1樓 2014-06-01 21:56:32

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

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

請(qǐng)問下晶體的生長(zhǎng) 已經(jīng)有5人回復(fù)
高效液相和氣相的維護(hù)問題已經(jīng)有13人回復(fù)
目前常用的求取液相擴(kuò)散系數(shù)的方法已經(jīng)有3人回復(fù)
求教液相和氣相色譜的標(biāo)準(zhǔn)品的一般濃度已經(jīng)有7人回復(fù)
請(qǐng)問氣相，液相怎么拆？已經(jīng)有10人回復(fù)
是選擇氣相還是液相已經(jīng)有7人回復(fù)
請(qǐng)問氣相色譜和液相色譜都分別適合哪些物質(zhì)的檢測(cè)？已經(jīng)有6人回復(fù)
湖泊中氨根離子，硝酸根離子，以及磷酸根離子的擴(kuò)散系數(shù)是多少，謝謝各位！已經(jīng)有3人回復(fù)
求助：液相擴(kuò)散法培養(yǎng)單晶已經(jīng)有7人回復(fù)
【求助】求UDF高手幫忙看看我的問題已經(jīng)有18人回復(fù)
【求助】求助：學(xué)習(xí)氣相液相看那些資料已經(jīng)有18人回復(fù)
【求助】請(qǐng)教：液相與氣相檢測(cè)的區(qū)別？？已經(jīng)有4人回復(fù)
【討論】如何培養(yǎng)單晶呀已經(jīng)有3人回復(fù)
【交流】氣相液相出峰有什么規(guī)律已經(jīng)有4人回復(fù)
【求助】培養(yǎng)單晶的大蝦來幫幫忙啊已經(jīng)有21人回復(fù)
【資源】液相色譜常見問題及處理方法已經(jīng)有43人回復(fù)
氣液反應(yīng)速率常數(shù)可查么？如果不能怎么能得到已經(jīng)有7人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

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