Exemplo n.º 1
0
        public static AirCoef_res alpha1(double vel, double za, int curve, GeometryInput geoInput, int hexType)
        {
            AirCoef_res res = new AirCoef_res();
            //空气侧换热系数计算
            double Pt         = geoInput.Pt;
            double Pr         = geoInput.Pr;
            double Do         = geoInput.Do;
            double Fthickness = geoInput.Fthickness;
            double P_f        = 25.4 / geoInput.FPI / 1000;
            int    Nrow       = geoInput.Nrow;

            double W_hex = Nrow * Pr;

            double rho_a_st = 1.2;                                        //kg/m3
            double cp_a     = (hexType == 0 ? 1.027 * 1000: 1.02 * 1000); //J/kg/K
            double k_a      = 0.025;                                      //W/K/m
            double mu_a     = 18.2 * Math.Pow(10, -6);                    //Pa*s

            double De_c = 4 * (P_f - Fthickness) * (Pt * Pr - Math.PI * Math.Pow(Do + 2 * Fthickness, 2.0) / 4) / (2 * (Pt * Pr - Math.PI * Math.Pow(Do + 2 * Fthickness, 2.0) / 4) + Math.PI * (Do + 2 * Fthickness) * (P_f - Fthickness));

            double V_f       = vel * (P_f * Pt * Pr) / ((P_f - Fthickness) * (Pt * Pr - Math.PI * Math.Pow(Do + 2 * Fthickness, 2.0) / 4));
            double Re_a      = De_c * V_f * rho_a_st / mu_a;
            double Pr_a      = mu_a * cp_a / k_a;
            double Nusselt_a = 2.1 * Math.Pow(Re_a * Pr_a * De_c / W_hex, 0.38);
            double ha        = Nusselt_a * k_a / De_c * za;

            //空气侧压降计算
            double f    = (0.43 + 35.1 * Math.Pow(Re_a * De_c / W_hex, -1.07)) * De_c / W_hex;
            double dP_a = 2 * f * W_hex * rho_a_st * Math.Pow(V_f, 2) / De_c;

            //ha = za * ha;

            return(res = new AirCoef_res {
                ha = ha, dP_a = dP_a
            });
        }
Exemplo n.º 2
0
 public static double[,] alpha_cal(double[,] ha, double[,] Va, double Va_ave, double Vel_ave, double za, int curve, GeometryInput geoInput_air, int hexType, int N_tube, int Nelement)
 {
     for (int i = 0; i < N_tube; i++)
     {
         for (int j = 0; j < Nelement; j++)
         {
             if (i == 0 && j == 0)
             {
                 ha[i, j] = AirHTC.alpha1(Va[i, j] * (Vel_ave / Va_ave), za, curve, geoInput_air, hexType).ha;
             }
             else
             {
                 for (int ii = 0; ii <= i; ii++)
                 {
                     for (int jj = 0; jj < Nelement; jj++)
                     {
                         //if (VaDistri.Va[i, j] == VaDistri.Va[ii, jj])
                         if (Math.Abs((Va[i, j] - Va[ii, jj]) / Va[i, j]) < 0.001)
                         {
                             ha[i, j] = ha[ii, jj];
                             break;
                         }
                     }
                 }
                 if (ha[i, j] == 0)
                 {
                     ha[i, j] = AirHTC.alpha1(Va[i, j] * (Vel_ave / Va_ave), za, curve, geoInput_air, hexType).ha;
                 }
             }
         }
     }
     return(ha);
 }