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

返回列表

mydearkaren

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 181.2
散金: 16
帖子: 73
在線: 26.2小時(shí)
蟲號(hào): 2795305
注冊(cè): 2013-11-12
性別: MM
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

[求助] 離散相邊界條件UDF求助

參考幫助文檔寫了一個(gè)壁面碰撞的UDF，可以加載到FLUENT中，但是似乎完全沒有起到作用。
我想實(shí)現(xiàn)的功能是這樣的，當(dāng)顆粒速度大于給定值時(shí)，顆粒與壁面碰撞后彈開（可能有滑移，可能無滑移），否則顆粒被壁面捕獲（trap）。
加載UDF程序后，不論怎么調(diào)整給定值，顆�？偸菑楅_，從未被捕獲。甚至我把整個(gè)彈開部分的程序都刪掉，只留下Trap_Particle(p)結(jié)果還是一樣。
不知道哪里出了問題，請(qǐng)大家?guī)兔纯囱絶

程序如下：

/*reflect boundary condition for inert particles*/
#include "udf.h"
DEFINE_DPM_BC(bc_reflect,p,t,f,f_normal,dim)
{
real alpha;/*angle of particle path with face normal*/
real vn=0.,vt=0.,u=0.3,ud=0.15;
real nor_coeff=0.8,tan_coeff=0.3;
real normal[3];
real x[2];
real NV_VEX(x);
int i,idim=dim;

/*dim is always 2 in 2D compilation. Need special treatment for 2D axisymmetric and swirl flows*/
if (rp_axi_swirl)
{
real R=sqrt(P_POS(p)[1]*P_POS(p)[1]+P_POS(p)[2]*P_POS(p)[2]);
if (R>1.e-20)
{
/*二維法向量轉(zhuǎn)換為三維法向量*/
idim=3;
normal[0]=f_normal[0];
normal[1]=(f_normal[1]*P_POS(p)[1])/R;
normal[2]=(f_normal[1]*P_POS(p)[2])/R;
}
else
{
for (i=0;i<idim;i++)
normal=f_normal;
}
}
if(NV_MAG(P_VEL(p))>0.02)
{
Reflect_Particle(p,f_normal,dim,f,t,t,f);

if(p->type==DPM_TYPE_INERT)
{
alpha=M_PI/2-acos( MAX(-1,MIN(1,NV_DOT(normal,P_VEL(p))/MAX(NV_MAG(P_VEL(p)),DPM_SMALL))));
/*M_PI是一個(gè)宏定義，表示pi，3.14159265*/
/*acos反余弦函數(shù)*/
/*NV_DOT(x,u)在3D中等同于x[0]*u[0]+x[1]*u[1]+x[2]*u[2],2D中為x[0]*u[0]+x[1]*u[1]*/
/*NV_MAG（x）求模運(yùn)算。相當(dāng)于計(jì)算sqrt(x[0]*u[0]+x[1]*u[1]+x[2]*u[2])*/

if((NNULLP(t))&&(THREAD_TYPE(t)==THREAD_F_WALL))
/*NNULLP(t）是什么？*/
F_CENTROID(x,f,t);
/*F_CENTROID是找到面心坐標(biāo)，儲(chǔ)存在x中，2D中即face整個(gè)面的坐標(biāo)*/

/*calculate the normal component, rescale its magnitude by the coefficient of restitution and substract the change*/

/*compute normal velocity*/
for(i=0;i<idim;i++)
vn+=P_VEL(p)*normal;

/*substract off normal velocity.*/
for(i=0;i<idim;i++)
P_VEL(p)-=vn*normal;
vt=NV_MAG(P_VEL(p));
if (vt>7/2*u*(1+nor_coeff)*vn)
{
/*顆粒和壁面有滑移碰撞*/
for (i=0;i<idim;i++)
P_VEL(p)-=ud*(1+nor_coeff)*vn*P_VEL(p)/vt;
vn=-1*nor_coeff*vn;
}
/*顆粒和壁面無滑移碰撞*/
else
{
for (i=0;i<idim;i++)
P_VEL(p)*=5/7;
vn=-1*nor_coeff*vn;
}
for (i=0;i<idim;i++)
P_VEL(p)+=vn*normal;

/*store new velocity in P_VEL0 of particle*/
for(i=0;i<idim;i++)
P_VEL0(p)=P_VEL(p);

return PATH_ACTIVE;
}
}
else

Trap_Particle(p);

}

回復(fù)此樓

» 猜你喜歡

一志愿武理材料305分求調(diào)劑已經(jīng)有5人回復(fù)
材料專碩英一數(shù)二306 已經(jīng)有3人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有7人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線已經(jīng)有7人回復(fù)
295求調(diào)劑已經(jīng)有5人回復(fù)
354求調(diào)劑已經(jīng)有5人回復(fù)
材料專業(yè)求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有7人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有3人回復(fù)

1樓 2015-06-10 22:03:06

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

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

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

FLUENT仿真微通道多相流模型設(shè)置已經(jīng)有13人回復(fù)
求助：關(guān)于fluent多相流分析中液相顯示問題已經(jīng)有7人回復(fù)
這種次相體積分?jǐn)?shù)方程該用什么多相流模型？已經(jīng)有6人回復(fù)
有關(guān)fluent UDF的問題已經(jīng)有15人回復(fù)
Fluent源項(xiàng)是離散的數(shù)據(jù)文件已經(jīng)有10人回復(fù)
關(guān)于DPM模型模擬液滴的若干問題已經(jīng)有16人回復(fù)
fluent中離散相模擬時(shí)，非穩(wěn)態(tài)顆粒跟蹤怎樣輸出每個(gè)時(shí)間步的顆粒位置等信息？已經(jīng)有21人回復(fù)
求助：FLUENT進(jìn)出口流量不守恒的問題已經(jīng)有19人回復(fù)
【求助】fluent UDF設(shè)置邊界條件已經(jīng)有13人回復(fù)
fluent模擬噴霧過程相關(guān)問題求助已經(jīng)有34人回復(fù)
離散相融化問題已經(jīng)有8人回復(fù)
UDF的編譯通不過已經(jīng)有8人回復(fù)
DPM模型邊界設(shè)置問題已經(jīng)有9人回復(fù)
fluent模擬激光照射鋁板問題。求助已經(jīng)有6人回復(fù)
邊界條件UDF求助已經(jīng)有4人回復(fù)
運(yùn)用離散相模型模擬氣固兩相中如何添加顆粒的相關(guān)受力已經(jīng)有16人回復(fù)
求助氣固兩相流的邊界條件問題已經(jīng)有8人回復(fù)
【求助】氣液分離兩相流動(dòng)模擬fluent gambit出入口的設(shè)置已經(jīng)有11人回復(fù)

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

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

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