fluent里的用戶自定義函數(shù)(文件)
【正文】 UDS求解過程（續(xù)） ? 設(shè)定 Unsteady Function none (default), userdefined function （ DEFINE_UDS_UNSTEADY) UDS求解過程（續(xù)） ? 設(shè)定邊界條件 Defineboundary conditions UDS求解過程（續(xù)） ? （如果有源項）流體面板 defineboundary conditionsfluid UDS求解過程（續(xù)） ? 設(shè)定求解參數(shù)，給定初始值，計算 ? 檢查求解結(jié)果 scalarn Diffusion Coef. Of Scalarn 。在 ntx86下再建 2d、 2d_host、 2d_node等 ? (2d）然后在 2d文件夾下 ,把 \\fluent6. \src\和 \fluent6. \src \ ,將 名為 makefile ? 接著修改 ,要改成以下格式 : SOURCES = $(SRC)(源程序名字 ) VERSION = 2d(是二維還是三維 ) PARALLEL_NODE = none(有沒有并行 ) ? 最后點(diǎn)開始 程序 附件 命令提示符 ,進(jìn)入 DOS界面 ,訪問\libudf\ntx86\2d,然后敲 nmake,這樣就 OK了 . ? 然后在 FLUENT里 load這個 libudf就行了 . Compiled UDFs－ 方法二 ? 在 Compiled UDF里 ,點(diǎn) add source codes,找到源程序 ,加入 , ? 然后 build,再 LOAD就行了。 ? Compiled UDFs ： 編譯后的 UDF由 C語言系統(tǒng)的編譯器編譯成本地目標(biāo)碼。 C_UDMI(c,t,0) = (temptmin)/(tmaxtmin)。 tavg += temp*volume。 /* get cell volume */ temp = C_T(c,t)。 Thread *t。 /* declare domain pointer since it is not passed a */ /* argument to DEFINE macro */ { real tavg = 0.。 /* loop over all cell threads in the domain */ thread_loop_c (t,domain) { /* loop over all cells */ begin_c_loop_all (c,t) { C_CENTROID(xc,c,t)。它可以作為解釋程序或者編譯后的 UDF在FLUENT中執(zhí)行。 ? name是你所指定的 UDF的名字。 ? DEFINE_INIT 完成和修補(bǔ)一樣的功能，只是它以另一種方式 —— UDF來完成。區(qū)域變量提供了存取網(wǎng)格中所有單元和表面的線程。name是你所指定的 UDF的名字。 ? 可以用 DEFINE_ADJUST來修改流動變量（如：速度，壓力）并計算積分。 /* write out kount to data file */ } DEFINE_RW_FILE(reader, fp) { printf(Reading UDF data from data file...\n)。 /* define global variable kount */ DEFINE_ADJUST(demo_calc, domain) { kount++。 ? 這三個函數(shù)一起工作如下 。 ? 在 demo_calc中 ， kount的值每次迭代后增加因為每次迭代調(diào)用DEFINE_ADJUST一次 。 y = x[1]。 DEFINE_PROFILE(inlet_x_velocity, thread, index) { real x[ND_ND]。 循環(huán)宏用來在區(qū)域中每個面上循環(huán)以創(chuàng)建型線或數(shù)據(jù)數(shù)組 。 dS[eqn] = WEIGHT。 C_CENTROID(x, cell, thread)。 函數(shù)采用自變量 ds（ 它是數(shù)組的名字 ） 并設(shè)置由eqn指定的元素為關(guān)于速度 （ w_vel） 導(dǎo)數(shù)的值 。amp。指針 rr指向的目標(biāo)是設(shè)置反應(yīng)速率 。 else mu_lam = 1.。 DEFINE_PROPERTY(cell_viscosity, cell, thread) { real mu_lam。 DEFINE_PROFILE(inlet_x_velocity, thread, index) 定義在 define DEFINE_PROFILE(name, t, i) void name(Thread *t, int i) 在編譯過程中延伸為 void inlet_x_velocity(Thread *thread, int index) 名字為 inlet_x_velocity的函數(shù)不返回值由于它被聲明為空的數(shù)據(jù)類型 。 一旦 UDF被編譯 ， 它的名字 （ 例如 ，inlet_x_velocity） 將在 FLUENT適當(dāng)?shù)倪吔鐥l件面板 （ 例如 ， Velocity Inlet面板 ） 的下拉列表中變?yōu)榭梢姷暮涂蛇x的 。 ? 一旦函數(shù)被編譯 （ 和連接 ） ， 你為你的 UDF選擇的名字在 FLUENT下拉列表中將變成可見的和可選的 。 Domain數(shù)據(jù)類型代表了 FLUENT中最高水平的數(shù)據(jù)結(jié)構(gòu)。 Thread數(shù)據(jù)類型是 FLUENT中的數(shù)據(jù)結(jié)構(gòu)。 } FLUENT求解過程中 UDFs的先后順序 非耦合求解器 耦合求解器 FLUENT 網(wǎng)格拓?fù)? ? 單元（ cell） 區(qū)域被分割成的控制容積 ? 單元中心（ cell center） FLUENT中場數(shù)據(jù)存儲的地方 ? 面（ face） 單元（ 2D or 3D）的邊界 ? 邊（ edge） 面（ 3D）的邊界 ? 節(jié)點(diǎn)（ node） 網(wǎng)格點(diǎn) ? 單元線索（ cell thread） 在其中分配了材料數(shù)據(jù)和源項的單元組 ? 面線索（ face thread） 在其中分配了邊界數(shù)據(jù)的面組 ? 節(jié)點(diǎn)線索（ node thread） 節(jié)點(diǎn)組 ? 區(qū)域（ domain
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號-1