一個簡單而完整的UDS例子

1、一個簡單而完整的UDS例子 gearboy標準方程形式： (默認)Fluent中各項的定義如下：名稱表達式各項在UDS中需要定義的量在UDS中對應的宏非穩(wěn)態(tài)項離散形式的非穩(wěn)態(tài)項(詳細見UDF幫助)DEFINE_UDS_UNSTEADY對流項通量 DEFINE_UDS_FLUX擴散項擴散率 DEFINE_DIFFUSIVITY假定需要求解方程為當，即如果計算達到穩(wěn)定時，此方程解析解為，其中c為常數(shù)（假定其為0），可由邊界條件確定。因此，當解穩(wěn)定時，等值線為一系列45度的斜直線。將此方程與標準方程對比發(fā)現(xiàn)：名稱表達式各項在UDS中需要定義的量需要在UDS中定義的宏非穩(wěn)態(tài)項離散形式的非穩(wěn)態(tài)項DEFI

2、NE_UDS_UNSTEADY對流項通量（單位向量）DEFINE_UDS_FLUX擴散項擴散率 DEFINE_DIFFUSIVITY源項22DEFINE_SOURCE邊界條件 (b代表邊界處)DEFINE_PROFILE*注：源項和擴散率為常數(shù)，可以不用定義宏，但為了說明問題，仍采用宏來定義。假定如圖所示的計算域，其中正方形邊長為10，原點在正方形中心，四邊均為wall。將如下源碼，寫入一個UDF文件，例如MyFistUDF.c，與cas文件放在同一目錄。#include "udf.h"DEFINE_UDS_UNSTEADY(MyUnsteady,c,t,i,apu,su)

3、real physical_dt, vol, phi_old;physical_dt = RP_Get_Real("physical-time-step");vol = C_VOLUME(c,t);*apu = -vol / physical_dt; /*implicit part*/phi_old = C_STORAGE_R(c,t,SV_UDSI_M1(i);*su = vol*phi_old/physical_dt; /*explicit part*/DEFINE_UDS_FLUX(MyFlux,f,t,i)real NV_VEC(unit_vec), NV_VEC(

4、A); /聲明矢量變量F_AREA(A, f, t);NV_DS(unit_vec, =, 1, 1, 1, *, 1); /單位矢量賦值return NV_DOT(unit_vec, A); /矢量點積DEFINE_DIFFUSIVITY(MyDiff,c,t,i)return 1.0;DEFINE_SOURCE(MySource,c,t,dS,eqn) dSeqn=0; return 2.0;DEFINE_PROFILE(MyProfile,thread,index)real xND_ND; /* this will hold the position vector */real xx,y

5、y;face_t f;begin_f_loop(f,thread) F_CENTROID(x,f,thread);xx = x0;yy=x1;F_PROFILE(f,thread,index) = (xx+yy); /end_f_loop(f,thread)編譯通過并載入啟動非穩(wěn)態(tài)求解器，并在Define-User-Defined-Scalars面板中添加一個UDS，選擇自定義的通量和非穩(wěn)態(tài)項在材料面板中，選擇自定義擴散率在boundary conditions的Fluid面板中，選擇自定義源項在boundary conditions的壁面邊界條件面板中，選擇自定義邊界條件在solve-controls-solution中禁止Flow計算，只計算自定義標量迭代過程中的殘差曲線，可以看出正逐漸趨向穩(wěn)定采用后處理等值線面板，觀察自定義標量的等值線