UDF-使用VOF計(jì)算泡沫的內(nèi)部循環(huán)

2017-01-19  by:CAE仿真在線  來(lái)源:互聯(lián)網(wǎng)



Problem:
It is important to get the internal circulation inside a bubble in many cases. Here the observer is sitting at the bubble itself.
1) Get the converged solution
2) The define on demand function below computes the velocity of the centroid of the bubble (and also the coordinates of the centroid).
3) The velocity vectors are modified depending upon this velocities.
Note: If you save this case file after executing the udf, then you will lose the actual velocity data. So make sure you have the original data securely saved.

 

Resolution:
**************************************************************************
* Udf to compute the velocity of the centroid of the bubble
* change the absolute velocity field to velocity field relative to
* the centroid of the bubble
* Written by Suman Basu of Fluent India
**************************************************************************/

#include "udf.h"

DEFINE_ON_DEMAND(move_with_bubble)
{
int phase_domain_index;
real vel_sum[ND_ND],co_ord1[ND_ND],co_ord[ND_ND],volume,average_velocity[ND_ND],x[ND_ND];
cell_t cell;
Thread *cell_thread;
Domain *subdomain,*mixture_domain=Get_Domain(1);
NV_S(vel_sum,=,0.);
NV_S(co_ord,=,0.);
volume=0.;
NV_S( average_velocity,=,0.);
NV_S(co_ord1,=,0.);
/* loop over all subdomains (phases) in the superdomain (mixture) */
sub_domain_loop(subdomain, mixture_domain, phase_domain_index)
{
/* loop if secondary phase */
if (DOMAIN_ID(subdomain) == 3)
/* loop over all cell threads in the secondary phase domain */
{
thread_loop_c (cell_thread,subdomain)
{
/* loop over all cells in secondary phase cell threads */
begin_c_loop(cell,cell_thread)
{
C_CENTROID(x,cell,cell_thread);
co_ord[0]+=C_VOF(cell,cell_thread)*C_VOLUME(cell,cell_thread)*x[0];
co_ord[1]+=C_VOF(cell,cell_thread)*C_VOLUME(cell,cell_thread)*x[1];
vel_sum[0]+= C_VOF(cell,cell_thread)*C_VOLUME(cell,cell_thread)*C_U(cell,cell_thread);
vel_sum[1]+= C_VOF(cell,cell_thread)*C_VOLUME(cell,cell_thread)*C_V(cell,cell_thread);
volume+=C_VOF(cell,cell_thread)*C_VOLUME(cell,cell_thread);
}
end_c_loop(cell,cell_thread)
}
}
}
NV_VS(co_ord1,=,co_ord,/,volume);
NV_VS(average_velocity,=,vel_sum,/,volume);
Message("Centroid[0]=%f\nCentroid[1]=%f\n",co_ord1[0],co_ord1[1]);
Message("V_Centroid[0]=%f\nV_Centroid[1]=%f\n",average_velocity[0],average_velocity[1]);

thread_loop_c (cell_thread,mixture_domain)
{
begin_c_loop(cell,cell_thread)
{
C_U(cell,cell_thread)=C_U(cell,cell_thread)-average_velocity[0];
C_V(cell,cell_thread)=C_V(cell,cell_thread)-average_velocity[1];
}
end_c_loop(cell,cell_thread)
}

}

開(kāi)放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:UDF-使用VOF計(jì)算泡沫的內(nèi)部循環(huán) Fluent培訓(xùn) Fluent流體培訓(xùn) Fluent軟件培訓(xùn) fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學(xué)反應(yīng) fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析 

編輯
在線報(bào)名:
  • 客服在線請(qǐng)直接聯(lián)系我們的客服,您也可以通過(guò)下面的方式進(jìn)行在線報(bào)名,我們會(huì)及時(shí)給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國(guó)服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號(hào)金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項(xiàng)目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com