| 查看: 501 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者李小白007将赠送您 30 个金币 | ||
[求助]
固体壁面吸附甲醛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客户端 |
回复此楼» 猜你喜欢
郑州大学光电器件方向招收一名博士
已经有20人回复
-
氧化硅薄膜玻璃 求购
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有216人回复
-
vasp极化强度
已经有0人回复
-
vasp计算的一些问题
已经有0人回复
-
国科温州研究院 项目聘用有意“二战”的优秀本科毕业生
已经有26人回复
-
SCAPS-1D学习交流
已经有0人回复
20