//EQUATIONS
//Units: SI, Temperatures in Celsius, pressures in bar
//Project file: D:\_classement\_Thopt\THERMOPTIM_Pro_282\proj\PBMR_IHX_N2_En.prj
//Date and Time: 2024-08-11 19:34:56

//Flow rate unit: -

//PROCESSES

//Process: IHX
//Equation: 1
m_dot_IHX = m_dot_HPregenerator // Upstream process - HP regenerator
// Comment = //Exchange process connected to a heat exchanger

//Process: IHX core
//Equation: 2
m_dot_IHXcore = m_dot_core // Upstream process - core
// Comment = //Exchange process connected to a heat exchanger

//Process: intercooler
//Equation: 3
m_dot_intercooler = m_dot_LPcompressor // Upstream process - LP compressor
// Comment = mDeltaH not set
//Equation: 4
h_10 =  calcH_TP("N2";T = T_10;P = p_10)// Enthalpy
//Equation: 5
Q_dot_intercooler/m_dot_intercooler = h_10 - h_9 // Upstream point - 9 - Downstream point - 10 - DeltaH/flow

//Process: HP regenerator
//Equation: 6
m_dot_HPregenerator = m_dot_HPcompressor // Upstream process - HP compressor
// Comment = //Exchange process connected to a heat exchanger

//Process: precooler
//Equation: 7
m_dot_precooler = m_dot_LPregenerator // Upstream process - LP regenerator
// Comment = mDeltaH not set
//Equation: 8
h_8 =  calcH_TP("N2";T = T_8;P = p_8)// Enthalpy
//Equation: 9
Q_dot_precooler/m_dot_precooler = h_8 - h_7 // Upstream point - 7 - Downstream point - 8 - DeltaH/flow

//Process: LP regenerator
//Equation: 10
m_dot_LPregenerator = m_dot_powerturbine // Upstream process - power turbine
// Comment = //Exchange process connected to a heat exchanger

//Process: core
// Comment = mDeltaH not set
//Equation: 11
h_coreoutlet =  calcH_TP("He";T = T_coreoutlet;P = p_coreoutlet)// Enthalpy
//Equation: 12
Q_dot_core/m_dot_core = h_coreoutlet - h_corecompressoroutlet // Upstream point - core compressor outlet - Downstream point - core outlet - DeltaH/flow

//Process: HP turbine
//Equation: 13
m_dot_HPturbine = m_dot_IHX // Upstream process - IHX
//Equation: 14
s_3 = calcS_PH("N2";P = p_3;H = h_3) // Upstream point - 3 - Downstream point - 4
// Comment = Polytropic reference
//Equation: 15
ds_4 = -(1 - etaT_HPturbine)*8.314/M_3*ln(p_4/p_3) // Upstream point - 3 - Downstream point - 4
//Equation: 16
s_4 = s_3 + ds_4 // Entropy - 4
//Equation: 17
M_3 = 28.0134 // Molar mass - 4
//Equation: 18
etaT_HPturbine = 0.9// Polytropic efficiency
//Equation: 19
h_4 = calcH_PS("N2";P = p_4;S = s_4) // Enthalpy
// Comment = Polytropic coefficient: k = -Math.log(aval.p/amont.p)/Math.log(aval.V/amont.V)
//Equation: 20
//h_4 = calcH_PS("N2";P = p_4;S = s_4) // Enthalpy
//Equation: 21
T_4 = calcT_PH("N2";P = p_4 ;H = h_4) // Downstream point - 4
// Comment = Given outlet pressure
//Equation: 22
p_4 = 55.908073701989174// Outlet pressure
//Equation: 23
W_dot_HPturbine = m_dot_HPturbine*(h_4 - h_3) // DeltaH 

//Process: LP compressor
//Equation: 24
m_dot_LPcompressor = m_dot_precooler // Upstream process - precooler
//Equation: 25
s_8 = calcS_PH("N2";P = p_8;H = h_8) // Upstream point - 8 - Downstream point - 9
// Comment = Polytropic reference
//Equation: 26
ds_9 = (1 - etaT_LPcompressor)/etaT_LPcompressor*8.314/M_8*ln(p_9/p_8) // Upstream point - 8 - Downstream point - 9
//Equation: 27
s_9 = s_8 + ds_9 // Entropy - 9
//Equation: 28
M_8 = 28.0134 // Molar mass - 9
//Equation: 29
etaT_LPcompressor = 0.9// Polytropic efficiency
// Comment = Polytropic coefficient: k = -Math.log(aval.p/amont.p)/Math.log(aval.V/amont.V)
//Equation: 30
h_9 = calcH_PS("N2";P = p_9;S = s_9) // Enthalpy
//Equation: 31
T_9 = calcT_PH("N2";P = p_9 ;H = h_9) // Downstream point - 9
// Comment = Given outlet pressure
//Equation: 32
p_9 = 42.9// Outlet pressure
//Equation: 33
W_dot_LPcompressor = m_dot_LPcompressor*(h_9 - h_8) // DeltaH 

//Process: HP compressor
//Equation: 34
m_dot_HPcompressor = m_dot_intercooler // Upstream process - intercooler
//Equation: 35
s_10 = calcS_PH("N2";P = p_10;H = h_10) // Upstream point - 10 - Downstream point - 1
// Comment = Polytropic reference
//Equation: 36
ds_1 = (1 - etaT_HPcompressor)/etaT_HPcompressor*8.314/M_10*ln(p_1/p_10) // Upstream point - 10 - Downstream point - 1
//Equation: 37
s_1 = s_10 + ds_1 // Entropy - 1
//Equation: 38
M_10 = 28.0134 // Molar mass - 1
//Equation: 39
etaT_HPcompressor = 0.9// Polytropic efficiency
// Comment = Polytropic coefficient: k = -Math.log(aval.p/amont.p)/Math.log(aval.V/amont.V)
//Equation: 40
h_1 = calcH_PS("N2";P = p_1;S = s_1) // Enthalpy
//Equation: 41
T_1 = calcT_PH("N2";P = p_1 ;H = h_1) // Downstream point - 1
// Comment = Given outlet pressure
//Equation: 42
p_1 = 70.7// Outlet pressure
//Equation: 43
W_dot_HPcompressor = m_dot_HPcompressor*(h_1 - h_10) // DeltaH 

//Process: LP turbine
//Equation: 44
m_dot_LPturbine = m_dot_HPturbine // Upstream process - HP turbine
//Equation: 45
//s_4 = calcS_PH("N2";P = p_4;H = h_4) // Upstream point - 4 - Downstream point - 5
// Comment = Polytropic reference
//Equation: 46
ds_5 = -(1 - etaT_LPturbine)*8.314/M_4*ln(p_5/p_4) // Upstream point - 4 - Downstream point - 5
//Equation: 47
s_5 = s_4 + ds_5 // Entropy - 5
//Equation: 48
M_4 = 28.0134 // Molar mass - 5
//Equation: 49
etaT_LPturbine = 0.9// Polytropic efficiency
//Equation: 50
h_5 = calcH_PS("N2";P = p_5;S = s_5) // Enthalpy
// Comment = Polytropic coefficient: k = -Math.log(aval.p/amont.p)/Math.log(aval.V/amont.V)
//Equation: 51
//h_5 = calcH_PS("N2";P = p_5;S = s_5) // Enthalpy
//Equation: 52
T_5 = calcT_PH("N2";P = p_5 ;H = h_5) // Downstream point - 5
// Comment = Given outlet pressure
//Equation: 53
p_5 = 45.773492624997495// Outlet pressure
//Equation: 54
W_dot_LPturbine = m_dot_LPturbine*(h_5 - h_4) // DeltaH 

//Process: power turbine
//Equation: 55
//s_5 = calcS_PH("N2";P = p_5;H = h_5) // Upstream point - 5 - Downstream point - 6
// Comment = Polytropic reference
//Equation: 56
ds_6 = -(1 - etaT_powerturbine)*8.314/M_5*ln(p_6/p_5) // Upstream point - 5 - Downstream point - 6
//Equation: 57
s_6 = s_5 + ds_6 // Entropy - 6
//Equation: 58
M_5 = 28.0134 // Molar mass - 6
//Equation: 59
etaT_powerturbine = 0.9// Polytropic efficiency
//Equation: 60
h_6 = calcH_PS("N2";P = p_6;S = s_6) // Enthalpy
// Comment = Polytropic coefficient: k = -Math.log(aval.p/amont.p)/Math.log(aval.V/amont.V)
//Equation: 61
//h_6 = calcH_PS("N2";P = p_6;S = s_6) // Enthalpy
//Equation: 62
T_6 = calcT_PH("N2";P = p_6 ;H = h_6) // Downstream point - 6
// Comment = Given outlet pressure
//Equation: 63
p_6 = 26.0// Outlet pressure
//Equation: 64
W_dot_powerturbine = m_dot_powerturbine*(h_6 - h_5) // DeltaH 

//Process: core compressor
//Equation: 65
m_dot_corecompressor = m_dot_IHXcore // Upstream process - IHX core
//Equation: 66
s_corecompressorinlet = calcS_PH("He";P = p_corecompressorinlet;H = h_corecompressorinlet) // Upstream point - core compressor inlet - Downstream point - core compressor outlet
// Comment = Isentropic reference
//Equation: 67
hs_corecompressoroutlet = calcH_PS("He";P = p_corecompressoroutlet;S = s_corecompressorinlet) // Downstream point - core compressor outlet
//Equation: 68
etaT_corecompressor = 0.9// Isentropic efficiency
//Equation: 69
h_corecompressoroutlet = h_corecompressorinlet + (hs_corecompressoroutlet - h_corecompressorinlet)/etaT_corecompressor // Upstream point - core compressor inlet - Downstream point - corecompressoroutlet
//Equation: 70
T_corecompressoroutlet = calcT_PH("He";P = p_corecompressoroutlet ;H = h_corecompressoroutlet) // Downstream point - core compressor outlet
// Comment = Given outlet pressure
//Equation: 71
p_corecompressoroutlet = 70.0// Outlet pressure
//Equation: 72
W_dot_corecompressor = m_dot_corecompressor*(h_corecompressoroutlet - h_corecompressorinlet) // DeltaH 

//NODES

//HEAT EXCHANGERS

//Heat exchanger: regenerator
//Equation: 73
mCp_HPregenerator = (h_2 - h_1)/(T_2 - T_1)*m_dot_HPregenerator // mCpf =deltaH/deltaT - HP regenerator
//Equation: 74
mCp_LPregenerator = (h_7 - h_6)*m_dot_LPregenerator/(T_7 - T_6) // mCpc =-deltaH/deltaT - LP regenerator
// Comment = mCpc>mCpf
//Equation: 75
UA_regenerator = NTU_regenerator *mCp_HPregenerator // Cold fluid - HP regenerator
//Equation: 76
R_regenerator = mCp_HPregenerator /mCp_LPregenerator // Hot fluid - LP regenerator - Cold fluid - HP regenerator
//Equation: 77
T_2 = T_1 + epsilon_regenerator*(T_6 - T_1) // Hot fluid outlet temperature
//Equation: 78
h_2 =  calcH_TP("N2";T = T_2;P = p_2)// Enthalpy
//Equation: 79
h_7 = - m_dot_HPregenerator /m_dot_LPregenerator*(h_2 - h_1) + h_6 // Hot fluid - LP regenerator - Cold fluid - HP regenerator
//Equation: 80
T_7 =  calcT_PH("N2";P = p_7;H = h_7)// Hot fluid outlet temperature
// Comment = epsilon given value
//Equation: 81
epsilon_regenerator = 0.9 // Given value
//Equation: 82
argLn_regenerator = (1 - epsilon_regenerator*R_regenerator)/(1 - epsilon_regenerator)
//Equation: 83
NTU_regenerator = 1/(1 - R_regenerator)*ln(argLn_regenerator) // Counterflow heat exchanger
//Equation: 84
Q_dot_LPregenerator = m_dot_LPregenerator*(h_7 - h_6) // DeltaH hot fluid
//Equation: 85
Q_dot_HPregenerator = m_dot_HPregenerator*(h_2 - h_1) // DeltaH cold fluid

//Heat exchanger: IHX core
//Equation: 86
mCp_IHX = (h_3 - h_2)/(T_3 - T_2)*m_dot_IHX // mCpf =deltaH/deltaT - IHX
//Equation: 87
mCp_IHXcore = (h_corecompressorinlet - h_coreoutlet)*m_dot_IHXcore/(T_corecompressorinlet - T_coreoutlet) // mCpc =-deltaH/deltaT - IHX core
// Comment = mCpc>mCpf
//Equation: 88
UA_IHXcore = NTU_IHXcore *mCp_IHX // Cold fluid - IHX
//Equation: 89
R_IHXcore = mCp_IHX /mCp_IHXcore // Hot fluid - IHX core - Cold fluid - IHX
//Equation: 90
T_3 = T_2 + epsilon_IHXcore*(T_coreoutlet - T_2) // Hot fluid outlet temperature
//Equation: 91
h_3 =  calcH_TP("N2";T = T_3;P = p_3)// Enthalpy
//Equation: 92
h_corecompressorinlet = - m_dot_IHX /m_dot_IHXcore*(h_3 - h_2) + h_coreoutlet // Hot fluid - IHX core - Cold fluid - IHX
//Equation: 93
T_corecompressorinlet =  calcT_PH("He";P = p_corecompressorinlet;H = h_corecompressorinlet)// Hot fluid outlet temperature
// Comment = epsilon given value
//Equation: 94
epsilon_IHXcore = 0.8462035134515719 // Given value
//Equation: 95
argLn_IHXcore = (1 - epsilon_IHXcore*R_IHXcore)/(1 - epsilon_IHXcore)
//Equation: 96
NTU_IHXcore = 1/(1 - R_IHXcore)*ln(argLn_IHXcore) // Counterflow heat exchanger
//Equation: 97
Q_dot_IHXcore = m_dot_IHXcore*(h_corecompressorinlet - h_coreoutlet) // DeltaH hot fluid
//Equation: 98
Q_dot_IHX = m_dot_IHX*(h_3 - h_2) // DeltaH cold fluid

//Number of equations: 98

//POINTS WITH SATURATION TEMPERATURE SET


//OTHER POINTS WITH PRESSURE SET


//SET FLOW RATES

//Equation: 99
m_dot_core = 180.0// Given flow 
//Equation: 100
m_dot_powerturbine = 700.0// Given flow 

//OVERALL BALANCE

//Equation: 101
useful_Energy = W_dot_HPturbine + W_dot_LPcompressor + W_dot_HPcompressor + W_dot_LPturbine + W_dot_powerturbine + W_dot_corecompressor
//Equation: 102
purchased_Energy = Q_dot_core
//Equation: 103
eta_global = abs(useful_Energy/purchased_Energy)
T_8 = 27.9
p_2 = 70
p_3 = 67.4
p_8 = 25.48
p_7 = 25.74
p_coreoutlet = 67.4
T_coreoutlet = 900
p_10 = 42.47
p_corecompressorinlet = 66.5
T_10 = 27.6