| 查看: 1611 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
在壁面添加质量源项 已有1人参与
|
|||
|
请问如何在壁面添加质量源项呢？我想加一个-KC，就是一级反应动力学里的那个，使到壁面的污染物以-kc减少，我看了很久，知道需要提取临近壁面的cell，可是，如何表示方程式和用壁面的污染物质量分数呢？我实在是憋不出来了，UDF最多写成这个样子，求助虫友们帮我完善UDF，谢谢，谢谢！！！其中一级反应动力学常数是0.274，壁面的ID是13，我用yk表示污染物的质量分数 #include "udf.h" #define K 0.274 DEFINE_SOURCE(face_source, cell, face, thread, yk, dS, eqn) { real source; fluid_ID = 13; Domain *domain; ?Thread *t = Lookup_Thread(domain,fluid_ID); ?Thread *t0 = THREAD_T0(t); cell_t c0 = F_C0(f,tf); tc0 = THREAD_T0(tf); source=yk*exp(-Kt) dS[eqn]=-K*yk return source } 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
085601 329分调剂
已经有3人回复
-
本2一志愿C9-333分，材料科学与工程，求调剂
已经有3人回复
-
一志愿大连理工大学材料求调剂
已经有5人回复
-
材料与化工调剂一志愿大连海事085600，349
已经有6人回复
-
一志愿中海洋320化学工程与技术学硕求调剂
已经有4人回复
-
材料科学与工程调剂
已经有11人回复
-
311求调剂一志愿合肥工业大学
已经有4人回复
-
一志愿食品科学与工程083200求调剂
已经有4人回复
-
材料科学与工程求调剂
已经有10人回复
-
一志愿北化085600材料专硕275|有文章专利｜求调剂
已经有11人回复
landykun
铁虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 2672.6
- 散金: 55
- 红花: 7
- 帖子: 643
- 在线: 89.9小时
- 虫号: 1772570
- 注册: 2012-04-23
- 性别: GG
- 专业: 流体力学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
1 源项都是单位时间内的源项，是个速率，所以 我觉得不用加，你再研究一下； 2 根据公式 公式怎么样就怎么样，肯定要转换啊，质量分数变成浓度； 3 你画网格的时候，壁面一层单独画成一个zone，如果是用icem画网格 就是单独一个block 如果是gambit或者其他的 我还真不完全清楚，然后在改zone里面 加载source；如果你想在整个计算区域都加载source 然后判断临近边界的区域 我觉得太麻烦了，可能还浪费计算时间。所以，如果网格好画的话，建议在壁面第一层，单独画个zone。如果实在不行 那你就采用判断方式吧。 |
4楼2019-06-03 12:38:21
landykun
铁虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 2672.6
- 散金: 55
- 红花: 7
- 帖子: 643
- 在线: 89.9小时
- 虫号: 1772570
- 注册: 2012-04-23
- 性别: GG
- 专业: 流体力学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
首先 不知道你憋了多久写的，不过写的确实不好 ，建议好好看看fluent udf帮助里面的source部分，很简单的。 1 DEFINE_SOURCE (name, c, t, dS, eqn) 这个宏，你怎么在这个括号里面随便继续加变量？这个变量是不能变的。 2 包括cell 以及thread，通过这个宏fluent已经自动给你传进来了，你还Thread *t = Lookup_Thread(domain,fluid_ID); Thread *t0 = THREAD_T0(t); cell_t c0 = F_C0(f,tf); tc0 = THREAD_T0(tf)，都是多余的。你只需要指定每个cell里面的消耗速率值即可，置于是在哪个cell 哪个thread，这个宏已经给你传进来了，你不用管了。 3你的消耗速率公式是什么？ 4 ds的数值 我写udf的时候 基本不反馈； 5 消耗速率计算计算公式里面还有时间吗，不太清楚。 5 我觉的应该这么写， #include "udf.h" DEFINE_SOURCE (kc_source, c, t, dS, eqn) { real source; real yk,k yk=每个cell的组分质量百分比宏函数，你找找udf帮助，我忘了； k=0.274; source=消耗速率计算公式； return source } |
» 本帖已获得的红花（最新10朵）
2楼2019-06-02 11:31:36
送红花一朵 |
谢谢给予回复，谢谢！ 1.我的反应公式就是一级反应动力学公式，dc/dt=-Kc，所以把它解析解之后，就带了时间，请问我这个公式需要那么变吗？因为看源项案例，都把他们解析解了 2.式中是浓度，但是入口条件只有写质量分数，我直接当他们相等了。可以直接用这个公式吗？还是需要换算一下？ 2.我看到别人说加壁面源项需要先把临近壁面的网格提出来，我加的不是整体的源项，是一个面的，所以加的那些指针，请问如果不加直接在壁面上的udf导入，不知道宏能不能找到壁面临近的cell呢？ 发自小木虫Android客户端 |
3楼2019-06-03 08:10:35
5楼2019-06-04 10:37:54
