| 查看: 759 | 回复: 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客户端 |
回复此楼» 猜你喜欢
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有80人回复
-
物理学I论文润色/翻译怎么收费?
已经有81人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
20