[求助] fluent多孔介质阻力系数是位置的函数，但不管阻力系数大还是小，速度压力都没明显区别

一个流体域里包含了一个多孔域，编了一个udf，想要多孔域实现类似格栅的效果，隔一段阻力系数大，孔隙率接近0，空气很难通过，隔一段阻力系数小，孔隙率接近1，空气顺利通过。 采用ergun公式算的阻力系数，  体积当量直径采用的是 6倍体积比表面积 采用的physical velocity 算的。 然后计算之后发现多孔域的速度压力云图，速度矢量图都没有体现出阻力系数大，阻力系数小的区别，但是却有压力降，似乎fluent自动给多孔域平均了各种阻力系数。 附上我的UDF:（有的时候孔隙率的设置能用UDF定义的，有的时候就float报错） 请大神指正，金币可以全部奉上。 #include "udf.h" DEFINE_PROFILE(porosity,t,i) {    real xc[ND_ND];    cell_t c;    begin_c_loop(c,t)       {      if (xc[0]>-44.1164 && xc[0]<-39.1164&&xc[0]>-54.1164 && xc[0]<-49.1164&&xc[0]>-64.1164 && xc[0]<-59.1164&&xc[0]>-74.1164 && xc[0]<-69.1164&&xc[0]>-84.1164 && xc[0]<-79.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)      {        F_PROFILE(c,t,i)=1;     }      else if (xc[0]>-39.1164 && xc[0]<-34.1164&&xc[0]>-49.1164 && xc[0]<-44.1164&&xc[0]>-59.1164 && xc[0]<-54.1164&&xc[0]>-69.1164 && xc[0]<-64.1164&&xc[0]>-79.1164 && xc[0]<-74.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)      {        F_PROFILE(c,t,i)=0;       }     } end_c_loop(c,t) } DEFINE_PROFILE(v_resistance,t,i) {    real xc[ND_ND];    cell_t c;    begin_c_loop(c,t)       {      if (xc[0]>-44.1164 && xc[0]<-39.1164&&xc[0]>-54.1164 && xc[0]<-49.1164&&xc[0]>-64.1164 && xc[0]<-59.1164&&xc[0]>-74.1164 && xc[0]<-69.1164&&xc[0]>-84.1164 && xc[0]<-79.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)      {        F_PROFILE(c,t,i)=0;     }      else if (xc[0]>-39.1164 && xc[0]<-34.1164&&xc[0]>-49.1164 && xc[0]<-44.1164&&xc[0]>-59.1164 && xc[0]<-54.1164&&xc[0]>-69.1164 && xc[0]<-64.1164&&xc[0]>-79.1164 && xc[0]<-74.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)      {        F_PROFILE(c,t,i)=664781714.9993;       }     } end_c_loop(c,t) } DEFINE_PROFILE(v_resist1,t,i) {    real xc[ND_ND];    cell_t c;    begin_c_loop(c,t)       {      if (xc[0]>-44.1164 && xc[0]<-39.1164&&xc[0]>-54.1164 && xc[0]<-49.1164&&xc[0]>-64.1164 && xc[0]<-59.1164&&xc[0]>-74.1164 && xc[0]<-69.1164&&xc[0]>-84.1164 && xc[0]<-79.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)      {        F_PROFILE(c,t,i)=0;     }      else if (xc[0]>-39.1164 && xc[0]<-34.1164&&xc[0]>-49.1164 && xc[0]<-44.1164&&xc[0]>-59.1164 && xc[0]<-54.1164&&xc[0]>-69.1164 && xc[0]<-64.1164&&xc[0]>-79.1164 && xc[0]<-74.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)      {        F_PROFILE(c,t,i)=664781714999.3;       }     } end_c_loop(c,t) } DEFINE_PROFILE(i_resistance,t,i) {    real xc[ND_ND];    cell_t c;    begin_c_loop(c,t)       {      if (xc[0]>-44.1164 && xc[0]<-39.1164&&xc[0]>-54.1164 && xc[0]<-49.1164&&xc[0]>-64.1164 && xc[0]<-59.1164&&xc[0]>-74.1164 && xc[0]<-69.1164&&xc[0]>-84.1164 && xc[0]<-79.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)      {        F_PROFILE(c,t,i)=0;     }      else if (xc[0]>-39.1164 && xc[0]<-34.1164&&xc[0]>-49.1164 && xc[0]<-44.1164&&xc[0]>-59.1164 && xc[0]<-54.1164&&xc[0]>-69.1164 && xc[0]<-64.1164&&xc[0]>-79.1164 && xc[0]<-74.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)      {        F_PROFILE(c,t,i)=7368209.3023;       }     } end_c_loop(c,t) } DEFINE_PROFILE(i_resist1,t,i) {    real xc[ND_ND];    cell_t c;    begin_c_loop(c,t)       {      if (xc[0]>-44.1164 && xc[0]<-39.1164&&xc[0]>-54.1164 && xc[0]<-49.1164&&xc[0]>-64.1164 && xc[0]<-59.1164&&xc[0]>-74.1164 && xc[0]<-69.1164&&xc[0]>-84.1164 && xc[0]<-79.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)      {        F_PROFILE(c,t,i)=0;     }      else if (xc[0]>-39.1164 && xc[0]<-34.1164&&xc[0]>-49.1164 && xc[0]<-44.1164&&xc[0]>-59.1164 && xc[0]<-54.1164&&xc[0]>-69.1164 && xc[0]<-64.1164&&xc[0]>-79.1164 && xc[0]<-74.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)      {        F_PROFILE(c,t,i)=7368209302.3;       }     } end_c_loop(c,t) }

回复此楼