| 查看: 477 | 回复: 0 | ||
[求助]
udf 质量流率 每个网格循环
|
|
#include"mem.h" #define p_pretotal 350000 #define pcr 185000 #define R 287.06 #define k 1.4 #define T 300 DEFINE_PROFILE(mass_flow_udf,t,i) { real p_sum = 0; real r = 0,tem; double V; real xc[ND_ND]; face_t f; Domain *d; d = Get_Domain(1); /*获取液态区域的网格 */ begin_f_loop(f,t) { F_CENTROID(xc,f,t); /* 将网格坐标存入数组xc中 */ p_sum = F_P(f, t); /* 将单元中的压力值赋予 p_sum中*/ tem= F_T(f, t); r=F_P(f, t)/(R*tem); if(p_sum<=p_pretotal) if(p_sum>=pcr) V=sqrt((2*k*R*T/(k-1))*(1-pow((p_sum/p_pretotal),(k-1)/k))); else V=sqrt(k*R*T*pow((pcr/p_pretotal),(k-1)/k)); else V=0; F_PROFILE(f,t,i)=0.0002*r*V; } end_f_loop(f,t); } 获取每个网格入口边界处的压力 然后有对应的速度公式 然后通过计算得到质量流率 再返回给这个网格，要循环每个网格 ，但是现在加载进去就给质量流率的边界条件返回了一个网格上的质量流率，求大神知道 哪里里写错了 |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有247人回复
-
散金币，求好运，祝面上顺利！
已经有78人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
