| 查看: 2410 | 回復(fù): 21 | |||
opporoppor新蟲 (初入文壇)
|
[求助]
一個(gè)入口速度函數(shù)的UDF編程,新手實(shí)在無(wú)力,望大神不吝賜教! 已有2人參與
|
|
入口速度是V=80*丨sin((pi*x)/0.6)丨, “丨”表示絕對(duì)值,“pi”表示圓周率。 fluent新手,對(duì)宏命令知之甚少,希望大神能幫我編寫一個(gè)函數(shù)。深表感激! |
|
#include "udf.h" DEFINE_PROFILE(velocity1,t,i) /* 定義屬性宏 三個(gè)參數(shù) 第一是名稱 第二個(gè)是線程指針 第三個(gè)是變量標(biāo)號(hào) */ { real x[ND_ND]; /* 存放2D或者3D位置參數(shù) */ real y; /* y值為模型中的y坐標(biāo)值 */ cell_t f; /* 創(chuàng)建網(wǎng)格體指針 */ begin_c_loop(f,t) /* 體線程t中所有網(wǎng)格體指針f的遍歷循環(huán) */ { C_CENTROID(x,f,t); /* 把體線程t中所有指針f中的網(wǎng)格體的中心點(diǎn)的坐標(biāo)值賦給數(shù)組x[ND_ND] */ y = x[1]; /* 把x[ND_ND]中的第二個(gè)元素集合賦給變量y */ if ( y > 0.5 ) /* 條件判斷,當(dāng)y值大于0.5時(shí),執(zhí)行下列語(yǔ)句 */ { C_PROFILE(f, t, i) = y^2; /* y^2是表示入口速度等于y坐標(biāo)值的平方,這里可以改成你要的關(guān)系式 */ } else /* 條件判斷 */ { C_PROFILE(f, t, i) = 1; /* 條件判斷,當(dāng)y值小于0.5時(shí),執(zhí)行語(yǔ)句 */ } } end_c_loop(f,t) /* 遍歷循環(huán)結(jié)束 */ } 最后要解釋方法導(dǎo)入到fluent里 再在速度入口條件下導(dǎo)入 |


新蟲 (初入文壇)
|
fluent udf 是采用c++編寫的 如果是一個(gè)整數(shù)的話 那么他的絕對(duì)是 abs() 長(zhǎng)整型式 labs() 雙精度數(shù)是fabs() 正弦函數(shù)是 sin() 所以你的可以是labs(sin(f(x))) f(x)就是你的表達(dá)式 自己多試試吧 y>0.5 就是當(dāng)y軸坐標(biāo)大于0.5時(shí) 讓你的速度等于那個(gè)等式 y<0.5 執(zhí)行另外一個(gè)表達(dá)式 就是分段函數(shù) 你這里是單一區(qū)間 不用這個(gè)條件判斷 直接寫就行了 |

新蟲 (初入文壇)

新蟲 (初入文壇)
金蟲 (正式寫手)
|
#include "udf.h" #include "math.h" #define pi 3.14159 DEFINE_PROFILE(inlet_x_velocity, thread, position) { real x[ND_ND]; /* this will hold the position vector */ real y; face_t f; begin_f_loop(f, thread) { F_CENTROID(x,f,thread); y = x[0]; F_PROFILE(f, thread, position) =80.*abs(sin(pi*y/.6)); } end_f_loop(f, thread) } 不保證對(duì),自己多想多試。 |
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 340求調(diào)劑 +4 | 話梅糖111 2026-03-24 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +7 | 天天好運(yùn)來(lái)上岸?/a> 2026-03-24 | 7/350 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +7 | holy J 2026-03-21 | 7/350 |
|
|
[考研]
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
[考研] 299求調(diào)劑 +7 | 某某某某位 2026-03-21 | 7/350 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +7 | vv迷 2026-03-22 | 7/350 |
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +6 |  ̄^ ̄゜汗 2026-03-19 | 9/450 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 9/450 |
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽(yáng)CF 2026-03-18 | 5/250 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|