版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

yangshuaibyx

金蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1740.5
散金: 106
帖子: 167
在線: 645.8小時
蟲號: 2513219
注冊: 2013-06-19
性別: GG
專業(yè): 機械測試理論與技術

[求助] 曳力模型阻力系數(shù)的UDF

氣固兩相流數(shù)值模擬時fluent 中曳力模型的自定義UDF程序文件中有兩個模型，阻力系數(shù)調用函數(shù)時需要調用這兩個曳力模型嗎？還是只調用一個就行
/* This is a routine for customizing default Syamlal drag law in Fluent 6.
The default drag law uses 0.8 (for void<=0.85) and 2.65 (void>0.85) for
bfac. This is for a min fluid vel of 25 cm/s. The current drag law has been
tuned for a min fluid vel of 8 cm/s and uses 0.28 and 9.07 for these
parameters. */

#include "udf.h"
#include "sg_mphase.h"

# define pi 4.*atan(1.)
#define diam2 3.e-4

DEFINE_EXCHANGE_PROPERTY(custom_drag_syam, cell, mix_thread, s_col, f_col)
{
Thread *thread_g, *thread_s;
real x_vel_g, x_vel_s, y_vel_g, y_vel_s, abs_v, slip_x, slip_y,
   rho_g, rho_s, mu_g, reyp, afac,
   bfac, void_g, vfac, fdrgs, taup, k_g_s;

/* find the threads for the gas (primary) and solids (secondary phases).
These phases appear in columns 2 and 1 in the Interphase panel respectively*/

thread_g = THREAD_SUB_THREAD(mix_thread, s_col);/*gas phase*/
thread_s = THREAD_SUB_THREAD(mix_thread, f_col);/* solid phase*/

/* find phase velocities and properties*/

x_vel_g = C_U(cell, thread_g);
y_vel_g = C_V(cell, thread_g);

x_vel_s = C_U(cell, thread_s);
y_vel_s = C_V(cell, thread_s);

slip_x = x_vel_g - x_vel_s;
slip_y = y_vel_g - y_vel_s;

rho_g = C_R(cell, thread_g);
rho_s = C_R(cell, thread_s);

mu_g = C_MU_L(cell, thread_g);

/*compute slip*/
abs_v = sqrt(slip_x*slip_x + slip_y*slip_y);

/*compute reynolds number*/

reyp = rho_g*abs_v*diam2/mu_g;

/* compute particle relaxation time */

taup = rho_s*diam2*diam2/18./mu_g;

void_g = C_VOF(cell, thread_g);/* gas vol frac*/

/*compute drag and return drag coeff, k_g_s*/

afac = pow(void_g,4.14);

if(void_g<=0.85)
  bfac = 0.281632*pow(void_g, 1.28);
else
  bfac = pow(void_g, 9.076960);

vfac = 0.5*(afac-0.06*reyp+sqrt(0.0036*reyp*reyp+0.12*reyp*(2.*bfac-
            afac)+afac*afac));
fdrgs = void_g*(pow((0.63*sqrt(reyp)/vfac+4.8*sqrt(vfac)/vfac),2))/24.0;

k_g_s = (1.-void_g)*rho_s*fdrgs/taup;

return k_g_s;

}

DEFINE_EXCHANGE_PROPERTY(custom_drag_ihme, cell, mix_thread, s_col, f_col)

{

Thread *thread_g, *thread_s;
real x_vel_g, x_vel_s, y_vel_g, y_vel_s, abs_v, slip_x, slip_y,
   rho_g, rho_s, mu_g, reyp, cd, eg,
   void_g, k_g_s;

/* find the threads for the gas (primary) and solids (secondary phases).
These phases appear in columns 2 and 1 in the Interphase panel respectively*/

thread_g = THREAD_SUB_THREAD(mix_thread, s_col);/*gas phase*/
thread_s = THREAD_SUB_THREAD(mix_thread, f_col);/* solid phase*/

/* find phase velocities and properties*/

x_vel_g = C_U(cell, thread_g);
y_vel_g = C_V(cell, thread_g);

x_vel_s = C_U(cell, thread_s);
y_vel_s = C_V(cell, thread_s);

slip_x = x_vel_g - x_vel_s;
slip_y = y_vel_g - y_vel_s;

rho_g = C_R(cell, thread_g);
rho_s = C_R(cell, thread_s);

mu_g = C_MU_L(cell, thread_g);

/*compute slip*/
abs_v = sqrt(slip_x*slip_x + slip_y*slip_y);

/*compute reynolds number*/

reyp = rho_g*abs_v*diam2/mu_g;

cd = (24./(reyp+SMALL)) + 5.48*pow((reyp+SMALL),-0.573) + 0.36;

void_g = C_VOF(cell, thread_g);/* gas vol frac*/

eg = pow(void_g,-2.65);

k_g_s = (3./4.)*(cd*(1.-void_g)*abs_v*rho_g*eg)/diam2;

return k_g_s;

}

回復此樓

» 猜你喜歡

281求調劑（0805）已經(jīng)有16人回復
304求調劑已經(jīng)有6人回復
材料工程專碩調劑已經(jīng)有6人回復
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復
085700資源與環(huán)境308求調劑已經(jīng)有3人回復
求材料調劑已經(jīng)有8人回復
294求調劑材料與化工專碩已經(jīng)有5人回復
一志愿華中科技大學，080502，354分求調劑已經(jīng)有4人回復
一志愿吉林大學材料學碩321求調劑已經(jīng)有6人回復
085410人工智能專碩317求調劑（0854都可以）已經(jīng)有3人回復

1樓 2014-05-14 21:25:57

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

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關的精華帖子，希望有用哦~

關于DPM模型模擬液滴的若干問題已經(jīng)有16人回復
Fluent中DPM模型的UDF中使用P_POS()指針，然后給我提示undeclared variable 怎么辦? 已經(jīng)有6人回復
Fluent的UDF是否可以自定義輻射模型或湍流模型已經(jīng)有11人回復
兩相流歐拉模型模擬氣液分層流動，第二相diameter如何確定已經(jīng)有10人回復
請問液固流化床內(nèi)固體物料的初始流化速度與物料的初始堆積高度是否有關？已經(jīng)有22人回復
udf在非host節(jié)點下并行計算message 已經(jīng)有11人回復
fluent 中二維模型中，材料屬性隨著溫度變化的udf源代碼？已經(jīng)有14人回復
單機多核并行計算下UDF的問題已經(jīng)有9人回復
Fluent模擬流化床兩相顆粒之間的曳力系數(shù)選擇已經(jīng)有11人回復
fluent計算多相流問題時，UDF能否提取單元內(nèi)分相速度，以及單元平均速度已經(jīng)有5人回復
fluent 的混合模型中granular如何設置已經(jīng)有5人回復
歐拉雙流體模型中如何模擬出布風板？已經(jīng)有10人回復
UDF 改為并行計算已經(jīng)有9人回復
vof模型算氣液兩相逆流使用udf加動量源項碰到問題請高手解答已經(jīng)有14人回復
運用離散相模型模擬氣固兩相中如何添加顆粒的相關受力已經(jīng)有16人回復
曳力函數(shù)的的UDF 已經(jīng)有15人回復
流化床曳力系數(shù)UDF 已經(jīng)有9人回復
udf中的平方怎么寫已經(jīng)有5人回復
fluent中如何編寫風剖面的UDF 已經(jīng)有6人回復
udf 重力方向變化已經(jīng)有17人回復
【求助】請問fluent能模擬超重力的情況嗎？已經(jīng)有4人回復
【求助】UDF氣體的密度已經(jīng)有10人回復
【求助】fluent計算中，udf編譯通過，初始化出錯已經(jīng)有11人回復
【求助】如何改變Fluent中的曳力模型（UDF）已經(jīng)有12人回復

點擊這里搜索更多相關資源

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

相關版塊跳轉我要訂閱樓主 yangshuaibyx 的主題更新