版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

dqyzgyttro

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1673.6
帖子: 101
在線: 72小時(shí)
蟲號(hào): 4389334
注冊(cè): 2016-01-31
專業(yè): 多相流熱物理學(xué)

[求助] 求助利用uds求解方程組，uds發(fā)散已有2人參與

在收斂的流場(chǎng)中設(shè)置四個(gè)uds，求解四個(gè)方程，算的是穩(wěn)態(tài)問(wèn)題，不計(jì)算流場(chǎng)，只計(jì)算四個(gè)標(biāo)量。
現(xiàn)在我按照uds的格式寫好各項(xiàng)，其中對(duì)流項(xiàng)和源項(xiàng)都有利用C_UDSI引用其他三個(gè)標(biāo)量，同時(shí)還用ON_DEMAND給各標(biāo)量賦了初值。
代碼可以通過(guò)編譯，問(wèn)題是開(kāi)始計(jì)算就會(huì)顯示uds發(fā)散，不知道問(wèn)題出在哪里，向各位高手求助。
謝謝解答！@wuming524

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有3人回復(fù)

1樓 2016-11-03 20:31:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

qc12345

至尊木蟲 (職業(yè)作家)

應(yīng)助: 146 (高中生)
金幣: 13419.3
紅花: 26
帖子: 3293
在線: 699.2小時(shí)
蟲號(hào): 3071487
注冊(cè): 2014-03-19
專業(yè): 園藝學(xué)與植物營(yíng)養(yǎng)學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

??UDS?????????????н??!

贊一下

回復(fù)此樓

聽(tīng)海

2樓2016-11-03 22:06:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dqyzgyttro

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1673.6
帖子: 101
在線: 72小時(shí)
蟲號(hào): 4389334
注冊(cè): 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;
}

贊一下

回復(fù)此樓

3樓2016-11-03 22:18:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

匿名

用戶注銷 (小有名氣)

應(yīng)助: 31 (小學(xué)生)
金幣: 1248.2
散金: 68
紅花: 7
帖子: 220
在線: 78.3小時(shí)
蟲號(hào): 0
注冊(cè): 2015-01-18
專業(yè): 工程熱物理與能源利用

感謝參與，應(yīng)助指數(shù) +1

本帖僅樓主可見(jiàn)

贊一下

4樓2016-11-04 07:14:56

已閱申請(qǐng)仿真EPI 回復(fù)此樓編輯查看我的主頁(yè)

dqyzgyttro

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1673.6
帖子: 101
在線: 72小時(shí)
蟲號(hào): 4389334
注冊(cè): 2016-01-31
專業(yè): 多相流熱物理學(xué)

引用回帖:

4樓: Originally posted by xpxp1991 at 2016-11-04 07:14:56
初值用define-init

用了，似乎問(wèn)題沒(méi)解決呢

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

5樓2016-11-04 10:28:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

匿名

用戶注銷 (小有名氣)

應(yīng)助: 31 (小學(xué)生)
金幣: 1248.2
散金: 68
紅花: 7
帖子: 220
在線: 78.3小時(shí)
蟲號(hào): 0
注冊(cè): 2015-01-18
專業(yè): 工程熱物理與能源利用

本帖僅樓主可見(jiàn)

贊一下

6樓2016-11-04 12:22:03

已閱申請(qǐng)仿真EPI 回復(fù)此樓編輯查看我的主頁(yè)

dqyzgyttro

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1673.6
帖子: 101
在線: 72小時(shí)
蟲號(hào): 4389334
注冊(cè): 2016-01-31
專業(yè): 多相流熱物理學(xué)

引用回帖:

6樓: Originally posted by xpxp1991 at 2016-11-04 12:22:03
Cell_t用法有問(wèn)題

是什么問(wèn)題呢，不太明白

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

7樓2016-11-04 12:41:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

蔡榮需努力

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 34.3
帖子: 21
在線: 5.6小時(shí)
蟲號(hào): 23281994
注冊(cè): 2020-07-29
專業(yè): 結(jié)構(gòu)工程

樓主，解決了么，我碰到相同的問(wèn)題了。希望可以交流指點(diǎn)！

贊一下(1人)

回復(fù)此樓

8樓2020-10-05 15:13:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 dqyzgyttro 的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽(yáng)CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無(wú)懈可擊111
[考研] 274求調(diào)劑 +5	S.H1 2026-03-18	5/250	2026-03-18 21:27 by guosr9609
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過(guò)，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 307求調(diào)劑 +3	冷笙123 2026-03-17	3/150	2026-03-18 09:55 by macy2011
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 283求調(diào)劑 +3	聽(tīng)風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321

24小時(shí)熱門版塊排行榜

[求助] 求助利用uds求解方程組，uds發(fā)散 已有2人參與

» 猜你喜歡

【答案】應(yīng)助回帖

[求助] 求助利用uds求解方程組，uds發(fā)散已有2人參與