| 查看: 1591 | 回复: 9 | ||
[求助]
如何在不同的壁面用不同的壁面函数 已有2人参与
|
|
一个球体，分左半球和右半球，左半球用一个壁面函数，右半球用另一个壁面函数，如何实现？UDF都比较简单，一个是HELP文件里的原文，另一个就是仅仅增加了一个系数2. 2个UDF分别为： UDF1： #include "udf.h" DEFINE_WALL_FUNCTIONS(USER,f,t,c0,t0,wf_ret,yPlus,Emod) { real wf_value switch (wf_ret) { case UPLUS_LAM: wf_value= yPlus; break; case UPLUS_TRB: wf_value= log(Emod*yPlus)/KAPPA; break; case DUPLUS_LAM: wf_value = 1.0; break; case DUPLUS_TRB: wf_value= 1./(KAPPA*yPlus); break; case D2UPLUS_TRB: wf_value = -1./(KAPPA*yPlus*yPlus); break; default: printf("Wall function return value unavailable\n"  ;} return wf_value; } UDF2： #include "udf.h" DEFINE_WALL_FUNCTIONS(USER,f,t,c0,t0,wf_ret,yPlus,Emod) { real wf_value switch (wf_ret) { case UPLUS_LAM: wf_value= 2*yPlus; break; case UPLUS_TRB: wf_value= log(2*Emod*yPlus)/KAPPA; break; case DUPLUS_LAM: wf_value = 2.0; break; case DUPLUS_TRB: wf_value= 1./(KAPPA*yPlus); break; case D2UPLUS_TRB: wf_value = -1./(KAPPA*yPlus*yPlus); break; default: printf("Wall function return value unavailable\n" ;} return wf_value; } |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 注册化工工程师须知
已经有23人回复
- 过饱和体系中的气核析出问题
已经有1人回复
2楼2015-08-06 12:25:11
3楼2015-08-06 23:52:13
4楼2015-08-07 09:13:03
5楼2015-08-07 10:12:16
6楼2015-08-07 10:27:30
7楼2015-08-07 11:18:29
8楼2015-08-07 14:38:28
ffengqiao
金虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 2224.6
- 散金: 213
- 红花: 5
- 帖子: 629
- 在线: 214.3小时
- 虫号: 1655125
- 注册: 2012-03-01
- 性别: GG
- 专业: 工程热物理相关交叉领域
9楼2015-08-07 22:21:32
10楼2015-12-23 21:31:52
