public static CalcResult ElementCal(string[] fluid, double[] composition, double dh, double l,
double Aa_fin, double Aa_tube, double A_r_cs, double Ar, double tai,
double tri, double pri, double hri, double mr, double g, double ma, double ha,
double eta_surface, double zh, double zdp, int hexType, double thickness, double conductivity, double Pwater)
{
var r = new Refrigerant.SATTTotalResult();
double href = 0;
double gg = 9.8;
double tsat;
int phase1 = 1;
int phase2 = 2;
tsat = Refrigerant.SATP(fluid, composition, pri, phase1).Temperature;
r = Refrigerant.SATTTotal(fluid, composition, tsat).SATTTotalResult;
double Vol_tubes = A_r_cs * l; //Tube volume, for charge calculation
double h_l = r.EnthalpyL;
double h_v = r.EnthalpyV;
double Tri_mod;
double alpha;
double M;
var res_element = new CalcResult();
// **********Superheated state**********
if (hri > h_v && fluid[0] != "Water")
{
if (hri < 1.02 * h_v)
{
Tri_mod = tri + 0.5; //"for Tri modification in the transition region"
}
else
{
Tri_mod = tri;
}
res_element = SPElement.ElementCalc(fluid, composition, dh, l, Aa_fin, Aa_tube, A_r_cs, Ar, tai, Tri_mod, pri, hri, mr, g, ma, ha, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater);
res_element.x_i = (hri - h_l) / (h_v - h_l);
res_element.x_o = (res_element.hro - h_l) / (h_v - h_l);
alpha = 1; //set void fraction to 1 to identify a superheated state
//{Equations are for charge calculation}
double T_avg = (tri + res_element.Tro) / 2; //Average temperature of the element
double P_avg = (pri + res_element.Pro) / 2; //Average pressure of the element
double rho = Refrigerant.TPFLSH(fluid, composition, T_avg + 273.15, P_avg).D *r.Wm; //density(ref$, T=T_avg, P=P_avg)
M = Vol_tubes * rho; //"Mass calculated"
}
// **********Twophase state**********"
if (hri <= h_v && hri >= h_l && fluid[0] != "Water")
{
res_element = TPElement.ElementCalc(fluid, composition, dh, l, Aa_fin, Aa_tube, A_r_cs, Ar, tai, tri, pri, hri, mr, g, ma, ha, eta_surface, zh, zdp, hexType, thickness, conductivity);
//x=x_o "outlet quality of the element"
double x_avg = (res_element.x_i + res_element.x_o) / 2; //Average quality of the element
if (x_avg < 0) //"If negative, set quality to inlet value"
{
x_avg = res_element.x_i;
}
double T_avg = (tri + res_element.Tro) / 2; //Average temperature of the element
//Call VOIDFRACTION_charge(ref$,x_avg, T_avg, G_r, Dh: alpha) "Average void fraction of the element"
alpha = 1;
//{Equations are for charge calculation}
double P_avg = (pri + res_element.Pro) / 2; //Average pressure of the element
double rho_l = Refrigerant.SATP(fluid, composition, P_avg, phase1).DensityL *r.Wm;
double rho_v = Refrigerant.SATP(fluid, composition, P_avg, phase1).DensityV *r.Wm;
//{Call VOIDFRACTION_pressure(ref$, x_avg, P_avg : alpha_p) "Baroczy void fraction model" }
M = Vol_tubes * (alpha * rho_v + (1 - alpha) * rho_l); //Mass calculated
}
//**********Subcooled state**********
if (hri < h_l || fluid[0] == "Water")
{
if (hri > 0.98 * h_l)
{
Tri_mod = tri - 0.5; //"for Tri modification in the transition region"
}
else
{
Tri_mod = tri;
}
if (fluid[0] == "Water")
{
Tri_mod = tri - 0.0001;
}
res_element = SPElement.ElementCalc(fluid, composition, dh, l, Aa_fin, Aa_tube, A_r_cs, Ar, tai, Tri_mod, pri, hri, mr, g, ma, ha, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater);
//Call SUBCOOLED(ref$, Dh, L, A_a, A_r, Tai, Tri_mod, Pri, hri, m_r, G_r, m_a, h_air, eta_surface: Tro, Pro, hro, Tao, Q, h_ref, R_1, R_1a, R_1r, DELTAP, Vel_r )
res_element.x_i = (hri - h_l) / (h_v - h_l);
res_element.x_o = (res_element.hro - h_l) / (h_v - h_l);
//x=x_o "outlet quality of the element"
//{x=-1 "set quality to -100 to identify a subcooled state"}
alpha = -1; //set void fraction to -100 to identify a subcooled state
//{Equations are for charge calculation}
double T_avg = (tri + res_element.Tro) / 2; //Average temperature of the element
double P_avg = (pri + res_element.Pro) / 2; //Average pressure of the element
double rho = Refrigerant.TPFLSH(fluid, composition, T_avg, P_avg).D *r.Wm; //density(ref$, T=T_avg, P=P_avg)
M = Vol_tubes * rho; //Mass calculated
}
return(res_element);
}
public static CalcResult ElementCal(string fluid, double l,
double Aa_fin, double Aa_tube, double A_r_cs, double Ar, Geometry geo, double tai, double RHi,
double tri, double pri, double hri, double mr, double g, double ma, double ha, double haw,
double eta_surface, double zh, double zdp, int hexType, double thickness, double conductivity, double Pwater, AbstractState coolprop, double[,] SourceTableData)
{
//AbstractState coolprop = AbstractState.factory("HEOS", fluid);
double tsat;
coolprop.update(input_pairs.PQ_INPUTS, pri * 1000, 0);
tsat = coolprop.T();
//tsat = CoolProp.PropsSI("T", "P", pri * 1000, "Q", 0, fluid);
double Vol_tubes = A_r_cs * l; //Tube volume, for charge calculation
coolprop.update(input_pairs.QT_INPUTS, 0, tsat);
double h_l = coolprop.hmass() / 1000;
//double h_l = CoolProp.PropsSI("H", "T", tsat, "Q", 0, fluid) / 1000 ;
coolprop.update(input_pairs.QT_INPUTS, 1, tsat);
double h_v = coolprop.hmass() / 1000;
//double h_v1 = CoolProp.PropsSI("H", "T", tsat, "Q", 1, fluid) / 1000 ;
double Tri_mod;
double alpha;
var res_element = new CalcResult();
// **********Superheated state**********
if (hri > h_v && fluid != "Water")
{
if (hri < 1.02 * h_v) //ignore modification for now, ruhao,20180226
{
Tri_mod = tri + 0.5 * 0; //"for Tri modification in the transition region"
}
else
{
Tri_mod = tri;
}
res_element = SPElement.ElementCalc3(fluid, l, Aa_fin, Aa_tube, A_r_cs, Ar, geo, tai, RHi, Tri_mod, pri, hri, mr, g, ma, ha, haw, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater, coolprop, SourceTableData);
if (res_element.Pro < 0)
{
return(res_element);
}
res_element.x_i = (hri - h_l) / (h_v - h_l);
res_element.x_o = (res_element.hro - h_l) / (h_v - h_l);
alpha = 1; //set void fraction to 1 to identify a superheated state
//{Equations are for charge calculation}
coolprop.update(input_pairs.HmassP_INPUTS, hri * 1000, pri * 1000);
double rho = coolprop.rhomass();
res_element.M = Vol_tubes * rho; //"Mass calculated"
}
// **********Twophase state**********"
if (hri <= h_v && hri >= h_l && fluid != "Water")
{
res_element = TPElement.ElementCalc1(fluid, l, Aa_fin, Aa_tube, A_r_cs, Ar, geo, tai, RHi, tri, pri, hri, mr, g, ma, ha, haw, eta_surface, zh, zdp, hexType, thickness, conductivity, coolprop, SourceTableData);
if (res_element.Pro < 0)
{
return(res_element);
}
//x=x_o "outlet quality of the element"
double x_avg = (res_element.x_i + res_element.x_o) / 2; //Average quality of the element
if (x_avg < 0) //"If negative, set quality to inlet value"
{
x_avg = res_element.x_i;
}
double T_avg = (tri + res_element.Tro) / 2; //Average temperature of the element
//Call VOIDFRACTION_charge(ref$,x_avg, T_avg, G_r, Dh: alpha) "Average void fraction of the element"
alpha = 1;
//{Equations are for charge calculation}
double P_avg = (pri + res_element.Pro) / 2; //Average pressure of the element
coolprop.update(input_pairs.PQ_INPUTS, P_avg * 1000, 0);
double rho_l = coolprop.rhomass();
//double rho_l = CoolProp.PropsSI("D", "P", P_avg * 1000, "Q", 0, fluid);
coolprop.update(input_pairs.PQ_INPUTS, P_avg * 1000, 1);
double rho_v = coolprop.rhomass();
//double rho_v = CoolProp.PropsSI("D", "P", P_avg * 1000, "Q", 1, fluid);
//{Call VOIDFRACTION_pressure(ref$, x_avg, P_avg : alpha_p) "Baroczy void fraction model" }
double alpha_homog = 1 / (1 + (1 - x_avg) / x_avg * (rho_v / rho_l)); // Homogeneous model, Intermittent flow void fraction
res_element.M = Vol_tubes * (alpha_homog * rho_v + (1 - alpha_homog) * rho_l); //Mass calculated
}
//**********Subcooled state**********
if (hri < h_l || fluid == "Water")
{
if (hri > 0.98 * h_l) //ignore modification for now, ruhao,20180226
{
Tri_mod = tri - 0.5 * 0; //"for Tri modification in the transition region"
}
else
{
Tri_mod = tri;
}
if (fluid == "Water")
{
Tri_mod = tri - 0.0001;
}
res_element = SPElement.ElementCalc3(fluid, l, Aa_fin, Aa_tube, A_r_cs, Ar, geo, tai, RHi, Tri_mod, pri, hri, mr, g, ma, ha, haw, eta_surface, zh, zdp, hexType, thickness, conductivity, Pwater, coolprop, SourceTableData);
if (res_element.Pro < 0)
{
return(res_element);
}
//Call SUBCOOLED(ref$, Dh, L, A_a, A_r, Tai, Tri_mod, Pri, hri, m_r, G_r, m_a, h_air, eta_surface: Tro, Pro, hro, Tao, Q, h_ref, R_1, R_1a, R_1r, DELTAP, Vel_r )
res_element.x_i = (hri - h_l) / (h_v - h_l);
res_element.x_o = (res_element.hro - h_l) / (h_v - h_l);
//x=x_o "outlet quality of the element"
//{x=-1 "set quality to -100 to identify a subcooled state"}
alpha = -1; //set void fraction to -100 to identify a subcooled state
//{Equations are for charge calculation}
coolprop.update(input_pairs.HmassP_INPUTS, hri * 1000, pri * 1000);
double rho = coolprop.rhomass();
res_element.M = Vol_tubes * rho; //Mass calculated
}
return(res_element);
}
