opporoppor

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105.8
帖子: 20
在線: 6.2小時
蟲號: 3862070
注冊: 2015-05-11
專業(yè): 流體力學(xué)

[求助] 一個入口速度函數(shù)的UDF編程，新手實在無力，望大神不吝賜教！已有2人參與

入口速度是V=80*丨sin（（pi*x）/0.6）丨， “丨”表示絕對值，“pi”表示圓周率。
fluent新手，對宏命令知之甚少，希望大神能幫我編寫一個函數(shù)。深表感激！

回復(fù)此樓

» 猜你喜歡

086003食品工程求調(diào)劑已經(jīng)有3人回復(fù)
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑已經(jīng)有12人回復(fù)
一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑已經(jīng)有12人回復(fù)
材料調(diào)劑已經(jīng)有5人回復(fù)
081700 調(diào)劑 267分已經(jīng)有8人回復(fù)
請教下大家 2026年國家基金申請是雙盲審嗎？已經(jīng)有5人回復(fù)
276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷已經(jīng)有10人回復(fù)
一志愿南航材料專317分求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)學(xué)碩求調(diào)劑已經(jīng)有5人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

fluent udf請教已經(jīng)有16人回復(fù)

1樓 2015-05-18 14:36:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liukan12

銀蟲 (小有名氣)

應(yīng)助: 48 (小學(xué)生)
金幣: 499.6
紅花: 15
帖子: 198
在線: 63小時
蟲號: 607312
注冊: 2008-09-19
專業(yè): 分離過程

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

#include "udf.h"
DEFINE_PROFILE(velocity1,t,i)  /* 定義屬性宏三個參數(shù) 第一是名稱第二個是線程指針第三個是變量標(biā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)格體的中心點的坐標(biāo)值賦給數(shù)組x[ND_ND] */
   y = x[1];       /* 把x[ND_ND]中的第二個元素集合賦給變量y */
   if ( y > 0.5 )    /* 條件判斷，當(dāng)y值大于0.5時，執(zhí)行下列語句 */
   {
   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時，執(zhí)行語句 */
   }

   }
end_c_loop(f,t)       /* 遍歷循環(huán)結(jié)束 */
}
最后要解釋方法導(dǎo)入到fluent里再在速度入口條件下導(dǎo)入

贊一下

回復(fù)此樓

ANSYSFluent流體仿真研究人員。

2樓2015-05-18 17:59:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liukan12

銀蟲 (小有名氣)

應(yīng)助: 48 (小學(xué)生)
金幣: 499.6
紅花: 15
帖子: 198
在線: 63小時
蟲號: 607312
注冊: 2008-09-19
專業(yè): 分離過程

如果只要對面遍歷的話把上程序的C改成F Cell 改成Face

贊一下

回復(fù)此樓

ANSYSFluent流體仿真研究人員。

3樓2015-05-18 18:00:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

opporoppor

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105.8
帖子: 20
在線: 6.2小時
蟲號: 3862070
注冊: 2015-05-11
專業(yè): 流體力學(xué)

引用回帖:

3樓: Originally posted by liukan12 at 2015-05-18 18:00:18
如果只要對面遍歷的話把上程序的C改成F Cell 改成Face

老師您好，我這個是3D建模分析，您發(fā)的程序里面的解釋語句，我只能理解一小部分，比如“/* 條件判斷，當(dāng)y值大于0.5時，執(zhí)行下列語句 */”我就理解不了，另外一個，這個程序涉及絕對值和sin（）函數(shù)，我不會呀。懇請老師能辛苦一下，幫我編一個不需要改的函數(shù)，需要什么參數(shù)，您說。
我的金幣都可以給您！orz

贊一下

回復(fù)此樓

4樓2015-05-18 19:28:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liukan12

銀蟲 (小有名氣)

應(yīng)助: 48 (小學(xué)生)
金幣: 499.6
紅花: 15
帖子: 198
在線: 63小時
蟲號: 607312
注冊: 2008-09-19
專業(yè): 分離過程

【答案】應(yīng)助回帖

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

贊一下

回復(fù)此樓

ANSYSFluent流體仿真研究人員。

5樓2015-05-18 20:09:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

opporoppor

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105.8
帖子: 20
在線: 6.2小時
蟲號: 3862070
注冊: 2015-05-11
專業(yè): 流體力學(xué)

引用回帖:

5樓: Originally posted by liukan12 at 2015-05-18 20:09:53
fluent udf 是采用c++編寫的
如果是一個整數(shù)的話那么他的絕對是 abs（）長整型式 labs() 雙精度數(shù)是fabs()
正弦函數(shù)是 sin（）
所以你的可以是labs(sin(f(x))) f(x)就是你的表達式
自己多試試吧
y>0.5 就 ...

感謝老師的幫助，這個“real x[ND_ND]; /* 存放2D或者3D位置參數(shù) */”，我是3D模型，這個應(yīng)該怎么寫呢？我的模型很簡單，就是一個直管，請問這個“模型中y的坐標(biāo)值”怎么寫呢。。。老師，我其實懂的真的不多，我也不好意思一會兒問您一個問題，真的希望您能幫我寫一個詳細的orz，感激不盡！

贊一下

回復(fù)此樓

6樓2015-05-18 20:40:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liukan12

銀蟲 (小有名氣)

應(yīng)助: 48 (小學(xué)生)
金幣: 499.6
紅花: 15
帖子: 198
在線: 63小時
蟲號: 607312
注冊: 2008-09-19
專業(yè): 分離過程

ND 就是No dimension 就是說如果你是2維你申明的就是一個二維數(shù)組 x[0]放的是x坐標(biāo) x[1]是y坐標(biāo) 如果你是3維那么就是一個三維數(shù)組 x[0]放的是x坐標(biāo) x[1]是y坐標(biāo) x[2]是z坐標(biāo)

贊一下

回復(fù)此樓

ANSYSFluent流體仿真研究人員。

7樓2015-05-18 21:01:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

opporoppor

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105.8
帖子: 20
在線: 6.2小時
蟲號: 3862070
注冊: 2015-05-11
專業(yè): 流體力學(xué)

老師，幫幫忙吧，給我寫一個orz

[ 發(fā)自手機版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

8樓2015-05-18 23:42:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

moons2

金蟲 (正式寫手)

應(yīng)助: 144 (高中生)
金幣: 1440.5
紅花: 8
帖子: 431
在線: 239.4小時
蟲號: 3784607
注冊: 2015-04-03
性別: GG
專業(yè): 流體力學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
opporoppor: 金幣+4, ★★★很有幫助 2015-05-19 18:53:02

#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)
}

不保證對，自己多想多試。

贊一下

回復(fù)此樓

9樓2015-05-19 09:44:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

opporoppor

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 105.8
帖子: 20
在線: 6.2小時
蟲號: 3862070
注冊: 2015-05-11
專業(yè): 流體力學(xué)

引用回帖:

9樓: Originally posted by moons2 at 2015-05-19 09:44:46
#include "udf.h"
#include "math.h"
#define pi 3.14159
DEFINE_PROFILE(inlet_x_velocity, thread, position)
{
real x; /* this will hold the position vector */ ...

非常感謝老師！請問這個程序?qū)?yīng)的我的模型，是不是管道軸向必須是X向呢

[ 發(fā)自手機版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

10樓2015-05-19 10:37:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 求調(diào)劑 +7	十三加油 2026-03-21	7/350	2026-03-23 23:48 by 熱情沙漠
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +3	17501029541 2026-03-23	4/200	2026-03-23 23:47 by Txy@872106
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +5	niko- 2026-03-22	5/250	2026-03-23 22:01 by fuyu_
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	9/450	2026-03-23 12:57 by yqdszhdap－
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋懷
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11	楓橋ZL 2026-03-18	13/650	2026-03-22 20:26 by edmund7
[考研] 308求調(diào)劑 +3	墨墨漠 2026-03-21	3/150	2026-03-22 16:54 by i_cooler
[考研] 生物學(xué)071000 329分求調(diào)劑 +5	我愛生物生物愛?/a> 2026-03-17	5/250	2026-03-22 16:42 by tcx007
[考研] 尋找調(diào)劑 +4	倔強芒? 2026-03-21	4/200	2026-03-22 16:14 by 木托莫露露
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +9	niko- 2026-03-19	10/500	2026-03-22 16:08 by ColorlessPI
[考研] 305分求調(diào)劑（食品工程） +4	Sxy112 2026-03-21	6/300	2026-03-22 15:26 by 無懈可擊111
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無際的草原
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3	李軼男003 2026-03-20	3/150	2026-03-21 13:28 by 搏擊518
[考研] 307求調(diào)劑 +3	wyyyqx 2026-03-17	3/150	2026-03-21 03:20 by JourneyLucky
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO