| 查看: 752 | 回复: 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客户端 |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有100人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有1人回复
-
0702一志愿吉大B区求调剂
已经有2人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
10