24小时热门版块排行榜    

查看: 2114  |  回复: 8

jinjiaoshou

新虫 (初入文坛)

[求助] 求助,多孔介质燃烧的udf编辑问题 已有2人参与

控制方程如上图,我现在有双能量方程的udf,但是燃烧中包含了组分变化的项,请问怎么用udf写入啊
DEFINE_UDS_UNSTEADY(myunsteady, c, t, i, apu, su)
{
real physical_dt, vol, phi_old;
physical_dt = RP_Get_Real("physical-time-step";
vol=C_VOLUME(c,t);
*apu=-vol/physical_dt;         
phi_old=C_STORAGE_R(c,t,SV_UDSI_M1(i));
*su=vol*phi_old/physical_dt;         
}

求助,多孔介质燃烧的udf编辑问题
控制方程.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinjiaoshou

新虫 (初入文坛)

→_→

[ 发自小木虫客户端 ]
2楼2015-05-18 16:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liukan12

银虫 (小有名气)

【答案】应助回帖

#include "udf.h"
DEFINE_PROFILE(velocity1,t,i)
{
    real x[ND_ND];    /* 存放2D或者3D位置参数 */
    real y;          /* y值为模型中的y坐标值 */
    cell_t f;        /* 创建网格体指针 */
    begin_c_loop(f,t)  /* 体线程t中所有网格体指针f的遍历循环 */
      {
       C_CENTROID(x,f,t);  /* 把体线程t中所有指针f中的网格体的中心点的坐标值赋给数组x[ND_ND] */
       y = x[1];          /* 把x[ND_ND]中的第二个元素集合赋给变量y,x[0]则是x坐标,x[2]是z坐标,根据你的需要 */
       if ( y > 0.5 )      /* 条件判断,当y值大于0.5时,执行下列语句 */
           {
                   C_PROFILE(f, t, i) = y^2;   /* y^2是表示入口速度等于y坐标值的平方,这里可以改成你要的关系式 */
       }
       else                              /* 条件判断 */
           {
                   C_PROFILE(f, t, i) = 1;      /* 条件判断,当y值小于0.5时,执行语句 */
       }
          
      }
    end_c_loop(f,t)          /* 遍历循环结束 */
}

你把上面的C换成F的话 就是从遍历体循环到面循环 如果是二维的话 基本上都是F
然后你还要通过interpret 跟fluent hook一下 然后再在velocity入口条件下 导入这个udf
ANSYSFluent流体仿真研究人员。
3楼2015-05-18 17:44:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liukan12

银虫 (小有名气)

不好意思 ,不小心回复错人了
ANSYSFluent流体仿真研究人员。
4楼2015-05-18 17:55:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinjiaoshou

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by liukan12 at 2015-05-18 17:55:14
不好意思 ,不小心回复错人了

大神你知道全局常数怎么写么,就是两段常数不一样的,我写的是这样的提示错误
#include "udf.h"
#define dp
real x[ND_ND];
real y;
cell_t c;
Thread *t;
C_CENTROID(x,c,t);
y=x[1];
if(y>0.06)
dp=0.06;
else
dp=0.03;
但是我把上面这段放在Define_property里,当作一个局部变量也会出问题,不知道为什么
5楼2015-05-18 18:08:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

5428liu

新虫 (小有名气)

引用回帖:
5楼: Originally posted by jinjiaoshou at 2015-05-18 18:08:35
大神你知道全局常数怎么写么,就是两段常数不一样的,我写的是这样的提示错误
#include "udf.h"
#define dp
real x;
real y;
cell_t c;
Thread *t;
C_CENTROID(x,c,t);
y=x;
if(y>0.06)
dp ...

楼主你好,我用的控制方程跟你用的一模一样,我在添加弥散项的时候(方程2),需要各组分的密度和比热,会计算出很多结果,请问你是怎么表示的呢,是需要矩阵吗?
6楼2015-09-18 14:12:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

5428liu

新虫 (小有名气)

【答案】应助回帖

请教楼主扩散速度编写udf需要什么宏呢?急求急求
7楼2015-09-18 14:39:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

need帮助

新虫 (初入文坛)

师兄你好,能分享一下你的双能量方程吗,有偿
来学习模拟的知识
8楼2020-06-22 15:59:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笑已泛黄

新虫 (小有名气)

引用回帖:
7楼: Originally posted by 5428liu at 2015-09-18 14:39:09
请教楼主扩散速度编写udf需要什么宏呢?急求急求

define_diffusity

发自小木虫Android客户端
9楼2020-06-22 20:35:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jinjiaoshou 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 18:59 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:50 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-15 15:50 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 14:38 by VMDqgnlKimtZ
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 14:14 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
信息提示
请填处理意见