24小时热门版块排行榜    

查看: 1641  |  回复: 0

hu9924

金虫 (小有名气)

[求助] 一个关于fluent自定义修改密度的udf

在模拟计算气体燃烧,单项,需要使用udf修改密度项,挂上后发现密度值不变,不知道错在哪里,大家一起讨论下。
关于这部分的udf代码如下:
DEFINE_PROPERTY(cell_density,c,t)
{
        real rho;
        real t_c1;
        real t_c2;
        real row_ding,row_iding1,row_iding2;/*density and intermediate variable */
        int j,j_ding,i,i_ding;
/*下面部分是根据两个自定义的标量值来查询每一点的密度值,这部分查询计算程序没有问题*/
        if (C_UDSI(c,t,0)==1)

        {
                row_ding=row[0][0];
       
        }
        else if (C_UDSI(c,t,0)==2)
        {
       
        row_ding=row[100][0];
        }
        else
        {


        t_c1=C_UDSI(c,t,0);
        t_c2=C_UDSI(c,t,1)/(C_UDSI(c,t,0)*(1-C_UDSI(c,t,0)));
                for (j=0;j<=100;j++)
            {

                if ((t_c2>=c2[j])&&(t_c2<c2[j+1]))
                {
                                j_ding=j;
                }
                else if (t_c2=c2[100])
                {
                                j_ding=100;
                }
                       
        }
                for (i=0;i<=100;i++)
                {

                        if ((t_c1=c1)&&(t_c1<c1[i+1]))
                        {
                                i_ding=i;
                       
                        }

                        else if (t_c1=c1[100])
                        {
                                i_ding=100;
                               
                        }
                       
                       
                }
       
               
                if((i_ding<100)&&(j_ding<100))
                    {
                               
                                row_iding1=(t_c2-c2[j_ding])*(row[i_ding][j_ding+1]-row[i_ding][j_ding])/(c2[j_ding+1]-c2[j_ding])+row[i_ding][j_ding];
                                row_iding2=(t_c2-c2[j_ding])*(row[i_ding+1][j_ding+1]-row[i_ding+1][j_ding])/(c2[j_ding+1]-c2[j_ding])+row[i_ding+1][j_ding];
                                row_ding=(t_c1-c1[i_ding])*(row_iding2-row_iding1)/(c1[i_ding+1]-c1[i_ding])+row_iding1;
               
                     }
               
               
                else if((i_ding=100)&&(j_ding<100))
                    {
                                row_ding=(t_c2-c2[j_ding])*(row[i_ding][j_ding+1]-row[i_ding][j_ding])/(c2[j_ding+1]-c2[j_ding])+row[i_ding][j_ding];
              
                     }
                else if ((i_ding<100)&&(j_ding=100))

                    {
                                row_ding=(t_c1-c1[i_ding])*(row[i_ding+1][j_ding]-row[i_ding][j_ding])/(c1[i_ding+1]-c1[i_ding])+row[i_ding][j_ding];
                    }
                else if ((i_ding=100)&&(j_ding=100))
                     {
           row_ding=row[i_ding][j_ding];
                     }
        }
rho=row_ding;
return rho;
}
回复此楼

» 猜你喜欢

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

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 hu9924 的主题更新
信息提示
请填处理意见