 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 712 | 回復(fù): 1 | ||
| 【懸賞金幣】回答本帖問(wèn)題，作者王大大天將贈(zèng)送您 15 個(gè)金幣 | ||
王大大天新蟲(chóng) (小有名氣)
|
[求助]
做CFD-PBM耦合的，不導(dǎo)入uDF的時(shí)候沒(méi)有問(wèn)題，導(dǎo)入以后出現(xiàn)了這樣的問(wèn)題，，萬(wàn)分感謝
|
|
|
導(dǎo)入以后初始化，開(kāi)始計(jì)算還沒(méi)計(jì)算一步，就出現(xiàn)了這個(gè)問(wèn)題 received a fatal signal（segmentation fault） 如果有大佬知道 請(qǐng)求大佬幫幫我 萬(wàn)分感謝: |
回復(fù)此樓» 猜你喜歡
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有98人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有5人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生
已經(jīng)有2人回復(fù)
王大大天
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 115
- 散金: 20
- 帖子: 55
- 在線: 15.7小時(shí)
- 蟲(chóng)號(hào): 18102829
- 注冊(cè): 2019-10-10
- 專業(yè): 無(wú)機(jī)化工
|
程序如下#include "udf.h" #include "math.h" #include "stdio.h" #include "mem.h" #include "global.h" #define total_N 31 #define delta_L 2.00E-06 #define Kvvv 0.866 #define Ka 5.196 #define Growth_Kg 21 #define Growth_Ea 41.6 #define Growth_R 0.008314 /* k J/mol */ #define Growth_g 1.9 #define minmod_theta 1.5 #define rho_c 1296 /* 僅在擬均相下使用 */ #define rho_l 789 #define M_paracetamol 151.17 /* g/mol */ #define M_ethanol 46.07 /* g/mol */ #define Diffu 7.09e-10 /* m2/s Stokes-Einstein assume constant T */ #define Mean_L 0.0001 /* use for the calculation of mass transfer Kd */ #define Dynam_visc 0.001201 /* Dynamic viscosity, N*s/m2 or Pa*s */ #define Nucleation_Kn 1.597e-10 /* primary nucleation kn #/(min*m3) */ #define Nucleation_n 2.276 /* primary nucleation n */ #define Nucleation_E 7e-20 /* can be regressed from later data(Ea,kg,g,e,gamma)*/ #define Nucleation_gamma 2.5e-3 /* surface tension J */ #define Diameter_molecular 7.18e-10 /* m */ #define k_Boltzmann 1.381e-23 /* J/K */ #define diffu_m 0 /* 擴(kuò)散系數(shù)，此時(shí)為粒子，設(shè)為0 */ real source; real G_j1; real G_j2; real S_eq; real delta_S; real delta_C; real relative_S; real Nucleation_rate1; /* primary nucleation rate */ real Nucleation_rate2; /* second nucleation rate */ real Kd; real Moment_2; /* uds kg/kg */ real crystal_L[total_N] = { 0, 2.00E-06, 4.00E-06, 6.00E-06, 8.00E-06, 1.00E-05, 1.20E-05, 1.40E-05, 1.60E-05, 1.80E-05, 2.00E-05, 2.20E-05, 2.40E-05, 2.60E-05, 2.80E-05, 3.00E-05, 3.20E-05, 3.40E-05, 3.60E-05, 3.80E-05, 4.00E-05, 4.20E-05, 4.40E-05, 4.60E-05, 4.80E-05, 5.00E-05, 5.20E-05, 5.40E-05, 5.60E-05, 5.80E-05, 6.00E-05 }; /* 本udf平均離散*/ real min22(real a,real b) { if (a>b) return b; else return a; } real max22(real a,real b) { if (a<b) return b; else return a; } real minmod(real minmod_0,real minmod_1,real minmod_2) { real c; if (minmod_0>0 && minmod_1>0 && minmod_2>0) { c=min22(minmod_0,minmod_1); return min22(c,minmod_2); } else if (minmod_0<0 && minmod_1<0 && minmod_2<0) { c=max22(minmod_0,minmod_1); return max22(c,minmod_2); } else return 0; } DEFINE_SOURCE(fn_source,c,t,dS,eqn) { int i; real fdensitya; real fdensityb; real fdensityc; real fdensityd; real fdensitye; int k; int j=eqn-570+1; S_eq=0.0002955*exp(0.02179*C_T(c,t)); relative_S=C_YI(c,t,0)/S_eq; delta_S=C_YI(c,t,0)-S_eq; delta_C=delta_S*rho_l/M_paracetamol; if (delta_C>0) { G_j1=Growth_Kg*exp(-Growth_Ea/Growth_R/C_T(c,t))*pow(delta_C,Growth_g); G_j2=G_j1; } else { Kd=Diffu/Mean_L*(2+0.8*pow(C_D(c,t)*Mean_L*Mean_L*Mean_L*Mean_L/Dynam_visc/Dynam_visc/Dynam_visc,0.2)*pow(Dynam_visc/Diffu/rho_l,1/3)); G_j1 = Kd*Ka*M_paracetamol/3/Kvvv/rho_c*delta_C; G_j2=G_j1; } if (j==0||j==1) { Moment_2=0; for (k=0;k<(total_N-1);k++) {Moment_2+=rho_l*C_UDSI(c,t,k)*4/rho_c/Kvvv/(pow(crystal_L[k+1],4)-pow(crystal_L[k],4))*pow(crystal_L[k],2)*delta_L; } if ((C_YI(c,t,0)-S_eq) > 0) { Nucleation_rate1 = Nucleation_Kn*pow((C_YI(c, t, 0)-S_eq),Nucleation_n)/60; Nucleation_rate2 = Nucleation_E*Ka*Moment_2*Diffu/pow(Diameter_molecular,4)*exp(-3.1415*pow(Nucleation_gamma*Diameter_molecular*Diameter_molecular/k_Boltzmann/C_T(c,t),2)/log(relative_S)); } else { Nucleation_rate2 = 0; /* #/(m3*s) */ Nucleation_rate1 = 0; /* #/(m3*s) */ } } else { Nucleation_rate2=0; Nucleation_rate1=0; } if (j>=2 && j<=total_N-4) { fdensitya=C_UDSI(c,t,j-2)*4/rho_c/Kvvv/(pow(crystal_L[j-1],4)-pow(crystal_L[j-2],4)); fdensityb=C_UDSI(c,t,j-1)*4/rho_c/Kvvv/(pow(crystal_L[j],4)-pow(crystal_L[j-1],4)); fdensityc=C_UDSI(c,t,j)*4/rho_c/Kvvv/(pow(crystal_L[j+1],4)-pow(crystal_L[j],4)); fdensityd=C_UDSI(c,t,j+1)*4/rho_c/Kvvv/(pow(crystal_L[j+2],4)-pow(crystal_L[j+1],4)); fdensitye=C_UDSI(c,t,j+2)*4/rho_c/Kvvv/(pow(crystal_L[j+3],4)-pow(crystal_L[j+2],4)); } else if (j==1) { fdensitya=0; fdensityb=C_UDSI(c,t,j-1)*4/rho_c/Kvvv/(pow(crystal_L[j],4)-pow(crystal_L[j-1],4)); fdensityc=C_UDSI(c,t,j)*4/rho_c/Kvvv/(pow(crystal_L[j+1],4)-pow(crystal_L[j],4)); fdensityd=C_UDSI(c,t,j+1)*4/rho_c/Kvvv/(pow(crystal_L[j+2],4)-pow(crystal_L[j+1],4)); fdensitye=C_UDSI(c,t,j+2)*4/rho_c/Kvvv/(pow(crystal_L[j+3],4)-pow(crystal_L[j+2],4)); } else if (j==0) { fdensitya=0; fdensityb=0; fdensityc=C_UDSI(c,t,j)*4/rho_c/Kvvv/(pow(crystal_L[j+1],4)-pow(crystal_L[j],4)); fdensityd=C_UDSI(c,t,j+1)*4/rho_c/Kvvv/(pow(crystal_L[j+2],4)-pow(crystal_L[j+1],4)); fdensitye=C_UDSI(c,t,j+2)*4/rho_c/Kvvv/(pow(crystal_L[j+3],4)-pow(crystal_L[j+2],4)); } else if (j==total_N-3) { fdensitya=C_UDSI(c,t,j-2)*4/rho_c/Kvvv/(pow(crystal_L[j-1],4)-pow(crystal_L[j-2],4)); fdensityb=C_UDSI(c,t,j-1)*4/rho_c/Kvvv/(pow(crystal_L[j],4)-pow(crystal_L[j-1],4)); fdensityc=C_UDSI(c,t,j)*4/rho_c/Kvvv/(pow(crystal_L[j+1],4)-pow(crystal_L[j],4)); fdensityd=C_UDSI(c,t,j+1)*4/rho_c/Kvvv/(pow(crystal_L[j+2],4)-pow(crystal_L[j+1],4)); fdensitye=fdensityc; } else if (j==total_N-2) { fdensitya=C_UDSI(c,t,j-2)*4/rho_c/Kvvv/(pow(crystal_L[j-1],4)-pow(crystal_L[j-2],4)); fdensityb=C_UDSI(c,t,j-1)*4/rho_c/Kvvv/(pow(crystal_L[j],4)-pow(crystal_L[j-1],4)); fdensityc=C_UDSI(c,t,j)*4/rho_c/Kvvv/(pow(crystal_L[j+1],4)-pow(crystal_L[j],4)); fdensityd=fdensityc; fdensitye=fdensityc; } if (G_j1>=0) source=rho_l*rho_c*Kvvv/4/delta_L*(pow(crystal_L[j+1],4)-pow(crystal_L[j],4))* (-G_j2*(fdensityc+0.5*delta_L*minmod(minmod_theta/delta_L*(fdensityc-fdensityb),0.5/delta_L*(fdensityd - fdensityb), minmod_theta/delta_L*(fdensityd - fdensityc)))+G_j1*(fdensityb+0.5*delta_L*minmod(minmod_theta/delta_L*(fdensityb-fdensitya),0.5/delta_L*(fdensityc-fdensitya), minmod_theta/delta_L*(fdensityc - fdensityb)))+Nucleation_rate2/rho_l/2+Nucleation_rate1/rho_l/2); else if (G_j1<0) source=rho_l*rho_c*Kvvv/4/delta_L*(pow(crystal_L[j+1],4)-pow(crystal_L[j],4))* (-G_j2*(fdensityd-0.5*delta_L*minmod(minmod_theta/delta_L*(fdensityd-fdensityc),0.5/delta_L*(fdensitye -fdensityc), minmod_theta/delta_L*(fdensitye - fdensityd)))+G_j1*(fdensityc-0.5*delta_L*minmod(minmod_theta/delta_L*(fdensityc-fdensityb),0.5/delta_L*(fdensityd-fdensityb), minmod_theta/delta_L*(fdensityd -fdensityc)))); dS[eqn] = 0; return source; |
2樓2020-07-23 10:59:29
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 王大大天 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過(guò) a 線，總分 281 求調(diào)劑 +4 | 香香 Lu 2026-03-03 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 22408-273求調(diào)劑-擔(dān)任3個(gè)項(xiàng)目的負(fù)責(zé)人-1篇國(guó)際期刊論文（一作）1篇核心期刊論文在投。 +3 | 沒(méi)想好取什么名 2026-03-03 | 3/150 |
|
|
|
[考研] 085600材料與化工（高分子）290分求調(diào)劑 +8 | wengyujian 2026-03-04 | 8/400 |
|
|
|
[考研] 0856材料與化工求調(diào)劑！ +3 | 化工考生111 2026-03-04 | 7/350 |
|
|
|
[考研] 0857求調(diào)劑 +4 | CDDD1013 2026-03-04 | 4/200 |
|
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
 5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
|
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +8 | yang婷 2026-03-02 | 10/500 |
|
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
|
[考研] 化工專碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
|
[考研]
材料270求調(diào)劑
6+6
|
Eiiiio 2026-03-01 | 11/550 |
|
|
|
[考研]
 計(jì)算機(jī)學(xué)碩 分?jǐn)?shù)285求調(diào)劑
+4
|
glwshine 2026-03-02 | 5/250 |
|
|
|
[考研] 284求調(diào)劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
|
[考研] 1 +7 | 黑！在干嘛 2026-02-28 | 8/400 |
|
|
|
[考研] 272求調(diào)劑 +9 | 材紫有化 2026-02-28 | 9/450 |
|
|
|
[考研] 一志愿東北大學(xué)材料專碩328，求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|
|
|
[論文投稿]
求助coordination chemistry reviews 的寫(xiě)作模板
10+3
|
ljplijiapeng 2026-02-27 | 4/200 |
|
|
|
[論文投稿]
Optics letters投稿被拒求助
30+3
|
luckyry 2026-02-26 | 4/200 |
|