public static CalcResult NinMout()
{
//string[] fluid = new string[] { "Water" };
string[] fluid = new string[] { "R410A.MIX" };
//string[] fluid = new string[] { "ISOBUTAN" };
double[] composition = new double[] { 1 };
CalcResult res = new CalcResult();
int Nrow = 2;
double[] FPI = new double[Nrow + 1];
//FPI = new double[] { 1.27, 1.27, 1.27, 1.27, 1.27, 1.27, 1.27, 2.6, 2.6, 2.6, 2.6, 2.6, 2.6, 2.6, 2.6, 5.2, 5.2, 5.2, 5.2, 5.2, 5.2 };
FPI = new double[] { 15, 15 };
double Pt = 1 * 25.4 * 0.001;
double Pr = 0.75 * 25.4 * 0.001;
double Di = 8.4074 * 0.001; //8 6.8944
double Do = 10.0584 * 0.001; //8.4 7.35
double Fthickness = 0.095 * 0.001;
double thickness = 0.5 * (Do - Di);
//double n_tubes = 10;
//double n_rows = 2;
//int[] Ntube = { 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 };
//int[] Ntube = { 2, 2, 2, 2 };
int[] Ntube = { 24, 24 };
int N_tube = Ntube[0];
double L = 914.4 * 0.001;
int Nelement = 1;
int[,] CirArrange;
//CirArrange = new int[,] { { 25, 26, 27, 28, 4, 3, 2, 1, 0, 0 }, { 29, 30, 31, 7, 6, 5, 0, 0, 0, 0 }, { 32, 33, 34, 35, 11, 10, 9, 8, 0, 0 },
//{ 36, 37, 38, 39, 40, 16, 15, 14, 13, 12 }, { 41, 42, 43, 19, 18, 17, 0, 0, 0, 0 }, { 44, 45, 46, 47, 48, 24, 23, 22, 21, 20 } };
//CircuitNumber CircuitInfo = new CircuitNumber();
//CircuitInfo.number = new int[] { 4, 2 }; //4in 2out
//CircuitInfo.TubeofCir = new int[] { 8, 6, 8, 10, 6, 10 }; //{ 4, 8 };
//CircuitInfo.UnequalCir = new int[] { 5, 5, 6, 6, 0, 0 };
//// [19 - 17 - 15 - 13 11 9 7 5 3 1] <====Air
//// [20 - 18 - 16 - 14 12 10 8 6 4 2] <====Air
//5in2out, test3
//CirArrange = new int[,] { { 25, 26, 27, 28, 4, 3, 2, 1, 0, 0 }, { 29, 30, 31, 7, 6, 5, 0, 0, 0, 0 }, { 43, 42, 41, 17, 18, 19, 0, 0, 0, 0 },
//{40, 39, 15, 16, 0, 0, 0, 0, 0, 0 }, { 44, 45, 46, 47, 48, 24, 23, 22, 21, 20 }, { 36, 37, 38, 14, 13, 12, 0, 0, 0, 0 },
//{32, 33, 34, 35, 11, 10, 9, 8, 0, 0 }};
//CirArrange = new int[,] { { 41, 42, 43, 19, 18, 17, 0, 0, 0, 0 }, { 44, 45, 46, 47, 48, 24, 23, 22, 21, 20 }, { 25, 26, 27, 28, 4, 3, 2, 1, 0, 0 },
//{29, 30, 31, 7, 6, 5, 0, 0, 0, 0 }, { 32, 33, 34, 35, 11, 10, 9, 8, 0, 0 }, { 36, 37, 38, 39, 40, 16, 15, 14, 13, 12 }
//};
//Paralle flow
CirArrange = new int[, ] {
{ 17, 18, 19, 43, 42, 41, 0, 0, 0, 0 }, { 20, 21, 22, 23, 24, 48, 47, 46, 45, 44 }, { 1, 2, 3, 4, 28, 27, 26, 25, 0, 0 },
{ 5, 6, 7, 31, 30, 29, 0, 0, 0, 0 }, { 8, 9, 10, 11, 35, 34, 33, 32, 0, 0 }, { 12, 13, 14, 15, 16, 40, 39, 38, 37, 36 }
};
CircuitNumber CircuitInfo = new CircuitNumber();
CircuitInfo.number = new int[] { 4, 2 }; //4in 2out
CircuitInfo.TubeofCir = new int[] { 6, 10, 8, 6, 8, 10 }; //{ 4, 8 };
CircuitInfo.UnequalCir = new int[] { -5, -6, 5, 5, 6, 6 }; //{ 3, 4, -3, -3, -4, -4 };
double mr = 0.076;
double Vel_a = 1.8; //m/s
double H = Pt * N_tube;
double Hx = L * H;
double rho_a_st = 1.2; //kg/m3
double Va = Vel_a * Hx;
double ma = Va * rho_a_st; //Va / 3600 * 1.2; //kg/s
int curve = 1; //
double za = 1; //Adjust factor
double ha = AirHTC.alpha(Vel_a, za, curve) * 1.5; //71.84;//36.44;
double eta_surface = 1;
double zh = 1;
double zdp = 1;
double tai = 26.67;
double tri = 7.2;
double te = tri;
double pe = Refrigerant.SATT(fluid, composition, te + 273.15, 1).Pressure;
double P_exv = 1842.28; //kpa
double T_exv = 24; //C
double conductivity = 386; //w/mK for Cu
double Pwater = 0;
int hexType = 0; //*********************************0 is evap, 1 is cond******************************************
//double densityL = Refrigerant.SATT(fluid, composition, te + 273.15, 1).DensityV;
//double hri = Refrigerant.ENTHAL(fluid, composition, tri + 273.15, densityL).Enthalpy ;
double wm = Refrigerant.WM(fluid, composition).Wm; //g/mol
//hri = hri / wm - 140;
double hri = Refrigerant.TPFLSH(fluid, composition, T_exv + 273.15, P_exv).h / wm - (fluid[0] == "Water" ? 0 : 140);
//double hri = 354.6;
//double xin = 0.57;
double[, ,] ta = new double[Nelement, N_tube, Nrow + 1];
//string AirDirection="DowntoUp";
string AirDirection = "Counter";
ta = InitialAirProperty.AirTemp(Nelement, Ntube, Nrow, tai, te, AirDirection);
GeometryResult geo = new GeometryResult();
//GeometryResult[,] geo_element = new GeometryResult[,] { };
GeometryResult[,] geo_element = new GeometryResult[N_tube, Nrow];
for (int k = 0; k < Nrow; k++)
{
for (int j = 0; j < N_tube; j++)
{
geo_element[j, k] = Areas.Geometry(L / Nelement, FPI[k], Do, Di, Pt, Pr, Fthickness);
//geo_element[i] = Areas.Geometry(L / element, FPI[i], Do, Di, Pt, Pr);
geo.Aa_tube += geo_element[j, k].Aa_tube;
geo.Aa_fin += geo_element[j, k].Aa_fin;
geo.A_a += geo_element[j, k].A_a;
geo.A_r += geo_element[j, k].A_r;
geo.A_r_cs += geo_element[j, k].A_r_cs;
//geo.A_ratio += geo_element[j,k].A_ratio;
}
}
geo.A_ratio = geo.A_r / geo.A_a;
res = Slab.SlabCalc(CirArrange, CircuitInfo, Nrow, Ntube, Nelement, fluid, composition, Di, L, geo_element, ta, te, pe, hri,
mr, ma, ha, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater);
//res = Slab.SlabCalc(Npass, N_tubes_pass, fluid, composition, Dh, L, geo.A_a, geo.A_r_cs, geo.A_r, tai, tri, pe, hri,
// mr, ma, ha, eta_surface, zh, zdp);
//Tsh_calc = res.Tro - (Refrigerant.SATP(fluid, composition, res.Pro, 1).Temperature - 273.15);
// res = Slab.SlabCalc(Npass, N_tubes_pass, fluid, composition, Dh, L, geo.A_a, geo.A_r_cs, geo.A_r, tai, tri, pe, hri,
// mr, ma, ha, eta_surface, zh, zdp);
// Tsh_calc = res.Tro - (Refrigerant.SATP(fluid, composition, res.Pro, 1).Temperature - 273.15);
return(res);
}
public static CalcResult Water_Midea9()
{
string[] fluid = new string[] { "Water" };
double[] composition = new double[] { 1 };
CalcResult res = new CalcResult();
int Nrow = 1;
double[] FPI = new double[Nrow + 1];
FPI = new double[] { 21 };
double Pt = 21 * 0.001;
double Pr = 18.19 * 0.001;
double Di = 6.8944 * 0.001; //8
double Do = 7.35 * 0.001; //8.4
double Fthickness = 0.095 * 0.001;
double thickness = 0.5 * (Do - Di);
int[] Ntube = { 12 };
int N_tube = Ntube[0];
double L = 410 * 0.001;
int Nelement = 1;
int[,] CirArrange;
CirArrange = new int[, ] {
{ 1, 2, 3, 4 }, { 5, 6, 7, 8 }, { 9, 10, 11, 12 }
};
//CirArrange = new int[,] { { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 } };
CircuitNumber CircuitInfo = new CircuitNumber();
CircuitInfo.number = new int[] { 3, 3 };
CircuitInfo.TubeofCir = new int[] { 4, 4, 4 }; //{ 4, 8 };
// [19 - 17 - 15 - 13 11 9 7 5 3 1] <====Air
// [20 - 18 - 16 - 14 12 10 8 6 4 2] <====Air
double mr = 9.99 / 60;
double Vel_a = 2; //m/s
double H = Pt * N_tube;
double Hx = L * H;
double rho_a_st = 1.188; //kg/m3
double Va = Vel_a * Hx;
double ma = Va * rho_a_st; //Va / 3600 * 1.2; //kg/s
int curve = 1; //
double za = 1; //Adjust factor
double zh = 1;
double zdp = 1;
double eta_surface = 0.8284;
double ha = AirHTC.alpha(Vel_a, za, curve) / 79 * 77.42;//71.84;//36.44;
double tai = 19.98;
double tri = 44.98;
double tc = tri;
double pc = Refrigerant.SATT(fluid, composition, tc + 273.15, 1).Pressure;
double Pwater = 395; //kpa
double conductivity = 386; //w/mK for Cu
int hexType = 1; //*********************************0 is evap, 1 is cond******************************************
double wm = Refrigerant.WM(fluid, composition).Wm; //g/mol
double hri = Refrigerant.TPFLSH(fluid, composition, tc + 273.15, Pwater).h / wm - 0.5 - (fluid[0] == "Water" ? 0 : 140);
double[, ,] ta = new double[Nelement, N_tube, Nrow + 1];
//string AirDirection="DowntoUp";
string AirDirection = "Counter";
ta = InitialAirProperty.AirTemp(Nelement, Ntube, Nrow, tai, tc, AirDirection);
GeometryResult geo = new GeometryResult();
//GeometryResult[,] geo_element = new GeometryResult[,] { };
GeometryResult[,] geo_element = new GeometryResult[N_tube, Nrow];
for (int k = 0; k < Nrow; k++)
{
for (int j = 0; j < N_tube; j++)
{
geo_element[j, k] = Areas.Geometry(L / Nelement, FPI[k], Do, Di, Pt, Pr, Fthickness);
geo.Aa_tube += geo_element[j, k].Aa_tube;
geo.Aa_fin += geo_element[j, k].Aa_fin;
geo.A_a += geo_element[j, k].A_a;
geo.A_r += geo_element[j, k].A_r;
geo.A_r_cs += geo_element[j, k].A_r_cs;
//geo.A_ratio += geo_element[j,k].A_ratio;
}
}
geo.A_ratio = geo.A_r / geo.A_a;
res = Slab.SlabCalc(CirArrange, CircuitInfo, Nrow, Ntube, Nelement, fluid, composition, Di, L, geo_element, ta, tc, pc, hri,
mr, ma, ha, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater);
return(res);
}
public static CalcResult Water_Heat_Midea9()
{
string[] fluid = new string[] { "Water" };
double[] composition = new double[] { 1 };
CalcResult res = new CalcResult();
int Nrow = 1;
double[] FPI = new double[Nrow + 1];
FPI = new double[] { 21 };
double Pt = 21 * 0.001;
double Pr = 18.19 * 0.001;
double Di = 6.8944 * 0.001; //8
double Do = 7.35 * 0.001; //8.4
double Fthickness = 0.095 * 0.001;
double thickness = 0.5 * (Do - Di);
int[] Ntube = { 12 };
int N_tube = Ntube[0];
double L = 410 * 0.001;
int Nelement = 1;
int[,] CirArrange;
CirArrange = new int[, ] {
{ 1, 2, 3, 4 }, { 5, 6, 7, 8 }, { 9, 10, 11, 12 }
};
CircuitNumber CircuitInfo = new CircuitNumber();
CircuitInfo.number = new int[] { 3, 3 };
CircuitInfo.TubeofCir = new int[] { 4, 4, 4 }; //{ 4, 8 };
int N = 28;
for (int i = 0; i < N; i++)
{
double[] mr = new double[] { 10.01, 14.01, 18.0, 21.0, 25.0, 29.01, 31.0, 9.99, 14.01, 18.01, 21.01, 25.0, 29.01, 31.0, 10.01, 14.01, 18.01, 21.01, 25.0, 29.01, 31.01, 9.99, 14.0, 18.0, 21.0, 24.99, 29.0, 31.01 }; // 60;
mr[i] = mr[i] / 60;
double[] Vel_a = new double[] { 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0 }; //m/s
double H = Pt * N_tube;
double Hx = L * H;
double rho_a_st = 1.188; //kg/m3
double[] Va = new double[N];
Va[i] = Vel_a[i] * Hx;
double[] ma = new double[N];
ma[i] = Va[i] * rho_a_st; //Va / 3600 * 1.2; //kg/s
int curve = 1; //
double za = 1; //Adjust factor
double zh = 1;
double zdp = 1;
double[] eta_surface = new double[] { 0.8764, 0.8764, 0.8764, 0.8764, 0.8764, 0.8764, 0.8764, 0.865, 0.865, 0.865, 0.865, 0.865, 0.865, 0.865, 0.855, 0.855, 0.855, 0.855, 0.855, 0.855, 0.855, 0.846, 0.846, 0.846, 0.846, 0.846, 0.846, 0.846 };
//double ha = AirHTC.alpha(Vel_a, za, curve) / 79 * 78.7;//71.84;//36.44;
//double[] ha = new double[] { 58.92, 58.92, 58.92, 58.92, 58.92, 58.92, 58.92, 65.61, 65.61, 65.61, 65.61, 65.61, 65.61, 65.61, 71.19, 71.19, 71.19, 71.19, 71.19, 71.19, 71.19, 75.95, 75.95, 75.95, 75.95, 75.95, 75.95, 75.95 };
double[] ha = new double[] { 58.51, 58.51, 58.51, 58.51, 58.51, 58.51, 58.51, 64.92, 64.92, 64.92, 64.92, 64.92, 64.92, 64.92, 70.41, 70.41, 70.41, 70.41, 70.41, 70.41, 70.41, 75.42, 75.42, 75.42, 75.42, 75.42, 75.42, 75.42 };
double[] tai = new double[] { 20.01, 20.0, 20.02, 20.0, 19.99, 20.02, 20.0, 19.98, 20.0, 19.99, 19.99, 19.99, 19.98, 20.0, 19.98, 20.01, 19.99, 20.0, 20.0, 19.99, 19.99, 19.98, 20.01, 20.0, 20.0, 20.01, 20.01, 20.02 };
double[] tri = new double[] { 45.01, 45.01, 44.99, 44.99, 45.0, 45.0, 44.99, 44.99, 45.01, 44.99, 45.01, 44.99, 45.01, 45.01, 45.0, 45.01, 44.99, 45.01, 45.01, 44.99, 45.01, 44.98, 45.01, 44.99, 45.0, 45.0, 45.01, 45.0 };
double[] tc = tri;
double[] pc = new double[N];
pc[i] = Refrigerant.SATT(fluid, composition, tc[i] + 273.15, 1).Pressure;
double Pwater = 395; //kpa
double conductivity = 386; //w/mK for Cu
int hexType = 1; //*********************************0 is evap, 1 is cond******************************************
double wm = Refrigerant.WM(fluid, composition).Wm; //g/mol
double[] hri = new double[N];
hri[i] = Refrigerant.TPFLSH(fluid, composition, tc[i] + 273.15, Pwater).h / wm - 0.5 - (fluid[0] == "Water" ? 0 : 140);
double[, ,] ta = new double[Nelement, N_tube, Nrow + 1];
string AirDirection = "Counter";
ta = InitialAirProperty.AirTemp(Nelement, Ntube, Nrow, tai[i], tc[i], AirDirection);
GeometryResult geo = new GeometryResult();
GeometryResult[,] geo_element = new GeometryResult[N_tube, Nrow];
for (int k = 0; k < Nrow; k++)
{
for (int j = 0; j < N_tube; j++)
{
geo_element[j, k] = Areas.Geometry(L / Nelement, FPI[k], Do, Di, Pt, Pr, Fthickness);
//geo_element[i] = Areas.Geometry(L / element, FPI[i], Do, Di, Pt, Pr);
geo.Aa_tube += geo_element[j, k].Aa_tube;
geo.Aa_fin += geo_element[j, k].Aa_fin;
geo.A_a += geo_element[j, k].A_a;
geo.A_r += geo_element[j, k].A_r;
geo.A_r_cs += geo_element[j, k].A_r_cs;
//geo.A_ratio += geo_element[j,k].A_ratio;
}
}
geo.A_ratio = geo.A_r / geo.A_a;
res = Slab.SlabCalc(CirArrange, CircuitInfo, Nrow, Ntube, Nelement, fluid, composition, Di, L, geo_element, ta, tc[i], pc[i], hri[i],
mr[i], ma[i], ha[i], eta_surface[i], zh, zdp, hexType, thickness, conductivity, Pwater);
//using (StreamWriter wr = File.AppendText(@"D:\Work\Simulation\Test\Midea9_heat.txt"))
//{
// wr.WriteLine("Q, {0}, DP, {1}, href, {2}, Ra_ratio, {3}, Tao, {4}, Tro, {5}", res.Q, res.DP, res.href, res.Ra_ratio, res.Tao, res.Tro);
//}
}
return(res);
}
public static CalcResult Water_Midea5()
{
string[] fluid = new string[] { "Water" };
//string[] fluid = new string[] { "R410A.MIX" };
//string[] fluid = new string[] { "ISOBUTAN" };
double[] composition = new double[] { 1 };
CalcResult res = new CalcResult();
int Nrow = 3;
double[] FPI = new double[Nrow + 1];
//FPI = new double[] { 1.27, 1.27, 1.27, 1.27, 1.27, 1.27, 1.27, 2.6, 2.6, 2.6, 2.6, 2.6, 2.6, 2.6, 2.6, 5.2, 5.2, 5.2, 5.2, 5.2, 5.2 };
FPI = new double[] { 17, 17, 17 };
double Pt = 21 * 0.001;
double Pr = 13.37 * 0.001;
double Di = 6.8944 * 0.001; //8
double Do = 7.35 * 0.001; //8.4
double Fthickness = 0.095 * 0.001;
double thickness = 0.5 * (Do - Di);
//double n_tubes = 10;
//double n_rows = 2;
//int[] Ntube = { 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 };
//int[] Ntube = { 2, 2, 2, 2 };
int[] Ntube = { 14, 14, 14 };
int N_tube = Ntube[0];
double L = 367 * 0.001;
int Nelement = 1;
int[,] CirArrange;
//CirArrange = new int[,] { { 8, 6, 4, 2, 1, 3, 5, 7 } };//actual, counter-paralle, Q=83.1
//CirArrange = new int[,] { { 20, 18, 16, 14, 12, 10, 8, 6, 4, 2, 1, 3, 5, 7, 9, 11, 13, 15, 17, 19 } };//actual, counter-paralle, Q=83.1
//CirArrange = new int[,] { { 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 1, 3, 5, 7, 9, 11, 13, 15, 17, 19 } }; //paralle-paralle, better Q=85.3
//CirArrange = new int[,] { { 20, 18, 16, 14, 12, 10, 8, 6, 4, 2, 19, 17, 15, 13, 11, 9, 7, 5, 3, 1 } };//counter-counter, Q=82.4
//CirArrange = new int[,] { { 14, 12, 10, 8, 6, 4, 2, 1, 3, 5, 7, 9, 11, 13 } };//actual, counter-paralle, Q=76
//CirArrange = new int[,] { {16, 14, 12, 10, 8, 6, 4, 2, 1, 3, 5, 7, 9, 11, 13, 15 } };//actual, counter-paralle, Q=79
//CirArrange = new int[,] { {42, 40, 38, 36, 34, 32, 30, 28, 26, 24, 22, 20, 18, 16, 14, 12, 10, 8, 6, 4, 2,
// 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33, 35, 37, 39, 41 } };//actual, counter-paralle, Q=79
//CirArrange = new int[,] { {32, 30, 28, 26, 24, 22, 20, 18, 16, 14, 12, 10, 8, 6, 4, 2,
// 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31 } };//actual, counter-paralle, Q=79
//CirArrange = new int[,] { { 7, 8, 2, 1, 5, 3 }, {5, 8, 9, 10, 11}, {5, 88, 5, 4, 3 } };
//CirArrange = new int[,] { { 7, 8, 2, 1, 0, 0, 0, 0 }, { 9, 10, 11, 12, 6, 5, 4, 3 } };
//CirArrange = new int[,] { { 7, 1, 0, 0, 0, 0 }, { 8, 9, 3, 2, 0, 0 }, { 10, 11, 12, 6, 5, 4 } };
//CirArrange = new int[,] { { 7, 8, 2, 1, 9, 10, 11, 12, 6, 5, 4, 3 } };
//List<string> productType = new List<string>();
//CirArrange = new int[,] { { 7, 8, 9, 3, 2, 1 }, { 12, 11, 10, 4, 5, 6 } };
//CirArrange = new int[,] { { 7, 8, 9, 4, 5, 6 }, { 12, 11, 10, 3, 2, 1 } };
//CirArrange = new int[,] { { 7, 8, 5, 6 }, { 10, 9, 2, 1 }, { 12, 11, 4, 3 } };
//CirArrange = new int[,] { { 7, 8, 2, 1 }, { 10, 9, 3, 4 }, { 11, 12, 6, 5 } };
CirArrange = new int[, ] {
{ 29, 30, 16, 15, 1, 2 }, { 32, 31, 17, 18, 4, 3 }, { 33, 34, 20, 19, 5, 6 }, { 36, 35, 21, 22, 8, 7 }, { 37, 38, 24, 23, 9, 10 }, { 40, 39, 25, 26, 12, 11 }, { 41, 42, 28, 27, 13, 14 }
};
CircuitNumber CircuitInfo = new CircuitNumber();
CircuitInfo.number = new int[] { 7, 7 };
CircuitInfo.TubeofCir = new int[] { 6, 6, 6, 6, 6, 6, 6 }; //{ 4, 8 };
// [19 - 17 - 15 - 13 11 9 7 5 3 1] <====Air
// [20 - 18 - 16 - 14 12 10 8 6 4 2] <====Air
// Ncir=1, 20in, 20->19 1out
//int total = 0;
//if (CirArrange != null)
//{
// foreach (var seg in CirArrange)
// {
// total += seg.
// }
//}
double mr = 23.0 / 60;
double Vel_a = 1.2; //m/s
double H = Pt * N_tube;
double Hx = L * H;
double rho_a_st = 1.2; //kg/m3
double Va = Vel_a * Hx;
double ma = Va * rho_a_st; //Va / 3600 * 1.2; //kg/s
int curve = 1; //
double za = 1; //Adjust factor
double zh = 1;
double zdp = 1;
double eta_surface = 0.89;
double ha = AirHTC.alpha(Vel_a, za, curve);//71.84;//36.44;
double tai = 20;
double tri = 45;
double tc = tri;
double pc = Refrigerant.SATT(fluid, composition, tc + 273.15, 1).Pressure;
double Pwater = 305; //kpa
double conductivity = 386; //w/mK for Cu
int hexType = 1; //*********************************0 is evap, 1 is cond******************************************
//double densityL = Refrigerant.SATT(fluid, composition, te + 273.15, 1).DensityV;
//double hri = Refrigerant.ENTHAL(fluid, composition, tri + 273.15, densityL).Enthalpy ;
double wm = Refrigerant.WM(fluid, composition).Wm; //g/mol
//hri = hri / wm - 140;
double hri = Refrigerant.TPFLSH(fluid, composition, tc + 273.15, Pwater).h / wm - 0.5 - (fluid[0] == "Water" ? 0 : 140);
//double hri = 354.6;
//double xin = 0.57;
double[, ,] ta = new double[Nelement, N_tube, Nrow + 1];
//string AirDirection="DowntoUp";
string AirDirection = "Counter";
ta = InitialAirProperty.AirTemp(Nelement, Ntube, Nrow, tai, tc, AirDirection);
GeometryResult geo = new GeometryResult();
//GeometryResult[,] geo_element = new GeometryResult[,] { };
GeometryResult[,] geo_element = new GeometryResult[N_tube, Nrow];
for (int k = 0; k < Nrow; k++)
{
for (int j = 0; j < N_tube; j++)
{
geo_element[j, k] = Areas.Geometry(L / Nelement, FPI[k], Do, Di, Pt, Pr, Fthickness);
//geo_element[i] = Areas.Geometry(L / element, FPI[i], Do, Di, Pt, Pr);
geo.Aa_tube += geo_element[j, k].Aa_tube;
geo.Aa_fin += geo_element[j, k].Aa_fin;
geo.A_a += geo_element[j, k].A_a;
geo.A_r += geo_element[j, k].A_r;
geo.A_r_cs += geo_element[j, k].A_r_cs;
//geo.A_ratio += geo_element[j,k].A_ratio;
}
}
geo.A_ratio = geo.A_r / geo.A_a;
res = Slab.SlabCalc(CirArrange, CircuitInfo, Nrow, Ntube, Nelement, fluid, composition, Di, L, geo_element, ta, tc, pc, hri,
mr, ma, ha, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater);
//res = Slab.SlabCalc(Npass, N_tubes_pass, fluid, composition, Dh, L, geo.A_a, geo.A_r_cs, geo.A_r, tai, tri, pe, hri,
// mr, ma, ha, eta_surface, zh, zdp);
//Tsh_calc = res.Tro - (Refrigerant.SATP(fluid, composition, res.Pro, 1).Temperature - 273.15);
// res = Slab.SlabCalc(Npass, N_tubes_pass, fluid, composition, Dh, L, geo.A_a, geo.A_r_cs, geo.A_r, tai, tri, pe, hri,
// mr, ma, ha, eta_surface, zh, zdp);
// Tsh_calc = res.Tro - (Refrigerant.SATP(fluid, composition, res.Pro, 1).Temperature - 273.15);
return(res);
}