| 查看: 516 | 回复: 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); } 获取每个网格入口边界处的压力 然后有对应的速度公式 然后通过计算得到质量流率 再返回给这个网格，要循环每个网格 ，但是现在加载进去就给质量流率的边界条件返回了一个网格上的质量流率，求大神知道 哪里里写错了 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有92人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有50人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
30