| 2 | 1/1 | 返回列表 |
| 查看: 2090 | 回復(fù): 1 | ||
[求助]
【已完結(jié)】關(guān)于fluent模擬波浪的問題
|
|
本人是fluent新手,想利用udf模擬波浪,下面是我編的源項程序,在運行的時候出現(xiàn)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 ] |
木蟲 (正式寫手)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311(085601)求調(diào)劑 +6 | liziyeyeye 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 085600,321分求調(diào)劑 +7 | 大饞小子 2026-03-31 | 7/350 |
|
|
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +6 | GENJIOW 2026-03-31 | 6/300 |
|
|
[考研] 環(huán)境工程 085701,267求調(diào)劑 +13 | minht 2026-03-29 | 13/650 |
|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +6 | 12.27 2026-03-27 | 6/300 |
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 13/650 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +10 | 傾____蕭 2026-03-27 | 11/550 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +8 | 哇呼哼呼哼 2026-03-27 | 10/500 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-30 | 6/300 |
|
|
[論文投稿]
chinese chemical letters英文版投稿求助
130+3
|
Yishengeryi 2026-03-30 | 3/150 |
|
|
[考研] 323分 食品與營養(yǎng)調(diào)劑 +3 | 嘿ooo 2026-03-31 | 3/150 |
|
|
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +5 | pipiver 2026-03-30 | 5/250 |
|
|
[考研] 317求調(diào)劑 +8 | 十閑wx 2026-03-24 | 8/400 |
|
|
[考研] 材料化工340求調(diào)劑 +3 | jhx777 2026-03-30 | 3/150 |
|
|
[考研] 材料專碩 085600求調(diào)劑 +7 | BBQ233 2026-03-30 | 7/350 |
|
|
[考研] 342求調(diào)劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 一志愿上海理工能源動力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|