FLUENT個人學(xué)習(xí)資料

用戶自定義函數(shù)UDF編寫1.邊界條件UDF1?1邊界條件與坐標(biāo)有關(guān)v=20-20X(-^-)2x 0.075進(jìn)口速度為拋物線的源代碼#include"udfh" /*C語言頭文件邊界條件，自動獲取DEFINE_PROFILE(inlet_x_velocity,thread,index){realx[ND_ND];realy;face_tf;/*variablestatement,positionvector數(shù)組*//*variablestatement,positionvalue*//*variablestatement,facethreadforloop*/begin_f_loop(f,thread)/*thread獲取Fluent中用戶定義的邊界線*/{F_CENTROID(x,f,thread); /*獲取面質(zhì)心,視維度將坐標(biāo)返回給x*/y=x[1]; /*x[1]表示y方向坐標(biāo)值*/F_PROFILE(f,thread,index)=20.-y*y/(.0745*.0745)*20.;廠end_f_loop(f,thread)/*thread建立的線索引進(jìn)行face遍歷*/}一1.2邊界條件與溫度有關(guān)T=T=300+100Xsin(兀J)X 0.005邊界溫度為正弦變化的函數(shù)#include"udf.h"/*C語言頭文件*/#definePI3.141592654 .十…丁邊界條件，自動獲取DEFINE_PROFILE(inlet_x_velocity,thread,index)/*variablestatement,positionvector數(shù)組/*variablestatement,positionvector數(shù)組*//*variablestatement,positionvalue*//*variablestatement,facethreadforloop*/realr[ND_ND];realx;face_tf;begin_f_loop(f,thread)/*thread獲取Fluent中用戶定義的邊界線*/{—F_CENTROID(r,f,thread); /*獲取面質(zhì)心,視維度將坐標(biāo)返回給x*/x=r[0]; /*r[0]表示x方向坐標(biāo)值*/F_PROFILE(f,thread,index)=300.+100.*sin(PI*x/0.005);廠end_f_loop(f,thread)/*thread建立的線索引進(jìn)行face遍歷*/}一1.3邊界條件與時間有關(guān)#include"udf.h" /*C語言頭文件*/#definePI3.141592654DEFINE_PROFILE(unsteady,thread,index){——face_tf; /*variablestatement,facethreadforloop*/begin_f_loop(f,thread)/*thread獲取Fluent中用戶定義的邊界線*/{—realt=RP_Get_Real(“flow-time”)；/*獲取計算過程中當(dāng)前時間*/F_PROFILE(f,thread,index)=20.+5.0*sin(10.*t);/*修改邊界條件*/廠end_f_loop(f,thread) /*thread建立的線索引進(jìn)行face遍歷*/}—2.內(nèi)部物性UDF2?1流體粘度與溫度有關(guān)#include"udf.h"DEFINE_PROPERTY(cell_viscosity,cell,thread)/*修改物性*/{ /*用戶不選擇邊界，函數(shù)自動遍歷整個計算區(qū)域*/realmu_lam; /*變量聲明，存儲粘度值*/realtemp=C_T(cell,thread); /*對計算區(qū)域cell遍歷，取溫度*/if(temp>288.)mu_lam=5.5e-3;elseif(temp>286.)mu_lam=143.2135-0.49725*temp;elsemu_lam=1.;returnmu_lam; /*返回給物性參數(shù)，流體粘度選項*/} — /*物性參數(shù)修改為返回型函數(shù)*/3?源項編寫UDF3.1多孔介質(zhì)區(qū)域的x維度源項方程式source=-0.5xC2py\v\v（其中C2是常數(shù)，取100）2 XX 2源項source中p可時時獲取，對定區(qū)域內(nèi)是常量。同理，質(zhì)心y是常量。此例中，源項的變量為耳，求解出dS/dq=dS[eqn]。XX???source=-A\v\v=dS/du=-2^\v\XXXX源項UDF必須給出源項對所求變量的直接導(dǎo)數(shù)，源項的導(dǎo)數(shù)用于對源項的線性化。下面是根據(jù)上述方程編寫的UDF的C源程序。

include"udf.h"defineC2100.0DEFINE_SOURCE(xmom_source,cell,thread,dS,eqn){ /*全局cell遍歷，x方向源項添加*/realx[ND_ND]; /*存儲質(zhì)心坐標(biāo)，源項中y值*/realcon,source;C_CENTROID(x,cell,thread); /*獲取質(zhì)心坐標(biāo)*/con=C2*0.5*C_R(cell,thread)*x[1]; /*得源項中的密度或A值*/source=-con*fabs(C_U(cell,thread))*C_U(cell,thread);/*源項式*/dS[eqn]=-2.*con*fabs(C_U(cell,thread)); /*源項倒數(shù)*/returnsource; /*返回源項*/} /*CU、CV、CW三維度下速度*/3.2管內(nèi)流動源項的編寫溫度為290K液態(tài)金屬進(jìn)入管道，管道上壁280K對液態(tài)金屬冷卻。當(dāng)液態(tài)金屬溫度低于288K時，就添加動量方程的源項。由于添加了源項，液態(tài)金屬速度會逐漸降為零。在本例中，僅以液態(tài)金屬的速度大小來表征是否凝固，如果速度等于零，則認(rèn)為液態(tài)金屬已經(jīng)凝固，而不是求解能量方程。x方向動量源項方程為：(式中，C=20)S=S=—Cv=込=—C=dS[eqn]dvx#include"udf.h"#defineCON20.0DEFINE_SOURCE(cell_x_source,cell,thread,dS,eqn){"""realsource;if(C_T(cell