24小时热门版块排行榜    

查看: 2078  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 yunziaaaaa 2026-03-01 4/200 2026-03-01 23:45 by L135790
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-01 5/250 2026-03-01 23:31 by L135790
[考研] 0854复试调剂 276 +3 wmm9 2026-03-01 3/150 2026-03-01 23:13 by 热情沙漠
[考研] 265分求调剂不调专业和学校有行学上就 +6 礼堂丁真258 2026-02-28 8/400 2026-03-01 22:50 by jian_
[考研] 274求调剂 +3 cgyzqwn 2026-03-01 6/300 2026-03-01 21:24 by cgyzqwn
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 12/600 2026-03-01 19:48 by 好好好1233
[考研] 306分材料调剂 +4 chuanzhu川烛 2026-03-01 5/250 2026-03-01 19:48 by 无际的草原
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[考研] 313求调剂 +3 水流年lc 2026-02-28 3/150 2026-03-01 16:01 by 新能源达人
[考研] 307求调剂 +5 wyyyqx 2026-03-01 5/250 2026-03-01 15:21 by Fff-1
[考研] 284求调剂 +6 天下熯 2026-02-28 6/300 2026-03-01 14:19 by Ducount.Y
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见