//EQUATIONS
//Units: SI, Temperatures in Celsius, pressures in bar
//Project file: D:\_classement\_Thopt\THERMOPTIM_Pro_282\proj\refrig_light_SC10.prj
//Date and Time: 2024-08-24 17:31:39

//Flow rate unit: -

//GAS COMPOSITIONS

//PROCESSES

//Process: refrigeration effect
//Equation: 1
m_dot_refrigerationeffect = m_dot_throttling // Upstream process - throttling
// Comment = mDeltaH not set
//Equation: 2
Q_dot_refrigerationeffect = m_dot_refrigerationeffect*(h_1 - h_4) // Upstream point - 4 - Downstream point - 1 - DeltaH/flow

//Process: condenser
//Equation: 3
m_dot_condenser = m_dot_desuperheating // Upstream process - desuperheating
// Comment = mDeltaH not set
//Equation: 4
Q_dot_condenser = m_dot_condenser*(h_3 - h_3a) // Upstream point - 3a - Downstream point - 3 - DeltaH/flow

//Process: desuperheating
//Equation: 5
m_dot_desuperheating = m_dot_compressor // Upstream process - compressor
// Comment = mDeltaH not set
//Equation: 6
Q_dot_desuperheating = m_dot_desuperheating*(h_3a - h_2) // Upstream point - 2 - Downstream point - 3a - DeltaH/flow

//Process: compressor
//Equation: 7
m_dot_compressor = m_dot_refrigerationeffect // Upstream process - refrigeration effect
//Equation: 8
s_1 = s_Ph("R134a",p_1,h_1) // Upstream point - 1 - Downstream point - 2
// Comment = Isentropic reference
//Equation: 9
hs_2 = h_Ps("R134a",p_2,s_1) // Downstream point - 2
//Equation: 10
etaT_compressor = 0.75// Isentropic efficiency
//Equation: 11
h_2 = h_1 + (hs_2 - h_1)/etaT_compressor // Upstream point - 1 - Downstream point - 2
//Equation: 12
T_2 = T_Ph("R134a",p_2,h_2) // Downstream point - 2
// Comment = Given outlet pressure
//Equation: 13
p_2 = 12.0// Outlet pressure
//Equation: 14
W_dot_compressor = m_dot_compressor*(h_2 - h_1) // DeltaH 

//Process: throttling
// Comment = Isenthalpic throttling
//Equation: 15
//p_4 = 1.78// Given outlet pressure
//Equation: 16
xl_4 = 0.// Saturated liquid quality
//Equation: 17
Tl_4 = T_4- 0.01// Saturated liquid temperature
//Equation: 18
xv_4 = 1.// Saturated vapor quality
//Equation: 19
Tv_4 = T_4+ 0.01// Saturated vapor temperature
//Equation: 20
hl_4 = hsat_Px("R134a",p_4,xl_4)// Saturated liquid enthalpy
//Equation: 21
hv_4 = hsat_Px("R134a",p_4,xv_4)// Saturated vapor enthalpy
//Equation: 22
x_4 = (h_3 - hl_4)/(hv_4 - hl_4)// Quality
//Equation: 23
T_4 = Tsat_P("R134a",p_4) // Downstream point - 4
//Equation: 24
h_4 = hsat_Px("R134a",p_4,x_4) // Enthalpy

//NODES

//HEAT EXCHANGERS

//Number of equations: 24

//POINTS WITH SATURATION TEMPERATURE SET


//Point 1

//Outlet point of process refrigeration effect
//Equation: 25
p_1 = 1.78// P (bar)
//Equation: 26
x_1 = 1.0// Quality
//Equation: 27
dTsat_1 = 1.0// Deviation from Tsat
//Equation: 28
T_1 = Tsat_P("R134a",p_1)+dTsat_1// set Tsat (Celsius)
//Equation: 29
h_1 = hsat_Px("R134a",p_1,x_1)// Enthalpy


//Point 3a

//Outlet point of process desuperheating
//Equation: 30
p_3a = 12.0// P (bar)
//Equation: 31
x_3a = 1.0// Quality
//Equation: 32
dTsat_3a = 0.0// Deviation from Tsat
//Equation: 33
T_3a = Tsat_P("R134a",p_3a)+dTsat_3a// set Tsat (Celsius)
//Equation: 34
h_3a = hsat_Px("R134a",p_3a,x_3a)// Enthalpy


//Point 3

//Outlet point of process condenser
//Equation: 35
p_3 = 12.0// P (bar)
//Equation: 36
x_3 = 0.0// Quality
//Equation: 37
dTsat_3 = -2.0// Deviation from Tsat
//Equation: 38
T_3 = Tsat_P("R134a",p_3)+dTsat_3// set Tsat (Celsius)
//Equation: 39
h_3 = hsat_Px("R134a",p_3,x_3)// Enthalpy


//Point 4
//Equation: 40
p_4 = 1.78// P (bar)
//Equation: 41
//x_4 = 0.328174791// Quality
//Equation: 42
dTsat_4 = 0.0// Deviation from Tsat
//Equation: 43
T_4 = Tsat_P("R134a",p_4)+dTsat_4// set Tsat (Celsius)
//Equation: 44
//h_4 = hsat_Px("R134a",p_4,x_4)// Enthalpy


//OTHER POINTS WITH PRESSURE SET


//SET FLOW RATES


//OVERALL BALANCE

//Equation: 45
useful_Energy = Q_dot_refrigerationeffect
//Equation: 46
purchased_Energy = W_dot_compressor
//Equation: 47
eta_global = abs(useful_Energy/purchased_Energy)
//Equation: 48
m_dot_throttling=1