| 24小時(shí)熱門(mén)版塊排行榜 |
| 2 | 1/1 | 返回列表 |
| 查看: 2089 | 回復(fù): 1 | ||
summer_cai新蟲(chóng) (初入文壇)
|
[求助]
【已完結(jié)】關(guān)于fluent模擬波浪的問(wèn)題
|
|
本人是fluent新手,想利用udf模擬波浪,下面是我編的源項(xiàng)程序,在運(yùn)行的時(shí)候出現(xiàn)overflow的錯(cuò)誤,運(yùn)行不下去,想請(qǐng)教各位大神看看是不是我編的程序錯(cuò)了,拜托了 #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 ] |
木蟲(chóng) (正式寫(xiě)手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|