| 查看: 785 | 回复: 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客户端 |
回复此楼» 猜你喜欢
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有205人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有3人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
5