UDF 如何调用气体分压？

weselyas

我的反应器里面的总压是固定的2.0MPa，其中有诸多的气体组份：CO，CO2，H20，CH4等，我想调用CO2的分压（CO2分压力随反应的进行而变化的），怎么调用的啊？非常感谢给位大侠！

gearboy

原帖由 weselyas 于 2011-11-28 07:52 发表

                               
登录/注册后可看大图

我的反应器里面的总压是固定的2.0MPa，其中有诸多的气体组份：CO，CO2，H20，CH4等，我想调用CO2的分压（CO2分压力随反应的进行而变化的），怎么调用的啊？非常感谢给位大侠！

                 int i;
        real co2_mole_fract,P_co2,P_h2o,P_total;
        real Mw[N_GAS],total_mole;
        Domain*domain;
        Material *mix_mat,*sp;
               
                    domain=Get_Domain(1);
                mix_mat=mixture_material(domain);
                total_mole=0;
                mixture_species_loop(mix_mat,sp,i)
                {
                        Mw=MATERIAL_PROP(sp,PROP_mwi);  //molecular weight of species i
                        total_mole+=C_YI(c,t,i)/Mw;
                }
                co2_mole_fract=(C_YI(c,t,Species_Panel_Num[CO2_NUM])/Mw[Species_Panel_Num[CO2_NUM]])/total_mole;
                P_total=ABS_P(C_P(c,t),op_pres);  // absolute total pressure,Pa
                P_co2=P_total*co2_mole_fract;   //co2 partial pressure,Pa