 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2467 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
dqyzgyttro金蟲 (小有名氣)
|
[求助]
求助利用uds求解方程組，uds發(fā)散 已有2人參與
|
|
|
在收斂的流場中設(shè)置四個uds，求解四個方程，算的是穩(wěn)態(tài)問題，不計算流場，只計算四個標(biāo)量。 現(xiàn)在我按照uds的格式寫好各項，其中對流項和源項都有利用C_UDSI引用其他三個標(biāo)量，同時還用ON_DEMAND給各標(biāo)量賦了初值。 代碼可以通過編譯，問題是開始計算就會顯示uds發(fā)散，不知道問題出在哪里，向各位高手求助。 謝謝解答！@wuming524 |
回復(fù)此樓» 猜你喜歡
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有3人回復(fù)
-
332求調(diào)劑
已經(jīng)有3人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有7人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有17人回復(fù)
-
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑
已經(jīng)有10人回復(fù)
-
0817調(diào)劑
已經(jīng)有3人回復(fù)
-
焦慮
已經(jīng)有13人回復(fù)
dqyzgyttro
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1673.6
- 帖子: 101
- 在線: 72小時
- 蟲號: 4389334
- 注冊: 2016-01-31
- 專業(yè): 多相流熱物理學(xué)
7樓2016-11-04 12:41:47
qc12345
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 146 (高中生)
- 金幣: 13666.3
- 紅花: 26
- 帖子: 3294
- 在線: 699.2小時
- 蟲號: 3071487
- 注冊: 2014-03-19
- 專業(yè): 園藝學(xué)與植物營養(yǎng)學(xué)
2樓2016-11-03 22:06:45
dqyzgyttro
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1673.6
- 帖子: 101
- 在線: 72小時
- 蟲號: 4389334
- 注冊: 2016-01-31
- 專業(yè): 多相流熱物理學(xué)
|
程序如下 #include "udf.h" #include "mem.h" #include "math.h" DEFINE_ON_DEMAND(demandbeginning) { Thread *t; cell_t c; face_t f; Domain *d; d=Get_Domain(1); thread_loop_c(t,d) { begin_c_loop(c,t) { C_UDSI(c,t,0)=0.000001; C_UDSI(c,t,1)=C_U(c,t); C_UDSI(c,t,2)=C_V(c,t); C_UDSI(c,t,3)=C_W(c,t); } end_c_loop(c,t); } } DEFINE_UDS_FLUX(flux0, f, t, i) { Thread *t0,*t1=NULL; cell_t c0,c1=-1; real NV_VEC(psi_vec),NV_VEC(A),flux; NV_D(psi_vec,=,0,0,0); c0=F_C0(f,t); t0=F_C0_THREAD(f,t); F_AREA(A,f,t); if(NULL==F_C1_THREAD(f,t)) { NV_DS(psi_vec,=,C_UDSI(c0,t0,1),C_UDSI(c0,t0,2),C_UDSI(c0,t0,3),*,998); flux=NV_DOT(psi_vec,A); } else { c1=F_C1(f,t); t1=F_C1_THREAD(f,t); NV_DS(psi_vec,=,C_UDSI(c0,t0,1),C_UDSI(c0,t0,2),C_UDSI(c0,t0,3),*,988); NV_DS(psi_vec,+=,C_UDSI(c1,t1,1),C_UDSI(c1,t1,2),C_UDSI(c1,t1,3),*,0.5); flux=NV_DOT(psi_vec,A); } return flux; } DEFINE_UDS_FLUX(flux1, f, t, i) { Thread *t0,*t1=NULL; cell_t c0,c1=-1; real NV_VEC(psi_vec),NV_VEC(A),flux1; NV_D(psi_vec,=,0,0,0); c0=F_C0(f,t); t0=F_C0_THREAD(f,t); F_AREA(A,f,t); if(NULL==F_C1_THREAD(f,t)) { NV_DS(psi_vec,=,C_UDSI(c0,t0,1),C_UDSI(c0,t0,2),C_UDSI(c0,t0,3),*,998); flux1=C_UDSI(c0,t0,0)*NV_DOT(psi_vec,A); } else { c1=F_C1(f,t); t1=F_C1_THREAD(f,t); NV_DS(psi_vec,=,C_UDSI(c0,t0,1),C_UDSI(c0,t0,2),C_UDSI(c0,t0,3),*,988); NV_DS(psi_vec,+=,C_UDSI(c1,t1,1),C_UDSI(c1,t1,2),C_UDSI(c1,t1,3),*,0.5); flux1=(C_UDSI(c0,t0,0)+C_UDSI(c1,t1,0))*0.5*NV_DOT(psi_vec,A); } return flux1; } DEFINE_SOURCE(myudssourceu,c,t,dS,eqn) { real x[ND_ND]; real CD,Re,force,k,source,a,b; Domain *d; d=Get_Domain(1); C_CENTROID(x,c,t); a=pow((C_U(c,t)*C_U(c,t)+C_V(c,t)*C_V(c,t)+C_W(c,t)*C_W(c,t)),0.5); b=pow((C_UDSI(c,t,1)*C_UDSI(c,t,1)+C_UDSI(c,t,2)*C_UDSI(c,t,2)+C_UDSI(c,t,3)*C_UDSI(c,t,3)),0.5); Re=998*pow((a*a-2*a*b+b*b),0.5)*0.00002/0.0000179; if (Re<=1000.0) { CD=24*(1+0.15*pow(Re,0.687))/Re; } else { CD=0.44; } force=CD*Re/24; k=18*0.0000179*force/998/0.00002/0.00002; source=998*C_UDSI(c,t,0)*k*(C_U(c,t)-C_UDSI(c,t,1)); dS[eqn]=-998*C_UDSI(c,t,0)*k; return source; } DEFINE_SOURCE(myudssourcev,c,t,dS,eqn) { real x[ND_ND]; real CD,Re,force,k,source,a,b; Domain *d; d=Get_Domain(1); C_CENTROID(x,c,t); a=pow((C_U(c,t)*C_U(c,t)+C_V(c,t)*C_V(c,t)+C_W(c,t)*C_W(c,t)),0.5); b=pow((C_UDSI(c,t,1)*C_UDSI(c,t,1)+C_UDSI(c,t,2)*C_UDSI(c,t,2)+C_UDSI(c,t,3)*C_UDSI(c,t,3)),0.5); Re=998*pow((a*a-2*a*b+b*b),0.5)*0.00002/0.0000179; if (Re<=1000.0) { CD=24*(1+0.15*pow(Re,0.687))/Re; } else { CD=0.44; } force=CD*Re/24; k=18*0.0000179*force/998/0.00002/0.00002; source=998*C_UDSI(c,t,0)*k*(C_V(c,t)-C_UDSI(c,t,1)); dS[eqn]=-998*C_UDSI(c,t,0)*k; return source; } DEFINE_SOURCE(myudssourcew,c,t,dS,eqn) { real x[ND_ND]; real CD,Re,force,k,source,a,b; Domain *d; d=Get_Domain(1); C_CENTROID(x,c,t); a=pow((C_U(c,t)*C_U(c,t)+C_V(c,t)*C_V(c,t)+C_W(c,t)*C_W(c,t)),0.5); b=pow((C_UDSI(c,t,1)*C_UDSI(c,t,1)+C_UDSI(c,t,2)*C_UDSI(c,t,2)+C_UDSI(c,t,3)*C_UDSI(c,t,3)),0.5); Re=998*pow((a*a-2*a*b+b*b),0.5)*0.00002/0.0000179; if (Re<=1000.0) { CD=24*(1+0.15*pow(Re,0.687))/Re; } else { CD=0.44; } force=CD*Re/24; k=18*0.0000179*force/998/0.00002/0.00002; source=998*C_UDSI(c,t,0)*k*(C_W(c,t)-C_UDSI(c,t,3)); dS[eqn]=-998*C_UDSI(c,t,0)*k; return source; } |
3樓2016-11-03 22:18:08
匿名
用戶注銷 (小有名氣)
- 應(yīng)助: 31 (小學(xué)生)
- 金幣: 1248.2
- 散金: 68
- 紅花: 7
- 帖子: 220
- 在線: 78.3小時
- 蟲號: 0
- 注冊: 2015-01-18
- 專業(yè): 工程熱物理與能源利用
4樓2016-11-04 07:14:56
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研]
 271材料工程求調(diào)劑
+5
|
.6lL 2026-03-18 | 5/250 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +3 | 何潤采123 2026-03-18 | 3/150 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 10/500 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研]
0854，計算機(jī)類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|