24小时热门版块排行榜    

查看: 1665  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料292调剂 +3 橘颂思美人 2026-03-23 3/150 2026-03-23 19:36 by pswait
[考研] 考研化学308分求调剂 +7 你好明天你好 2026-03-23 8/400 2026-03-23 18:39 by macy2011
[考研] 一志愿211 初试270分 求调剂 +3 谷雨上岸 2026-03-23 4/200 2026-03-23 15:53 by ACS Nano——
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +5 vv迷 2026-03-22 5/250 2026-03-23 07:36 by Iveryant
[考研] 311求调剂 +6 冬十三 2026-03-18 6/300 2026-03-22 20:18 by edmund7
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
信息提示
请填处理意见