 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 713 | 回復: 1 | |||
| 【懸賞金幣】回答本帖問題，作者王大大天將贈送您 15 個金幣 | |||
[求助]
做CFD-PBM耦合的，不導入uDF的時候沒有問題，導入以后出現(xiàn)了這樣的問題，，萬分感謝
|
|||
|
導入以后初始化，開始計算還沒計算一步，就出現(xiàn)了這個問題 received a fatal signal（segmentation fault） 如果有大佬知道 請求大佬幫幫我 萬分感謝: |
回復此樓» 猜你喜歡
重慶交通大學光子學微結構與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復
-
一志愿鄭大材料學碩298分，求調(diào)劑
已經(jīng)有6人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有258人回復
-
尋合作：應力腐蝕多尺度模擬
已經(jīng)有6人回復
-
考研交流
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
幫我的英語口語老師找學生
已經(jīng)有3人回復
|
程序如下#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 /* 擴散系數(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
相關版塊跳轉(zhuǎn)
我要訂閱樓主 王大大天 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +9 | 一個學習者 2026-03-04 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +4 | ounce. 2026-03-04 | 5/250 |
|
|
|
[考研] 中科大材料299求調(diào)劑 +4 | DAIjiayo 2026-03-05 | 6/300 |
|
|
|
[基金申請]
 成果系統(tǒng)訪問量大，請15分鐘后再嘗試。由此給您造成的不便，敬請諒解。
+15
|
xhuama 2026-03-02 | 17/850 |
|
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 4/200 |
|
|
|
[考研] 376求調(diào)劑 +3 | 王浩然sam 2026-03-04 | 3/150 |
|
|
|
[考研] 一志愿武漢理工大學-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +4 | 熬夜的貓頭鷹 2026-03-02 | 4/200 |
|
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
 5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
|
[考研] 一志愿東北大學化學314分求調(diào)劑 +5 | lr1212.. 2026-03-02 | 5/250 |
|
|
|
[考研] 一志愿中科大能動297求調(diào)劑，本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
|
[考研] 338求調(diào)劑 +5 | 18162027187 2026-03-02 | 6/300 |
|
|
|
[考研] 0854復試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
|
[考研] 化工專碩342，一志愿大連理工大學，求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
|
[考研] 化工299分求調(diào)劑 一志愿985落榜 +5 | 嘻嘻(*^ω^*) 2026-03-01 | 5/250 |
|
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|