版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

dqyzgyttro

金蟲 (小有名氣)

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

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

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

回復此樓

» 猜你喜歡

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

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

dqyzgyttro

金蟲 (小有名氣)

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

程序如下
#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

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

查看全部 8 個回答

qc12345

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

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

【答案】應助回帖

感謝參與，應助指數 +1

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

贊一下

回復此樓

聽海

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

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

匿名

用戶注銷 (小有名氣)

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

感謝參與，應助指數 +1

本帖僅樓主可見

贊一下

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

已閱申請仿真EPI 回復此樓編輯查看我的主頁

dqyzgyttro

金蟲 (小有名氣)

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

引用回帖:

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

用了，似乎問題沒解決呢

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

贊一下

回復此樓

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

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

查看全部 8 個回答

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

[考研] 281求調劑（0805） +6	煙汐憶海 2026-03-16	16/800	2026-03-19 00:59 by 煙汐憶海
[考研] 一志愿華中科技大學，080502，354分求調劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 一志愿吉林大學材料學碩321求調劑 +4	Ymlll 2026-03-18	6/300	2026-03-18 22:15 by li123456789.
[考研] 085600材料與化工調劑 324分 +8	llllkkkhh 2026-03-18	8/400	2026-03-18 21:01 by Catalysis25
[考研] 295求調劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 材料專碩274一志愿陜西師范大學求調劑 +6	薛云鵬 2026-03-13	6/300	2026-03-18 14:14 by 脫穎而出
[考研] 312求調劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 有沒有道鐵/土木的想調劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 材料與化工專碩調劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 275求調劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 0703 物理化學調劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 326求調劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 294求調劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 中科大材料與化工319求調劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 255求調劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 復試調劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 308 085701 四六級已過求調劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 266求調劑 +4	學員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 工科278分求調劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 0817化學工程與技術考研312分調劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

dqyzgyttro

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

» 猜你喜歡

dqyzgyttro

qc12345

【答案】應助回帖

匿名

dqyzgyttro

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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