| 查看: 534 | 回复: 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); } 获取每个网格入口边界处的压力 然后有对应的速度公式 然后通过计算得到质量流率 再返回给这个网格，要循环每个网格 ，但是现在加载进去就给质量流率的边界条件返回了一个网格上的质量流率，求大神知道 哪里里写错了 |
回复此楼» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有6人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有69人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复