| 查看: 2082 | 回复: 1 | ||
[求助]
【已完结】关于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 ] |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
我在想关于fluent的一个问题。
已经有3人回复
- 圆柱绕流问题求助
已经有15人回复
- 用fluent模拟激波管的问题
已经有7人回复
- fluent模拟后处理
已经有4人回复
- fluent中patch区域问题
已经有21人回复
- fluent做太阳能电池模拟
已经有8人回复
- 求助：FLUENT进出口流量不守恒的问题
已经有19人回复
- fluent模拟喷雾过程相关问题求助
已经有34人回复
- fluent模拟动画慢动作
已经有7人回复
- fluent新手问题
已经有10人回复
- fluent 非预混燃烧模拟的问题
已经有7人回复
- 关于FLUENT仿真风力发电系统的问题
已经有25人回复
- 请教一个关于闭合回路环流的fluent模拟问题
已经有17人回复
- fluent遇到的问题
已经有4人回复
- 【求助】Fluent中关于设置稳态流动还是非稳态流动的问题
已经有6人回复
- 【求助】Fluent中的对称问题
已经有10人回复
lx478404227
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1651.5
- 散金: 212
- 帖子: 399
- 在线: 188.8小时
- 虫号: 4144218
- 注册: 2015-10-14
- 性别: GG
- 专业: 工程热物理相关交叉领域
2楼2015-12-17 23:18:27
6
