| 查看: 463 | 回复: 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); } 获取每个网格入口边界处的压力 然后有对应的速度公式 然后通过计算得到质量流率 再返回给这个网格，要循环每个网格 ，但是现在加载进去就给质量流率的边界条件返回了一个网格上的质量流率，求大神知道 哪里里写错了 |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有85人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有0人回复