博研聯(lián)盟_一個簡單而完整的UDS例子

1、標(biāo)準(zhǔn)方程形式:獨+_L，_rWL】=s巾（默認(rèn)爐;）ctj（cXjFluent中各項的定義如下:名稱表達(dá)式各項在UDS中需要定義的量在UDS中對應(yīng)的宏非穩(wěn)態(tài)項cP*ct離散形式的非穩(wěn)態(tài)項AtAt,.-3Papuiisu（詳細(xì)見UDF幫助）DEFINE_UDS_UNSTEADY對流項a通因通量聯(lián)DEFINE_UDS_FLUX擴散項生12cXj擴散率rDEFINE_DIFFUSIVITY假定需要求解方程為+2-=2日因I班）22=2,此方程解析解為巾=x+y+c,出I火；其中c為常數(shù) （假定其為0） ,可由邊界條件確定。 因此， 當(dāng)解穩(wěn)定時，*-x+y,等值線為一系列45度的斜直線。將此方程與標(biāo)準(zhǔn)

2、方程對比發(fā)現(xiàn):名稱表達(dá)式各項在UDS中需要定義的量需要在UDS中定義的宏非穩(wěn)態(tài)項那a離散形式的非穩(wěn)態(tài)項_也口+義巾nAt1AtDEFINE_UDS_UNSTEADY對流項因通量爐=I-A（單位向量I=1,1,1）DEFINE_UDS_FLUX擴散項C2*2CXj擴散率1=1DEFINE_DIFFUSIVITY源項22DEFINE_SOURCE邊界條件c=0旬=Xb+yb（b代表邊界處）DEFINE_PROFILE*注：源項和擴散率為常數(shù)，可以不用定義宏，但為了說明問題，仍采用宏來定義。一個簡單而完整的 UDS 例子gearboy當(dāng)tT8，即如果計算達(dá)到穩(wěn)定時, i； i ii將如下源碼，寫入一

3、個UDF文件，例如MyFistUDF.c,與cas文件放在同一目錄。#includeudf.hDEFINE_UDS_UNSTEADY(MyUnsteady,c,t,i,apu,su)realphysical_dt,vol,phi_old;physical_dt=RP_Get_Real(physical-time-step);vol=C_VOLUME(c,t);*apu=-vol/physical_dt;implicitpart*/phi_old=C_STORAGE_R(c,t,SV_UDSI_M1(i);*su=vol*phi_old/physical_dt;explicitpart*/DEF

4、INE_UDS_FLUX(MyFlux,f,t,i)realNV_VEC(unit_vec),NV_VEC(A);/聲明矢量變量F_AREA(A,f,t);NV_DS(unit_vec,=,1,1,1,*,1);單位矢量賦值returnNV_DOT(unit_vec,A);矢量點積DEFINE_DIFFUSIVITY(MyDiff,c,t,i)return1.0;DEFINE_SOURCE(MySource,c,t,dS,eqn)dSeqn=0;return2.0;假定如圖所示的計算域，其中正方形邊長為10,原點在正方形中心，四邊均為wall。DEFINE_PROFILE(MyProfile,

5、thread,index)(realxND_ND;/*thiswillholdthepositionvector*/realxx,yy;face_tf;begin_f_loop(f,thread)(F_CENTROID(x,f,thread);xx=x0;yy=x1;F_PROFILE(f,thread,index)=(xx+yy);4b=xb十yb)end_f_loop(f,thread)啟動非穩(wěn)態(tài)求解器，并在Define-User-Defined-Scalars面板中添加一個UDS,選擇自定義的通量和非穩(wěn)態(tài)項User-DefinedScalarNumberofUser-DefinedSca

6、larsOK|Cancel|Help|在材料面板中，選擇自定義擴散率FluxFunctionMyFlux：libudtUnsteadyFuinctionMyUnsteady:lib編譯通過并載入在boundaryconditions的Fluid面板中，選擇自定義源項MsmeMaterialTypeairfluid二lalerialsChemicalFormulaFluentFluidMaterialsairMixtureUser-DenonePropertiesDensity(kg/m3constantEdit1|Uier-DefinedFunctionsconstant1.225Viscosity(kg/m-s)1.7894e-05MyUnsteady;:libudfMyDiff;libudfUDSDiffusivity(kg/m-sMySaurce：libudfMyProfile：libudfjCanccljHelpChange/CreateDeleteCloseHelpFluei在boundaryconditions的壁面邊界條件面板中，選擇自定義邊界條件在solve-controls-solution中禁止Flow計算，只計算自定義標(biāo)量迭代過程中的殘差曲線，可以看出正逐漸