24小时热门版块排行榜    

查看: 2063  |  回复: 1

summer_cai

新虫 (初入文坛)

[求助] 【已完结】关于fluent模拟波浪的问题

本人是fluent新手,想利用udf模拟波浪,下面是我编的源项程序,在运行的时候出现overflow的错误,运行不下去,想请教各位大神看看是不是我编的程序错了,拜托了
#include<math.h>
#include "udf.h"
#define k 2*3.14159265/8.0
#define H 0.8
#define w 2.617994
#define h 4.0
#define T 2*3.14159265/2.617994
#define deltax 0.1
DEFINE_SOURCE(zbc_source,c,t,dS,eqn)
{
real x[ND_ND];
real q,source;
real y;
real dt;
dt=RP_Get_Real("flow-time" );
C_CENTROID(x,c,t);
y=x[1];
q=C_R(c,t)*H*w*(cosh(k*y+k*h)*cos(w*dt)/sinh(k*h)+3/8*H*k*cosh(2*k*y+2*k*h)*cos(2*w*dt)/(sinh(k*h)*sinh(k*h)*sinh(k*h)*sinh(k*h)))/deltax;
source=q;
dS[eqn]=0;
return source;
}
DEFINE_SOURCE(zbx_source,c,t,dS,eqn)
{
real x[ND_ND];
real q,source;
real y;
real dt;
dt=RP_Get_Real("flow-time" ) ;
C_CENTROID(x,c,t) ;
y=x[1];
q=C_R(c,t)*H*w*(cosh(k*y+k*h)*cos(w*dt)/sinh(k*h)+3/8*H*k*cosh(2*k*y+2*k*h)*cos(2*w*dt)/(sinh(k*h)*sinh(k*h)*sinh(k*h)*sinh(k*h)))/deltax;
source=C_R(c,t)*C_U(c,t)*q;
dS[eqn]=C_R(c,t)*q;
return source;
}
DEFINE_SOURCE(zby_source,c,t,dS,eqn)
{
real x[ND_ND];
real q,source,q1;
real y;
real dt;
dt=RP_Get_Real("flow-time" ) ;
C_CENTROID(x,c,t);
y=x[1];
q=C_R(c,t)*H*w*(cosh(k*y+k*h)*cos(w*dt)/sinh(k*h)+3/8*H*k*cosh(2*k*y+2*k*h)*cos(2*w*dt)/(sinh(k*h)*sinh(k*h)*sinh(k*h)*sinh(k*h)))/deltax;
q1=C_R(c,t)*H*w*(k*sinh(k*y+k*h)*cos(w*dt)/sinh(k*h)+3/4*H*k*k*sinh(2*k*y+2*k*h)*cos(2*w*dt)/(sinh(k*h)*sinh(k*h)*sinh(k*h)*sinh(k*h)))/deltax;
source=C_R(c,t)*C_V(c,t)*q+1/3*C_R(c,t)*C_MU_EFF(c,t)*q1;
dS[eqn]=C_R(c,t)*q;
return source;
}

[ Last edited by 1592203609 on 2014-3-20 at 17:59 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

lx478404227

木虫 (正式写手)

同求啊,我要模拟波浪对管道的作用
2楼2015-12-17 23:18:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 summer_cai 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:47 by vs90ilomwc
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 5/250 2026-02-08 08:32 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 4/200 2026-02-08 08:27 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:26 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:12 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:07 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 08:06 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:52 by vs90ilomwc
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:46 by vs90ilomwc
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:32 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:27 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:26 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:07 by vs90ilomwc
[硕博家园] 博士延得我,科研能力直往上蹿 +8 偏振片 2026-02-02 8/400 2026-02-08 06:52 by liyeqik
[教师之家] 有院领导为了换新车,用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人) +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
信息提示
请填处理意见