public static CalcResult SlabCalc(int[,] CirArrange, CircuitNumber CircuitInfo, int Nrow, int[] Ntube, int Nelement, string fluid, //double Npass, int[] N_tubes_pass,
double l, Geometry geo, double[, ,] ta, double[, ,] RH,
double te, double pe, double hri, double mr, double[,] ma, double[,] ha, double[,] haw,
double eta_surface, double zh, double zdp, int hexType, double thickness, double conductivity, double Pwater, string Airdirection, CapiliaryInput cap_inlet, CapiliaryInput cap_outlet, AbstractState coolprop, double[,] SourceTableData)
{
//------->
// R2 R1
// [11 1] <====
// [12 2] <====
// <==== Air
// [13 3] <====
// [14 4] <====
// [15 5] <====
// [16 6] <====
// [17 7] <====
// [18 8] <====
// [19 9] <====
// [20 10] <====
// Ncir=1, 11in, 20->10 1out
// [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
// CirArrange = new int[,] { { 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1 } };
// Nrow=2
// Ncir=2
//AbstractState coolprop = AbstractState.factory("HEOS", fluid);
double[] Q1 = new double[50];
double tri = te; //Refrigerant.SATP(fluid, composition, pri, 1).Temperature - 273.15;
double pri = pe;
int Nciri = CircuitInfo.number[0];
int Nciro = CircuitInfo.number[1];
int Ncir = (Nciri == Nciro ? Nciri : Nciri + Nciro);
int N_tube = Ntube[0];
int N_tube_total = 0;
int iRow = 0;
int iTube_o = 0;
int iTube_n = 0;
int index_o = 0;
int index_n = 0;
double te_calc_org = 0;
CirArr[] cirArr = new CirArr[Nrow * N_tube];
CirArrforAir cirArrforAir = new CirArrforAir();
cirArrforAir = CirArrangement.ReadCirArr(CirArrange, CircuitInfo, Nrow, Ntube, 0);
cirArr = cirArrforAir.CirArr;
CalcResult res_slab = new CalcResult();
double[] pri_cir = new double[Ncir]; //[element, tube, row]
double[] hri_cir = new double[Ncir];
double[] tri_cir = new double[Ncir];
double[] mr_ciri = new double[Nciri];
List <double[]> mr_ciri_base = new List <double[]>();
double[] mr_ciro = new double[Nciro];
int[] Ngroupin = new int[Nciro];
int index = 0;
int restartDP_index = 0;
int N_tube2 = 0;
int[] index_cir = new int[Ncir];
int index_mr_ciri_base = 0;
CalcResult[] r = new CalcResult[Ncir];
CalcResult[] r1 = new CalcResult[Ncir];
CalcResult[] r2 = new CalcResult[Ncir]; //for NinMout only
CalcResult[] res_cir2 = new CalcResult[Nciro + 1];
CalcResult[] res_type = new CalcResult[Nciri + 1];
double[,] Q_detail = new double[N_tube, Nrow];//detail output
double[,] DP_detail = new double[N_tube, Nrow];
double[,] Tri_detail = new double[N_tube, Nrow];
double[,] Pri_detail = new double[N_tube, Nrow];
double[,] hri_detail = new double[N_tube, Nrow];
double[,] Tro_detail = new double[N_tube, Nrow];
double[,] Pro_detail = new double[N_tube, Nrow];
double[,] hro_detail = new double[N_tube, Nrow];
double[,] href_detail = new double[N_tube, Nrow];
double[,] mr_detail = new double[N_tube, Nrow];
double[,] charge_detail = new double[N_tube, Nrow];
int flag_ciro = 0;
int Ncir_forDP = 0;
double[] mr_forDP = new double[Nciri];
int k;
double te_calc = 0;
CheckAir airConverge = new CheckAir();
CheckDP dPconverge = new CheckDP();
CheckPri priconverge = new CheckPri();
for (int i = 0; i < Nrow; i++)
{
N_tube_total += Ntube[i];
}
for (int i = 0; i < Nciro; i++)
{
mr_ciro[i] = mr / Nciro;
}
bool index_outbig;
if (CircuitInfo.UnequalCir == null || CircuitInfo.UnequalCir[0] > 0)
{
index_outbig = false;
}
else
{
index_outbig = true;
}
if (CircuitInfo.UnequalCir != null)
{
for (int j = 0; j < Nciro; j++)
{
for (int i = 0; i < Ncir; i++)
{
if (CircuitInfo.UnequalCir[i] == Nciri + 1 + j)
{
Ngroupin[j]++;
}
}
//for (int i = 0; i < Nciri; i++) mr_ciri[i] = mr_ciro[j] / Ngroupin[j];
}
}
int iterforAir = 0;
int iterforDP = 0;
int iterforPri = 0;
double pri_1 = 0;
double pri_2 = 0;
double pri_3 = 0;
double pri_4 = 0;
double tri1 = 0;
double te1 = 0;
double iterforDP_ciri = 0;
double iterforDP_ciro = 0;
//Starting properties
iterforDP_ciro = 0;
double[] mr_forDP_o_4 = new double[Nciri];
double[] mr_forDP_o_3 = new double[Nciri];
double[] mr_forDP_o_2 = new double[Nciri];
double[] mr_forDP_o_1 = new double[Nciri];
do
{
#region //AirConverge
do
{
//iterforDP = 0;
r = new CalcResult[Ncir];
r1 = new CalcResult[Ncir];
res_cir2 = new CalcResult[Nciro + 1];
flag_ciro = (index_outbig ? 1 : 0);
//tri = tri;
//制冷制热模块计算切换
if (hexType == 0)
{
coolprop.update(input_pairs.PQ_INPUTS, pri * 1000, 0);
tri = coolprop.T() - 273.15;
tri = CoolProp.PropsSI("T", "P", pri * 1000, "Q", 0, fluid) - 273.15;
int a = 1;
}
else
{
coolprop.update(input_pairs.PQ_INPUTS, pri * 1000, 0);
te = coolprop.T() - 273.15;
te = CoolProp.PropsSI("T", "P", pri * 1000, "Q", 0, fluid) - 273.15;
int a = 1;
}
for (int j = 0; j < (flag_ciro == 1 ? (index_outbig ? Nciri + 1 : 1) : Nciro + 1); j++)
{
if (j >= Nciro)
{
j = j - Nciro; //for Nciro
flag_ciro = (index_outbig ? 0 : 1);
}
if (j == 1 && index_outbig && index == 0)
{
j = j - 1;
flag_ciro = 0;
}
index = 0;
iterforDP_ciri = 0;
double[] mr_forDP_4 = new double[Nciri];
double[] mr_forDP_3 = new double[Nciri];
double[] mr_forDP_2 = new double[Nciri];
double[] mr_forDP_1 = new double[Nciri];
double [] mr0_forDP = new double[5];
#region //DPConverge
do
{
res_type = new CalcResult[Nciri + 1];
k = 0;
if (!index_outbig)
{
for (int i = 0; i < (flag_ciro == 1 ? Nciro : (Nciri == Nciro ? Ncir : Ncir - Nciro)); i++)
{
if (flag_ciro == 1)
{
pri_cir[i + Ncir - Nciro] = res_cir2[i].Pro;
hri_cir[i + Ncir - Nciro] = res_cir2[i].hro;
tri_cir[i + Ncir - Nciro] = res_cir2[i].Tro;
}
else
{
pri_cir[i] = pri;
hri_cir[i] = hri;
tri_cir[i] = tri;
}
}
}
else
{
for (int i = 0; i < (flag_ciro == 1 ? Nciro : Ncir); i++)
{
if (flag_ciro == 1)
{
pri_cir[i] = pri;
hri_cir[i] = hri;
tri_cir[i] = tri;
}
else
{
if (CircuitInfo.UnequalCir[i] == Nciri + 1 + j)
{
pri_cir[i] = r2[j].Pro;
hri_cir[i] = r2[j].hro;
tri_cir[i] = r2[j].Tro;
}
}
}
}
for (int i = 0; i < Ncir; i++)
{
if (flag_ciro == 1)
{
//汇管计算
if (CircuitInfo.UnequalCir[i] <= 0)
{
//for (int i = 0; i < Ncir; i++)
//r[i] = Circuit.CircuitCalc(i, cirArr, CircuitInfo, Nrow, Ntube, Nelement, fluid, dh, l, geo.element, ta, RH,
// tri_cir[i], pri_cir[i], hri_cir[i], mr_ciro[k], ma, ha, haw, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater, Airdirection, d_cap, lenth_cap);
r[i] = Circuit.CircuitCalc(i, cirArr, CircuitInfo, Nrow, Ntube, Nelement, fluid, l, geo, ta, RH,
tri_cir[i], pri_cir[i], hri_cir[i], mr_ciro[k], ma, ha, haw, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater, cap_inlet, cap_outlet, coolprop, SourceTableData);
if (r[i].Pro < 0)
{
res_slab.Pro = -10000000; return(res_slab);
}
r1[k] = r[i].ShallowCopy();
r2[k] = r[i].ShallowCopy();
if (!index_outbig)
{
r1[k].DP += res_cir2[k].DP;
}
index_cir[k] = i;
k++;
Ncir_forDP = Nciro;
mr_forDP = (double[])mr_ciro.Clone(); // mr_forDP = mr_ciro
if (k == Nciro)
{
break;
}
}
}
else if (Nciri == Nciro || CircuitInfo.UnequalCir[i] == Nciri + 1 + j)
{
//均匀流路计算和不均匀流路开始部分（独立管）计算
if (index == 0)
{
if (Nciri == Nciro && iterforPri == 0)
{
mr_ciro.CopyTo(mr_ciri, 0);
}
else if (Nciri != Nciro)
{
if (restartDP_index == 1 || !priconverge.flag)
{
var mm = mr_ciri_base[j].Sum();
//foreach (var item in mr_ciri_base[j])
//mm += item;
mr_ciri[k] = mr_ciri_base[j][k] * mr_ciro[j] / mm;//(mr / Nciro);
}
else
{
mr_ciri[k] = mr_ciro[j] / Ngroupin[j];
}
}
}
//else mr_ciri_base.CopyTo(mr_ciri[k], 0);
//for (int i = 0; i < Ncir; i++)
//首次流路计算
if (CircuitInfo.CirType != null && CircuitInfo.CirType.flag == true)
{
bool circuit_cap = false;
//if (d_cap[index] == 0 && lenth_cap[index] == 0) circuit_cap = true;
if ((CircuitInfo.CirType.flag == true) && (CircuitInfo.CirType.type[i, 0] == 0) && (res_type[CircuitInfo.CirType.type[i, 1]] != null && circuit_cap))
{
//r[i] = res_type[CircuitInfo.CirType.type[i, 1]];
r[i] = new CalcResult();
r[i].DP = res_type[CircuitInfo.CirType.type[i, 1]].DP;
r[i].DP_cap = res_type[CircuitInfo.CirType.type[i, 1]].DP_cap;
r[i].href = res_type[CircuitInfo.CirType.type[i, 1]].href;
r[i].hro = res_type[CircuitInfo.CirType.type[i, 1]].hro;
r[i].M = res_type[CircuitInfo.CirType.type[i, 1]].M;
r[i].Pro = res_type[CircuitInfo.CirType.type[i, 1]].Pro;
r[i].Q = res_type[CircuitInfo.CirType.type[i, 1]].Q;
r[i].R_1 = res_type[CircuitInfo.CirType.type[i, 1]].R_1;
r[i].R_1a = res_type[CircuitInfo.CirType.type[i, 1]].R_1a;
r[i].R_1r = res_type[CircuitInfo.CirType.type[i, 1]].R_1r;
r[i].Ra_ratio = res_type[CircuitInfo.CirType.type[i, 1]].Ra_ratio;
r[i].RHout = res_type[CircuitInfo.CirType.type[i, 1]].RHout;
r[i].Tao = res_type[CircuitInfo.CirType.type[i, 1]].Tao;
r[i].Tri = res_type[CircuitInfo.CirType.type[i, 1]].Tri;
r[i].Tro = res_type[CircuitInfo.CirType.type[i, 1]].Tro;
r[i].x_i = res_type[CircuitInfo.CirType.type[i, 1]].x_i;
r[i].x_o = res_type[CircuitInfo.CirType.type[i, 1]].x_o;
r[i].Vel_r = res_type[CircuitInfo.CirType.type[i, 1]].Vel_r;
r[i].mr = res_type[CircuitInfo.CirType.type[i, 1]].mr;
r[i].Q_detail = new double[N_tube, Nrow];//
r[i].DP_detail = new double[N_tube, Nrow];
r[i].Tri_detail = new double[N_tube, Nrow];
r[i].Pri_detail = new double[N_tube, Nrow];
r[i].hri_detail = new double[N_tube, Nrow];
r[i].Tro_detail = new double[N_tube, Nrow];
r[i].Pro_detail = new double[N_tube, Nrow];
r[i].hro_detail = new double[N_tube, Nrow];
r[i].href_detail = new double[N_tube, Nrow];
r[i].mr_detail = new double[N_tube, Nrow];
r[i].Tao_Detail = new double[Nelement, N_tube, Nrow];
r[i].RHo_Detail = new double[Nelement, N_tube, Nrow];
r[i].charge_detail = new double[N_tube, Nrow];
for (int m = 0; m < CircuitInfo.TubeofCir[i]; m++)
{
index_o = 0;
index_n = 0;
if (i == 0)
{
index_n = 0;
}
else
{
for (int n = 1; n <= i; n++)
{
index_n += CircuitInfo.TubeofCir[n - 1];
}
}
if (res_type[CircuitInfo.CirType.type[i, 1]].index == 0)
{
index_o = 0;
}
else
{
for (int n = 1; n <= res_type[CircuitInfo.CirType.type[i, 1]].index; n++)
{
index_o += CircuitInfo.TubeofCir[n - 1];
}
}
iRow = cirArr[m + index_o].iRow;
iTube_o = cirArr[m + index_o].iTube;
iTube_n = cirArr[m + index_n].iTube;
r[i].Q_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].Q_detail[iTube_o, iRow];//
r[i].DP_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].DP_detail[iTube_o, iRow];
r[i].Tri_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].Tri_detail[iTube_o, iRow];
r[i].Pri_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].Pri_detail[iTube_o, iRow];
r[i].hri_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].hri_detail[iTube_o, iRow];
r[i].Tro_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].Tro_detail[iTube_o, iRow];
r[i].Pro_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].Pro_detail[iTube_o, iRow];
r[i].hro_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].hro_detail[iTube_o, iRow];
r[i].href_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].href_detail[iTube_o, iRow];
r[i].mr_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].mr_detail[iTube_o, iRow];
r[i].charge_detail[iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].charge_detail[iTube_o, iRow];
for (int p = 0; p < Nelement; p++)
{
//ta[p, iTube_n, iRow + 1] = res_type[CircuitInfo.CirType.type[i, 1]].Tao_Detail[p, iTube_o, iRow];
//RH[p, iTube_n, iRow + 1] = res_type[CircuitInfo.CirType.type[i, 1]].RHo_Detail[p, iTube_o, iRow];
r[i].Tao_Detail[p, iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].Tao_Detail[p, iTube_o, iRow];
r[i].RHo_Detail[p, iTube_n, iRow] = res_type[CircuitInfo.CirType.type[i, 1]].RHo_Detail[p, iTube_o, iRow];
}
}
//r[i].Tao_Detail = ta;
//r[i].RHo_Detail = RH;
}
else
{
//r[i] = Circuit.CircuitCalc(i, cirArr, CircuitInfo, Nrow, Ntube, Nelement, fluid, dh, l, geo.element, ta, RH,
//tri_cir[i], pri_cir[i], hri_cir[i], mr_ciri[k], ma, ha, haw, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater, Airdirection, d_cap, lenth_cap);
r[i] = Circuit.CircuitCalc(i, cirArr, CircuitInfo, Nrow, Ntube, Nelement, fluid, l, geo, ta, RH,
tri_cir[i], pri_cir[i], hri_cir[i], mr_ciri[k], ma, ha, haw, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater, cap_inlet, cap_outlet, coolprop, SourceTableData);
if (r[i].Pro < 0)
{
res_slab.Pro = -10000000; return(res_slab);
}
if (CircuitInfo.CirType.type[i, 0] == 0)
{
res_type[CircuitInfo.CirType.type[i, 1]] = r[i];
res_type[CircuitInfo.CirType.type[i, 1]].index = i;
}
}
}
else
{
//r[i] = Circuit.CircuitCalc(i, cirArr, CircuitInfo, Nrow, Ntube, Nelement, fluid, dh, l, geo.element, ta, RH,
//tri_cir[i], pri_cir[i], hri_cir[i], mr_ciri[k], ma, ha, haw, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater, Airdirection, d_cap, lenth_cap);
r[i] = Circuit.CircuitCalc(i, cirArr, CircuitInfo, Nrow, Ntube, Nelement, fluid, l, geo, ta, RH,
tri_cir[i], pri_cir[i], hri_cir[i], mr_ciri[k], ma, ha, haw, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater, cap_inlet, cap_outlet, coolprop, SourceTableData);
if (r[i].Pro < 0)
{
res_slab.Pro = -10000000; return(res_slab);
}
}
r1[k] = r[i].ShallowCopy();
index_cir[k] = i;//不均匀流路的输出才会用到
k++;
if (k == (Nciri == Nciro ? Ncir : Ngroupin[j]))
{
Ncir_forDP = (Nciri == Nciro ? Ncir : Ngroupin[j]);
mr_forDP = (double[])mr_ciri.Clone();
break;
}
}
}
if (index_outbig && flag_ciro == 1)
{
break;
}
index++;
//dPconverge = CheckDPforCircuits.CheckDPConverge(mr, mr_ciri, r, Ncir);
dPconverge = CheckDPforCircuits.CheckDPConverge(hexType, res_cir2, iterforPri, flag_ciro, mr_forDP, r1, Ncir_forDP);
if (flag_ciro == 0)
{
iterforDP_ciri++;
if (iterforDP_ciri >= 5)
{
mr_forDP_4 = mr_forDP_3;
mr_forDP_3 = mr_forDP_2;
mr_forDP_2 = mr_forDP_1;
mr_forDP_1 = mr_forDP;
try
{
if (mr_forDP_1[0] < mr_forDP_2[0] && (Math.Abs(mr_forDP_1[0] - mr_forDP_3[0]) / mr_forDP_1[0] < 0.0001) ||
Math.Abs(mr_forDP_1[0] - mr_forDP_4[0]) / mr_forDP_1[0] < 0.0001)
{
dPconverge.flag = true;
}
}
catch (Exception ex)
{
continue;
}
}
restartDP_index = 0;
if (!dPconverge.flag)
{
dPconverge.mr.CopyTo(mr_ciri, 0); //mr_ciri = dPconverge.mr;
}
}
else //(flag_ciro == 1)
{
iterforDP_ciro++;
if (iterforDP_ciro >= 5)
{
mr_forDP_o_4 = mr_forDP_o_3;
mr_forDP_o_3 = mr_forDP_o_2;
mr_forDP_o_2 = mr_forDP_o_1;
mr_forDP_o_1 = mr_forDP;
try
{
if (mr_forDP_o_1[0] < mr_forDP_o_2[0] && (Math.Abs(mr_forDP_o_1[0] - mr_forDP_o_3[0]) / mr_forDP_o_1[0] < 0.0001) ||
Math.Abs(mr_forDP_o_1[0] - mr_forDP_o_4[0]) / mr_forDP_o_1[0] < 0.0001)
{
dPconverge.flag = true;
}
}
catch (Exception ex)
{
continue;
}
}
if (dPconverge.flag)
{
restartDP_index = 0;
}
else
{
restartDP_index = 1;
dPconverge.mr.CopyTo(mr_ciro, 0); //mr_ciro = dPconverge.mr;
break;
}
}
iterforDP++;
N_tube2 = 0;
#region //Result print out
if (dPconverge.flag)
{
if (Nciri == Nciro)
{
//te_calc = Refrigerant.SATP(fluid, composition, r[j].Pro, 1).Temperature;
coolprop.update(input_pairs.PQ_INPUTS, r[j].Pro * 1000, 0);
te_calc = coolprop.T();
te_calc = CoolProp.PropsSI("T", "P", r[j].Pro * 1000, "Q", 0, fluid);
int a = 1;
}
else
{
if (mr_ciri_base.Count == Nciro && flag_ciro == 0)
{
mr_ciri_base.RemoveAt(index_mr_ciri_base);
mr_ciri_base.Insert(index_mr_ciri_base, mr_forDP);
index_mr_ciri_base++;
index_mr_ciri_base %= mr_ciri_base.Count;
}
if (mr_ciri_base.Count < Nciro)
{
mr_ciri_base.Add(mr_forDP); //keep original mr ratio for fast iter
}
j = (flag_ciro == 1 ? j + Nciro : j);
res_cir2[j] = new CalcResult();
for (int i = 0; i < (flag_ciro == 1 ? Nciro : Ngroupin[j]); i++)
{
res_cir2[j].Q += r1[i].Q;
res_cir2[j].M += r1[i].M;
res_cir2[j].hro += (flag_ciro == 1 ? mr_ciro[i] : mr_ciri[i]) * r1[i].hro;
if (fluid == "Water")
{
res_cir2[j].Tro += (flag_ciro == 1 ? mr_ciro[i] : mr_ciri[i]) * r1[i].Tro;
}
res_cir2[j].Vel_r = r1[i].Vel_r;
res_cir2[j].href += r1[i].href * CircuitInfo.TubeofCir[index_cir[i]];
res_cir2[j].R_1 += r1[i].R_1 * CircuitInfo.TubeofCir[index_cir[i]];
res_cir2[j].R_1a += r1[i].R_1a * CircuitInfo.TubeofCir[index_cir[i]];
res_cir2[j].R_1r += r1[i].R_1r * CircuitInfo.TubeofCir[index_cir[i]];
N_tube2 += CircuitInfo.TubeofCir[index_cir[i]];
}
res_cir2[j].DP = r1[(flag_ciro == 1 ? Nciro : Ngroupin[j]) - 1].DP;
res_cir2[j].Tao_Detail = ta;
res_cir2[j].Pro = r1[(flag_ciro == 1 ? Nciro : Ngroupin[j]) - 1].Pro;
res_cir2[j].hro = res_cir2[j].hro / (flag_ciro == 1 ? mr : mr_ciro[j]);
res_cir2[j].href = res_cir2[j].href / N_tube2;
res_cir2[j].R_1 = res_cir2[j].R_1 / N_tube2;
res_cir2[j].R_1a = res_cir2[j].R_1a / N_tube2;
res_cir2[j].R_1r = res_cir2[j].R_1r / N_tube2;
res_cir2[j].Tri = tri;
//te_calc = Refrigerant.SATP(fluid, composition, res_cir2[j].Pro, 1).Temperature;
coolprop.update(input_pairs.PQ_INPUTS, res_cir2[j].Pro * 1000, 0);
te_calc = coolprop.T();
te_calc = CoolProp.PropsSI("T", "P", res_cir2[j].Pro * 1000, "Q", 0, fluid);
if (fluid == "Water")
{
res_cir2[j].Tro = res_cir2[j].Tro / (flag_ciro == 1 ? mr : mr_ciro[j]) - 273.15;
}
else
{
res_cir2[j].Tro = CoolProp.PropsSI("T", "P", res_cir2[j].Pro * 1000, "H", res_cir2[j].hro * 1000, fluid) - 273.15;
}
}
}
#endregion
#endregion
} while (!dPconverge.flag && iterforDP < 200);
if (Nciri == Nciro)
{
break;
}
if (index_outbig && (j == Nciro - 1) && (res_cir2[0] != null))
{
for (int i = 0; i < Nciro; i++)
{
r2[i].DP += res_cir2[i].DP;
}
flag_ciro = 1;
Ncir_forDP = Nciro;
mr_forDP = (double[])mr_ciro.Clone(); // mr_forDP = mr_ciro
dPconverge = CheckDPforCircuits.CheckDPConverge(hexType, res_cir2, iterforPri, flag_ciro, mr_forDP, r2, Ncir_forDP);
iterforDP_ciro++;
if (iterforDP_ciro >= 5)
{
mr_forDP_o_4 = mr_forDP_o_3;
mr_forDP_o_3 = mr_forDP_o_2;
mr_forDP_o_2 = mr_forDP_o_1;
mr_forDP_o_1 = mr_forDP;
try
{
if (mr_forDP_o_1[0] < mr_forDP_o_2[0] && (Math.Abs(mr_forDP_o_1[0] - mr_forDP_o_3[0]) / mr_forDP_o_1[0] < 0.0001) ||
Math.Abs(mr_forDP_o_1[0] - mr_forDP_o_4[0]) / mr_forDP_o_1[0] < 0.0001)
{
dPconverge.flag = true;
}
}
catch (Exception ex)
{
continue;
}
}
if (!dPconverge.flag)
{
restartDP_index = 1;
dPconverge.mr.CopyTo(mr_ciro, 0); //mr_ciro = dPconverge.mr;
}
break;
}
}
if (Airdirection == "顺流")
{
airConverge.flag = true;
for (int ii = 0; ii < Ncir; ii++)
{
for (int i = 0; i < Nrow; i++)
{
for (int j = 0; j < N_tube; j++)
{
for (int kk = 0; kk < Nelement; kk++)
{
if (r[ii].Tao_Detail[kk, j, i] != 0)
{
ta[kk, j, i + 1] = r[ii].Tao_Detail[kk, j, i];
}
if (r[ii].RHo_Detail[kk, j, i] != 0)
{
RH[kk, j, i + 1] = r[ii].RHo_Detail[kk, j, i];
}
}
}
}
}
}
else//Counter
{
airConverge = CheckAirConvergeforCircuits.CheckAirConverge(cirArrforAir.TotalDirection, Nrow, N_tube, Nelement, ta, RH, r); //taout_calc, RHout_calc
ta = airConverge.ta;
RH = airConverge.RH;
iterforAir++;
}
//Add Q converge criterion to avoid results oscillation, ruhao 20180426
if (Airdirection != "顺流") //No airConverge iter for Parallel
{
for (int i = 0; i < Ncir; i++)
{
Q1[iterforAir - 1] += r[i].Q;
}
try
{
if (Q1[iterforAir - 1] < Q1[iterforAir - 2] && (Math.Abs(Q1[iterforAir - 1] - Q1[iterforAir - 3]) / Q1[iterforAir - 1] < 0.0001) ||
Math.Abs(Q1[iterforAir - 1] - Q1[iterforAir - 4]) / Q1[iterforAir - 1] < 0.0001)
{
airConverge.flag = true;
}
}
catch (Exception ex)
{
continue;
}
}
} while (!airConverge.flag && iterforAir < 50);
#endregion
//using (StreamWriter wr = File.AppendText(@"D:\Work\Simulation\Test\MinNout.txt"))
//{
//for (int i = 0; i < Ncir; i++)
//{
//wr.WriteLine("Q, {0}, DP, {1}, href, {2}, Ra_ratio, {3}, Tao, {4}, Tro, {5}, mr, {6}", r[i].Q, r[i].DP, r[i].href, r[i].Ra_ratio, r[i].Tao, r[i].Tro, r[i].mr);
//}
//}
if (restartDP_index == 1)
{
priconverge.flag = false;
}
else if (hexType == 0 && (fluid != "Water"))
{
priconverge = CheckPin.CheckPriConverge(te, te_calc - 273.15, te_calc_org - 273.15, pri, pe, r[Ncir - 1].Pro); //res_slab.Pro
iterforPri++;
if (iterforPri >= 20)
{
pri_4 = pri_3;
pri_3 = pri_2;
pri_2 = pri_1;
pri_1 = pri;
try
{
if (pri_1 < pri_2 && (Math.Abs(pri_1 - pri_3) / pri_1 < 1e-5) ||
Math.Abs(pri_1 - pri_4) / pri_1 < 1e-5)
{
priconverge.flag = true;
}
}
catch (Exception ex)
{
continue;
}
}
pri = priconverge.pri;
te_calc_org = te_calc;
if (priconverge.flag && iterforPri == 1 && iterforDP == 1)
{
priconverge.flag = false; //to avoid not even iterate but converge by chance
}
}
else
{
priconverge.flag = true;
}
} while (!priconverge.flag && iterforPri < 100);
// if (iterforDP >= 200)
//{
// return res_slab;
// throw new Exception("iter for DPConverge > 100.");
//}
//if (iterforPri >= 50)
//{
// return res_slab;
// throw new Exception("iter for Pri > 50.");
//}
#region //Result print out
for (int i = 0; i < Ncir; i++)
{
res_slab.Q += r[i].Q;
res_slab.M += r[i].M;
if (Nciri == Nciro)
{
res_slab.hro += mr_ciri[i] * r[i].hro;
}
res_slab.href += r[i].href * CircuitInfo.TubeofCir[i];
res_slab.R_1 += r[i].R_1 * CircuitInfo.TubeofCir[i];
res_slab.R_1a += r[i].R_1a * CircuitInfo.TubeofCir[i];
res_slab.R_1r += r[i].R_1r * CircuitInfo.TubeofCir[i];
for (int j = 0; j < N_tube; j++)//detail output
{
for (k = 0; k < Nrow; k++)
{
if (r[i].Q_detail[j, k] != 0)
{
Q_detail[j, k] = r[i].Q_detail[j, k];
}
if (r[i].DP_detail[j, k] != 0)
{
DP_detail[j, k] = r[i].DP_detail[j, k];
}
if (r[i].Tri_detail[j, k] != 0)
{
Tri_detail[j, k] = r[i].Tri_detail[j, k];
}
if (r[i].Pri_detail[j, k] != 0)
{
Pri_detail[j, k] = r[i].Pri_detail[j, k];
}
if (r[i].hri_detail[j, k] != 0)
{
hri_detail[j, k] = r[i].hri_detail[j, k];
}
if (r[i].Tro_detail[j, k] != 0)
{
Tro_detail[j, k] = r[i].Tro_detail[j, k];
}
if (r[i].Pro_detail[j, k] != 0)
{
Pro_detail[j, k] = r[i].Pro_detail[j, k]; //
}
if (r[i].hro_detail[j, k] != 0)
{
hro_detail[j, k] = r[i].hro_detail[j, k];
}
if (r[i].href_detail[j, k] != 0)
{
href_detail[j, k] = r[i].href_detail[j, k];
}
if (r[i].mr_detail[j, k] != 0)
{
mr_detail[j, k] = r[i].mr_detail[j, k];
}
if (r[i].charge_detail[j, k] != 0)
{
charge_detail[j, k] = r[i].charge_detail[j, k];
}
}
}
}
if (Nciri == Nciro)
{
res_slab.hro = res_slab.hro / mr;
res_slab.Pro = r[Ncir - 1].Pro;
res_slab.DP_cap = r[Ncir - 1].DP_cap;
res_slab.Vel_r = r[Ncir - 1].Vel_r;
}
else if (!index_outbig)
{
res_slab.hro = res_cir2[Nciro].hro;
res_slab.Pro = res_cir2[Nciro].Pro;
res_slab.DP_cap = res_cir2[Nciro].DP_cap;
res_slab.Vel_r = res_cir2[Nciro].Vel_r;
}
else
{
res_slab.hro = res_cir2[Nciro - 1].hro;
res_slab.Pro = res_cir2[Nciro - 1].Pro;
res_slab.DP_cap = res_cir2[Nciro - 1].DP_cap;
res_slab.Vel_r = res_cir2[Nciro - 1].Vel_r;
}
res_slab.Pri = pri;
res_slab.Tri = tri;
res_slab.hri = hri;
res_slab.mr = mr;
res_slab.DP = pri - res_slab.Pro;
res_slab.Tao_Detail = ta;
res_slab.RHo_Detail = RH;
res_slab.href = res_slab.href / N_tube_total;
for (int i = 0; i < ha.GetLength(0); i++)
{
for (int j = 0; j < ha.GetLength(1); j++)
{
res_slab.ha += ha[i, j];
}
}
res_slab.ha = res_slab.ha / ha.Length;
res_slab.R_1 = res_slab.R_1 / N_tube_total;
res_slab.R_1a = res_slab.R_1a / N_tube_total;
res_slab.R_1r = res_slab.R_1r / N_tube_total;
coolprop.update(input_pairs.PQ_INPUTS, res_slab.Pro * 1000, 0);
te_calc = coolprop.T() - 273.15;
te_calc = CoolProp.PropsSI("T", "P", res_slab.Pro * 1000, "Q", 0, fluid) - 273.15;
coolprop.update(input_pairs.QT_INPUTS, 0, te_calc + 273.15);
double densityLo = coolprop.rhomass();
//double densityLo = CoolProp.PropsSI("D", "T", te_calc + 273.15, "Q", 0, fluid);
coolprop.update(input_pairs.QT_INPUTS, 1, te_calc + 273.15);
double densityVo = coolprop.rhomass();
//double densityVo = CoolProp.PropsSI("D", "T", te_calc + 273.15, "Q", 1, fluid);
coolprop.update(input_pairs.DmassT_INPUTS, densityLo, te_calc + 273.15);
double hlo = coolprop.hmass() / 1000;
//double hlo = CoolProp.PropsSI("H", "T", te_calc + 273.15, "D", densityLo, fluid) / 1000 ;
coolprop.update(input_pairs.DmassT_INPUTS, densityVo, te_calc + 273.15);
double hvo = coolprop.hmass() / 1000;
//double hvo = CoolProp.PropsSI("H", "T", te_calc + 273.15, "D", densityVo, fluid) / 1000 ;
res_slab.x_o = (res_slab.hro - hlo) / (hvo - hlo);
coolprop.update(input_pairs.PQ_INPUTS, pri * 1000, 0);
double hli = coolprop.hmass() / 1000;
//double hli = CoolProp.PropsSI("H", "P", pri * 1000, "Q", 0, fluid) / 1000 ;
coolprop.update(input_pairs.PQ_INPUTS, pri * 1000, 1);
double hvi = coolprop.hmass() / 1000;
//double hvi = CoolProp.PropsSI("H", "P", pri * 1000, "Q", 1, fluid) / 1000 ;
res_slab.x_i = (res_slab.hri - hli) / (hvi - hli);
coolprop.update(input_pairs.HmassP_INPUTS, res_slab.hro * 1000, res_slab.Pro * 1000);
res_slab.Tro = coolprop.T() - 273.15;
//res_slab.Tro = CoolProp.PropsSI("T", "P", res_slab.Pro * 1000, "H", res_slab.hro * 1000, fluid) - 273.15;
double h = res_slab.hro;
for (int j = 0; j < N_tube; j++)
{
for (int i = 0; i < Nelement; i++)
{
res_slab.Tao += res_slab.Tao_Detail[i, j, Nrow];
res_slab.RHout += res_slab.RHo_Detail[i, j, Nrow];
}
}
res_slab.Tao = res_slab.Tao / (N_tube * Nelement);
res_slab.RHout = res_slab.RHout / (N_tube * Nelement);
res_slab.Ra_ratio = res_slab.R_1a / res_slab.R_1;
for (int i = 0; i < ma.GetLength(0); i++)
{
for (int j = 0; j < ma.GetLength(1); j++)
{
res_slab.ma += ma[i, j];
}
}
res_slab.Va = res_slab.ma / 1.2 * 3600;
res_slab.Q_detail = Q_detail;//detail output
res_slab.DP_detail = DP_detail;
res_slab.Tri_detail = Tri_detail;
res_slab.Pri_detail = Pri_detail;
res_slab.hri_detail = hri_detail;
res_slab.Tro_detail = Tro_detail;
res_slab.Pro_detail = Pro_detail;
res_slab.hro_detail = hro_detail;
res_slab.href_detail = href_detail;
res_slab.mr_detail = mr_detail;
res_slab.Aa = geo.TotalArea.A_a;
res_slab.Ar = geo.TotalArea.A_r;
res_slab.AHx = geo.TotalArea.A_hx;
res_slab.N_row = Nrow;
res_slab.tube_row = N_tube;
res_slab.charge_detail = charge_detail;
return(res_slab);
#endregion
}
public static CalcResult SlabCalc(int[,] CirArrange, CircuitNumber CircuitInfo, int Nrow, int[] Ntube, int Nelement, string[] fluid, double[] composition, //double Npass, int[] N_tubes_pass,
double dh, double l, GeometryResult[,] geo, double[, ,] ta,
double te, double pe, double hri, double mr, double ma, double ha,
double eta_surface, double zh, double zdp, int hexType, double thickness, double conductivity, double Pwater)
{
//------->
// R2 R1
// [11 1] <====
// [12 2] <====
// <==== Air
// [13 3] <====
// [14 4] <====
// [15 5] <====
// [16 6] <====
// [17 7] <====
// [18 8] <====
// [19 9] <====
// [20 10] <====
// Ncir=1, 11in, 20->10 1out
// [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
// CirArrange = new int[,] { { 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1 } };
// Nrow=2
// Ncir=2
double wm = Refrigerant.WM(fluid, composition).Wm;
double tri = te; //Refrigerant.SATP(fluid, composition, pri, 1).Temperature - 273.15;
double pri = pe;
int Nciri = CircuitInfo.number[0];
int Nciro = CircuitInfo.number[1];
int Ncir = (Nciri == Nciro ? Nciri : Nciri + Nciro);
int N_tube = Ntube[0];
int N_tube_total = 0;
CirArr[] cirArr = new CirArr[Nrow * N_tube];
cirArr = CirArrangement.ReadCirArr(CirArrange, CircuitInfo, Nrow, Ntube);
CalcResult res_slab = new CalcResult();
double[] pri_cir = new double[Ncir]; //[element, tube, row]
double[] hri_cir = new double[Ncir];
double[] tri_cir = new double[Ncir];
double[] mr_ciri = new double[Nciri];
List <double[]> mr_ciri_base = new List <double[]>();
double[] mr_ciro = new double[Nciro];
int[] Ngroupin = new int[Nciro];
int index = 0;
int restartDP_index = 0;
int N_tube2 = 0;
int[] index_cir = new int[Ncir];
CalcResult[] r = new CalcResult[Ncir];
CalcResult[] r1 = new CalcResult[Ncir];
CalcResult[] r2 = new CalcResult[Ncir]; //for NinMout only
CalcResult[] res_cir2 = new CalcResult[Nciro + 1];
int flag_ciro = 0;
int Ncir_forDP = 0;
double[] mr_forDP = new double[Nciri];
int k;
double te_calc = 0;
CheckDP dPconverge = new CheckDP();
CheckPri priconverge = new CheckPri();
for (int i = 0; i < Nrow; i++)
{
N_tube_total += Ntube[i];
}
for (int i = 0; i < Nciro; i++)
{
mr_ciro[i] = mr / Nciro;
}
bool index_outbig;
if (CircuitInfo.UnequalCir == null || CircuitInfo.UnequalCir[0] > 0)
{
index_outbig = false;
}
else
{
index_outbig = true;
}
if (CircuitInfo.UnequalCir != null)
{
for (int j = 0; j < Nciro; j++)
{
for (int i = 0; i < Ncir; i++)
{
if (CircuitInfo.UnequalCir[i] == Nciri + 1 + j)
{
Ngroupin[j]++;
}
}
//for (int i = 0; i < Nciri; i++) mr_ciri[i] = mr_ciro[j] / Ngroupin[j];
}
}
int iterforDP = 0;
int iterforPri = 0;
//Starting properties
do
{
r = new CalcResult[Ncir];
r1 = new CalcResult[Ncir];
res_cir2 = new CalcResult[Nciro + 1];
flag_ciro = (index_outbig ? 1 : 0);
tri = Refrigerant.SATP(fluid, composition, pri, 1).Temperature - 273.15;
for (int j = 0; j < (flag_ciro == 1 ? (index_outbig ? Nciri + 1 : 1) : Nciro + 1); j++)
{
if (j >= Nciro)
{
j = j - Nciro; //for Nciro
flag_ciro = (index_outbig ? 0 : 1);
}
if (j == 1 && index_outbig && index == 0)
{
j = j - 1;
flag_ciro = 0;
}
index = 0;
do
{
k = 0;
if (!index_outbig)
{
for (int i = 0; i < (flag_ciro == 1 ? Nciro : (Nciri == Nciro ? Ncir : Ncir - Nciro)); i++)
{
if (flag_ciro == 1)
{
pri_cir[i + Ncir - Nciro] = res_cir2[i].Pro;
hri_cir[i + Ncir - Nciro] = res_cir2[i].hro;
tri_cir[i + Ncir - Nciro] = res_cir2[i].Tro;
}
else
{
pri_cir[i] = pri;
hri_cir[i] = hri;
tri_cir[i] = tri;
}
}
}
else
{
for (int i = 0; i < (flag_ciro == 1 ? Nciro : Ncir); i++)
{
if (flag_ciro == 1)
{
pri_cir[i] = pri;
hri_cir[i] = hri;
tri_cir[i] = tri;
}
else
{
if (CircuitInfo.UnequalCir[i] == Nciri + 1 + j)
{
pri_cir[i] = r2[j].Pro;
hri_cir[i] = r2[j].hro;
tri_cir[i] = r2[j].Tro;
}
}
}
}
for (int i = 0; i < Ncir; i++)
{
if (flag_ciro == 1)
{
if (CircuitInfo.UnequalCir[i] <= 0)
{
//for (int i = 0; i < Ncir; i++)
r[i] = Circuit.CircuitCalc(i, cirArr, CircuitInfo, Nrow, Ntube, Nelement, fluid, composition, dh, l, geo, ta,
tri_cir[i], pri_cir[i], hri_cir[i], mr_ciro[k], ma, ha, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater);
r1[k] = r[i].ShallowCopy();
r2[k] = r[i].ShallowCopy();
if (!index_outbig)
{
r1[k].DP += res_cir2[k].DP;
}
index_cir[k] = i;
k++;
Ncir_forDP = Nciro;
mr_forDP = (double[])mr_ciro.Clone(); // mr_forDP = mr_ciro
if (k == Nciro)
{
break;
}
}
}
else if (Nciri == Nciro || CircuitInfo.UnequalCir[i] == Nciri + 1 + j)
{
if (index == 0)
{
if (Nciri == Nciro)
{
mr_ciro.CopyTo(mr_ciri, 0);
}
else
{
if (restartDP_index == 1 || !priconverge.flag)
{
mr_ciri[k] = mr_ciri_base[j][k] * mr_ciro[j] / (mr / Nciro);
}
else
{
mr_ciri[k] = mr_ciro[j] / Ngroupin[j];
}
}
}
//else mr_ciri_base.CopyTo(mr_ciri[k], 0);
//for (int i = 0; i < Ncir; i++)
r[i] = Circuit.CircuitCalc(i, cirArr, CircuitInfo, Nrow, Ntube, Nelement, fluid, composition, dh, l, geo, ta,
tri_cir[i], pri_cir[i], hri_cir[i], mr_ciri[k], ma, ha, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater);
r1[k] = r[i].ShallowCopy();
index_cir[k] = i;
k++;
if (k == (Nciri == Nciro ? Ncir : Ngroupin[j]))
{
Ncir_forDP = (Nciri == Nciro ? Ncir : Ngroupin[j]);
mr_forDP = (double[])mr_ciri.Clone();
break;
}
}
}
if (index_outbig && flag_ciro == 1)
{
break;
}
index++;
//dPconverge = CheckDPforCircuits.CheckDPConverge(mr, mr_ciri, r, Ncir);
dPconverge = CheckDPforCircuits.CheckDPConverge(flag_ciro, mr_forDP, r1, Ncir_forDP);
if (flag_ciro == 0)
{
restartDP_index = 0;
if (!dPconverge.flag)
{
dPconverge.mr.CopyTo(mr_ciri, 0); //mr_ciri = dPconverge.mr;
}
}
else //(flag_ciro == 1)
{
if (dPconverge.flag)
{
restartDP_index = 0;
}
else
{
restartDP_index = 1;
dPconverge.mr.CopyTo(mr_ciro, 0); //mr_ciro = dPconverge.mr;
break;
}
}
iterforDP++;
N_tube2 = 0;
#region //Result print out
if (dPconverge.flag)
{
if (Nciri == Nciro)
{
te_calc = Refrigerant.SATP(fluid, composition, r[j].Pro, 1).Temperature;
}
else
{
if (mr_ciri_base.Count < Nciro)
{
mr_ciri_base.Add(mr_forDP); //keep original mr ratio for fast iter
}
j = (flag_ciro == 1 ? j + Nciro : j);
res_cir2[j] = new CalcResult();
for (int i = 0; i < (flag_ciro == 1 ? Nciro : Ngroupin[j]); i++)
{
res_cir2[j].Q += r1[i].Q;
res_cir2[j].M += r1[i].M;
res_cir2[j].hro += (flag_ciro == 1 ? mr_ciro[i] : mr_ciri[i]) * r1[i].hro;
if (fluid[0] == "Water")
{
res_cir2[j].Tro += (flag_ciro == 1 ? mr_ciro[i] : mr_ciri[i]) * r1[i].Tro;
}
res_cir2[j].Vel_r = r1[i].Vel_r;
res_cir2[j].href += r1[i].href * CircuitInfo.TubeofCir[index_cir[i]];
res_cir2[j].R_1 += r1[i].R_1 * CircuitInfo.TubeofCir[index_cir[i]];
res_cir2[j].R_1a += r1[i].R_1a * CircuitInfo.TubeofCir[index_cir[i]];
res_cir2[j].R_1r += r1[i].R_1r * CircuitInfo.TubeofCir[index_cir[i]];
N_tube2 += CircuitInfo.TubeofCir[index_cir[i]];
}
res_cir2[j].DP = r1[(flag_ciro == 1 ? Nciro : Ngroupin[j]) - 1].DP;
res_cir2[j].Tao_Detail = ta;
res_cir2[j].Pro = r1[(flag_ciro == 1 ? Nciro : Ngroupin[j]) - 1].Pro;
res_cir2[j].hro = res_cir2[j].hro / (flag_ciro == 1 ? mr : mr_ciro[j]);
res_cir2[j].href = res_cir2[j].href / N_tube2;
res_cir2[j].R_1 = res_cir2[j].R_1 / N_tube2;
res_cir2[j].R_1a = res_cir2[j].R_1a / N_tube2;
res_cir2[j].R_1r = res_cir2[j].R_1r / N_tube2;
te_calc = Refrigerant.SATP(fluid, composition, res_cir2[j].Pro, 1).Temperature;
if (fluid[0] == "Water")
{
res_cir2[j].Tro = res_cir2[j].Tro / (flag_ciro == 1 ? mr : mr_ciro[j]);
}
else
{
res_cir2[j].Tro = Refrigerant.PHFLSH(fluid, composition, res_cir2[j].Pro, (res_cir2[j].hro + 140) * wm).t - 273.15;
}
}
}
#endregion
} while (!dPconverge.flag && iterforDP < 100);
if (Nciri == Nciro)
{
break;
}
if (index_outbig && j == Nciro - 1)
{
for (int i = 0; i < Nciro; i++)
{
r2[i].DP += res_cir2[i].DP;
}
flag_ciro = 1;
Ncir_forDP = Nciro;
mr_forDP = (double[])mr_ciro.Clone(); // mr_forDP = mr_ciro
dPconverge = CheckDPforCircuits.CheckDPConverge(flag_ciro, mr_forDP, r2, Ncir_forDP);
if (!dPconverge.flag)
{
restartDP_index = 1;
dPconverge.mr.CopyTo(mr_ciro, 0); //mr_ciro = dPconverge.mr;
}
break;
}
}
using (StreamWriter wr = File.AppendText(@"D:\Work\Simulation\Test\MinNout.txt"))
{
for (int i = 0; i < Ncir; i++)
{
wr.WriteLine("Q, {0}, DP, {1}, href, {2}, Ra_ratio, {3}, Tao, {4}, Tro, {5}, mr, {6}", r[i].Q, r[i].DP, r[i].href, r[i].Ra_ratio, r[i].Tao, r[i].Tro, r[i].mr);
}
}
if (restartDP_index == 1)
{
priconverge.flag = false;
}
else if (hexType == 0 && (fluid[0] != "Water"))
{
priconverge = CheckPin.CheckPriConverge(te, te_calc - 273.15, pri, pe, r[Ncir - 1].Pro); //res_slab.Pro
iterforPri++;
pri = priconverge.pri;
if (priconverge.flag && iterforPri == 1 && iterforDP == 1)
{
priconverge.flag = false; //to avoid not even iterate but converge by chance
}
}
else
{
priconverge.flag = true;
}
} while (!priconverge.flag && iterforPri < 20);
if (iterforDP >= 100)
{
throw new Exception("iter for DPConverge > 100.");
}
if (iterforPri >= 20)
{
throw new Exception("iter for DPPri > 20.");
}
#region //Result print out
for (int i = 0; i < Ncir; i++)
{
res_slab.Q += r[i].Q;
res_slab.M += r[i].M;
if (Nciri == Nciro)
{
res_slab.hro += mr_ciri[i] * r[i].hro;
}
res_slab.href += r[i].href * CircuitInfo.TubeofCir[i];
res_slab.R_1 += r[i].R_1 * CircuitInfo.TubeofCir[i];
res_slab.R_1a += r[i].R_1a * CircuitInfo.TubeofCir[i];
res_slab.R_1r += r[i].R_1r * CircuitInfo.TubeofCir[i];
}
if (Nciri == Nciro)
{
res_slab.hro = res_slab.hro / mr;
res_slab.Pro = r[Ncir - 1].Pro;
res_slab.Vel_r = r[Ncir - 1].Vel_r;
}
else if (!index_outbig)
{
res_slab.hro = res_cir2[Nciro].hro;
res_slab.Pro = res_cir2[Nciro].Pro;
res_slab.Vel_r = res_cir2[Nciro].Vel_r;
}
else
{
res_slab.hro = res_cir2[Nciro - 1].hro;
res_slab.Pro = res_cir2[Nciro - 1].Pro;
res_slab.Vel_r = res_cir2[Nciro - 1].Vel_r;
}
res_slab.Pri = pri;
res_slab.Tri = tri;
res_slab.hri = hri;
res_slab.mr = mr;
res_slab.DP = pri - res_slab.Pro;
res_slab.Tao_Detail = ta;
res_slab.href = res_slab.href / N_tube_total;
res_slab.ha = ha;
res_slab.R_1 = res_slab.R_1 / N_tube_total;
res_slab.R_1a = res_slab.R_1a / N_tube_total;
res_slab.R_1r = res_slab.R_1r / N_tube_total;
te_calc = Refrigerant.SATP(fluid, composition, res_slab.Pro, 1).Temperature;
double densityLo = Refrigerant.SATT(fluid, composition, te_calc, 1).DensityL; //mol/L
double densityVo = Refrigerant.SATT(fluid, composition, te_calc, 2).DensityV; //mol/L
//double wm = Refrigerant.WM(fluid, composition).Wm;
double hlo = Refrigerant.ENTHAL(fluid, composition, te_calc, densityLo).Enthalpy / wm - (fluid[0] == "Water" ? 0 : 140);
double hvo = Refrigerant.ENTHAL(fluid, composition, te_calc, densityVo).Enthalpy / wm - (fluid[0] == "Water" ? 0 : 140);
res_slab.x_o = (res_slab.hro - hlo) / (hvo - hlo);
double densityLi = Refrigerant.SATT(fluid, composition, tri + 273.15, 1).DensityL; //mol/L
double densityVi = Refrigerant.SATT(fluid, composition, tri + 273.15, 2).DensityV; //mol/L
//double wm = Refrigerant.WM(fluid, composition).Wm;
double hli = Refrigerant.ENTHAL(fluid, composition, tri + 273.15, densityLi).Enthalpy / wm - (fluid[0] == "Water" ? 0 : 140);
double hvi = Refrigerant.ENTHAL(fluid, composition, tri + 273.15, densityVi).Enthalpy / wm - (fluid[0] == "Water" ? 0 : 140);
res_slab.x_i = (res_slab.hri - hli) / (hvi - hli);
res_slab.Tro = Refrigerant.PHFLSH(fluid, composition, res_slab.Pro, (res_slab.hro + (fluid[0] == "Water" ? 0 : 140)) * wm).t - 273.15;
for (int j = 0; j < N_tube; j++)
{
for (int i = 0; i < Nelement; i++)
{
res_slab.Tao += res_slab.Tao_Detail[i, j, Nrow];
}
}
res_slab.Tao = res_slab.Tao / N_tube;
res_slab.Ra_ratio = res_slab.R_1a / res_slab.R_1;
res_slab.ma = ma;
res_slab.Va = ma / 1.2 * 3600;
return(res_slab);
#endregion
}
//Start
public static CalcResult SlabCalc(int[,] CirArrange, CircuitNumber CircuitInfo, int Nrow, int[] Ntube, int Nelement, string fluid, //double Npass, int[] N_tubes_pass,
double l, Geometry geo, double[, ,] ta, double[, ,] RH, double te, double pe, double hri, double mr, double[,] ma, double[,] ha, double[,] haw,
double eta_surface, double zh, double zdp, int hexType, double thickness, double conductivity, double Pwater, int AirFlowDirection, List <NodeInfo> Nodes, int N_Node, CapiliaryInput cap_inlet, CapiliaryInput cap_outlet, AbstractState coolprop, double[,] SourceTableData)
{
double tri = te; //Refrigerant.SATP(fluid, composition, pri, 1).Temperature - 273.15;
double pri = pe;
double te_calc_org = 0;
int N_tube = Ntube[0];
CirArr[] cirArr = new CirArr[Nrow * N_tube];
CirArrforAir cirArrforAir = new CirArrforAir();
cirArrforAir = CirArrangement.ReadCirArr(CirArrange, CircuitInfo, Nrow, Ntube, AirFlowDirection);
cirArr = cirArrforAir.CirArr;
CheckAir airConverge = new CheckAir();
CheckDP dPconverge = new CheckDP();
CheckPri priconverge = new CheckPri();
int iterforAir = 0;
int iterforDP = 0;
int iterforPri = 0;
CalcResult r = new CalcResult();
CalcResult res_slab = new CalcResult();
int N_cir = CircuitInfo.TubeofCir.Length;
int N_tube_total = CircuitInfo.TubeofCir.Sum();
CalcResult[] res_cir = new CalcResult[N_cir];
double[, ,] taout_calc = new double[Nelement, N_tube, Nrow];
double[, ,] RHout_calc = new double[Nelement, N_tube, Nrow];
double[,] Q_detail = new double[N_tube, Nrow];//detail output
double[,] DP_detail = new double[N_tube, Nrow];
double[,] Tri_detail = new double[N_tube, Nrow];
double[,] Pri_detail = new double[N_tube, Nrow];
double[,] hri_detail = new double[N_tube, Nrow];
double[,] Tro_detail = new double[N_tube, Nrow];
double[,] Pro_detail = new double[N_tube, Nrow];
double[,] hro_detail = new double[N_tube, Nrow];//
double[,] href_detail = new double[N_tube, Nrow];
double[,] mr_detail = new double[N_tube, Nrow];
double[,] charge_detail = new double[N_tube, Nrow];
int index_Node = 0;
int index_couple = 0;
int index_last_Node = 0;
bool index_end = false;
int[] index_status = new int[N_Node];
int[] index_FullStatus = new int[N_Node];
int[] index_DP = new int[N_Node];//
double mri_cal = 0;
double pri_cal = 0;
double hri_cal = 0;
double tri_cal = 0;
int index_cir = 0;
double te_calc;
do //Pri iteration
{
if (hexType == 0) // need to be under pri iteration
{
tri = CoolProp.PropsSI("T", "P", pri * 1000, "Q", 0, fluid) - 273.15;
}
else
{
te = CoolProp.PropsSI("T", "P", pri * 1000, "Q", 0, fluid) - 273.15;
}
int iiii = 0;
do//Air iteration
{
//Node[] Nodes=new Node[N_node];
//Nodes[0]=First Node;
for (int i = 0; i < N_Node; i++)
{
index_FullStatus[i] = Math.Max(Nodes[i].N_in, Nodes[i].N_out);
index_status[i] = 0;
}
index_Node = SearchNode.FindNextNode(-1, 0, Nodes, N_Node); //Find the first node
index_status[index_Node] = 0; //initial status
//index_FullStatus[index_Node]=Nodes[index_Node].N_out;
//index_DP[index_Node] = 0;
index_end = false;
for (int i = 0; i < 1000; i++) //for DP converge
{
if (Nodes[index_Node].inlet[0] == -1) //First Node(diverse)
{
Nodes[index_Node].mri[0] = mr;
Nodes[index_Node].pri[0] = pri;
Nodes[index_Node].hri[0] = hri;
Nodes[index_Node].tri[0] = tri;
//index_cir=Node[index].out[i];//?
}
if (Nodes[index_Node].type == 'D' || Nodes[index_Node].type == 'S') //Diverse Node Distribution
{
if (index_DP[index_Node] == 0) //first time calculation
{
Nodes[index_Node].mro[index_status[index_Node]] = Nodes[index_Node].mri[0] / Nodes[index_Node].N_out;
}
else
{
Nodes[index_Node].mro[index_status[index_Node]] = Nodes[index_Node].mri[0] * Nodes[index_Node].mr_ratio[index_status[index_Node]];
}
Nodes[index_Node].pro[index_status[index_Node]] = Nodes[index_Node].pri[0];
Nodes[index_Node].tro[index_status[index_Node]] = Nodes[index_Node].tri[0];
Nodes[index_Node].hro[index_status[index_Node]] = Nodes[index_Node].hri[0];
if (Nodes[index_Node].outType[index_status[index_Node]] == 0)//0:out is Node,1:out is Circuit
{
index_last_Node = index_Node;
index_Node = SearchNode.FindNextNode(index_Node, index_status[index_Node], Nodes, N_Node);//find node
Nodes[index_Node].pri[0] = Nodes[index_last_Node].pro[index_status[index_last_Node]];
Nodes[index_Node].tri[0] = Nodes[index_last_Node].tro[index_status[index_last_Node]];
Nodes[index_Node].hri[0] = Nodes[index_last_Node].hro[index_status[index_last_Node]];
Nodes[index_Node].mri[0] = Nodes[index_last_Node].mro[index_status[index_last_Node]];
index_status[index_Node] = 0;
//i_end=Nodes[index_Node].N_out;
continue;
}
else//out is Circuit
{
mri_cal = Nodes[index_Node].mro[index_status[index_Node]];
pri_cal = Nodes[index_Node].pro[index_status[index_Node]];
tri_cal = Nodes[index_Node].tro[index_status[index_Node]];
hri_cal = Nodes[index_Node].hro[index_status[index_Node]];
index_cir = Nodes[index_Node].outlet[index_status[index_Node]];
}
}
else if (Nodes[index_Node].type == 'C')//Converge Node
{
mri_cal = Nodes[index_Node].mro[0];
pri_cal = Nodes[index_Node].pro[0];
hri_cal = Nodes[index_Node].hro[0];
tri_cal = Nodes[index_Node].tro[0];
index_cir = Nodes[index_Node].outlet[0];
}
//index_status[index_Node]=i;//status
r = Circuit.CircuitCalc(index_cir, cirArr, CircuitInfo, Nrow, Ntube, Nelement, fluid, l, geo, ta, RH,
tri_cal, pri_cal, hri_cal, mri_cal, ma, ha, haw, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater, cap_inlet, cap_outlet, coolprop, SourceTableData);
res_cir[index_cir] = r;
if (r.Pro < 0)
{
res_slab.Pro = -10000000; return(res_slab);
}
for (int aa = 0; aa < Nrow; aa++)// detail result print
{
for (int bb = 0; bb < N_tube; bb++)
{
Q_detail[bb, aa] = r.Q_detail[bb, aa] == 0 ? Q_detail[bb, aa] : r.Q_detail[bb, aa];
DP_detail[bb, aa] = r.DP_detail[bb, aa] == 0 ? DP_detail[bb, aa] : r.DP_detail[bb, aa];
Pri_detail[bb, aa] = r.Pri_detail[bb, aa] == 0 ? Pri_detail[bb, aa] : r.Pri_detail[bb, aa];
Tri_detail[bb, aa] = r.Tri_detail[bb, aa] == 0 ? Tri_detail[bb, aa] : r.Tri_detail[bb, aa];
hri_detail[bb, aa] = r.hri_detail[bb, aa] == 0 ? hri_detail[bb, aa] : r.hri_detail[bb, aa];
Pro_detail[bb, aa] = r.Pro_detail[bb, aa] == 0 ? Pro_detail[bb, aa] : r.Pro_detail[bb, aa];
Tro_detail[bb, aa] = r.Tro_detail[bb, aa] == 0 ? Tro_detail[bb, aa] : r.Tro_detail[bb, aa];
hro_detail[bb, aa] = r.hro_detail[bb, aa] == 0 ? hro_detail[bb, aa] : r.hro_detail[bb, aa];
href_detail[bb, aa] = r.href_detail[bb, aa] == 0 ? href_detail[bb, aa] : r.href_detail[bb, aa];
mr_detail[bb, aa] = r.mr_detail[bb, aa] == 0 ? mr_detail[bb, aa] : r.mr_detail[bb, aa];
charge_detail[bb, aa] = r.charge_detail[bb, aa] == 0 ? charge_detail[bb, aa] : r.charge_detail[bb, aa];
for (int cc = 0; cc < Nelement; cc++)
{
taout_calc[cc, bb, aa] = r.Tao_Detail[cc, bb, aa] == 0 ? taout_calc[cc, bb, aa] : r.Tao_Detail[cc, bb, aa];
RHout_calc[cc, bb, aa] = r.RHo_Detail[cc, bb, aa] == 0 ? RHout_calc[cc, bb, aa] : r.RHo_Detail[cc, bb, aa];
}
}
}
index_Node = SearchNode.FindNextNode(index_Node, index_status[index_Node], Nodes, N_Node); //Find Next Node
if (Nodes[index_Node].type == 'D') //Diverse Node cal
{
index_status[index_Node] = 0;
//i_end=Nodes[index_Node].N_out;
Nodes[index_Node].pri[0] = r.Pro;
Nodes[index_Node].tri[0] = r.Tro;
Nodes[index_Node].hri[0] = r.hro;
Nodes[index_Node].mri[0] = r.mr;
continue;
}
else if (Nodes[index_Node].type == 'C')//Converge Node cal
{
for (int ii = 0; ii < N_Node; ii++)
{
if (Nodes[ii].couple == Nodes[index_Node].couple && ii != index_Node)
{
index_couple = ii;
break;
}
}//Find the Couple Node
Nodes[index_Node].pri[index_status[index_couple]] = r.Pro;
Nodes[index_Node].tri[index_status[index_couple]] = r.Tro;
Nodes[index_Node].hri[index_status[index_couple]] = r.hro;
Nodes[index_Node].mri[index_status[index_couple]] = r.mr;
for (int k = 0; k < 100; k++)
{
if (index_status[index_couple] < index_FullStatus[index_couple] - 1)//continue
{
index_Node = index_couple;
index_status[index_Node]++;
//i_end=Nodes[index_couple].N_out;
break;
}
else if (index_status[index_couple] == index_FullStatus[index_couple] - 1)//dp converge calculation
{
//DPconverge=DPConverge(Nodes[index_Node].mri,Nodes[index_Node].pri);
dPconverge = CheckDPforCircuits.CheckDPConverge2(hexType, iterforPri, Nodes[index_Node].mri, Nodes[index_Node].pri, Nodes[index_couple].pri[0], index_FullStatus[index_couple]);
iterforDP++;
if (dPconverge.flag == false)//need to be modified
{
Nodes[index_couple].mro = dPconverge.mr;
Nodes[index_couple].mr_ratio = dPconverge.mr_ratio;
index_DP[index_couple]++;
index_Node = index_couple;
index_status[index_Node] = 0;
//i_end=Nodes[index_Node].N_out;
break;
}
else if (dPconverge.flag == true)
{
index_DP[index_couple] = 1;
Nodes[index_couple].mr_ratio = dPconverge.mr_ratio;
double mr_sum = 0;
double tro_ave = 0;
double hro_ave = 0;
for (int ii = 0; ii < index_FullStatus[index_couple]; ii++)
{
mr_sum += Nodes[index_Node].mri[ii];
}
for (int ii = 0; ii < index_FullStatus[index_couple]; ii++)
{
tro_ave += Nodes[index_Node].mri[ii] * Nodes[index_Node].tri[ii];
hro_ave += Nodes[index_Node].mri[ii] * Nodes[index_Node].hri[ii];
}
Nodes[index_Node].mro[0] = mr_sum;
Nodes[index_Node].tro[0] = tro_ave / mr_sum;
Nodes[index_Node].hro[0] = hro_ave / mr_sum;
Nodes[index_Node].pro[0] = Nodes[index_Node].pri[0];
if (Nodes[index_Node].outType[0] == 0)//0:out is Node
{
index_last_Node = index_Node;
index_Node = SearchNode.FindNextNode(index_Node, 0, Nodes, N_Node);
if (Nodes[index_Node].type == 'C')
{
for (int ii = 0; ii < N_Node; ii++)
{
if (Nodes[ii].couple == Nodes[index_Node].couple && ii != index_Node)
{
index_couple = ii;
break;
}
}//Find the Couple Node
Nodes[index_Node].pri[index_status[index_couple]] = Nodes[index_last_Node].pro[0];
Nodes[index_Node].tri[index_status[index_couple]] = Nodes[index_last_Node].tro[0];
Nodes[index_Node].hri[index_status[index_couple]] = Nodes[index_last_Node].hro[0];
Nodes[index_Node].mri[index_status[index_couple]] = Nodes[index_last_Node].mro[0];
continue;
}
else if (Nodes[index_Node].type == 'D')
{
index_status[index_Node] = 0;
Nodes[index_Node].mri[0] = Nodes[index_last_Node].mro[0];
Nodes[index_Node].tri[0] = Nodes[index_last_Node].tro[0];
Nodes[index_Node].hri[0] = Nodes[index_last_Node].hro[0];
Nodes[index_Node].pri[0] = Nodes[index_last_Node].pro[0];
break;
}
}
else if (Nodes[index_Node].outType[0] == -1)
{
index_end = true;
break;
}
else//out is Circuit
{
break;
}
continue;
}
} //end if
} //end for
if (index_end == true)
{
break;
}
}
else if (Nodes[index_Node].type == 'E')//for 1 out case
{
Nodes[index_Node].pri[0] = r.Pro;
Nodes[index_Node].tri[0] = r.Tro;
Nodes[index_Node].hri[0] = r.hro;
Nodes[index_Node].mro[0] = r.mr;
Nodes[index_Node].pro[0] = r.Pro;
Nodes[index_Node].tro[0] = r.Tro;
Nodes[index_Node].hro[0] = r.hro;
Nodes[index_Node].mro[0] = r.mr;
break;
}
}//end out for
airConverge = CheckAirConvergeforCircuits.CheckAirConverge2(cirArrforAir.TotalDirection, Nrow, N_tube, Nelement, ta, RH, taout_calc, RHout_calc); //taout_calc, RHout_calc
ta = airConverge.ta;
RH = airConverge.RH;
iterforAir++;
} while (airConverge.flag == false && iterforAir < 50);
if (hexType == 0)
{
te_calc = CoolProp.PropsSI("T", "P", Nodes[index_Node].pro[0] * 1000, "Q", 0, fluid);
priconverge = CheckPin.CheckPriConverge(te, te_calc - 273.15, te_calc_org - 273.15, pri, pe, Nodes[index_Node].pro[0]); //res_slab.Pro
iterforPri++;
pri = priconverge.pri;
te_calc_org = te_calc;
if (priconverge.flag && iterforPri == 1)
{
priconverge.flag = false; //to avoid not even iterate but converge by chance
}
}
} while (priconverge.flag == false && iterforPri < 20);
//result print
if (iterforDP >= 200)
{
return(res_slab);
throw new Exception("iter for DPConverge > 200.");
}
if (iterforPri >= 50)
{
return(res_slab);
throw new Exception("iter for Pri > 50.");
}
#region result print
for (int i = 0; i < N_cir; i++)
{
res_slab.R_1 += res_cir[i].R_1 * CircuitInfo.TubeofCir[i];
res_slab.R_1a += res_cir[i].R_1a * CircuitInfo.TubeofCir[i];
res_slab.R_1r += res_cir[i].R_1r * CircuitInfo.TubeofCir[i];
}
for (int j = 0; j < N_tube; j++)//detail output
{
for (int k = 0; k < Nrow; k++)
{
res_slab.Q += Q_detail[j, k];
res_slab.M += charge_detail[j, k];
res_slab.href += href_detail[j, k];
}
}
res_slab.hro = Nodes[index_Node].hro[0];
res_slab.Pro = Nodes[index_Node].pro[0];
res_slab.Pri = pri;
res_slab.Tri = tri;
res_slab.hri = hri;
res_slab.mr = mr;
res_slab.DP = pri - res_slab.Pro;
res_slab.Tao_Detail = ta;
res_slab.RHo_Detail = RH;
res_slab.href = res_slab.href / N_tube_total;
for (int i = 0; i < ha.GetLength(0); i++)
{
for (int j = 0; j < ha.GetLength(1); j++)
{
res_slab.ha += ha[i, j];
}
}
res_slab.ha = res_slab.ha / ha.Length;
res_slab.R_1 = res_slab.R_1 / N_tube_total;
res_slab.R_1a = res_slab.R_1a / N_tube_total;
res_slab.R_1r = res_slab.R_1r / N_tube_total;
te_calc = CoolProp.PropsSI("T", "P", res_slab.Pro * 1000, "Q", 0, fluid) - 273.15;
double densityLo = CoolProp.PropsSI("D", "T", te_calc + 273.15, "Q", 0, fluid);
double densityVo = CoolProp.PropsSI("D", "T", te_calc + 273.15, "Q", 1, fluid);
double hlo = CoolProp.PropsSI("H", "T", te_calc + 273.15, "D", densityLo, fluid) / 1000;
double hvo = CoolProp.PropsSI("H", "T", te_calc + 273.15, "D", densityVo, fluid) / 1000;
res_slab.x_o = (res_slab.hro - hlo) / (hvo - hlo);
double hli = CoolProp.PropsSI("H", "P", pri * 1000, "Q", 0, fluid) / 1000;
double hvi = CoolProp.PropsSI("H", "P", pri * 1000, "Q", 1, fluid) / 1000;
res_slab.x_i = (res_slab.hri - hli) / (hvi - hli);
res_slab.Tro = CoolProp.PropsSI("T", "P", res_slab.Pro * 1000, "H", res_slab.hro * 1000, fluid) - 273.15;
for (int j = 0; j < N_tube; j++)
{
for (int i = 0; i < Nelement; i++)
{
res_slab.Tao += res_slab.Tao_Detail[i, j, Nrow];
res_slab.RHout += res_slab.RHo_Detail[i, j, Nrow];
}
}
res_slab.Tao = res_slab.Tao / (N_tube * Nelement);
res_slab.RHout = res_slab.RHout / (N_tube * Nelement);
res_slab.Ra_ratio = res_slab.R_1a / res_slab.R_1;
for (int i = 0; i < ma.GetLength(0); i++)
{
for (int j = 0; j < ma.GetLength(1); j++)
{
res_slab.ma += ma[i, j];
}
}
res_slab.Va = res_slab.ma / 1.2 * 3600;
res_slab.Q_detail = Q_detail;//detail output
res_slab.DP_detail = DP_detail;
res_slab.Tri_detail = Tri_detail;
res_slab.Pri_detail = Pri_detail;
res_slab.hri_detail = hri_detail;
res_slab.Tro_detail = Tro_detail;
res_slab.Pro_detail = Pro_detail;
res_slab.hro_detail = hro_detail;
res_slab.href_detail = href_detail;
res_slab.mr_detail = mr_detail;
res_slab.Aa = geo.TotalArea.A_a;
res_slab.Ar = geo.TotalArea.A_r;
res_slab.AHx = geo.TotalArea.A_hx;
res_slab.N_row = Nrow;
res_slab.tube_row = N_tube;
res_slab.charge_detail = charge_detail;
#endregion
return(res_slab);
} //end function