24小时热门版块排行榜

返回列表

summer_cai

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 671
帖子: 9
在线: 27.6小时
虫号: 2881761
注册: 2013-12-19
性别: MM

[求助] 【已完结】关于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 ]

回复此楼

» 猜你喜欢

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

1楼 2013-12-28 11:51:33

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

lx478404227

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 1651.5
散金: 212
帖子: 399
在线: 188.8小时
虫号: 4144218
注册: 2015-10-14
性别: GG
专业: 工程热物理相关交叉领域

同求啊，我要模拟波浪对管道的作用

赞一下

回复此楼

2楼2015-12-17 23:18:27

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

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

返回列表