| 24小時(shí)熱門版塊排行榜 |
| 3 | 1/1 | 返回列表 |
| 查看: 1736 | 回復(fù): 2 | |||
| 【懸賞金幣】回答本帖問(wèn)題,作者無(wú)敵小張!將贈(zèng)送您 5 個(gè)金幣 | |||
無(wú)敵小張!新蟲(chóng) (初入文壇)
|
[求助]
fluent的udf
|
||
|
大神們,有段fluent的udf,但是看不懂,有沒(méi)有大神能給解答一下是什么意思 研一萌新沒(méi)啥金幣,直接分享,大家有需要就拿去用 #include "udf.h" #define U 18 #define I 200 #define PI 3.1415928 #define R 0.003 #define n 0.8 #define v 0.003 #define L 0.005 #define u0 1.26e-6 DEFINE_SOURCE(S_X,cell,thread,dS,eqn) { real xc[ND_ND]; real t; real source; real x,y,z; real r; C_CENTROID(xc,cell,thread); x=xc[0]; y=xc[1]; z=xc[2]; t=RP_Get_Real("flow-time" ;r=sqrt((x-v*t)*(x-v*t)+y*y); if(r<R) { source=-1*u0*I*I/(4*PI*PI*R*R*r)*exp(-1*r*r/(2*R*R)) *(1-exp(-1*r*r/(2*R*R)))*(1-z/L)*(1-z/L)*x/L; dS[eqn]=0; } else source=dS[eqn]=0; return source; } DEFINE_SOURCE(S_Y,cell,thread,dS,eqn) { real xc[ND_ND]; real t; real source; real x,y,z; real r; C_CENTROID(xc,cell,thread); x=xc[0]; y=xc[1]; z=xc[2]; t=RP_Get_Real("flow-time" ;r=sqrt((x-v*t)*(x-v*t)+y*y); if(r<R) { source=-1*u0*I*I/(4*PI*PI*R*R*r)*exp(-1*r*r/(2*R*R)) *(1-exp(-1*r*r/(2*R*R)))*(1-z/L)*(1-z/L)*y/L; dS[eqn]=0; } else source=dS[eqn]=0; return source; } DEFINE_SOURCE(S_Z,cell,thread,dS,eqn) { real xc[ND_ND]; real t; real source; real x,y,z; real r; C_CENTROID(xc,cell,thread); x=xc[0]; y=xc[1]; z=xc[2]; t=RP_Get_Real("flow-time" ;r=sqrt((x-v*t)*(x-v*t)+y*y); if(r<R) { source=-1*u0*I*I/(4*PI*PI*L*r*r) *(1-exp(-1*r*r/(2*R*R)))*(1-z/L)*(1-z/L); dS[eqn]=0; } else source=dS[eqn]=0; return source; } DEFINE_PROFILE(heat_flux,thread,index) { real xc[ND_ND]; face_t f; real t; real heat,T0=300; real x,y,z; real r; real hc = 80, sigmab = 5.67e-8, epsilon = 0.06; t=RP_Get_Real("flow-time" ;begin_f_loop(f,thread) { F_CENTROID(xc,f,thread); x=xc[0]; y=xc[1]; z=xc[2]; r=sqrt((x-v*t)*(x-v*t)+y*y); if(r<sqrt(6)*R) heat=n*I*U/(2*PI*R*R)*exp(-1*r*r/(2*R*R)); else heat=0; F_PROFILE(f, thread, index) = heat - hc*(F_T(f,thread), - T0) - sigmab*epsilon*(pow(F_T(f,thread), 4.) - pow(T0, 4.)); } end_f_loop(f,thread) } |
新蟲(chóng) (初入文壇)
新蟲(chóng) (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程085601調(diào)劑求老師收留 +8 | 強(qiáng)木木木 2026-03-07 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿河海大學(xué)085900土木水利專碩279求調(diào)劑不挑專業(yè) +3 | SunWwWwWw 2026-03-10 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +5 | 邱gl 2026-03-10 | 8/400 |
|
|
[考研] 工科0856專碩化學(xué)工程269能調(diào)劑嗎 +6 | 我想讀研11 2026-03-10 | 6/300 |
|
|
[考研] 調(diào)劑 +5 | 調(diào)劑的考研學(xué)生 2026-03-09 | 5/250 |
|
|
[考博] 26申博求助 +3 | 跳躍餅干 2026-03-10 | 4/200 |
|
|
[考研]
一志愿天津大學(xué)材料與化工275求調(diào)劑
10+5
|
穿只靴子 2026-03-07 | 19/950 |
|
|
[考研] 311求調(diào)劑 +3 | 牛乳糖的卡卡 2026-03-10 | 3/150 |
|
|
[考研] 材料工程專碩調(diào)劑自薦信,初試323 +6 | 李白26 2026-03-07 | 6/300 |
|
|
[考研] 328,0703考生求調(diào)劑,一志愿為東北師范大學(xué) +3 | 觀素律 2026-03-09 | 4/200 |
|
|
[考研] 334求調(diào)劑 +5 | Trying] 2026-03-07 | 5/250 |
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
[考研] 考研求調(diào)劑,工科,總分304 +4 | 矮子不想長(zhǎng)高了 2026-03-08 | 7/350 |
|
|
[考研] 一志愿211 085600 280數(shù)二英二求調(diào)劑 +3 | 月山斜 2026-03-06 | 3/150 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過(guò)來(lái) +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
[考研] 275求調(diào)劑 +4 | 大爆炸難民 2026-03-06 | 5/250 |
|
|
[考研] 材料085601一志愿哈工大317 +4 | 壓迫感行 2026-03-04 | 4/200 |
|
|
[考研] 274環(huán)境工程求調(diào)劑 +6 | 扶柳盈江 2026-03-05 | 6/300 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-04 | 3/150 |
|