| 查看: 733 | 回复: 0 | ||
[求助]
固体壁面吸附甲醛udf编写源项
|
|
各位：我的问题是壁面吸附甲醛。 壁面吸附函数为S=(C(x,t)-C(x,t-1))*Ks/delta time，这里的浓度是用隐式格式，即此时和上一时刻浓度的差值。需要指出的是，只有在近壁面的时候才发生吸附作用，其他区域S=0.。 我想通过组分输运方程附加源项来实现，以下是我的udf，望各位不吝赐教，本人不甚感激 #include"udf.h" #define Ks=1000.0 /*质源表达式*/ DEFINE_SOURCE(mass_source,f,t,dS,eqn) { real physical_dt,phi_old,con,source=0.0; int zone_ID=22;/*壁面ID*/ cell_t c0; face_t f; Thread *t0; Thread *inletstrip=Lookup_Thread(Get_Domain(1),zone_ID); c0=F_C0(f,inletstrip); t0=THREAD_T0(inletstrip); /*壁面指针*/ /*begin_f_loop(f,t0) {*/ physical_dt=RP_Get_Real("physical-time-step"  ;/*求出delta time*/ phi_old=C_STORAGE_R(f,t0,SV_Yi_M1(i)); /*上一时刻的组分分数*/ con=0.5*Ks*(F_Yi(f,t0)/*-phi_old*/); source=2.0*con/physical_dt; /*} end_f_loop(f,t0) */ dS[eqn]=0.0; return source; } @月只蓝 @beefly 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有101人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有0人回复