24小时热门版块排行榜    

查看: 1660  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856调剂 +8 刘梦微 2026-02-28 8/400 2026-03-02 19:23 by caszguilin
[考研] 290分材料工程085601求调剂 数二英一 +4 llx0610 2026-03-02 4/200 2026-03-02 19:19 by caszguilin
[考研] 求调剂 +10 yunziaaaaa 2026-03-01 11/550 2026-03-02 19:17 by caszguilin
[考研] 085600 英一数二272求调剂 5+6 vida_a 2026-03-01 16/800 2026-03-02 19:13 by zhukairuo
[考研] 化工京区271求调剂 +6 11ing 2026-03-02 6/300 2026-03-02 18:52 by caszguilin
[考研] 环境调剂 +3 chenhanheng 2026-03-02 3/150 2026-03-02 17:37 by houyaoxu
[考研] 一志愿华中科技大学,化学专业344分,求调剂 +3 邢xing1 2026-03-02 3/150 2026-03-02 17:32 by houyaoxu
[考研] 一志愿东北大学材料专硕328,求调剂 +3 shs1083 2026-03-02 3/150 2026-03-02 17:27 by houyaoxu
[考研] 0703 总分319求调剂 +3 Xinyuu 2026-03-02 3/150 2026-03-02 17:20 by houyaoxu
[考研] 考研复试调剂,过国家线的同学都可报名 +4 黑!在干嘛 2026-02-28 5/250 2026-03-02 16:48 by chock1337
[考研] 化工专硕348,一志愿985求调剂 +6 弗格个 2026-02-28 9/450 2026-03-02 14:09 by liyongv
[考博] 诚招农业博士 +3 心欣向荣 2026-02-28 3/150 2026-03-02 13:33 by 时间不狗
[考研] 272求调剂 +7 材紫有化 2026-02-28 7/350 2026-03-02 12:48 by 无际的草原
[考研] 26考研报考西工大材料308分求调剂 +4 weizhong123 2026-03-01 4/200 2026-03-02 12:46 by 无际的草原
[考研] 295求调剂 +8 19171856320 2026-02-28 8/400 2026-03-02 11:19 by yuchj
[考研] 0856材料与化工,270求调剂 +8 YXCT 2026-03-01 9/450 2026-03-02 11:01 by 无际的草原
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 307求调剂 +5 wyyyqx 2026-03-01 5/250 2026-03-01 15:21 by Fff-1
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
信息提示
请填处理意见