from math import log # ln

n = 1 # arbitraire
R = 8.314
g = 1.4 # air = diatomique
Cp = n*R*g/(g-1)
Cv = n*R/(g-1)

# Lenoir moteur
VA = 1 # arbitraire
PA = 1 # arbitraie
TA = None ### A REMPLACER

VB = VA/2 # enonce
PB = PA # enonce
TB = None ### A REMPLACER

TC = TA # enonce
VC = VB # enonce non utile
#PC = n*R*TC/VC # calcul non utile

QAB = None # isoP
QBC = None # isoV
QCA = None #isoT

eta = None # rendement en fonction de QAB, QBC et QCA
# Verification
print("n moteur Lenoir")
#print(round(eta,6))
print(round(1-g/(1+2*(g-1)*log(2)),6)) # apres simplification



# Stirling pac
alpha = 2 # arbitraire
beta = 3 # arbitraire

V1 = 1 # arbitraire
V2 = alpha*V1 # enonce
T1 = 300 # arbitraire
T2 = beta*T1 # enonce

QAB = None # isoV
QBC = None # isoT
QCD = None # isoV
QDA = None # isoT
Qf = None # analyse a exprimer en fonction de QAB/QBC/QCA/QDA (2 a identifier)
Qc = None # analyse a exprimer en fonction de QAB/QBC/QCA/QDA (2 a identifier)

e = None # efficacite en fonction de Qf et Qc
# Verification
print("e pac Stirling")
#print(round(e,6))
print(round(1/(1+((beta-1)+(g-1)*log(alpha))/((g-1)*beta*log(1/alpha)+(1-beta))),6))



# Atkinson frigo
alpha = 2 # arbitraire
beta = 3 # arbitraire

VA = 1
VB = VA
VD = VA/alpha
VC = VD*beta

TA = 300 # arbitraire
PA = n*R*TA/VA # calcul
TB = 100 # arbitraire
PB = n*R*TB/VB # calcul
# suite de calculs eventuels pour des essais avec valeurs numeriques

print("e frigo Atkinson")
print() # vos calculs pour obtenir e ici
print(round(1/((1/g)*((alpha**g-beta**g)/(alpha-beta))-1),6))