版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

steven_tan

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 153.5
帖子: 48
在线: 16.7小时
虫号: 3137296
注册: 2014-04-15
专业: 工程管理

[求助] 模拟波浪已有1人参与

我现在在用边界造波法模拟波浪，但一直没做出正确的模拟，下面是我的UDF文件，麻烦大神看看哪里不正确，我的模型长L0=100m,水深H=1.5m,水面以上H0-H=1.5m,坐标原点在自由界面处
#include "udf.h"
real AA=0.3; /*波高*/
real TT=2.0;  /*周期*/
real LL=6.2;  /*波长*/
real FL=100;  /*水槽长度*/

real kk=2.0*3.1415925/6.2;
real ww=2.0*3.1415925/2.0;
real h=1.5; /*水深*/
// x direction  velocity  //
DEFINE_PROFILE(x_velocity,thread,index)
{
  real x[ND_ND];
  real xx;
  real y;
  face_t f;
  real u = 0;
  real  t = RP_Get_Real("flow-time"

;
  begin_f_loop(f,thread)
{
   F_CENTROID(x,f,thread);

   xx = x[0];

   y = x[1];
   if(y<(AA*cos(kk*xx-ww*t)))
   u=3.14*2.0*AA*cosh(kk*y)*cos(kk*xx-ww*t)/TT/sin(kk*h);
   else u = 0.0;
F_PROFILE(f,thread,index) = u;
}

  end_f_loop(f,thread)
}
// y direction  velocity  //
DEFINE_PROFILE(y_velocity,thread,index)
{
  real x[ND_ND];
  real xx;
  real y;
  face_t f;
real  v = 0;
  real  t = RP_Get_Real("flow-time"

;

  begin_f_loop(f,thread)
{
   F_CENTROID(x,f,thread);
   xx = x[0];

   y = x[1];
if(y<(AA*cos(kk*xx-ww*t)))

   v=3.14*2.0*AA*sinh(kk*y)*sin(kk*xx-ww*t)/TT/sinh(kk*h);
else v = 0.0;
F_PROFILE(f,thread,index) = v;

}

  end_f_loop(f,thread)
}

//入口体积分数//
DEFINE_PROFILE(voffactor,thread,index)

{
  real x[ND_ND];
  real xx;
  real y;
  face_t f;
  real  t = RP_Get_Real("flow-time"

;

  begin_f_loop(f,thread)
{
   F_CENTROID(x,f,thread);
   xx=x[0];
   y = x[1];

   if(y<(AA*cos(kk*xx-ww*t)))

   F_PROFILE(f,thread,index) = 1.0;

   else F_PROFILE(f,thread,index) = 0.0;

}

  end_f_loop(f,thread)
}

// 消波x direction  momentum  //
DEFINE_SOURCE(x_mom_source,c,t,dS,eqn)
{
real x[ND_ND];
real x_source;
real y;
real chan;
real xishu;
real time=RP_Get_Real("flow-time"

;
y=x[1];
C_CENTROID(x,c,t);
if(x[0]>=FL-LL&&x[0]<=FL)
{
if(y>=0-AA&&y<=0)
{
xishu=8*(x[0]-(FL-LL))/LL;
x_source=-C_U(c,t)*998.2*xishu;
dS[eqn]=-998.2*xishu;
}
}
return x_source;
}
// y direction momentum //
DEFINE_SOURCE(y_mom_source,c,t,dS,eqn)
{
real x[ND_ND];
real y_source;
real y;
real chan;
real xishu;
real time=RP_Get_Real("flow-time"

;
y=x[1];
C_CENTROID(x,c,t);
if(x[0]>=FL-LL&&x[0]<=FL)
{
if(y>=0-AA&&y<=0)
{
xishu=8*(x[0]-(FL-LL))/LL;
y_source=-C_V(c,t)*998.2*xishu;
dS[eqn]=-998.2*xishu;
}
}
return y_source;
}
DEFINE_PROFILE(out_pressure,thread,position)
{
real x[ND_ND];
real y;
face_t f;
begin_f_loop(f,thread)
{
F_CENTROID(x,f,thread);
y=x[1];
if(y<=0)
F_PROFILE(f,thread,position)=998.2*9.81*(-y);
else
F_PROFILE(f,thread,position)=0;
}
end_f_loop(f,thread)
}
//出口体积分数//
DEFINE_PROFILE(outlet_voffactorc,thread,position)
{
real x[ND_ND];
real y;
face_t f;
begin_f_loop(f,thread)
{
F_CENTROID(x,f,thread);
y=x[1];
if(y<=0)
F_PROFILE(f,thread,position)=1;
else
F_PROFILE(f,thread,position)=0;
}
end_f_loop(f,thread)
}

模型.png

回复此楼

» 猜你喜欢

一志愿中科大材料与化工，353分还有调剂学校吗已经有13人回复
085600材料与化工调剂已经有6人回复
320求调剂已经有6人回复
期刊推荐已经有5人回复
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业已经有7人回复
又一批高校组建人工智能学院师资行吗不是骗人吗已经有4人回复
297，工科调剂?河南农业大学本科已经有15人回复
337求调剂已经有3人回复
申博/考博已经有5人回复
申博已经有3人回复

1楼 2014-10-19 23:12:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

死天

专家顾问 (著名写手)

专家经验: +60
应助: 399 (硕士)
金币: 3648.9
红花: 89
帖子: 1336
在线: 443.8小时
虫号: 2643057
注册: 2013-09-09
专业: 流体力学
管辖: 仿真模拟

【答案】应助回帖

感谢参与，应助指数 +1

源项里面 if条件后的其他条件呢？？？

赞一下

回复此楼

2楼2014-10-20 08:08:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

steven_tan

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 153.5
帖子: 48
在线: 16.7小时
虫号: 3137296
注册: 2014-04-15
专业: 工程管理

引用回帖:

2楼: Originally posted by 死天 at 2014-10-20 08:08:26
源项里面 if条件后的其他条件呢？？？

if(x[0]>=FL-LL&&x[0]<=FL)
{
if(y>=0-AA&&y<=0)
{
xishu=8*(x[0]-(FL-LL))/LL;
y_source=-C_V(c,t)*998.2*xishu;
dS[eqn]=-998.2*xishu;
}
}
还有其他条件吗？？编写源项我只知道这么多，您能告诉我还需要什么吗，拜谢

赞一下

回复此楼

3楼2014-10-20 09:06:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

死天

专家顾问 (著名写手)

专家经验: +60
应助: 399 (硕士)
金币: 3648.9
红花: 89
帖子: 1336
在线: 443.8小时
虫号: 2643057
注册: 2013-09-09
专业: 流体力学
管辖: 仿真模拟

引用回帖:

3楼: Originally posted by steven_tan at 2014-10-20 09:06:27
if(x>=FL-LL&&x<=FL)
{
if(y>=0-AA&&y<=0)
{
xishu=8*(x-(FL-LL))/LL;
y_source=-C_V(c,t)*998.2*xishu;
dS=-998.2*xishu;
}
}
还有其他条件吗？？编写源项我只知道 ...

那ELSE呢？IF只是满足这个条件的时候，那么其他条件呢？至少先个C的条件语句的书看看啊

赞一下

回复此楼

4楼2014-10-20 09:14:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

steven_tan

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 153.5
帖子: 48
在线: 16.7小时
虫号: 3137296
注册: 2014-04-15
专业: 工程管理

引用回帖:

4楼: Originally posted by 死天 at 2014-10-20 09:14:33
那ELSE呢？IF只是满足这个条件的时候，那么其他条件呢？至少先个C的条件语句的书看看啊...

if没有必要一定要跟else的吧，1、if{}；2、if{}……else；3、if{}……else if；这些都可以的
我写这个UDF的意思就是满足if条件，就调用程序；不满足，就不调用程序

赞一下

回复此楼

5楼2014-10-20 09:47:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

steven_tan

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 153.5
帖子: 48
在线: 16.7小时
虫号: 3137296
注册: 2014-04-15
专业: 工程管理

引用回帖:

4楼: Originally posted by 死天 at 2014-10-20 09:14:33
那ELSE呢？IF只是满足这个条件的时候，那么其他条件呢？至少先个C的条件语句的书看看啊...

您那有没有可以用的源项文件，我这个虽然编译没有问题，但是感觉在FLUENT里没有起到应有的效果

赞一下

回复此楼

6楼2014-10-20 09:56:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

死天

专家顾问 (著名写手)

专家经验: +60
应助: 399 (硕士)
金币: 3648.9
红花: 89
帖子: 1336
在线: 443.8小时
虫号: 2643057
注册: 2013-09-09
专业: 流体力学
管辖: 仿真模拟

引用回帖:

5楼: Originally posted by steven_tan at 2014-10-20 09:47:02
if没有必要一定要跟else的吧，1、if{}；2、if{}……else；3、if{}……else if；这些都可以的
我写这个UDF的意思就是满足if条件，就调用程序；不满足，就不调用程序...

那不满足条件的返回值是什么？源项有2个返回值，source 和ds[eqn].少了一半

赞一下

回复此楼

7楼2014-10-20 10:02:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

死天

专家顾问 (著名写手)

专家经验: +60
应助: 399 (硕士)
金币: 3648.9
红花: 89
帖子: 1336
在线: 443.8小时
虫号: 2643057
注册: 2013-09-09
专业: 流体力学
管辖: 仿真模拟

【答案】应助回帖

引用回帖:

6楼: Originally posted by steven_tan at 2014-10-20 09:56:49
您那有没有可以用的源项文件，我这个虽然编译没有问题，但是感觉在FLUENT里没有起到应有的效果

...

看我楼上的发言，你自己的实验，你自己去分析各个条件。实验不一样，UDF能一样吗?

赞一下

回复此楼

8楼2014-10-20 10:03:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

steven_tan

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 153.5
帖子: 48
在线: 16.7小时
虫号: 3137296
注册: 2014-04-15
专业: 工程管理

引用回帖:

8楼: Originally posted by 死天 at 2014-10-20 10:03:24
看我楼上的发言，你自己的实验，你自己去分析各个条件。实验不一样，UDF能一样吗?...

我不是专攻这个方向的，单纯就是做模拟，希望能做出来波形，我现在的情况是波形做不出，周围没人做这个，自己捣鼓很长时间了，也没解决问题，所以希望您能给一个正确的UDF文件，学习一下，拜托啦

赞一下

回复此楼

9楼2014-10-20 10:26:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

死天

专家顾问 (著名写手)

专家经验: +60
应助: 399 (硕士)
金币: 3648.9
红花: 89
帖子: 1336
在线: 443.8小时
虫号: 2643057
注册: 2013-09-09
专业: 流体力学
管辖: 仿真模拟

引用回帖:

9楼: Originally posted by steven_tan at 2014-10-20 10:26:39
我不是专攻这个方向的，单纯就是做模拟，希望能做出来波形，我现在的情况是波形做不出，周围没人做这个，自己捣鼓很长时间了，也没解决问题，所以希望您能给一个正确的UDF文件，学习一下，拜托啦

...

我也不是搞这个。我怎么会知道

赞一下

回复此楼

10楼2014-10-20 11:01:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 steven_tan 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +11	否极泰来2026 2026-04-15	13/650	2026-04-20 22:31 by Equinoxhua
[考研] 085600材料与化工调剂 5+3	孜孜不倦2002 2026-04-19	6/300	2026-04-20 21:25 by babero
[考研] 320求调剂 +6	深郊akm 2026-04-17	6/300	2026-04-20 18:57 by fs26jie
[论文投稿] 期刊推荐 +3	材料研究生 2026-04-15	5/250	2026-04-20 16:02 by 豆豆7758
[考研] 337求调剂 +3	jyz04 2026-04-18	3/150	2026-04-20 12:24 by 研可安
[考博] 申博/考博 +4	啃面包的小书虫 2026-04-17	5/250	2026-04-20 10:47 by YuY66
[考博] 申博 +3	Xyyx. 2026-04-18	3/150	2026-04-20 10:44 by YuY66
[考研] 通信工程求调剂！！！ +7	zlb770521 2026-04-14	7/350	2026-04-19 20:56 by Equinoxhua
[考研] 297，工科调剂? +11	河南农业大学-能 2026-04-14	11/550	2026-04-19 20:07 by Equinoxhua
[考研] 291求调剂 +11	关忆北. 2026-04-14	11/550	2026-04-19 17:16 by 中豫男
[考研] 304求调剂 +8	castLight 2026-04-16	8/400	2026-04-19 17:14 by 中豫男
[考研] 求调剂推荐 +9	小聂爱学习 2026-04-14	9/450	2026-04-19 17:03 by 中豫男
[考研] 297，工科调剂? +5	河南农业大学-能 2026-04-14	5/250	2026-04-18 15:17 by Equinoxhua
[考研] 22408 312求调剂 +24	门路摸摸 2026-04-14	26/1300	2026-04-18 13:04 by wunaiy88
[考研] 收到复试调剂但是去不了 +8	小蜗牛* 2026-04-16	8/400	2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4	Zyt1314520.. 2026-04-17	5/250	2026-04-18 08:28 by babysonlkd
[考研] 322求调剂 +6	tekuzu 2026-04-17	6/300	2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9，生物学326求调剂 +9	刘墨墨 2026-04-15	9/450	2026-04-16 17:14 by 崔崔崔cccc
[基金申请] RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一 +7	zju2000 2026-04-14	18/900	2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 一志愿A区211，22408 321求调剂 +6	随心所欲☆ 2026-04-15	7/350	2026-04-15 21:45 by lbsjt

24小时热门版块排行榜

[求助] 模拟波浪 已有1人参与

» 猜你喜欢

【答案】应助回帖

【答案】应助回帖

[求助] 模拟波浪已有1人参与