| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1277 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
豬的主人新蟲(chóng) (初入文壇)
|
[求助]
udf致命錯(cuò)誤 已有1人參與
|
|
|
我仿真的模型是表層土壤淺埋換熱管,研究非飽和土壤中熱濕耦合 遷移,現(xiàn)階段只考慮液態(tài)水的遷移。模型有混合mixture模型,porous模型,利用udf編寫(xiě)方程,程序如下: #include"udf.h" #include<math.h> DEFINE_UDS_UNSTEADY(myunsteady,c,t,i,apu,su) { real physical_dt,vol,phi_old; physical_dt=RP_Get_Real("physical_time_step" ;vol=C_VOLUME(c,t); *apu=-vol/physical_dt; phi_old=C_STORAGE_R(c,t,SV_UDSI_M1(i)); *su=vol*phi_old/physical_dt; } DEFINE_SOURCE(source_moving,cell,sec_th,dS,eqn) { real A[ND_ND]; real source; real temp0,temp1,temp2; real vol; real vof; real a0[ND_ND],a1[ND_ND],a2[ND_ND]; real diff; real xv[ND_ND],yv[ND_ND],zv[ND_ND]; Domain*d; Thread*mix_th; face_t f; d=Get_Domain(2); mix_th=THREAD_SUPER_THREAD(sec_th); //指針指向混合項(xiàng)的第二項(xiàng)即液態(tài)水 temp0=C_T_G(cell,mix_th)[0]; temp1=C_T_G(cell,mix_th)[1]; temp2=C_T_G(cell,mix_th)[2]; vol=C_VOLUME(cell,mix_th); //控制體體積; vof=C_UDSI(cell,mix_th,0); //控制體 內(nèi)含濕量; diff=-2.3495*pow(10,-5)*pow(vof,6)+2.0370*pow(10,-5)*pow(vof,5)-6.6946*pow(10,-6)*pow(vof,4)+1.0438*pow(10,-6)*pow(vof,3)-7.8699*pow(10,-8)*pow(vof,2)+2.7350*pow(10,-9)*vof+5.4039*pow(10,-12); //等溫?cái)U(kuò)散系數(shù) xv[0]=1;xv[1]=0;xv[2]=0; yv[0]=0;yv[1]=1;yv[2]=0; zv[0]=0;zv[1]=0;zv[2]=1; begin_f_loop(f,mix_th) //歷遍混合相的 所有面; { F_AREA(A,f,mix_th); a0[ND_ND]+=temp0*NV_DOT(A,xv); //x方向的溫度梯度與控制體x方向 投影面積的乘積; a1[ND_ND]+=temp1*NV_DOT(A,yv); a2[ND_ND]+=temp2*NV_DOT(A,zv); } end_f_loop(f,mix_th) source=diff*(a0[0]+a1[1]+a2[2])/vol; dS[eqn]=0; return source; } DEFINE_DIFFUSIVITY(MyDiff,c,t,i) { real vof; real con_zl; real con_nd; real MyDiff; vof=C_UDSI(c,t,0); con_zl=0.075636667-0.000114214*C_T(c,t)-0.000000279*C_T(c,t)*C_T(c,t); con_nd=0.00041212*exp(-C_T(c,t)/12.37666)+0.00119916*exp(-C_T(c,t)/41.70239)+0.00018058; MyDiff=9.15003*pow(vof,3.613655)*con_zl*0.0014953452/(con_nd*0.073240827); return MyDiff; } 在進(jìn)行計(jì)算的時(shí)候顯示致命錯(cuò)誤,求大神指點(diǎn) @wuming524 |
新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +3 | zzz,,r 2026-04-02 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 0832食品科學(xué)與工程學(xué)碩282調(diào)劑 +4 | 魚(yú)在水中游a 2026-04-02 | 7/350 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-04-02 | 4/200 |
|
|
[考研] 266分,一志愿電氣工程,本科材料,求材料專(zhuān)業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-04-01 | 11/550 |
|
|
[考研] 一志愿北交大材料工程總分358 +8 | cs0106 2026-04-01 | 9/450 |
|
|
[考研] 07生物學(xué)求調(diào)劑 一志愿同濟(jì)大學(xué)359分 +3 | LAMC. 2026-03-30 | 3/150 |
|
|
[考研] 311求調(diào)劑 +16 | zchqwer 2026-04-01 | 18/900 |
|
|
[考研] 303求調(diào)劑 +8 | DLkz1314. 2026-03-30 | 8/400 |
|
|
[考研] 301求調(diào)劑 +13 | A_JiXing 2026-04-01 | 13/650 |
|
|
[考研] 一志愿安徽大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)碩,331分求調(diào)劑 +5 | 蔣昌鵬qtj 2026-04-01 | 5/250 |
|
|
[考研] 0856材料化工調(diào)劑 總分330 +18 | zhubinhao 2026-03-27 | 18/900 |
|
|
[考研] 求化學(xué)調(diào)劑 +12 | wulanna 2026-03-28 | 12/600 |
|
|
[考研]
|
小羊36 2026-03-30 | 4/200 |
|
|
[考研] 一志愿中海洋材料357 +4 | 麥恩莉. 2026-03-30 | 4/200 |
|
|
[考研] 一志愿西電085401數(shù)一英一299求調(diào)劑 六級(jí)521 +4 | 愛(ài)吃大鴨梨 2026-03-31 | 4/200 |
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專(zhuān)碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
[考研]
|
鐘llll 2026-03-26 | 5/250 |
|
|
[考研] 一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|