24小时热门版块排行榜    

查看: 1794  |  回复: 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的回帖

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的回帖
查看全部 9 个回答

jinjiaoshou

新虫 (初入文坛)

→_→

[ 发自小木虫客户端 ]
2楼2015-05-18 16:40:05
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见